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.     

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.     

Variations of the critical frequency <Emphasis Type="Italic">foF</Emphasis>2 at middle latitudes during strong disturbance in the electric field observed in the <Emphasis Type="Italic">F</Emphasis> region over saint-santin

Morphological analysis of variations of the critical frequency foF2 in the midlatitude ionosphere at various sectors of local time is carried out on the basis of data from ground-based stations of vertical sounding of the ionosphere in the period when during use of the incoherent scatter radar at Saint-Santin an anomalously strong increase in the electric field was observed at heights of the ionospheric F region in the period of enhanced geomagnetic activity (4+ < Kp < 6−). The obtained picture of the space-time distribution of disturbances in foF2 makes it possible to assume that they could be caused by penetration to middle latitudes of the large-scale electric field of the magnetospheric convection directed westward in the nighttime and morning hours and eastward in the noon and evening sectors.     

Spatial and time variations in the <Emphasis Type="Italic">foF</Emphasis>2 (night)/<Emphasis Type="Italic">foF</Emphasis>2 (day) ratio: Specification of some effects

A series of important aspects of the time and spatial variations in the critical frequency ratio fo(night)/fo(day) considered by one of the authors [Danilov, 2007, 2008] is specified. The list of the Eastern-Hemisphere stations, for which an analysis of the above indicated ratio was performed, is completed. The available stations of the Western Hemisphere are considered. It is shown that the character of the variations in the fo(night)/fo(day) ratio is independent of limitations imposed on the Ap index of magnetic activity for the analyzed days. The consideration of the long-term variations in geomagnetic activity using the annual mean value of the Ap index does not influence the principal conclusions of the study, although in some cases changes statistical characteristics of the obtained behavior of the fo(night)/fo(day) ratio after 1980. A comparison of the results, based on the interpretation of the trends of the fo(night)/fo(day) value with the changes in the stratosphere-ionosphere coupling from the 1980s to the 1990s obtained by the authors earlier, confirms the assumption that there occurs a systematic change in the zonal wind in the upper atmosphere.     

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.     

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.     

<Emphasis Type="Italic">R/S</Emphasis> analysis of the <Emphasis Type="Italic">Dst</Emphasis> index

The paper addresses estimation of the Hurst exponent for time series of the hourly values of the Dst index for the period from 1957 to 2011. It is found that the Hurst exponent is 0.79–0.94 for yearly intervals and 0.8–1.0 for monthly intervals. Based on R/S graphs, the Dst cycles are identified; they range from 3–4 months to 2.2 years and from 8.5 to 22 years in length. It is shown that a Dst time series can be quite satisfactorily described by an α-stable Levy process.     

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.     

Long-term trends in the ratio of the daytime and nighttime values of <Emphasis Type="Italic">foF</Emphasis>2

The consideration of the relation between the daytime and nighttime values of the critical frequency F2, foF2 of the ionospheric F2 layer, started in the previous publication of the authors, is continued. The main regularities in variations in the correlation coefficient R(foF2) characterizing this relation are confirmed using larger statistical material (more ionospheric stations and longer observational series). Long-term trends in the R(foF2) value are found: at all stations the negative value of R(foF2) increases with time after 1980.     

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.     

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.     

Multi-scale Tomography for Crustal <Emphasis Type="Italic">P</Emphasis> and <Emphasis Type="Italic">S</Emphasis> Velocities in Southern California

Seismic tomography is a viable tool in building depth-velocity models in the presence of strong lateral velocity variations. In this study 3-D P- and S-velocity models for the crust of southern California are constrained using more than 1,000,000 P-wave first arrivals and 130,000 S-wave arrivals from local earthquakes. To cope with the uneven distribution of raypaths, a multi-scale tomography is applied with overlapping model cells of different sizes. Within the 300 × 480 × 39 km³ model volume, the smallest cell size is 10 × 10 × 3 km³. During the iterations of velocity updating, earthquake hypocenters are determined using both P and S arrivals, and full 3-D ray tracing is implemented. Except near the edges and in the lower crust, the resultant models are robust according to various tests on the effects of reference models, resolution and signal-to-noise ratio. The tomographic velocities at shallow depths correlate very well with the regional geology of southern California. In the upper crust the P-wave and S-wave models exhibit slow velocities in major sedimentary basins and fast velocities in areas of crystalline rocks. Mid-crustal low velocity zones are present under the Coso Range, San Gabriel Mountains, and a large portion of the Mojave Desert. P- and S-velocity patterns maintain their similarity in the lower crust though the models are less reliable there.     

