Issue 
J. Space Weather Space Clim.
Volume 3, 2013



Article Number  A30  
Number of page(s)  28  
DOI  https://doi.org/10.1051/swsc/2013052  
Published online  15 October 2013 
Research Article
Model of traveling ionospheric disturbances
Kharkiv V. N. Karazin National University, Kharkiv, Ukraine
^{*} Corresponding author: FedorenkoYP@mail.ru
Received:
19
September
2011
Accepted:
5
September
2013
A multiscale semiempirical model of traveling ionospheric disturbances (TIDs) is developed. The model is based on the following assumptions: (1) TIDs are generated by acousticgravity waves (AGWs) and propagate as pressure waves; (2) time intervals between adjacent extrema of atmospheric pressure oscillations in a disturbance source are constant; (3) the pressure extrema propagate from the source up to ~14 000 km at a constant horizontal velocity; (4) the velocity of each extremum is determined only by its number in a TID train. The model was validated using literature data on disturbances generated by about 20 surface and highaltitude nuclear explosions, two volcano explosions, one earthquake and by energetic proton precipitation events in the magnetospheric cusp of the northern hemisphere. Model tests using literature data show that the spatial and temporal TID periods may be predicted with an accuracy of 12%. Adequacy of the model was also confirmed by our observations collected using transionospheric sounding. The following TID parameters: amplitudes, horizontal spatial periods, and a TID front inclination angle in a vertical plane are increasing as the distance between an AGW and the excitation source is increasing. Diurnal and seasonal variability of the TID occurrence, defined as ratio of TID events to the total number of observations for the corresponding period, is not observed. However, the TID occurrence was growing from ~50% in 1987 to ~98% in 2010. The results of other studies asserting that the TID occurrence does not depend on the number of sunspots and magnetic activity are confirmed. The TID occurrence has doubled over the period from 1987 to 2010 indicating increasing solar activity which is not associated with sunspot numbers. The dynamics of spatial horizontal periods was studied in a range of 150–35 000 km.
Key words: semiempirical model / multiscale traveling ionospheric disturbances / horizontal velocity of extremum / nuclear explosion / volcano explosive eruption / energetic proton precipitation / magnetospheric cusp / transionospheric sounding / disturbance occurrence / sunspot number / magnetic activity
© Y.P. Fedorenko et al., Published by EDP Sciences 2013
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
Traveling ionospheric disturbances (TIDs) are known to be generated by acousticgravity waves (AGWs) propagating in the neutral atmosphere. The TIDtrack, i.e., the spatial quasiperiodic modulation of electron concentration N along a surface concentric with the earth’s surface, leads to modification of radio wave parameters on propagation paths crossing this track. Therefore, a predictive TID model, updatable by current measurements, may be useful for practical applications (e.g., for reduction of navigation measurement errors and for radio wave propagation forecasts).
AGWs of natural origin play an important role in energy and momentum transfer from high latitudes to mid and low latitudes during precipitation of solar wind protons into the magnetospheric cusp (RTO Technical Report 2009; Karpachev et al. 2010) and from equatorial latitudes to higher latitudes during tropical storms (Vadas & Crowley 2010). Therefore, a study of AGWTID phenomena is of great scientific value, allowing deeper understanding of the dominant energy redistribution mechanisms in the atmosphere.
Depending on the wavelength (spatial period) L, TIDs are usually separated into two classes: largescale (LS) and mediumscale traveling ionospheric irregularities (MS TID). The LS TIDs are characterized by wavelengths L > 1000 km, temporal periods T ≈ 0.5–4 h, and horizontal velocities V ≈ 0.4–1000 km s^{−1} (Hocke & Schlegel 1996; RTO Technical Report 2009). The corresponding parameters for MS TIDs are L ≈ 100–600 km, T ≈ 0.25–1 h, and V < 0.25 km s^{−1}. It was predicted theoretically (Hocke & Schlegel 1996) and proven experimentally (Fedorenko et al. 2008, 2010a, 2010b) that minimal wavelengths L for MS TID do not exceed 100–150 km. Research carried out within the last six decades has revealed that MS TIDs occur much more frequently than LS TIDs. With rare exceptions, LS TIDs are moving from polar regions to the equator, they are observed often, but not always, after auroral disturbances during magnetic storms (Francis 1975; Hocke & Schlegel 1996). Some scientists assume that the abovementioned TID classification is defined by the nature of their origin. However, our longterm experimental studies, using satellite beacons, allow us to conclude that both LS and MS traveling ionospheric disturbances are of the same nature and really present different phases of disturbance dynamics. Nevertheless, for convenience we will use the LS and MS terminology in this paper.
Due to their large horizontal scales, LS TIDs can only propagate in the earthionosphere waveguide. MS TIDs can move parallel to the earth surface in guided wavemode (e.g., Gossard & Hooke 1975) as well as in free wavemode (e.g., Francis 1974). In the latter case, TID waves are reflected by the Earth’s surface. Until now, there exists no generally accepted opinion concerning MS TID propagation and an effective layer thickness where disturbances are observed. For example, according to Iyer (1983) and Senior et al. (2006) the layer thickness varies in the range 100–120 km, but our research has shown that the thickness was only about 50 km in most cases, and values ~100 km were observed only occasionally (Dorohov et al. 2004; Fedorenko et al. 2008, 2010a, 2010b, 2011).
Now, after 66 years since TIDs discovery (Munro 1948), the theory of interaction between AGWs and TIDs is far from being complete, some predictions and deductions have not been confirmed experimentally and have not been understood completely yet. In particular, the growth of wavelength L in a TID train with increasing distance from an AGW source, predicted by existing theoretical models (Francis 1974; Bristow et al. 1994), has not been observed experimentally. According to the abovementioned theoretical approaches, the amplitude A _{ d } of a TID train as well as L grow when the distance from a disturbance source increases. Here A _{ d } = ∆N/N where ∆N is the difference between background and perturbed value of N. Increase of the amplitude A _{ d } during TID propagation was demonstrated in recent experimental research (Tsugawa et al. 2007a), however, the authors could not explain the result. These effects have been clearly discovered and analyzed in our numerous experiments carried out under specific heliogeophysical conditions (see Sect. 5). Due to limited measurement accuracy, there exists a misconception among researchers (see, for example, Evans et al. 1983; Ogawa et al. 1988; Tsugawa et al. 2007a, 2007b; Kotake et al. 2007; Song et al. 2011; Otsuka et al. 2013), concerning diurnal and seasonal occurrence of TIDs (see Sect. 5). According to our observations, there are no minima and maxima in diurnal and seasonal TID occurrence. Further, our observations, carried out during 1987–2010, revealed that TID occurrence does not depend on the phases of the solar activity cycle.
The invariability of the propagation velocity of the wave train extrema was deduced from theoretical treatments and was confirmed experimentally 45 years ago for TIDs generated by nuclear tests and highpower volcanic explosions (Row 1967; Breitling & Kupferman 1967; Kanellakos 1967; Albee & Kanellakos 1968; Lomax & Nielson 1968; Roberts et al. 1982a, 1982b; Igarashi et al. 1994). However, some scientists can hardly accept the fact that this velocity is practically constant up to 14 000 km. Therefore, up to now this phenomenon was not explained. At the same time, this effect plays a key role in the unified semiempirical model of AGWTID, presented in this article.
During about 20 years after TID discovery in 1946 and the first publication (Munro 1948), which described these experiments, a theoretical AGWTID model has been developed (Row 1967; Cole & Greifinger 1969; Liu & Yeh 1971; Yeh & Liu 1974; Francis 1974; Gossard & Hooke 1975). Many authors have created this model, or rather its varieties, and throughout decades geophysicists (e.g., Bristow et al. 1994) repeated its main features.
The model is based on linearized equations of geophysical hydrodynamics for an isothermal nondissipative atmosphere and a point perturbation source. Winds were taken into account only if TIDs propagate in waveguide mode (WGM) (Gossard & Hooke 1975). The wind velocity and direction were also taken into consideration in a transformation process from freewave mode (FWM) to waveguide mode. However, almost all these assumptions are violated in processes of AGW excitation and propagation. For example, in nuclear tests, highpower volcano explosive eruptions and earthquakes the overall dimensions of AGW sources may vary from one to ten kilometers. Energetic proton precipitation into the magnetosphere at 60° magnetic latitude can form a cusp whose length may vary from ~8000 to 10 000 km, while a cusp length in a night sector is extended along Lshell up to 10–12 time zones (Halperin et al. 1990). Therefore, it is no surprise that a TID front, generated by a highlatitude source, has a length that exceeds 2000 kilometers along a parallel in a zonal direction (Fedorenko et al. 1997; Tsugawa et al. 2007a, 2007b; Kotake et al. 2007; Song et al. 2011; Otsuka et al. 2013). A prior guess, asserting that relative changes of atmospheric density at distances of hundreds and thousands of kilometers from the source of local energy release, is much less than unity is quite courageous. Numerous experiments indicate that such approach is not valid, since the TID amplitude A _{ d } may vary from ten to several ten percent at large distances from the source and often reaches 70–80% (Breitling & Kupferman 1967; Kanellakos 1967; Albee & Kanellakos 1968; Lomax & Nielson 1968; Taran et al. 1979; Afraimovich et al. 2005; Tyrnov et al. 2005; Yakovets et al. 2008; Fedorenko et al. 2010a, 2010b, 2011, 2012). Thus, these approximations may be applicable and reasonable only at the first stage of research.
Surprisingly, the past 60 years of investigations have proven that such a simplified TID model, grounded on simplifying and often not fulfilled assumptions, was generally confirmed by experimental data. Although several of its predictions, for example, spatiotemporal dynamics of the AGW train, have not been confirmed experimentally. But some results, reliably confirmed by experiments, have been simply ignored by investigators, probably due to the surprising nature of these results. In particular, the abovementioned facts, concerning propagation velocity and amplitudes of the AGWtrain, have not been analyzed. TID properties are essentially very similar to properties of solitary waves, generated by nonlinear processes (e.g., Lonngren & Scott 1978; Bullough & Caudrey 1980). A single maximum of solitary waves moves with constant velocity, while each of the AGW extrema moves with its own constant velocity. It is commonly accepted that both AGWTIDs and solitary waves propagate over large distances without essential modification of their form and amplitude. A solitarylike TID behavior was once recorded when one TID transited through another, and both remained unchanged (Heisler 1958). Therefore, we suppose that the theory of AGWTID origin and propagation may be based on the soliton theory. Summing up this short critical review, we may conclude that a complete AGWTID theory does not exist. There exist only two reasonable hypotheses, but their application to particular TID observations sometimes does not give good results.
The nature of regularly observed TIDs, generated by highlatitude natural energy sources, has not been disclosed yet. For years scientists all over the world have been trying to solve this problem. It is worth mentioning the campaigns for simultaneous TID observations at the high, middle, and low latitudes in both hemispheres by all existing diagnostic means (Hocke & Schlegel 1996). Over the years of observation, there appeared many theories, describing generation mechanisms of regularly observed AGWs. Nonetheless morphological investigations conducted over the last 15 years have facilitated identification of disturbance sources. For example, only the most probable source, namely energetic particle precipitation into the magnetospheric cusp, was mentioned in RTO Technical Report (2009), but the Joule heating by particle precipitation had been mentioned by investigators many times before this report was published. We have also made the same conclusion. It was reported by Karpachev et al. (2010) that the disturbance source, responsible for the LS AGW, recorded by various methods during a substorm observed in 1998, appears to be located in the cusp region.
In our model we combine wellknown parameters and relations, obtained in the past by various researchers. Our model links different results from theoretical and experimental morphological studies and can predict spatiotemporal dynamics of TID parameters. Certainly, our model is not yet complete. The model is described in Section 2, selection and justification of model parameters are given in Section 3. Sections 4 and 5 are devoted to model verification using experimental data, taken from available literature and our experimental results, obtained by transionospheric sounding by satellite signals.
Thus, the present article aims at the development of a semiempirical model of AGWTID generation and propagation and its comprehensive experimental verification.
2. Description of the AGWTID model
Our semiempirical model is based on the following assumptions:

atmospheric density oscillations are arising under the influence of an external stimulus;

time intervals between adjacent pressure extrema are constant;

each pressure extremum propagates from a source as a wave with a horizontal velocity which is constant up to ~14 000 km;

