Next Article in Journal
A Rapid, Low-Cost, and High-Precision Multifrequency Electrical Impedance Tomography Data Acquisition System for Plant Phenotyping
Next Article in Special Issue
Phase Centre Corrections of GNSS Antennas and Their Consistency with ATX Catalogues
Previous Article in Journal
Empirical Models for Estimating Air Temperature Using MODIS Land Surface Temperature (and Spatiotemporal Variables) in the Hurd Peninsula of Livingston Island, Antarctica, between 2000 and 2016
Previous Article in Special Issue
Retrieval of Soil Moisture Content Based on Multisatellite Dual-Frequency Combination Multipath Errors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Comparative Study of Predominantly Daytime and Nighttime Lightning Occurrences and Their Impact on Ionospheric Disturbances

1
Department of Land Surveying and Geo-Informatics, The Hong Kong Polytechnic University, Kowloon, Hong Kong
2
Research Institute of Land and Space, The Hong Kong Polytechnic University, Kowloon, Hong Kong
3
School of Civil and Transportation Engineering, Guangdong University of Technology, Guangzhou 510006, China
4
Department of Building Environment and Energy Engineering, The Hong Kong Polytechnic University, Kowloon, Hong Kong
5
Department of Geomatic Engineering, Kwame Nkrumah University of Science and Technology, Kumasi AK-448-4746, Ghana
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(13), 3209; https://doi.org/10.3390/rs14133209
Submission received: 20 May 2022 / Revised: 28 June 2022 / Accepted: 29 June 2022 / Published: 4 July 2022
(This article belongs to the Special Issue GNSS, Space Weather and TEC Special Features)

Abstract

:
Space weather events adversely impact the operations of Global Navigation Satellite Systems (GNSS). Understanding space weather mechanisms, interactions in the atmosphere, and the extent of their impact are useful in developing prediction and mitigation models. In this study, the hourly lightning occurrence and its impact on ionospheric disturbances, quantified using the Rate of Total electron content Index (ROTI), were assessed. The linear correlation between diurnal lightning activity and ROTI in the coastal region of southern China where lightning predominates in the daytime was initially negative contrary to a positive correlation in southern Africa where lighting predominates in the evening. After appreciating and applying the physical processes of gravity waves, electromagnetic waves and the Trimpi effect arising from lightning activity, and the time delay impact they have on the ionosphere, the negative correlation was overturned to a positive one using cross-correlation. GNSS has demonstrated its capability of revealing the impact lightning has on the ionosphere at various times of the day.

Graphical Abstract

1. Introduction

The wide spatio-temporal advantages of Global Navigation Satellite Systems (GNSS) have made it an effective tool not only in positioning and navigation but also in studying and monitoring weather events. Geomagnetic and solar radiation, which are the key drivers of ionospheric plasma variations have been extensively observed using GNSS [1,2,3,4]. Tropospheric weather events such as cyclones [5,6,7,8], earthquakes [9,10,11,12,13], and/or man-made events such as rocket launches [14] and explosive bursts [15] have also been observed to cause plasma variations through the use of GNSS.
Thunderstorms/lightning, another troposphere weather event, has also garnered interest in the GNSS scientific community for adopting the advantages of GNSS to gain more insight into its activities. In recent GNSS lightning studies, Osei-Poku et al. [16] evaluated common total electron content (TEC) detrending techniques during lightning events. Rahmani et al. [17] probed the vertical coupling effect of a thunderstorm on the lower ionosphere. The ionospheric response to thunderstorms in West Africa has been reported by Ogunsua et al. [18], whereas Tang et al. [19] and Liu et al. [20] presented gravity waves resulting from thunderstorms. Blanc et al. [21] also showed gravity wave measurements in the ionosphere following thunderstorms in Europe. Lay [22] reported acoustic wave activity above thunderstorms. Amin [23] also studied the hourly lightning activity and its effect on the ionosphere using GNSS in southern Africa.
The quest to harness GNSS for more insightful studies about lightning continues. As a recommendation in their work, Amin [23] suggested the use of more long-term data to study the correlation between hourly lightning activity and irregular ionospheric events as their study had a limited dataset with few case studies. The use of long-term data would provide more consistent, reliable, and insightful information In this regard, this current study uses four-year long-term data to study the relationship between the hourly lightning activity and the rate of the TEC index (ROTI) in the low-latitude coastal sea region of southern China. Similar techniques used in Amin [23] are deployed for an effective comparison. The results are compared to Amin [23] and those in other geographical regions. The initial results in this present study showed a negative linear correlation between hourly lightning activity and ROTI as opposed to the positive linear correlation in Amin [23]. One reason could be attributed to lightning interactions with the ionosphere in the different geographical regions. In the subsequent sections, the data and methods used are described. Further investigations are made to find and discuss the physical mechanisms underlying the differences in the results. The derived conclusions are then presented.

2. Data and Methods

2.1. Lightning Data

The lightning data is from a local lightning network in the coastal region of southern China. The network consists of about 17 Vaisala LS8000 sensors that provide the geolocation and source peak currents of lightning in the area [24]. The number of times a current is recorded is positively correlated with lightning activity [25]; a day with a lightning count greater than 10,000 is deemed a “lightning day” [16]. The data span is from 2014 to 2017.

2.2. GNSS Data

