Next Article in Journal
Visualization of Prediction Methods for Wildfire Modeling Using CiteSpace: A Bibliometric Analysis
Next Article in Special Issue
Atmospheric and Ionospheric Effects of La Palma Volcano 2021 Eruption
Previous Article in Journal
Comparison of Lichens and Mosses as Biomonitors of Airborne Microplastics
Previous Article in Special Issue
Propagation of Perturbations in the Lower and Upper Atmosphere over the Central Mediterranean, Driven by the 15 January 2022 Hunga Tonga-Hunga Ha’apai Volcano Explosion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Role of Neutral Wind Velocity and Its Vertical Component on Predictability of Formation and Localization of Sporadic E (Es)

by
Goderdzi G. Didebulidze
*,
Giorgi Dalakishvili
,
Maya Todua
and
Lekso Toriashvili
Space Research Center, School of Natural Sciences and Medicine, Ilia State University, Kakutsa Cholokashvili Ave. 3/5, 0162 Tbilisi, Georgia
*
Author to whom correspondence should be addressed.
Atmosphere 2023, 14(6), 1008; https://doi.org/10.3390/atmos14061008
Submission received: 11 April 2023 / Revised: 11 May 2023 / Accepted: 8 June 2023 / Published: 10 June 2023
(This article belongs to the Special Issue State-of-the-Art in Gravity Waves and Atmospheric-Ionospheric Physics)

Abstract

:
The formation and localization of sporadic E (Es) layers predicted by the ion vertical drift velocity and its vertical change in the lower thermosphere during nighttime are shown analytically and by numerical simulations. The consideration of the existence of a minimum negative value of the vertical change of the ion vertical drift velocity as a necessary condition of formation of the Es layer and determining ion convergence rate into this layer is extended in case of the presence of vertical wind. Upward vertical wind can shift the convergence regions upward, while downward vertical wind shifts them downward, unlike the cases of the presence of only meridional and zonal winds. It also changes the ion convergence rate compared to the one with just a horizontal wind. It is shown for the first time that the upward constant wind also causes the convergence of ions with the maximum rate in the region where the ion-neutral collision frequency is equal to their cyclotron frequency. While demonstrating the presented theory by numerical simulations, HWM14 data is used for the meridional and zonal wind velocity profiles and the presence of its vertical component is assumed. In this case, in addition to the estimated ion drift velocity and its vertical changes, their initial distribution and ambipolar diffusion also determine the development of ion convergence/divergence processes. For a small magnitude of vertical wind velocity, its significant influence on the ions/electrons behavior is demonstrated, which shows the importance of both tidal wind and wind changes caused by the propagation of atmospheric gravity waves on the formation of the Es layer. In this theoretical research, it is shown that the realistic profile of the wind velocity, which takes into account the vertical component along with its zonal and meridional ones, is important for the prediction of the Es layers formation, as well as regions of ion depletion.

1. Introduction

Sporadic E (Es) is a high-density narrow and irregular ionosphere layer [1,2], and its existence was noticed in the ionograms of 1933–1935 [3]. In this layer, the ion/electron density can be comparable to or greater than those of the regular ionosphere E, F1 and F2 layers. Therefore, Es may influence the quality of radio communications [4,5].
Es layers occur during nighttime and more frequently during daytime in the lower thermosphere, at altitudes between 90–150 km [3,6,7]. This fact indicates that in this range of altitudes, the conditions for ion/electron convergence into a narrow layer are provided [8,9,10,11,12,13,14]. This convergence and, consequently, the formation of the Es layer can be caused by the horizontal wind with vertical shear (windshear theory) [8,9,15,16,17,18,19,20], by an electric field and its variations [21,22] and/or by the combined action of all these effects [23,24,25]. By using the set of magnetohydrodynamic equations for the mid-latitude lower thermosphere, the ion vertical drift velocity can be estimated by the wind velocity, the vertical distribution of this region’s dominant neutral particles and the geomagnetic field. By this drift velocity and its vertical changes, the height regions can be estimated where the ion/electron convergence or depletion occurs [12,14,26]. The loss rate of molecular ions (e.g., NO+ and O2+) increases during ion/electron convergence into a thin dense layer, which in turn causes the dominance of long-lived heavy metallic ions (e.g., Fe+ and Mg+), which is a phenomenon that has been observed in nighttime Es layers [27,28].
In recent studies [12,26] it was shown that in addition to earlier considerations [2,3,15,29,30], which took into account the “windshear” effect, the wind magnitude and direction were important for the development of ion/electron convergence (or divergence/depletion) in certain height regions. It was also shown that Es layer density is controlled by ion ambipolar diffusion as well.
In the present study, the suggested theoretical mechanism of the Es layer formation takes into account an arbitrary wind height profile. This approach is different from the earlier considerations [12,26]. In this study, the recently obtained results of the case of horizontal wind are also generalized. This makes the formation of the Es layers more predictable. These processes can also be described by an analytical approach. The behavior of the vertical distribution of ion density N i ( z , t ) in an ambipolar approach and correspondingly the formation of ionosphere layers, similar to [12,26,31] (with some extension, including gravity effect as well), can be described by the following equation:
N i t = w i v z N i w i v ( z ) N i z + z D a ( z ) N i z + N i H i + Q L
Here, in the expression of ion vertical drift velocity w i , the w i v part caused by neutral wind and the part D a ( z ) N i z + N i H i caused by ions ambipolar diffusion ( D a ) are separated. This separation makes further analysis more convenient. D a is the coefficient of ion ambipolar diffusion. Q and L are the ion production and loss rates, respectively. H i is the ion scale height. The effect of gravity in the lower thermosphere (h < 140 km) is relatively small ( N i z > > N i H i ). Here and hereafter, the vertical coordinate z represents the residual between some actual h and initial heights h o ( z = h h o ).
Despite the fact that ionosphere plasma is weakly ionized and the influence of neutrals on the ion/electron drift velocity and correspondingly on their density behavior is important, the formation of the regular ionosphere E, F1 and F2 layers is still possible, even in the absence of neutral wind ( V = 0 and w i v = 0 ). For example, in the lower thermosphere, the formation of the E layer by the Chapman mechanism is mainly determined by the photoionization rate ( Q ) of molecules (e.g., O2 and N2) with vertically decreasing density [5].
For the formation of the F1 layer during daytime (mostly in the summer season, at a height of about 180–250 km), the photoionization of atomic oxygen is important, which decreases in density with height, since the electron production is dominated by extreme ultraviolet (UV, 10–100 nm) radiation ionizing atomic oxygen. During the night, the photoionization is absent ( Q = 0 ), so the F1 layer vanishes, and the F2 becomes a dominant layer. This layer regularly exists (with peak electron density at about 300–400 km height) and it is maintained by the increase in ion ambipolar diffusion displacement ( D a ( z ) N i z + N i H i ) above the electron density peak height [5]), as well as by the increase in ion recombination below this height. Here, these processes are slower than in the E and F1 regions, allowing the F2 layer to exist during both daytime, when photoionization increases ( Q > 0 ), and nighttime, when ionization is absent ( Q = 0 ). Although the neutral wind is not necessary for the formation of the E, F1 and F2 layers, it can significantly influence ion/electron vertical drift described by w i v , as in Equation (1), and therefore affect the behavior of these layers.
In the suggested study, unlike the regular ionosphere E, F1 and F2 layers, the formation of the sporadic E (Es) layer at mid-latitude lower thermosphere is mainly determined by the neutral wind. The ion vertical drift velocity, caused by the combined effect of vertically changing ion-neutral collision and Lorentz forcing, changes during the existence of homogeneous or inhomogeneous wind ( w i v z 0 ). At some height regions, this drift velocity is zero ( w i v = 0 ), depending on the direction of the wind. These height regions w i v = 0 will be noted as a node of ion drift velocity. When the ion upstream and downstream flows ( ( d i v w w i v z ) min < 0 ) towards this region occur, then their convergence into a thin Es-type layer is possible at this node [26]. The processes of ions convergence into a thin layer should take place until the increase in the density due to the convergence exceeds its decrease due to diffusive displacement ( D a ), and total loss ( Q L < 0 ) caused by their production ( Q ) and loss ( L ) rates.
For those wind directions when the ions downstream and upstream flows occur ( ( d i v w w i v z ) max > 0 ), the opposite phenomenon takes place in these height regions—the ion/electron depletion. Further, we will show that both horizontal and vertical components of the wind velocity are important for the estimation of the ion/electron convergence/divergence regions and prediction of the development of corresponding processes.
To show the suggested theoretical mechanism of sporadic E formation, we will solve Equation (1) using an analytic approach when
w i v z C = c o n s t a n t
and D a = c o n s t . In this case, the solution of Equation (1) has the following form:
N i ( z , t ) = e C ( t t o ) 2 [ π D a ( t t o ) ] 1 / 2 + e [ z w i v ( t t o ) z ] 2 4 D a ( t t o ) N i o ( z ) d z
Here,
G i ( z , t ) = e C ( t t o ) 2 [ π D a ( t t o ) ] 1 / 2 e [ z w i v ( t t o ) ] 2 4 D a ( t t o )
is the analogy of the source function [32], which corresponds to the solution of Equation (1) in the case of Q L = 0 . When C > 0 , then G i ( z , t ) increases with time, and when C < 0 , it decreases. Equation (3) also shows that for ion vertical drift velocity w i v > 0 , the initial peak height z m ( t = t o ) = 0 of G i ( z , t ) shifts upward and for w i v < 0 it descends. We note that the upward shift of the ion’s initial layer, w i v , is not noticeable when ion diffusion is neglected.
In Equation (3), the value + e [ z w i v ( t t o ) z ] 2 4 D a ( t t o ) N i o ( z ) d z indicates that, in addition to the diffusive displacement of the ion’s initial distribution N i o ( z ) , their redistribution also occurs in the height region of its drift velocity node, w i v = 0 .
When the increase in ions density ( e C ( t t o ) > 1 ), caused by vertical changes of ion vertical drift velocity ( C > 0 ), exceeds the diffusive displacement of their initial density ( + e [ z w i v ( t t o ) z ] 2 4 D a ( t t o ) N i o ( z ) d z ), then the density in the layers formed in the nodes increases; see Equation (3). Equation (3) also shows that when C < 0 and, correspondingly, e C ( t t o ) < 1 , then in this region the opposite phenomenon takes place—its density decreases (depletes).
Actually, C = w i v z changes with height and their influence on ion/electron density variations will differ across various height regions. If a certain height C has a maximal positive value ( C max > 0 ), then, considering ion/electron density behavior approximation by Equation (3), ions/electrons will be accumulated into a thin layer in the vicinity of this height, and the Es-type layer formation is expectable.
Here and hereafter, the superscript prime denotes the vertical derivative of the corresponding function, and the subscript max or min denotes its minimum and maximum values, respectively.
Using the approximated value of w i v ( z ) at the region z = z o , it can be shown that
w i v w i v ( z = z o ) + ( z z o ) w i v z z = z o + 1 2 ( z z o ) 2 2 w i v z 2 z = z o
where C = w i v z z = z o is maximal ( C max > 0 ) or minimal ( C min < 0 ), then 2 w i v z 2 z = z o = 0 .
Equations (2)–(5) show that when C max > 0 , then the ions flow upstream N i w i v z < z o ( z z o ) w i v z z = z o > 0 and downstream N i w i v z > z o ( z z o ) w i v z z = z o < 0 to the region z = z o , where their density increases maximally ( exp [ C ( t t o ) ] > 1 ), as in Equations (3) and (4). So, the formation of a high-density Es-type narrow layer occurs at the region of the ion vertical drift velocity node z = z o ( w i v = 0 ).
We note that this process can be developing along with the vertical drift of ions, during which the formed layer moves or is localized in the region where w i v 0 . If the height region with w i v = 0 is shifted from the region with C > 0 , then, according to Equations (3) and (4), the ion vertical convergence at this region is expectable during ions’ upward ( w i v ( z = z o > z c ) > 0 ) or downward ( w i v ( z = z o > z c ) < 0 ) drift. Equation (4) also shows that when C < 0 , then at the z = z o region, the opposite effect occurs: the ions downstream N i w i v z < z o ( z z o ) w i v z z = z o < 0 and upstream N i w i v z > z o ( z z o ) w i v z z = z o > 0 flows increase, resulting in the depletion of ions density at regions with C min < 0 , as in Equations (4) and (5).
We note that, for certain directions of the horizontal wind, vertical changes of the ion vertical drift velocity already have minimal negative ( C max > 0 ) or/and maximal positive ( C min < 0 ) values, even if the wind is homogeneous [12]. Therefore, at height ranges of 90–150 km, the processes of ions convergence into a thin layer or their depletion are possible. The formation of the Es-type layer (at a region with w i v = 0 ) occurs until its density increase is balanced by the increased ions diffusive outflow from this layer ( e C ( t t o ) = 1 2 [ π D a ( t t o ) ] 1 / 2 + e [ z w i v ( t t o ) z ] 2 4 D a ( t t o ) N i o ( z ) d z ).
It is also important to note that heavy ions are mainly distributed between about 90–150 km [33,34] and correspondingly the increase in density of the Es layer continues during its upward w i v + 0 (towards the topside) or downward w i v 0 (towards the bottom) motions where vertical changes of the vertical drift velocity vanish ( w i v z 0 ).
The wind in the lower thermosphere can be caused by vertically propagating tidal motion, planetary waves and/or atmospheric gravity waves (AGWs). In this case, the ion vertical drift velocity can have several nodes where the ions’ convergence (when C max > 0 ) or divergence (when C min < 0 ) is possible. So, during the presence of more than one node with the condition C max > 0 , the formation of multilayered sporadic E is expectable [26].
So, on the basis of the realistic height profile of the neutral wind velocity and the ions’ initial distribution, the prediction of the formation of the Es-layer in the mid-latitude lower thermosphere is possible. Under the same conditions, the estimation of the ions depletion region is also possible. In these cases, both horizontal and vertical components of the wind profile are important, which will be shown in the suggested study.
Hereafter, for demonstration of the above-derived analytical theory of the formation of sporadic E and ions density depletion under the influence of neutral wind, we will solve Equation (1) numerically. In the equation for ions’ vertical drift velocity, w i the neutral wind effect ( w i v ) will be distinguished from its drift by ions’ ambipolar diffusion. The convergence/divergence rate parameter w i v z C = c o n s t a n t will be generalized for a variable wind velocity profile and the equations describing the role of (1) its magnitude and direction, (2) the vertical changes of the horizontal components, and (3) the vertical changes of the vertical component will be separated. Unlike earlier considerations [12], the influence of vertical wind on both the convergence of ions and their depletion in a region of a certain height of the lower thermosphere will be shown.
In the first stage, we show only the effect of vertical wind on the convergence/divergence of ions. Further, we will consider its existence in the horizontal wind, whose meridional and zonal velocity profiles are determined by HWM14 data [35]. As an example, in numerical simulations, we use the HWM14 data for the mid-latitude (45° ± 2° N, 45° ± 2° E) lower thermosphere region at nighttime (from LT 23:00) around the spring equinox (e.g., day 102).
We will choose a small magnitude of vertical wind velocity, which can be caused by both tidal motions and AGWs often observed in the lower thermosphere [36,37]. We will show its influence on the location of the Es layer (e.g., the location with w i v = 0 or w i v 0 ), as well as the presence of convergence (with w i v / z min < 0 ) and divergence (with w i v / z max > 0 ) height regions, and the development of corresponding ion/electrons density convergence/depletion processes in these regions. For the considered wind velocity profiles, the location of the height regions of maximum convergence and divergence rates, as well as the parameters caused by the magnitude, direction and vertical change of the wind, will be demonstrated. Theoretical and corresponding numerical simulations show that in the mid-latitude regions the formation of the sporadic E (Es), along with the significant effect of vertical shear of the horizontal wind (shear theory) [38,39], the factor of the magnitude and direction of the three-dimensional wind is important.
Sporadic E contains mainly metallic ions and, consequently, it is important to know the realistic wind velocity height profile, to study their global and regional distributions.