the speed of each extremum is determined by its index in a TID train.
2.1. The source parameters
The AGW source is defined as a volume whose characteristic dimension is equal to ~5 km. Atmospheric pressure oscillates in this volume due to local release of energy and momentum. Oscillations start with an intense decrease of the air pressure in the volume. The initial disturbance, generated by energetic proton precipitation in the magnetospheric cusp, is described as an extended source, consisting of the abovementioned volume sources, aligned along geomagnetic latitudes in the vicinity of 60–75° at an altitude of 100–200 km. The time interval between an adjacent pressure minimum and maximum, i.e., the temporal halfperiod, is τ _{0} = 3.95 min. Each extremum generates an atmospheric pressure wave propagating along the Earth’s surface with velocity V _{ i }. The velocity of each extremum does not vary over distances up to 14 000 km. For sphereshaped sources the wavefront section in a horizontal plane is a circle. If the source is proton precipitation in the magnetospheric cusp, a section is a line or an arc along a magnetic parallel, oriented along the cusp longitudinal axis. In both cases influence of winds, regional irregularities of temperature and pressure were neglected.
The velocities V _{ i } and wave amplitudes A _{ i }, of the first 16 extrema, normalized by peak value, are given in Table 1. We postulate that A _{ i } and V _{ i } do not vary during TID propagation.
Amplitudes A _{ i } and velocities V _{ i } for extremum index i.
2.2. AGWTID parameters
All AGW parameters do not depend upon the geomagnetic field, but TID amplitudes depend on the angle between the AGW wave vector and a plane through the magnetic meridian, and upon the magnetic dip angle in the vicinity of the measuring station. Albee & Kanellakos (1968), Lomax & Nielson (1968) have ascertained that if an AGW source is located in the southern magnetic hemisphere and TID measurements are performed southward of the AGW source or northward of the magnetic equator, AGW train maxima (minima) give rise to electron concentration maxima (minima), i.e., N _{ m } increases (decreases) relative to the background value. Thus, an altitudinal profile N(z) is shrunk and z _{ m } is decreased.
If TID observations are performed to the north of a source, but to the south of the magnetic equator, the AGW train maxima (minima) give rise to electron concentration minima (maxima) and the altitude profile of N(z) stretches out. In all these cases, the total electron content N _{0} is practically unchanged.
Oscillations of ionized atmospheric components, entrained by an AGW in a direction transverse to magnetic field lines, do not lead to the generation of TIDs. The Lorentz force makes charged particles move only along the geomagnetic field line. In the vicinity of the North geographic pole the magnetic field is directed toward the Earth surface, while in the southern hemisphere it is directed upward. According to our theoretical analysis, carried out without involvement of a model of the Earth’s electric field, air pressure maxima (minima) in an AGW correspond to maxima (minima) of N _{ m } in a TID, but application of this model will, probably, explain all observable effects. At the same time, the whole range of the abovementioned experimental observations has not been explained yet. The experimental data collected in the same observation cycle in both hemispheres for the same initiating source are in fact only from the four US nuclear explosions. They were produced in 1962 in the Pacific Ocean near the magnetic equator. From these data, the following empirical relationships are deduced.
If a disturbance source is located in the northern hemisphere, the sign of ΔN (the difference between TID extremum N _{ m } and background N values) is defined by the sign of the following empirical relationship (Albee & Kanellakos 1968):(1) where θ is the magnetic dip angle which is positive in the northern hemisphere and negative in the southern hemisphere, φ _{ TID } is the angle between the AGW wave vector and the magnetic north direction at the observation point. Positive values of J _{ΔN } are proportional to ΔN/N _{ m }. For negative J _{ΔN } the value (ΔN/N _{ m } + 1) is proportional to exp(J _{ΔN }). According to our estimates, the formula for J _{ΔN } was derived using the third, fourth, and fifth TID extrema values given in Table 1. Albee & Kanellakos (1968) assumed that an AGW is a transverse wave and therefore oscillations in N are accomplished only along a local vertical. Since such extrema generate LS AGWs, we can ascertain that this approximation is valid.
At distance R from a source the time interval t _{ i,i+1} = t _{ i+1} − t _{ i } between adjacent extrema in the AGW train, i.e., the temporal halfperiod τ _{ i }, is calculated using a linear dependence:(2)
The sum of τ _{ i } and τ _{ i+1} is the temporal quasiperiod T _{ i }. The curves τ _{ i }(R) for the first 14 quasiperiods are given in Figure 1, panel 1. The reference point for the index i coincides with the first extremum. The curves τ _{ i }(R) for i from 11 to 14, presenting the longest periods, differ from each other essentially, since they are created by a lowvelocity part of an AGW train. The first 10 halfperiods τ _{ i } differ from each other less than the other halfperiods. The dependence of τ _{ i } on the extrema index i is not always a monotonic function (see Fig. 1, panel 2).
Fig. 1.
Temporal halfperiods τ _{ i } versus distance R to the source, predicted by the model: panel 1 for i = 1–14, panel 2 for i = 1–10. 
Excitation of an AGW train starts by generating a minimum whose amplitude is A _{0} = 0.48A _{3} and velocity exceeds 11 km s^{−1}. Since the first extremum has a small amplitude and large propagation velocity, it may be difficult to detect it experimentally. An initial disturbance onset results in an air pressure decrease, which, in some cases, especially during a solar activity maximum, may be masked by subsequent extrema having larger amplitudes (up to ~2 times). Therefore, in our model the zerotime reference is related to the first extremum (maximum) whose amplitude A _{1} = 0.83A _{3} and velocity V _{1} = 3 km s^{−1}, although the disturbance onset occurs 3.95 min before this event. When the second extremum (minimum) is initiated, the first AGW spatial halfperiod moves away from the source area. When the third extremum (second maximum) is originated, the wave train represents the full quasiperiod of the horizontal wave. After origination of the fourth and fifth extrema the wave train contains the two full spatial quasiperiods. Formation of each two subsequent extrema increases the number of quasiperiods in a wave train by unity.
Let a horizontal distance between adjacent AGW train extrema, i.e., a horizontal spatial halfperiod, be l _{ i }. Then the sum of adjacent halfperiods l _{ i } and l _{ i+1 } is a spatial quasiperiod L _{ i } . According to our assumptions, l _{ i } is an increasing function of time:(3)where t′ = t − iτ _{0}. The halfperiod l _{ i } arises at time t′ = 0. Figure 2 shows l _{ i }(t) for the first 14 quasiperiods, calculated using equation (3) and data presented in Table 1. The minimum lengths of l _{ i } and L _{ i } are observed at the moment of their initiation. The data for the first fourteen halfperiods and seven quasiperiods are given in Table 2. For example, the fifth halfperiod, whose length is 79.1 km, originates at the 40th minute after disturbance offset. In time, the lengths of halfperiod l _{ i } and quasiperiod L _{ i } increase linearly from the lowest values observed at origination.
Fig. 2.
Spatial halfperiods l _{ i } versus time t, predicted by the model. Zerotime reference corresponds to the time of disturbance initiation. 
Minimal values of l _{ i } and L _{ i }.
Further, we assume a TID train consisting of the first fourteen halfperiods. When an AGW train moves away from a disturbance source, its halfperiods l _{ i } and quasiperiods L _{ i } are linearly increasing. Let R _{1} be a distance between a disturbance source and the first extremum in a TID train. The l _{ i } and L _{ i } for the first extrema at distances R _{1} = 1500 and 3000 km are given in Table 3 and Table 4. These values of R _{1} are observed 8.3 and 16.7 min after generation of the first extremum. The selected values of R _{1} are typical for the observations (see Sects. 4 and 5).
l _{ i } and L _{ i } for R _{1} = 1500 km.
l _{ i } and L _{ i } for R _{1} = 3000 km.
Apparently, at distances 1500–3000 km the first period of the wave train presents an LS TID and the next periods are MS TIDs. The smallest value of L _{ i } for the first extremum at such distances is more than 300 km.
2.3. The layer height and thickness
LS AWGTIDs are observed in the Earthionosphere waveguide, whose effective thickness is 2H = 160 km and a maximum TID amplitude is observed at the height z _{ m } of the F2 layer maximum. MS AGWs propagate in freewave mode and after reflection from the Earth’s surface at different angles and generate at ionospheric heights MS TIDs which are observed in a layer with effective thickness 2H = 50 km. The incidence and reflection angles are the greater the farther the reflection point from the disturbance source is.
2.4. Spatial scans of halfperiods
For a given value of R _{1,} the distance R _{ i } between the disturbance source and the nearest end of halfperiod l _{ i } in an AGW train is given by the expression:(4)
The spatiotemporal dynamics of halfperiods l _{ i } for different times t _{1} are shown in Figure 3, where positions of halfperiods R _{ i } and the first extremum R _{1} are marked by markers on rectilinear segments, parallel to the axis R. This plot is known as the spatial scan of halfperiods l _{ i }(R).
Fig. 3.
Spacetime dynamics of AGW train, consisting of 14 halfperiods. 
Panel 1 of this figure clearly shows the strong difference between the lengths of the first four halfperiods of the TID train, which are generated by the first five highspeed extrema (recall that V _{0} is not included in the calculations). The halfperiods l _{ i }, generated by lowspeed extrema of the TID train, are shown in Figure 3, panel 2. All fourteen halfperiods are formed at time t _{1} = 0.93 h, and only for t _{1} ≥ 0.97 h the left boundary of the fourteenth halfperiod and the whole TID train starts to move away from the excitation source. Thus, the differences between adjacent l _{ i } values in the lowspeed TID train part are much smaller than those in the highspeed part.
Figure 4 shows the plots L _{ i }(R) for distances R _{1} = 500–8000 km. Here R is the distance from the source to the middle of each quasiperiods. As can be seen L _{ i }(R) is an almost linear function in the range R = 1000–1500 km for the train part nearest to the disturbance source (corresponding to MS AGWs).
Fig. 4.
Spatial halfperiods L versus distances to a disturbance source R for various values of R _{1}. 
To analyze spacetime dynamics of AGW quasiperiods we performed a thought experiment. Let there exist a hypothetical observer capable of measuring visible AGW quasiperiods at predefined moments. The observer sequentially settles at measuring stations located at distances R _{ S } from the source. As soon as a left boundary of a nearest to the source quasiperiod crosses the vertical line at the observer’s site, the observer starts measuring the next quasiperiod, which is excited nearer to the source than the previous quasiperiod. To allow the observer to perform the measurements, we calculated the L _{ i }(R) scans, i.e., a set of quasiperiods for the set of t _{1}.
Figure 5 shows the graphs L(t) plotted with increments ΔR _{ S } = 500 km for the range R _{ S } = 500–3500 km where according to estimates the cusp can be located in our measurements (see Sect. 5). At first glance, sawtooth pieces on the curve L(t) most notably observed until a minimum value of L followed by almost monotonic increase indicate model inadequacy. In fact, such behavior should be observed at nearly continuous measurements of AGW spatial periods over a measuring station for different distances from a disturbance source.
Fig. 5.
Quasiperiod L versus time t for various distances R _{ S } between the measuring station and TID sources. Zerotime reference corresponds to disturbance initiation time. 
To explain the results obtained by the observer at various distances R _{ S } we will analyze the L(t) dynamics at R _{ S } = 500 km (Fig. 5, panel 1). The first value of L = 1482 km corresponds to a time delay t = 0.138 h (8.3 min) relative to the time of the first extremum excitation. At that moment the first extremum is at distance of 1500 km from a source, and the first quasiperiod consists of the first and second halfperiods. According to the measurements of our hypothetical observer the first quasiperiod as well as the next quasiperiods increases linearly as they propagate over the measuring station. This phase of the AGW movement over the measuring station corresponds to the first, linearly increasing part of the plot (up to ~3027 km at t = 0.324 h). As soon as the initial part of the second quasiperiod, i.e., the right boundary, crosses the observer’s vertical, the measured period will consist of the second and third halfperiods, either of the two being much shorter than the first quasiperiod. As a result this period, consisting of the second and third halfperiods, is considerably shorter. In our example its length is about 973 km. This new period is also linearly increasing up to ~1200 km, and it changes in a stepwise manner to a considerably smaller quasiperiod which consists of the third and the fourth halfperiods and is about 611 km long.
Thus, while the AGW is moving over the measuring station, the observed quasiperiod and its temporal gradient in each sawtoothed fragment of the plot L(t) are becoming smaller and smaller, since differences between adjacent values of V _{ i } in the train are reduced. This process continues until a minimum quasiperiod L _{ min } = 200 km is reached. Then almost monotonic increase L(t) is observed, since the differences between the adjacent extrema with indexes i > 8 are practically constant. The experimental observation of the sawtooth effect is illustrated in Section 5.
2.5. Slope of the AGW front in a vertical plane
The inclination of the TID front, defined by an angle ξ, measured in a vertical plane, is determined by:(5) (6)where R _{ E } is the Earth’s radius, z _{0} is the source altitude over the Earth surface, and the distance R is determined at the F2layer maximum height. The plots ξ(R) for both propagation modes are shown in Figure 6. Calculations of ξ(R) for the freewave mode were made using following parameters: the TID height is z _{ m } = 300 km and disturbance source heights are z _{0} = 100 and 200 km (curves 1 and 2, respectively). The distance R refers to a height z _{ m } = 300 km.
Fig. 6.
Inclination angles ξ versus R: 1 – MS TID, z _{0} = 100 km; 2 – MS TID, z _{0} = 200 km; 3 – LS TID. 
2.6. Spatial distribution of electron density
If a TID propagates in a direction defined by an azimuthal angle α, a spatial modulation of electron density N in a horizontal plane is described by:(7) where φ, λ are geographic latitude and longitude of the measuring station, respectively; ξ _{50} is the inclination angle at φ = φ _{0}; A _{ d } = A _{ d50} − (φ − φ _{0}) k _{ Ad }, N _{ n } is the undisturbed spatial distribution of electron density at a given moment, defined by a global ionospheric model, A _{ d50} is the relative amplitude of a TID at φ = φ _{0}, Φ is the initial phase shift of the corresponding disturbance harmonic, k _{ Ad } and k _{ ξ } are proportionality coefficients, and z is the height over the Earth’s surface. Plus or minus signs in the formulas for Δφ and Δλ are selected depending upon the relative position of the TID source and m. If the source is to the north or west of the measuring station, minus signs are used. Otherwise, plus signs are used. The functions F(z) and F _{1}(φ, λ) describe modifications of TID amplitude and phase depending on z, φ and λ. The global ionospheric model may be improved using navigation satellite signals (Tyrnov et al. 2005).
3. The choice of the model parameters
The model is based on the following postulates:

extrema periodicity;

spacetime constancy of V _{ i } and A _{ i };

invariability of τ _{0} during excitation of 16 train extrema.
The key model parameters are:

extrema velocities V _{ i };

relative amplitudes of train extrema A _{ i };

halfperiod of atmospheric pressure oscillations τ _{0} in an excitation source;
Due to the high velocities of the first extrema (see Table 1) accurate experimental values can be obtained only by continuous monitoring of ionospheric parameters, in particular, the total electron content or N _{ m }. Such approach was implemented by Roberts et al. (1982a, 1982b). In the first article the total electron content was monitored by measuring the Faraday rotation angle of linearly polarized VHF signals radiated by geostationary satellites. The experiment detected a disturbance initiated by a volcanic explosion (eruption of the volcano St. Helens, on May 16, 1980). These data were used to adapt the model to an isothermal atmosphere disturbance by a point source (Cole & Greifinger 1969; Liu & Yeh 1971). As a result, parameters of a simplified global equivalent isothermal atmosphere have been determined, namely: the speed of sound C = 486 m/s, adiabatic index γ = 1.38; the freefall acceleration due to gravity averaged over height range from 0 to 350 km g = 9.33 m/s^{2}.
The aforementioned model, equivalent isothermal atmosphere adapted, as it is named by Roberts et al. (1982a) was successfully tested by Roberts et al. (1982b), using experimental data obtained after the explosion of the volcano St. Helens. For the diagnosis of disturbances the Faraday effect of pulsar PSR 0329 +54 signals, recorded by a highly sensitive radio interferometer, was used.
Using data presented in the latter article, we calculated V _{ i } values for i from 1 to 11. The V _{ i } values for i from 12 to 15 are derived as linear extrapolation of V _{ i } values for i from 9 to 11. V _{0} was found via an extrapolation of ninthorder regression for V _{ i } values for i from 1 to 11. In what follows all V _{ i } values for i from 0 to 15 will be called the V _{ i } model (the label Model in Figure 7). The data thus obtained are given in Table 1. The spacetime constancy of V _{ i } values was experimentally confirmed by Breitling & Kupferman (1967), Row (1967), Kanellakos (1967), Albee & Kanellakos (1968), Lomax & Nielson (1968), Roberts et al. (1982a, 1982b), Igarashi et al. (1994) for distances from the disturbance source up to ~14 000 km. A _{ i } values are taken from Bristow et al. (1994), calculated for a delay relative to the onset of AGW excitation equal to 3.333 h. This delay is close to the average value used in our experimental data.
Fig. 7.
(a) Extrema velocities of a TID train versus its index for surface nuclear explosions (panels 1, 2, and 3) and volcano eruption (panel 4). (b) Extrema velocities of TID train versus its index for highaltitude nuclear explosions. 
We postulated invariability of τ _{0} during generation of all 16 AGW extrema. Following Albee & Kanellakos (1968), we have used the formula τ _{0}= 2πC/(γ g) (Row 1967). We substitute the parameters from an equivalent isothermal atmosphere adopted model (Roberts et al. 1982a) and obtained τ _{0} = 3.95 min.
The slope of the AGWTID front in the vertical plane for the MS TID (5) and LS TID (6) was obtained from Fedorenko et al. (2010b). The relation (5) was obtained for propagation of omnidirectional waves in a spherically symmetric medium without refraction. The relation (6), valid for the waveguide mode of AGW propagation, was obtained from Kunitsyn et al. (2007).
The values of 2H for the MS and LS TID, the height of the waveguide longitudinal axis for the LS TID, and the layer in which the MS TID are observed, were inferred from transionospheric sounding using signals of Russian loworbital navigation satellites (Cicada and CicadaM), and from incoherent scatter radar data (Fedorenko et al. 2008, 2010a, 2010b, 2011, 2012).
4. Model validation using literature data
4.1. Velocities of extrema
An estimate of the onset velocity V _{0} ≈ 10 km s^{−1} (in the model V _{0} ≈ 11.3 km s^{−1}) is given only in one single publication (Taran et al. 1985). The V _{0} value was estimated from incoherent scatter radar observations at a distance of 3500 km from the source of a surface explosion of yield Q ≈ 0.24 Mt TNT, detonated in the USSR on 28/11/1981 at 02:31 UT. Since the onset speed is difficult to measure experimentally, we could observe only the consequences of the initial disturbance passage. Therefore we could determine the extremum velocities V _{ i }, starting from V _{1}.
In Figure 7a, panel 1, the values V _{ i } for i from 1 to 15 given in Table 1 (solid line) are compared with our calculations based on the experimental data presented by Adushkin & Gorelyi (1994) (dashed line). In this paper the averaged data, obtained after 27 nuclear explosions with equivalent yields Q from 5 to10 and from 20 to 30 Mt TNT, are given. In the figure, only data for the second range of Q are presented. All the explosions were made at Novaya Zemlya during mass testing of nuclear weapons in the USSR. Explosion heights were unknown. However, we suppose the heights did not exceed ~5 km. Averaging was carried out by the method of epoch superposition. The critical frequencies of the F2layer (f _{0}F2) were recorded using ionosondes at two measuring stations (R _{ S } = 1000 and R _{ S } = 4000 km) in increments of 15 min, starting from the moment of the explosion. The arrival times of the first three extrema for R _{ S } = 1000 km in accordance with our V _{ i } model should be 5.6, 14.4, and 24.4 min. Therefore, we believe that these extrema could not be detected by this station. The V _{ i } values for extrema from 1 to 3 were determined using data obtained by the second station, and for extrema from 10 to 15 using data collected by the first station. The values for extrema from 4 to 9 were obtained by averaging data from both stations.
Data for the first five, the fastest, extrema, have the greatest measurement errors. The confidence intervals for the means were obtained using 15min intervals of f _{0}F2 measurements. The confidence intervals for extrema from 1 to 3 were calculated from the data of the second station, and for all others extrema from the first station. The ionosonde measurement time is about 1 min, the data were read from a small graph which was obtained by data averaging over an unknown number of explosions. As noted, data from 27 explosions were analyzed, but how many of them were in the range from 20 to 30 Mt are not specified in the article. We believe that the confidence intervals calculated for the 15min sampling intervals are only a part of the uncertainty interval. The error of the experimental data relative to the model values is 25.8%. On the whole, these data confirm the validity of the V _{ i } model.
On panel 2, the first 11 extrema of the model are compared with those obtained after a nuclear explosion (Q = 58 Mt TNT) over Novaya Zemlya, detonated by the Soviet Union, at 30/10/1961 08:33 UT at altitude z _{ expl } = 3.7 km (Row 1967). In this experiment, the critical frequencies f _{0}F2 were measured by ionosondes at several stations in 15 min intervals. The stations were located at ranges 1250–4400 km from the explosions. Despite the enormous nuclear yield, the 15min sampling rate allowed to record only the extrema four to eleven. The minima, i.e., extrema with numbers 6, 8, and 10, were determined by interpolation using a ninth order regression. Omission of the first three extrema can be easily explained by the following considerations. In accordance with the V _{ i } model, the first extremum should be detected by the ionosondes closest to the source Kiruna (R _{ S } = 1250 km) and Lycksele (R _{ S } = 1650 km), only in 6.9 min and 9.2 min, respectively, but it could not be recorded since the first measurement was made only 12 min after the explosion. The next two extrema (the duration of their existence over the measuring station was probably less than a minute) have also been missed owing to the large sampling period. Therefore, the first registered disturbance is actually a fourth extremum (minimum). Taking into account these limitations and the measurement errors we may conclude that the model and experimental curves are in good agreement, and the average error is 13.9%. Without such correction the difference between the experimental data and the V _{ i } model could not be explained by any model of Earth’s atmosphere. The confidence intervals were calculated for the Uppsala station (R _{ S } = 2123 km) which is about the average distance from the explosion center to measuring stations.
Panel 3 shows the data collected by a network of groundbased ionosondes (located in Pacific Islands) after a megaton lowaltitude US nuclear explosion in Johnston Island area on 30/10/1962 at 16:02 UT (Albee & Kanellakos 1968). Ground stations were located at R _{ S } = 1000–6500 km.
The first disturbance recorded in this experiment corresponds to a transition through zero of the AGW train and before the second maximum (third extremum), predicted by the model. According to our analysis, missing the first two extrema in this experiment cannot be explained by the measurement resolution (equal to 5 min), but by orientation of the wave vector of the first AGW halfperiod relative to the geomagnetic field lines at the measuring stations. The first halfperiod is classified as an LS AGW whose wave vector is parallel to the Earth’s surface. Magnetic inclination angles for all measuring stations were in the range 30–40°, and vertical oscillations in transverse AGW are practically orthogonal to the magnetic field lines. Only in subsequent train extrema, propagating in freewave mode, the wave vector after reflection from the Earth’s surface is oriented so as to allow plasma oscillations along the lines of force. Taking these considerations into account the experimental data error is 5.3%.
The first TID maximum following the volcano Pinatubo explosive eruption, recorded 15/06/1991 near Japan (R _{ S } = 3000 km) (Igarashi et al. 1994), corresponds to the seventh extremum predicted by the model. Since the velocities of this and subsequent extrema are small (see Table 1), the measurement in 15min intervals of the apparent reflection height of the F2layer provided data for comparison with the model which was proven to be in error by 4.4% (see panel 4 in Fig. 7a).
Comparison of data predicted by the V _{ i } model and experimental data collected after five highaltitude nuclear explosions over Johnston Island in 1962 (z _{ expl } varied from several dozens to 400 km) is presented in Figure 7b. The experiment involved more than 60 ionospheric stations (R _{ S } was in the range 1000–14 000 km).
In panels 1–3 the data from the studies of Lomax & Nielson (1968) and Breitling & Kupferman (1967) are combined. These data were collected for the same nuclear explosions at ranges of ~1000–5000 km and 4000–14 000 km with different temporal discretization equal to 1 min in the first work and 1 h in the second. Confidence intervals for data were estimated at R _{ S } = 2000 km and R _{ S } = 8000 km, which correspond to the mean distances for each measuring station involved in the two studies. The considerable spread of the confidence intervals can be explained by the different time sampling for the data.
Panel 1 shows the results obtained after an explosion with Q = 1.4 Mt TNT on 09/07/1962 at 09:00 UT at the height z _{ expl } = 399 km. The data errors for the minute and hour discretization are 34% and 79%, respectively. The maximum R _{ S } was approximately 14 000 km. The results from a nuclear test (Q = 0.02 Mt TNT) performed 20/10/1962 at 08:30 UT at the height z _{ expl } = 147 km are shown in panel 2. Here, the data errors are 14% and 56%, respectively. The dependence of the extremum velocity on its number in the TID train for a third nuclear test (Q = 0.41 Mt TNT) on 01/11/1962 at 12:10 UT at z _{ expl } = 97.5 km is illustrated in panel 3. The data errors are 25% and 34%, respectively. In the last two experiments, the maximum value of R _{ S } was ~12 500 km.
Panels 4 and 5 show the results calculated using the data in Breitling & Kupferman (1967) for nuclear explosions performed 26/10/1962 at 10:00 UT and 04/11/1962 at 07:30 UT at the height of a few tens of kilometers. The nuclear yields were ≈0.2–1 Mt TNT and 0.02 Mt TNT, respectively. The data errors of the experimental results for these explosions are 28% and 20%, respectively. The confidence intervals were calculated for R _{ S }=8000 km, and the maximal value of R _{ S } was about 12 500 km.
The experimental data errors relative to the data predicted by the V _{ i } model for nearsurface and highaltitude nuclear explosions are approximately the same. Errors assigned to the data from Breitling & Kupferman (1967) are not taken into account in this comparison because to their large errors due to the 1h discretization. Our conclusions were obtained using only data with approximately the same measurement error (measurement interval equal to 15 min and less). In general, we claim that the results of the above experiments confirm the V _{ i } model adequately.
4.2. The dynamics of horizontal spatial halfperiods
A spatial train position, consisting of four halfperiods defined by zero values of the disturbance, was given in the work by Albee & Kanellakos (1968). Observational data were obtained after a US megaton nuclear detonation near the Earth’s surface on 30/10/1962. We used these data to test the model validity. Figure 8 shows the spatial position of the observed (1) and modeled (2) trains of spatial halfperiods. As noted above, the first two extrema were missed in these nuclear tests, and the measurements start at zero perturbation (the second in the train), located in front of the third extremum (second maximum). Thus, the experimental train starts from the second halfperiod.
Fig. 8.
Scans of AGW train halfperiods. 1 – experimental values (Albee & Kanellakos 1968), 2 – values predicted by the model. 
Table 1 shows propagation velocities V _{ i }(i) of the first 15 extrema (i is extremum index). We obtained the velocity of AGW train fragments located between each pair of adjacent extrema and the corresponding unperturbed values of air pressure fluctuations in the train using tenth regression of V _{ i }(i). Table 5 shows the values of phase velocities of the train fragment V _{0i } (i _{0} is the index of undisturbed area in the train). Calculation of spatial halfperiods was carried out using relation (3), where the velocities V _{ i } were substituted by V _{0i } (see Table 5).
Velocities V _{0i } versus extremum index i.
The second and third halfperiods predicted by the model coincide practically with corresponding experimental data as far as the values and positions relative to the disturbance source are concerned. The halfperiod indexing goes from right to left and is performed using the model halfperiod marked by the markers from 1 to 7 in upper part of Figure 8. Experimental extrema observed by Albee & Kanellakos (1968) are connected by quasivertical dashed lines with extrema predicted by our model. According to the figures with experimental data ΔN/N(t) given in the article by Albee & Kanellakos (1968), only the first three extrema, corresponding to the third, fourth, and fifth extrema predicted by the model, are clearly marked, since their amplitudes A _{ d } reach 50–100%, while A _{ d } for the subsequent extrema are ~15 times smaller and are close to natural background fluctuations.
Therefore, only the second and third experimental halfperiods are indubitable. We suppose that the last two experimental halfperiods shown in the Figure 8 are really combination of the fourth and fifth, sixth and seventh quasiperiods. This is evidenced by the amplitude and spatial location of these TID fragments relative to the source. The last experimental halfperiod could also be the sum of the 6th, 7^{th}, and 8th halfperiods, since the 8th model halfperiod is very small and is compared with marker dimensions.
These assumptions are confirmed in Figure 9, which presents the experimental (curve 1) and model (curve 2) dependences of the halfperiod lengths upon the distance to the disturbance center. Experimental and model values for the second and third halfperiods are practically the same, while those for the fourth and fifth halfperiods are not as close as the others, due to the growth of experimental errors. The model error for the second and third halfperiods is about 6.5%.
Fig. 9.
Lengths of halfperiods versus distance between their centers and a source. 1 – experimental values (Albee & Kanellakos 1968); 2 – values predicted by the model. 
4.3. The dependence of the TID temporal period upon the distance from the source
The relation between the first measured temporal period and the distance to the explosion, carried out by the US 10/30/1962, can be found in Albee & Kanellakos (1968) (see curve 2 in our Fig. 10). As stated above, the first and second extrema were missed in this experiment, and the first recorded period corresponds to the sum of the third and fourth halfperiods. The starting point of the first halfperiod is the transition through zero ΔN _{ m } between the second and third model extrema. Therefore, we compared the experimental curve with the model curve, which is the sum of the third and fourth halfperiods (curve 1). The average error of the model curve does not exceed 3.6%. For the same nuclear test we compare in Figure 11 the model (curve 1) and experimental (curves 2) spatial dependencies for the third (panel a) and fourth (panel b) halfperiods. The experimental curves are the firstorder regression of the experimental values shown in the graph by square markers. The average errors of the model dependence relative to the regression of experimental values for the third and fourth halfperiods are 7% and 12%, respectively.
Fig. 10.
Temporal periods versus a distance R to the source. 1 – values predicted by the model; 2 – experimental values (Albee & Kanellakos 1968). 
Fig. 11.
The third (panel a) and fourth (panel b) temporal halfperiods versus a distance R to the source. 1 – values predicted by the model; 2 – experimental values (Albee & Kanellakos 1968). 
In the works by Row (1967) and Roberts et al. (1982a) the experimental data obtained from a nuclear test in Novaya Zemlya 30/10/1961 and from the St. Helens volcano explosion 18/05/1980 are presented. In these studies, they reported that 10 to 11 extrema were recorded in the TID trains for different distances from the disturbance sources. In Figure 12 we compare the spatial dependence of the average time period revealed in these experiments and two times the model halfperiods averaged over the first ten halfperiods (curve 1). Curve 2 presents the dynamics of TID quasiperiods registered after the nuclear explosion performed by the USSR at Novaya Zemlya (Row 1967).
Fig. 12.
Temporal period T versus distance R to the source 1 – values predicted by the model; 2 – Row (1967); 3 – Roberts et al. (1982a). 
Curve 3 is the average temporal period in the AGW train, obtained by Roberts et al. (1982a) after the St. Helens volcano explosion. The average error of the model (curve 1) with respect to the experimental curves (curves 2 and 3) is 10%.
4.4. Electron density decrease observed in connection with the propagation of a TID train
Lomax & Nielson (1968), analyzing the effects of three highaltitude nuclear tests, discovered an unexplained effect of marked decrease in electron density along the propagation path of a TID train within an area of 2500–3000 km radius. The effect was detected using vertical and oblique sounding with 1min resolution, which significantly improved the data quality, as compared to earlier experiments conducted with sampling periods of 1 h, 15 min, and rarely 5 min. The measurement of total electron content by satellite beacons, made 12 and 15 h after the explosion, revealed a TEC reduction up to 30% as compared to TEC measured the day before. The nuclear tests were performed at night on 09/07/1962 (Q = 1.4 Mt TNT, z _{ expl } = 399 km), 20/10/1962 (Q = 0.02 Mt TNT, z _{ expl } = 147 km), and 01/11/1962 (Q = 0.41 Mt TNT, z _{ expl } = 97.5 km).
According to our AGWTID model the propagation velocity of the zero extremum (minimum) in the train is more than 11.3 km s^{−1}, its amplitude A _{0} = 0.48. Therefore, at a distance of 3,390 km from the explosion center the atmospheric pressure will be much lower than the natural background value 5 min after the explosion. Since the lifetime of disturbances in the ionospheric Fregion can be more than 10 h (Tyrnov 2005), a long living area of reduced electron density values with a spatial extent of about 7000 km is formed. As noted above, the second extremum with amplitude A _{ i } = 0.83 (first peak), propagating at 3 km s^{−1}, was not recorded in these experiments since air oscillations in AGW were orthogonal to the magnetic field lines. For the St. Helens volcano eruption (Roberts et al. 1982a), the registration conditions were more favorable in this respect, and the first extremum as well as the onset and subsequent extrema were registered.
We now return to the results of highaltitude nuclear tests. The extrema of AGW trains having smaller velocities and propagating in an area with reduced density develop the next slower TID train extrema. Thus, at distances ranging from 1000 to 5000 km the measurement results are a superposition of longlived disturbance phases and oscillations with temporal periods from 30 to 120 min (see also Fig. 12). We believe that the abovementioned unexplained effect indicates that the onset has already passed and gave rise to a region of lower concentration (also called a hole), which may be long living due to low solar activity. Such regions are formed both during nuclear tests (e.g. Lomax & Nielson 1968) and energetic particles precipitation in the magnetospheric cusp after sudden commencement of magnetic storm (Shashunkina 1966, 1968a, 1968b, 1972).
A similar signature of a low electron density region was observed after a nuclear explosion at Novaya Zemlya on 30/10/1961 (Q = 58 Mt TNT, z _{ expl } = 3.7 km). Adushkin & Gorelyi (1994) give results in graphical form of the F2layer critical frequency (f _{0}F2) measurement, made in the period from 27/09/1961 to 05/10/1961 at the following ionospheric stations: midlatitude AlmaAta (USSR, R _{ S } = 3500 km), southern hemisphere Mundaring (Australia, R _{ S } = 12 000 km), and highlatitude Mirny (Antarctic, R _{ S } = 16 000 km). The measurements began 100 h before the explosion and continued for up to 200 h after the explosion. The critical frequencies f _{0}F2 above AlmaAta were reduced for 44 h, and then began to increase, while above the Southern Hemisphere sites they were enhanced for about 50 h at Mundaring and about 59 h at Mirny and then began to decrease. The maximal amplitude (about 50%) of this longperiod wave was observed over Mundaring. Data presented in the paper do not allow us to calculate the wave velocity due to small plot scale, but we can assert that the concentration began decreasing over AlmaAta 17 h after explosions, and began increasing after 21 h over Mundaring and after 35 h over Mirny. The experimental results presented in that paper indicate that if an AGW propagates from north to south in the northern hemisphere, AGW air pressure minima generate TID minima, and in the southern hemisphere, AWG minima generate TID maxima. Moreover, the duration of TID halfperiods is increasing while the TID moves away from the explosion epicenter. These two effects confirm that our model is valid at distances up to 16 000 km from the disturbance source.
4.5. Global observation of AGW and TID
Disturbances generated by the nuclear tests on Novaya Zemlya on 23/10/1961 (Q = 25 Mt TNT) and 30/10/1961 (Q = 58 Mt TNT, z _{ expl } = 3.7 km) are the most studied globally, since atmospheric pressure waves propagated across the globe. For the second event pressure disturbances were observed by all ground stations. Isolines of disturbance delay time in both hemispheres were similar to concentric circles centered at the explosion epicenter (Wexler & Hass 1962). Using highquality microbarograms recorded in Washington at R _{ S } = 7000 km (Wexler & Hass 1962), we estimated that the velocities of pressure waves at the Earth’s surface for these two events were equal to 0.33 and 0.31 km s^{−1}. The velocity of the pressure wave generated by the second event, averaged over the globe as reported by Wexler & Hass (1962) and Obayashi (1962), was 0.31 km s^{−1}.
However, ionospheric disturbances were not observed by all stations. Ionospheric disturbances generated by the two events were observed in Europe, Asia, India, and Southeast Asia. Small disturbances were observed in Japan, and only very minor disturbances or no disturbances at all were observed in North America (Obayashi 1962), despite the fact that R _{ S } for stations in North America are comparable with R _{ S } for the European stations. These dead zones of the ionospheric disturbances confirm the validity of our model, which postulates that the ionospheric manifestation of atmospheric waves is determined by the angle between the wave vector of AGW and the geomagnetic meridian plane, and by the projection of the ionized component velocity on the geomagnetic field vector. Possible distortions of the AGW front, caused by regional wind velocities and by pressure and temperature heterogeneities, observed, for example, by Wexler & Hass (1962) should be also taken into account. As it turned out, the front curvature of the direct wave can reach 40°.
5. Model validation by ionosphere sounding using signals of loworbit satellites
5.1. The research methodology
Spatial distributions of the electron density were obtained from radio signals radiated by the navigational satellites Cicada and CicadaM (Russia) orbiting at an altitude of approximately 1000 km at an inclination of 83°. The onboard transmitter continuously radiates coherently at 150 and 400 MHz. A ground station located at the Radiophysical Observatory of V.N. Karazin Kharkiv National University (φ = 49.63° N, λ = 36.32° E) receives the signals and determines the ionospheric component of the Doppler shift δf _{ i } for 12–16 min every 1.5–2 h. Registrations of δf _{ i }(t) are indicators of the state of the ionosphere, allowing to investigate the regular and perturbed parameters of the spatial distribution of electron density for a region in the latitude range 33–66° N, corresponding to ~3600 km at the Earth’s surface. During a single satellite pass, the ionosphere parameters can be inferred over a latitude range of 24° (about 2700 km). Details of the method are described by Tyrnov et al. (2005) and Fedorenko et al. (2010b, 2012).
5.2. TID occurrence
We have analyzed the dependency of δf _{ i }(t), recorded in 1987–2010 under specific heliogeophysical conditions. The amplitude of the effects caused by the TIDs is ten to hundred times greater than the instrumental error δ = 0.005 Hz. The high measurement accuracy and huge amount of data, accumulated over two cycles of solar activity, allow us to disprove conclusions of other researchers, who conducted experiments with the loworbit navigation satellite NNS (Transit) and NNSS (USA) but with insufficient measurement accuracy. According to our observations, the generally accepted view concerning the presence of latitudinal, seasonal, and diurnal occurrence of MS TID is not true. For example, Evans et al. (1983) and Ogawa et al. (1988) state that MS TIDs are observed least often in the night, and in the summer there is a minimum of their occurrence. We have analyzed the data provided in the abovecited papers, and concluded that the results may be explained by the lack of the measurement precision. There exists a threshold value of minimum electron density N _{ m }, below which the measurements are impossible. It should be noted that both the abovecited papers make similar statements. Evans et al. (1983) claim that the poor measurement accuracy influences only the amplitude of detected dependencies and does not affect their validity: “... which we believe are attributable, in part, to the diurnal and seasonal variation in total electron content, together with a minimum size detectable disturbance”. As far as diurnal occurrence is concerned, Ogawa et al. (1988) state: “... which may result from easier detectability of TID’s during daytime hours”. Observations described in these studies were carried out in different hemispheres, so when Evans et al. (1983) have observed a maximum in the winter, Ogawa et al. (1988) have observed a summer minimum at the same time (and vice versa). According to our hypothesis, most TIDs of natural origin are generated by energetic particle precipitation in the magnetospheric cusp, therefore, the TID occurrence should not depend upon the seasons. Moreover, even a single act of AGW excitation generates TIDs observed almost for 24 h. Therefore, maxima and minima in the diurnal and seasonal dependence of the TID occurrence would be observed if the measurement accuracy is low.
Tsugawa et al. (2007a, 2007b); Kotake et al. (2007); Song et al. (2011); Otsuka et al. (2013) and Vlasov et al. (2011) after almost three decades of research performed by Evans et al. (1983) and Ogawa et al. (1988) have also studied the diurnal and seasonal dependence of TID occurrence. However, the dependences obtained by these authors contradict each other and previous results. For example, Tsugawa et al. (2007b) observed minimum occurrence in the summer 2006 while Vlasov et al. (2011) detected maxima in the summers 2007–2008. Vlasov deduced that a seasonal dependence of TID amplitude was absent and that MS TIDs are not observed in near midnight hours. These findings are consistent with our results. Differing conclusions regarding diurnal and seasonal occurrence may be explained by measurement errors of AGWTID parameters, which were determined by measuring the TEC, electron density, ion velocity, electron and ion temperatures. A possible cause of such dependence is explained by Vlasov et al. (2011) as follows: “The seasonal variation in AGWTID occurrence achieved in this study follows closely the seasonal variation of ionospheric electron content at the observing site. It remains still as a topic for future studies whether the observed seasonal variation is natural or caused by seasonal variation in the performance of our observation method”; “The AGW signatures are clearest in ion velocity data. One has to remember that if we do not see the signatures in the other ISR parameters it does not mean their total absence. For example, in electron density the AGW signatures can be masked by some other background processes”.
Continuous registration of δf _{ i }(t) allowed us to diagnose the ionosphere for all seasons on a 24h basis for different levels of solar and geomagnetic activity. We have determined that TID occurrence does not depend on the time of day. Table 6 shows the seasonal TID occurrence during the years of moderate and low solar activity. Following the paper of Evans et al. (1983), we classify measurements made from October 27 to February 4, from April 30 to August 2, from February 5 to April 29, and from August 3 to October 26 as winter, summer, and equinox measurements, respectively. Figure 13 illustrates the annual average of TID occurrence (circles) and its linear regression (curve 1), the yearly average of sunspot number (curve 2), and annual average of TID occurrence with linear trend removed (curve 3) for the period from 1987 to 2010.
Fig. 13.
TID occurrence and solar spot numbers in 1987–2010: 1 – linear regression of TID occurrence; 2 – W/100, W is yearly average Wolf number; 3 – TID occurrence without a linear trend, 4 – the occurrence of TIDs propagating in waveguide and freewave modes and their linear regressions for the period from 1997 to 2010. 
TID occurrence.
Based on the above analysis we conclude that during the 23year period the TID occurrence does not depend on the season and on sunspot number. The TID occurrence, defined as ratio of TID events to the total number of observations for the corresponding period, increased from ~50% in 1987 to ~98% in 2010, and the average was 87%. This conclusion may explain the scatter of occurrence values observed during the 60year period. At the same time, we support the conclusion of most researchers that MS TID occurrence does not depend on the magnetic activity (e.g., Evans et al. 1983; Ogawa et al. 1988; Hocke & Schlegel 1996; Tyrnov et al. 2005; Tsugama et al. 2007a). In the years of minimum solar activity, about 70% of LS TIDs are observed in the midnight time and about 30% in the morning.
In addition, we managed to explain latitudinal dependence of MS TID occurrence and reasons why MS TIDs are simultaneously detected both on the northern and southern branches of the δf _{ i }(t) registrations. It turns out that TIDs exist in the entire observation area and are detected at one or both δf _{ i }(t) branches, depending on the inclination angle of the TID front in a vertical plane (Fedorenko et al. 2010b). This effect will be described in detail below.
We use this effect to classify TIDs by their propagation mode, namely, the waveguide mode and freewave mode. According to our model, the LS TIDs propagate only in waveguide mode, while MS TIDs can propagate in both modes. If the spatial horizontal period of MS TID does not depend on the longitude, we assume that the TID propagates in the waveguide mode. If MS TIDs were observed only in the northern or southern part of the area where δf _{ i }(φ) was recorded or their amplitudes are dominated in one of these parts, we assume that a freewave mode was present. The distribution of the TID with respect to the propagation mode was carried out according to the prevailing influence of one or another TID type on the registration of δf _{ i }(t). Therefore, if an MS TID and an LS TID are simultaneously traced in a δf _{ i }(t) registration under the conditions that LS TID only slightly deforms registration obtained under undisturbed conditions, the TID is classified as propagating in freewave mode. We analyzed only the registrations in which such classification can be reliably made. As a result of this selection we retained for periods 1977–1999 and 2009–2010 only about 50% and 30%, respectively, of the available registrations. The analysis results are shown in Table 7.
Seasonal occurrence of TIDs, propagating in waveguide (WGM) and freewave mode (FWM).
As can be seen, the occurrence of TIDs propagating in a freewave mode (MS TIDs only) is from 3 to 18 times greater than that for TIDs propagating in a waveguide mode. We did not take into account the summer data for the period from 2009 to 2010 when the ratio WGM/FWG was equal to 0.1, for yet unexplained reasons, which are anomalous as far as the propagation mechanism of natural TID is concerned. Perhaps it may be explained by small statistics for the selected period. Overall, the results confirm the wellknown conclusion that LS TIDs are observed much less frequently than MS TID. We will show below that an LS TID may exist for several hours and MS TIDs of different scales can propagate against their background. These data are consistent with the dynamics of TID and its horizontal dimensions as described by our model, which predicts that an initial disturbance generates highspeed LS TID, followed by MS TID propagating for a much longer time period with a much lower velocity.
At the bottom panel of Figure 13 the occurrence of TIDs propagating in waveguide and freewave modes and their linear regressions for the period from 1997 to 2010 are shown. No dependence on sunspot number was detected. The regression lines are nearly horizontal, suggesting no epoch or solar cycle dependence. According to the data for this period the TID occurrence propagating in freewave and waveguide mode is 76% and 24%, respectively. In fact, the statistics reflects the observation time of MS and LS TID. Thus, MS TIDs are observed three times more often than LS TIDs.
5.3. The minimum values of TIDs horizontal scales
As stated above, L _{ min } may vary in the range 200–747 km for R _{ S } between 500 and 3500 km. Using cusp models, presented by Köhnlein & Raitt (1977) and Deminov et al. (1992), we calculated that in our experiments registrations of δf _{ i }(t) start at distances R ~ 500 km from the cusp. Therefore, the minimum value of L, calculated according to our model, should not exceed 200 km. With a decrease of R the value of L _{ min } is also reduced down to the minimum value of the selected train period, observed at its onset. We now estimate the value of L _{ min } for R ~ 500 km in a slightly different way. According to our model, during initiation the AGWTID spatial quasiperiods have the minimum values decreasing with their index in the TID train (see Table 2). We consider the quasiperiods which generate the MS AGWTIDs at the given R. For example, the initial value for the seventh quasiperiod L _{ min } = 47 km linearly increases as the distance from the source grows. At a distance R = 480 km the value is L = 245 km, and at R = 1022 km it is L = 412 km. The first extrema are now at R _{1} = 1200 km and 20 000 km, respectively. For the fourth quasiperiod, which was recorded much more frequently in our experiments, L _{ min } = 104 km at R = 512 km and R _{1} = 7000 km, L predicted by the model is equal to 438 km. At R = 1006 km and R _{1} = 11 000 km the model value is L = 689 km. Let us examine the panel 2 in Figure 4, which shows L(R) for the trains in which the first extremum was at distances R _{1} in the range from 4500 km to 8000 km, corresponding to an elapsed time of 25–44 min after the initiation of the disturbance. As can be seen from Figure 4, at a distance of R = 500 km the L values for R _{1} are in the range of 400 km to 600 km. Thus, if MS TIDs dominate the ionosphere above a measuring station, the TID with periods less than 100–150 km should not be detected in the δf _{ i }(t) registrations.
We now turn to our observational data obtained under specific heliogeophysical conditions. As a result of our data analysis, we found that the observed minimum values of L are greater than ~150 km. Similar conclusions were obtained by Francis (1975) and Fedorenko et al. (2010a, 2010b). The experimental verification of the lower bound for L ~ 200 km validates our model adequacy.
5.4. Precipitation of solar wind energetic particles into the magnetospheric cusp
We suppose that the source of most disturbances observed during the whole period of our experimental studies is associated with the precipitation of solar wind energetic particles into the magnetospheric cusp of the northern hemisphere.
This statement may be confirmed by the following observations. The horizontal length L, the amplitude A _{ d }, and the modulus of the inclination angle ξ are increasing with decreasing latitude φ. Since ξ is negative, the front TID in the vertical plane is tilted toward low φ i.e., in the direction of AWG propagation. Now we will give separate experimental results supporting our conclusion.
Figure 14 shows δf _{ i }(φ) and L(φ), obtained in winter, summer, and autumn, under average solar and geomagnetic activity (W = 52–97, Kp = 3 … 5+) around noon local time.
Fig. 14.
The results of MS TID observation for ξ = 42 – 75°. Parameters for the most southern irregularity are given in lower part of panels 1, 2, and 3; parameters for the rest of irregularities are given in upper part of these panels. 
The L values were determined by measuring the distance between adjacent minima or maxima δf _{ i }(φ), at various latitudes φ. Therefore, the number of recorded values is almost twice the number of TID quasiperiods appearing on the δf _{ i }(φ) registrations. The experiment date, the time of the satellite traverse (t _{0}), the satellite zenith angle at the traverse (γ _{0}), L, A _{ d }, ξ, and the unperturbed N _{ m } and z _{ m } over the measuring station at the time of measurement are shown in the plot. Hereinafter, we use the Standard Moscow Time (SMT), offset is UT+2.
The abovementioned parameters were obtained under the assumption that the TID wave vector is directed along a geographic meridian. The MS TIDs are clearly identified on the plots δf _{ i }(φ). It should be kept in mind that when a radio beam approaches the satellite traverse and the angle γ _{0} is in a range ~0°–20°, the TID amplitude inferred from the measured δf _{ i }(φ) is significantly reduced even if the TID amplitude A _{ d } stays unchanged in the entire region of radio sounding (Tyrnov et al. 2005; Fedorenko et al. 2010b). Therefore, near the traverse our methodology is less suitable for TID detection. In addition, for large values of ξ (>40°) TIDs are observed only at that part of the δf _{ i }(φ) registration chain (before or after the traverse), which is further from the source (Fedorenko et al. 2010b). This effect specifies the direction of TID movement (south or north) and the approximate location of the source (to the south or north relative to a measuring station).
The halfwidth H of the layer which guides the TID was 25 km. We have found that for most irregularities the parameters are in the ranges: A _{ d } = 0.07–0.17, L = 250 km, ξ = 42°–75°. For the most southern irregularities the corresponding values are: A _{ d } = 0.13–0.5, L = 415–600 km, ξ = 59–75°. We calculated the above parameters using the model of perturbed spatial distribution N (6). The parameters N, N _{ m }, and z _{ m } of the unperturbed model were determined by adapting the global model of N by the method described by Tyrnov et al. (2005). All δf _{ i }(φ) registrations indicate that TIDs always appear to the south of the latitude traverse, clearly indicating that the disturbance source was located northward of the measuring station.
The correlation coefficients of the linear regression L(φ) were in the range from 0.78 to 0.99. Figure 4 shows the approximately linear increase in spatial quasiperiods L of a TID train as predicted by our model when the wave packet moves away from the source. Here we present the model dependence of quasiperiods L upon the distance between their midpoints and the source for various distances R _{1}.
Figure 15 shows the data accumulated at midnight summer time during minimal and small solar and magnetic activity. In this experiment the MS TIDs were observed in a layer with halfthickness H = 50 km, though such values in our experiments were seldom observed. The TID amplitude for most irregularities registered in this satellite pass was A _{ d } = 0.08. If a TID wave propagated in the northsouth direction, the most intense spectral component of the irregularity period was L = 425 km. As in previous experiments we observed that the ξ was negative and the function L(φ) may be approximated by a linear regression (q = 0.63). It turned out that functions L(φ) are decreasing and TIDs were observed to the south of the measuring station. The above observations indicate that the excitation source of registered TIDs was located at high latitudes.
Fig. 15.
The results of MS TID observation for ξ = –(50 – 62)°. 
Table 8 shows the estimated values of distance R between a TID source and the measuring station calculated using formula (5). It was assumed that the energy of the precipitating particles is released at an altitude of 200 km (see also Fig. 6). As can be seen, the distance from the source was in the range 650–1900 km.
The distances between the MS TID and measuring station, shown in Figures 14 and 15.
To verify the above calculations we estimated distances to the disturbance source by another method. As noted above, the TID quasiperiods, defined by our model, are increasing as distance from the disturbance source increases. As mentioned above, the specific distribution of spatial quasiperiods in the TID train is observed for selected values of R _{1} (see Fig. 4). Therefore, if the experimentally observed distribution is approximated by the function L(R) predicted by the model for the chosen R _{1} values, the distances between the disturbance source and each quasiperiod can be specified. To illustrate the calculations we choose the three values of quasiperiods shown in Figure 14 (panel 2), since these data have the highest correlation coefficient of linear regression. The approximation for the function L(R) was made using a sampling increment for R _{1} equal to 500 km. The results are illustrated in Figure 16. The model error relative to the experimental data is equal to 5.8%. Figure 14 (panels 1 and 2) indicates that the quasiperiod L closest to the measuring station (latitude 49.63° N) has the smallest length. As can be seen from Figure 16, the distance to the middle of this quasiperiod is 750 km, which is close to the distance between the source and the measuring station, R = 650 km, estimated by ξ values, averaged over the satellite pass (see Table 8).
Fig. 16.
1 – the experimental curve R(L); 2 – the curve R(L), predicted by the model. 
5.5. Earthquake in South Sumatra
Figure 17 shows the data recorded September 30, 2009 and October 01, 2009 around midnight during a period of low magnetic and low solar activity. The local time of observation, the magnetic and solar activities were similar to those in the experiment whose results are shown in Figure 15. However, in contrast to those data and the data illustrated in Figure 14, the linear regression of L(φ) presented in Figure 17 is a decreasing function of latitude φ. The TIDs are observed in the northern part of the δf _{ i }(φ) registration chain. According to our model, this suggests that the source of these TIDs is located to the south of the southern boundary of radio visibility (35° N), and that the disturbance moves from the south to the north.
Fig. 17.
The plots δf _{ i }(φ) and L(φ), registered after South Sumatra Earthquake. 
We found that the same propagation direction was also observed in other registration cycles during magnetic storms. In these rare cases, the location of the magnetospheric cusp was located south of or above the measuring station. These conclusions were made by calculation of the cusp location using empirical models by Köhnlein & Raitt (1977) and Deminov et al. (1992), which defines the cusp position depending on the longitude and geomagnetic activity. The experimental data presented in Figure 17, correspond to low magnetic activity, therefore the magnetospheric cusp was excluded as AGW source.
We took into account the fact that the abovementioned data were registered after the earthquake in South Sumatra (Indonesia) at 30/09/2009 13:16:08.6 MST (10:16:08.6 UT). The earthquake magnitude was 7–7.5, the depth was 90 km, the epicenter coordinates were 0.84° S, 99.8° E. The azimuthal angle to the epicenter measured from north was ζ = 111.3° and the distance between the epicenter and the measuring station R _{ quake } = 8200 km.
The registrations of δf _{ i }(t), obtained after the earthquake and their approximation by the model under unperturbed and perturbed conditions, are presented in Figure 18 on the left and right panels, respectively. Adaptation of the global spacetime regular model of N was carried out using the method described by Tyrnov et al. (2005). The values of N _{ m } and z _{ m } shown in the figure were calculated from model values corrected for the measuring station coordinates. The TID parameters were estimated by an iterative technique, using the possible range of the parameters defined during the direct problem solution. We have not always achieved exact phase matching between the model and experimental curves, since our main objective was to estimate the spatial scale and amplitudes of the disturbances. To avoid cluttering, the tick marks at the time axis origin indicate the hours and minutes and further along the axis only minutes of a current hour. We found that when an MS AGW, formed by extrema numbers 7, 8, and 9 (inferred from the model), propagates from the epicenter to the measuring station, it transforms into an LS AGW and, in addition, the azimuthal angle of the wave vector changes from 111.3 to ~150°. We believe that such effect may be explained by regional winds on the propagation path. A less probable cause may be regional anomalies of pressure and temperature. As can be seen, in all three sounding intervals the MS TID with L = 2134 km has the largest amplitude A _{ d } = 0.6. All parameters of this TID: L, A _{ d }, H, ξ, and Ф did not change for at least 3.3 h. In the years close to solar activity maximum, the lifetime of the AWG track in the ionosphere can be more than 5 h. Therefore, the δf _{ i }(t) registrations shown on panels 2 and 3 indicate that in this experiment the wind velocity component along the AGW wave vector was negligible.
Fig. 18.
The recordings of δfi(t) obtained after South Sumatra Earthquake (solid lines) and their model approximations for undisturbed conditions (left panel) and for disturbed conditions (right panel). “SN” – satellite moves from the south to the north; “NS” – satellite moves from the north to the south. 
The ξ value for the TID of any scale is 80°. The positive sign and large value of ξ indicate that the disturbance source is located thousands of kilometers to the south of the measuring station. Above, we have identified the source that gave rise to the LS TID and now we consider an MS TID with L = 350 km and A _{ d } = 0.2 (see panel 1). The amplitude of this irregularity stayed constant for about 30 min while the period L was increasing to 500 km (panel 2). Then the period did not change within 2.8 h, but the amplitude increased to A _{ d } = 0.3. Undoubtedly, these MS TIDs are also generated by a southern source, but in our opinion, they cannot be related to the earthquake because lowspeed extrema, observed as an MS TID, should be observed later, but not simultaneously with an LS TID.
A series of earthquakes with magnitude from 6 to 7.5 occurred prior to the abovementioned earthquake and the shocks nearest in time were observed at 8:25:46 and 11:24:57 MST (05:25:46 and 8:24:57 UT), respectively. The first shock took place in Tonga Islands (15.56° S, 173.49° W, R _{ quake } = 15 368 km, ζ = 134), and the second in Samoa Island (16.72° S, 172.87° W, R _{ quake } = 15 513 km, ζ = 134°) at a depth of 10 km. These shocks occurred 2 and 5 h ahead of the AGW, generated by the South Sumatra earthquake, and their epicenter azimuths are approximately equal to that of the AGW source. A superposition of the extrema, generated by these earthquakes, could give rise to the MS TID in the sounding region. According to our estimates, the extrema of each of these earthquakes at delays equal to ~11.5 and 14.5 h could generate over the measuring station only LS TID (L = 1500 km). The distance between the earthquake epicenters was equal to ~15 000 km, and it does not contradict this hypothesis, since disturbances, arising from highaltitude nuclear explosions with an energy in the range from 0.02 to 1.4 Kt TNT, were observed at similar distances (see Sect. 4).
5.6. The time dynamics of L and A _{ d }
According to our TID model, initial values of velocity extrema V _{ i } differ very much. Thus, even adjacent initial spatial halfperiods in a TID train differ essentially in magnitude. This effect has always been observed experimentally in the first hours after the initial disturbance. Registrations of δf _{ i }(t) of MS TID are illustrated in Figure 19 where all identified disturbances were moving from the north to south, i.e., from the cusp to the equator. The direction of TID movement is unambiguously defined by the ξ sign. The excitation time of the TIDs, detectable in registrations presented in panels 1 and 2, was about 19:00 and 21:00, respectively. Our estimates show that location and energy of the TID sources were almost identical in both events, considering the very similar manifestation of disturbances (see panels 1 and 2) and on almost equal delays relative to the excitation time (about 2 h). The TID excitation time inferred from the registrations δf _{ i }(t) is a characteristic indicator for highenergy proton precipitation into the magnetospheric cusp.
Fig. 19.
Examples of MS TID with A _{ d } = 0.18 … 0.95 and ξ = –45 … –80°. The recordings of δfi(t) obtained after precipitation energetic protons near the magnetospheric cusp (solid lines) and their model approximations for undisturbed conditions (left panel) and for disturbed conditions (right panel). “SN” – satellite moves from the south to the north; “NS” – satellite moves from the north to the south. 
The dependencies noted above show that in the first hours after the AGW initiation L and A _{ d } are substantially increasing with increasing distance from the excitation source. Tsugawa et al. (2007a) observed the effect of TID amplitude increase but did not explain it. However, the effect is predicted by our model. The effect was also observed in our study. Thus, the data presented in the first and second panels show that the periods are increasing from 400 to 800 km and from 320 to 1000 km, respectively (increase from 2 to 3 times), during the second hour after excitation. The observations recorded 14 h after TID excitation and presented in the third panel show the period increasing from 400 to 650 km. This suggests that these disturbances are generated by subsequent slower quasiperiods of the AGW train for which, according to the model, the period increases much less. Increasing of L and A _{ d } for MS TID moving to the equator over a region with a typical extent of about 3000 km gives further evidence that the source of this TID is at northern latitudes. Using a large data set and another methodology, Fedorenko et al. (2010a) have ascertained that the L values are increasing if the latitude is decreasing.
Using δf _{ i }(t) registrations, we found that when there exist only the MS TIDs or MS TIDs and LS TID simultaneously, a general rule can be formulated as: “the larger is L, the greater is the relative amplitude A _{ d }”. For LS and MS disturbances the ratio between A _{ d } values may be up to 10. This rule is predicted by our model and is confirmed in most of our experiments. A similar correlation has been established by Fedorenko et al. (2010a) using statistical data analysis for MS TID.
Additional confirmations of our findings are shown in Figures 20a and 20b, where for low solar activity (W = 18–25) the diurnal variations of L, A _{ d }, N _{ m }, and z _{ m } above the measuring station for the period between 02/03/10 and 05/03/10 are shown. The diurnal data L(t) and A _{ d }(t) were obtained by δf _{ i }(t) registrations using the method proposed by Fedorenko et al. (2012), which allows to identify the TID having a spatial scale in the range 200–35 000 km and relative amplitudes from ~0.02 to 1 over an observation region with a radius of about 1500 km. This method may also be used to analyze MS TIDs and LS TID simultaneously observed over the same region.
Fig. 20a.
The diurnal dependencies of L, A _{ d }, z _{ m }, and N _{ m }, registered on 2 and 3 of March 2010. 
Fig. 20b.
The diurnal dependencies of L, A _{ d }, z _{ m } and N _{ m }, registered on 4 and 5 of March 2010. 
The plots L(t) and A _{ d }(t) for LS TID and MS TID are presented in Figures 20a and 20b, respectively. If LS and MS TIDs were observed in the southern and northern parts of the registration station chain and have different L values, then their average values are indicated in the figure. We have estimated that for the LS TID with L ≈ 1500–35 000 km the amplitudes A _{ d } are in a range 0.4–1, and for MS TID with L ≈ 150–800 km in the range A _{ d } ≈ 0.02–0.5. We believe that the LS TIDs, presented in Figures 20a and 20b, are generated near midnight.
The sawtooth curves of L(t) may be explained by large differences in adjacent spatial periods at the beginning of the TID train. Such dynamics are readily explained by our semiempirical model.
An important feature of the data obtained during minimum solar activity is that when a disturbance starts, an LS TID is generated, and a significant decrease in N, as compared to the unperturbed background values, is observed. Such volume of reduced density is known as a hole. According to our model, oscillations in a disturbance source start with lowering air pressure up to two times (see Table 1).
This area of reduced air pressure propagates as a wave in the horizontal plane in all directions at 11 km s^{−1}. In our experiment, an area of reduced atmospheric pressure was identified by monitoring N. The latitude of minimum values of N within such volume is moving in the range from ~60° to 40° N. A relative amplitude A _{ low } of holes in this latitude range is reduced between 1 and ~0.3 (A _{ low } = ΔN _{ low }/N _{ n }, where ΔN _{ low } is the absolute decrease of N _{ m } in a hole and N _{ n } is the unperturbed N _{ m }). These changes were recorded during the daily observation periods (lasting for about 15–20 h). Thus, during the day (in fact, from one event of particle participation to another) we observed a significant decrease in the background N levels during the minimum of solar activity (W = 18–25). As noted in Section 4 the significant decrease of N was also observed after nuclear experiments in 1962, at low solar activity (W = 11–28).
According to our observations, TIDs are generated frequently at high latitudes (in the magnetospheric cusp area) and can be observed in the midnight sector. The TIDs may be excited in the morning cusp and traverse the observation point after passing over the North Pole. Here we present only the observation results without analysis of the process in detail. Often the event (precipitation of energetic protons) occurs once every 2 or 3 days, but sometimes it may appear two or three times per day (see, Fig. 19). In the latter case, L and A _{ d } have no time to decrease to values normally observed at the end of the day with a single disturbance source (see, for example, Fig. 20a for 02/03/2010). Rather, a spatial AGW trace existing in the ionosphere may be superseded by a new trail having larger amplitudes and quasiperiods. If the time between consecutive precipitation events is more than a day, the values of L and A _{ d } are reduced to their minima reaching ~150 km for L and up to several percent for A _{ d }.
From the data shown in Figure 20, we may conclude that the values of L and A _{ d } decay very slowly over about 5 h. This effect may be explained by our model, which states that only the first two periods of an AGW train have maximum A _{ d } and subsequent periods are decreasing linearly with time. Due to differences between adjacent train velocities V _{ i }, decreasing from one to another halfperiod, the corresponding values of L are also reduced, at first very sharply and then only little. Consequently, the third, fourth, and subsequent halfperiods of an AGW train could not erase the previous longliving AGW trace in the ionosphere and replace it by lower values of L and A _{ d }. Therefore, this largescale and large amplitude trace stays in the ionosphere without being influenced by subsequent TID phases. Only quasiperiods of a train following much later can establish a new TID with lower values of L and A _{ d }.
During daytime and, particularly around noon when background electron density is increased as compared to nighttime values, the MS AGW train can be observed to be superposed over an existing longliving LS TID. If it propagates with significantly lower values of L and A _{ d }, the original LS TID track becomes modulated by the spatial quasiharmonics of the MS TID.
Figure 21 is similar to Figure 20, but the data were collected during maximum solar activity (from 31/07/2002 (W = 148) to 01/08/2002 (W = 137)). Sharp increases of L on 31/0/2002 at ~19:00 and 22:40, and on 01/08/2002 at 12:27 are attributed to separate precipitation events of energetic protons. Sawtooth pieces are present in the L(t) curves registered both for the minimum and maximum solar activity, but these effects are manifested in different ways.
Fig. 21.
Dynamics of parameters for TIDs and F2 maximum, observed on 31 July and 1 August 2010. 
The maximum values of periods L _{ max } reach 30 000–35 000 km and about 3200–4000 km during daily observations in periods of solar minimum and maximum activity, respectively. The nearly tenfold reduction of the maximum values of L _{ max }, for solar maximum, is explained by a significant decrease of the disturbance lifetime, since the highestspeed LS TIDs could not be recorded because of the limited cadence of our measurements which was equal to ~1.5–2 h. Such a large difference of L _{ max } for minimum and maximum solar activity can be partially explained by the different delays of the first extremum, which is 1–2 h and 11 h for minimum and maximum solar activity, respectively.
We found that regions with reduced electron density were practically never observed under conditions of maximum solar activity, but areas with high density relative to unperturbed background values were detected during the whole measurement period. On 01/08/02 the disturbance amplitude A _{ high } in this region changed in the range 0.3–0.6 for about 20 h according to the quasiharmonic law. Here A _{ high } = ΔN _{ high } /N _{ n }; ΔN _{ high } = N − N _{ n }; N and N _{ n } are disturbed and undisturbed electron concentrations, respectively. The latitude of the center of this region fluctuated in antiphase with A _{ high } around ~45° N. The amplitude of these fluctuations was in the range 3°–5°, and the period was increasing almost linearly from 100 to 200 min during the time interval ~13:15–16:30. Only one LS TID was recorded on 31/07/02 at ~19:30. The TID parameters were found to be L = 5000 km, A _{ d } = 0.9, and A _{ high } = 2.5.
5.7. Spatial dependence of TID inclination angle ξ
From the analysis in Section 2 we deduced that the further a TID quasiperiod is away from the excitation source, the greater is the absolute value of ξ. The effect is clearly visible in Figure 19 (panels 2 and 3) and is summarized in Table 9. The parameters: L, A _{ d }, and the absolute value of ξ increase while the TID is moving away from the excitation source. The front of an MS TID is tilted in the direction of disturbance propagation.
This tendency also holds for LS TIDs. However, if L > 900–1000 km, ξ can be determined with an error not less than ±30°. In many cases, one can only determine a sign of ξ and conclude that ξ is within an expected range (Fedorenko et al. 2012). See, for example, five cases in the article Fedorenko et al. (2012) for L = 1800–2500 km and A _{ d } = 0.5–0.75, obtained at night under solar minimum conditions. For L > 2000–3000 km the absolute value of this angle, in some cases, can be determined only with accuracy δξ ~50–70°. These restrictions are imposed by the geometry of radio sounding for large L.
6. Summary
In this article we have discussed the following observations and their implications:

