Mechanisms of stratification of the <Emphasis Type="Italic">F</Emphasis>2 layer and formation of the <Emphasis Type="Italic">F</Emphasis>3 and <Emphasis Type="Italic">G</Emphasis> layers in the equatorial ionosphere

The paper is dedicated to the studies of formation mechanisms of additional layers in the equatorial ionosphere carried out using numerical simulations with use of the Global Self-Consistent Model of the Thermosphere, Ionosphere, and Protonosphere (GSM TIP) modified in the part of the solution of the electric field equation in the Earth’s ionosphere. Calculations were preformed for quiet geomagnetic conditions using the MSIS-90 model for the calculation of thermospheric parameters. The obtained spatio-temporal pattern of thermospheric circulation and the variations in the dynamo electric field obtained on its basis make it possible to reproduce the stratification effect of the F2 layer and the appearance of the F3 layer in the equatorial ionosphere due to the action of the nonuniform in height zonal electric field at the geomagnetic equator. On the basis of the earlier presented results of calculations using the modified GSM TIP model, the appearance of a maximum in the vertical profile of the electron density at a height of ∼1000 km formed by H+ ions, which we called the G layer, has been predicted. Numerical simulations showed that this layer is formed by the meridional component of the thermospheric wind and is related to the formation of the nighttime midlatitude maximum at heights of the ionospheric F region.     

Behavior of <Emphasis Type="Italic">foF</Emphasis>2 and <Emphasis Type="Italic">hmF</Emphasis>2 after sunset

The time behavior of the foF2 and hmF2 values at the time moment T(ss + 2 h) 2 h after sunset is considered. It is assumed that at this moment, the horizontal winds in the thermosphere in the strongest way influence hmF2 and, therefore, foF2. It is found that a fairly well pronounced and statistically significant change (trend) is observed for the foF2(ss + 2)/foF2(14) ratio, the sign of the change being different for different stations and even different seasons at the same station. A similar picture is obtained for the value of hmF2(ss + 2). It is shown that a positive correlation between the trends of these two values is observed. This confirms the initial concept of the paper that the foF2 and hmF2 trends are caused by long-term trends in the thermospheric dynamics.     

Variability of parameters of the <Emphasis Type="Italic">F</Emphasis>2-layer maximum in the quiet midlatitude ionosphere under low solar activity: 1. Statistical properties

Results of statistical analysis of the properties of variability of F2-layer maximum parameters (critical frequency foF2 and the height hmF2) in quiet midlatitude ionosphere under low solar activity in the daytime (1000–1500 LT) and nighttime (2200–0300 LT) hours are presented on the basis of Irkutsk station data for 2007–2008. It is found that the distribution density of δfoF2 could be presented as consisting of two distinctly different normal laws of this distribution, one of which corresponds to weak (|δfoF2| < 10%) fluctuations in foF2 and the other corresponds to strong (30% > |δfoF2| > 10%) fluctuations. Weak fluctuations in foF2 to a substantial degree are related to ionospheric variability at times less of than 1–3 h and determine the δfoF2 variability in the daytime hours. Strong fluctuations in foF2 are mainly related to day-to-day variability of the ionosphere at a fixed local time, the variability increasing by approximately a factor of 3 during the transition from day to night and determining the δfoF2 variability in the nighttime hours. The distribution density of ΔhmF2 is close to the normal distribution law. An interpretation of the different character of the distribution densities of δfoF2 and ΔhmF2 is given.     

Numerical simulation of the <Emphasis Type="Italic">F</Emphasis>2-layer stratification and appearance of the <Emphasis Type="Italic">F</Emphasis>3 and <Emphasis Type="Italic">G</Emphasis> layers in the equatorial ionosphere: The morphology of the phenomena

This work is devoted to a numerical simulation of the equatorial ionosphere, performed using the GSM TIP model completed with a new block for calculating the electric field. It has been indicated that the usage of the wind system calculated according to the MSIS-90 model makes it possible to reproduce the electromagnetic drift velocities at the equator, the effect of the F2-layer stratification, and the appearance of the F3 layer in the equatorial ionosphere. The calculations performed using the modified GSM TIP model made it possible to detect a maximum in the electron density vertical profile at an altitude of ∼1000 km, formed by H⁺ ions, which we called the G layer. If this layer actually exists, it can be observed during sounding the low-latitude ionosphere from satellites during dark time of day.     