2. Parameters Describing the Formation of Sporadic E under the Influence of Neutral Wind

In this section, we consider the behavior of ions/electrons and the possibility of sporadic E formation in the mid-latitude lower thermosphere under the influence of arbitrarily directed neutral wind.
In this case, we use Equation (1) to study the behavior of heavy metallic ions (Fe+) and electrons under the assumption of quasi-neutrality ( N i = N e ), where the ions’ vertical drift velocity w i v (determined by the wind) and their vertical changes w i v z vary in time. Here, along with the effect of wind shear, the influence of its magnitude and direction on the process of ions/electrons convergence into a narrow Es-type layer (convergence instability) and the importance of balancing this convergence by ambipolar diffusion will be taken into account. We also solve the ion/electron continuity equation with an analytical approach, which, similar to Equations (2)–(4), shows the formation of the Es layer and the tendency for its localization in the region of nodes of ion drift velocity w i v . We use the analytical profiles of ion/electron density to set initial conditions, and then solve their basic equation, which takes into account the presence of wind with arbitrary direction, using numerical methods in the following cases: (1) in the presence of vertical wind only; (2) in the presence of horizontal wind whose meridional and zonal components are determined by HWM14 data; and (3) in the presence of wind determined by the above-mentioned vertical and zonal components. In these cases, the process of formation and localization of the Es layer throughout the night (covering 12-h intervals) will be shown. Assuming the presence of a vertical wind component, we demonstrate its significant influence on the formation and localization of the Es layer. The demonstrated numerical results align with its analytically predicted behavior. Additionally, we will show the impact of vertical wind on the ion/electron depletion process.
The neutral particles with velocity V ( V x , V y   V z ) in the lower thermosphere influence the ion vertical drift via the combined action of the Lorentz forcing and ion-neutral collision. We take a right-handed set of coordinates (x, y, z) with x directed to the magnetic north, y to the west and z vertically upward. In this case, the dependence of the ions’ vertical drift velocity w i on the neutral wind velocity, taking into account their ambipolar diffusion ( D a ), can be described by the equation [22,29,30]:
w i = w i v D a ( z ) 1 N i N i z
where the influence of neutral wind velocity V ( V x , V y   V z ) on the ion vertical drift velocity w i is considered in the value of w i v , which has the following form [2,23,30]:
w i v ( z , t ) = C x ( z ) V x ( z , t ) C y ( z ) V y ( z , t ) + C z ( z ) V z ( z , t )
Here, considerations of the recent studies [12,26] are extended and Equation (7) for w i v takes into account not only the horizontal V x (meridional) and V y (zonal) components of neutral wind velocity, but also the influence of its vertical V z components on ion vertical drift velocity is included.
In Equations (6) and (7):
C x ( z ) = 1 1 + κ 2 ( z ) sin I cos I
C y ( z ) = κ ( z ) 1 + κ 2 ( z ) cos I
C z ( z ) = sin 2 I + κ 2 ( z ) κ 2 ( z ) + 1
D a ( z ) = κ 2 ( z ) + sin 2 I 1 + κ 2 ( z ) 2 k T m i ν i n ( z )
Here, κ ( z ) = ν i n ( z ) / ω i , ν i n is the ion-neutral collision frequency, whose height profile will be used similarly to the one applied in [12,40]. ω i = e B / m i is the ion gyrofrequency; B ( B x = B cos I ,   0 ,   B z = B sin I ) is the Earth’s magnetic field vector, and I is the magnetic dip angle. T = ( T e + T i ) / 2 is the mean plasma temperature; T i and T e are the ion and electron temperatures, respectively; k B is the Boltzmann constant; m i is the ion mass.
To investigate the behavior of the height profile of electron/ion density N e ( z , t ) (described in the previous section) in more detail, assuming quasi-neutrality N e = N i , neglecting production Q and loss rates L, and using ions vertical drift velocity, as in Equations (6) and (7), we rewrite the continuity Equation (1) as follows:
N i t = C t N i w i v N i z + z D a N i z
where
C t ( z , t ) w i v z = C v ( z , t ) + C s h ( z , t ) + C v z ( z , t )
Here, in the value of ions/electrons total convergence C t ( z , t ) > 0 ( w i v z < 0 ) or divergence C t ( z , t ) < 0 ( w i v z > 0 ) rate, Equation (13), we have separated the terms depending on: (1) the effects of wind velocity direction (with meridional V x ( z , t ) , zonal V y ( z , t ) and vertical V z ( z , t ) components) and magnitude described by C v ( z , t ) ; (2) the effects of meridional and zonal wind velocity vertical shear ( V x ( z , t ) z and V y ( z , t ) z ) described by C s h ( z , t ) ; and (3) the effects of vertical changes of wind velocity vertical component ( V z ( z , t ) z ) described by C v z ( z , t ) . In this case, see Equation (7), as C v ( z , t ) , C s h ( z , t ) and C v z ( z , t ) have the following forms:
C v ( z , t ) = V x ( z , t ) C x ( z ) z + V y ( z , t ) C y ( z ) z V z ( z , t ) C z ( z ) z
C s h ( z , t ) = C x ( z ) V x ( z , t ) z + C y ( z ) V y ( z , t ) z
C v z ( z , t ) = C z ( z ) V z ( z , t ) z
Equations (14)–(16) show that the neutral meridional V x ( z , t ) , zonal V y ( z , t ) and vertical V z ( z , t ) wind velocity factors C x , C y , C z of ion vertical drift velocity w i v , Equation (6), as well as their vertical changes C x ( z ) z , C y ( z ) z , C z ( z ) z significantly determine the condition of the ion/electrons convergence ( ( w i v z ) min < 0 ) or divergence ( ( w i v z ) max > 0 ) caused by the magnitude and direction of velocity V and its vertical changes ( V x ( z , t ) z , V y ( z , t ) z and V z ( z , t ) z ). The role of ambipolar diffusion of ions ( D a N i z ), as in Equations (6) and (11), is also important in the development of these processes.
In this case, the tendency of ion/electron convergence into a thin layer when the minimal negative value in w i v z (maximal value of C t o t ( z , t ) ) is present, also can be seen and analyzed from the analytic approach of solution of Equation (12), which has the following form:
N e ( z , t ) N o m exp 2 D a H i c 2 + C t ( z , t ) ( t t o ) z [ z o m + w i v ( t t o ) ] H i c 2
Here, the slight changes in C t ( z , t ) over time t t o < < H i c 2 / 2 D a are assumed.
H i c is the characteristic scale height of ions, which at some initial time t = t o determines the main ion/electron layer thickness ( 2 H i c ) and the height region z z o m = ± H i c , where their density decreases by a factor of e. H i c 2 2 D a is the characteristic time of the ion density decrease by their ambipolar diffusion. The behavior of ion/electron density described by Equation (17) is similar to the one described by Equation (3) under the assumption of Gaussian-type distribution, at the initial time t = t o , with the maximal density (peak density) N o m at the height z = z o m (initial peak height). When C t ( z , t = t o ) > 2 D a H i c 2 , then Equation (17) shows a tendency for the formation of ions/electrons high-density sporadic Es-type layer ( N m ( t > t o ) N o m > 1 ). The peak height of this layer varies in time as follows:
z = z o m + w i v ( t t o )
i.e., moves upward ( w i v > 0 ) or downward ( w i v < 0 ) by ions’ vertical drift velocity w i v and it could be localized at the height region with w i v = 0 (the drift velocity node). When C t ( z , t = t o ) < 0 occurs at the region of the ion drift velocity node, ion/electron density experiences depletion ( N m ( t > t o ) N o m < 1 ), as predicted by Equation (17).
Thus, the behavior of the ion/electron density profile described by Equation (17), which is similar to [12,26], shows the tendency of Es layer formation when considering both horizontal and vertical winds. This solution takes into account the effect of ions’ ambipolar diffusion on the density of their layer, as well as its vertical drift, as seen in Equation (18). In this case, the velocity factor C z ( z ) and its vertical change profile C z ( z ) z , as seen in Equation (15), significantly determine the effect of the vertical wind on the ions’/electrons’ vertical drift, as in Equation (7), and on their convergence into an Es type layer. These quantities, as well as the meridional C x ( z ) and zonal C y ( z ) velocity factors, their vertical changes C x ( z ) z and C y ( z ) z , the meridional and zonal wind vertical shears, and the vertical changes of vertical wind, determine the minimal negative value of vertical changes of the ion vertical drift velocity, as in Equation (13), which is necessary for the formation of the Es-type layer, as in Equation (17). Note that, actually, Equation (17) gives continuous spatial and temporal variations of the ion/electron density height profile for the time of t t o < < H i c 2 / 2 D a . However, we use them in accordance with the space and time resolutions determined by the neutral atmosphere (e.g., NRLMSISE-00) and horizontal wind data (e.g., HWM14).
Figure 1 demonstrates the features of the height profiles of C x ( z ) , C y ( z ) and C z ( z ) (panel a) and, correspondingly, (panel b): C x ( z ) z with a maximal value of about 2.41 × 10−5 s−1 at height ≈ 121 km, C y ( z ) z with a maximal value of about 1.72 × 10−5 s−1 at height ≈ 115 km and minimal value −0.94 × 10−5 s−1 at height ≈ 131 km, C z ( z ) z with maximal value 1.23 × 105 s−1 at height ≈ 121 km (panel b). These features lead to the presence of a minimal negative value of vertical changes of ions’ vertical drift velocity at lower thermospheres, even in cases wherein neutral wind velocity is constant. For example, in the case of only vertical wind ( V x = 0 ,   V y = 0 ,   V z 0 ), the vertical changes of the ion vertical drift velocity ( w i v z V z C z z ) have minimal negative value at the height of about 121 km, and the formation of the Es-type layer for the upward vertical wind ( V z > 0 ) is possible; see Equations (14), (17) and (18).
We note that for the isothermal atmosphere, using a barometric distribution of the neutral particle densities exp ( z / H ) and, correspondingly, ion-neutral collision frequency ν i n ( z ) exp ( z / H ) , for the minimal values of vertical changes of the vertical wind velocity factor C z , we get
( C z z ) max C z max = cos 2 I 2 H
The above analysis and Equations (17)–(19) show that the upward wind V z > 0 can cause additional convergence at the height region ν i n = ω i in case of the presence of northward wind V x > 0 . It also will influence the expectable development of ion/electron depletion in this region caused by southward wind. It is also expectable the influence of vertical wind on Es-type layer location due to its significant effect on ion/electron vertical drift velocity. For example, Equations (6), (9) and (14) show that, to the bottom of the lower thermosphere h < 110 km where the ion-neutral collisional frequency exceeds the ion cyclotron frequency ( κ > > 1 ), the meridional C x ( z ) and zonal C y ( z ) wind velocity factors vanish, but the vertical wind velocity factor C z ( z ) 1 . In this case, w i v ( z , t ) V z ( z , t ) , the initial layer descends ( w i v < 0 ) or moves upward ( w i v > 0 ) by vertical wind velocity V z . Additionally, for the upper altitudes h > 130 km, where the zonal wind factor disappears, and in the absence of meridional wind, the vertical drift of ions/electrons is determined by the vertical wind speed w i v ( z , t ) sin 2 I V z ( z , t ) .
In the next section, we will show numerically the behavior of ion/electron density described by Equation (17) under the influence of neutral wind. In the first stage, we will assume the presence of vertical wind only, then we will show its effect in case of the presence of the meridional and zonal winds with velocities specified by the data of the HWM14 [35]. In these cases, the formation of an Es-type layer (when w i v z max C t ( z , t = t o ) max > 2 D a H i c 2 ) and the localization of ions in the region of their vertical drift velocity nodes (for example, at height with w i v = 0 ) or regions where w i v 0 , predicted by Equations (17) and (18), will be shown throughout the night ( t t o 12   h ). Equations (14)–(16) will be used to estimate the height regions of total maximum convergence ( C t max ( z , t ) = C v ( z , t ) + C s h ( z , t ) max > 0 ) and divergence ( C t min ( z , t ) = C v ( z , t ) + C s h ( z , t ) min < 0 ) of ions and, accordingly, the parameters C v ( z , t ) and C s h ( z , t ) , which depend on the magnitude, direction and vertical shear of the neutral wind speed, respectively, and determine the above-mentioned regions.