The local GNSS data were obtained from the Hong Kong Satellite Reference (HK SatRef). The GNSS receivers have a sampling rate of 30 s. Information on the network is given in the works of Ji et al. [26] and Kumar et al. [27]. More information on HK SatRef is available at the website of the Hong Kong Survey Department (https://www.geodetic.gov.hk/en/rinex/downv.aspx, (accessed on 14 June 2019)).

2.3. Method

2.3.1. ROTI

ROTI defined by Pi et al. [28] as the root mean square of the TEC rate is used to characterize ionospheric irregularity. ROTI could be used as a proxy for the scintillation [29,30]. ROTI is computed from the GNSS data as follows. First, a 15° elevation cut-off angle was used to eliminate the multipath effect [6]. A geometry-free linear combination of pseudo- and carrier-phase signals was then used to compute slant TEC (STEC) at 30 s sampling intervals. STEC was converted to vertical TEC (VTEC) by applying a mapping function using Equation (1) below, where Re is the earth’s radius, θ is the elevation angle at the ionospheric pierce point (IPP) of the signal–receiver path, and hi is the ionospheric single layer, approximated at 350 km.
VTEC = 1 ( R e cos θ R e + h i ) 2 STEC
Finally, at five-minute intervals of the rate of TEC (ROT), ROTI was computed according to Equation (2), where ROT and ROTI are in TEC units (TECu: 1 TECu = 1016 e/m2) and the notation · is the averaging operation [31].
ROTI = < ROT 2 > < ROT > 2
ROTI as shown in Equation (2) is usually for a single satellite-receiver pair. ROTI values exceeding 0.2 TECu are used to indicate that ionospheric scintillation has happened [32]. However, ROTI average (ROTIavg) is the average value of ROTI over 30 min for all satellites received by a single station; following Oladipo et al. [33], a scintillation is deemed to occur for a ROTIavg value exceeding a threshold of 0.8 TECu. To effectively compare this study to that of Amin [23], ROTIavg is adopted.

2.3.2. Hourly Occurrence

For the diurnal hourly occurrence, the number of times lightning occurred and ROTIavg was greater than the threshold were accumulated within hourly intervals for all days in each year. A correlation is then looked for between the lightning and ROTIavg hourly occurrences.

2.3.3. Selection Criteria

To avoid the geomagnetic storm and solar radiation effects, only days with disturbance storm time (Dst) greater than −30 nT [34] and solar flux index (F10.7 index) less than 150 sfu (solar flux units, 1 sfu = 10−22 watt per square meter-hertz) [35] were selected. Table 1 shows the total number of days for the years 2014–2017 that had lightning activity only.

3. Results

According to Tang et al. [19], the stations of HK Satref are quite close, hence their observations are similar. Only the observations from one station (HKOH) are presented.
Figure 1 shows the ROTIavg from 2014 to 2017. In Figure 1, the highest values are mostly in the nighttime and between 0.07–0.2 TECu; this indicates that a scintillation is often a nighttime event. This is coherent with Tang et al. [31] and Ji et al. [26] who have shown that a greater percentage of ROTI lies between 0.02 and 0.05 TECu in the Hong Kong region. In this present study, the ROTIavg threshold is set at 0.075 TECu instead of 0.8 TECu, unlike Oladipo et al. [33] but the same as Nishioka et al. [36], who suggested a scintillation threshold of 0.075 TECu in the Asian region. More so, the years 2014–2017 are at the declining phase of the 24th solar cycle [37] where scintillation values are low compared to the high solar active years of 2002 and 2012 in Oladipo et al. [33] and Amin [23], respectively. Also, Jacobsen [38] and Liu et al. [39] have demonstrated that GNSS receiver types, configurations, and sampling rates influence ROTI values. These differences in ROTI values (thresholds) arising from the technique used, GNSS receiver configurations, and geographic location should be important factors to consider when developing regional and global models.
On the diurnal hourly occurrence, the number of times lightning occurred and the ROTIavg was greater than 0.075 TECu were accumulated within hourly intervals for all days in each year. Panels a to d in Figure 2 show the annual hourly diurnal ROTIavg occurrence against that of the lightning occurrence. All years showed a similar trend. The trend reveals that the lighting occurrence peaked before that of ROTIavg. This resulted in negative linear correlation values of −0.364, −0.41, −0.371, and −0.421 with significant values (p-value) of 0.05, 0.04, 0.05, and 0.04 at a confidence interval of 95% (α = 0.05) in chronological order from 2014 to 2017 as seen in panels e to h in the right column in Figure 2. The p-values indicate that despite being negative, the linear correlation is statistically significant. This implies that as lightning increases, ROTI decreases and vice versa. These are contrary to observations in Amin [23] where both lightning and ROTIavg peaked at the same time in the evening resulting in positive linear correlation values. As the objective of this study is to find the relationship between lightning and ROTI, reasons are discussed to explain the differences between this observation and that of [23].

4. Discussions

The objective of this work is to appreciate and understand the relationship between hourly lightning activity and ionospheric irregularities (ROTI) through the use of long-term data. The initial results seen in Figure 2 show a negative linear correlation contrary to the positive linear correlation reported in Amin [23]. One reason for the observed differences in linear correlation could be attributed to the interaction between local variations in mesoscale and the topographic nature of the geographic region [40,41] resulting in different diurnal trends of lightning. That is, different regions tend to respond to lightning differently. From the observed data and results in other studies, there seems to be a common trend in the diurnal pattern of lightning. Places in the northern hemisphere seem to have a common pattern. Likewise, places in the southern hemisphere. For instance, Figure 5 in Williams et al. [42] shows a similar pattern of diurnal lightning as observed in this current study, and their Figure 6 is similar to that of Figure 5.13 of Amin [23]. In Figure 5 of Williams et al. [42], Australia, which is in the southern hemisphere, shows a broad peak of lightning from noon, which lasts until about 23:00 local time (LT). This pattern is similar to that of South Africa, also in the southern hemisphere. Brazil also shows a similar trend as seen in Pinto et al. [40]. In Figure 6 of Williams et al. [42], French Guyana, which is in the northern hemisphere and in South America, exhibits a lightning peak of around 15 LT, which lasts for about 3 h (till 18 LT). This pattern is similar to that of Hong Kong. Other similar trends are observed in the Indian region [43] and the United States [41,44], all in the northern hemisphere. The broad width of the lightning peak observed in the southern hemisphere regions should cover their peak time of ROTIavg, which happens in the evening (mostly between 20 and 23 h). Statistically, lightning and ROTIavg having similar trends and with their peaks coinciding at the same time stands a high chance of giving a positive linear correlation value, whereas peaks that occur far apart at different times would give a negative linear correlation value.
Correlation may not necessarily mean causation. Possible mechanisms influencing the diurnal correlations in this work and that of Amin [23] are the lightning activity interactions with the ionosphere through lightning-related current and electrical discharges [45,46] and processes of gravity and electromagnetic wave [24,47,48,49].
A thunderstorm with lightning activity may affect the electron density in the lower ionosphere (particularly the D-layer) through two basic physical processes: the gravity wave (GW) due to the convection and thunder [47] and the electromagnetic wave (EM) due to the lightning [24]. The GW will generate a disturbance of the mesosphere (60–90 km) and the lower thermosphere (above 90 km), thus of the neutral atmosphere, which drags electrons from the ionosphere. Although the GW is a transversal wave, any changes in its amplitude generated from a lightning discharge will perform as a longitudinal wave that propagates upwards at the speed of the order of sound in air. The higher the air density, the stronger and faster the GW. Therefore, the GW due to the lightning and thunder and the thunderstorm convection disturbance would be stronger and faster as it propagates in the lower (higher air density) D-layer during the daytime than it does in the higher (lower air density) D-layer during the nighttime. The GW due to the lightning and thunder and the thunderstorm convection disturbance on ROTI lags the lightning activity due to its very low propagation speed at a high altitude. Theoretical studies have shown that the GW needs about one to several hours to propagate to higher altitudes depending on speed and period [50]. Using the wave model of Row [51] and Francis [52], which assumes that the GW propagates upwards through the atmosphere, Taylor et al. [53] found that the GW from a thunderstorm took about six hours before reaching the airglow layers (80 km above ground surface), which implies a vertical velocity of about 5 ms−1. In Figure 2, the time delay between lightning and ROTIavg is about 5–7 h, similar to that observed by [53]. At an altitude of 350 km and a time delay of 5–7 h, the vertical velocity of the GW would be approximately 10–15 ms−1. These observations show that the average vertical velocity of the GW is about 5–15 ms−1 confirming the low propagation speed of the GW. Also, the GW is predominant in the horizontal component which takes a relatively longer time to reach higher atmosphere heights [19,20,54]. This could best account for the difference in peak times of lightning and the ROTI pattern in this study, where lightning is predominantly in the afternoon. Again, Lay [22] made an observation where lightning occurring predominantly in the daytime (14–16 LT) had its related ionospheric disturbances around 00–02 LT similar to the observations in the left panels in Figure 2. Thus, lightning that occurs predominantly in the afternoon sees the ionospheric disturbance some hours later leading to an initial negative linear correlation.
Lightning-generated electromagnetic waves (EM) on the other hand propagate at the speed of light and reach the lower ionosphere almost immediately after the lightning activity has taken place. The strength of the impact of the EM on ROTI depends on the density of the electrons present. The higher the electron density, the stronger the impact. Theoretical simulations have shown that the EM could lead to a reduction in electron density at lower altitudes and an increase in electron density at higher altitudes in an ionospheric D-layer due to Joule heating effects [55,56,57], which have also been proven by lightning sferics observations [24,58,59]. Lightning sferics also produce very-low-frequency (VLF) discharges that contribute to the Trimpi effect. The Trimpi effect refers to transient perturbations caused by electron precipitation on sub ionospheric propagating waves [60]. The perturbations associated with sferics take only about 0.6 s to begin. Also, the Trimpi effect is only observed in nighttime ionospheric conditions and not daytime [61,62]. Therefore, the lightning-generated EM and the Trimpi effect would have a stronger impact on ROTI during the nighttime than during the daytime. The appearance of the impact of the EM and Trimpi effect on ROTI is almost on par with the time the lightning activity happened. Thus, lightning that occurs predominantly in the evening sees an ionospheric disturbance almost instantly leading to a positive linear correlation.
To further illustrate the time of lightning dominance and the time lags of its associated ionospheric disturbance, two diurnal zone divisions are made. That is, 07–17 and 18–06 LT to represent the daytime and evening time zones, respectively [30]. In Figure 2, 2014–2017 have similar lightning and ROTIavg trends. Some days are selected to illustrate this further as most of the individual days have similar observations. The assessment is done by finding the cross-correlation between lightning and ROTI instead of the linear correlation. Cross-correlation studies different variables to identify their similarity and draws characteristics relative to each other based on time to derive new information. Pseudorandom noise codes (PRN) from the individual satellite-receiver pairs available at the time ROTI lags the lightning are presented to provide extra evidence of ionospheric disturbances. The TEC of the PRN is detrended using the Savitzky–Golay filter of order 6 and window length of 120 min [16]. Detrended TEC (DTEC) is filtered with a bandpass of frequencies between 1 and 2.8 MHz and 4.2 and 8.2 MHz to derive ionospheric gravity (IGW) and ionospheric acoustic (IAW) waves, respectively, from lightning [17,63]. A filtered DTEC amplitude above 0.08 TECu and 0.025 TECu [63] and ROTI greater than 0.2 TECu signify the presence of IGW, IAW, and ionospheric disturbances, respectively. The selected days for the daytime and evening are presented in Figure 3 and Figure 4, respectively. The left panels of the top rows for each day in Figure 3 and Figure 4 show the ROTIavg and lightning counts. The right panels of the top rows show the cross-correlation between the ROTIavg and lightning counts from the left panel. The bottom rows show the DTEC, IGW, IAW, and ROTI of the PRN, which was available at the time ROTI lags the lightning. The magenta lines are the thresholds of IGW, IAW, and ROTI.
Figure 3 shows the days with lightning predominantly in the daytime. In panels b, h, and n in Figure 3, ROTIavg is seen to have a negative time lag of about 4 h to the lightning activity but is strongly positively correlated. The time of the disturbances as seen in the PRNs in the bottom rows confirms the delay. The negative time lags; however, positive coefficients indicate that the disturbances are associated with lightning, which can be attributed to the delayed impact of the gravity wave mechanism from the predominantly daytime lightning on the ionosphere.
Figure 4 shows the days with lightning predominantly in the evening. From panels b, h, and n in Figure 4, ROTIavg is seen to have a time lag of zero to a few minutes to the lightning activity and with strong positive correlation coefficients. The time of the disturbances as seen in the PRNs in the bottom rows is on par with the lightning activity. The time lag of zero to a few minutes and the positive coefficients indicate that the disturbances are associated with lightning, which can be explained by the almost immediate impact of the electromagnetic wave mechanism and the Trimpi effect resulting from the predominantly nighttime lightning on the ionosphere.
Following the mechanisms of the gravity wave, electromagnetic wave, and the Trimpi effect as discussed above, cross-correlation of the lightning and ROTIavg counts from the left panels in Figure 2 yielded positive correlation values of 0.89, 0.82, 0.72, and 0.61 for 2014 to 2017, respectively as shown in Figure 5. The values are similar to 0.86 reported by Amin [23]. This shows that the diurnal hourly lightning activity and associated ROTIavg have a positive correlation.
The observations in Figure 3, Figure 4 and Figure 5 confirm that ionospheric disturbances from lightning occurring in the afternoon take a relatively longer time to manifest, whereas those due to lightning occurring in the evening are almost immediate. Using cross-correlation provides a better means of accessing the relationship between lightning activity and its associated ROTI at various times of the day. The cross-correlation coefficient is statistically significant when its absolute value is greater than the expected cross-correlation coefficient (ECCF), which is computed as [64].
ECCF = 2 n | k |
where n is the total number of observations and k is the number of observations between zero and the point of the maximum cross-correlation value on the lag axis. Negative and positive k are to the left and right of zero on the lag axis, respectively.
Table 2 is a summary of the cross-correlation coefficients of days and years in Figure 3, Figure 4 and Figure 5 and their ECCF. The observations of the days in Figure 3 and Figure 4 are at 30 min intervals, whereas those of the years in Figure 5 are at hourly intervals. The cross-correlation coefficients being greater than the ECCFs show that lightning activity is significantly positively correlated with its associated ROTI regardless of the time of day the lightning was predominant.

5. Conclusions

In this work, a correlational study between hourly lightning activity and ROTI was carried out using long-term data. The data were from 2014 to 2017 with the study area being Hong Kong. The linear correlation between the hourly lightning activity and ROTI was negative compared to a similar study in southern Africa, where the correlation was positive. The initial differences in linear correlation values between this present study and that in southern Africa are attributed to the time differences at which lightning and ROTI peaked. Although lightning peaked in the daytime in Hong Kong, its associated ROTI peaked in the evening. For southern Africa, both lightning and ROTI peaked at the same time in the evening. Probing further, this current study found that the ROTI resulting from lightning predominantly in the daytime lagged due to slow propagation of gravity wave mechanisms but did not lag when lightning was predominantly in the evening due to electromagnetic waves and the Trimpi mechanisms. The temporal discrepancy could be explained by either or both mechanisms and it would require simulations to confirm their likelihood. Using datasets obtained from radio occultation missions such as COSMIC would be effective for quantifying these mechanisms and their perturbations at the various layers of the ionosphere. Unlike the linear correlation that showed that predominantly daytime lightning is negatively correlated with its ROTI, cross-correlation offered a better means to access the lightning activity–ROTI relationship. Cross-correlation revealed the time delays and showed that the lightning activity and its associated ROTI are positively correlated, which can be explained by the physical mechanisms. The results and observations have shown that GNSS can reveal the impact lightning activity has on the ionosphere at various times of the day. The models for ionospheric scintillation simulations, prediction, and forecasting purposes based on lightning activities should consider and incorporate this observation when being developed.

Author Contributions

Conceptualization and research design were done by L.O.-P. and W.C. L.O.-P., L.T., M.C. and A.A.A. contributed to data processing and analysis. W.C. and M.C. supervised the work. Writing and editing of the manuscript were done by all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the University Grants Committee of Hong Kong under the scheme Research Impact Fund on the project R5009-21, and the Research Institute of Land and Space, Hong Kong Polytechnic University, and the National Natural Science Foundation of China (Grant number 41804021).

Data Availability Statement

The Dst-index data (IAGA 2002-like format) can be obtained from the Data Analysis Center for Geomagnetism and Space Magnetism, Kyoto University, operating WDC for Geomagnetism, Kyoto (http://wdc.kugi.kyoto-u.ac.jp/dstae/index.html). The Bz component of the interplanetary magnetic field was obtained from the GSFC/SPDF OMNIWeb interface at http://omniweb.gsfc.nasa.gov. The GNSS data can be obtained from Hong Kong SatRef of the Lands Department of the Hong Kong Government (https://www.geodetic.gov.hk/en/rinex/downv.aspx). All data center websites were accessed on 14 June 2019.

Acknowledgments

The research institutes and government departments where data was obtained for this study are duly acknowledged for their continuous support of scientific studies. The authors also thank the anonymous reviewers for their useful and insightful comments for improving the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cherniak, I.; Zakharenkova, I. Large-Scale Traveling Ionospheric Disturbances Origin and Propagation: Case Study of the December 2015 Geomagnetic Storm. Space Weather 2018, 16, 1377–1395. [Google Scholar] [CrossRef]
  2. Lu, Y.; Wang, Z.; Ji, S.; Chen, W. Assessing the Positioning Performance under the Effects of Strong Ionospheric Anomalies with Multi-GNSS in Hong Kong. Radio Sci. 2020, 55, 1–18. [Google Scholar] [CrossRef]
  3. Afraimovich, E.L.; Astafyeva, E.I.; Demyanov, V.V.; Edemskiy, I.K.; Gavrilyuk, N.S.; Ishin, A.B.; Kosogorov, E.A.; Leonovich, L.A.; Lesyuta, O.S.; Palamartchouk, K.S.; et al. A review of GPS/GLONASS studies of the ionospheric response to natural and anthropogenic processes and phenomena. J. Space Weather Space Clim. 2013, 3, A27. [Google Scholar] [CrossRef] [Green Version]
  4. Lay, E.H.; Shao, X.M.; Carrano, C.S. Variation in total electron content above large thunderstorms. Geophys. Res. Lett. 2013, 40, 1945–1949. [Google Scholar] [CrossRef]
  5. Chen, G.; Zhou, C.; Liu, Y.; Zhao, J.; Tang, Q.; Wang, X.; Zhao, Z. A statistical analysis of medium-scale traveling ionospheric disturbances during 2014–2017 using the Hong Kong CORS network. Earth Planets Space 2019, 71, 52. [Google Scholar] [CrossRef]
  6. Yang, Z.; Liu, Z. Observational study of ionospheric irregularities and GPS scintillations associated with the 2012 tropical cyclone Tembin passing Hong Kong. J. Geophys. Res. Space Phys. 2016, 121, 4705–4717. [Google Scholar] [CrossRef]
  7. Polyakova, A.S.; Perevalova, N.P. Investigation into impact of tropical cyclones on the ionosphere using GPS sounding and NCEP/NCAR Reanalysis data. Adv. Space Res. 2011, 48, 1196–1210. [Google Scholar] [CrossRef]
  8. Mao, T.; Wang, J.; Yang, G.; Yu, T.; Ping, J.; Suo, Y. Effects of typhoon Matsa on ionospheric TEC. Chin. Sci. Bull. 2010, 55, 712–717. [Google Scholar] [CrossRef]
  9. Kong, J.; Yao, Y.; Zhou, C.; Liu, Y.; Zhai, C.; Wang, Z.; Liu, L. Tridimensional reconstruction of the Co-Seismic Ionospheric Disturbance around the time of 2015 Nepal earthquake. J. Geod. 2018, 92, 1255–1266. [Google Scholar] [CrossRef]
  10. Miyazaki, S.I.; Larson, K.M. Coseismic and early postseismic slip for the 2003 Tokachi-oki earthquake sequence inferred from GPS data. Geophys. Res. Lett. 2008, 35, L04302. [Google Scholar] [CrossRef] [Green Version]
  11. Craig, T.J.; Calais, E. Strain accumulation in the New Madrid and Wabash Valley seismic zones from 14 years of continuous GPS observation. J. Geophys. Res. Solid Earth 2014, 119, 9110–9129. [Google Scholar] [CrossRef]
  12. Calais, E.; Minster, J.B. GPS, earthquakes, the ionosphere, and the Space Shuttle. Phys. Earth Planet. Inter. 1998, 105, 167–181. [Google Scholar] [CrossRef]
  13. Afraimovich, E.L.; Ding, F.; Kiryushkin, V.V.; Astafyeva, E.I.; Jin, S.; Sankov, V.A. TEC response to the 2008 Wenchuan Earthquake in comparison with other strong earthquakes. Int. J. Remote Sens. 2010, 31, 3601–3613. [Google Scholar] [CrossRef] [Green Version]
  14. Mendillo, M. The effect of rocket launches on the ionosphere. Adv. Space Res. 1981, 1, 275–290. [Google Scholar] [CrossRef]
  15. Kundu, B.; Senapati, B.; Matsushita, A.; Heki, K. Atmospheric wave energy of the 2020 August 4 explosion in Beirut, Lebanon, from ionospheric disturbances. Sci. Rep. 2021, 11, 2793. [Google Scholar] [CrossRef]
  16. Osei-Poku, L.; Tang, L.; Chen, W.; Mingli, C. Evaluating Total Electron Content (TEC) Detrending Techniques in Determining Ionospheric Disturbances during Lightning Events in A Low Latitude Region. Remote Sens. 2021, 13, 4753. [Google Scholar] [CrossRef]
  17. Rahmani, Y.; Alizadeh, M.M.; Schuh, H.; Wickert, J.; Tsai, L.-C. Probing vertical coupling effects of thunderstorms on lower ionosphere using GNSS data. Adv. Space Res. 2020, 66, 1967–1976. [Google Scholar] [CrossRef]
  18. Ogunsua, B.O.; Srivastava, A.; Bian, J.; Qie, X.; Wang, D.; Jiang, R.; Yang, J. Significant Day-time Ionospheric Perturbation by Thunderstorms along the West African and Congo Sector of Equatorial Region. Sci. Rep. 2020, 10, 8466. [Google Scholar] [CrossRef]
  19. Tang, L.; Chen, W.; Chen, M.; Osei-Poku, L. Statistical observation of thunderstorm-induced ionospheric Gravity waves above low-latitude areas in the northern hemisphere. Remote Sens. 2019, 11, 2732. [Google Scholar] [CrossRef] [Green Version]
  20. Liu, T.; Yu, Z.; Ding, Z.; Nie, W.; Xu, G. Observation of Ionospheric Gravity Waves Introduced by Thunderstorms in Low Latitudes China by GNSS. Remote Sens. 2021, 13, 4131. [Google Scholar] [CrossRef]
  21. Blanc, E.; Ceranna, L.; Hauchecorne, A.; Charlton-Perez, A.; Marchetti, E.; Evers, L.G.; Kvaerna, T.; Lastovicka, J.; Eliasson, L.; Crosby, N.B.; et al. Toward an Improved Representation of Middle Atmospheric Dynamics Thanks to the ARISE Project. Surv. Geophys. 2018, 39, 171–225. [Google Scholar] [CrossRef] [Green Version]
  22. Lay, E.H. Ionospheric Irregularities and Acoustic/Gravity Wave Activity Above Low-Latitude Thunderstorms. Geophys. Res. Lett. 2018, 45, 90–97. [Google Scholar] [CrossRef] [Green Version]
  23. Amin, M.M. Influence of Lightning on Electron Density Variation in the Ionosphere Using WWLLN Lightning Data and GPS Data. Master’s Thesis, University of Cape Town, Cape Town, South Africa, 2015. [Google Scholar]
  24. Qin, Z.; Chen, M.; Lyu, F.; Cummer, S.A.; Gao, Y.; Liu, F.; Zhu, B.; Du, Y.P.; Usmani, A.; Qiu, Z. Prima Facie Evidence of the Fast Impact of a Lightning Stroke on the Lower Ionosphere. Geophys. Res. Lett. 2020, 47, e2020GL090274. [Google Scholar] [CrossRef]
  25. Salem, M.A.; Liu, N.; Rassoul, H.K. Modification of the lower ionospheric conductivity by thunderstorm electrostatic fields. Geophys. Res. Lett. 2016, 43, 5–12. [Google Scholar] [CrossRef]
  26. Ji, S.; Chen, W.; Wang, Z.; Xu, Y.; Weng, D.; Wan, J.; Fan, Y.; Huang, B.; Fan, S.; Sun, G.; et al. A study of occurrence characteristics of plasma bubbles over Hong Kong area. Adv. Space Res. 2013, 52, 1949–1958. [Google Scholar] [CrossRef]
  27. Kumar, S.; Chen, W.; Chen, M.; Liu, Z.; Singh, R.P. Thunderstorm-/lightning-induced ionospheric perturbation: An observation from equatorial and low-latitude stations around Hong Kong. J. Geophys. Res. Space Phys. 2017, 122, 9032–9044. [Google Scholar] [CrossRef]
  28. Pi, X.; Mannucci, A.J.; Lindqwister, U.J.; Ho, C.M. Monitoring of global ionospheric irregularities using the worldwide GPS network. Geophys. Res. Lett. 1997, 24, 2283–2286. [Google Scholar] [CrossRef]
  29. Yang, Z.; Liu, Z. Correlation between ROTI and Ionospheric Scintillation Indices using Hong Kong low-latitude GPS data. GPS Solut. 2016, 20, 815–824. [Google Scholar] [CrossRef]
  30. Yang, Z.; Liu, Z. Low-Latitude Ionospheric Density Irregularities and Associated Scintillations Investigated by Combining COSMIC RO and Ground-Based Global Positioning System Observations Over a Solar Active Period. J. Geophys. Res. Space Phys. 2018, 123, 3998–4014. [Google Scholar] [CrossRef]
  31. Tang, L.; Louis, O.-P.; Chen, W.; Chen, M. A ROTI-Aided Equatorial Plasma Bubbles Detection Method. Remote Sens. 2021, 13, 4356. [Google Scholar] [CrossRef]
  32. Wu, Q.; Chou, M.-Y.; Schreiner, W.; Braun, J.; Pedatella, N.; Cherniak, I. COSMIC observation of stratospheric gravity wave and ionospheric scintillation correlation. J. Atmos. Sol. Terr. Phys. 2021, 217, 105598. [Google Scholar] [CrossRef]
  33. Oladipo, O.A.; Schüler, T. Equatorial ionospheric irregularities using GPS TEC derived index. J. Atmos. Sol. Terr. Phys. 2013, 92, 78–82. [Google Scholar] [CrossRef]
  34. Loewe, C.A. Classification and mean behavior of magnetic storms. J. Geophys. Res. A Space Phys. 1997, 102, 14209–14213. [Google Scholar] [CrossRef]
  35. Yu, S.; Liu, Z. The ionospheric condition and GPS positioning performance during the 2013 tropical cyclone Usagi event in the Hong Kong region. Earth Planets Space 2021, 73, 66. [Google Scholar] [CrossRef]
  36. Nishioka, M.; Saito, A.; Tsugawa, T. Occurrence characteristics of plasma bubble derived from global ground-based GPS receiver networks. J. Geophys. Res. Space Phys. 2008, 113. [Google Scholar] [CrossRef]
  37. Patel, K.; Singh, A.; Singh, S.B.; Singh, A.K. Causes responsible for intense and severe storms during the declining phase of Solar Cycle 24. J. Astrophys. Astron. 2019, 40, 4. [Google Scholar] [CrossRef]
  38. Jacobsen, K.S. The impact of different sampling rates and calculation time intervals on ROTI values. J. Space Weather Space Clim. 2014, 4, A33. [Google Scholar] [CrossRef] [Green Version]
  39. Liu, Z.; Yang, Z.; Xu, D.; Morton, Y.J. On Inconsistent ROTI Derived From Multiconstellation GNSS Measurements of Globally Distributed GNSS Receivers for Ionospheric Irregularities Characterization. Radio Sci. 2019, 54, 215–232. [Google Scholar] [CrossRef]
  40. Pinto, I.R.C.A.; Pinto, O.; Rocha, R.M.L.; Diniz, J.H.; Carvalho, A.M.; Filho, A.C. Cloud-to-ground lightning in southeastern Brazil in 1993: 2. Time variations and flash characteristics. J. Geophys. Res. Atmos. 1999, 104, 31381–31387. [Google Scholar] [CrossRef] [Green Version]
  41. Lopez, R.E.; Holle, R.L. Diurnal and Spatial Variability of Lightning Activity in Northeastern Colorado and Central Florida during the Summer. Mon. Weather Rev. 1986, 114, 1288–1312. [Google Scholar] [CrossRef] [Green Version]
  42. Williams, E.R.; Heckman, S.J. The local diurnal variation of cloud electrification and the global diurnal variation of negative charge on the Earth. J. Geophys. Res. Atmos. 1993, 98, 5221–5234. [Google Scholar] [CrossRef]
  43. Kandalgaonkar, S.S.; Tinmaker, M.I.R.; Kulkarni, J.R.; Nath, A. Diurnal variation of lightning activity over the Indian region. Geophys. Res. Lett. 2003, 30. [Google Scholar] [CrossRef]
  44. Easterling, D.R.; Robinson, P.J. The Diurnal Variation of Thunderstorm Activity in the United States. J. Appl. Meteorol. Climatol. 1985, 24, 1048–1058. [Google Scholar] [CrossRef]
  45. Kuo, C.L.; Lee, L.C. Ionospheric plasma dynamics and instability caused by upward currents above thunderstorms. J. Geophys. Res. Space Phys. 2015, 120, 3240–3253. [Google Scholar] [CrossRef] [Green Version]
  46. Collier, A.B.; Hughes, A.R.W.; Lichtenberger, J.; Steinbach, P. Seasonal and diurnal variation of lightning activity over southern Africa and correlation with European whistler observations. Ann. Geophys. 2006, 24, 529–542. [Google Scholar] [CrossRef] [Green Version]
  47. Alexander, M.J.; Holton, J.R.; Durran, D.R. The Gravity Wave Response above Deep Convection in a Squall Line Simulation. J. Atmos. Sci. 1995, 52, 2212–2226. [Google Scholar] [CrossRef] [Green Version]
  48. Vadas, S.L.; Liu, H.-I. Generation of large-scale gravity waves and neutral winds in the thermosphere from the dissipation of convectively generated gravity waves. J. Geophys. Res. Space Phys. 2009, 114, A10310. [Google Scholar] [CrossRef]
  49. Pasko, V.P.; Inan, U.S.; Bell, T.F.; Reising, S.C. Mechanism of ELF radiation from sprites. Geophys. Res. Lett. 1998, 25, 3493–3496. [Google Scholar] [CrossRef] [Green Version]
  50. Occhipinti, G.; Rolland, L.; Lognonné, P.; Watada, S. From Sumatra 2004 to Tohoku-Oki 2011: The systematic GPS detection of the ionospheric signature induced by tsunamigenic earthquakes. J. Geophys. Res. Space Phys. 2013, 118, 3626–3636. [Google Scholar] [CrossRef]
  51. Row, R.V. Acoustic-gravity waves in the upper atmosphere due to a nuclear detonation and an earthquake. J. Geophys. Res. 1967, 72, 1599–1610. [Google Scholar] [CrossRef]
  52. Francis, S.H. Global propagation of atmospheric gravity waves: A review. J. Atmos. Terr. Phys. 1975, 37, 1011–1054. [Google Scholar] [CrossRef]
  53. Taylor, M.J.; Hapgood, M.A. Identification of a thunderstorm as a source of short period gravity waves in the upper atmospheric nightglow emissions. Planet. Space Sci. 1988, 36, 975–985. [Google Scholar] [CrossRef]
  54. Chimonas, G. Enhancement of sporadicEby horizontal transport within the layer. J. Geophys. Res. 1971, 76, 4578–4586. [Google Scholar] [CrossRef]
  55. Taranenko, Y.N.; Inan, U.S.; Bell, T.F. Interaction with the lower ionosphere of electromagnetic pulses from lightning: Heating, attachment, and ionization. Geophys. Res. Lett. 1993, 20, 1539–1542. [Google Scholar] [CrossRef]
  56. Marshall, R.A.; Inan, U.S.; Chevalier, T.W. Early VLF perturbations caused by lightning EMP-driven dissociative attachment. Geophys. Res. Lett. 2008, 35. [Google Scholar] [CrossRef]
  57. Marshall, R.A. Effect of self-absorption on attenuation of lightning and transmitter signals in the lower ionosphere. J. Geophys. Res. Space Phys. 2014, 119, 4062–4076. [Google Scholar] [CrossRef]
  58. Cheng, Z.; Cummer, S.A.; Baker, D.N.; Kanekal, S.G. Nighttime D region electron density profiles and variabilities inferred from broadband measurements using VLF radio emissions from lightning. J. Geophys. Res. Space Phys. 2006, 111. [Google Scholar] [CrossRef] [Green Version]
  59. Shao, X.M.; Lay, E.H.; Jacobson, A.R. Reduction of electron density in the night-time lower ionosphere in response to a thunderstorm. Nat. Geosci. 2012, 6, 29–33. [Google Scholar] [CrossRef]
  60. Armstrong, W.C. Recent advance from studies of the Trimpi effect. Antarct. Jour. 1983, 18, 281–290. [Google Scholar]
  61. Leyser, T.B.; Inan, U.S.; Carpenter, D.L.; Trimpi, M.L. Diurnal variation of burst precipitation effects on subionospheric VLF/LF signal propagation near L = 2. J. Geophys. Res. Space Phys. 1984, 89, 9139–9143. [Google Scholar] [CrossRef]
  62. Rodger, C.J. Subionospheric VLF perturbations associated with lightning discharges. J. Atmos. Sol. Terr. Phys. 2003, 65, 591–606. [Google Scholar] [CrossRef]
  63. Lay, E.H.; Shao, X.M.; Kendrick, A.K.; Carrano, C.S. Ionospheric acoustic and gravity waves associated with midlatitude thunderstorms. J. Geophys. Res. A Space Phys. 2015, 120, 6010–6020. [Google Scholar] [CrossRef]
  64. Interpret all Statistics and Graphs for Cross Correlation. Available online: https://support.minitab.com/en-us/minitab/18/help-and-how-to/modeling-statistics/time-series/how-to/cross-correlation/interpret-the-results/all-statistics-and-graphs/ (accessed on 10 June 2022).
Figure 1. ROTIavg for the years 2014–2017. The abscissa axis shows the local time (LT: UT + 8) in hours. The ordinate axis shows the months of the year. The color bar shows the ROTIavg values in TECu.
Figure 1. ROTIavg for the years 2014–2017. The abscissa axis shows the local time (LT: UT + 8) in hours. The ordinate axis shows the months of the year. The color bar shows the ROTIavg values in TECu.
Remotesensing 14 03209 g001
Figure 2. Panels in the left column (ad) show the diurnal hourly occurrence of ROTIavg and lightning. Panels in the right column (eh) are scatter diagrams showing the linear correlation between the diurnal hourly occurrence of lightning and ROTIavg for the respective years.
Figure 2. Panels in the left column (ad) show the diurnal hourly occurrence of ROTIavg and lightning. Panels in the right column (eh) are scatter diagrams showing the linear correlation between the diurnal hourly occurrence of lightning and ROTIavg for the respective years.
Remotesensing 14 03209 g002
Figure 3. Some selected days (21 November 2015, 27 July 2015, 22 July 2015) where lightning was predominantly in the daytime (07–17 local time) and its cross-correlation with ROTIavg (top rows of each day). The associated ionospheric disturbances are also indicated by DTEC, IGW, IAW, and ROTI (bottom rows of each day) by satellites that passed some hours after the lightning activity. The magenta lines show the threshold of ±0.08 TECu, ±0.025 TECu and 0.2 TECu for IGW, IAW and ROTI respectively.
Figure 3. Some selected days (21 November 2015, 27 July 2015, 22 July 2015) where lightning was predominantly in the daytime (07–17 local time) and its cross-correlation with ROTIavg (top rows of each day). The associated ionospheric disturbances are also indicated by DTEC, IGW, IAW, and ROTI (bottom rows of each day) by satellites that passed some hours after the lightning activity. The magenta lines show the threshold of ±0.08 TECu, ±0.025 TECu and 0.2 TECu for IGW, IAW and ROTI respectively.
Remotesensing 14 03209 g003
Figure 4. Some selected days (9 August 2015, 11 July 2016, 2 July 2015) where lightning was predominantly in the evening (18–06 local time) and its cross-correlation with ROTIavg (top rows of each day). The associated ionospheric disturbances indicated by DTEC, IGW, IAW, and ROTI (bottom rows of each day) by satellites that passed at the time of lightning activity. The magenta lines show the threshold of ±0.08 TECu, ±0.025 TECu and 0.2 TECu for IGW, IAW and ROTI respectively.
Figure 4. Some selected days (9 August 2015, 11 July 2016, 2 July 2015) where lightning was predominantly in the evening (18–06 local time) and its cross-correlation with ROTIavg (top rows of each day). The associated ionospheric disturbances indicated by DTEC, IGW, IAW, and ROTI (bottom rows of each day) by satellites that passed at the time of lightning activity. The magenta lines show the threshold of ±0.08 TECu, ±0.025 TECu and 0.2 TECu for IGW, IAW and ROTI respectively.
Remotesensing 14 03209 g004
Figure 5. Cross-correlation of lightning and ROTIavg counts from left panels in Figure 2.
Figure 5. Cross-correlation of lightning and ROTIavg counts from left panels in Figure 2.
Remotesensing 14 03209 g005
Table 1. Number of lightning days in each year void of geomagnetic storm and solar radiation effects.
Table 1. Number of lightning days in each year void of geomagnetic storm and solar radiation effects.
YearNumber of Lightning Days
201485
201581
2016108
2017102
Table 2. Cross-correlation coefficient of days and years in Figure 3, Figure 4 and Figure 5 and their ECCF.
Table 2. Cross-correlation coefficient of days and years in Figure 3, Figure 4 and Figure 5 and their ECCF.
Cross-Correlation CoefficientnkECCF
Days in Figure 321 November 20150.7648−90.32
27 July 20150.7648−60.31
22 July 20150.848−110.33
Days in Figure 49 August 20150.764800.3
2 July 20150.8948−10.29
11 July 20160.9248−10.29
Years in Figure 520140.8924−80.5
20150.8224−80.5
20160.7124−80.5
20170.624−70.4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Osei-Poku, L.; Tang, L.; Chen, W.; Chen, M.; Acheampong, A.A. Comparative Study of Predominantly Daytime and Nighttime Lightning Occurrences and Their Impact on Ionospheric Disturbances. Remote Sens. 2022, 14, 3209. https://doi.org/10.3390/rs14133209

AMA Style

Osei-Poku L, Tang L, Chen W, Chen M, Acheampong AA. Comparative Study of Predominantly Daytime and Nighttime Lightning Occurrences and Their Impact on Ionospheric Disturbances. Remote Sensing. 2022; 14(13):3209. https://doi.org/10.3390/rs14133209

Chicago/Turabian Style

Osei-Poku, Louis, Long Tang, Wu Chen, Mingli Chen, and Akwasi Afrifa Acheampong. 2022. "Comparative Study of Predominantly Daytime and Nighttime Lightning Occurrences and Their Impact on Ionospheric Disturbances" Remote Sensing 14, no. 13: 3209. https://doi.org/10.3390/rs14133209

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