Next Article in Journal
Enabling Fine Sample Rate Settings in DSOs with Time-Interleaved ADCs
Previous Article in Journal
Evaluation of Electrochromic Properties of Polypyrrole/Poly(Methylene Blue) Layer Doped by Polysaccharides
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Long-Term Observation of the Quasi-3-Hour Large-Scale Traveling Ionospheric Disturbances by the Oblique-Incidence Ionosonde Network in North China

Electronic Information School, Wuhan University, Wuhan 430072, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(1), 233; https://doi.org/10.3390/s22010233
Submission received: 1 December 2021 / Revised: 25 December 2021 / Accepted: 28 December 2021 / Published: 29 December 2021
(This article belongs to the Section Remote Sensors)

Abstract

:
The oblique-incidence ionosonde network in North China is a very unique system for regional ionospheric observation. It contains 5 transmitters and 20 receivers, and it has 99 ionospheric observation points between 22.40° N and 33.19° N geomagnetic latitudes. The data of the ionosonde network were used to investigate the statistical characteristics of the quasi-3-h large-scale traveling ionospheric disturbances (LSTIDs). From September 2009 to August 2011, 157 cases of the quiet-time LSTIDs were recorded; 110 cases traveled southward, 46 cases traveled southwestward and only 1 case traveled southeastward. The LSTIDs mainly appeared between 10:00 and 19:00 LT in the months from September to the following May. We compared the data of the Beijing, Mohe and Yakutsk digisondes and found that the LSTIDs are most likely to come from the northern auroral region. These LSTIDs may be induced by the atmospheric gravity waves (AGWs) and presented obvious seasonal and diurnal varying features, indicating that the thermospheric wind field has played an important role.

1. Introduction

Large-scale traveling ionospheric disturbances (LSTIDs) are generally considered as the manifestation of the atmospheric gravity waves (AGWs) in the ionosphere and are usually observed while the geomagnetic storms and auroral substorms occur [1,2,3,4,5]. Joule heating and Lorentz forces caused by the enhancement of the auroral electrojet or intense precipitation of charged particles are the main excitation sources of the LSTIDs, and they have typical horizontal velocities between 300 and 1000 m/s, horizontal wavelengths of more than 1000 km and period in the range of 30 min to 5 h. These periodic waves carry huge energy from the lower neutral atmosphere to the upper ionosphere and induce periodic electron density fluctuations in the ionosphere [4,6]. Almost all of the reported individual LSTID events occur during geomagnetic storms or auroral substorms (e.g., [7,8,9,10]). There are also a number of statistical analyses of the LSTIDs. The LSTIDs were found to occur in both the magnetically quiet and disturbed days according to the total electron content (TEC) data of the GPS observation network in Japan [11]. The statistical study of the TEC data in Europe shows that the amplitude of the LSTIDs has an obvious positive correlation with the AE index [12]. The occurrence rate of the LSTIDs in North America presents peaks at 12:00 and 19:00 LT [13]. The occurrence rate of the LSTIDs in China displays a maximum in winter and a minimum in summer [5]. In brief, there are many reports about the LSTIDs, and they come to different conclusions. However, very few reports specifically targeted the magnetically quiet LSTIDs.
A variety of equipment types, such as ionosondes [14,15,16], incoherent scatter radars [17,18,19,20], satellite beacons [21], HF Doppler radars [22] and airglow imagers [23,24,25], can be used to observe LSTIDs. LSTIDs span thousands of kilometers with very long wavelength. One piece of radio equipment can only observe one spot of the ionosphere. The airglow imager can observe a large area of the ionosphere but only operates at night. Nowadays, the high-density observations of large area ionosphere are mainly performed by the GPS receiver network. The TEC recorded by a GPS receiver is the total number of the electrons along the satellite–ground link. A vertical-incidence (VI) ionosonde operating in the high-frequency (HF) band can measure more parameters of the ionosphere, such as the maximum electron density and peak height of each ionospheric layer. Theoretically, the peak electron density variation gained through the swept-frequency ionogram can accurately reflect many kinds of ionospheric disturbance. Considering the cost of construction, it is impossible to set up a dense detection network composed of the VI ionosondes. However, the oblique-incidence (OI) ionosonde is a kind of low-cost system for network observation. We used the data of the OI ionosonde network in North China to study the mid-latitude LSTIDs. In Section 2, the detection network and the data analysis method are introduced. The typical examples of the LSTIDs and the statistical characteristics for two-year observation are presented in Section 3. In Section 4, the possible source and the reason for the diurnal and seasonal dependence of the recorded LSTIDs are discussed. A conclusion is made in the last section.

2. Typical Observation Cases

