 Letter
 Open Access
 Published:
Statistical study of mediumscale traveling ionospheric disturbances using SuperDARN Hokkaido ground backscatter data for 2011
Earth, Planets and Space volume 67, Article number: 22 (2015)
Abstract
We describe an automated technique to determine parameters of traveling ionospheric disturbances (TIDs) using the Super Dual Auroral Radar Network (SuperDARN) high frequency (HF) radar data. The technique is based on the analysis of minimum ground backscatter range variations corresponding to different radar beams. Using this technique, we processed the SuperDARN Hokkaido radar data for 2011 and revealed statistical distributions of mediumscale TID (MSTID) azimuth and apparent horizontal velocity. We found four peaks with a distinct diurnal and seasonal dependence in the MSTID azimuth occurrence rate distributions. Northeast MSTID azimuths (20° to 50°) are typical of the summer and equinox morning hours; southeast azimuths (100° to 140°) prevail in the winter daytime; southwest azimuths (190° to 220°) are typical mostly in the summer and equinox nighttime and in the equinox evening; northwest azimuths (280° to 320°) are typical of the summer daytime and evening. The apparent horizontal velocities are generally within the 100 to 160 m/s range. The obtained results agree well with earlier studies by other researchers. However, there are also certain differences. The summer daytime northwestward MSTIDs are not indicated in the earlier studies. The nighttime horizontal velocities are 1.5 to 2 times higher than those in the daytime. Furthermore, winter velocity values are about 1.5 times higher than those in other seasons. These differences might be associated with the peculiarities of the data recorded by different facilities, or the features of the processing techniques, and require further investigation for their interpretation.
Correspondence/findings
Introduction
By traveling ionospheric disturbances (TIDs), we understand wavelike electron density disturbances caused by atmospheric gravity waves (AGWs) propagating at ionospheric heights. One should distinguish between mediumscale and largescale TIDs (Hunsucker 1982). Mediumscale TIDs (MSTIDs) have horizontal wavelengths of about several hundreds of kilometers and periods of 15 min up to 1 h. Largescale TIDs (LSTIDs) have horizontal wavelengths of more than 1,000 km and periods of 30 min to 3 h.
TIDs have been studied for decades (Hocke and Schlegel 1996; Hunsucker 1982), including studies that have used the Super Dual Auroral Radar Network (SuperDARN) high frequency (HF) radar data (Bristow and Greenwald 1995; Bristow et al. 1994; Chisham et al. 2007; Karhunen et al. 2006; Samson et al. 1990; Ogawa et al. 2012). A TID passage within the radar fieldofview leads to HF wave focusing/defocusing effects. The latter manifests themselves in periodic variations of the received ground backscatter power (Samson et al. 1990). One of the ground backscatter characteristics most sensitive to TIDs and measured by SuperDARN radar is the minimum range corresponding to the skip distance (Arnold et al. 1998; Stocker et al. 2000). Karhunen et al. (2006) proposed a technique involving a special multifrequency mode (METSO). This technique uses skip distance measurements to build a twodimensional phase map of the disturbance. From this map, one can obtain the azimuth, horizontal velocity, and horizontal wavelength as a function of the disturbance period.
In contrast to different case studies examining the properties of individual TIDs, statistical studies, using wide datasets, enable us to reveal the general pattern of AGW propagation. A number of papers presenting statistical studies of TIDs have been published in the last 15 years (Afraimovich et al. 1999, 2000; Ding et al. 2008; Kotake et al. 2007; Medvedev et al. 2013; Otsuka et al. 2011; Shiokawa et al. 2003; Tsugawa et al. 2004). There are also statistical studies that utilized SuperDARN radar data (Bristow et al. 1996; He et al. 2004; Ichihara et al. 2013; Ishida et al. 2008; Grocott et al. 2013). Apparently, a relatively small number of such studies are related to the complexity of extensive dataset processing, as well as the issues concerning an automated interpretation of HF radar data.
In this paper, we present a statistical analysis of MSTID azimuth and apparent horizontal velocity calculated using the SuperDARN Hokkaido HF radar data obtained in 2011. We used an automated technique where a crosscorrelation analysis of minimum ground backscatter range variations is used to determine TID parameters. We did not place limitations on the spatial scale or temporal periodicity of the observed disturbances, but we analyzed the minimum range variations as they were observed. A statistical sampling for the crosscorrelation function calculation is varied within the 1.5 to 7 h interval. Therefore, both MSTIDs and LSTIDs are included in the analysis, but we suppose the LSTID effect is negligible in our statistical results because LSTIDs are observed about an order of magnitude less frequently. In addition, we focused on investigating the seasonal and diurnal dependence of the above parameters to compare them with other recent statistical studies.
Instrumentation and analysis method
We used the 2011 SuperDARN Hokkaido HF radar (43.53° N, 143.61° E) ground backscatter data in this study (http://center.stelab.nagoyau.ac.jp/hokkaido/index.html). The phasedarray antenna system enables the radiation pattern to be focused in one of 16 azimuth directions (narrow beam with approximately 3.5° width). The radar measures echoes from ranges of up to approximately 3,500 km. Figure 1 shows a map of the radar location, as well as all 16 beam directions that form the approximately 50degreewide fieldofview. In the standard operational mode, the radar continuously scans the entire fieldofview. Scanning takes either 1 or 2 min, then the process is repeated depending upon the integration time along each beam, which is 3 or 6 s, respectively. The operating frequency is selected from the allowed frequency range and remains relatively constant throughout a day. We analyzed the data obtained at the approximately 11 MHz operating frequency because it was mostly used in 2011. The SuperDARN Hokkaido radar measures rangetime dependencies of the following echo parameters at each beam: power, Doppler velocity, Doppler spectral width, and elevation angle.
Ground backscatter is one of two echo types recorded by the radar. In addition to ground backscatter echoes, there are also echoes backscattered directly by ionospheric irregularities. Separation between the two echo types is performed in the standard way (He et al. 2004). The ground backscatter echoes are supposed to have a Doppler velocity and Doppler spectral width of typically less than 50 and 30 m/s, respectively. Figure 2a illustrates an example of the SuperDARN Hokkaido data, showing the ground backscatter echo power received by the radar along Beam 0 versus slant range and Japanese standard time (JST = UT + 9 h) on November 24, 2011. The oblique lines with enhanced power clearly seen in Figure 2a manifest the MSTID passage through the radar fieldofview.
Data preprocessing
As shown by Milan et al. (1997), SuperDARN radars can receive ground backscatter echoes not only from the front but also from behind the radar. In the case of the Hokkaido radar, echoes from behind might be observed more frequently during daytime at low solar activity, when the critical frequency is reduced so that there is no ionospheric reflection in front of the radar. However, HF waves may be reflected from behind the radar from northern slope of the equatorial anomaly, where the critical frequency is greater. To exclude the echoes coming from the rear fieldofview, additional interferometer data (elevation angle) processing is provided, according to Milan et al. (1997).
We are interested in MSTIDs propagating at the heights of the F2 ionospheric layer. Except for the regular onehop propagation ground backscatter echoes reflected by the F2 ionospheric layer, one may also record other echoes. Among them, there are regular echoes reflected by the E layer, as well as twohop propagation echoes reflected by the E and F2 layers. In Figure 2a, a regular F2 onehop propagation echo is seen at 6:00 through 18:00 JST within the 700 to 1,700 km range, simultaneously with an F2 twohop propagation echo (above approximately 1,800 km range), and with other (irregular) echoes at nighttime. To correctly exclude the echoes not related to the F2 onehop propagation, we developed a model of ground backscatter characteristics. It describes diurnal and seasonal variations of the F2 onehop propagation ground backscatter characteristics for all beams and for the main operating frequencies of the SuperDARN Hokkaido radar. The model is based on the ground backscatter simulation using the original software developed within the waveguide approach and the International Reference Ionosphere model (Bilitza and Reinisch 2008). The simulation was performed for the 15th day of each month. Figure 3 shows, for example, diurnal and seasonal dependencies of four model parameters: minimum range, skip distance, elevation angle, and ratio between the total range to ground scatter and the range to the ionospheric reflection point (the ray trajectory apogee). Figure 3 is a plot for Beam 0, with an operating frequency of 11 MHz, for December 2010 through January 2012. The range ratio (see Figure 3d) characterizes the asymmetry of the ray trajectory caused by the regular ionospheric gradients. As we can see for poleward trajectories, this ratio is greater than two. In addition, the model does not predict ground backscatter echoes at nighttime in October to April. This is because of the low nighttime critical frequency during these months, resulting in the HF wave passing through the ionosphere without reflection. By linear interpolation, the model allows us to calculate a mask in which the regular ground backscatter echo should be for each day of the year.
Further preprocessing is illustrated in Figure 2b, where the red and black lines show the calculated model minimum range and mask, respectively, for Beam 0 on November 24, 2011. The squares show the ranges where the maxima of the received ground backscatter power were found within the mask at each moment in time. The squares represent the required dependence of the slant range versus time. We fit the dependencies for all beams by penalized regression splines (not shown in Figure 2b) which is found to be good because of adaptive control over smoothing (Ruppert et al. 2003). Then, using splines, we recalculate the dependencies on an equidistant time grid common for all beams. As a final step, we remove the diurnal trend from the dependencies to obtain the minimum range variations corresponding to the TID effects.
Technique for TID parameter estimation
The technique for TID parameter estimation (Oinats et al. 2012, 2013) is based on three assumptions. The first two assumptions, concerning the TID’s model, are that a disturbance has a plane phase front within the radar fieldofview, and the phase velocity vertical component is equal to zero. The third assumption is that the ionosphere mainly influences the ground backscatter at the apogee of the propagation trajectory, i.e., where reflection occurs. Using the third assumption, we can consider the minimum range variations caused by ionosphere variations that occur directly in the reflection point. Figure 1 shows the approximate location of the ionospheric reflection points for the SuperDARN Hokkaido radar (red circles).
The disturbance’s phase is a function of observational point radius vector, \( \overrightarrow{r} \), and of time, t. According to the first assumption, we have:
where \( \overrightarrow{k} \) and ω are the TID wave vector and cyclic frequency, respectively, and ϕ _{0} is an initial phase. The radar fieldofview is about 50° in the azimuth direction. This corresponds to about 500 km between reflection points for Beams 0 and 15. Therefore, we should write Expression 1 in terms of spherical coordinates with regard to the radar location. For the case of zero vertical phase velocity, we can write:
where Φ is the TID azimuth (clockwise from north), \( k=\left\overrightarrow{k}\right=2\pi /\lambda \), λ is the TID spatial wavelength, (θ′, φ ') are the latitude and longitude of the radar location, and r _{0} is the Earth’s radius. Thus, the time lag between variations corresponding to the ith and jth radar beams is given by expression (Oinats et al. 2013):
where υ = ω/k is the TID phase velocity, and (θ _{i}, φ _{i}) and (θ _{j}, φ _{i}) are the coordinates (latitude and longitude) of the reflection points for the ith and jth beams. Excluding the i = j case, we can obtain 240 equations of type Equation 3 with two unknown parameters for all 16 radar beams. Certainly, one needs only two combined equations to determine two unknown parameters. This corresponds to measuring the variations at three reflection points. However, we should keep in mind the possible inaccuracy of data processing. For example, an ambiguity might arise when minimum range variations are extracted from the radar data. The coordinates of the reflection points are also determined approximately. Therefore, the solution for the equation set is found through twodimensional leastsquares fitting. The discrepancy between the observed and the theoretical time lags is given by:
The observed time lag, \( \Delta {t}_{ji}^{\prime } \), for each pair of variations is determined as a time shift corresponding to the maximum of the crosscorrelation function. The statistical sampling size corresponds to the prevailing period in the variation’s spectrum within the 0.8 to 5 h interval.
To determine the coordinates of the reflection points, we first calculate the ground distance to the reflection point using the slant range and the wellknown formulas based on Breit and Tuve’s theorem (Davies 1990). Two parameters are needed: the slant range to the ionospheric reflection point and either the effective reflection height or the elevation angle. Unfortunately, the elevation angle measurements made by the Hokkaido radar have not as yet been calibrated. Therefore, we assume the effective height of the ionospheric reflection is constant at 250 km for all ranges (we found that this value was appropriate for the 2011 Hokkaido radar data). Such an approach is quite reasonable and has been used earlier by other researchers (Bristow et al. 1994; Grocott et al. 2013). In addition, as shown by Karhunen et al. (2006), even when sounding is done at different operating frequencies, the ionospheric reflection height remains approximately the same.
The slant range to the reflection point is equal to half of the slant range to the ground scatter in the case of a spherically symmetric ionosphere. However, in general, there might be an asymmetry of the trajectory due to ionospheric gradients. Our estimations show that the difference for the ground distance to the reflection point can reach up to 15% because of the regular ionospheric gradients for the 2011 Hokkaido radar data. This leads to the emergence of corresponding errors in the TID horizontal velocity estimations, especially for diurnal intervals, when the solar terminator passes through the radar fieldofview. To minimize the error, we calculate the slant range to the reflection point using the model range ratio (see Figure 3d). Finally, we calculate the appropriate latitude and longitude with reference to the beam azimuth and the radar coordinates.
In Equation 4, we take into account only those time lags, for which the correlation coefficient exceeds the 0.85 threshold. Minimization of the discrepancy is performed using the LevenbergMarquardt algorithm for solving nonlinear least squares problems. If the number of the time lags satisfying the mentioned threshold is less than three, minimization is not performed.
The above processing was carried out for each local time. In Figure 4a, the calculated diurnal variations of the azimuth and horizontal velocity are shown for November 24, 2011 (see Figure 2). As one can see, the azimuth and horizontal velocity vary during the day. There are ‘jumps’ that can be connected with uncertainties in time lag calculations because of the presence of components with smaller periods in the variations. However, the common diurnal trends are seen clearly enough. The azimuth in the morning hours is relatively stable, at about 210°, then changes rapidly to approximately 120° at 10:00 JST, then increases again to approximately 220°. The disturbance propagation direction changes clockwise from southeast to southwest during the daytime. The velocity increases from approximately 150 m/s in the morning to approximately 350 m/s in the evening. Figure 4b,c shows the corresponding distributions of the azimuths and velocities, respectively. Because the parameters vary during the day, there are no single maxima in the distributions, but the ranges of parameter change are clear enough.
Results and discussion
All the SuperDARN Hokkaido radar data for 2011 were processed in the manner described above. Figure 5 shows the relative occurrence rate of azimuths and of apparent horizontal velocities for the entire year. Each bar in the distribution is proportional to the number of cases for which the calculated azimuth (or velocity) is within the given 10degreewidth interval (or 20 m/s). Only the variations for which the average root mean square (RMS) over all beams exceeded 60 km were taken into account (thus, the variation amplitude exceeds the radar range gate). As seen in Figure 5a, there are four dominant MSTID propagation directions: northeast (approximately 35°), southeast (approximately 120°), southwest (approximately 205°), and northwest (approximately 295°). The relative occurrence rate for these directions is approximately 1.5 to 2 times higher than the background rate. The most frequent horizontal velocities are within the 100 to 160 m/s range.
To study the diurnal dependence, we calculated the relative occurrence rate over each hour of the diurnal interval as was done by Kotake et al. (2007) and Otsuka et al. (2011), who used the GPS total electron content (TEC) maps over Japan (GEONET) to determine MSTID parameters. Figures 6a and 7a show the results for the whole year for azimuth and velocity, respectively. As we can see, there is a distinct diurnal dependence of the four dominant MSTID groups. Azimuths corresponding to the northeast peak (20° to 50°) are typical of the morning hours (4 to 6 JST). Southeast azimuths (100° to 140°) prevail in the daytime (8 to 13 JST). Southwest azimuths (190° to 220°) are typical mostly at night and in the evening (16 to 2 JST). Northwest azimuths (280° to 320°) are typical of daytime and in the evening (12 to 17 JST).
We calculated similar distributions separately for three seasons in 2011: equinox (March, April, September, and October), summer (May, June, July, and August) and winter (November, December, January, and February). As seen in Figure 6d, the daytime southward MSTIDs are most typical of the winter. Furthermore, the propagation direction of such disturbances changes clockwise during the day from the southeast (approximately 100°) at 9 JST to the southwest (approximately 210°) at 17 JST. Such behavior of daytime MSTIDs was indicated earlier and is well explained by the filtering properties of the neutral wind (Afraimovich et al. 1999; Kotake et al. 2007). However, Otsuka et al. (2011) did not find such a distinct dependence for the Japan region, indicating only the southsoutheast direction for the winter daytime MSTIDs. Ishida et al. (2008) also indicated the dominant southsoutheast direction for the winter daytime MSTIDs.
As seen in Figure 6b,c, the nighttime southwestward MSTIDs are typical of the summer and equinoxes. Other authors (Afraimovich et al. 1999; Kotake et al. 2007; Otsuka et al. 2011) found that such MSTIDs were also typical in winter. Ichihara et al. (2013) noted that approximately two thirds of the nighttime MSTIDs propagate southwestward, regardless of the season. Southwestward propagation of most of the observed nighttime MSTIDs was indicated by Shiokawa et al. (2003), who used the 630nm airglow images at two stations in Japan from 1998 till 2000. The nighttime southwestward MSTIDs were also mentioned by Medvedev et al. (2013), who analyzed the Irkutsk incoherent scatter radar and the ionosonde DPS4 data using the technique for full TID wave vector determination. The occurrence rate is seen to be weak in the winter nighttime (see Figure 6d) because the regular ground backscatter echoes are almost absent during this period. Morning northeastward MSTIDs are observed mainly in the summer and equinoxes (see Figure 6b,c, respectively). The evening northwestward MSTIDs are typical of summer (see Figure 6c). This is also consistent with Kotake et al. (2007) and Otsuka et al. (2011).
Thus, our results agree well with previous studies. However, there are also several differences. According to Figure 6c, the northwestward MSTIDs are typical of the summer daytime. Other studies do not indicate this. As clearly seen in Figure 7a, the horizontal velocity nighttime values are 1.5 to 2 times higher than those for the daytime. According to Medvedev et al. (2013), the absolute velocity nighttime values are approximately two times higher than those for the daytime. However, there is no diurnal dependence found for the apparent horizontal velocity. It is also interesting that, according to Figure 7d, the horizontal velocities in winter are about 1.5 times greater than those in other seasons. Explanation for the mentioned differences requires further investigation. Most likely, these differences can be related to the peculiarities of the data recorded by different facilities or the features of the processing techniques. In particular, one of the issues with the Hokkaido radar is an absence of calibrated elevation angle measurements. This leads to difficulties in accurate identification of the overlapping E and F2 ground backscatter echoes, which is the case during the summer daytime. Another related issue is the rough mapping of the ionospheric reflection points. The latter mainly affects the accuracy of the horizontal velocity determination as discussed by Grocott et al. (2013).
Conclusion
We developed an original automated technique for SuperDARN HF radar data processing to determine the main TID parameters. The technique is based on crosscorrelation analysis of the ground backscatter minimum range variations at different radar beams. The 2011 SuperDARN Hokkaido radar dataset was processed using this technique. We calculated the occurrence rates of the MSTID azimuth and of MSTID apparent horizontal velocity for the entire year, as well as for different seasons, depending on the local time.
There are four peaks with a distinct diurnal and seasonal dependence in the MSTID azimuth occurrence rate distributions. Northeast MSTID azimuths (20° to 50°) are typical of the summer and equinox morning hours; southeast azimuths (100° to 140°) prevail in the winter daytimes; southwest azimuths (190° to 220°) are typical mostly in the summer, and equinox nighttime and in the equinox evening; northwest azimuths (280° to 320°) are typical of the summer daytime and evening. The apparent horizontal velocities are generally within the 100° to 160 m/s range. The results agree well with earlier studies by other researchers. However, there are several differences whose explanation requires further investigation.
The technique described here can be applied to the extensive SuperDARN Hokkaido radar dataset, as well as to the data of other SuperDARN radars. Such a study, if performed, could be useful for better understanding of AGW propagation and generation mechanisms.
Abbreviations
 AGW:

atmospheric gravity waves
 HF:

high frequency
 LSTID:

largescale traveling ionospheric disturbance
 MSTID:

mediumscale traveling ionospheric disturbance
 SuperDARN:

Super Dual Auroral Radar Network
 TID:

traveling ionospheric disturbance
References
Afraimovich EL, Boitman ON, Zhovty EI, Kalikhman AD, Pirog TG (1999) Dynamics and anisotropy of traveling ionospheric disturbances as deduced from transionospheric sounding data. Radio Sci 34:477–487
Afraimovich EL, Lipko YV, Vugmeister B (2000) Determining dynamic parameters of differentscale ionospheric irregularities over northern Siberia. J Atmos SolTerr Phys 62:133–140
Arnold NF, Jones T, Robinson T (1998) Validation of the CUTLASS HF radar gravity wave observing capability using EISCAT CP1 data. Ann Geophys 16:1392–1399
Bilitza D, Reinisch BW (2008) International reference ionosphere 2007: improvements and new parameters. Adv Space Res 42:599–609, doi:10.1016/j.asr.2007.07.048
Bristow WA, Greenwald RA (1995) Estimating gravity wave parameters from oblique highfrequency backscatter: modeling and analysis. J Geophys Res 100(A3):3639–3648
Bristow WA, Greenwald RA, Samson JC (1994) Identification of highlatitude acoustic gravity wave sources using the Goose Bay HF radar. J Geophys Res 99(A1):319–331
Bristow WA, Greenwald RA, Villain JP (1996) On the seasonal dependence of mediumscale atmospheric gravity waves in the upper atmosphere at high latitudes. J Geophys Res 101(A7):15685–15699
Chisham G, Lester M, Milan SE, Freeman MP, Bristow WA, Grocott A, McWilliams KA, Ruohoniemi JM, Yeoman TK, Dyson PL, Greenwald RA, Kikuchi T, Pinnock M, Rash JPS, Sato N, Sofko GJ, Villain JP, Walker ADM (2007) A decade of the Super Dual Auroral Radar Network (SuperDARN): scientific achievements, new techniques and future directions. Surv Geophys 28:33–109, doi:10.1007/s1071200790178
Davies K (1990) Ionospheric radio. Peter Peregrinus Ltd, London
Ding F, Wan W, Liu L, Afraimovich EL, Voeykov SV, Perevalova NP (2008) A statistical study of largescale traveling ionospheric disturbances observed by GPS TEC during major magnetic storms over the years 20032005. J Geophys Res. doi:10.1029/2008JA013037
Grocott A, Hosokawa K, Ishida T, Lester M, Milan SE, Freeman MP, Sato N, Yukimatu AS (2013) Characteristics of mediumscale traveling ionospheric disturbances observed near the Antarctic Peninsula by HF radar. J Geophys Res Space Phys 118:5830–5841, doi:10.1002/jgra.50515
He LS, Dyson P, Parkinson ML, Wan W (2004) Studies of medium scale travelling ionospheric disturbances using TIGER SuperDARN radar sea echo observations. Ann Geophys 22:4077–4088
Hocke K, Schlegel K (1996) A review of atmospheric gravity waves and travelling ionospheric disturbances: 19821995. Ann Geophys 14:917–940
Hunsucker RD (1982) Atmospheric gravity waves generated in the highlatitude ionosphere: a review. Rev Geophys Space Phys 20(2):293–315, doi:10.1029/RG020i002p00293
Ichihara A, Nishitani N, Ogawa T, Tsugawa T (2013) Northwardpropagating nighttime mediumscale traveling ionospheric disturbances observed with SuperDARN Hokkaido HF radar and GEONET. Adv Polar Sci 24:42–49, doi:10.3724/SP.J.1085.2013.00042
Ishida T, Hosokawa K, Shibata T, Suzuki S, Nishitani N, Ogawa T (2008) SuperDARN observations of daytime MSTIDs in the auroral and midlatitudes: possibility of longdistance propagation. Geophys Res Lett 35:L13102, doi:10.1029/2008GL034623
Karhunen TJT, Robinson TR, Arnold NF, Lester M (2006) Determination of the parameters of travelling ionospheric disturbances in the highlatitude ionosphere using CUTLASS coherent scatter radars. J Atmos SolTerr Phys 68:558–567, doi:10.1016/j.jastp.2005.03.021
Kotake N, Otsuka Y, Ogawa T, Tsugawa T, Saito A (2007) Statistical study of mediumscale traveling ionospheric disturbances observed with the GPS networks in Southern California. Earth Planets Space 59:95–102
Medvedev AV, Ratovsky KG, Tolstikov MV, Alsatkin SS, Scherbakov AA (2013) Studying of the spatialtemporal structure of wavelike ionospheric disturbances on the base of Irkutsk incoherent scatter radar and Digisonde data. J Atmos SolTerr Phys 105–106:350–357, doi:10.1016/j.jastp.2013.09.001
Milan SE, Jones TB, Robinson TR, Thomas EC, Yeoman TK (1997) Interferometric evidence for the observation of ground backscatter originating behind the CUTLASS coherent HF radars. Ann Geophys 15:29–39
Ogawa T, Nishitani N, Tsugawa T, Shiokawa K (2012) Giant ionospheric disturbances observed with the SuperDARN Hokkaido HF radar and GPS network after the 2011 Tohoku earthquake. Earth Planets Space 64(12):1295–1307
Oinats AV, Kurkin VI, Kutelev KA, Nishitani N (2012) The outlook of SuperDARN radars application for monitoring of the ionospheric dynamics in Russia. Physical Bases of Instrumentation 1(3):3–18 (in Russian)
Oinats AV, Kurkin VI, Nishitani N, Saito A (2013) On the determination of traveling ionospheric disturbances parameters using SuperDARN radar data. Electromagnetic Waves and Electronic Systems 18(8):30–39 (in Russian)
Otsuka Y, Kotake N, Shiokawa K, Ogawa T, Tsugawa T, Saito A (2011) Statistical study of mediumscale traveling ionospheric disturbances observed with a GPS receiver network in Japan. In: Abdu M, Pancheva D, Bhattacharyya A (eds) Aeronomy of the Earth’s atmosphere and ionosphere. Springer Netherlands, Dordrecht
Ruppert D, Wand MP, Carroll RJ (2003) Semiparametric regression. Cambridge University Press, New York
Samson JC, Greenwald RAA, Ruohoniemi JM, Frey A, Baker KB (1990) Goose Bay radar observations of Earthreflected, atmospheric gravity waves in the highlatitude ionosphere. J Geophys Res 95(A6):7693–7709
Shiokawa K, Ihara C, Otsuka Y, Ogawa T (2003) Statistical study of nighttime mediumscale traveling ionospheric disturbances using midlatitude airglow images. J Geophys Res 108(A1):1052, doi:10.1029/2002JA009491
Stocker AJ, Arnold NF, Jones TB (2000) The synthesis of travelling ionospheric disturbance (TID) signatures in HF radar observations using ray tracing. Ann Geophys 18:56–64, doi:10.1007/s0058500000564
Tsugawa T, Saito A, Otsuka Y (2004) A statistical study of largescale traveling ionospheric disturbances using the GPS network in Japan. J Geophys Res 109:A06302, doi:10.1029/2003JA010302
Acknowledgements
We would like to thank Dr. K. G. Ratovsky and Dr. A. V. Medvedev at the Institute of SolarTerrestrial Physics of the Siberian Branch of the Russian Academy of Sciences for useful discussions and helpful comments. We are very grateful to the referees for their appropriate and constructive suggestions. This study was done with financial support from the Russian Foundation for Basic Research (under grants No. 140500259, No. 140500588, and No. 120500865a).
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
AO elaborated the technique of TID parameter estimation, provided the processing and calculations, and drafted the manuscript. VK proposed the idea of the study and supervised the work of AO at ISTP. NN provided the Hokkaido SuperDARN radar data, supervised the work of AO at STEL, and participated in the drafting the manuscript. All the authors participated in the discussion and interpretation of the obtained statistical results. All authors read and approved the final manuscript.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.
The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.
To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Oinats, A.V., Kurkin, V.I. & Nishitani, N. Statistical study of mediumscale traveling ionospheric disturbances using SuperDARN Hokkaido ground backscatter data for 2011. Earth Planets Space 67, 22 (2015). https://doi.org/10.1186/s4062301501924
Received:
Accepted:
Published:
Keywords
 Mediumscale traveling ionospheric disturbances
 Atmospheric gravity waves
 SuperDARN Hokkaido radar