Largeand mediumscale TIDs are different stages of an initial perturbation generated by the same source. The spatial and temporal dynamics of the TID parameters may be predicted by our semiempirical model. The model was verified via analysis of observations published in the scientific literature. Our results show that the spatial and temporal periods are predicted with accuracy not greater than 12%.

The following parameters of TIDs: the spatial period, amplitude, and front inclination angle in a vertical plain increase with decreasing latitude in almost all observations made during the two solar cycles in the range of latitudes ~33–66° N. This indicates that the sources of most TIDs of natural origin are located at high latitudes.

Spatiotemporal dynamics of the TIDs, generated presumably by the precipitation of energetic protons in the magnetospheric cusp of the northern hemisphere, are described by our empirical model.

The semiempirical model predicts the sawtooth pieces on time dynamics of spatial TID periods observed in our transionospheric radiosounding experiments.

Directly proportional relationship between the amplitude and horizontal spatial period of TIDs was observed during the minimum and maximum of solar activity.

No systematic diurnal and seasonal variability of TID occurrence was found, but the TID occurrence was growing from~50% in 1987 to ~98% in 2010. The results of other studies suggesting that TID occurrence does not depend on sunspot number and magnetic activity are confirmed.

For the period 1987–2010 the occurrence of TIDs, propagating in the freewave mode and waveguide mode, is 76% and 24%, respectively.