As shown in Figure 1, the OI ionosonde network in North China and three digisondes were used to investigate the large-scale oscillations in the ionosphere. The three digisondes are located in Beijing (40.3° N, 116.2° E), Mohe (53.5° N, 122.37° E) and Yakutsk (62.0° N, 129.6° E). Their geomagnetic latitudes are 30.2° N, 43.4° N, and 52.4° N, respectively, and are indicated as red triangles in Figure 1. They record ionograms from 0.5 to 12 MHz with 50 kHz step every 15 min. The ionospheric parameters as well as the electron density profile of the ionospheric F2 layer can be derived from the manually scaled trace of the ordinary waves by the SAO-Explorer [26]. The blue dots display the ionospheric observation points of the OI ionosonde network in North China.
Different from the traditional VI ionosonde, the transmitter and receiver of an OI ionosonde are separated, and it can also measure the ionospheric parameters such as the virtual and peak heights, as well as the critical frequency, like a traditional VI ionosonde [27,28]. As shown in Figure 2a, the ionospheric observation point is in the middle between the transmitter and receiver of an OI ionosonde. The OI ionosonde also operates in swept-frequency mode and records the echoes of each operating frequency. A typical swept-frequency OI ionogram is displayed in Figure 2b, and the echo trace of the ordinary waves is manually scaled and indicated by the gray curve. The distance between the transmitter and receiver is 1152.78 km. There is a certain relationship between the oblique-incidence ionogram and the vertical-incidence ionogram. The oblique echo trace is converted to the vertical ordinary trace by Martyn’s equivalent path theorem with the assumption of a spherically stratified ionosphere and no magnetic field [29]. The converted echo trace is shown in Figure 2c, and the blue curve displays a clear vertical echo trace of ionospheric F1 and F2 layers. The critical frequency of the F2 layer (foF2) used in the following is obtained by this echo trace conversion method, and the long period observation verified the accuracy of this method with the comparison of the VI observation [30,31]. Chen et al. [31] made an experimental comparison between oblique and vertical detections. They believed that the deviations of foF2 between the vertical and oblique observations never exceed 1 MHz and are generally less than 0.5 MHz.
In a network consisting of many OI ionosondes, the signal of one transmitter can be received by several receivers, and one receiver can receive the signal of several transmitters; thus, fewer radio systems can be used to obtain more ionospheric observation points. In North China, the OI ionosonde network composed of 5 transmitters and 20 receivers was constructed to observe the regional ionosphere with 99 observation points [32]. The blue circles in the map of Figure 2 display the observation points of the ionosonde network in North China. The observation area covers the latitudes between 32.34° N and 43.03° N, the longitudes between 111.35° E and 124.43° E and the geomagnetic latitudes between 22.40° N and 33.19° N. The ionosonde network records 99 OI ionograms every 30 min. Through the OI echo trace inversion method, we have estimated the foF2 at all the observation points of the ionosonde network. In order to eliminate the day-to-day variations and so highlight the periodic disturbances in the F2 layer, we estimate the difference foF2 (∆foF2) compared with the monthly value,
foF2 = (foF2 − foF2m)/foF2m
where foF2m represents the monthly mean of the foF2 during the magnetically quiet days.
The data of the OI ionosonde network from 1 September 2009 to 21 August 2011 are analyzed. Except for the maintenance time, the network has been working for 684 days. Excluding the impact of geomagnetic activities, there are 595 magnetically quiet days with the maximum Kp index ≤ 3. During these days, 157 cases of the LSTIDs from the north have been recorded. The observed LSTIDs mainly traveled to the south and southwest. We introduce two typical cases of them firstly.
One case of the LSTID occurring on 24 November 2010 is shown in Figure 3. The foF2 depletion firstly appeared at 14:30 LT in the top right corner of Figure 3b and then extended southwestward to cover all the observation area at 16:30 LT. At 17:00 LT, the foF2 enhancement emerged in the top right corner of Figure 3g and also moved to the southwest. At 20:00 LT in Figure 3m, the regional enhancement had moved out of the observation region. Soon after, the second circle of the foF2 enhancement occurred in the northeast in Figure 3n. The regional enhancement traveled southwest through the observation region and moved out at 23:00 LT. The other typical case was recorded on 22 January 2011 and is displayed in Figure 4. The regional foF2 enhancement emerged in the north of the map in Figure 4b at 8:30 LT and then moved southward. At 10:30 LT in Figure 4f, it propagated out of the observation area. The second and third circles of the southward moving foF2 enhancement appeared at 12:30 LT in Figure 4j and at 16:00 LT in Figure 4q, respectively. Both cases present the obvious periodic oscillations traveling across the North China ionosphere.
The latitude and longitude dependences of the two cases of ionospheric oscillations are drawn in the top two rows of Figure 5. The oscillations recorded on 24 November 2010 and 22 January 2011 are displayed in the left and right columns of the figure, respectively. The first row shows the latitude–time dependence of the ∆foF2 values recorded by the 99 observation points. The disturbances of several periods are clearly displayed. They lasted for more than 10 h with periods of 3.9 h in Figure 5a and 4.3 h in Figure 5f. The wave peaks and troughs slanting to the right in both plots indicate the southward moving tendency. The longitude–time dependence of the ∆foF2 values is displayed in the second row of Figure 5, and the obvious fluctuation features are also present in these two plots. The peaks and troughs in Figure 5b move to the smaller longitudes, indicating the westward zonal velocity. In Figure 5g, the peaks and troughs at different longitudes appeared almost simultaneously, implying that the waves were propagating to almost due south. According to the slopes of wave peak and trough in the latitude–time and longitude–time plots, the horizontal speed and the azimuth of the oscillation in the left column are estimated as 169.03 ± 13.92 m/s and 221.57 ± 0.56°, and the same parameters in the right column are 218.24 ± 12.88 m/s and 182.18 ± 0.31°. The computing method for the wave speed and azimuth is provided in the Appendix A. This algorithm can accurately calculate the parameters we need.
To further study the characteristics of the two LSTIDs, the observations of the three digisondes at different latitudes are analyzed. The third row of Figure 5 displays the electron density profiles estimated from the ionograms of the Beijing digisonde. The black curves display the variations of the peak height of the F2 layer (hmF2). The variations of the plasma frequency at the heights from 190 to 230 km with 10 km step are presented in the fourth row of Figure 5. The displayed plasma frequency curves have passed through the band-pass filter of 1–6 h. The periodic fluctuations appear on all the displayed heights and they are consistent with the wave phase in the top two rows. The solid and hollow circles indicate the data below and above the hmF2, respectively. We connect the wave peaks and troughs at different heights and find that the wave phase velocities are downward. The echoes of an ionosonde are reflected below hmF2, so the solid circles in Figure 5d,i are applied to estimate the vertical phase velocities. According to the slope of the connecting lines, the vertical phase velocities of the two LSTIDs are estimated as 16.81 ± 4.69 and 32.74 ± 8.93 m/s downward, respectively. Usually, most of the large-scale disturbances in the ionosphere stem from the lower atmosphere, and the downward phase velocities indicate that they are very likely the upward gravity waves [33,34,35]. The filtered foF2 curves recorded by the digisondes in Beijing, Mohe and Yakutsk are shown in the bottom row of Figure 5. The band-pass filter of 1–6 h is also applied. The wave peaks and troughs are connected by the slope lines, which indicate the southward propagating velocities of the LSTIDs. Due to the very long distance between the three digisondes, the traveling speed of the waves cannot be accurately estimated, but it clearly illuminates that the recorded LSTIDs come from the auroral region north of Yakutsk.

3. Statistical Analysis