Winter anomaly of the <Emphasis Type="Italic">E</Emphasis>-layer critical frequency in the nighttime auroral zone

Analysis of the annual variation of the E-layer critical frequency median foE in the nighttime (22?02 LT) auroral zone by the data of several stations of the Northern Hemisphere has shown the median maximum in winter and minimum in summer, even though the summer contribution of solar radiation to foE is greater. Thus, a new phenomenon was discovered—an foE median winter anomaly in the nighttime auroral zone. Its amplitude (ratio of winter to summer foE figures) can reach 10–15%; however, this anomaly was weakly expressed and statistically insignificant at particular stations located in the auroral zone. The winter anomaly is more distinct for foE _avr, the median of the E-layer critical frequency foE caused by the auroral source of atmospheric ionization, i.e., excluding the solar radiation contribution to foE. For foE _avr, the amplitude of the winter anomaly can reach 15–20%. Based on the qualitative analysis, it has been found that foE winter anomaly is stipulated by the winter/summer asymmetry of energy flow of accelerated electrons, which induce discrete aurorae in the nighttime auroral zone.     

Variations in the fractal dimension of fluctuations of the ionospheric <Emphasis Type="Italic">F</Emphasis>2 layer critical frequency

The 40-year period of observations of short-term variations (with characteristic times of up to 1–2 days) in the critical frequency of the ionospheric F2 layer (foF2) is analyzed. The continuous (with a step of 1 h) series of fluctuations (F) of the foF2 critical frequency (with eliminated daily variations) has been calculated using the hourly variations in foF2 at Moscow stations. The fractal dimension (FRH) of the fluctuations, characterizing short-term variations in foF2, has been determined and analyzed on a 30-day interval, using the Higuchi method. It has been established that FRH estimates substantially change in time. The 11-year cycle, which is in antiphase with the solar cycle, and the total annual and semiannual variations, similar to the variations observed in the normalized critical frequency of the E region and in the electron density of the D region, are clearly defined in these changes. Thus, the parameters of fast variations in the ionospheric F2 layer are affected by the phase of the 11-year solar cycle and by the position of the Earth in the orbit or seasonal variations in the atmosphere.     

Vertical structure of the midlatitude winter <Emphasis Type="Italic">F</Emphasis> region of the ionosphere during postmidnight enhancements in <Emphasis Type="Italic">NmF</Emphasis>2

The behavior of the vertical structure of the ionospheric F2 layer, including the variations in the heights of the maximum and bottom of the layer, its half-thickness, and electron content at some fixed heights during postmidnight enhancements in the electron density at the F2 layer maximum (NmF2), has been studied based on the data of the ionospheric vertical sounding, conducted in Alma-Ata (76°55′ E, 43°15′ N) in 2005–2006. The analysis of the amplitude and phase relationships between the measured parameters of the layer made it possible to qualitatively complete the existing concepts of the mechanisms by which the discussed effect is maintained. It is shown that the accelerated decrease in the electron density of the layer within a short time interval preceding the beginning of the postmidnight increase in NmF2 is governed not only by recombination processes but also by the plasma redistribution over the increasing thickness of the layer. The regularly observed effect of the delay in the moment of reversal in the motion direction of the layer bottom relative to the corresponding moment for the layer maximum made it possible to conclude that the meridional wind asynchronously reverses its direction from the poleward daytime to the equatorward nighttime in the entire layer: the direction changes later with decreasing height.     

Combined Method of <Emphasis Type="Italic">F</Emphasis>-, <Emphasis Type="Italic">S</Emphasis>-, and <Emphasis Type="Italic">R</Emphasis>-Approximations of Increased Dimensionality in Solving the Problems of Geophysics and Geomorphology