The minimum values of the horizontal scale of TIDs do not exceed ~150–200 km, predicted by our model and confirmed by observations.
Acknowledgments
We are grateful to all researchers for their contribution to theoretical and experimental studies of TID since their first observation by Munro (1948). We also express our special sincere gratitude to participants of nuclear experiments for their devotion, and to scientists, who have used these events for geophysical researches and have published their results, frequently with frank remarks and doubts. Their contribution to development of our model is invaluable. We are grateful to the editor and two anonymous reviewers for useful recommendations and questions.
References
 Adushkin, V.V., and K.I. Gorelyi, The response of the ionospheric Fregion to atmospheric nuclear explosions, Dynamic Processes in Geospheres, Geophysics of Strong Disturbances. Collected Scientific Transactions (in Russian), Moscow, 239–248, 1994. [Google Scholar]
 Afraimovich, E.L., S.V. Voyenkov, K.G. Ratovsky, and P.V. Tatarinov, GPSdetection of the single internal gravitational waves generated during powerful magnetic storms. In: Trudy 21 vserossiiskoi nauchnoi konferentsi po rasprostraneniyu radiovoln (Proc. 21th AllRussia Sci. Conference on Radiowave Propagation), Ioshkar Ola, 1, 124–128, 2005. [Google Scholar]
 Albee, P.R., and D.P. Kanellakos, A spatial model of the Fregion ionospheric traveling disturbance following a lowaltitude nuclear explosion, J. Geophys. Res., 73 (3), 1039–1053, 1968. [CrossRef] [Google Scholar]
 Breitling, W.J., and R.A. Kupferman, Traveling ionospheric disturbances associated with nuclear detonations, J. Geophys. Res., 72 (1), 307–315, 1967. [CrossRef] [Google Scholar]
 Bristow, W.A., R.A. Greenwald, and J.C. Samson, Identification of highlatitude acoustic gravity wave sources using the Goose Bay HF radar, J. Geophys. Res., 99 (A1), 319–331, 1994. [CrossRef] [Google Scholar]
 Bullough, R.K., and P.J. Caudrey, Solitons, SpringerVerlag, Berlin, Heidelberg, New York, 1980. [CrossRef] [Google Scholar]
 Cole, J.D., and C. Greifinger, Acousticgravity waves from an energy source at the ground in the isothermal atmosphere, J. Geophys. Res., 74 (14), 3693–3703, 1969. [CrossRef] [Google Scholar]
 Deminov, M.G., A.T. Karpachev, V.V. Afonin, and J. Shmilauer, Change of a position of the main ionospheric trough depending on a longitude and geomagnetic activity, Geomagn. Aeron., 32 (5), 185–188, 1992. [Google Scholar]
 Dorohov, V.L., V.G. Somov, V.N. Fedorenko, and Yu.P. Fedorenko, Diagnostics of medium scale traveling ionospheric disturbances using radio sounding of the ionosphere by signals from loworbiting navigation satellites, Visn. KharkivKarazin Nat. Univ., Ser.: Radiofiz. Elektron. (in Russian), 2 (646), 224–229, 2004. [Google Scholar]
 Evans, J.V., J.M. Holt, and R.H. Wand, A differential doppler study of traveling ionospheric disturbances from Millstone Hill, Radio Sci., 18 (3), 435–451, 1983. [CrossRef] [Google Scholar]
 Fedorenko, V.N., Yu.P. Fedorenko, and I.I. Shagimuratov, Results of the ionosphere study by means of diversity reception of radio signals of loworbiting navigation satellites, Geomagn. Aeron., 37 (3), 346–349, 1997. [Google Scholar]
 Fedorenko, Yu.P., O.F. Tyrnov, and V.N. Fedorenko, Parametres of empirical model of traveling ionospheric disturbances, Elektromagnitnye volny i elektronnye sistemy, 13, 21–46, 2008. [Google Scholar]
 Fedorenko, Yu.P., O.F. Tyrnov, and V.N. Fedorenko, Parameters of traveling ionospheric disturbances estimated from satellite beacon observations in low Earth orbit, Geomagn. Aeron., 50, 489–503, 2010a. [CrossRef] [Google Scholar]
 Fedorenko, Yu.P., V.N. Fedorenko, and V.L. Dorohov, Front inclination in a vertical plane of MediumScale Traveling Ionospheric Disturbances and an effective thickness of a layer of their detection, Visn. Kharkiv Karazin Nat. Univ., Ser.: Radiofiz. Elektron., (in Russian), 17 (927), 109–120, 2010b. [Google Scholar]
 Fedorenko, Yu.P., V.N. Fedorenko, and V.N. Lysenko, Parameters of the mediumscale traveling ionospheric disturbances model deduced from measurements, Geomagn. Aeron., 51 (1), 88–104, 2011. [CrossRef] [Google Scholar]
 Fedorenko, Yu.P., V.N. Fedorenko, and V.L. Dorohov, Diagnostics of parametres of largescale traveling ionospheric disturbances with the help radioscopy ionosphere by signals loworbital navigating satellites, Visn. Kharkiv Karazin Nat. Univ., Ser.: Radiofiz. Elektron., (in Russian), 20 (1010), 97–112, 2012. [Google Scholar]
 Francis, S.H., A theory of mediumscale traveling ionospheric disturbances, J. Geophys. Res., 79 (34), 5245–5260, 1974. [CrossRef] [Google Scholar]
 Francis, S.H., Global propagation of atmospheric gravity waves: A review, J. Atmos. Terr. Phys., 37 (6–7), 1011–1054, 1975. [CrossRef] [Google Scholar]
 Gossard, E.E., and W.H. Hooke, Waves in the Atmosphere, Elsevier, Amsterdam, Netherlands, 456 p., 1975. [Google Scholar]
 Halperin, Yu.I., L.D. Sivtseva, V.M. Filippov, and V.L. Khalipov, The SubAuroral Upper Ionosphere, Nauka, Novosybyrsk (in Russian), 191, 1990. [Google Scholar]
 Heisler, L.H., Anomalies in ionosonde records due to traveling ionospheric disturbances, Aust. J. Phis., 11, 79–90, 1958. [CrossRef] [Google Scholar]
 Hocke, K., and K. Schlegel, A review of atmospheric gravity waves and traveling ionospheric disturbances: 1982–1995, Ann. Geophys., 14, 917–940, 1996. [Google Scholar]
 Igarashi, K., S. Kainuma, I. Nishimuta, S. Okamoto, H. Kuroiwa, T. Tanaka, and T. Ogawa, Ionospheric and atmospheric disturbances around Japan caused by the eruption of Mount Pinatubo on 15 June 1991, J. Atmos. Terr. Phys., 56 (9), 1227–1234, 1994. [CrossRef] [Google Scholar]
 Iyer, K.N., Effect of traveling ionospheric disturbances on HF phasepath measurements, Indian J. Radio Space Phys., 12, 47–49, 1983. [Google Scholar]
 Kanellakos, D.P., Response of the ionosphere to the passage of acousticgravity waves generated be lowaltitude nuclear explosions, J. Geophys. Res., 72 (17), 4559–4576, 1967. [CrossRef] [Google Scholar]
 Karpachev, A.T., N. Beloff, T.D. Carozzi, P.F. Denisenko, T.J.T. Karhunem, and M. Lester, Detection of large scale TIDs associated with the dayside cusp using SuperDARN Data, J. Atmos. Sol. Terr. Phys., 72 (9–10), 653–661, 2010. [CrossRef] [Google Scholar]
 Köhnlein, W., and W.J. Raitt, Position of the midlatitude trough in the topside ionosphere as deduced from ESPO 4 observations, Planet. Space Sci., 25, 600–602, 1977. [CrossRef] [Google Scholar]
 Kotake, N., Y. Otsuka, T. Ogawa, T. Tsugawa, and A. Saito, Statistical study of mediumscale traveling ionospheric disturbances observed with the GPS networks in Southern California, Earth Planets Space, 59, 95–102, 2007. [Google Scholar]
 Kunitsyn, W.E., S.N. Suraev, and R.R. Akhmedov, Modeling propagation of acoustic gravity waves in the atmosphere for different surface sources, Vestn. Mosk. Univ., Sert. 3. Fiz. Astron., 2, 59–63, 2007. [Google Scholar]
 Liu, C.H., and K.C. Yeh, Exciting of acousticgravity waves in an isothermal atmosphere, Tellus, 23, 150–163, 1971. [CrossRef] [Google Scholar]
 Lomax, J.B., and D.L. Nielson, Observation of acousticgravity wave effects showing geomagnetic field dependence, J. Atmos. Terr. Phys., 30, 1033–1050, 1968. [CrossRef] [Google Scholar]
 Lonngren, K., and A. Scott, Soliton in Actions, Academic Press, New York, San FranciscoLondon, 1978. [Google Scholar]
 Munro, G.H., Shortperiod changes in the F region of the ionosphere, Nature, 162, 886–887, 1948. [CrossRef] [Google Scholar]
 Obayashi, T., Widespread ionospheric disturbances due to nuclear explosions during October 1961, Nature, 196 (4849), 24–27, 1962. [CrossRef] [Google Scholar]
 Ogawa, T., K. Igarashi, K. Aikyo, and H. Maeno, Satellite observation of medium scale traveling ionospheric disturbances over Syowa Station, Proc NIPR Symposium Upper Atmosphere, 1, 192–198, 1988. [Google Scholar]
 Otsuka, Y., K. Suzuki, S. Nakagawa, M. Nishioka, K. Shiokawa, and T. Tsugawa, GPS observations of mediumscale traveling ionospheric disturbances over Europe, Ann. Geophys., 31, 163–172, 2013. [CrossRef] [Google Scholar]
 Roberts, D.H., J.A. Klobuchar, P.F. Fougere, and D.H. Hendrickson, A largeamplitude traveling ionospheric disturbance produced by the May 18, 1980, explosion of Mount St. Helens, J. Geophys. Res., 87, 6291–6301, 1982a. [CrossRef] [Google Scholar]
 Roberts, D.H., A.E.E. Rogers, B.R. Allen, C.L. Bennett, B.F. Bukker, P.E. Greenfield, C.R. Lawrence, and T.A. Clark, Radio interferometric detection of a traveling ionospheric disturbance excited by Explosion of Mount St. Helens, J. Geophys. Res., 87 (A8), 6302–6306, 1982b. [CrossRef] [Google Scholar]
 Row, R.V., Acousticgravity waves in the upper atmosphere due to a nuclear detonation and an earthquake, J. Geophys. Res., 72 (5), 1599–1610, 1967. [CrossRef] [Google Scholar]
 RTO Technical Report TRIST051, Characterising the Ionosphere. Final Report of Task Group IST051, Published January 2009. 1.1.1.2 (Hawlitschka Stefan). 1.3.8 (Prikryl Paul). [Google Scholar]
 Senior, A., M.J. Kosch, T.K. Yeoman, M.T. Rietveld, and I.W. McCrea, Effects of high latitude atmospheric gravity wave disturbances on artificial HF radar backscatter, Ann. Geophys., 24, 2347–2361, 2006. [CrossRef] [Google Scholar]
 Shashunkina, V.M., Ionospheric effect of the sudden beginning of a magnetic storm on July, 15th, 1959, Geomagn. Aeron. (in Russian), 6 (1), 146–149, 1966. [Google Scholar]
 Shashunkina, V.M., Study of planetary ionospheric effect of the sudden beginning of a magnetic storm, Ionosf. Issled. (in Russian), 16, 91–96, 1968a. [Google Scholar]
 Shashunkina, V.M., Variations of ionospheric parametres during effect SC depending on solar and magnetic activity according to the Lindaudata, Geomagn. Aeron. (in Russian), 8 (2), 352–354, 1968b. [Google Scholar]
 Shashunkina, V.M., Results of study of ionospheric effect of the sudden beginning of a magnetic storm, Ionosf. Issled. (in Russian), 20, 154–165, 1972. [Google Scholar]
 Song, Q., F. Ding, W.X. Wan, L.B. Liu, and B.Q. Ning, Monitoring nighttime mediumscale traveling ionospheric disturbances using the GPS network over North America, J. Geophys., 54, 162–168, 2011. [Google Scholar]
 Taran, V.I., Yu.I. Pod’yachii, V.I. Golovin, V.I. Vashchenko, and I.D. Arkad’ev, Traveling ionospheric disturbances detected using the incoherent scatter method, Ionos. Issled., 27, 102–110, 1979. [Google Scholar]
 Taran, V.I., Yu.I. Pod’yachii, A.N. Smirnov, and L.J. Gerstein, Disturbances of the ionosphere after a ground level burst on supervision by a method of incoherent scatter, Izvest. akadem. nauk SSSR, Phisica Zemly (in Russian), 11, 75–79, 1985. [Google Scholar]
 Tsugawa, T., N. Kotake, Y. Otsuka, and A. Saito, Mediumscale traveling ionospheric disturbances observed by GPS receiver network in Japan, A short review GPS Solutions, 11 (2), 139–144, DOI: 10.1007/s1029100600455 Key: citeulike:1207528, 2007a. [CrossRef] [Google Scholar]
 Tsugawa, T., Y. Otsuka, A.J. Coster, and A. Saito, New characteristics of mediumscale traveling ionospheric disturbances detected with dense widecoverage TEC maps over North America, Geophys.Res. Lett., 34, 22–101, DOI: 10.1029/2007GL031663, 2007b. [CrossRef] [Google Scholar]
 Tyrnov, O.F., Yu.P. Fedorenko, and L.F. Chernogor, Studying wavelike disturbances of electron density using radio sounding of the ionosphere by coherent signals of navigation satellites, Usp. Sovrem. Radioelektron., 1, 36–80, 2005. [Google Scholar]
 Vadas, S.L., and G. Crowley, Sources of the traveling ionospheric disturbaces observed by the ionospheric TIDDBIT sounder near Wallops Island on 30 October 2007, J. Geophys. Res., 115 (A07324), 1–24, 2010. [CrossRef] [Google Scholar]
 Vlasov, A., K. Kauristie, M. van de Kamp, J.P. Luntama, and A. Pogoreltsev, A study of traveling ionospheric disturbances and atmospheric gravity waves using EISCAT Svalbard Radar IPYdata, Ann. Geophys., 29, 2101–2116, 2011. [CrossRef] [Google Scholar]
 Wexler, H., and W.A. Hass, Global atmospheric pressure effects of the October 30, 1961, explosion, J. Geophys. Res., 67 (10), 3875–3879, 1962. [CrossRef] [Google Scholar]
 Yakovets, A.F., V.V. Vodyannikov, G.I. Gordienko, J.F. Ashkaliev, Yu.G. Litvinov, and S.B. Akasov, The response of a night middlelatitude ionosphere to passage of an atmospheric gravitational wave, Geomagn. Aeron., 48 (4), 534–541, 2008. [CrossRef] [Google Scholar]
 Yeh, K.C., and C.H. Liu, Acousticgravity waves in the upper atmosphere, Rev. Geophys. Space Phys., 12 (2), 193–616, 1974. [CrossRef] [Google Scholar]