3. Numerical Results and Discussion

In order to assess the role of vertical wind in the process of ion convergence, as the first step, we present (numerically) the behavior of ions/electrons, described by Equation (12), in the nighttime conditions, in the cases of constant vertical wind directed upward ( V z > 0 ) or downward ( V z < 0 ). Similarly, we consider the behavior of ions/electrons density and the possibility of Es layer formation under the influence of horizontal wind, with the meridional V x and zonal V y components of the velocity determined by HWM14 data. Further, in order to illustrate the influence of the vertical component of the three-dimensional (arbitrary direction) wind on the ions/electrons density and accordingly, on the formation of the Es layer, we will incorporate the vertical component in the model horizontal wind. As an example, we will choose the vertical wind velocity magnitude of V z = ± 4   m / s . This value in the lower thermosphere is about an order of magnitude lower than the typical values of meridional and zonal wind speeds caused by tidal motions (about 30–80 m/s) [41]. The consideration of the effect of even such a small value in the vertical wind velocity should demonstrate its importance in ions’ vertical drift and their convergence into the Es-type layer. Neutral wind velocity vertical components of this magnitude and higher can be caused by both tidal motions and the propagation of AGWs frequently observed in the lower thermosphere and mesopause regions [36,37].
We will perform numerical simulations for the location of the initial layer of ions at: (1) h o m = 115   km , where their maximal density in the E region is mostly observed, and (2) h o m =   100 km , where the high density of metallic ions (e.g., Fe+) is expected [33]. These initial ion/electron distributions N i ( z , t = t o ) = N e ( z , t = t o ) N o ( z ) correspond to the Gaussian-type distribution defined at the initial moment t = t o by Equation (17), where the ion scale is chosen to be approximately twice as large ( H i c = 16   km ) as the atmospheric scale (H = 8–10 km). In the presented simulation, the neutral particle densities of the lower thermosphere are used from the NRLMSISE-00 model [42], for the mid-latitude regions 45° ± 2° N, 45° ± 2° E, and I = 61° ± 2°, for a time around the solar activity minimum phase (e.g., 12 April 1998). The geomagnetic field and its inclination for the given region of the mid-latitude lower thermosphere are chosen in accordance with the International Geomagnetic Reference Field (IGRF) model (https://geomag.bgs.ac.uk/data_service/models_compass/igrf_calc.html (accessed on 3 April 2023)). In the numerical simulation, the height and time resolutions are 0.2 km and 1.4 10 4   h , respectively.
The height profile of the ion vertical drift velocity w i ( z , t ) , its vertical changes w i ( z , t ) / z and, correspondingly, the ion/electron density N i ( z , t ) = N e ( z , t ) , as in Equation (11), will be solved numerically similarly to [12,43,44,45] for the whole night t t o 12   h .
Figure 2 shows that the initial layer of ion/electron density N i ( z , t = t o ) = N e ( z , t = t o ) N o ( z ) (panel a) undergoes a slight decrease caused by the ions’ ambipolar diffusion. The downward neutral wind velocity V z = 4   m / s causes ions downward flux with velocity w i v ( z , t ) = C z ( z ) V z accompanied by a maximal decrease in density at height h ≈ 121 km (panel b) by divergence rate C v min ( V x = 0 ,   V = 0 ,   V z = 4   m / s ) = 4.9 × 10 5   s 1 . The upward neutral wind with V z = 4   m / s (panel c) causes ion/electron upward flux, accompanied by their increased diffusive displacement for upper heights and convergence with the maximal rate C v max ( V x = 0 ,   V = 0 ,   V z = 4 m / s ) = 4.9 × 10 5   s 1 at h ≈ 121 km resulting in density increase (up to about 1.09) for time t t o 0.9   h .
Similar behavior of the ion/electron density occurs for its initial distribution with peak height at h o m = 100   km (Figure 2d) and undergoes negligibly small changes due to ambipolar diffusion. In this case, downward neutral wind V z = 4   m / s causes the most downward flux of the initial layer by ion drift velocity w i v ( z , t ) = C z ( z ) V z , which is about the same as neutrals velocity ( w i v V z ) for lower heights about h < 110 km ( C z 1 , Figure 1a). In this case, at ions divergence rate peak height h = 121 km for a time interval t t o 3   h , a relatively significant decrease ( N e ( h = 121   km , t t o = 3   h ) / N o ( h = 121   km ,   t = t o ) 0.003 ) occurs compared to the upper location of the initial layer ( N e ( h = 121   km ,   t t o = 3   h ) / N o ( h = 121   km ,   t = t o ) 0.02 ). For this lower location of the initial layer, the upward neutral wind V z = 4   m / s causes upward drift of this layer, which converges ( N e ( h = 121   km , t t o ) / N o ( h = 121   km ,   t = t o ) > 1 ) at about h = 121 km, with a slightly smaller density ( N e ( h = 126   km , t t o = 0.9   h ) / N o m ( h o m = 100   km ,   t t o = 1.9 h ) 1.07 ) compared to the one at the upper location of the initial peak ( N e ( h 126   km , t t o = 0.9   h ) / N o m ( h o m = 115   km ,   t t o = 1.9     h ) 1.09 ), and peak value at about height region h = 126 km occurs a little later t t o = 1.9   h than in the case when the initial peak is located at upper height.
So, the vertical wind can influence the dynamics of ion/electron convergence processes in the lower thermosphere. In the considered cases of constant wind velocity V z (Figure 2), the maximal/minimal value of convergence/divergence rate is determined by the peak value of C z z , as in Equation (19), and located at the same height at about h = 121 km (with ν i n = ω i ). This location is the same (Figure 1) as one of the maximal values of vertical changes of the meridional wind velocity factor C x z .
The ion/electron convergence/divergence rate due to the vertical neutral wind always occurs in the lower thermosphere, and its influence on their density, in combination with horizontal wind, is expected to be significant. In this case the atmospheric waves, i.e., AGWs, can also create regions (spaced with distances about vertical wavelength) with maximal ( C v z max > 0 )/minimal ( C v z min < 0 ) values C v z , as in Equation (16); thus, the development of additional convergence/divergence processes and formation of multilayered sporadic E is expectable, similarly to the cases given in [26].
The presence of neutral wind with meridional and zonal components will have a combined effect on ion/electron convergence/divergence processes. In this case, the existence of the region with minimized downstream and upstream flows (e.g., w i v = 0 ), as in Equation (7), is possible and additional convergence (or divergence) of ions at height (h = 121 km) C z max = cos 2 I 2 H can be more important for ions convergence close to this height region, compared to the case of absence of meridional and zonal components, i.e., compared to the case of only vertical wind (Figure 2).
Thus, the vertical wind can have a significant effect on the vertical movement of the ion/electron layer, as well as on their convergence/depletion process. Note that for constant vertical wind ( V z = ± 4   m / s ), the height profile of the ion drift velocity w i v ( z , t ) = C z ( z ) V z is proportional to C z (Figure 1a), as in Equation (7). Therefore, the height profile of its vertical change w i v z = V z C z z , determining ion convergence rate, as in Equation (14), is proportional to C z z (Figure 1b). For the sake of brevity, these profiles are not demonstrated. In this case, vertical changes of the ion-neutral collision frequency ν i n ( z ) , which result in maximal values in convergence ( V z > 0 )/divergence ( V z < 0 ) rate factor C z z (Figure 1b), as in Equation (19), and, correspondingly, the ions convergence into Es-type layer or their depletion at the region with ν i n = ω i , are possible even in case of constant vertical wind (Figure 2). Similarly, maximal positive (or minimal negative) values in the convergence (or divergence) factors C x z and C y z (Figure 1b) show the possibility of the formation of the Es layer (or ions depletion) due to constant meridional and zonal wind [12].
We consider the only horizontal wind whose meridional V x ( z , t ) and zonal V y ( z , t ) velocities are determined by the HWM14 data [35], as well as the cases when the vertical component ( V z = ± 4   m / s ) is incorporated. In these cases, we will estimate the profiles of the ion vertical drift velocity and its vertical change. Additionally, we will identify the regions of maximum convergence ( C t max > 0 ) and divergence ( C t min < 0 ) rates of ions/electrons, as well as regions (with w i v = 0 ) where the formation of the Es-type layer is possible. In the same cases, the behavior of ions/electrons during the night will be presented and processes of the expected Es-type layer, along with processes of ion/electron density depletion will be shown.
Figure 3 demonstrates nightly variations of the typical height profiles of (a) the meridional wind velocity V x ( z , t ) and (b) the zonal wind velocity V y ( z , t ) , described by HWM14 data (Day 102, t o = LT   23 : 00 ). To determine the height profile of the ion vertical drift velocity w i v , as in Equations (7)–(10), and its vertical changes w i v z , as in Equations (13)–(16), we will use the above-mentioned meridional V x ( z , t ) and zonal V y ( z , t ) wind velocity data (Figure 2), as well as the data where the vertical wind with V z = 4   m / s and V z = 4   m / s is incorporated.
Figure 4a–c show that in the lower thermosphere, in case of the existence of neutral wind, the ion vertical drift velocity w i v height profile is characterized by regions with w i v = 0 or w i v 0 . In these cases, vertical changes of the ion vertical drift velocity (Figure 4d–f) also have minimal negative values w i v z min < 0 ( C t max > 0 ). According to Equations (17) and (18), in the regions where the above-mentioned conditions are satisfied the Es-type layers can be formed. These layers are localized at regions with w i v = 0 or descent to the regions where w i v 0 .
In the case of the presence of wind with a downward vertical velocity component V z = 4   m / s (Figure 4b,e) the regions w i v = 0 are shifted downward compared to the case of only horizontal wind (Figure 4a,d). For example, the height regions with w i v = 0 (Figure 4a) and the highest minimal negative value in w i v z min < 0 ( C t max > 0 ) (about −0.003 s−1), located at about 128–129 km (Figure 4d), for a time interval 1.5   h t t o 4.5   h (Figure 4a), are shifted downward to about 126–127 km (Figure 4b).
In the case of the upward component of the neutral wind velocity V z = 4   m / s (Figure 4c,f) the regions of upper nodes are displaced upward. For example, the height regions w i v = 0 located at about 128–129 km, for time intervals 1.5   h t t o 4.5   h , are displaced upward to about 129–130 km (Figure 2c).
At lower heights h < 110, in the case of only horizontal wind ( V z = 0 ), the ion drift velocity vanishes w i v 0 (Figure 4a), as in Equations (7)–(10). When the vertical wind component is directed upward ( V z > 0 ) or downward ( V z < 0 ), then for these regions of the lower thermosphere (about h < 110 km) the ion vertical drift velocity is directed downward with w i v 4   m / s or upward with w i v 4   m / s , respectively.
Figure 4 also shows the presence of the regions with w i v z max > 0 ( C t min < 0 ), where ion/electron depletion is expectable. The altitude of these regions varies in the height interval from about 116 km ( t = t o ) to 112 km ( t t o = 12   h ). The noticeable minimal values C t min 0.0021       s 1 occur for time t t o 1.7   h . The presence of the constant downward V z = 4   m / s (panel e) and upward V z = 4   m / s (Figure 4f) vertical components of neutral wind do not give significant changes in these parameters compared to those in the case of only horizontal wind (Figure 4d).
We will demonstrate the behavior of the ion/electron density height profile N e ( z , t ) / N o m and formation of the Es-type high-density thin layer (with peak density N e m / N o m > 1 ), H i c ( t = t o ) > H i c ( t > t o ) and/or ion/electron density depletion ( N e / N o m < < 1 ). These phenomena can be predicted by analyses of ions drift velocity w i ( z , t ) (Figure 4a–c) and its vertical changes w i v z (Figure 4d–f), which are estimated using HWM14 data of meridional V x ( z , t ) and zonal wind velocities V y ( z , t ) (Figure 3a,b).
Neutral wind leads to the existence of the regions where the condition w i v z min < 0 i.e., C t max > 0 is satisfied (Figure 3 and Figure 4). Figure 5 shows that close to the above-mentioned regions, the formation of the Es layers and their localization occur at regions with w i = 0 (all panels) and w i v 0 (Figure 5a,d). For the used horizontal wind profile (Figure 3) the high-density Es layer N m 1 / N o m = 5.4 (with H i c = 1.8   km ) is formed for a time around t t o 3 5   h (Figure 5a,d) at the location of the ion drift velocity node h m 1 = 128   km ( w i = 0 ). Additionally, the secondary Es layer N m 2 / N o m = 1.02 is formed at a time t t o = 12   h at lower heights ( h m 2 < 100   k m ) and it descends towards the height about h m 2 = 93.4   km , where w i v 0 . When the density peak height of the ion/electron initial distribution is relatively higher h o m = 115   km (close to the height region of the upper maximum of convergence rate C t max > 0 and node w i = 0 ), the density in the upper Es layer (with H i c = 1.8   km ) is higher N m 1 / N o m = 5.4 (Figure 5a) than when the initial distribution density peak height is lower h o m = 100   k m , ( N m 1 / N o m = 2.8 , H i c = 1.8   km ) (Figure 5d). In these cases (Figure 5a,d), the lower location (close to the lower maximum of convergence rate C t max > 0 around h = 99 km–104 km) of initial peak height (panel d, h o m = 100   km ) causes higher density N m 2 / N o m = 2.7 in the lower Es layer ( h m 2 = 93.4   km ) (Figure 5d) than in case of higher location h o m = 115   km of initial peak height (Figure 5a).
Figure 5 also shows the existence of a divergence region C t min < 0 , where the ion/electron density continuously decreases/depletes ( N e ( z , t ) / N o m < 1 ). The location of this region varies in the height interval from 116 km ( t = t o ) to 112 km ( t t o =   12   h ). For example, at the time t t o = 3   h , for the higher location of the initial peak height N e ( h = 116   km ,   h o m = 115   km ,   t t o = 3   h ) / N o ( h = 116   km ,     t = t o ) 0.02 (Figure 5a), while for the lower location N e ( h = 100   km ,   h o m = 115   km ,   t t o = 3   h ) / N o ( h = 116   km ,     t = t o ) 0.09 (Figure 5d).
Below the height regions of ions depletion (e.g., h < 110 km) where the ions drift velocity w i v 0 , at around h = 105 km the ion/electron total convergence rate is maximal C t max ( 2 ± 9 ) × 10 5   s 1 and their ions converge into descending ( w i 0 ) Es-type layer ( N e ( z , t ) / N o m > 1 ), with the continuous increase in densities and with the peak height located at about h m 2 = 93.4   km at t t o = 12   h . Here, the peak density of the formed Es N m 2 / N o m = 1.02 (Figure 5a) is smaller for the higher location of the ion/electron initial layer ( h m o = 115   km ) compared to those for the lower location (close to the region with C t max > 0 ) of the initial layer ( h m o = 100   km ) when N m 2 / N o m = 2.7 (Figure 5d).
So, in the case of the presence of horizontal wind with only meridional and zonal components of velocities (Figure 3), two regions occur with condition C t max > 0 (Figure 4). This leads to the formation of sporadic E with a double peak at about h m 1 = 128   km and h m 2 = 93.4 km (Figure 5a,d). In this case, the upper Es layer is located at the height region with w i v = 0 ( h m 1 = 128   km ) while the lower one descends to the height region where ions’ vertical drift velocity vanishes w i v 0 .
Figure 5b,e show that the downward direction of the vertical component of wind velocity V z < 0 , which causes an increase in ion/electron downward flux, in case of V z = 4   m / s , leads to the lower location h m 1 ( V z = 4   m / s ) = 126.4   km < h m 1 ( V z = 0 ) = 128   km of the ion vertical drift velocity node ( w i = 0 ) and, therefore, the Es layer formation occurs at about h m 1 = 126.8   km . Here, the higher location of the initial density peak height ( h m o = 115   km ) and increase in ion/electron flux from upper regions to the region of ion drift velocity node h ≈ 126.4 km lead to the higher density N m 1 / N o m = 5.6 ( H i c = 1.7 km ) (Figure 5b) in the formed Es layer for a time around t t o = 3.4   h . Here, the relatively small value in the ion’s diffusive displacement to the heights above its convergence region also causes an increase in their peak density. In the case of the lower location, below the divergence region (about height 116 km) of the initial layer ( h m o = 100 km ), the increased downward flux of ions, caused by the downward vertical wind V z = 4   m / s , leads to the formation of sporadic E with relatively lower density, N m 1 / N o m = 1.3 ( H i c = 1.7 km ) and at the lower location 126.4 km (Figure 5e), than in the case of higher location of initial density peak (Figure 5b).
Figure 5b,e show that the downward wind with additional depletion around the h ≈ 121 km region (Figure 2b,e) causes a relatively higher decrease in ion electron density, compared to the case of the wind with only horizontal components of velocity.
In the case of downward vertical wind velocity, secondary high-density Es layers are not formed at the lower heights of h < 110 km. This is because, unlike the case of wind with only horizontal velocity components, the region w i v 0 does not occur and the ions/electrons layer descends with a velocity close to neutral vertical wind velocity w i v 4   m / s ( w i v V z ) without significant changes in its density (Figure 5b,e). During the downward drift with velocity w i v 4   m / s , the ion/electrons cross the region of their convergence rate maximum (around h = 99–104 km) within a short time (about 1 h). For the values, C t max 2 10 5 s 1 this does not lead to a significant increase in the ions/electrons density, neither in the case of higher ( h m o = 115 km ) nor lower ( h m o = 100 km ) location of the initial layer. As a result, the formation of the secondary Es-type layer ( N e ( z , t ) / N o m > 1 ) at lower heights <110 km does not occur. Figure 5c,f show that the wind velocity component directed upward V z > 0 causes the additional upward drift of ions/electrons, i.e., reduction of their downward drift velocity, as in Equation (7). In this case, the value V z = 4   m / s leads to the relatively higher location of the ions drift velocity node w i ( h 129.4 km ) = 0 (compared to h m 1 = 128   km in the case of V z = 0 ) and the Es layer formation occurs at about h m 1 = 129.4   km . In this case of the higher location of the ion/electron initial layer ( h m o = 115 km ), the additional ion/electron upward drift and increase in ambipolar diffusion also lead to a lower density N m 1 / N o m = 1.9 (Figure 5c) in the formed Es layer than it was in cases of V z = 0 (Figure 5a) and V z = 4   m / s (Figure 5b).
The lower location of the initial peak height ( h m o = 100 km ) and upward vertical wind ( V z = 4   m / s ) support the increase in ion/electron flux from lower to upper heights and their additional convergence at the height region of about 121 km (Figure 2f) leads to the relatively high density N m 1 / N o m = 4.7 in this Es layer formed at about 129.4 km height compared to the densities in Es layers formed for cases of V z = 0 ( N m 1 / N o m = 2.7 , h m 1 = 128 km Figure 5d) and V z = 4   m / s ( N m 1 / N o m = 1.3 , h m 1 = 126.8 km Figure 5e). In this case, in spite of the relatively lower value of the convergence rate, due to upward vertical wind ( 4.9 × 10 5 s 1 ), the combined effect of horizontal wind direction, its shear and vertical wind gives high density ( N m 1 ( V z = 4   m / s ) / N 1 m ( V z = 0 ) 1.7 ) in the Es-type layer at about 129.4 km height.
The upward neutral wind and corresponding increased upward ion/electrons flux also result in their depletion in the divergence region (around 116 km). In both cases of the ion/electron initial layer being located at higher or lower altitudes, their depletion ( N e ( z = 116   km ,     t t o > 0 ) / N e ( z = 116   km , t = t o ) < 1 ) in the divergence region with C t min < 0 for time of about t t o > 2.5   h is bigger (Figure 5b,c,e,f), compared to the one during the absence of vertical wind (Figure 5a,d). However, for a time of about t t o < 2.5   h the relatively small depletion at the height region around 116 km, due to ion/electron upward flux to this region by vertical wind velocity V z > 0 , occurs at the lower location ( z o m = 100   km ) of the initial layer (Figure 5f).
Figure 5c,f also show that the increased ions/electrons upward flux at a lower height (h < 110 km) due to upward vertical wind V z > 0 also results in their density rarefication at lower height regions with maximal values of convergence rate (99–104 km) and the formation of the secondary Es-type layer at these height regions do not occur.
The upward/downward wind causes the vanishing of the region w i v 0 at lower heights h < 110 km ( w i v V z ) and the formation of an Es-type layer in this region is less expectable. So, in this case of vertical wind velocity, the necessary condition C t max > 0 for Es layer formation at regions about 99–104 km occurs, but it is not sufficient for its formation. It is important to note that the presence of condition C t max > 0 for the formation of the Es-type is necessary but not sufficient. In the considered cases (Figure 5), the condition C t max > 0 occurs for the whole night, but the development of the processes of ions/electrons convergence into a thin dense Es-type layer also depends on the initial distribution of these charged particles. For instance, if we use the same distribution for ion/electron density (Figure 2a,d) for different moments of the night, then the convergence of ions/electrons into the Es-type layer and its localization will occur at the regions where conditions C t max > 0 , w i v = 0 and w i v 0 are satisfied at a particular moment of night mentioned above. This is not demonstrated for the sake of brevity.
Note that we have demonstrated theoretically the formation of the Es layer for one typical wind profile, including the HWM14 data, for the meridional and zonal wind velocities at mid-latitudes (45° N; 45° S). This methodology is not limited to HWM14 and can be applied to other model data, as well as measured wind profiles, including the data close to real time.
We can also apply our model to different initial distributions of charged particles, which fits both theory and measurements. In this case, by theoretically estimating the height profile of ions’ vertical drift velocity w i v , as in Equation (7), and w i v z , the regions with conditions w i v z min < 0 ( C t max > 0 ) necessary for the development of ion/electron convergence into narrow dense Es-type layers also can be identified. This layer is located at the altitudes with minimal value of ions upstream and downstream flows (e.g., w i v = 0 ) or it moves to the region with w i v 0 and C t max > 0 . In this case, the condition C t max = ( C v + C s h + C v z ) max > 2 D a H i c 2 necessary for sporadic E formation, for a particular mid-latitude region (with given B and I), is determined not only by the horizontal ‘windshear’ convergence rate C s h , as in Equation (15), but also by the combined effect of rates caused by the realistic wind (three-dimensional) velocity magnitude, direction C v , as in Equation (14), and the vertical changes of vertical wind C v z , as in Equation (16). We note that in the lower thermosphere, the value of ions ambipolar diffusion, as in Equation (10), does not exceed 0.00001 m2/s. As a result, the condition C t max > 2 D a H i c 2 at the initial time ( H i c = 10 ÷ 16 km ) is satisfied for the small values (e.g., ± 4   m / s ) of wind velocity. However, the effect of diffusion control on the Es layer is important during its formation (e.g., when H i c = 1 ÷ 3 km ).
Below, we demonstrate the values of C t max > 0 and C t min < 0 (Figure 4 and Figure 5) and corresponding values of parameters C v ( z , t ) and C s h ( z , t ) . These values are determined by w i v z profile, which in turn is estimated using the HWM14 data (Figure 3).
Figure 6 shows that the condition C t max > 0 ( w i v z min < 0 ) necessary for ion/electron convergence into Es-type thin layer (panels a and c) and the condition C t min < 0 ( w i v z max > 0 ) necessary for their divergence/depletion (panel b) are determined by the combined effect of convergence/divergence rates C v ( z , t ) , determined by wind velocity direction and magnitude, and C s h ( z , t ) , determined by horizontal wind velocity vertical shear. For instance, in the case of only horizontal wind, the convergence rate in the upper Es layer, formed for time tto = 3–5 h at heights about 126–130 km (Figure 5), C t max 0.0027   s 1 (Figure 4 and Figure 5) is the sum of C v 0.00099   s 1 and C s h 0.0017   s 1 . The latter shows that in the upper Es layer (Figure 5) about 60% of ion/electron convergence into this layer is attributed to the ‘shear effect’ (by the presence of horizontal wind with vertical shear).
At lower heights (h < 110 km), the higher value C t max 9 10 5   s 1 in the Es layer (Figure 5a,d) is primarily caused by the convergence rate C v , which is determined by wind value and direction. In this region, the shear effect is negligible and even negative C s h < 0 .
At the altitudes of approximately = 116 km, in the time interval tto < 4 h, the magnitude of C t min < 0   w i v z max > 0 ) reaches 0.002   1   s 1 , with approximately 30% ( C s h 0.0006   s 1 ) of this value caused by the horizontal wind shear effect, while the remaining 70% ( C v 0.0015   s 1 ) is determined by the wind magnitude and direction.
At the heights of the peaks of the upper Es layer (h = 126 km–130 km), the divergence (≈116 km) and the lower Es layer (93–104 km), the component of the wind velocity has negligible contribution in the values of C t max > 0 and C t min < 0 . However, the ions, due to their vertical drift caused by the vertical component of the wind velocity, pass through the region where only upward/downward vertical wind ( V z = ± 4   m / s ) causes the maximum convergence/divergence effect (height h ≈ 121 km); therefore, the vertical wind component influences the densities of ions in Es layers, as well as in the depletion regions.
The presented analytical and numerical study shows that the realistic (three-dimensional) profile of the neutral wind velocity in the lower thermosphere (h = 90 km–150 km) significantly determines the behavior of ion/electrons density in this region and therefore, the region of sporadic E (Es) formation and localization. Here, vertical profiles of the ions and vertical drift velocity and its vertical change, estimated by taking into account the vertical component of the wind velocity along with its horizontal component, influence the Es layer formation processes, since the above-mentioned profiles determine the condition, w i z min < 0 (or C t max > 0 ), necessary for ion/electron convergence into a narrow dense layer. Additionally, by the condition w i z max > 0 (or C t min < 0 ) the location of ion/electrons depletion regions is identified. It is important to take into account the presence of vertical wind because it affects both the location of regions of ions/electrons convergence into a narrow dense layer, as well as the dynamics of the development of these processes.
The presented theory can be extended by taking into account the production and loss rates (Q and L) of charged particles (e.g., NO+ and O2+) in Equation (1), as well as considering the effects of the electric field and the ion drift caused by this field, in Equation (6). This extended study is planned for the near future.
The wind velocity profile given by the HWM14 is mainly determined by tidal motions and planetary waves [35]. The presented theoretical model allows its application for relatively short-term wind variations caused by AGWs and/or other dynamical processes. As a result, it is possible to study the regional features of the origin of the discussed ionospheric irregular structures caused by the above-mentioned processes. To study the feature of ions’ global and regional distributions, it is important to consider molecular ions (e.g., NO+ and O2+), along with metal ions (e.g., Fe+, Be+) in the three-dimensional wind. This is also the subject of future studies.

