A&A 386, 653-657 (2002)
DOI: 10.1051/0004-6361:20020279
V. I. Zhukov
Central Astronomical Observatory of Russian Academy of Sciences, 65 Pulkovo, 196140 St. Petersburg, Russia
Received 29 August 2001 / Accepted 29 January 2002
Abstract
We have calculated the spectrum of eigenmodes of umbral
oscillations. The eigenmodes of oscillations with periods from
several tens of minutes (g-modes) to several tens of seconds
(p-modes) have been found. It is shown that the 3-min umbral
oscillations are p-modes modified by the magnetic field.
Key words: Sun: sunspots - Sun: oscillations - magnetohydrodynamics (MHD)
In the beginning of the 70s the discovery of sunspot umbral oscillations (Bhatnagar et al. 1972; Giovanelli 1972; Beckers & Schultz 1972; Bhatnagar & Tanaka 1972) attracted attention of many theorists, as it seemed that the umbral oscillations could provide unique information about the sub-photospheric structure of sunspots.
For years, two basic approaches have been used to study oscillations in sunspots (see the reviews by Staude 1999 and Bogdan 2000). Starting from the study of Uchida & Sakurai (1975), the spectrum of eigenmodes of umbral oscillations has been calculated (Scheuer & Thomas 1981; Thomas & Scheuer 1982; Zhukov et al. 1987; Wood 1990, 1997; Hasan 1991; Cally & Bogdan 1993; Cally et al. 1994; Bogdan & Cally 1997; Gore 1997, 1998; Lites et al. 1998). In the second approach, Zugzda et al. (1983, 1984, 1987); Settele et al. (1999) (see also Gurman & Leibacher 1984) calculated the coefficient of the transmission for a longitudinal wave mode excited by broad-band noise from the deep layers of the convective zone. However, it follows from the paper of Zhukov & Efremov (1988) that these two approaches are essentially equivalent, provided the sunspot umbra contains a resonator from which wave energy can leak. It should be noted, however, that two factors affect the calculated transmission coefficient of the longitudinal wave mode: the partial reflection of the wave and the partial transformation of the longitudinal wave into transverse wave mode. These factors may result in the situation where the transmission coefficient of the longitudinal wave is a monotonous function of frequency (see Fig. 1 in the study of Zhukov & Efremov 1988, and also the study of Settele et al. 1999). In this case, the presence of a resonator for the transverse wave mode in the sunspot umbra can be indicated only by the frequency dependence of the coefficient of transformation of the longitudinal wave mode into the transverse mode (Fig. 4 in the study of Zhukov & Efremov 1988). Thereby, the resonance frequencies cannot always be determined only by calculation of the transmission coefficient of the longitudinal wave mode. That is why in the present study we use the first approach, which while being substantially simpler, unfortunately, does not provide detailed information about linear transformation of waves in a sunspot umbra.
The first and the most consistent attempt to give physical interpretation of the observed oscillations was made by Scheuer & Thomas (1981). They tried to calculate the spectrum of umbral eigenoscillations in sunspots by solving the complete linearized set of equations of magnetohydrodynamics. However, as the authors themselves noted, they were not able to find asymptotic solutions for this set of equations for . In their calculations they used some conditions for sufficiently deep layers of the umbra which in fact did not take into consideration the presence of the convective zone. As a result, Scheuer & Thomas (1981) did not succeed in finding any eigenmodes trapped in the convective zone.
Here we assume that in sufficiently deep layers of umbra, the Alfvén speed is considerably less than the sound speed. In the weak field approximation ( ), the full set of equations of magnetohydrodynamics proved to be divided into two independent subsystems whose asymptotes can be easily obtained (Sect. 3, see also Hasan & Christensen-Dalsgaard 1992 (for the isothermal atmosphere) and Cally & Bogdan 1993; Cally et al. 1994; Bogdan & Cally 1997 (for the polytropic atmosphere)) and, hence, the problem of calculation of the eigenmodes of the umbral oscillations can be completely defined.
The calculations made under the conditions on the infinity (for ) obtained in Sect. 3 have shown that the eigenmodes of the umbral oscillations display periods from several tens of minutes (g-modes) to several tens of seconds (p-modes), and that both the 5-min and 3-min oscillations are eigenmodes of the umbral oscillations.
Let us assume that the magnetic field is homogeneous and vertical
(
). Then for an ideal gas, assuming that
all small perturbations have the form
,
after linearization the basic set of equations of
magnetohydrodynamics is reduced to two independent subsystems, one
of which describes the propagation of Alfvén waves (here it will
not be considered), and the other the propagation of
magneto-atmospheric waves. It has the following form in a
cartesian system of coordinates with the Z-axis directed upwards
(Ferraro & Plumpton 1958)
The system (1)-(2) is a rather complicated set of equations; its solutions have been studied in detail only for an isothermal atmosphere (Ferraro & Plumpton 1958; Lerroy & Schwartz 1982; Zhugzhda & Dzhalilov 1982).
In our opinion, at the present stage, until the origin of the umbral oscillations has been ultimately determined, the problem should not be complicated with attempts to describe the structure of the atmosphere of a sunspot umbra with all possible accuracy (for example, with the use of semi-empirical models). Therefore, in the present study we solve the set of Eqs. (1)-(2) numerically for the model of umbra similar to that used by Scheuer & Thomas (1981). We assume that our umbra model is a two-layer atmosphere. The sound speed and density in each layer are given as follows:
Layer 1 (Chromosphere - corona, z>0)
Layer 2 (Convective zone, z<0)
In most studies, zero boundary conditions were used in the calculation of the spectrum of eigenmodes of umbral oscillations. In fact, however, (see Zhukov & Efremov 1988; Cally & Bogdan 1993; Cally et al. 1994; Bogdan & Cally 1997; Lites et al. 1998; Settele at al. 1999), in order to correctly carry out the numerical integration of the set of Eqs. (1)-(2), it is necessary to know the asymptote of its solutions both in the corona (at ) and convective zone (at ).
As it was mentioned above, the solution of the system of Eqs. (1)-(2) is well known for the isothermal
atmosphere; in particular, the asymptotic solutions of this system
at
which are necessary for our purposes have the
form (Leroy & Schwartz 1982):
Let us find asymptotic solutions of the system
(1)-(2) in the deep layers of umbra (
). As the sound speed and density increase with depth, it
is obvious that in sufficiently deep layers of the umbra, the
Alfvén speed becomes far smaller than the sound speed.
Therefore, let us introduce a small parameter
,
where
is the depth, at
which
.
Then seeking the solutions of this system
in the form of the series:
In the second case, the following equation for
is
obtained:
Now, for the problem to be formulated completely it is necessary
to find a dispersion relation which specifies the spectrum of
umbra eigenoscillations. For this purpose, conditions on the
boundary of discontinuity at z = 0 must be found. The conditions
on the boundary of the discontinuity in a medium with a vertical
magnetic field has been considered in a number of studies (see,
for example, Leroy & Schwartz 1982); they have the form
Taking into account the fact that in the corona (at
)
the exponential decaying solution (v_{z1t}) and the
solution corresponding to upwards propagated slow magnetoacoustic
waves (v_{z2t}) should be taken, for layer 1 (z>0) we have
In deep layers of the convective zone it is necessary to take the
exponential decaying solution of the Eq. (6) (v_{z1b})
and the solution for downwards propagated slow magnetoacoustic
waves (5) (v_{z2b}) and for the layer 2 (z<0) we
have
Thus, from the boundary conditions (7) the following set
of equations for the definition of arbitrary constants
T 1, T _
2, B 1 and B 2 is obtained:
Figure 1: Variation of the sound speed with height in the umbra. | |
Open with DEXTER |
Figure 2: Variation of the density with height in the umbra. | |
Open with DEXTER |
To compare our results with observations, we adopt the following
numerical values for the parameters in our umbral models:
,
,
,
,
,
,
.
The variation of the sound
speed and density with height for the chosen values of the
parameters and for
is presented in
Figs. 1 and 2.
Figure 3: Diagnostic diagram for the sunspot umbral oscillations. | |
Open with DEXTER |
Figure 4: Diagnostic diagram for the umbral oscillations without magnetic field. | |
Open with DEXTER |
The diagnostic diagram for the umbral oscillations presented in
Fig. 3 was calculated for
Figure 5: Variations of frequencies of the f and p_{n}- modes () with the magnetic field for . | |
Open with DEXTER |
The calculated eigenperiods of the umbral oscillations (see
Fig. 3) for two values of the wave number (i.e., in fact,
for two sunspots of different size) and the observed periods of
the oscillations taken from the study of Thomas et al.
1984) are presented in Table 1. We can see from
Table 1 that despite the rather arbitrary choice of some parameters
for the sunspot umbra in our calculations (e.g., the parameter
,
which essentially specifies the height of the transition
region, or the magnitude z_{2} which, as our calculations have
shown, strongly influences the spectrum of oscillations), the
calculated values of the periods are generally close to those
observed. It thereby follows from the obtained results that both
the 5-min and 3-min oscillations in the sunspot umbra are the
eigenmodes of umbral oscillations. It is necessary to note,
however, that the 3-min umbral oscillations are high-frequencies
p-modes (with n> 4) appreciably modified by the magnetic
field.
mode | Obs. | ||
(s) | (s) | (s) | |
g_{2} | 51.5(min) | 44.0(min) | - |
g_{1} | 28.4(min) | 24.9(min) | - |
f | 392 | 352 | 366 |
p_{1} | 246 | 225 | 301 |
p_{2} | 201 | 195 | 270 |
p_{3} | 186 | 172 | 197 |
p_{4} | 168 | 154 | 171 |
p_{5} | 164 | 139 | 155 |
It follows from the asymptotic solutions of the system (1)-(2) (Sect. 3) that energy can leak from the sunspot umbra to the corona and deep layers of a sunspot in the form of slow MHD waves. The radiation of the waves to the corona and deep layers of a sunspot implies that the eigenfrequencies of the umbral oscillations must be complex ( ) (see Zhukov & Efremov 1988; Cally & Bogdan 1993; Bogdan & Cally 1997), and therefore the results obtained in the present paper are reliable provided that , i.e. the leakage of the wave energy from the umbral resonator is sufficiently small. In this case, as shown for example by Zhukov (2001), the leakage of the wave energy from the umbral resonator ( ) should result in minor change of , which as a first approximation can be neglected. Moreover, the initial assumption accepted while calculating eigenfrequencies for g, f and p-modes can be shown to be actually valid, taking into account conditions at infinity ( ) derived in Sect. 3. However, it follows from the study of Zhukov & Efremov (1988), in which the spectrum of eigenoscillations of the sunspot umbra was calculated (in the approximation of an incompressible medium, which, generally speaking, is applicable only to the f-mode), that the radiation of the wave energy from sunspots can be substantial (see also Cally & Bogdan 1993; Bogdan & Cally 1997) and therefore it is very important to take into account the leakage of wave energy from the sunspot umbra in calculations. In addition to that, in order to achieve the best consistency between calculations and observations, it is necessary to calculate the spectrum of the umbral oscillations for the various z_{2} which dramatically affect both f and several first p-modes.
We note in conclusion, that since open boundary conditions are acquired for the umbra, the umbral oscillations (at least within the 3-min band) may be excited by the mechanism suggested by Zugzda et al. (1983). We consider, however, that generation of oscillations in sunspots constitutes a rather complicated problem, which will require substantial effort in order to be resolved.
Acknowledgements
I would like to thank the referees for useful suggestions that helped improve the manuscript.