The interrelation between different variants of the method of linear integral representations in the spaces of an arbitrary dimension is considered. The combined approximations of the topography and geopotential fields allows the selection of the optimal parameters of the method in solving a wide range of inverse problems in geophysics and geomorphology, as well as a most thorough use of the a priori information about the elevations and elements of anomalous fields. A method for numerically solving an inverse problem on finding the equivalent, in terms of the external field, mass distributions in the ordinary three-dimensional (3D) space and in the four-dimensional (4D) space is described.     

Combined method of <Emphasis Type="Italic">F</Emphasis>-, <Emphasis Type="Italic">S</Emphasis>-, and <Emphasis Type="Italic">R</Emphasis>-approximations in solving the problems of geophysics and geomorphology

The interrelation between different modifications of the method of linear integral representation is studied. Combined approximations of the topography and geopotential fields enable more refined tuning of the method in solving inverse problems of geophysics and geomorphology and provide a more complete allowance for the a priori information about the surface elevation data and elements of anomalous fields. A technique for finding the numerical solution for the inverse problem for determining the mass distributions equivalent in terms of the external field is presented. The results of the mathematical experiment are discussed.     

Mechanism of formation of <Emphasis Type="Italic">Q</Emphasis>-disturbances in the <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript> region of the equatorial ionosphere

Using model simulations, the morphological picture (revealed earlier) of the disturbances in the F ₂ region of the equatorial ionosphere under quiet geomagnetic conditions (Q-disturbances) is interpreted. It is shown that the observed variations in the velocity of the vertical E × B plasma drift, related to the zonal E _y component of the electric field, are responsible for the formation of Q-disturbances. The plasma recombination at altitudes of the lower part of the F ₂ region and the dependence of the rate of this process on heliogeophysical conditions compose the mechanism of Q-disturbance formation at night. The daytime positive Q-disturbances are caused exclusively by a decrease in the upward E × B drift, and this type of disturbances could be related to the known phenomenon of counter electrojet. Possible causes of formation of the daytime negative Q-disturbances are discussed.     

Linear transformations of anomalous potential fields using <Emphasis Type="Italic">S</Emphasis> and <Emphasis Type="Italic">R</Emphasis> approximations

Methods for linear transformations of anomalous physical fields based on R and S approximations of the elements of the initial field are described. The results of the numerical experiment and analytical extension of the gravity field measured in a detailed gravity survey in two regions of Russia are presented.     

Attenuative Dispersion of <Emphasis Type="Italic">P</Emphasis> Waves and Crustal <Emphasis Type="Italic">Q</Emphasis> in Turkey and Germany

We have measured group delays of the spectral components of high-frequency P-waves along two portions of the North Anatolian Fault Zone (NAFZ) in Turkey and in a region of southern Germany. Assuming that the observed dispersion is associated with attenuation in the crust and that it can be described by a continuous relaxation model, we obtained Q and the high-frequency relaxation times for those waves for each of the three regions. Individual P-wave Q values exhibit large scatter, but mean values in the NAFZ increase from about 25 to 60 over the distance range 5–90 km. Mean Q values are somewhat higher in the eastern portion of the NAFZ than in the western portion for measurements made at distances between 10 and 30 km. P-wave Q values in Germany range between about 50 and 300 over the hypocentral distance range 20–130 km. In that region we separated the effects of Q for basement rock (2–10 km depth) from that of the overlying sediment (0–2 km depth) using a least-squares method. Q varies between 100 and 500 in the upper 8–10 km of basement, with mean values for most of the distance range being about 250. Q in the overlying sediments ranges between 6 and 10. Because of large scatter in the Q determinations we investigated possible effects that variations of the source-time function of the earthquakes and truncation of the waveform may have on Q determinations. All of our studies indicate that measurement errors are relatively large and suggest that useful application of the method requires many observations, and that the method will be most useful in regions where the number of oscillations following the initial P pulse is minimized. Even though there is large scatter in our Q determinations, the mean values that we obtained in Turkey are consistent with those found in earlier studies. Our conclusions that Q is significantly higher in the basement rock of Germany than in the basement rock of Turkey and that Q is lower in western Turkey than in eastern Turkey are also consistent with results of Q studies using Lg coda.