4. Conclusions

(1) The behavior of the ion/electron concentration and the possibility of the formation of an Es-type layer, as well as their depletion processes under the influence of three-dimensional, which takes into account its vertical V z component along with the meridional V x and zonal V y components of the velocity V , are shown. In this case, the presence of a minimum negative value of the vertical change of the vertical drift velocity of ions w i v z min < 0 in the lower thermosphere, as the main condition for predicting the formation of the Es layer [26], was extended analytically and with appropriate numerical methods, in case of presence of the vertical component of the wind velocity V z .
(2) Using the analytical approach for ion/electron concentration with consideration of ambipolar diffusion (17) shows that their narrow-layer total convergence rate C t = C v + C s h + C v z ( N e t w i v z = C t > 0 ) is determined by both horizontal and vertical wind velocity magnitude, direction ( C v ) and vertical change ( C s h and C v z ); see Equations (14)–(16). The vertical movement of the ion/electron layer and, accordingly, the localization of the Es-type layer in the regions of the minimal reduction (or with w i v = 0 ) of their lower and upper fluxes or regions where these fluxes vanish ( w i v 0 ), is determined by the profile of the ions vertical drift velocity w i v , as in Equation (7). In a given region (with particular B and I), it is possible to identify the regions of the development of ion convergence C t max > 2 D a H i c 2 (Es-type layer formation) and their depletion conditions C t min < 0 . The convergence ( C v z > 0 )/divergence ( C v z < 0 ) factor caused by the vertical wind is also influenced by the neutral particles’ wind velocity profile.
(3) Downward/upward vertical wind will shift the convergence regions (with w i v = 0 or w i v 0 ) down/up, respectively, from their locations that existed in the case of just meridional and zonal winds. In the case of the presence of the wind velocity vertical component, the rates of convergence and divergence of the ions are different from their values in the case of just a horizontal wind.
Along with the profiles of the meridional ( C x , C x / z ) and zonal ( C y , C y / z ) factors of the wind velocity, the ones of the vertical wind factor ( C z ) and its vertical change ( C z / z ) are important in the height profile of the ion drift velocity w i v and its vertical change w i v z (see Figure 1). In this case, it is shown for the first time that the upward constant wind ( V z > 0 ), like the northward wind ( V x > 0 ), causes the convergence of ions with the maximum rate in the region of the maximum of C z / z (h ≈ 121 km), where the collision frequency of ions-neutrals is equal to their cyclotron frequency ν i n = ω i (Figure 1). In this height region, a downward wind ( V z < 0 ), like a southward one ( V x < 0 ), causes ion/electron depletion (Figure 2).
(4) The significant influence of horizontal wind, ions diffusion and their initial distribution, as well as the vertical wind on the behavior of ions/electrons and, therefore, the formation of Es-layer (as well as ion density depletion) has been demonstrated by numerical simulations (Figure 2 and Figure 5).
Accordingly, at first, the behavior of ions/electrons in the presence of only vertical wind was demonstrated (Figure 2), then the case of only horizontal wind determined by typical HWM14 data was considered (Figure 3) and, finally, the case of three-dimensional wind with the horizontal and vertical components was considered.
The vertically upward/downward component of the three-dimensional wind velocity will shift the regions of ion convergence (e.g., w i v = 0 or w i v 0 ) upward or downward, respectively, and change both the location of the formed Es layer and its density relative to values of corresponding parameters estimated in the case of only horizontal wind (Figure 5). Here, the increase in diffusive displacement towards the upper regions (e.g., h > 120 km) also reduces the density of the upper Es-layer and increases it at the lower location.
The location of the initial ion/electron layer relative to the expected regions of convergence, as well as the ions drift caused by the vertical wind, influences the conditions for their convergence into a narrow layer in these regions (Figure 4 and Figure 5). For example, in the lower regions (h < 110 km), due to the downward/upward drift of ions ( w i v V z ), caused by downward/upward wind, respectively, the ions localization does not occur in the regions with w i v 0 (Figure 5b,c,e,f), which appeared in cases of only horizontal wind (Figure 5a,d). So, the presence of a maximum rate of convergence ( C t max > 0 ) in the lower regions may not be sufficient for the formation of a high-density Es layer during the vertical wind.
In the case of three-dimensional wind with an upward component (Figure 5f), the formation of an Es layer with a density higher than in the case of only horizontal wind was shown (Figure 5d). In this case, the additional convergence caused by the vertically upward wind in the region increases the convergence caused by the horizontal wind in the regions close to the ion vertical drift velocity node.
(5) In the region of maximum ion/electron divergence w i z max > 0 (or C t min < 0 ) (e.g., at a height of about 116 km), during the presence of vertical upward/downward wind the depletion of ion/electron concentration is higher than in case of the horizontal wind alone (Figure 5).
(6) The convergence/divergence factors caused by the values of both the horizontal wind vertical shear ( C s h ) and magnitude, and the direction of the wind velocity ( C v z ), which determine the maximum rates of convergence ( C t max > 0 ) and divergence ( C t min < 0 ), have been shown for the whole considered nighttime interval (Figure 6).
(7) The vertical wind velocity used in the numerical simulations is much smaller in the lower thermosphere than its observed [41] and model horizontal wind averages [35]. This can be caused by the tidal motions, as well as the propagation of AGWs. In the presented theoretical model, the three-dimensional profile of the wind velocity, which takes into account the vertical component along with its zonal and meridional ones, is important for the prediction of the regions of the ion convergence and formation of Es layers, and also the regions of ion depletion. Therefore, consideration of the significant effect of the wind shear [3,38,39], as well as the wind velocity magnitude and direction [12,26], including the vertical wind described above, is important to study the regional and global distributions of the metal ions [34,46,47,48,49,50] that make up the lower thermosphere.