We develop four criteria to accurately recognize the LSTIDs from other ionospheric disturbances: (1) The amplitude of the enhancement or depletion exceeds the diurnal average by at least 5%. (2) At least 60% of the ionospheric observation points of the ionosonde network as well as the Beijing digisonde have recorded the same periodic disturbance. (3) The duration of the wave peak and trough exceeds 30 min. (4) More than two circles of the disturbances are recorded. In the magnetically quiet days from 1 September 2009 to 21 August 2011, 157 cases of the LSTIDs from the north have been recorded by the OI ionosonde network. The total occurrence rate of the quiet-time LSTIDs is 26.39%, which is defined as a fraction of the days the quiet-time LSTIDs were observed against the total number of magnetically quiet days. The four seasons in North China are defined as spring (March–April), summer (May–August), autumn (September–October), and Winter (November–February). The occurrence number and occurrence rate of the quiet-time LSTIDs in four seasons are presented in Table 1. The LSTIDs mainly occur in winter and can be observed almost every two days. In summer, only two cases have been recorded. The occurrence rate in autumn is 10.7% more than that in spring. This is consistent with the previous observations of Tang et al. [36]. The periods of the LSTIDs recorded in the four seasons all tend to have Gaussian distribution between 2–5 h, as shown in Figure 6a, and they are centered at 3.44 h. Figure 6b shows the distribution of the horizontal velocities of the southward and southwestward traveling disturbances, respectively. The velocities of the southwestward waves vary mainly between 106.67 and 256.01 m/s, and those of the southward waves are mainly between 153.12 and 312.52 m/s. The median azimuths and the standard deviations of the southward and southwestward traveling LSTIDs are 181.58 ± 14.43° and 224.98 ± 17.43°, respectively.
The occurrence rates of the LSTIDs varying with local time in different months are shown in Figure 7 and present obvious diurnal and seasonal dependence. The occurrence rate is defined as a fraction of the number of days of the quiet-time LSTIDs observed against the total number of magnetically quiet days of a given month, the same method used by Otsuka et al. [37]. The mid-latitude LSTIDs appear mainly between 9:00 and 19:00 LT from September to March. The statistical study of Tsugawa et al. [11] also showed that under the quiet conditions of Kp ≤ 3, the occurrence number of the LSTIDs was much larger in the daytime than in the nighttime.
Two-year azimuth variations of the mid-latitude LSTIDs are presented as the blue dots in Figure 8. Except for a few days of downtime as indicated by the gray vertical bars, the ionosonde network was well operated from 1 September 2009 to 21 August 2011. As shown in Figure 8, the recorded LSTIDs come from about three fixed directions. Only one southeastward LSTID emerged on 17 October 2010. There is an obvious boundary between the southwestward and southward traveling LSTIDs as shown by the red horizontal line at 200° azimuth. There are 46 cases above the red horizontal line traveling southwestward with the azimuth between 207.55 and 242.40°. Their average azimuth is 218.66° and is indicated by the dotted blue horizontal line above. There are 110 cases below the red horizontal line propagating southward with the azimuth between 167.14 and 196.01°. The dotted blue horizontal line below is used to indicate their average azimuth of 182.95°. It is very likely that the two kinds of LSTIDs stemmed from two different sources. In the studies of the ionospheric disturbances induced by earthquakes, thunderstorms and so on [35,38,39,40], the waves traveled from the source to all directions, and the measured wave direction depends on the orientation of the wave source relative to the observer. Both the southwestward and southward LSTIDs present the same seasonal characteristics. They appeared mainly from September to the following May and disappeared from June to August. It is very likely that the same background factor dominates the propagation of the LSTIDs in two different directions.

4. Discussion

Almost all the case studies of the LSTIDs focus on the disturbances associated with storms (e.g., [10,41,42,43]). Some statistical studies also considered that LSTIDs are strongly linked to the geomagnetic activity at high latitudes (e.g., [13,36,44]). However, there are a few investigations about the LSTIDs occurring during the magnetically quiet days, and much fewer investigations concern the LSTIDs with periods larger 2 h (e.g., [5,45]). We attempted to study the LSTIDs with periods larger than 1 h using the OI ionosonde network in North China. According to the two-year observations, we have found interesting LSTIDs with periods concentrated around 3.44 h. Most of the LSTIDs came from north and northeast, and their downward phase velocities between 190 and 230 km altitudes indicated that they may travel upward from the lower atmosphere and be related to AGWs. Compared with the observations of the three digisondes in different latitudes from high to mid-latitude, we find that both the southward and southwestward propagating LSTIDs stemmed from the north auroral region. However, these quiet-time LSTIDs have a larger period than the disturbed-time LSTIDs reported by Tsugawa et al. [11]. Quiet-time LSTIDs are often considered to be excited by the AGWs moving equatorward from a source in the dawn sector of the auroral zone [46]. The momentum and energy disturbed by the AWG are transmitted to the charged particles through the collision between neutral particles and charged particles. Thus, the periodic disturbances travel into the ionosphere from the neutral atmosphere to generate the TID along its propagation path.
Most AGWs are generated in the lower atmosphere and travel upward. Previous studies by many scholars have shown that the propagation of AGWs is very sensitive to the filtering effect of the background wind. Therefore, the propagation direction of gravity waves exhibits obvious seasonal characteristics. Medeiros et al. [47] investigated the dominant quasi-monochromatic gravity waves in the upper mesosphere from September 1998 to October 1999 using the all-sky imager at Cachoeira Paulista, Brazil (23° S, 45° W), and found the waves exhibiting an obvious seasonal dependence on the horizontal propagation direction, propagating toward the southeast during the summer months and toward the northwest during the winter, which were attributed to a strong filtering of the waves in the middle atmosphere by the stratospheric winds. Takeo et al. [48] analyzed the 557.7 nm airglow images obtained at Shigaraki MU Observatory (34.8° N, 136.1° E) and found the propagation direction of gravity waves presented clear seasonal variation in the mesopause region. The propagation direction of the gravity waves in the mesopause region was mainly controlled by the wind filtering. The global gravity wave variations in the thermosphere were investigated by Miyoshi et al. [49]. The gravity waves were found propagating upward from the lower atmosphere, and the wave energy at the thermospheric altitudes reveals obvious seasonal dependence, indicating the filtering effect of the background wind. The TIDs are suggested as the plasma density manifestations of the gravity waves at ionospheric F-layer heights [50]. When a gravity wave propagates upward in the atmosphere, its wave amplitude becomes larger and larger due to the lower and lower atmospheric density with the altitude increasing. Most of the upward gravity waves will break before they pass through the mesopause, and some of them can reach the heights of the thermosphere and modulate the plasma to form the periodic ionospheric disturbances. The downward phase velocities of the recorded LSTIDs indicate that they may be induced by the AGWs traveling upward. The statistics of the quiet-time LSTIDs from both north and northeast present similar seasonal dependence, namely the highest occurrence rate in winter and the lowest in summer, which implies the possible wind filtering effect in the atmosphere.
According to the law of energy conservation, the amplitude of a gravity wave increases as the atmospheric density decreases with height. When reaching the MLT region, the wave amplitude is large enough to generate instabilities, which would lead to wave breaking when the temperature lapse rate of the waves plus the mean flow becomes super-adiabatic [51,52,53,54,55]. Therefore, the background wind field at these heights and above becomes very important for the upward traveling waves. In the thermosphere, the neutral winds blow poleward during daytime and equatorward during nighttime due to the pressure gradients caused by the solar EUV heating of the thermosphere. Regarding the seasonal variation of the thermospheric winds, the wind blows from the summer hemisphere to the winter hemisphere. Consequently, the poleward wind during the daytime is stronger in winter than in summer. In Figure 7 and Figure 8, it can be seen that the diurnal and seasonal dependences of the LSTIDs are closely related to the thermospheric meridian wind field. The 2–5 h LSTIDs can be observed at mid-latitudes when they are traveling in the opposite direction of the mean meridian winds. When the mean meridian winds turn equatorward, both the southward and southwestward LSTIDs disappear at mid-latitudes. As is well known, the wind field in the same direction of the phase velocity of an AGW will stretch the AGW rays in the propagation direction and increase the energy attenuation of the waves in the atmosphere, thus hindering the AGW propagation. Conversely, when the AGW encounters a reverse wind field, the AGW rays will shrink to shorten the propagation path of the AGW, and thus the energy attenuation of the wave will be reduced to form a favorable environment for the wave propagation [56,57]. According to the dispersion relation for gravity waves, the vertical wavelength of a gravity wave becomes small when the background wind blows in the same direction as the gravity wave propagation, and thus the gravity waves tend to be dissipated due to high viscosity in the thermosphere. As a result, the gravity waves in winter can survive. This suggests that the atmospheric background wind field could play an important role in the LSTIDs’ appearance.
The seasonal and local time variation of the quiet-time LSTIDs is similar to that of the medium-scale traveling ionospheric disturbances (MSTIDs) propagating equatorward (e.g., [37,58,59]). Both of them occur frequently during the daytime in winter. The equatorward traveling MSTIDs are also suggested to be induced by AWGs. We speculate that the background wind field has a significant influence on both the LSTIDs and MSTIDs, and so there is a similarity between them when they are traveling in the same direction.