Variations in <Emphasis Type="Italic">foF</Emphasis>2 and <Emphasis Type="Italic">hmF</Emphasis>2 at the end of the 1990s and the beginning of the 2000s

This article considers sparse available data on variations in the main parameters of the ionospheric F2 layer foF2(ss + 2) and hmF2(ss + 2) at the end of the 1990s and the beginning of 2000s. It is shown that the monotonous behavior of hmF2(ss + 2) obtained in earlier publications for the period after 1980 is violated. The hmF2(ss + 2) behavior obtains a more complicated nature by time with a tendency towards a decrease in hmF2(ss + 2) at the beginning of a new century. A statistically significant relationship between foF2(ss + 2) and hmF2(ss + 2) is discovered confirming the Rishbeth statement that during the first hours after sunset, the critical frequency foF2 is governed by dynamical processes via changes in the F2-layer height. It is found that at the end of the interval in question, there is a tendency towards deviations from the above-mentioned dependence. The latter could be a manifestation of the fact that changes in the aeronomical parameters caused by the general cooling and contraction of the thermosphere begin influencing the foF2 value. It is found that in the summer months, the foF2(ss + 2) value demonstrated a systematic decline tendency from the “boundary date” towards the beginning of the 2000s.     

Relationship between variations in the characteristics of the high-latitude ionosphere and magnetic <Emphasis Type="Italic">PC</Emphasis> indices

It is shown in a joint analysis of ionospheric vertical sounding data at the arctic Heiss Island and antarctic Vostok stations and the geomagnetic PC index, which characterizes the geoefficient component of the interplanetary magnetic field, that, during a disturbed geomagnetic period when PC > 2 in years of solar activity (SA) maxima in the winter season, positive phases of ionospheric disturbances are predominantly observed. In the nighttime hours, an increase in the critical foF2 frequencies by a factor of 2–3 can occur. In a disturbed geomagnetic period at the PC > 1.5 level in the summer season, negative phases of ionospheric disturbances are mainly observed. In years of maximum and moderate SA, the decrease in foF2, as compared to their median values, happens at night (∼30%). In years of low SA, the decrease value is much lower. At a substantial decrease in the PC index level, in the region of the geomagnetic pole at the Vostok station, in some cases, a substantial increase in the electron density level in the F region occurs with a delay of 0.5 h. At the same time, a significant correlation (r = −0.57) is observed between variations in the PC index and foF2.     

Long-term changes in the relation between <Emphasis Type="Italic">foF</Emphasis>2 and <Emphasis Type="Italic">hmF</Emphasis>2

The change in the dependence of the F2-layer critical frequency on its height hmF2 is considered based on two sources of initial data used earlier by the authors. It is found that the slope k of the foF2 dependence on hmF2 systematically decreases from the earlier (“etalon”) period, 1958–1980, to the later periods of 1988–2010, 1998–2010, and 1998–2014. Since the foF2 value depends on the atomic oxygen concentration in the F region much more strongly than hmF2, the found decrease in k confirms the concept of a decrease in the atomic oxygen concentration in the thermosphere with time previously formulated by the authors.     

Regular changes in the critical frequency of the <Emphasis Type="Italic">F</Emphasis>2 layer of the quiet midlatitude ionosphere

A method for constructing the empirical model of the F2 layer critical frequency (foF2) under magnetically quiet conditions, aimed at analyzing disturbances of any nature, is proposed. This method has been analyzed, and typical features of regular changes in foF2 of the quiet ionosphere and day-to-day foF2 variability are analyzed using the data from Irkutsk and Slough stations as an example. In particular, it has been obtained that this model differs from the international IRI model, and this difference is mainly caused by the fact that the foF2 values in the IRI model do not correspond to quiet conditions. Therefore, this model gives a larger amplitude of the annual and semiannual variations in foF2 than the IRI model. In addition, this model more accurately reproduces the rate of foF2 annual variations at a fixed local time, especially in equinoxes, when foF2 variations can exceed 1 MHz within one month.     

Variations in statistical parameters of the <Emphasis Type="Italic">NmF</Emphasis>2 winter anomaly with latitude and solar activity