Cite this article as: Fedorenko Y P, Tyrnov O, Fedorenko V & Dorohov V: Model of traveling ionospheric disturbances. J. Space Weather Space Clim., 2013, 3, A30.
All Tables
Seasonal occurrence of TIDs, propagating in waveguide (WGM) and freewave mode (FWM).
The distances between the MS TID and measuring station, shown in Figures 14 and 15.
All Figures
Fig. 1.
Temporal halfperiods τ _{ i } versus distance R to the source, predicted by the model: panel 1 for i = 1–14, panel 2 for i = 1–10. 

In the text 
Fig. 2.
Spatial halfperiods l _{ i } versus time t, predicted by the model. Zerotime reference corresponds to the time of disturbance initiation. 

In the text 
Fig. 3.
Spacetime dynamics of AGW train, consisting of 14 halfperiods. 

In the text 
Fig. 4.
Spatial halfperiods L versus distances to a disturbance source R for various values of R _{1}. 

In the text 
Fig. 5.
Quasiperiod L versus time t for various distances R _{ S } between the measuring station and TID sources. Zerotime reference corresponds to disturbance initiation time. 

In the text 
Fig. 6.
Inclination angles ξ versus R: 1 – MS TID, z _{0} = 100 km; 2 – MS TID, z _{0} = 200 km; 3 – LS TID. 