Author Contributions

G.G.D.: conceptualization, investigation, methodology, and writing—original draft; G.D.: formal analyses, funding acquisition, methodology, and validation; M.T.: investigation, validation, visualization, writing—review and editing; L.T.: visualization and supplementary simulations. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by Georgian Shota Rustaveli National Science Foundation Grant no. FR-21-2282 5.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Acknowledgments

We used the NRLMSISE-00 model and data are available at https://ccmc.gsfc.nasa.gov/modelweb/models/nrlmsise00.php, accessed on 3 April 2023. We used a public domain database, The Open Science Framework, for dissemination (https://osf.io/search/?q=Dalakishvili&page=1 (accessed on 3 April 2023)) of our numerical code.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Whitehead, J.D. Recent work on mid-latitude and equatorial sporadic-E. J. Atmos. Terr. Phys. 1989, 51, 401–424. [Google Scholar] [CrossRef]
  2. Mathews, J.D. Sporadic E: Current views and recent progress. J. Atmos. Sol.-Terr. Phys. 1998, 60, 413–435. [Google Scholar] [CrossRef]
  3. Haldoupis, C. Midlatitude sporadic E. A typical paradigm of atmosphere-ionosphere coupling. Space Sci. Rev. 2012, 168, 441–461. [Google Scholar] [CrossRef]
  4. Maruyama, T.; Kato, H.; Nakamura, M. Meteor-induced transient sporadic E as inferred from rapid-runionosonde observations at midlatitudes. J. Geophys. Res. Space Phys. 2008, 113, A09308. [Google Scholar] [CrossRef] [Green Version]
  5. Kelley, C.M. The Earth’s Ionosphere; Elsevier: Amsterdam, The Netherlands, 2009; Section 6; pp. 267–342. [Google Scholar]
  6. Garcia-Fernandez, M.; Tsuda, T. A global distribution of sporadic E events revealed by means of CHAMP-GPS occultations. Earth Planets Space 2006, 58, 33–36. [Google Scholar] [CrossRef] [Green Version]
  7. Yu, B.; Xue, X.; Yue, X.; Yang, C.; Yu, C.; Dou, X.; Ning, B.; Hu, L. The global climatology of the intensity of the ionospheric sporadic E layer. Atmos. Chem. Phys. 2019, 19, 4139–4151. [Google Scholar] [CrossRef] [Green Version]
  8. Whitehead, J.D. The formation of the sporadic E layer in the temperate zone. J. Atmos. Terr. Phys. 1961, 51, 20–49. [Google Scholar] [CrossRef]
  9. Axford, W.I. The formation and vertical movement of dense ionized layers in the ionosphere. J. Geophys. Res. 1963, 68, 769–779. [Google Scholar] [CrossRef]
  10. Cosgrove, R.B.; Tsunoda, R.T. A direction-dependent instability of sporadic-E layers in the nighttime midlatitude ionosphere. Geophys. Res. Lett. 2002, 29, 1864. [Google Scholar] [CrossRef]
  11. Bishop, R.L.; Earle, G.D.; Larsen, M.F.; Swenson, C.M.; Carlson, C.G.; Roddy, P.A.; Fish, C.; Bullett, T.W. Sequential observations of the local neutral wind field structure associated with E region plasma layers. J. Geophys. Res. 2005, 110, A04309. [Google Scholar] [CrossRef]
  12. Dalakishvili, G.; Didebulidze, G.G.; Todua, M. Formation of sporadic E (Es) layer by homogeneous and inhomogeneous horizontal winds. J. Atmos. Sol.-Terr. Phys. 2020, 209, 105403. [Google Scholar] [CrossRef]
  13. Collinson, G.A.; McFadden, J.; Grebowsky, J.; Mitchell, D.; Lillis, R.; Withers, P.; Vogt, M.F.; Benna, M.; Espley, J.; Jakosky, B. Constantly forming sporadic E-like layers and rifts in the Martian ionosphere and their implications for Earth. Nat. Astron. 2020, 4, 486–491. [Google Scholar] [CrossRef]
  14. Shinagawa, H.; Tao, C.; Jin, H.; Miyoshi, Y.; Fujiwara, H. Numerical prediction of sporadic E layer occurrence using GAIA. Earth Planets Space 2021, 73, 28. [Google Scholar] [CrossRef]
  15. Haldoupis, C.; Pancheva, D. Planetary waves and midlatitude sporadic E layers: Strong experimental evidence for a close relationship. J. Geophys. Res. 2002, 107, 1078. [Google Scholar] [CrossRef]
  16. Yeh, W.-H.; Liu, J.-Y.; Huang, C.-Y.; Chen, S.-P. Explanation of the sporadic-E layer formation by comparing FORMOSAT-3/COSMIC data with meteor and wind shear information. J. Geophys. Res. Atmos. 2014, 119, 4568–4579. [Google Scholar] [CrossRef]
  17. Didebulidze, G.G.; Lomidze, L. The formation of sporadic E layers by a vortical perturbation excited in the horizontal wind shear flow. Ann. Geophys. 2008, 26, 1741–1749. [Google Scholar] [CrossRef] [Green Version]
  18. Resende, L.C.A.; Batista, I.S.; Denardini, C.M.; Batista, P.P.; Carrasco, A.J.; Andrioli, V.F.; Moro, J. Simulations of blanketing sporadic E-layer over the Brazilian sector driven by tidal winds. J. Atmos. Sol.-Terr. Phys. 2017, 154, 104–114. [Google Scholar] [CrossRef]
  19. Liu, Y.; Zhou, C.; Tang, Q.; Li, Z.; Song, Y.; Qing, H.; Ni, B.; Zhao, Z. The seasonal distribution of sporadic E layers observed from radio occultation measurements and its relation with wind shear measured by TIMED/TIDI. Adv. Space Res. 2018, 62, 426–439. [Google Scholar] [CrossRef]
  20. Oikonomou, C.; Haralambous, H.; Leontiou, T.; Tsagouri, I.; Buresova, D.; Mosna, Z. Intermediate descending layer and sporadic E tidelike variability observed over three mid-latitude ionospheric stations. Adv. Space Res. 2022, 69, 96–110. [Google Scholar] [CrossRef]
  21. Abdu, M.A.; MacDougall, J.W.; Batista, I.S.; Sobral, J.H.A.; Jayachandran, P.T. Equatorial evening prereversal electric field enhancement and sporadic E layer disruption: A manifestation of E and F region coupling. J. Geophys. Res. 2003, 108, 1254. [Google Scholar] [CrossRef] [Green Version]
  22. Yokoyama, T.; Yamamoto, M.; Fukao, S. Computer simulation of polarization electric fields as a source of midlatitude field-aligned irregularities. J. Geophys. Res. 2003, 108, 1054. [Google Scholar] [CrossRef]
  23. Nygrén, T.; Jalonen, L.; Oksman, J.; Turunen, T. The role of electric field and neutral wind direction in the formation of sporadic E-layers. J. Atmos. Terr. Phys. 1984, 46, 373–381. [Google Scholar] [CrossRef]
  24. Wakabayashi, M.; Ono, T. Multi-layer structure of mid-latitude sporadic-E observed during the SEEK-2 campaign. Ann. Geophys. 2005, 23, 2347–2355. [Google Scholar] [CrossRef] [Green Version]
  25. Hysell, D.L.; Munk, J.; McCarrick, M. Sporadic E ionization layers observed with radar imaging and ionospheric modification. Geophys. Res. Lett. 2014, 41, 6987–6993. [Google Scholar] [CrossRef]
  26. Didebulidze, G.G.; Dalakishvili, G.; Todua, M. Formation of Multilayered Sporadic E under an Influence of Atmospheric Gravity Waves (AGWs). Atmosphere 2020, 11, 653. [Google Scholar] [CrossRef]
  27. Roddy, P.A.; Earle, G.D.; Swenson, C.M.; Carlson, C.G.; Bullett, T.W. The composition and horizontal homogeneity of E region plasma layers. J. Geophys. Res. 2007, 112, A06312. [Google Scholar] [CrossRef]
  28. Woodcock, K.R.S.; Vondrak, T.; Meech, S.R.; Plane, J.M.C. A kinetic study of the reactions FeO+ + O, Fe+∙N2 + O, Fe+∙O2 + O and FeO+ + CO: Implications for sporadic E layers in the upper atmosphere. Phys. Chem. Chem. Phys. 2006, 8, 1812–1821. [Google Scholar] [CrossRef]
  29. Didebulidze, G.G.; Lomidze, L.N. Double atmospheric gravity wave frequency oscillations of sporadic E formed in a horizontal shear flow. Phys. Lett. A 2010, 374, 952–959. [Google Scholar] [CrossRef]
  30. Didebulidze, G.G.; Dalakishvili, G.; Lomidze, L.; Matiashvili, G. Formation of sporadic-E (Es) layers under the influence of AGWs evolving in a horizontal shear flow. J. Atmos. Sol.-Terr. Phys. 2015, 136, 163–173. [Google Scholar] [CrossRef]
  31. Lin, Y.C.; Chu, Y.H. Model simulations of ion and electron density profiles in ionospheric E and F regions. J. Geophys. Res. Space Phys. 2017, 122, 2505–2529. [Google Scholar] [CrossRef]
  32. Tikhonov, A.N.; Samarsky, A.A. Equations of Mathematical Physics, 5th ed.; Nauka: Moscow, Russia, 1977; Section VI; pp. 452–488. [Google Scholar]
  33. Feng, W.; Marsh, D.R.; Chipperfield, M.P.; Janches, D.; Höffner, J.; Yi, F.; Plane, J.M.C. A global atmospheric model of meteoric iron. J. Geophys. Res. Atmos. 2013, 118, 9456–9474. [Google Scholar] [CrossRef]
  34. Huba, J.D.; Krall, J.; Drob, D. Global ionospheric metal ion transport with SAMI3. Geophys. Res. Lett. 2019, 46, 7937–7944. [Google Scholar] [CrossRef]
  35. Drob, D.P.; Emmert, J.T.; Meriwether, J.W.; Makela, J.; Doornbos, E.; Conde, M.; Hernandez, G.; Noto, J.; Zawdie, K.A.; McDonald, S.E.; et al. An update to the Horizontal Wind Model (HWM): The quiet time thermosphere. Earth Space Sci. 2015, 2, 301–319. [Google Scholar] [CrossRef]
  36. Yiğit, E.; Koucká Knízová, P.; Georgieva, K.; Ward, W. A review of vertical coupling in the Atmosphere–Ionosphere system: Effects of waves, sudden stratospheric warmings, space weather, and of solar activity. J. Atmos. Sol.-Terr. Phys. 2016, 141, 1–12. [Google Scholar] [CrossRef]
  37. Sedlak, R.; Zuhr, A.; Schmidt, C.; Wüst, S.; Bittner, M.; Didebulidze, G.G.; Price, C. Intra-annual variations of spectrally resolved gravity wave activity in the upper mesosphere/lower thermosphere (UMLT) region. Atmos. Meas. Tech. 2020, 13, 5117–5128. [Google Scholar] [CrossRef]
  38. Qiu, L.; Zuo, X.; Yu, T.; Sun, Y.-Y.; Qi, Y. Comparison of global morphologies of vertical ion convergence and sporadic E occurrence rate. Adv. Space Res. 2019, 63, 3606–3611. [Google Scholar] [CrossRef]
  39. Qiu, L.; Yu, T.; Yan, X.; Sun, Y.-Y.; Zuo, X.; Yang, N.; Wang, J.; Qi, Y. Altitudinal and latitudinal variations in ionospheric sporadic-E layer obtained from FORMOSAT-3/ COSMIC radio occultation. J. Geophys. Res. Space Phys. 2021, 126, e2021JA029454. [Google Scholar] [CrossRef]
  40. Banks, P.M.; Kockarts, G. Aeronomy; Academic: New York, NY, USA, 1973; Part A; pp. 217–232. [Google Scholar]
  41. Larsen, M.F. Winds and shears in the mesosphere and lower thermosphere: Results from four decades of chemical release wind measurements. J. Geophys. Res. 2002, 107, 1215. [Google Scholar] [CrossRef]
  42. Picone, J.M.; Hedin, A.E.; Drob, D.P.; Aikin, A.C. NRLMSISE00 empirical model of the atmosphere: Statistical comparisons and scientific issues. J. Geophys. Res. 2002, 107, 1468. [Google Scholar] [CrossRef]
  43. Du Fort, E.C.; Frankel, S.P. Stability conditions in the numerical treatment of parabolic differential equations. MTAC 1953, 7, 135–152. [Google Scholar] [CrossRef] [Green Version]
  44. Lanser, D.; Verwer, G.J. Analysis of operator splitting for advection-diffusion-reaction problems from air pollution modelling. J. Com. Appl. Math. 1999, 111, 201–216. [Google Scholar] [CrossRef] [Green Version]
  45. Hundsdorfer, W.; Verwer, G.J. Numerical Solution of Time-Dependent Advection-Diffusion-Reaction Equations; Springer: Heidelberg, Germany, 2003; pp. 325–417. [Google Scholar]
  46. Dou, X.K.; Xue, X.H.; Li, T.; Chen, T.D.; Chen, C.; Qiu, S.C. Possible relations between meteors, enhanced electron density layers, and sporadic sodium layers. J. Geophys. Res. 2012, 115, A06311. [Google Scholar] [CrossRef] [Green Version]
  47. Delgado, R.; Friedman, J.S.; Fentzke, J.T.; Raizada, S.; Tepley, C.A.; Zhou, Q. Sporadic metal atom and ion layers and their connection to chemistry and thermal structure in the mesopause region at Arecibo. J. Atmos. Sol.-Terr. Phys. 2012, 74, 11–23. [Google Scholar] [CrossRef]
  48. Yuan, T.; Wang, J.; Cai, X.; Sojka, J.; Rice, D.; Oberheide, J.; Criddle, N. Investigation of the seasonal and local time variations of the high-altitude sporadic Na layer (Na) formation and the associated midlatitude descending E layer (E) in lower E region. J. Geophys. Res. Space Phys. 2014, 119, 5985–5999. [Google Scholar] [CrossRef]
  49. Cai, X.; Yuan, T.; Eccles, J.V.; Raizada, S. Investigation on the distinct nocturnal secondary sodium layer behavior above 95 km in winter and summer over Logan, UT (41.7° N, 112° W) and Arecibo Observatory, PR (18.3° N, 67° W). J. Geophys. Res. Space Phys. 2019, 124, 9610–9625. [Google Scholar] [CrossRef]
  50. Yu, B.; Xue, X.; Scott, C.J.; Wu, J.; Yue, X.; Feng, W.; Chi, Y.; Marsh, D.R.; Liu, H.; Dou, X.; et al. Interhemispheric transport of metallic ions within ionospheric sporadic E layers by the lower thermospheric meridional circulation. Atmos. Chem. Phys. 2021, 21, 4219–4230. [Google Scholar] [CrossRef]