The maximal R ratios of the winter-to-summer NmF2 values of each ionosonde are calculated for a specified UT under daytime quiet geomagnetic conditions and at approximately equal levels of solar activity, based on foF2 measurement data of 98 ionosondes at mid- and low geomagnetic latitudes of the Northern and Southern hemispheres for 1957–2009. The P(R > 1) conditional probability of NmF2 winter anomaly observations, as well as the most probable RMP and average <R> of R values are calculated for low, moderate, and high solar activity on the base of foF2 measurements during the periods December 22 ± 30 days and June 21 ± 30 days. Variations in P(R > 1), RMP, and 〈R〉 with latitude and solar activity are analyzed.     

Variations in <Emphasis Type="Italic">foF</Emphasis>2 at the end of the 1990s and the beginning of the 2000s according to the SPIDR system data

Variations in the critical frequency of the F2 layer at 22 ionospheric stations within the period 1990–2010 according to the SPIDR system data are considered. A confirmation of the negative trends in foF2 for seven stations considered by one of the authors earlier on the basis of median data is obtained. It is found that both negative (a decrease in foF2 with time) and positive (growth in foF2 with time) trends of the critical frequency are observed. During the later part of the considered period (after 1997), negative trends dominate. This fact manifests itself, apparently, in an increase in the role of the decrease in the thermospheric neutral temperature in the formation of foF2 trends.     

Gradient drift instability as a source of spread <Emphasis Type="Italic">F</Emphasis> in the region of large-scale irregularities in the low-latitude topside ionosphere

The expression for the increment of instability and decrement of diffusion damping of gradient drift waves for ionospheric altitudes above the F ₂ layer maximum is obtained. The gradient drift instability is used to interpret the observations of spread F in the region of large-scale horizontal irregularities of the electron density. Two types of such irregularities observed on board the Intercosmos-19 (IC-19) satellite in the region of low latitudes (a peak of the density in the dusk ionosphere and a trough of the density in the dawn ionosphere) are considered. It is shown that the observed gradients of the density and electric field values in the dawn and dusk ionospheric sectors are quite sufficient for the instability development criterion to be satisfied in both considered cases.     

Relationship between the daytime critical frequencies of the <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript> layer and thermospheric characteristics

The relationship between the critical frequency of the F ₂ layer and the atmospheric characteristics has been obtained in a general form. It has been shown that this relation makes it possible to sufficiently accurately describe the daytime values of foF2 while comparing them with the observed monthly median values. Such comparisons were performed, first, for the data of measurements in Irkutsk using the DPS-4 digital ionosonde in 2003–2006 and, second, based on the annual variations in the noon foF2 values at 24 stations of the Northern Hemisphere in 1984. The calculations were performed using the MSIS-86 thermospheric model [Hedin, 1987].     

Reaction of the critical frequency of the <Emphasis Type="Italic">F</Emphasis>2 layer to a sharp depletion in atmospheric pressure

Changes in the critical frequencies of the F2 layer at several midlatitude stations of ionospheric vertical sounding during a sharp depletion in atmospheric pressure under quiet solar and geomagnetic conditions are analyzed. It is shown that in such periods, the observed foF2 values differ from the mean values by approximately 10–15% and the deviations from the mean could be both negative (in the daytime hours) and positive (at night). Such variations in foF2 could be referred to the known class of ionospheric disturbances observed under a quiet geomagnetic situation, that is, to the so-called “Q-disturbances.” Analysis of wavelet spectra of foF2 variations shows the presence in the F region of oscillations of various periods (from 0.5 to 10 days). The decrease in the amplitude of daily variations during pressure depletion is found. Presumably, the observed effect is caused by the dynamic impact of waves formed in the lower atmosphere on the ionospheric F2 layer.     

Properties of the <Emphasis Type="Italic">F</Emphasis>2-layer critical frequency median in the nocturnal subauroral ionosphere during low and moderate solar activity

Based on an analysis of data from the European ionospheric stations at subauroral latitudes, it has been found that the main ionospheric trough (MIT) is not characteristic for the monthly median of the F2-layer critical frequency (foF2), at least for low and moderate solar activity. In order to explain this effect, the properties of foF2 in the nocturnal subauroral ionosphere have been additionally studied for low geomagnetic activity, when the MIT localization is known quite reliably. It has been found that at low and moderate solar activity during night hours in winter, the foF2 data from ionospheric stations are often absent in the MIT area. For this reason, a model of the foF2 monthly median, which was constructed from the remaining data of these stations, contains no MIT or a very weakly pronounced MIT.     

Spatial and time variations in critical frequencies of the ionospheric <Emphasis Type="Italic">F</Emphasis> region above the zone of equatorial earthquake preparation