In the text 
Fig. 7.
(a) Extrema velocities of a TID train versus its index for surface nuclear explosions (panels 1, 2, and 3) and volcano eruption (panel 4). (b) Extrema velocities of TID train versus its index for highaltitude nuclear explosions. 

In the text 
Fig. 8.
Scans of AGW train halfperiods. 1 – experimental values (Albee & Kanellakos 1968), 2 – values predicted by the model. 

In the text 
Fig. 9.
Lengths of halfperiods versus distance between their centers and a source. 1 – experimental values (Albee & Kanellakos 1968); 2 – values predicted by the model. 

In the text 
Fig. 10.
Temporal periods versus a distance R to the source. 1 – values predicted by the model; 2 – experimental values (Albee & Kanellakos 1968). 

In the text 
Fig. 11.
The third (panel a) and fourth (panel b) temporal halfperiods versus a distance R to the source. 1 – values predicted by the model; 2 – experimental values (Albee & Kanellakos 1968). 

In the text 
Fig. 12.
Temporal period T versus distance R to the source 1 – values predicted by the model; 2 – Row (1967); 3 – Roberts et al. (1982a). 

In the text 
Fig. 13.
TID occurrence and solar spot numbers in 1987–2010: 1 – linear regression of TID occurrence; 2 – W/100, W is yearly average Wolf number; 3 – TID occurrence without a linear trend, 4 – the occurrence of TIDs propagating in waveguide and freewave modes and their linear regressions for the period from 1997 to 2010. 

