remote sensing
Letter
Statistical Observation of Thunderstorm-Induced Ionospheric Gravity Waves above Low-Latitude Areas in the Northern Hemisphere
Long Tang 1,2, Wu Chen 2,*, Mingli Chen 3 and Osei-Poku Louis 2 1 School of Civil and Transportation Engineering, Guangdong University of Technology, Guangzhou 510006, China 2 Department of Land Surveying and Geo-Informatics, The Hong Kong Polytechnic University, Hong Kong, China 3 Department of Building Services Engineering, The Hong Kong Polytechnic University, Hong Kong, China * Correspondence: wu.chen@polyu.edu.hk
Received: 30 September 2019; Accepted: 19 November 2019; Published: 21 November 2019
Abstract: Gravity waves (GWs) generated in the lower atmosphere can propagate upwards to ionospheric height. In this study, we investigated the correlation between ionospheric GWs detected by Global Navigation Satellite System (GNSS)-derived total electron content data and thunderstorm events recorded by a local lightning-detection network in the low-latitude region of Southern China during a four-year period, from 2014 to 2017. Ionospheric GWs were detected on both thunderstorm and non-thunderstorm days. Daytime ionospheric GW activity on high-thunderstorm days showed a similar convex-function-like diurnal variation to thunderstorm activity, which is diﬀerent to the concave-function-like pattern on non-thunderstorm days. Daytime ionospheric GW activity on low-thunderstorm days showed an approximately linear rising trend and was of a larger magnitude than that of high-thunderstorm days, suggesting it may be mixed by non-thunderstorm origins. Night-time enhancement of ionospheric GW activity was observed on thunderstorm days but not on non-thunderstorm days. Furthermore, ionospheric GW activity on thunderstorm days showed a positive correlation to solar activity. These ﬁndings can eﬀectively distinguish thunderstorm-related ionospheric GWs from those of non-thunderstorm origins and provide more comprehensive knowledge of thunderstorm–ionosphere coupling in low-latitude areas.
Keywords: GNSS; thunderstorm; ionosphere; gravity waves; low-latitude area

1. Introduction
Gravity waves (GWs) are a crucial part of the dynamics process of the atmosphere [1]. GWs inﬂuence the distribution of atmospheric winds and temperature by transferring energy and momentum between diﬀerent atmospheric layers [2,3]. When they propagate upwards to ionospheric height and interact with ionospheric plasma, the plasma density is perturbed (electron density is also changed on the basis of the charge conservation law) and generates so-called ionospheric disturbances [4,5]. Here, we call the GW signals of ionospheric disturbances ionospheric GWs. Studies have shown that various origins, including solid earth events, such as earthquakes and tsunamis [6–9], and tropospheric weather events, such as typhoons, penetrative convections, and jet streams [10–12], can trigger GWs and propagate upwards to ionospheric height.
Thunderstorms, as a kind of commonly occurring tropospheric weather event, have attracted the focus of scientiﬁc communities associated with the generation and distribution of ionospheric GWs. Pierce and Coroniti [13] ﬁrst associated GWs with thunderstorm–ionosphere coupling in the 1960s. In

Remote Sens. 2019, 11, 2732; doi:10.3390/rs11232732

www.mdpi.com/journal/remotesensing

Remote Sens. 2019, 11, 2732

2 of 11