Spatial and time variations in the critical frequencies foF2 before the strong earthquake of August 15, 1963, with the magnitude M = 7.75 are analyzed. The epicenter of the earthquake was located in the vicinity of the magnetic equator in the American longitudinal sector. The data of the topside ionosonde on board Alouette-1 and of the series of ground-based ionosondes has been used for this purpose. The ground-based ionosondes made it possible to detect an insignificant anomalous decrease in foF2 within the zone of earthquake preparation a few days prior to the earthquake. This result confirms the conclusion drawn earlier on the basis of the satellite data. The modification of the ionosphere at the F-region level is more evident in the satellite than in the ground-based data. It is also noted that the character of the time variations in foF2 a day before the earthquake is similar to the so called “quiet time” Q-disturbances in the ionosphere, when the electron concentration at the F region maximum differs from the median values by more than 20% under undisturbed geophysical conditions.     

Critical frequencies <Emphasis Type="Italic">foF</Emphasis>2 as an indicator of trends in thermospheric dynamics

The time variations in three parameters during the last decades are considered. R(foF2) is the correlation coefficient between the nighttime and daytime values of foF2 for the same day. Stable trends are found for the minimum (R(foF2)(max)) and maximum (R(foF2)(min)) values of R(foF2) during a year. The foF2(night)/foF2(day) ratio demonstrates both, negative and positive trends, and the trend sign depends on the inclination I and declination D of the magnetic field. The correlation coefficient r(h, fo) between foF2 and the 100 hP level in the stratosphere demonstrates a decrease (in the years of maximum and minimum solar activity) from the 1980s to the 1990s. The trends in all three groups of data are considered under the assumption of long-term changes in the circulation in the upper atmosphere.     

Dispersion Measurements of <Emphasis Type="Italic">P</Emphasis> Waves and their Implications for Mantle <Emphasis Type="Italic">Q</Emphasis><Subscript><Emphasis Type="Italic">p</Emphasis></Subscript>

We analyze the anelasticity of the earth using group delays of P-body waves of deep (>200 km) events in the period range 4–32 s for epicentral distances of 5–85 degrees. We show that Time Frequency Analysis (TFA), which is usually applied to very dispersive surface waves, can be applied to the much less dispersive P-body waves to measure frequency-dependent group delays with respect to arrival times predicted from the CMT centroid location and PREM reference model. We find that the measured dispersion is due to: (1) anelasticity (described by the P-wave quality factor Q _p ), (2) ambient noise, which results in randomly distributed noise in the dispersion measurements, (3) interference with other phases (triplications, crustal reverberations, conversions at deep mantle boundaries), for which the total dispersion depends on the amplitude and time separation between the different phases, and (4) the source time function, which is dispersive when the wavelet is asymmetrical or contains subevents. These mechanisms yield dispersion ranging in the order of one to 10 seconds with anelasticity responsible for the more modest dispersion. We select 150 seismograms which all have small coda amplitudes extending to ten percent of the main arrival, minimizing the effect of interference. The main P waves have short durations, minimizing effects of the source. We construct a two-layer model of Q _p with an interface at 660 km depth and take Q _p constant with period. Our data set is too small to solve for a possible frequency dependence of Q _p . The upper mantle Q ₁ is 476 [299–1176] and the lower mantle Q ₂ is 794 [633–1064] (the bracketed numbers indicate the 68 percent confidence range of Q _p ^–1). These values are in-between the AK135 model (Kennett et al., 1995) and the PREM model (Dziewonski and Anderson, 1981) for the lower mantle and confirm results of Warren and Shearer (2000) that the upper mantle is less attenuating than PREM and AK135.     

Behavior of critical frequency <Emphasis Type="Italic">foF2</Emphasis> at various moments of the day: 1. Dependence on solar activity

Variation with solar activity level of the ratios of the critical frequencies of the F2 layer, foF2, is considered for various pairs of local time moments T1 and T2 for two seasons (winter and summer). It is found that, as a rule, the foF2(T1)/foF2(T2) behavior with the solar activity index F10.7 is different in winter and summer. The variations of foF2(T1)/foF2(T2) with F10.7 obtained from experimental data are considered from the viewpoint of the current theory of F2-layer formation. It is shown that the majority of observed characteristics of these variations are explained using the above-indicated theory. However, special cases of foF2(T1)/foF2(T2) behavior with solar activity are found for which there is still no physical explanation.