In the text 
Fig. 14.
The results of MS TID observation for ξ = 42 – 75°. Parameters for the most southern irregularity are given in lower part of panels 1, 2, and 3; parameters for the rest of irregularities are given in upper part of these panels. 

In the text 
Fig. 15.
The results of MS TID observation for ξ = –(50 – 62)°. 

In the text 
Fig. 16.
1 – the experimental curve R(L); 2 – the curve R(L), predicted by the model. 

In the text 
Fig. 17.
The plots δf _{ i }(φ) and L(φ), registered after South Sumatra Earthquake. 

In the text 
Fig. 18.
The recordings of δfi(t) obtained after South Sumatra Earthquake (solid lines) and their model approximations for undisturbed conditions (left panel) and for disturbed conditions (right panel). “SN” – satellite moves from the south to the north; “NS” – satellite moves from the north to the south. 

In the text 
Fig. 19.
Examples of MS TID with A _{ d } = 0.18 … 0.95 and ξ = –45 … –80°. The recordings of δfi(t) obtained after precipitation energetic protons near the magnetospheric cusp (solid lines) and their model approximations for undisturbed conditions (left panel) and for disturbed conditions (right panel). “SN” – satellite moves from the south to the north; “NS” – satellite moves from the north to the south. 

In the text 
Fig. 20a.
The diurnal dependencies of L, A _{ d }, z _{ m }, and N _{ m }, registered on 2 and 3 of March 2010. 

In the text 
Fig. 20b.
The diurnal dependencies of L, A _{ d }, z _{ m } and N _{ m }, registered on 4 and 5 of March 2010. 

In the text 
Fig. 21.
Dynamics of parameters for TIDs and F2 maximum, observed on 31 July and 1 August 2010. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.