5. Conclusions

Aside from the GPS network, the dense OI ionosonde network offers another possibility to continuously observe the horizontal propagation of TIDs by both day and night. It is an economical solution using multi-input and multi-output technology. The inversion parameters of the network can present the propagation characteristics of the oscillations on a horizontal plane through the peak plasma density variations at each observation point. Based on the observations of the OI ionosonde network in North China and the digisondes in Beijing, Mohe and Yakutsk, a statistical study of the quiet-time LSTIDs was conducted. An interesting phenomenon of the mid-latitude LSTIDs of quasi-3-h period has been reported. It highlights the observation particularity of the OI ionosonde network, and the long-term observations of this network will deepen our understanding of the atmosphere–ionosphere coupling as well as the influence of momentum from polar regions on the middle and low latitudes.

Author Contributions

Conceptualization, R.Z. and G.C.; data curation, G.C., Y.L. and W.G.; formal analysis, R.Z.; funding acquisition, G.C.; investigation, R.Z. and G.C.; methodology, G.C.; project administration, G.C.; resources, G.C.; software, R.Z.; supervision, R.Z., G.C. and Z.H.; validation, R.Z., G.C. and S.Z.; visualization, R.Z. and M.Z.; writing—original draft, R.Z.; writing—review and editing, R.Z. and G.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 41722404.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The Beijing and Mohe digisonde data can be obtained from the Chinese Meridian Project (http://159.226.22.74/ (accessed on 28 December 2021)). The data of the oblique-incidence ionosonde network can be obtained from the China seismo-ionospheric ground-based monitoring network (http://www.sign.ac.cn/ (accessed on 8 April 2021)). The data of the Yakutsk digisonde can be obtained from the Lowell GIRO Data Center (http://giro.uml.edu/ (accessed on 28 December 2021)). The Kp indexes can be obtained from the World Data Center for Geomagnetism Web interface (http://wdc.kugi.kyoto-u.ac.jp (accessed on 28 December 2021)).

Acknowledgments

The authors acknowledge the use of the Beijing and Mohe digisonde data from the Chinese Meridian Project (http://159.226.22.74/ (accessed on 28 December 2021)). The data of the oblique-incidence ionosonde network can be obtained from the China seismo-ionospheric ground-based monitoring network (http://www.sign.ac.cn/ (accessed on 8 April 2021)). The data of the Yakutsk digisonde come from the Lowell GIRO Data Center (http://giro.uml.edu/ (accessed on 28 December 2021)). The Kp indexes are obtained from the World Data Center for Geomagnetism Web interface (http://wdc.kugi.kyoto-u.ac.jp (accessed on 28 December 2021)). The authors appreciate the patience and helpful suggestions of the two reviewers.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Computing method for the wave speed and azimuth.
Given a periodic disturbance traveling with the speed Vo and azimuth α as indicated in Figure A1, the zonal velocity Vz and meridian velocity Vm of the disturbance are
Vz = Vo sin α
Vm = Vo cos α
Figure A1. Schematic diagram of the propagating parameters of the periodic disturbances in the horizontal plane. The azimuth α is clockwise from north. The letters N, E, S, and W are used to indicate the directions of north, east, south and west, respectively.
Figure A1. Schematic diagram of the propagating parameters of the periodic disturbances in the horizontal plane. The azimuth α is clockwise from north. The letters N, E, S, and W are used to indicate the directions of north, east, south and west, respectively.
Sensors 22 00233 g0a1
The projection speeds of the disturbance on the horizontal and vertical axes are Vx and Vy. It must be pointed out that a disturbance is not a point target, and thus VyVm and VxVz. Vx/Vy can be estimated by the slope of the wavefront in the longitude/latitude–time plot. When Vx and Vy are known, the wave speed Vo and wave azimuth α can be calculated as follows:
V x = V z + V m cot α = V z + V m 2 V z = V z 2 + V m 2 V z = V o 2 V z = V o sin α
V y = V m + V z tan α = V m + V z 2 V m = V m 2 + V z 2 V m = V o 2 V m = V o cos α
tan α = V z V m
By the Equations (A3) and (A4), we know
V x sin α = V y cos α tan α = V y V x
So,
α = arctan V y V x
V o = V y cos α = V x sin α
Given that the slope of the wave peak and wave trough in the latitude–time and longitude–time plots are lm and lz (°/s), respectively:
V y = l m R o
V x = l z R o cos 37 °
where Ro represents the distance of 1° latitude (or 1° longitude at the equator) on the Earth’s surface and the mean latitude of the ionosonde network is 37°. So, the wave azimuth α and the wave horizontal speed Vo can be also expressed as
α = arctan ( l m l z cos 37 ° )
V o = l m R o cos α

References

  1. Ding, F.; Wan, W.X.; Ning, B.Q.; Zhao, B.Q.; Li, Q.; Zhang, R.; Xiong, B.; Song, Q. Two-dimensional imaging of large-scale traveling ionospheric disturbances over China based on GPS data. J. Geophys. Res. Space Phys. 2012, 117, A08318. [Google Scholar] [CrossRef]
  2. Hocke, K.; Schlegel, K. A review of atmospheric gravity waves and travelling ionospheric disturbances: 1982–1995. Ann. Geophys. 1996, 14, 917–940. [Google Scholar] [CrossRef]
  3. Hunsucker, R.D. Atmospheric gravity waves generated in the high-latitude ionosphere: A review. Rev. Geophys. Space Phys. 1982, 20, 293–315. [Google Scholar] [CrossRef]
  4. Lei, J.H.; Burns, A.G.; Tsugawa, T.; Wang, W.B.; Solomon, S.C.; Wiltberger, M. Observations and simulations of quasiperiodic ionospheric oscillations and large-scale traveling ionospheric disturbances during the December 2006 geomagnetic storm. J. Geophys. Res. Space Phys. 2008, 113, A06310. [Google Scholar] [CrossRef] [Green Version]
  5. Song, Q.; Ding, F.; Wan, W.X.; Ning, B.Q.; Liu, L.B.; Zhao, B.Q.; Li, Q.; Zhang, R. Statistical study of large-scale traveling ionospheric disturbances generated by the solar terminator over China. J. Geophys. Res. Space Phys. 2013, 118, 4583–4593. [Google Scholar] [CrossRef]
  6. Fujiwara, H.; Miyoshi, Y. Characteristics of the large-scale traveling atmospheric disturbances during geomagnetically quiet and disturbed periods simulated by a whole atmosphere general circulation model. Geophys. Res. Lett. 2006, 33, L20108. [Google Scholar] [CrossRef]
  7. Figueiredo, C.A.O.B.; Wrasse, C.M.; Takahashi, H.; Otsuka, Y.; Shiokawa, K.; Barros, D. Large-scale traveling ionospheric disturbances observed by GPS dTEC maps over North and South America on Saint Patrick’s Day storm in 2015. J. Geophys. Res. Space Phys. 2017, 122, 4755–4763. [Google Scholar] [CrossRef]
  8. Horvath, I.; Lovell, B.C. Large-scale traveling ionospheric disturbances impacting equatorial ionization anomaly development in the local morning hours of the Halloween Superstorms on 29–30 October 2003. J. Geophys. Res. Space Phys. 2010, 115, A04302. [Google Scholar] [CrossRef] [Green Version]
  9. Idrus, I.I.; Abdullah, M.; Hasbi, A.M.; Husin, A.; Yatim, B. Large-scale traveling ionospheric disturbances observed using GPS receivers over high-latitude and equatorial regions. J. Atmos. Sol.-Terr. Phys. 2013, 102, 321–328. [Google Scholar] [CrossRef]
  10. Tsugawa, T.; Saito, A.; Otsuka, Y.; Yamamoto, M. Damping of large-scale traveling ionospheric disturbances detected with GPS networks during the geomagnetic storm. J. Geophys. Res. 2003, 108, 1127. [Google Scholar] [CrossRef]
  11. Tsugawa, T.; Saito, A.; Otsuka, Y. A statistical study of large-scale traveling ionospheric disturbances using the GPS network in Japan. J. Geophys. Res. 2004, 109, A06302. [Google Scholar] [CrossRef]
  12. Borries, C.; Jakowski, N.; Kauristie, K.; Amm, O.; Mielich, J.; Kouba, D. On the dynamics of large-scale traveling ionospheric disturbances over Europe on 20 November 2003. J. Geophys. Res. Space Phys. 2017, 122, 1199–1211. [Google Scholar] [CrossRef] [Green Version]
  13. Ding, F.; Wan, W.X.; Liu, L.B.; Afraimovich, E.L.; Voeykov, S.V.; Perevalova, N.P. A statistical study of large-scale traveling ionospheric disturbances observed by GPS TEC during major magnetic storms over the years 2003–2005. J. Geophys. Res. Space Phys. 2008, 113, A00A01. [Google Scholar] [CrossRef]
  14. Hajkowicz, L.A. A global study of large scale travelling ionospheric disturbances (TIDS) following a step-like onset of auroral substorms in both hemispheres. Planet. Space Sci. 1990, 38, 913–923. [Google Scholar] [CrossRef]
  15. Hajkowicz, L.A. Global onset and propagation of large scale travelling ionospheric disturbances as a result of the great storm of 13 March 1989. Planet. Space Sci. 1991, 39, 583–593. [Google Scholar] [CrossRef]
  16. Maeda, S.; Handa, S. Transmission of large-scale TIDs in the ionospheric F2-region. J. Atmos. Terr. Phys. 1980, 42, 853–859. [Google Scholar] [CrossRef]
  17. Nicolls, M.J.; Kelley, M.C. Strong evidence for gravity wave seeding of an ionospheric plasma instability. Geophys. Res. Lett. 2005, 32, L05108. [Google Scholar] [CrossRef]
  18. Panasenko, S.V.; Goncharenko, L.P.; Erickson, P.J.; Aksonova, K.D.; Domnin, I.F. Traveling ionospheric disturbances observed by Kharkiv and Millstone Hill incoherent scatter radars near vernal equinox and summer solstice. J. Atmos. Sol.-Terr. Phys. 2018, 172, 10–23. [Google Scholar] [CrossRef]
  19. Sheen, D.R.; Liu, C.H. Modeling F region gravity waves observed during the WAGS campaign: 1. Special event. Radio Sci. 1988, 23, 879–893. [Google Scholar] [CrossRef]
  20. Vadas, S.L.; Nicolls, M.J. Temporal evolution of neutral, thermospheric winds and plasma response using PFISR measurements of gravity waves. J. Atmos. Sol.-Terr. Phys. 2009, 71, 744–770. [Google Scholar] [CrossRef]
  21. Jacobson, A.R.; Carlos, R.C.; Massey, R.S.; Wu, G.H. Observations of traveling ionospheric disturbances with a satellite-beacon radio interferometer: Seasonal and local time behavior. J. Geophys. Res. 1995, 100, 1653–1665. [Google Scholar] [CrossRef]
  22. Georges, T.M. HF Doppler studies of traveling ionospheric disturbances. J. Atmos. Terr. Phys. 1968, 30, 735–746. [Google Scholar] [CrossRef]
  23. Mendillo, M.; Baumgardner, J.; Nottingham, D.; Aarons, J.; Reinisch, B.; Scali, J.; Kelley, M. Investigations of thermospheric-ionospheric dynamics with 6300-Å images from the Arecibo Observatory. J. Geophys. Res. Space Phys. 1997, 102, 7331–7343. [Google Scholar] [CrossRef] [Green Version]
  24. Shiokawa, K. Thermospheric wind during a storm-time large-scale traveling ionospheric disturbance. J. Geophys. Res. 2003, 108, 1423. [Google Scholar] [CrossRef] [Green Version]
  25. Shiokawa, K. Geomagnetic conjugate observation of nighttime medium-scale and large-scale traveling ionospheric disturbances: FRONT3 campaign. J. Geophys. Res. 2005, 110, A05303. [Google Scholar] [CrossRef]
  26. Huang, X.Q.; Reinisch, B.W. Vertical electron density profiles from the Digisonde network. Adv. Space Res. 1996, 18, 121–129. [Google Scholar] [CrossRef]
  27. Chen, G.; Zhao, Z.Y.; Zhu, G.Q.; Shi, S.Z. The Wuhan Ionospheric Sounding Systems. IEEE Geosci. Remote Sens. Lett. 2009, 6, 748–751. [Google Scholar] [CrossRef]
  28. Huang, X.Q.; Reinisch, B.W.; Kuklinski, W.S. Mid-point electron density profiles from oblique ionograms. Ann. Geophys. 1996, 39, 757–761. [Google Scholar] [CrossRef]
  29. Reilly, M.H.; Kolesar, J.D. A method for real height analysis of oblique ionograms. Radio Sci. 1989, 24, 575–583. [Google Scholar] [CrossRef]
  30. Phanivong, B.; Chen, J.; Dyson, P.L.; Bennett, J.A. Inversion of oblique ionograms including the Earth’s magnetic field. J. Atmos. Terr. Phys. 1995, 57, 1715–1721. [Google Scholar] [CrossRef]
  31. Chen, G.; Zhao, Z.Y.; Zhang, Y.N.; Yang, G.B.; Zhou, C.; Huang, S.; Li, T.; Li, N.; Sun, H.Q. Application of the oblique ionogram as vertical ionogram. Sci. China Technol. Sci. 2012, 55, 1240–1244. [Google Scholar] [CrossRef]
  32. Chen, G.; Wang, J.; Zhang, S.D.; Deng, Z.X.; Zhong, D.K.; Wu, C.; Jin, H.; Li, Y.X. Opposite latitudinal dependence of the premidnight and postmidnight oscillations in the electron density of midlatitude F layer. J. Geophys. Res. Space Phys. 2018, 123, 796–807. [Google Scholar] [CrossRef]
  33. Hines, C.O. Gravity waves in the atmosphere. Nature 1972, 239, 73–78. [Google Scholar] [CrossRef]
  34. Chen, G.; Wu, C.; Zhao, Z.Y.; Zhong, D.K.; Qi, H.; Jin, H. Daytime E region field-aligned irregularities observed during a solar eclipse. J. Geophys. Res. Space Phys. 2014, 119, 10633–10640. [Google Scholar] [CrossRef]
  35. Chen, G.; Wang, J.; Wu, C.; Huang, X.Q.; Zhong, D.K.; Qi, H.; Huang, L.; Li, Y.X. Multisite Remote Sensing for Tsunami-Induced Waves. IEEE Trans. Geosci. Remote Sens. 2016, 54, 7177–7184. [Google Scholar] [CrossRef]
  36. Tang, Q.L.; Wan, W.X.; Ning, B.Q.; Yuan, H. Properties of large-scale TIDs observed in central China. Sci. China Ser. A Math. 2002, 45, 156–160. [Google Scholar] [CrossRef]
  37. Otsuka, Y.; Suzuki, K.; Nakagawa, S.; Nishioka, M.; Shiokawa, K.; Tsugawa, T. GPS observations of medium-scale traveling ionospheric disturbances over Europe. Ann. Geophys. 2013, 31, 163–172. [Google Scholar] [CrossRef]
  38. Jin, S.G.; Jin, R.; Li, J.H. Pattern and evolution of seismo-ionospheric disturbances following the 2011 Tohoku earthquakes from GPS observations. J. Geophys. Res. Space Phys. 2014, 119, 7914–7927. [Google Scholar] [CrossRef] [Green Version]
  39. Liu, J.Y.; Chen, C.H.; Lin, C.H.; Tsai, H.F.; Chen, C.H.; Kamogawa, M. Ionospheric disturbances triggered by the 11 March 2011M9.0 Tohoku earthquake. J. Geophys. Res. Space Phys. 2011, 116, A06319. [Google Scholar] [CrossRef]
  40. Xu, J.Y.; Li, Q.Z.; Yue, J.; Hoffmann, L.; Straka, W.C.; Wang, C.M.; Liu, M.H.; Yuan, W.; Han, S.; Miller, S.D.; et al. Concentric gravity waves over northern China observed by an airglow imager network and satellites. J. Geophys. Res. Atmos. 2015, 120, 11058–11078. [Google Scholar] [CrossRef] [Green Version]
  41. Ding, F.; Wan, W.; Ning, B.; Wang, M. Large-scale traveling ionospheric disturbances observed by GPS total electron content during the magnetic storm of 29–30 October 2003. J. Geophys. Res. Space Phys. 2007, 112, A06309. [Google Scholar] [CrossRef]
  42. Ding, F.; Wan, W.; Ning, B.; Zhao, B.; Li, Q.; Wang, Y.; Hu, L.; Zhang, R.; Xiong, B. Observations of poleward-propagating large-scale traveling ionospheric disturbances in southern China. Ann. Geophys. 2013, 31, 377–385. [Google Scholar] [CrossRef] [Green Version]
  43. Zakharenkova, I.; Astafyeva, E.; Cherniak, I. GPS and GLONASS observations of large-scale traveling ionospheric disturbances during the 2015 St. Patrick’s Day storm. J. Geophys. Res. Space Phys. 2016, 121, 12138–12156. [Google Scholar] [CrossRef]
  44. Katamzi, Z.T.; Smith, N.D.; Mitchell, C.N.; Spalla, P.; Materassi, M. Statistical analysis of travelling ionospheric disturbances using TEC observations from geostationary satellites. J. Atmos. Sol.-Terr. Phys. 2012, 74, 64–80. [Google Scholar] [CrossRef]
  45. Ding, F.; Wan, W.X.; Li, Q.; Zhang, R.; Song, Q.; Ning, B.Q.; Liu, L.B.; Zhao, B.Q.; Xiong, B. Comparative climatological study of large-scale traveling ionospheric disturbances over North America and China in 2011–2012. J. Geophys. Res. Space Phys. 2014, 119, 519–529. [Google Scholar] [CrossRef] [Green Version]
  46. Millward, G.H.; Moffett, R.J.; Quegan, S.; Fuller-Rowell, T.J. Effects of an atmospheric gravity wave on the midlatitude ionospheric F layer. J. Geophys. Res. Space Phys. 1993, 98, 19173–19179. [Google Scholar] [CrossRef]
  47. Medeiros, A.F.; Taylor, M.J.; Takahashi, H.; Batista, P.P.; Gobbi, D. An investigation of gravity wave activity in the low-latitude upper mesosphere: Propagation direction and wind filtering. J. Geophys. Res. 2003, 108, 4411. [Google Scholar] [CrossRef] [Green Version]
  48. Takeo, D.; Shiokawa, K.; Fujinami, H.; Otsuka, Y.; Matsuda, T.S.; Ejiri, M.K.; Nakamura, T.; Yamamoto, M. Sixteen year variation of horizontal phase velocity and propagation direction of mesospheric and thermospheric waves in airglow images at Shigaraki, Japan. J. Geophys. Res. Space Phys. 2017, 122, 8770–8780. [Google Scholar] [CrossRef]
  49. Miyoshi, Y.; Fujiwara, H.; Jin, H.; Shinagawa, H. A global view of gravity waves in the thermosphere simulated by a general circulation model. J. Geophys. Res. Space Phys. 2014, 119, 5807–5820. [Google Scholar] [CrossRef]
  50. Hines, C.O. Internal atmospheric gravity waves at ionospheric heights. Can. J. Phys. 1960, 38, 1441–1481. [Google Scholar] [CrossRef]
  51. Lindzen, R.S. Thermally driven diurnal tide in the atmosphere. Q. J. R. Meteorol. Soc. 1967, 93, 18–42. [Google Scholar] [CrossRef]
  52. Lindzen, R.S. Lower atmospheric energy sources for the upper atmosphere. In Meteorological Monographs, Meteorological Investigations of the Upper Atmosphere; American Meteorological Society: Boston, MA, USA, 1968; pp. 37–46. [Google Scholar] [CrossRef]
  53. Hodges, R.R., Jr. Generation of turbulence in the upper atmosphere by internal gravity waves. J. Geophys. Res. 1967, 72, 3455–3458. [Google Scholar] [CrossRef]
  54. Hodges, R.R., Jr. Eddy diffusion coefficients due to instabilities in internal gravity waves. J. Geophys. Res. 1969, 74, 4087–4090. [Google Scholar] [CrossRef]
  55. Pitteway, M.L.V.; Hines, C.O. The viscous damping of atmospheric gravity waves. Can. J. Phys. 1963, 41, 1935–1948. [Google Scholar] [CrossRef]
  56. Ding, F.; Wan, W.; Yuan, H. The influence of background winds and attenuation on the propagation of atmospheric gravity waves. J. Atmos. Sol.-Terr. Phys. 2003, 65, 857–869. [Google Scholar] [CrossRef]
  57. Yeh, K.C.; Webb, H.D.; Cowling, D.H. Evidence of directional filtering of travelling ionospheric disturbance. Nat. Phys. Sci. 1972, 235, 131–132. [Google Scholar] [CrossRef]
  58. Kotake, N.; Otsuka, Y.; Ogawa, T.; Tsugawa, T.; Saito, A. Statistical study of medium-scale traveling ionospheric disturbances observed with the GPS networks in Southern California. Earth Planets Space 2007, 59, 95–102. [Google Scholar] [CrossRef] [Green Version]
  59. Otsuka, Y.; Shinbori, A.; Tsugawa, T.; Nishioka, M. Solar activity dependence of medium-scale traveling ionospheric disturbances using GPS receivers in Japan. Earth Planets Space 2021, 73, 22. [Google Scholar] [CrossRef]
Figure 1. Geographical locations of the observation points in Eastern Asia. The 99 blue circles indicate the ionospheric observation points of the oblique-incidence ionosonde network. The three red triangles show the locations of the Beijing, Mohe and Yakutsk digisondes.
Figure 1. Geographical locations of the observation points in Eastern Asia. The 99 blue circles indicate the ionospheric observation points of the oblique-incidence ionosonde network. The three red triangles show the locations of the Beijing, Mohe and Yakutsk digisondes.
Sensors 22 00233 g001
Figure 2. (a) Schematic diagram for ionospheric oblique-incidence detection. (b) A typical oblique-incidence ionogram with ordinary and extraordinary waves. The gray curve is the echo trace of the F1 and F2 layer ordinary wave. (c) The converted vertical-incidence echo trace through the gray curve in plot (b).
Figure 2. (a) Schematic diagram for ionospheric oblique-incidence detection. (b) A typical oblique-incidence ionogram with ordinary and extraordinary waves. The gray curve is the echo trace of the F1 and F2 layer ordinary wave. (c) The converted vertical-incidence echo trace through the gray curve in plot (b).
Sensors 22 00233 g002
Figure 3. Maps showing the ΔfoF2 variations at all the observation points of the ionosonde network from 14:00 to 23:30 LT with half an hour step on 24 November 2010.
Figure 3. Maps showing the ΔfoF2 variations at all the observation points of the ionosonde network from 14:00 to 23:30 LT with half an hour step on 24 November 2010.
Sensors 22 00233 g003
Figure 4. Maps showing the ΔfoF2 variations at all the observation points of the ionosonde network from 8:00 to 17:30 LT with half an hour step on 22 January 2011.
Figure 4. Maps showing the ΔfoF2 variations at all the observation points of the ionosonde network from 8:00 to 17:30 LT with half an hour step on 22 January 2011.
Sensors 22 00233 g004
Figure 5. Two typical examples of the LSTIDs that occurred between 14:00 and 24:00 on 24 November 2010 (ae) and between 8:00 and 18:00 on 22 January 2011 (fj). The first row shows the latitude dependence of the LSTIDs. The second row shows the longitude dependence of the LSTIDs. The third row shows the ionospheric electron density profiles recorded by the Beijing digisonde, and the black curves indicate the peak height of the F2 layer (hmF2). The fourth row shows the filtered plasma frequency variations at five selected heights. The red and blue circles represent the peak and trough of the observed fluctuations, and the solid lines connect the circles at different heights. The hollow circles indicate the data above the peak height of the F2 layer. The bottom row displays the filtered F2-layer critical frequency (foF2) recorded by the Yakutsk (YAKT), Mohe (MH) and Beijing (BJ) digisondes. The red and blue circles represent the peak and trough on the foF2 curves, and the solid lines connect the circles to show the wave propagating tendency.
Figure 5. Two typical examples of the LSTIDs that occurred between 14:00 and 24:00 on 24 November 2010 (ae) and between 8:00 and 18:00 on 22 January 2011 (fj). The first row shows the latitude dependence of the LSTIDs. The second row shows the longitude dependence of the LSTIDs. The third row shows the ionospheric electron density profiles recorded by the Beijing digisonde, and the black curves indicate the peak height of the F2 layer (hmF2). The fourth row shows the filtered plasma frequency variations at five selected heights. The red and blue circles represent the peak and trough of the observed fluctuations, and the solid lines connect the circles at different heights. The hollow circles indicate the data above the peak height of the F2 layer. The bottom row displays the filtered F2-layer critical frequency (foF2) recorded by the Yakutsk (YAKT), Mohe (MH) and Beijing (BJ) digisondes. The red and blue circles represent the peak and trough on the foF2 curves, and the solid lines connect the circles to show the wave propagating tendency.
Sensors 22 00233 g005
Figure 6. (a) Period distribution of the LSTIDs from the north in different seasons between 1 September 2009 and 21 August 2011. The solid curves are the Gaussian fitting of the histograms. (b) Horizontal velocity distribution of the LSTIDs propagating southward and southwestward.
Figure 6. (a) Period distribution of the LSTIDs from the north in different seasons between 1 September 2009 and 21 August 2011. The solid curves are the Gaussian fitting of the histograms. (b) Horizontal velocity distribution of the LSTIDs propagating southward and southwestward.
Sensors 22 00233 g006
Figure 7. Time–month map of the occurrence rate of the LSTIDs. LT = UT + 8 h.
Figure 7. Time–month map of the occurrence rate of the LSTIDs. LT = UT + 8 h.
Sensors 22 00233 g007
Figure 8. Azimuth variation of the recorded LSTIDs between September 2009 and August 2011. The azimuths are measured clockwise from north. The red horizontal line shows the boundary between the LSTIDs propagating southwestward and southward. The two dotted blue horizontal lines show the mean azimuths of the southwestward and southward traveling LSTIDs, respectively. The vertical gray bars are used to indicate the unoperated period of the ionosonde network.
Figure 8. Azimuth variation of the recorded LSTIDs between September 2009 and August 2011. The azimuths are measured clockwise from north. The red horizontal line shows the boundary between the LSTIDs propagating southwestward and southward. The two dotted blue horizontal lines show the mean azimuths of the southwestward and southward traveling LSTIDs, respectively. The vertical gray bars are used to indicate the unoperated period of the ionosonde network.
Sensors 22 00233 g008
Table 1. The occurrence number and occurrence rate of the quiet-time LSTIDs in four seasons from 1 September 2009 to 21 August 2011.
Table 1. The occurrence number and occurrence rate of the quiet-time LSTIDs in four seasons from 1 September 2009 to 21 August 2011.
SeasonAutumnWinterSpringSummer
Number of events29110162
Occurrence rate27.0%56.7%16.3%1.0%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, R.; Chen, G.; Li, Y.; Zhang, S.; Gong, W.; He, Z.; Zhang, M. Long-Term Observation of the Quasi-3-Hour Large-Scale Traveling Ionospheric Disturbances by the Oblique-Incidence Ionosonde Network in North China. Sensors 2022, 22, 233. https://doi.org/10.3390/s22010233

AMA Style

Zhang R, Chen G, Li Y, Zhang S, Gong W, He Z, Zhang M. Long-Term Observation of the Quasi-3-Hour Large-Scale Traveling Ionospheric Disturbances by the Oblique-Incidence Ionosonde Network in North China. Sensors. 2022; 22(1):233. https://doi.org/10.3390/s22010233

Chicago/Turabian Style

Zhang, Ruijiao, Gang Chen, Yaxian Li, Shaodong Zhang, Wanlin Gong, Zhiqiu He, and Min Zhang. 2022. "Long-Term Observation of the Quasi-3-Hour Large-Scale Traveling Ionospheric Disturbances by the Oblique-Incidence Ionosonde Network in North China" Sensors 22, no. 1: 233. https://doi.org/10.3390/s22010233

APA Style

Zhang, R., Chen, G., Li, Y., Zhang, S., Gong, W., He, Z., & Zhang, M. (2022). Long-Term Observation of the Quasi-3-Hour Large-Scale Traveling Ionospheric Disturbances by the Oblique-Incidence Ionosonde Network in North China. Sensors, 22(1), 233. https://doi.org/10.3390/s22010233

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

Article Metrics

Back to TopTop