recent decades, with the development of ionospheric detection technologies, especially the high spatial and temporal resolution of the Global Navigation Satellite System’s (GNSS) total electron content (TEC) data, many thunderstorm-related ionospheric GWs have been reported. Thunderstorm-related ionospheric GWs were reported at diﬀerent altitudes in the D [14–16] and F regions [17–19], and at diﬀerent latitudes, including both middle-latitude [20,21] and low-latitude [22] areas.
Most reported studies on thunderstorm-related ionospheric GWs are based on single or several cases. In order to better understand the correlation between ionospheric GWs and thunderstorm events, a statistical study using large amounts of data is necessary. Lay et al. [23] performed three months of statistical research comparing ionospheric GWs observed by GNSS TEC with Next Generation Weather Radar thunderstorm measurements in the middle-latitude Great Plains of the United States. They indicated that the disturbed coverage and amplitude of ionospheric GWs scaled with thunderstorm activity, both high and low activities. In addition, Lay [24] employed a similar method to a low-latitude area in South America and concluded that thunderstorms can exert a noticeable eﬀect on the coverage and magnitude of ionospheric GWs. Statistical studies with several months of data by Lay et al. [23] and Lay [24] showed a correlation between ionospheric GWs and thunderstorm events. However, more long-term data, not only on thunderstorm days (high thunderstorms and low thunderstorms) but also on non-thunderstorm days, should be analyzed to strengthen this correlation, especially in low-latitude areas, due to complex plasma dynamics. In addition, the study area in the low latitudes of the northern hemisphere was not involved.
In this study, we investigated the correlation between ionospheric GWs and thunderstorm events in the low latitudes of the northern hemisphere using four-year GNSS TEC data and thunderstorm data from Southern China. According to a local lightning-detection network, this area has tremendous thunderstorm events in spring and summer.
2. Data and Methods
2.1. Data
Thunderstorm and GNSS data from 2014 to 2017 were employed in the study. According to the solar-ﬂux index, 2014 was a high activity year and 2017 was a low activity year (https: //www.swpc.noaa.gov/). Thunderstorm data were detected by a local lightning-detection network that we developed and established in Southern China (magenta points in Figure 1a). The lightning-detection network can provide the location and magnitude of thunderstorm activity. The GNSS observation data with a 30 s sampling were collected from the public website of the Hong Kong Geodetic Survey Services (https://www.geodetic.gov.hk/; red points in Figure 1b). Here, multi-system GNSS data, including Global Positioning System (GPS) and GLObal NAvigation Satellite System (GLONASS), were employed to increase data coverage. BeiDou System (BDS) data have been available since 2016 and were not used in this study to maintain data-coverage consistency across all years [25,26]. Considering that the distance between adjacent stations was only about 8–15 km, only three stations, HKKT, HKNP, and HKOH, located in the north, west, and east of the network, were selected. The study area is located in the low latitudes of the northern hemisphere, from 106◦ to 122◦ East, and 14◦ to 30◦ North (red rectangle in Figure 1a), which was selected according to data coverage.

Remote Sens. 2019, 11, 2732 Remote Sens. 2019, 11, x FOR PEER REVIEW
(a) 40°N

(b) 36'

3 of 11 3 of 11

30°N

30' HKKT
24'

20°N 10°N

22°N 18.00'
12'

HKNP

HKOH

0°90°E

100°E

110°E

0 120°E

km

500 1000 1500

130°E

140°E

6'

54'

114°E

6'

km

0 4 8 12

12'

18'

24'

Figure 1. (a) Locations of the lightning network (magenta points) and Global Navigation Satellite
SFyisgtuemre(1G.N(aS)S)Lnoectawtioornks(bolfatchkerelcigtahntngilne)g. Tnheetwwohrikte(lminaegsehnotwa sptohientms)agannedtiGc eloqbuaaltoNr.aTvhigearteidonreScatatnelgliltee inSydsicteamtes(tGhNe SstSu)dnyeatwreoar.k(b()bLlaacrkger-escctaalnegmlea).pTfohre lowchaittieonlsinoef sthheowGsNtShSenmetawgonrekti(creedqupaotionrt.s)T. hTehrreeed srteactitoannsg—le iHndKiKcaTt,eHs tKhNe sPt,uadnydaHreKaO. (Hb), Lwaerrgeee‐smcaplleomyeadpifnotrhlioscsattuiodnys. of the GNSS network (red points).
2.2. MTehtrheoedsstations— HKKT, HKNP, and HKOH, were employed in this study.

2.2. MTheethioodnsospheric TEC series was calculated by the geometry-free combination of dual-frequency

GNSS carrier phases for each satellite–receiver pair [27]. The cut-oﬀ elevation angle was set as 20◦

for eaTchhesaiotenlolistep–hreercieciTveErCpsaeirr.ieAs wsiansglcea-llcauyleartemdobdyelthine gScehoameert[r2y8‐]frweeitchoamhbeiingahtitoonf o3f50dukaml‐fwreaqsuuesnecdy

tGo NcoSmS pcaurtreietrhpe hioanseosspfohreeraicchpiseartceellpitoe–inretcpeoivsietriopna.irW[2e7e].mTphleocyuetd‐otfhfeelseevcaotniodn-oarndgelrecwenatsesredt iaﬀse2r0e°nfcoer

meaecthhosdatteolleitliem–riencaeteivtehrepbaaicrk. Agrosuinngdlet‐rleanyderinmtohdeeTlEinC Sscehriaeesra[n2d8]ewxtirtahctaedhethigehitoonfo3sp50hekrmic dwisatsuurbsaendcteo

sceormiesp[u2t9e].the ionospheric pierce point position. We employed the second‐order center difference

wSwdmehithestettieruhnTreTrogehbτhdτeτaernrae=taemcpompe5rprepsemeslelliseiitriemtnuniunetdisdtisnssee[ta2tshothou9eeffe]ia.tdadtaiwhbwiﬀfelafeeaevrvbrefeeoeannwcrcwkcetiehgitsthreshtotepdepupepeenr[trd[3ieio30co0dt]dti].ro.TeTTTnnhihisdoseemfmaaiimounmunllptptotilihplsipitelptluiuieheddTdedeeErbibrCrcyayatGt4i4isosoWesiirfnfnaisa4e4c,c((stbπtπooeττrarc//ncTacTadua))nnsaaeefsfstxteeteethrrrrrvveatteehchataeameessdddpaaiilﬀfibtbftheeaaurernnededdneni-co‐coppenefaaoppsassrssrpowoﬁfhcicalleeettvseresseirsrc,..,

wTidwSiioosEaienttinCTthotatohEsiuaanspCaptngphthpiueehτteerrne(riir=TicoiibtcodE5Ge(dGTCgmsWWEicsUniocCsnn)osp;aUipeai1ncse)cc6g;Tcso.61o6uEar.–r6dniTCdt2–idaEUn0i2bnCg0emlg=nUemtidotfn1o=io0noninr11sof6io0tiesahsi1eams6eaelcee/mlphmdlel/elpvemai2vﬁtelr.eie2clefH.clid(ateH(giergd(oeererennear(nm,reetoeaira,oirfatnsaialsfionlloaaiylncfynamalgtmoeolcpessitrpnslpsoiisltrhstituhttheilhudsaareadinrelcdngaet0eriGh0ﬀt.gr1h.rWe1eterTrhreTsesEatshnE,hnhCobcCaoleU1endUl)cdp).a1.)orT.uo)NfoN.hfsc0oTee0eo.th1s.ettu15seeh5nt.etThuihTtaEnaaEotmiCttfCttUhpoThUeflEeiwT2twC2uτaEτadasCtstneiicmmcdahohneofedo-ds‐apTspeeedEnwenrTCrifaioEfooviodCsrder

wdaastagAAartimdmthdeeetethhdbooedodgniddnteoenssiccanrrgii0bb.ae5end◦ddb×beyyn0LdL.5aa◦oyyfleeoettanacaglh.li.[ta2u[23rdc3]e]wa–wrlaeasatmsietmueismdspielpnolggyoreiyindded,tihnaientnhddteihfstfehteasertteiatsnittmciiscetaeiplciraanonltcaaeelrnsyvsasa.liysl.swTisha. seTsshteuetdstytouad3r0yeaamwreinaas.

WTcGrgTlohieWreghuucihudnuonstndtrs,dnsdedet,eidhednertdehrsgtsnt.heootnsenoWrtthtmrrteomheoaﬀatkceuatoecccsvcetoe0ieclevvdl.olr5leewau°wrtggnhaaer×atgesissade,e0eneﬀcf.tﬀda5fhoeene°ecuiccdotnnttliniiovtvvoitsneoehes.gtnpao[giFo2thtrruu4seiecp]ddrrep;tihelcihrlcfe–eeGorlaswraiueWcmgtannirtstotGiucsdrtoedWhetfvce,eofeeetclrhccgrloaoteersgivpviuvmedeererr,erdpaaas[uaxage2gsirnen4emsi.d]nte;Iudgdtfihtmufeathehraeagiaecncrhmoatgigidmv3rmpi0eecldperiaemtlalcuiclgtihiedunnselutde.l3etreiI0orspmfvutfahamaeresrlpaspeigeacnswerdhshrisdaiot3eeoislmdﬀddclieesgewotcelhhlfttepiteisvrnoetueeoiarnnrrimogpoeg3srcadp0pisosdlthsrirwdtmeeocudereekididlrnecl3ee..

fateomharrcpieholsinhteuofofdsledpecshtoieovfrfeiiotcghnGreoiWdsepﬀcsheecdeltlruiivfcroeiGrnggiWroiendsao,ccsthehplelh3sn0edrtmiuhcraiGnitncWtgiemleslaedwcuhpares3irn0eiogfmfdeeicanwtcitahvisme3.a0eFlsupmoeritrrnheioectdoirmmrwdeoeerpdreee,,rcatihonomeddmptwhuaaetxesaidmavl(eusFromiagrgueaecrmemopr2adl)xi.etiudImnd, uatehnmoedf

ftohlleowavienrga,gweemuasxeim“cuomveraamgep”liatundde“saomf pthlietuedfefe”cttoivseimgrpildy creelplsredseunritntgheeraecchor3d0emd ipnartaimmeetpeersr.iod were

comIpnuttheeds(tFuidgyuraere2a)., tIhnethpelafsomllaowbuinbgbl,ewoeccuusres “frceoqvueernatgley,”easnpdec“iaalmlypaltitnuidgeh”t,taonsdimcopullydrseipgrneisﬁecnatnttlhye

crheacnogrdeetdhepaarmampleittuerdse. and distribution days with a plasma bubble. In addition,

coverage of ionospheric GWs [31]. We therefore eliminated days with high geomagnetic activity (kp > 4) were excluded.

Remote Sens. 2019, 11, 2732 Remote Sens. 2019, 11, x FOR PEER REVIEW

4 of 11 4 of 11

Figure 2. Diagram for the computation of statistical parameters.

In the study area, the plasma bubble occurs frequently, especially at night, and could
significantly change the amplitude and distribution coverage of ionospheric GWs [31]. We therefore
eliminated days with a plasma bubble. In addition, days with high geomagnetic activity (kp > 4) were
excluded.
FFiigguurree 22.. DDiiaaggrraamm ffoorr tthhee ccoommppuuttaattiioonn ooff ssttaattiissttiiccaall ppaarraammeetteerrss..
3.3R. Reseuslutslts yseelsiihxleggiAylcaomihneglrwsutaishinIFfnrdnnsihtiiscanhFnegaiathidniogduetnngahw.drFgtunsdeeliineyrdmngest3asuewciiimytth3lnrpuawoesarirdprlewnFo3karyrsi,grereigekatsteshrnuhauneeatutrrdlnhasnsetenustedapiasd)lo3,.altneina,samoDrt)soneshnt.maphstxDoeyoepulaasirsxahmtnpmpbayuwedhlmupsdraieeeibtlsewcprrvhembsileadictlaeotnneoaThdfpt.droEsTImlatfbanChdmEupstuiaCmhbessaldnaetuvibsarnsddrenliemilebbeinerdyturuisitaeesosoobtrotsbinccbsoobmcuco,tlrnsuouedbbmaerrrsobcairmrsenoyrlveserlvsdhetvyfsroedriweotfgrerrrdoakbhiqoohtcgeyumhbkgicegeyeuetchhhnoMoorticertfmguhgoleaihyeGneudoar,octGggnnNshmnfetoeNrSsposetsoaSopmtSpepgmirSecshSneacteearhseagiMprctatotnhaiititlcucoieaeltovyitmnrrGauoiict,ccnrrWbhyatHea,ietrcwvcHsrKeto,tiocne[tiKrKw3voyridSeg1KTirhewtde]h,nyTid.pleiteoeW,n,(dtrtbkieein2paoypemb0nnln2tyo>1tbodho0ht4e4ts1tetere)hpp4(rc,doewehllo(ooft.weouhoelctArtoltrrheaihdeeeccsrliedalerl. 3Gtr.iiiWRgoognneseooswrusseppledthhrsebeeyrrdiiccentoGeGncWW-ttehsdsutrinwnigdeagerlerelsrmtedodoremtbneytchotnrseoi.dgnTi‐nihtnshi.suaImnnlldetemharensostfnoottrlhhlmoaswt.otiTrhniheggirisenansam.arIleenyafstnrihessqe, dutfhoaealynltsotlwotyhfioenercageccuahanrryraeielnyafgsrriewisoq,enudroeaesnydptslihyovefidroeieccadccGuhiW rnyrtieosnagr hiwghe-rFteihgduuinrvdeide3resdptorienrmsteon(hHtsigSah,n‐dtaheiuxlyanmdmeparlsxetiomorfmutmh(HunnSud,medrbaseitloryromf easﬀxtiermoctkuivemecgonurindmtscbeeplrlesor≥fh2eof0fu)e,rcltoriewvce-otgrhdruiednddcebelrylsstto≥hr2em0)lo,(LcloaSwl, ‐ ldigatihlyutnnmidnaegxrnismteoturwmmor(nkLuSam,ndbdaeirolyonfmoesﬀapxehciemtirviucemdgTrniEudCmcesbleelrsrie<osf2o0ebfafsenecdrtvi>ve0ed)g,brayindtdhcneloGlsnN-<t2hS0Sunasndtadetri>so0tno),,rHmanK(dNKnTSo,,ntihn‐teh2nu0u1n4md(eboretsrhtoeorfm yeﬀe(aeNrcsSti,hvteahdgersniidummcielballresrrieossfzueelfrtfsoe)).cdtDiavayeyssg. rwTiadibthcleeal1lsplliiasstszsmetrhaoeb)eudﬀbaebyclste.ivoTeranbhuliegmh1bgleiesrotosmfthaHgeSne,efLfteiScc,taaicvntedivnNiutySmwdbaeyrseoinfnoHetaSpc,lhoLtSyt,eeaadnr..d AsNSshdoawysninineaFcihguyreear3., thunderstorm events mainly occurred from March to September, while ionospheric GWs were detected in all mon(at)hThsu.ndTerhstiosrmmsteroakne s that there are frequently occurring ionospheric GWs trigg1e00r0e0d by non‐thunderstorm origins. In the following analysis, days of each year
were divided into high‐thunderstorm (HS, daily maximum number of effective grid cells ≥20), low‐
thunderstorm (LS, daily maximum number of effective grid cells <20 and >0), and non‐thunderstorm (NS, the number of effe5c00t0ive grid cells is zero) days. Table 1 lists the effective number of HS, LS, and
NS days in each year.

Counts

0 10000 0
6 5000
3

(a)Thunderstorm stroke

2

4

6

8

10

12

Month

(b)Ionospheric dTEC series

Counts TECU

0

0

0-3

2

4

6

8

Month

10

12

-6 60

(b)Ionospheric dTEC series

2

4

6

8

10

12

Month

3

Figure 3. (a) Thunderstorm stroke counts recorded by a local lightning network in 2014; (b) ionospheric dTEC series observed0 by the GNSS station, HKKT, in 2014.

TECU

-3

-6

0

2

4

6

8

10

12

Month

Remote Sens. 2019, 11, 2732

5 of 11

Table 1. Eﬀective numbers of high-thunderstorm (HS), low-thunderstorm (LS), and non-thunderstorm (NS) days in each year.

Year

HS

LS

NS

2014

88

54

105

2015

77

75

82

2016

105

92

99

2017

86

86

146

We compared the diurnal averages of thunderstorm activity (coverage) on HS and LS days to the diurnal averages of ionospheric GW activity (amplitude and coverage) on HS, LS, and NS days for each year. Diurnally averaged results are shown in Figure 4. Thunderstorm activity and ionospheric GW activity in diﬀerent circumstances were all plotted in one ﬁgure for quantitative comparison. The left vertical axis represents thunderstorm activity (blue line with points for HS and dark green line for LS); the right vertical axis indicates ionospheric GW activity (bold red, cyan, and magenta lines for HS, LS, and NS days, respectively) in each panel. The top two panels (Figure 4a,b) show the diurnally averaged thunderstorm coverage (Tarea) versus ionospheric GW amplitude (Gamp), and thunderstorm coverage versus ionospheric GW coverage (Garea), in 2014, respectively. The next two panels (Figure 4c,d) show the same plots in 2015. The next two panels (Figure 4e,f) show the same plots in 2016. The bottom two panels (Figure 4g,h) show the same plots in 2017. The average coverage of ionospheric GWs on all days during the ﬁrst and last 30 min time periods was signiﬁcantly smaller than that of adjacent time periods (Figure 4b,d,f,h). This is due to data missing for 20 min by the second-order center diﬀerence method.
As shown in Figure 4, thunderstorm activity presented an obvious diurnal variation on HS days for each year. Thunderstorm coverage slowly increased at the beginning, began to accelerate at about 11:00 LT (local time), and peaked at around 15:00–17:00 LT. From there, it began to decline with an approximately symmetric pattern. The diurnal variation curve is like a convex function. Thunderstorm activity on LS days showed a similar diurnal pattern to that of HS days, although coverage was very low. It was apparent that ionospheric GW activity also presented a diurnal variation on all days for each year.
The amplitude of ionospheric GWs showed diurnal variation with two signiﬁcant crests on HS days (four left panels in Figure 4). It presented a similar convex-function-like pattern to thunderstorm coverage in the rising phase at daytime. The ﬁrst crest was at around 18:00–19:00 LT, which had an approximate 2–3 h delay compared to that of thunderstorm coverage. At about 21:00–22:00 LT, the ionospheric GW amplitude rose again and reached the second crest at around 00:00–01:00 LT. After the time of the crest, it began to decline, with a trough at dawn. Coverage of ionospheric GWs showed a similar diurnal variation pattern to that of amplitude on HS days (four right panels in Figure 4).
The diurnal variation pattern of ionospheric GW activity on NS days, including amplitude (four left panels in Figure 4) and coverage (four right panels in Figure 4), was diﬀerent to that of thunderstorm activity. It increased quickly at the beginning and slowed down at about 11:00 LT, with an unobvious crest at around 16:00–18:00 LT, which was like a concave function. At night-time, ionospheric GW activity on NS days did not show night-time enhancement and kept declining until dawn. Due to a lack of night-time enhancement, the magnitude of ionospheric GW activity on NS days was signiﬁcantly smaller than that of HS and LS days.

RReemmoottee SSeennss.. 22001199,, 1111,, x27F3O2R PEER REVIEW

66ooff 1111

Garea/Counts

Gamp/TECU Tarea/Counts

Tarea/Counts

HS
40 (a)2014 30 20 10

LS

HS GWs

LS GWs

40 (b)2014

0.3

30

20

0.2

10

NS GWs 80 60 40 20

Garea/Counts

Tarea/Counts

Gamp/TECU Tarea/Counts

40 (c)2015

40 (d)2015

80

30

0.3

30

60

20

20

40

10

0.2

10

20

40 (e)2016

40 (f)2016

80

30

0.3

30

60

20

20

40

10

0.2

10

20

Garea/Counts

Gamp/TECU Tarea/Counts

Tarea/Counts

Garea/Counts

Gamp/TECU Tarea/Counts

Tarea/Counts

40 (g)2017

40 (h)2017

80

30

0.3

30

60

20

20

40

10

0.2

10

20

08 14 20 02 08 LT/h

08 14 20 02 08 LT/h

Figure 4. Diurnal averages of thunderstorm activity and ionospheric gravity wave (GW) activity from F2i0g1u4rteo42.0D17iu. rTnhael alivneerwagitehs pofotinhtusnrdeperrsetsoernmtsatchtuivnidtyerasntodrmionaoctsipvhiteyri(cblgureavfoitryHwSaavned(GdaWrk) agcrteievnityfofrroLmS). 2T0h1e4btold20l1in7.eTrheperleinseenwtsitihonpoosipnhtseriecpGreWsenatcstitvhiutynd(reerdstfoorrmHaSc,tcivyiatyn (fbolruLeSf,oarnHdSmaangdendtaarkfogrrNeeSn).f(oar) LTSh)u. nTdheersbtoolrdmlicnoevreerpargeesvenertsuiosnioonsposhpehriecriGcWGWacatimviptylit(uredde fionr2H01S4,;c(yba)nthfuorndLeSr,satnordmmcaogveenrtaagefovreNrsSu)s. (iao)noTshpuhnedriecrGstWormcovcoervaegreagine 2v0e1r4s;u(sc)itohnuonsdpehresrtiocrmGWcovaemrapgleituvedresuins i2o0n1o4s;p(hbe)ritchGunWdearmstpolritmudceoivne2ra0g15e; v(de)rstuhus nidoenrosstoprhmercicovGerWagecvoevresruasgieonionsp2h0e1r4i;c G(cW) tchouvnedraegrestionrm201c5o; v(ee)rathguendverstuosrmiocnoovseprhagereicveGrsWus
aiomnpolsiptuhdeericinGW201a5m; p(ldit)udtheuind2e0r1s6to; r(mf) thcouvnedreargsetorvmerscuosveiroangoesvpehresruics iGonWospchovereircaGgeWincov20er1a5g; e(ein)
t2h0u1n6;d(egr)sthourmndceorvsteorramgecvoevresruasgeiovneorspuhs eiorincoGspWhearmic pGlWituadme pinli2tu0d16e;i(nf)2t0h1u7;n(dhe)rtshtournmdecrostvoerrmagceoversaugse
ivoenrosuspshioenriocsGphWercicovGeWragcoeviner2a0g1e6i;n(2g0) 1t7h.understorm coverage versus ionospheric GW amplitude in
2017; (h) thunderstorm coverage versus ionospheric GW coverage in 2017. Ionospheric GW activity on LS days did not present a convex-function-like pattern at daytime
but hAasdsahoswimnilianrFvigauriraeti4o,nthautnndiegrhstt-otrimea,cctoivmitpyaprreedsetnoteiodnaonspohbevriiocuGs dWiuarcntaivl ivtayrioantioHnSodnaHysS.dTahyse fdoaryetaimcheyreisairn.gThpuatntderenrsotonrLmScdoavyesratogoeksltohwe laypipnrcorxeiamseadteasththape ebeogfiannliinneg,,ebsepgeacniatlloyafcocrecleorvaetreaagtea(bfoutr 1ri1g:0h0t pLaTn(elloscianl Ftiimguer)e, a4n).dApmeapklietdudaet ianro2u0n14da1n5d:0200–157:w00asLaT.mFirxotmurethoefrceo, nitcabveegaanndtocodnevcleixnefuwnictthioanns a(Fpipgruorxei4maa,ct)e.lyGesnyemramllye,tirticseepmaettdertno.beTahemidxeiudrenﬀael ctv, abryiatthioenconcvuerxv-efunisctiloikne-likae pcaotntevrenxoffuHnSctdioayns. Tanhdunthdeercsotnocramvea-cfutinvcityiono-nlikLeSvdaariyastisohnopwaettderan soifmNilSardadyisu.rInnatlerpeasttienrgnlyt,othtehqatuaonf tHitaStidveayvsa,luaeltshoofutghhe cionvoersapgheerwicaGs Wvearcytilvoiwty. cIotuwldafsuartphpearrseunpt ptohrattthioeneoxsipstheenrcice oGfWa maicxteivditeyﬀeaclst.oApsreshseonwtendinaFdigiuurrnea4l, vthaerimataiognniotundaellodfaioynsofsoprheearcihc GyeWara. mplitude and coverage on NS days was larger than that of HS days at daTyhtiemaem. pTlhiteumdeagonf itoundoespohf eiorincoGspWhsersihcoGwWedadmiuprlintauldveaarniadtioconvweriathgetwono sLiSgndiafiycsanwtacsrestssenontiaHllSy days (four left panels in Figure 4). It presented a similar convex‐function‐like pattern to thunderstorm coverage in the rising phase at daytime. The first crest was at around 18:00–19:00 LT, which had an

Remote Sens. 2019, 11, 2732

7 of 11

between the magnitude of ionospheric GW activity on NS days and the magnitude of ionospheric GW

on HS days. At night-time, the amplitude of ionospheric GWs on LS days had a similar enhancement

to that of HS days, although the magnitude was smaller. Coverage on LS days in 2014 and 2015 did

not have an observable crest, but the magnitude was enhanced during this time according to the

variation trend.

To compare the annual variation of thunderstorm activity and ionospheric GW activity, we

replotted Figure 4 and put each activity during the four years into one ﬁgure. The corresponding

results are shown in Figure 5. Diﬀerent colors in all panels distinguish the years. The top two panels

(Figure 5a,b) show the diurnally averaged thunderstorm coverage on HS and LS days, respectively.

Figure 5c presents the diurnally averaged amplitude of ionospheric GWs on HS days, and Figure 5d

presents the diurnally averaged coverage of ionospheric GWs on HS days. The next two panels

(Figure 5e,f) show the same plots for ionospheric GWs on LS days. The bottom two panels (Figure 5g,h)

shoRwemottheeSesnasm. 20e1p9,l1o1t,sxfFoOrRioPnEoEsRpRhEeVrIiEcWGWs on NS days.

8 of 11

40 (a)HS

2014

2015

2016 20 (b)LS

2017

20

10

Tarea/Counts

Tarea/Counts

Gamp/TECU

0 (c)HS
0.3
0.2
(e)LS 0.3

Garea/Counts

0 (d)HS
60
40
20
0 (f)LS
60

Garea/Counts

Gamp/TECU

0.2 (g)NS
0.3
0.2

40
20
0 (h)NS
60
40
20

Garea/Counts

Gamp/TECU

08 14 20 02 08 LT/h

0 08 14 20 02 08
LT/h

Figure 5. Annual variation of diurnal averages of thunderstorm activity and ionospheric GW activity. (a) FTihguunred5er. sAtonrnmuaclovvaerriaagtieononofHdSiudranyasl;a(vbe)rtahguens doef rtshtuonrmdecrsotvoerrmagaectoivnitLySadnadyiso;n(co)spiohneorsicpGheWricacGtiWvity. ma(gan) iTtuhduendoenrsHtoSrmdacyosv; e(dra)gioenoonspHhSerdiacyGs;W(bc)othvuernadgeersotnorHmScdoavyersa; g(ee)oinonLoSspdhaeyrsi;c(cG)WionmosapghneitruicdGe W on mLSagdnaiytusd; e(fo) nioHnoSsdpahyesr;ic(dG) Wioncoosvpehreargice GonWLcSovdearyasg;e(go)niHonSodspayhse;r(iec)GioWnomspahgenriitcuGdWe omn aNgSnidtuadyes;on (h)LioSndoasypsh;e(rfi)c iGonWoscpohveerriacgGe Won cNoSvedraaygse. on LS days; (g) ionospheric GW magnitude on NS days; (h)
ionospheric GW coverage on NS days.

4. Discussion
Daytime ionospheric GW activity on HS days showed a similar convex‐function‐like diurnal variation pattern to thunderstorm activity, with a 2–3 hour delay in crest time. Theoretical research

Remote Sens. 2019, 11, 2732

8 of 11

As seen in Figure 5, the magnitude of thunderstorm activity and ionospheric GW activity varied each year. The magnitude of thunderstorm activity on HS days, especially during the daytime, is positively related to solar activity (Figure 5a). This correlation is unapparent on LS days (Figure 5b), which can be attributed to daytime classiﬁcation. The magnitude of ionospheric GW activity on HS days has a signiﬁcant positive correlation with solar activity during nearly all times of the day (Figure 5c). For example, the magnitude of ionospheric GWs at night-time in 2016 (low solar activity) was lower than that of 2014 and 2015 (high solar activity), although the magnitude of thunderstorm activity reverses (Figure 5a). In addition, the magnitude of ionospheric GW activity on LS days had an approximately positive correlation to solar activity, although there was no apparent correlation to thunderstorm activity (Figure 5e,f). This solar activity correlation was also suitable for ionospheric GW activity on NS days (Figure 5g,h).
4. Discussion
Daytime ionospheric GW activity on HS days showed a similar convex-function-like diurnal variation pattern to thunderstorm activity, with a 2–3 h delay in crest time. Theoretical research shows that GWs generated in the lower atmosphere need one to several hours (depending on period and speed) to propagate upwards to an ionospheric height of 350 km [32]. As such, a 2–3 h delay between ionospheric GW activity and thunderstorm activity is reasonable. In addition, ionospheric GW activity on NS days showed a concave-function-like variation pattern, especially in the rising phase, compared to that of thunderstorm activity. The above results provide evidence that daytime ionospheric GW activity on HS days is triggered by thunderstorm activity.
Daytime ionospheric GW activity on LS days presents an approximately linear diurnal variation before the time of the crest, which is diﬀerent to that of HS days. Furthermore, the magnitude of ionospheric GW activity (amplitude and coverage) on LS days is generally larger than that of HS days. The diﬀerence in background TEC is unlikely to cause abnormal phenomena when LS days have more signiﬁcant ionospheric GW activity than HS days in the daytime. Both HS activity and LS activity occur in spring and summer (Figure 3a), so background TEC is not signiﬁcantly diﬀerent. In addition, normal results at night-time can support this explanation. Considering the magnitude and diurnal variation patterns of ionospheric GW activity on NS days, it is very likely that multiple origins, including thunderstorm activity and non-thunderstorm activity, trigger ionospheric GWs on LS days. Thunderstorm-related GWs are mixed with non-thunderstorm-related GWs, and the original variation pattern is changed.
The diurnal variation of daytime ionospheric GW activity on LS days in this study area is diﬀerent to that of South America in the southern hemisphere, as reported by Lay [24]. According to Lay [24], the diurnal variation pattern on LS days is similar to that on HS days, and the magnitude of the former is always smaller than that of the latter. These distinct results may be attributed to the diﬀerent study areas and employed data. Ionospheric GW activity may show diﬀerences based on longitudinal diﬀerences between the southern region of China and South America, even though both are in low latitudes. In addition, only 11 days were used by Lay to calculate the diurnal averages of thunderstorm activity on LS days [24], while over 50 days in each year were used in this study (Table 1). Non-thunderstorm events during the 11 LS days in Lay’s study [24] could have been unnoticeable, so the original variation pattern was preserved.
Night-time enhancement of ionospheric GW activity on both HS and LS days seems to be related to thunderstorm activity because it only occurs on thunderstorm days and is correlated with the magnitude of thunderstorm activity. The night-time enhancement of ionospheric GW activity on thunderstorm days was also reported in other areas, including low-latitude areas in South America [24] and the Great Plains in the United States, which is a middle-latitude area of the northern hemisphere [18]. These observed results are conducive to intensifying the relationship between night-time enhancement of ionospheric GW activity and thunderstorm activity. The mechanism for night-time enhancement is still unclear at present. One possible cause is thunderstorm-induced Perkins instability [33], a common

Remote Sens. 2019, 11, 2732

9 of 11

factor in triggering medium-scale ionospheric disturbances at night. In addition, electric changes produced by lightning discharges could be a factor in the anomaly [34].
The variation of ionospheric GW activity on both HS and LS days showed a positive correlation to solar activity. This is in line with the fact that ionospheric TEC is high in a high solar activity year and vice versa, and the magnitude of ionospheric GWs is proportional to background TEC [5]. Furthermore, thunderstorm activity had a positive correlation with solar activity. Thus, it is the combined action of solar and thunderstorm activity that drives the variations in ionospheric GWs.
5. Conclusions
In this study, we investigated the amplitude and coverage of thunderstorm-related ionospheric GWs above the low latitude of the northern hemisphere using GNSS ionospheric TEC and lightning data observed in Southern China from the high solar activity year of 2014 to the low solar activity year of 2017. Compared to previous studies, the statistical analysis of long-term data in this study provided support for the correlation between thunderstorm events and ionospheric GWs. Furthermore, the observed results show that the diurnal variations of daytime ionospheric GW activity on LS days in this study area are diﬀerent to those on HS days. It is very likely that thunderstorm-induced ionospheric GWs are mixed with other non-thunderstorm activity on LS days. Ionospheric GW activity on thunderstorm days shows a positive correlation to solar activity, which can be attributed to the combined action of solar and thunderstorm activity. These ﬁndings provide more comprehensive knowledge about thunderstorm–ionosphere coupling in low-latitude areas. In the future, more work is needed to discuss the night-time enhancement mechanism of thunderstorm-related ionospheric GWs and the mixed eﬀect of ionospheric GWs on LS days in other study areas. In addition, more GNSS satellites, including BDS and Galileo, could be employed to further increase data coverage for more reliable results.
Author Contributions: L.T. led the writing of the manuscript and contributed to data analysis and research design. W.C. supervised this study, contributed to the research design, and served as the corresponding author. M.C. contributed to data processing. O.-P.L. contributed to manuscript editing.
Funding: This research was funded by National Natural Science Foundation of China (Grant number 41804021), Natural Science Foundation of Guangdong Province (Grant number 2017A030310242), National Key Research and Development Program of China (grant number No. 2016YFB0501803) and the Hong Kong RGC Joint Research Scheme (E-PolyU501/16).
Acknowledgments: The GNSS data used in this study are provided by the Geodetic Survey Section of Survey and Mapping Oﬃce in Hong Kong.
Conﬂicts of Interest: The authors declare no conﬂict of interest.
References
1. Nappo, C.J. An Introduction to Atmospheric Gravity Waves; Academic Press: San Diego, CA, USA, 2002. 2. 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 ﬁltering. J. Geophys. Res. Atmos. 2003, 108. [CrossRef] 3. Nakamura, T.; Aono, T.; Tsuda, T.; Admiranto, A.G.; Achmad, E. Mesospheric gravity waves over a tropical convective region observed by OH airglow imaging in Indonesia. Geophys. Res. Lett. 2003, 30. [CrossRef] 4. Hines, C.O. Internal atmospheric gravity waves at ionospheric heights. Can. J. Phys. 1960, 38, 1441–1481. [CrossRef] 5. Hooke, W.H. Ionospheric irregularities produced by internal atmospheric gravity waves. J. Atmos. Terr. Phys. 1968, 30, 795–823. [CrossRef] 6. Rolland, L.M.; Lognonné, P.; Astafyeva, E.; Kherani, E.A.; Kobayashi, N.; Mann, M.; Munekane, H. The resonant response of the ionosphere imaged after the 2011 oﬀ the Paciﬁc coast of Tohoku Earthquake. Earth Planets Space 2011, 63, 853–857. [CrossRef]

Remote Sens. 2019, 11, 2732

10 of 11

7. Savastano, G.; Komjathy, A.; Verkhoglyadova, O.; Mazzoni, A.; Crespi, M.; Wei, Y.; Mannucci, A.J. Real-Time Detection of Tsunami Ionospheric Disturbances with a Stand-Alone GNSS Receiver: A Preliminary Feasibility Demonstration. Sci. Rep.-UK 2017, 7, 46607. [CrossRef]
8. Rakoto, V.; Lognonné, P.; Roll, L.; Coïsson, P. Tsunami wave height estimation from GPS-derived ionospheric data. J. Geophys. Res. Space Phys. 2018, 123, 4329–4348. [CrossRef]
9. Tang, L.; Li, Z.; Zhou, B. Large-area tsunami signatures in ionosphere observed by GPS TEC after the 2011 Tohoku earthquake. GPS Solut. 2018, 22, 93. [CrossRef]
10. Chou, M.Y.; Lin, C.H.; Yue, J.; Tsai, H.F.; Sun, Y.Y.; Liu, J.Y.; Chen, C.H. Concentric traveling ionosphere disturbances triggered by Super Typhoon Meranti (2016). Geophys. Res. Lett. 2017, 44, 1219–1226. [CrossRef]
11. Stull, R.B. Internal gravity waves generated by penetrative convection. J. Atmos. Sci. 1976, 33, 1279–1286. [CrossRef]
12. Plougonven, R.; Zhang, F. Internal gravity waves from atmospheric jets and fronts. Rev. Geophys. 2014, 52, 33–76. [CrossRef]
13. Pierce, A.D.; Coroniti, S.C. A mechanism for the generation of acoustic-gravity waves during thunderstorm formation. Nature 1966, 210, 1209. [CrossRef]
14. Yue, J.; Vadas, S.L.; She, C.Y.; Nakamura, T.; Reising, S.C.; Liu, H.L.; Stamus, P.; Krueger, D.A.; Lyons, W.; Li, T. Concentric gravity waves in the mesosphere generated by deep convective plumes in the lower atmosphere near Fort Collins, Colorado. J. Geophys. Res. Atmos. 2009, 114. [CrossRef]
15. Lay, E.H.; Shao, X.M. Multi-station probing of thunderstorm-generated D-layer ﬂuctuations by using time-domain lightning waveforms. Geophys. Res. Lett. 2011, 38. [CrossRef]
16. Chou, C.C.; Dai, J.; Kuo, C.L.; Huang, T.Y. Simultaneous observations of storm-generated sprite and gravity wave over Bangladesh. J. Geophys. Res. Space Phys. 2016, 121, 9222–9233. [CrossRef]
17. Vadas, S.L.; Liu, H.L. Numerical modeling of the large-scale neutral and plasma responses to the body forces created by the dissipation of gravity waves from 6 h of deep convection in Brazil. J. Geophys. Res. Space Phys. 2013, 118, 2593–2617. [CrossRef]
18. Lay, E.H.; Shao, X.M.; Carrano, C.S. Variation in total electron content above large thunderstorms. Geophys. Res. Lett. 2013, 40, 1945–1949. [CrossRef]
19. Azeem, I.; Yue, J.; Hoﬀmann, L.; Miller, S.D.; Straka, W.C.; Crowley, G. Multisensor proﬁling of a concentric gravity wave event propagating from the troposphere to the ionosphere. Geophys. Res. Lett. 2015, 42, 7874–7880. [CrossRef]
20. Nishioka, M.; Tsugawa, T.; Kubota, M.; Ishii, M. Concentric waves and short-period oscillations observed in the ionosphere after the 2013 Moore EF5 tornado. Geophys. Res. Lett. 2013, 40, 5581–5586. [CrossRef]
21. Yue, J.; Hoﬀmann, L.; Joan Alexander, M. Simultaneous observations of convective gravity waves from a ground-based airglow imager and the AIRS satellite experiment. J. Geophys. Res. Atmos. 2013, 118, 3178–3191. [CrossRef]
22. 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. [CrossRef]
23. Lay, E.H.; Shao, X.M.; Kendrick, A.K.; Carrano, C.S. Ionospheric acoustic and gravity waves associated with midlatitude thunderstorms. J. Geophys. Res. Space Phys. 2015, 120, 6010–6020. [CrossRef]
24. Lay, E.H. Ionospheric Irregularities and Acoustic/Gravity Wave Activity Above Low-Latitude Thunderstorms. Geophys. Res. Lett. 2018, 45, 90–97. [CrossRef]
25. Padokhin, A.M.; Tereshin, N.A.; Yasyukevich, Y.V.; Andreeva, E.S.; Nazarenko, M.O.; Yasyukevich, A.S.; Kozlovtseva, E.A.; Kurbatov, G.A. Application of BDS-GEO for studying TEC variability in equatorial ionosphere on diﬀerent time scales. Adv. Space Res. 2019, 63, 257–269. [CrossRef]
26. Savastano, G.; Komjathy, A.; Shume, E.; Vergados, P.; Ravanelli, M.; Verkhoglyadova, O.; Meng, X.; Crespi, M. Advantages of Geostationary Satellites for Ionospheric Anomaly Studies: Ionospheric Plasma Depletion Following a Rocket Launch. Remote Sens. 2019, 11, 1734. [CrossRef]
27. Zhang, X.; Tang, L. Detection of ionospheric disturbances driven by the 2014 Chile tsunami using GPS total electron content in New Zealand. J. Geophys. Res. Space Phys. 2015, 120, 7918–7925. [CrossRef]
28. Schaer, S. Mapping and Predicting the Earth’s Ionosphere Using the Global Positioning System; Geod.-Geophys. Arb. Schweiz: Zurich, Switzerland, 1999.

Remote Sens. 2019, 11, 2732

11 of 11

29. Tang, L.; Zhang, X.; Li, Z. Observation of ionospheric disturbances induced by the 2011 Tohoku tsunami using far-ﬁeld GPS data in Hawaii. Earth Planets Space 2015, 67, 88. [CrossRef]
30. Zhang, X.; Tang, L. Traveling ionospheric disturbances triggered by the 2009 North Korean underground nuclear explosion. Ann. Geophys.-Ger. 2015, 33, 137–142. [CrossRef]
31. Ji, S.; Chen, W.; Wang, Z.; Xu, Y.; Weng, D.; Wan, J.; Fan, Y.; Huang, B.; Fan, S.; Sun, G. A study of occurrence characteristics of plasma bubbles over Hong Kong area. Adv. Space Res. 2013, 52, 1949–1958. [CrossRef]
32. 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. [CrossRef]
33. Perkins, F. Spread F and ionospheric currents. J. Geophys. Res. 1973, 78, 218–226. [CrossRef] 34. 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. 2013, 6, 29–33. [CrossRef]
© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

