Wave Structures in the Upper Ionospheric Plasma in Spherical Coordinates ()
1. Introduction
The increasing interest in large-scale (planetary) wave disturbances is caused by the fact that many ionospheric phenomena, such as super rotation of the Earth atmosphere [1], ionospheric precursors of natural processes [2] [3] and the ionospheric response to anthropogenic activity [4] [5] fall in the range of such waves. It was assumed that in natural conditions planetary waves generate in the ionosphere-stratosphere and propagate upward to ionospheric altitudes. However, the theoretical study of wave processes as the basis for energy transfers from lower to the upper atmosphere indicates that the system of stable zonal winds reliably screens (especially in summer) the upper atmosphere from the action of large-scale planetary waves generated in the troposphere-stratosphere [6]. Moreover, conditions most favorable for an upward propagation of only very long planetary waves (with a wavenumber of 1 - 2) are created during equinoxes, when the zonal winds change their direction [7].
Nevertheless, a large body of observational data has been accumulated by now [8] [9] [10], which clearly indicate that fast large-scale electromagnetic disturbances are observed at the ionospheric regions [11] [12] [13]. It is also observed that the phase velocities of these disturbances differ by almost an order of magnitude at daily and night conditions. High phase velocities and their considerable variation from day to night make it impossible to identify these disturbances with magnetohydrodynamic waves.
Under natural conditions, these disturbances are observed as background oscillations. The observations indicate that induced oscillations of this type are generated at the impulse excitation of the ionosphere during magnetic storms [2] or under the action on the ionosphere from below as a result of earthquakes, volcanic eruptions or artificial explosions [3] [14] [15] in the latter case; the disturbances appear as solitary wave structure.
In this connection, it becomes extremely necessary to study whether planetary electromagnetic waves can be generated and changed in the ionosphere itself. It is necessary to determine the factors responsible for the generation of such disturbances to interpret adequately the processes proceeding in the ionosphere during wave formation and propagation.
2. Initial Equations
To study the generation of the large-scale low-frequency processes in the ionosphere the magneto-hydrodynamic equations are used for the ionosphere taking into account the fact, that large-scale flows do not perturb density and concentration of the medium particles [16]. The small-scale acoustic gravity and magnetohydrodynamic (Alfven, magnetosonic and helicon) waves are to be altered [17] [18]. The Coriolis force becomes a determining factor for atmospheric processes with a characteristic period of several hours and larger. As far as inhomogeneities of the Earth’s rotation and the magnetic fields are permanently acting fundamental factors for the planetary-scale processes, it would be more appropriate to study the problem of generation of the planetary-scale ultra-low-frequency (ULF) electromagnetic wave structures in F-layer of the ionosphere in a spherical coordinate system. Let
denote an addition to the latitude
,
is a longitude; r—distance from the earth center. The velocity component is assumed to be zero according to r axis—
;
is velocity component along meridian (it is positive if velocity is directed to the north).
is velocity component along the parallels (it is positive if velocity is directed to the east). For simplicity, the geomagnetic field is assumed to have only vertical component
, i.e. the moderate and high latitudes are considered; where
is the geomagnetic field strength near the pole. The angular rotational velocity component
normal to the earth surface is important for the dynamical processes.
The zonal wind determined by the experimental expression
, where
is a constant angular velocity of the atmosphere zonal circulation (so-called circulation index) [11] [12] [13]. The value of
varies seasonably:
—in winter,
—in summer. The stream function can be presented as
. Then, ignoring
in the continuity equation along with incompressibility term we determine the components of the velocity by stream function:
,
, (1)
. (2)
Taking into account above discussions and defined quantities, the close system of dynamical equations for the large-scale wave structures—the equation of motion of the medium particles and induction equation for dissipative ionosphere [16] can be reduced to the form:
, (3)
. (4)
Here, we introduce designations:
(5)
where N is a concentration of the charged particles;
—a neutral concentration; M unit mass of ions and molecules; c is a light speed; R—the earth radius, e—elementary charge. In our task r is a parameter and therefore we shall replace r by the earth radius R in the (3)-(5) taking into account thickness of the atmosphere layer;
is constant coefficient of surface friction between atmospheric layers and is of order of 10−5 s−1 for the ionospheric heights. Further for investigation the large-scale (planetary) zonal flows’ dynamics in F-region of the ionosphere we will use the system of nonlinear Equations (3) and (4).
3. Linear Approximation
We begin an investigation of (3) and (4) from the analysis of motion at the initial stage when the amplitudes of the perturbations are small enough the nonlinear terms to be neglected in (3)-(4). Therefore, we shall seek the solution of these equations in the plane wave form
in linear approximation, where m is a whole number;
-frequency of perturbations and
is an unknown function of
.
Substituting the linear perturbed quantities into (3)-(4) we get a new equation for f:
(6)
The relation (6) presents the equation for the associated Lagrange polynoms. This equation has unique bounded solution if the sum of the second and third terms in the brackets equal to
(where n is a whole number). The latter condition defines a dispersion equation:
, (7)
Submitting eigen-frequency
and damping decrement
;
,
, from (6) the dispersion relation for eigen frequencies:
, (8)
and expression for decrement:
, (9)
will be obtained, where
,
,
.
Let us estimate the order of the coefficients in Equation (8) for the characteristic values of the parameters in F-region of the ionosphere:
,
,
,
,
,
. Thus,
,
,
. It follows, that
,
. Taking all these into account, from (8) we get a spectrum of oscillations:
,
. (10)
We introduce the angular velocity of the wave propagation
, using (10) and determine the linear phase velocity of the wave motion along the latitudinal circles (parallels)—
. Thus, we get the expression for the phase velocity of the waves:
. (11)
This wave is the fast and represents itself a new mode of eigen-frequencies of the oscillations in F-region of the ionosphere. The wave (10), (11) is weakly damped with decrement:
and can propagate along the parallel to the west as well as to the east. The wavelength is 103 km or great.
From induction Equation (4) amplitude of the geomagnetic pulsations generated by these waves,
, where
is a transversal displacement of the charged particles, can be estimated. At displacement
equals 0.1 and 1 km for
,
is equal to 8 and 80 nT respectively.
4. Nonlinear Evolution of the Initial Equations
Investigating the system (3)-(4) at the further stage, when the amplitudes of the perturbations are sufficiently grown, the nonlinearity becomes important In determining the dynamical processes. So, for further analysis of the ionospheric energy transfer processes the close system of nonlinear Equations ((3), (4)) should be considered. For simplicity we consider the motion at any fixed latitude
. It is convenient to introduce the new latitudinal and longitudinal coordinates
,
, the latter transformation causes the freezing of variable coefficients in the system of Equations (3), (4) in non-dissipative approximation (
). We get a closed system of nonlinear equations with constant coefficients solution of which seek in a form of the stationary regular waves
,
, spreading along the latitudes x with velocity U, without changing the shape,
[18]:
, (12)
where F is an arbitrary chosen differential function of its argument,
,
, (13)
where
is the characteristic phase velocity of ordinary Rossby waves.
Following [19] we use the polar coordinates along the earth surface
,
and the circle with the radius a. Further, the quantities
and
is requested to be twice continuously differentiable (including the circle
) along their arguments and to vanish exponentially at
. Then (12) will have the following solutions:
, (14)
(15)
where
is nth degree Bessel function and
—McDonald function; p and
parameters are related by the dispersion relation:
,
, (16)
According the dispersion Equation (16), solution (14) has two free parameters U and a. The disturbed solution obtained from (13) equals to
at
. This indicates that the wave is localized on the earth surface
. This type structures represent a pair of oppositely circulating vortices (cyclone-anticyclone) of equal intensity, spreading along the parallels on the background of purely zonal mean west to east flow. The nonlinear vortex structures move with velocity
. The characteristic scale is
.
The above wave structures will be generated first at the action on the ionosphere from above (e.g. by geomagnetic sub storms) or from below (by earthquakes, major artificial explosions, etc.).
5. Discussion and Conclusions
It has been established that large-scale (planetary) linear and nonlinear electromagnetic wave structures, caused by the permanently acting factor (latitude variation in the geomagnetic field), can be generated in the ionospheric dynamo region. Considering large-scale ionospheric flows with the characteristic extend of 103 km and larger, we cannot ignore nonuniformity of the geomagnetic field and the Earth rotation angular velocity, since the ordinary Rossby parameter
and its magnetic analog
are of the same order of magnitude in the F-region:
,
. The waves have a frequency of the order of
and a wavelength of about 103 km and longer and are almost nondamped
. According to the Maxwell equation for the current
, fast waves generate substantial micro pulsations of the geomagnetic field, which’s order of magnitude (where
is the wavelength of fast waves). Thus,
.
Based on the analytical solutions of the system of nonlinear dynamic Equations (3), (4), it has been shown that the large-scale electromagnetic waves are self-localized in the nondissipative ionospheric F-region in the form of nonlinear solitary vortex structures moving at the constant velocity along latitudinal circles. Nonlinear structures are composed of the pair of oppositely rotating dipole vortices (cyclone-anticyclone) and carry entrapped particles of the medium. The electromagnetic perturbations of slow
and fast
vortices are not less than 103 and 104 km respectively.
According to the equations of transport (3), (47) in the dissipative ionosphere, the energy and enstrophy of large-scale vortices are much greater than the dissipative term, and, correspondingly, such vortices relax very slowly.
Thus, nonlinear large-scale electromagnetic vortex structures are long-lived formations in the ionospheric dynamo region and, therefore, can play a substantial role in the processes of matter, heat and energy transfer and in strong turbulation of the medium [20].
The considered wave structures are the natural degrees of freedom of the ionospheric F-region. Therefore, if the ionosphere is affected from above or below magnetic storms, earthquakes, artificial explosions, etc., the wave structures will first of all be generated and/or enhanced on these models [21]. Solitary vortices will be generated at a certain power of a source [18], which is corroborated by observations [4] [22] [23].
Acknowledgements
The research leading to these results has received funding from the Shota Rustaveli National Science Foundation project N FR17_279.