Figure 1. Height profiles of (a) the neutral wind factors of ions vertical drift velocity meridional C x (full line), zonal C y (dashed line) and vertical C z (dotted line) and (b) corresponding ions vertical convergence/divergence factors C x C x z (full line), C y C y z (dashed line) and C z C z z (dotted line). The horizontal dashed–dotted lines denote the heights of maximal values of ( C x ) max ( C z ) max (h = 121.4 km), minimal negative values of ( C y ) min < 0 (115 km) and maximal value of ( C y ) max > 0 (h = 131 km), and it is located at the height where the ion-neutral collisional frequency is equal to the ion cyclotron frequency ν i n = ω i (≈121 km). This value for the mid-latitude region ( I = 45 ° ) is the same as for meridional wind ( C x z C x ) max = sin I cos I 2 H [10]. In this case, the vertical changes of the zonal wind velocity factor C y z have its maximal (for westward wind V y > 0 ) ( C y z C y ) max = cos I 4 H or minimal (for eastward wind V y < 0 ) value ( C y z C y ) min = cos I 4 H . These maximal and minimal values are located at the heights ≈115 km and ≈131 km, which are below and above by about 0.9 H (H = 8–11 km) from the height region (≈121 km), where ν i n = ω i [12].
Figure 1. Height profiles of (a) the neutral wind factors of ions vertical drift velocity meridional C x (full line), zonal C y (dashed line) and vertical C z (dotted line) and (b) corresponding ions vertical convergence/divergence factors C x C x z (full line), C y C y z (dashed line) and C z C z z (dotted line). The horizontal dashed–dotted lines denote the heights of maximal values of ( C x ) max ( C z ) max (h = 121.4 km), minimal negative values of ( C y ) min < 0 (115 km) and maximal value of ( C y ) max > 0 (h = 131 km), and it is located at the height where the ion-neutral collisional frequency is equal to the ion cyclotron frequency ν i n = ω i (≈121 km). This value for the mid-latitude region ( I = 45 ° ) is the same as for meridional wind ( C x z C x ) max = sin I cos I 2 H [10]. In this case, the vertical changes of the zonal wind velocity factor C y z have its maximal (for westward wind V y > 0 ) ( C y z C y ) max = cos I 4 H or minimal (for eastward wind V y < 0 ) value ( C y z C y ) min = cos I 4 H . These maximal and minimal values are located at the heights ≈115 km and ≈131 km, which are below and above by about 0.9 H (H = 8–11 km) from the height region (≈121 km), where ν i n = ω i [12].
Atmosphere 14 01008 g001
Figure 2. The behavior of the electron density N e ( z , t ) / N o m in case of absence is neutral ( V ( V x ,     V y   , V z ) = 0 ) wind for its initial peak height h o m = 115     km (panel (a)), for the cases of influence of vertical wind with downward V z = 4   m / s (panel (b)), and upward component V z = 4   m / s (panel (c)), the one for h o m = 100     km (panel (d)), the same downward (panel (e)) and upward (panel f) velocities, respectively. The height h = 121 km corresponds to the maximal divergence rate C v min ( V x = 0 ,   V = 0 ,   V z = 4   m / s ) = 4.9 × 10 5   s 1 for downward wind V z = 4   m / s (blue dashed lines on panels (b,e)) and maximal convergence rate C v max ( V x = 0 ,   V = 0 ,   V z = 4   m / s ) = 4.9 × 10 5   s 1 for upward wind V z = 4   m / s (yellow dashed lines on, panels c,f).
Figure 2. The behavior of the electron density N e ( z , t ) / N o m in case of absence is neutral ( V ( V x ,     V y   , V z ) = 0 ) wind for its initial peak height h o m = 115     km (panel (a)), for the cases of influence of vertical wind with downward V z = 4   m / s (panel (b)), and upward component V z = 4   m / s (panel (c)), the one for h o m = 100     km (panel (d)), the same downward (panel (e)) and upward (panel f) velocities, respectively. The height h = 121 km corresponds to the maximal divergence rate C v min ( V x = 0 ,   V = 0 ,   V z = 4   m / s ) = 4.9 × 10 5   s 1 for downward wind V z = 4   m / s (blue dashed lines on panels (b,e)) and maximal convergence rate C v max ( V x = 0 ,   V = 0 ,   V z = 4   m / s ) = 4.9 × 10 5   s 1 for upward wind V z = 4   m / s (yellow dashed lines on, panels c,f).
Atmosphere 14 01008 g002
Figure 3. The HWM14 data for mid-latitude nightly varying (45° N; 45° S, day 102, t o = LT   23 : 00 ) height profiles of (a) the meridional wind velocity V x ( z , t ) and (b) the zonal wind velocity V y ( z , t ) .
Figure 3. The HWM14 data for mid-latitude nightly varying (45° N; 45° S, day 102, t o = LT   23 : 00 ) height profiles of (a) the meridional wind velocity V x ( z , t ) and (b) the zonal wind velocity V y ( z , t ) .
Atmosphere 14 01008 g003
Figure 4. (ac) The nightly varying ion vertical drift velocity w i v and (df) its vertical changes w i v z in case of presence (a,d) only meridional and zonal wind velocities determined by HWM14 data (see Figure 1, V z = 0 ), (panels (b,e)) the presence of vertical wind with V z = 4   m / s and (panels (c,f)) V z = 4   m / s . Yellow lines correspond to the height region with w i v = 0 (full lines) and C t max > 0 (dotted lines). Dashed blue lines correspond to the region C t min < 0 .
Figure 4. (ac) The nightly varying ion vertical drift velocity w i v and (df) its vertical changes w i v z in case of presence (a,d) only meridional and zonal wind velocities determined by HWM14 data (see Figure 1, V z = 0 ), (panels (b,e)) the presence of vertical wind with V z = 4   m / s and (panels (c,f)) V z = 4   m / s . Yellow lines correspond to the height region with w i v = 0 (full lines) and C t max > 0 (dotted lines). Dashed blue lines correspond to the region C t min < 0 .
Atmosphere 14 01008 g004
Figure 5. The behavior of the electron density N e ( z , t ) / N o m and formation of Es layers under the influence of neutral wind with velocity V ( V x ,     V y   , V z ) for ions initial distribution, Equation (17), with H i c = 16   km , peak height h o m = 115   km (panels ac) and h o m = 100   km (panels df) are demonstrated in the following cases: (1) only horizontal wind ( V z = 0 ) with meridional V x ( z , t ) (Figure 3a) and zonal V y ( z , t ) (Figure 3b) components determined by HWM14; (2) wind velocity vertical component directed downward with V z = 4   m / s (panels b,e); and (3) wind velocity vertical component directed upward with V z = 4   m / s (panels c,f). The full and dashed yellow lines correspond to the altitudes w i = 0 (Figure 4a–c) and ion convergence regions C t max > 0 (Figure 4d–f), respectively. The blue dashed lines correspond to the divergence regions C t min < 0 .
Figure 5. The behavior of the electron density N e ( z , t ) / N o m and formation of Es layers under the influence of neutral wind with velocity V ( V x ,     V y   , V z ) for ions initial distribution, Equation (17), with H i c = 16   km , peak height h o m = 115   km (panels ac) and h o m = 100   km (panels df) are demonstrated in the following cases: (1) only horizontal wind ( V z = 0 ) with meridional V x ( z , t ) (Figure 3a) and zonal V y ( z , t ) (Figure 3b) components determined by HWM14; (2) wind velocity vertical component directed downward with V z = 4   m / s (panels b,e); and (3) wind velocity vertical component directed upward with V z = 4   m / s (panels c,f). The full and dashed yellow lines correspond to the altitudes w i = 0 (Figure 4a–c) and ion convergence regions C t max > 0 (Figure 4d–f), respectively. The blue dashed lines correspond to the divergence regions C t min < 0 .
Atmosphere 14 01008 g005
Figure 6. The nightly variations of (a) maximal total convergence rates C t max > 0 (yellow dashed line), corresponding C v (dashed line) and C s h (full line) for the upper Es-type layer; (b) maximal total divergence rate C t min < 0 (blue dashed line), corresponding C v and C s h ; and (c) the convergence rates for lower Es-type layers (see Figure 5).
Figure 6. The nightly variations of (a) maximal total convergence rates C t max > 0 (yellow dashed line), corresponding C v (dashed line) and C s h (full line) for the upper Es-type layer; (b) maximal total divergence rate C t min < 0 (blue dashed line), corresponding C v and C s h ; and (c) the convergence rates for lower Es-type layers (see Figure 5).
Atmosphere 14 01008 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Didebulidze, G.G.; Dalakishvili, G.; Todua, M.; Toriashvili, L. The Role of Neutral Wind Velocity and Its Vertical Component on Predictability of Formation and Localization of Sporadic E (Es). Atmosphere 2023, 14, 1008. https://doi.org/10.3390/atmos14061008

AMA Style

Didebulidze GG, Dalakishvili G, Todua M, Toriashvili L. The Role of Neutral Wind Velocity and Its Vertical Component on Predictability of Formation and Localization of Sporadic E (Es). Atmosphere. 2023; 14(6):1008. https://doi.org/10.3390/atmos14061008

Chicago/Turabian Style

Didebulidze, Goderdzi G., Giorgi Dalakishvili, Maya Todua, and Lekso Toriashvili. 2023. "The Role of Neutral Wind Velocity and Its Vertical Component on Predictability of Formation and Localization of Sporadic E (Es)" Atmosphere 14, no. 6: 1008. https://doi.org/10.3390/atmos14061008

APA Style

Didebulidze, G. G., Dalakishvili, G., Todua, M., & Toriashvili, L. (2023). The Role of Neutral Wind Velocity and Its Vertical Component on Predictability of Formation and Localization of Sporadic E (Es). Atmosphere, 14(6), 1008. https://doi.org/10.3390/atmos14061008

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop