<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Phys.</journal-id>
<journal-title>Frontiers in Physics</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Phys.</abbrev-journal-title>
<issn pub-type="epub">2296-424X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fphy.2020.00348</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Physics</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Short-Term Lightning Response to Ground Level Enhancements</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Wu</surname> <given-names>Qiong</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Li</surname> <given-names>Hui</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Wang</surname> <given-names>Chi</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>State Key Laboratory of Space Weather, National Space Science Center, Chinese Academy of Sciences</institution>, <addr-line>Beijing</addr-line>, <country>China</country></aff>
<aff id="aff2"><sup>2</sup><institution>University of Chinese Academy of Sciences</institution>, <addr-line>Beijing</addr-line>, <country>China</country></aff>
<aff id="aff3"><sup>3</sup><institution>School of Astronomy and Space Science, University of Chinese Academy of Sciences</institution>, <addr-line>Beijing</addr-line>, <country>China</country></aff>
<aff id="aff4"><sup>4</sup><institution>College of Earth and Planetary Sciences, University of Chinese Academy of Sciences</institution>, <addr-line>Beijing</addr-line>, <country>China</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Ioannis A. Daglis, National and Kapodistrian University of Athens, Greece</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Athanasios Papaioannou, National Observatory of Athens, Greece; Ilya Usoskin, University of Oulu, Finland; Konstantinos Kourtidis, Democritus University of Thrace, Greece</p></fn>
<corresp id="c001">&#x0002A;Correspondence: Hui Li <email>hli&#x00040;spaceweather.ac.cn</email></corresp>
<fn fn-type="other" id="fn001"><p>This article was submitted to Space Physics, a section of the journal Frontiers in Physics</p></fn></author-notes>
<pub-date pub-type="epub">
<day>04</day>
<month>09</month>
<year>2020</year>
</pub-date>
<pub-date pub-type="collection">
<year>2020</year>
</pub-date>
<volume>8</volume>
<elocation-id>348</elocation-id>
<history>
<date date-type="received">
<day>10</day>
<month>05</month>
<year>2020</year>
</date>
<date date-type="accepted">
<day>23</day>
<month>07</month>
<year>2020</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2020 Wu, Li and Wang.</copyright-statement>
<copyright-year>2020</copyright-year>
<copyright-holder>Wu, Li and Wang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract><p>Cosmic rays (CRs) are considered the primary energetic particle source of atmospheric ionization on Earth. Under the modulation of severe solar eruption events, CR variations are further speculated to impact the Earth&#x00027;s lightning activities. Previous researches show that CR intensity and lightning incidence are positively correlated on the time scale of several days to decades. However, to our knowledge, the global lightning response to short-term CR variation has not been studied in the literature. Ground level enhancements (GLEs) provide the opportunity to study such a possible link. As a small fraction of solar energetic particle events that could reach the energy level of several GeVs, GLEs can thus generate atmospheric cascades that could be recorded by ground-based neutron monitors. Furthermore, as GLEs generally take place within several 10 min to an hour, the lightning variations caused by potential meteorological factors could be maximally diminished in such a short time. During the operational period of the World Wide Lightning Location Network (Aug 2004 to now), three typical GLEs with the intensity &#x0003E;15% are analyzed from the International GLE Database, namely &#x00023;69 (Jan 20, 2005), &#x00023;70 (Dec 13, 2006), and &#x00023;71 (May 17, 2012). For each GLE event, the global lightning incidence presents a positive response to GLE (i.e., a significant enhancement within 20 min right after the GLE onset). Meanwhile, the relative amplitude of lightning response seems to be in direct proportion to GLE intensity (i.e., the more intensive the GLE is, the more obvious the increase in the lightning incidence is), which is further verified to be statistically significant by Monte Carlo test. By comparing lightning responses in different latitudinal zones, we find that more intensive lightning responses to GLEs seem to be at higher latitudes.</p></abstract>
<kwd-group>
<kwd>lightning incidence</kwd>
<kwd>ground level enhancement (GLE)</kwd>
<kwd>solar energetic particle (SEP)</kwd>
<kwd>cosmic ray</kwd>
<kwd>short-term</kwd>
<kwd>latitude effects</kwd>
</kwd-group>
<contract-num rid="cn001">41574159</contract-num>
<contract-num rid="cn001">41574169</contract-num>
<contract-num rid="cn001">41731070</contract-num>
<contract-num rid="cn001">41874203</contract-num>
<contract-sponsor id="cn001">National Natural Science Foundation of China<named-content content-type="fundref-id">10.13039/501100001809</named-content></contract-sponsor>
<contract-sponsor id="cn002">China Association for Science and Technology<named-content content-type="fundref-id">10.13039/100010097</named-content></contract-sponsor>
<contract-sponsor id="cn003">Youth Innovation Promotion Association of the Chinese Academy of Sciences<named-content content-type="fundref-id">10.13039/501100004739</named-content></contract-sponsor>
<counts>
<fig-count count="4"/>
<table-count count="2"/>
<equation-count count="0"/>
<ref-count count="31"/>
<page-count count="7"/>
<word-count count="4264"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>Introduction</title>
<p>The influence of space weather variation on the Earth&#x00027;s atmospheric electrical environment has consistently captured attention in recent years. Cosmic rays (CRs), including Galactic Cosmic Rays (GCRs) and Solar Energetic Particles (SEPs), are considered the primary source of energetic particle originating from the extraterrestrial space and impinging upon the Earth&#x00027;s atmosphere [<xref ref-type="bibr" rid="B1">1</xref>]. GCRs, modulated by solar activities, are found to be the main source of atmospheric ionization, which continually interact with atomic nuclei in the air and generate secondary particles [<xref ref-type="bibr" rid="B2">2</xref>]. A small fraction of SEP events, which are identified as Ground Level Enhancements (GLEs), tend to affect atmospheric ionization by sporadic penetrations into the atmosphere during explosive energy releases of severe solar activities.</p>
<p>A GLE event generally possesses the energy of several GeV and generates atmospheric cascades which could be recorded by ground-based neutron monitor (NM) stations [<xref ref-type="bibr" rid="B3">3</xref>]. A series of studies by Mishev et al. [<xref ref-type="bibr" rid="B4">4</xref>&#x02013;<xref ref-type="bibr" rid="B7">7</xref>] demonstrate that GLEs, especially the most extreme ones, have a significantly positive effect on the short-time variation of the atmospheric ionization at high-altitudes of around 12&#x02013;18 km in polar and sub-polar regions. Moderate or even negative responses are found at middle- and low-latitudes due to the accompanying Forbush Decreases (FDs). Similar results are found by Usoskin et al. [<xref ref-type="bibr" rid="B8">8</xref>] and Velinov et al. [<xref ref-type="bibr" rid="B9">9</xref>]. In addition to atmospheric ionization, the major GLEs are also claimed to result in a larger air conductivity and further create significant variations in the strength of the electric field, especially in its vertical component [<xref ref-type="bibr" rid="B10">10</xref>&#x02013;<xref ref-type="bibr" rid="B14">14</xref>]. Recently, Golubenko et al. [<xref ref-type="bibr" rid="B15">15</xref>] find that an extreme SEP event could lead to a large increase in fair-weather downward current density on a global scale. Such variations do not only exist in the upper atmosphere above thunderstorms, but also extend down to the ground level. Changes in atmospheric ionization and electric field resulting from GLEs may eventually apply to the global electric circuit by inferred effects on cloud microphysics, thunderstorms, and other climate parameters, etc. [<xref ref-type="bibr" rid="B16">16</xref>&#x02013;<xref ref-type="bibr" rid="B20">20</xref>].</p>
<p>Lightning is a proxy for the variation in the Earth&#x00027;s atmospheric electrical environment. Although primarily maintained by meteorological factors, lightning is also speculated to be positively modulated by GCRs [<xref ref-type="bibr" rid="B21">21</xref>&#x02013;<xref ref-type="bibr" rid="B24">24</xref>]. Most studies mainly focused on the timescale of several days to several decades, especially on FDs [<xref ref-type="bibr" rid="B25">25</xref>, <xref ref-type="bibr" rid="B26">26</xref>]. However, to our knowledge, lightning response to short-term CR variation has not been studied in the literature. GLEs provide a good opportunity to focus on such a possible link between CRs and atmospheric processes. On one hand, GLEs present the increasing phases of CR variation, which could be in contrast with the decreasing phases presented by FDs. On the other hand, as GLEs generally occur within several 10 min to an hour, the potential effects of meteorological factors on lightning could be maximally diminished within such a short time.</p>
<p>In this study, we focus on the responses of lightning incidence to SEP intensity during GLE events. The GLE events are selected from the International GLE database and the lightning data are obtained from the World Wide Lightning Location Network (WWLLN). Section 2 introduces the data sets and methodology. Section 3 gives the results of the lightning responses to GLEs and Monte Carlo tests, including some brief discussions. Section conclusion summarizes our work.</p>
</sec>
<sec id="s2">
<title>Data Sets and Methodology</title>
<p>The WWLLN<xref ref-type="fn" rid="fn0001"><sup>1</sup></xref> has been continuously providing global lightning detection at ground level since August 2004. Being operated by the University of Washington (USA) and the University of Otago (New Zealand), it currently has over 70 sensors which are generally uniformly distributed at distances over thousands of kilometers around the globe. These sensors identify lightning transients and process sferic wave packets in the very-low-frequency band (100 Hz&#x02212;24 kHz) [<xref ref-type="bibr" rid="B27">27</xref>]. Using Time of Group Arrival (TOGA) analysis, each lightning stroke location requires the TOGA from at least five sensors [<xref ref-type="bibr" rid="B28">28</xref>]. Typically, only about 15&#x02013;30% of lightning strokes detected by one sensor are detected by five or more [<xref ref-type="bibr" rid="B29">29</xref>]. In addition, the detection efficiency of the WWLLN is in direct proportion to the peak current intensity of lighting return strokes [<xref ref-type="bibr" rid="B30">30</xref>], which is now up to 80% of the strokes above 50 kA peak current. Thus, the following results are mainly valid for the lightning strokes with the peak current above 50 kA.</p>
<p>The International GLE database<xref ref-type="fn" rid="fn0002"><sup>2</sup></xref> collects and archives neutron counting rates concerning GLEs from the worldwide NM network [<xref ref-type="bibr" rid="B31">31</xref>]. The NM network has been successively monitoring GLEs since the 1950&#x00027;s. The first and strongest event detected on Feb 23, 1956 was recorded as GLE number 5 (&#x00023;5). The following events are numbered consecutively. There are five identified GLEs during the operational period of the WWLLN, namely &#x00023;68 (Jan 17, 2005), &#x00023;69 (Jan 20, 2005), &#x00023;70 (Dec 13, 2006), &#x00023;71 (May 17, 2012), and &#x00023;72 (Sep 10, 2017). In this study, only the GLEs satisfying the following criteria are analyzed:</p>
<list list-type="order">
<list-item><p>The largest enhancement in neutron counting rate observed by some NM station during the GLE event is required to be &#x0003E;15%.</p></list-item>
<list-item><p>The increasing phase of the GLE event is evident, and its rising time is required to be &#x0003C;30 min.</p></list-item>
</list>
<p>Thus, three GLEs, &#x00023;69 (Jan 20, 2005), &#x00023;70 (Dec 13, 2006), and &#x00023;71 (May 17, 2012), are selected and analyzed in this study. The largest enhancement in neutron counting rate is 4808.95% (&#x00023;69, observed by SOPB), 92.1% (&#x00023;70, observed by OULU), and 17.3% (&#x00023;71, observed by SOPO), respectively. The corresponding rising time is 10, 20, and 20 min, respectively.</p>
<p>To present the general variation characteristics of a GLE event, we further use the average neutron counting rate of several NM stations which are close to the areas with a high incidence of lightning during each GLE. The GLE intensity, defined as the average enhancement in neutron counting rate observed by these selected NM stations, is 40.0% (&#x00023;69), 9.5% (&#x00023;70), and 3.9% (&#x00023;71), respectively. Since both the lightning distribution and NM stations in operation change over time, the stations used for each GLE event are different, as listed in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption><p>NM stations used to determine the average SEP variation for the three GLEs.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th valign="top" align="left"><bold>No</bold>.</th>
<th valign="top" align="left"><bold>Date</bold></th>
<th valign="top" align="left"><bold>Station</bold></th>
<th valign="top" align="center"><bold>Latitude (<sup><bold>&#x000B0;</bold></sup>)</bold></th>
<th valign="top" align="center"><bold>Longitude (<sup><bold>&#x000B0;</bold></sup>)</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">69</td>
<td valign="top" align="left">Jan 20, 2005</td>
<td valign="top" align="left">ATHN</td>
<td valign="top" align="center">37.97</td>
<td valign="top" align="center">23.72</td>
</tr>
<tr>
<td/>
<td/>
<td valign="top" align="left">KGSN</td>
<td valign="top" align="center">&#x02212;42.99</td>
<td valign="top" align="center">147.29</td>
</tr>
<tr>
<td/>
<td/>
<td valign="top" align="left">SNTG</td>
<td valign="top" align="center">&#x02212;33.49</td>
<td valign="top" align="center">&#x02212;70.72</td>
</tr>
<tr>
<td valign="top" align="left">70</td>
<td valign="top" align="left">Dec 13, 2006</td>
<td valign="top" align="left">KING</td>
<td valign="top" align="center">&#x02212;42.99</td>
<td valign="top" align="center">147.29</td>
</tr>
<tr>
<td/>
<td/>
<td valign="top" align="left">PWNK</td>
<td valign="top" align="center">55.00</td>
<td valign="top" align="center">&#x02212;85.00</td>
</tr>
<tr>
<td/>
<td/>
<td valign="top" align="left">ROME</td>
<td valign="top" align="center">41.90</td>
<td valign="top" align="center">12.52</td>
</tr>
<tr>
<td/>
<td/>
<td valign="top" align="left">STGO</td>
<td valign="top" align="center">&#x02212;33.49</td>
<td valign="top" align="center">&#x02212;70.72</td>
</tr>
<tr>
<td valign="top" align="left">71</td>
<td valign="top" align="left">May 17, 2012</td>
<td valign="top" align="left">APTY</td>
<td valign="top" align="center">67.55</td>
<td valign="top" align="center">33.33</td>
</tr>
<tr>
<td/>
<td/>
<td valign="top" align="left">FSMT</td>
<td valign="top" align="center">60.00</td>
<td valign="top" align="center">&#x02212;112.00</td>
</tr>
<tr>
<td/>
<td/>
<td valign="top" align="left">IRK2</td>
<td valign="top" align="center">52.28</td>
<td valign="top" align="center">104.02</td>
</tr>
<tr>
<td/>
<td/>
<td valign="top" align="left">KGSN</td>
<td valign="top" align="center">&#x02212;42.99</td>
<td valign="top" align="center">147.29</td>
</tr>
<tr>
<td/>
<td/>
<td valign="top" align="left">MXCO</td>
<td valign="top" align="center">19.33</td>
<td valign="top" align="center">&#x02212;99.18</td>
</tr>
<tr>
<td/>
<td/>
<td valign="top" align="left">NWRK</td>
<td valign="top" align="center">39.68</td>
<td valign="top" align="center">&#x02212;75.75</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Relative variations of lightning incidence in 1.5 h for these three GLEs are investigated, from 20 min before the GLE onset to 70 min later than that. During this period, relative variations of lightning incidence are obtained by dividing the number of lightning strokes in every 10 min by the value in the first 10-min after the GLE onset. Similarly, relative variations of the neutron counting rate are obtained with a temporal resolution of 5 min. To maximally reduce the potential effects of meteorological factors, we focus on the lightning response within 20 min after the GLE onset and define the value in the second 10-min after the GLE onset as the lightning response intensity.</p>
</sec>
<sec sec-type="results" id="s3">
<title>Results</title>
<sec>
<title>Case Analysis</title>
<p><xref ref-type="fig" rid="F1">Figure 1A</xref> shows the results during GLE &#x00023;69. The relative variation in SEP, represented by the red line, rapidly increases to its peak value within 15 min after the GLE onset at 06:55, giving the GLE intensity of 40.0%. After the GLE onset, the lightning incidence histogram also presents a significant enhancement within 20 min, giving the lightning response intensity of 1.37. During the SEP recovery, the lightning incidence presents some relatively minor fluctuations. This indicates that short-term SEP intensity enhancement during a GLE event may lead to a sharp and timely increase in global lightning incidence.</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p><bold>(A)</bold> Lightning response during GLE &#x00023;69 (Jan 20, 2005). The red line represents the relative change of 5-min neutron counting rates, the gray bars represent the relative variation of 10-min lightning incidence. <bold>(B)</bold> Results of the Monte Carlo test in 2005. The total number of tests is 9,000. The blue solid line represents the mean value, the blue dotted line represents the value at 95% confidence level, and the red dashed line represents the observed lightning response intensity.</p></caption>
<graphic xlink:href="fphy-08-00348-g0001.tif"/>
</fig>
<p>To further verify the reliability of the above lightning response to GLE &#x00023;69, a Monte Carlo test is conducted as follows. Firstly, we randomly select a 20-min interval in 2005. Secondly, we determine the relative change of global lightning in this interval. Thirdly, we repeat such processes by 9,000 times, getting the histogram distribution of relative change of lightning incidence in 2005, as shown in <xref ref-type="fig" rid="F1">Figure 1B</xref>. The relative change in lightning incidence is generally between 0.6 and 1.4, with both the peak and mean value close to unity, presenting a normal distribution. The lightning response intensity during GLE &#x00023;69 is 1.37, satisfying the requirement of the 95% confidence level (1.25) to a fairly large extent. Thus, such a positive response of lightning incidence to GLE &#x00023;69 is statistically significant.</p>
<p>Lightning variations during GLE &#x00023;70 and &#x00023;71 are analyzed in the same way. As shown in <xref ref-type="fig" rid="F2">Figure 2A</xref>, the corresponding intensity of GLE &#x00023;70 is about 9.5%. The onset time is 03:05, with an increasing phase of about 25 min. During this period, the lightning incidence also increases simultaneously until 03:35, with a more rapid growth within the first 20 min, giving the lightning response intensity of 1.11. After that, the lightning incidence generally presents a similar decreasing trend following the recovery of SEP despite some fluctuations. Overall, there is a relatively evident positive response of lightning during GLE &#x00023;70. Similarly, a Monte Carlo test is performed for this event and further demonstrates the statistical significance of this relation. As shown in <xref ref-type="fig" rid="F2">Figure 2B</xref>, the lightning response intensity of GLE &#x00023;70 satisfies the requirement of the 85% confidence level.</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p><bold>(A)</bold> Lightning response during GLE &#x00023;70 (Dec 13, 2006). <bold>(B)</bold> Results of the Monte Carlo test in 2006. The blue dotted line represents the value at 85% confidence level.</p></caption>
<graphic xlink:href="fphy-08-00348-g0002.tif"/>
</fig>
<p>As shown in <xref ref-type="fig" rid="F3">Figure 3A</xref>, the onset time of GLE &#x00023;71 is 01:55. During this event, the relative change in SEP reaches its maximum of 3.9% within 25 min after the GLE onset. Global lightning incidence keeps increasing until 02:35, giving the response intensity of 1.03. It then presents minor fluctuations during GLE recovery. Similarly, the positive response of lightning is further verified by the Monte Carlo test. As shown in <xref ref-type="fig" rid="F3">Figure 3B</xref>, the lightning response intensity of GLE &#x00023;71 also satisfies the requirement of the 80% confidence level.</p>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p><bold>(A)</bold> Lightning response during GLE &#x00023;71 (May 17, 2012). <bold>(B)</bold> Results of the Monte Carlo test in 2012. The blue dotted line represents the value at 80% confidence level.</p></caption>
<graphic xlink:href="fphy-08-00348-g0003.tif"/>
</fig>
<p>According to the case analysis for the three GLEs, global lightning incidence presents a positive response to SEP rapid enhancement during a GLE event. Meanwhile, as the GLE intensity varies from 40.0% (&#x00023;69) to 9.5% (&#x00023;70), and then to 3.9% (&#x00023;71), the corresponding lightning response intensity decreases from 1.37 to 1.11, and then to 1.03, so does the satisfied confidence level of each Monte Carlo test (95, 85, and 80%). This indicates that the lightning response intensity seems to be in direct proportion to the GLE intensity. A plausible physical mechanism is speculated as follows. The atmospheric ionization rate during GLEs increases through the avalanche of particles with high energies [e.g., [<xref ref-type="bibr" rid="B4">4</xref>, <xref ref-type="bibr" rid="B5">5</xref>]], which further intensifies the air conductivity and the atmospheric electric field [e.g., [<xref ref-type="bibr" rid="B10">10</xref>, <xref ref-type="bibr" rid="B12">12</xref>, <xref ref-type="bibr" rid="B14">14</xref>, <xref ref-type="bibr" rid="B15">15</xref>]]. Such variations in the atmosphere might eventually apply to the global electric circuit [e.g., [<xref ref-type="bibr" rid="B18">18</xref>&#x02013;<xref ref-type="bibr" rid="B20">20</xref>]]. As an essential part of the electric circuit, lightning strokes are thus more likely to be triggered during GLEs.</p>
</sec>
<sec>
<title>Latitudinal Effect</title>
<p>To investigate the latitudinal effect, we study the lightning response intensities in different latitudinal zones. For the three GLE events, most of the lightning strokes (94.0&#x02013;99.5%) locate between 40&#x000B0;S and 40&#x000B0;N. To make the comparison of sufficient statistical significance, we roughly divide the global lightning activity into two latitudinal zones. One is at low latitudes (20&#x000B0;S-20&#x000B0;N), the other is at the middle to high latitudes (&#x0003E;20&#x000B0;S or &#x0003E;20&#x000B0;N).</p>
<p><xref ref-type="table" rid="T2">Table 2</xref> gives the results of lightning response intensity in different latitudinal zones. For GLE &#x00023;69, the lightning response intensity is 1.14 at low latitudes and increases to 1.86 at higher latitudes. For GLE &#x00023;70 and GLE &#x00023;71, the values are 1.05 vs. 1.32 and 1.02 vs. 1.05, respectively. The lightning response intensity is greater at higher latitudes.</p>
<table-wrap position="float" id="T2">
<label>Table 2</label>
<caption><p>Lightning response intensities in different latitudinal zones for the three GLEs.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th valign="top" align="left"><bold>GLE</bold></th>
<th valign="top" align="left" style="border-bottom: thin solid #000000;"><bold>No</bold>.</th>
<th valign="top" align="center" style="border-bottom: thin solid #000000;"><bold>69</bold></th>
<th valign="top" align="center" style="border-bottom: thin solid #000000;"><bold>70</bold></th>
<th valign="top" align="center" style="border-bottom: thin solid #000000;"><bold>71</bold></th>
</tr>
<tr>
<th/>
<th valign="top" align="left" style="border-bottom: thin solid #000000;"><bold>Date</bold></th>
<th valign="top" align="center" style="border-bottom: thin solid #000000;"><bold>Jan 20, 2005</bold></th>
<th valign="top" align="center" style="border-bottom: thin solid #000000;"><bold>Dec 13, 2006</bold></th>
<th valign="top" align="center" style="border-bottom: thin solid #000000;"><bold>May 17, 2012</bold></th>
</tr>
<tr>
<th/>
<th valign="top" align="left"><bold>Intensity</bold></th>
<th valign="top" align="center"><bold>40.0%</bold></th>
<th valign="top" align="center"><bold>9.5%</bold></th>
<th valign="top" align="center"><bold>3.9%</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Lightning response<break/> intensity</td>
<td valign="top" align="left">Global</td>
<td valign="top" align="center">1.37</td>
<td valign="top" align="center">1.11</td>
<td valign="top" align="center">1.03</td>
</tr>
<tr>
<td/>
<td valign="top" align="left">20&#x000B0;S-20&#x000B0;N</td>
<td valign="top" align="center">1.14</td>
<td valign="top" align="center">1.05</td>
<td valign="top" align="center">1.02</td>
</tr>
<tr>
<td/>
<td valign="top" align="left">&#x0003E;20&#x000B0;N or &#x0003E;20&#x000B0;S</td>
<td valign="top" align="center">1.86</td>
<td valign="top" align="center">1.32</td>
<td valign="top" align="center">1.05</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>SEP flux during GLEs is subject to the cutoff rigidity varying in phase with the Earth&#x00027;s geomagnetic latitude. The larger SEP flux at higher latitudes is more likely to affect the lightning activity accordingly. As shown in <xref ref-type="fig" rid="F4">Figure 4</xref>, the GLE intensity for the three GLEs indeed generally increases with the geomagnetic latitude of NM stations. Thus, it is reasonable to find lightning responses more intensively at higher latitudes during GLEs. As for the low latitudes, we speculate that GLEs might indirectly affect the lightning activity through the global electric circuit, which still needs to be verified in the future.</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p>Relation between geomagnetic latitude and GLE intensity for <bold>(A)</bold> GLE &#x00023;69, <bold>(B)</bold> GLE &#x00023;70 and <bold>(C)</bold> GLE &#x00023;71. The black points represent all the NM stations which observed the corresponding GLE event. The black dashed line represents the linear fitting relation between geomagnetic latitude and GLE intensity for each event.</p></caption>
<graphic xlink:href="fphy-08-00348-g0004.tif"/>
</fig>
</sec>
</sec>
<sec sec-type="conclusions" id="s4">
<title>Conclusion</title>
<p>GLEs generally take place within several 10 min to an hour. During such a short time, the interferences of potential meteorological factors could be maximally diminished when studying the effect of CRs on global lightning activities. To our knowledge, no previous works have been done in the literature focusing on the short-term lightning response to GLEs. In this study, we investigate the lightning responses to three typical GLEs (&#x00023;69, &#x00023;70, and &#x00023;71), aiming to provide some new observational evidence for the possible link between CRs and atmospheric processes. The major results are summarized as follows:</p>
<list list-type="order">
<list-item><p>Global lightning incidence presents a significant enhancement within 20 min right after the GLE onset.</p></list-item>
<list-item><p>The relative amplitude of lightning response seems to be in direct proportion to GLE intensity, which is further verified to be statistically significant by the Monte Carlo test.</p></list-item>
<list-item><p>Lightning response to GLE seems to be more intensive at higher latitudes.</p></list-item>
</list>
<p>As the detection efficiency of the WWLLN is up to 80% of the lightning strokes above 50 kA peak current [<xref ref-type="bibr" rid="B30">30</xref>], our results are mainly valid for the lightning strokes with a peak current above 50 kA. Meanwhile, only three GLE events with the relative intensity &#x0003E;15% and the rising time &#x0003C;30 min are studied during the operational period of the WWLLN. To further verify the effect of GLEs on the global lightning activities in the timescale of several 10 min, a comprehensive statistical survey is needed in the future. Besides, regional lightning responses during GLEs are worthy of specialized investigations as well.</p>
</sec>
<sec sec-type="data-availability-statement" id="s5">
<title>Data Availability Statement</title>
<p>Publicly available datasets were analyzed in this study. This data can be found here: 1. World Wide Lightning Location Network (WWLLN): <ext-link ext-link-type="uri" xlink:href="http://wwlln.net/">http://wwlln.net/</ext-link> 2. The International GLE Database: <ext-link ext-link-type="uri" xlink:href="https://gle.oulu.fi/&#x00023;/">https://gle.oulu.fi/&#x00023;/</ext-link>.</p>
</sec>
<sec id="s6">
<title>Author Contributions</title>
<p>QW contributed to the data processing, data analysis, data interpreting, and manuscript writing. HL contributed to the study designing, data analysis, data interpreting, and manuscript revising. CW contributed to the data interpreting and manuscript revising. All authors contributed to the article and approved the submitted version.</p>
</sec>
<sec id="s7">
<title>Conflict of Interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
</body>
<back>
<ack><p>We thank the University of Oulu/Sodankyla Geophysical Observatory for providing the International GLE Database. We thank the World Wide Lightning Location Network for providing the lightning location data.</p>
</ack>
<ref-list>
<title>References</title>
<ref id="B1">
<label>1.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bazilevskaya</surname> <given-names>GA</given-names></name> <name><surname>Usoskin</surname> <given-names>IG</given-names></name> <name><surname>Fl&#x000FC;ckiger</surname> <given-names>EO</given-names></name> <name><surname>Harrison</surname> <given-names>RG</given-names></name> <name><surname>Desorgher</surname> <given-names>L</given-names></name> <name><surname>B&#x000FC;tikofer</surname> <given-names>R</given-names></name> <etal/></person-group>. <article-title>Cosmic ray induced ion production in the atmosphere</article-title>. <source>Space Sci Rev.</source> (<year>2008</year>) <volume>137</volume>:<fpage>149</fpage>&#x02013;<lpage>73</lpage>. <pub-id pub-id-type="doi">10.1007/s11214-008-9339-y</pub-id></citation></ref>
<ref id="B2">
<label>2.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Usoskin</surname> <given-names>IG</given-names></name> <name><surname>Desorgher</surname> <given-names>L</given-names></name> <name><surname>Velinov</surname> <given-names>P</given-names></name> <name><surname>Storini</surname> <given-names>M</given-names></name> <name><surname>Fl&#x000FC;ckiger</surname> <given-names>EO</given-names></name> <name><surname>B&#x000FC;tikofer</surname> <given-names>R</given-names></name> <etal/></person-group>. <article-title>Ionization of the earth&#x00027;s atmosphere by solar and galactic cosmic rays</article-title>. <source>Acta Geophys.</source> (<year>2009</year>) <volume>57</volume>:<fpage>88</fpage>&#x02013;<lpage>101</lpage>. <pub-id pub-id-type="doi">10.2478/s11600-008-0019-9</pub-id></citation></ref>
<ref id="B3">
<label>3.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gopalswamy</surname> <given-names>N</given-names></name> <name><surname>Xie</surname> <given-names>H</given-names></name> <name><surname>Yashiro</surname> <given-names>S</given-names></name> <name><surname>Akiyama</surname> <given-names>S</given-names></name> <name><surname>Makela</surname> <given-names>P</given-names></name> <name><surname>Usoskin</surname> <given-names>IG</given-names></name></person-group>. <article-title>Properties of ground level enhancement events and the associated solar eruptions during solar cycle 23</article-title>. <source>Space Sci Rev.</source> (<year>2012</year>) <volume>171</volume>:<fpage>23</fpage>&#x02013;<lpage>60</lpage>. <pub-id pub-id-type="doi">10.1007/s11214-012-9890-4</pub-id></citation></ref>
<ref id="B4">
<label>4.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mishev</surname> <given-names>AL</given-names></name> <name><surname>Velinov</surname> <given-names>PIY</given-names></name> <name><surname>Mateev</surname> <given-names>L</given-names></name> <name><surname>Tassev</surname> <given-names>Y</given-names></name></person-group>. <article-title>Ionization effect of solar protons in the Earth atmosphere&#x02014;Case study of the 20 January 2005 SEP event</article-title>. <source>Adv Space Res.</source> (<year>2011</year>) <volume>48</volume>:<fpage>1232</fpage>&#x02013;<lpage>7</lpage>. <pub-id pub-id-type="doi">10.1016/j.asr.2011.06.004</pub-id></citation></ref>
<ref id="B5">
<label>5.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mishev</surname> <given-names>AL</given-names></name> <name><surname>Velinov</surname> <given-names>PIY</given-names></name> <name><surname>Mateev</surname> <given-names>L</given-names></name> <name><surname>Tassev</surname> <given-names>Y</given-names></name></person-group>. <article-title>Ionization effect of nuclei with solar and galactic origin in the Earth atmosphere during GLE 69 on 20 January 2005</article-title>. <source>J Atmos Solar Terrestrial Phys.</source> (<year>2012</year>) <volume>89</volume>:<fpage>1</fpage>&#x02013;<lpage>7</lpage>. <pub-id pub-id-type="doi">10.1016/j.jastp.2012.06.012</pub-id></citation></ref>
<ref id="B6">
<label>6.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mishev</surname> <given-names>AL</given-names></name> <name><surname>Velinov</surname> <given-names>PIY</given-names></name></person-group>. <article-title>Time evolution of ionization effect due to cosmic rays in terrestrial atmosphere during GLE 70</article-title>. <source>J Atmos Solar Terrestrial Phys.</source> (<year>2015</year>) <volume>129</volume>:<fpage>78</fpage>&#x02013;<lpage>86</lpage>. <pub-id pub-id-type="doi">10.1016/j.jastp.2015.04.016</pub-id></citation></ref>
<ref id="B7">
<label>7.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mishev</surname> <given-names>AL</given-names></name> <name><surname>Velinov</surname> <given-names>PIY</given-names></name></person-group>. <article-title>Ion production and ionization effect in the atmosphere during the Bastille day GLE 59 due to high energy SEPs</article-title>. <source>Adv Space Res.</source> (<year>2018</year>) <volume>61</volume>:<fpage>316</fpage>&#x02013;<lpage>25</lpage>. <pub-id pub-id-type="doi">10.1016/j.asr.2017.10.023</pub-id></citation></ref>
<ref id="B8">
<label>8.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Usoskin</surname> <given-names>IG</given-names></name> <name><surname>Kovaltsov</surname> <given-names>GA</given-names></name> <name><surname>Mironova</surname> <given-names>IA</given-names></name> <name><surname>Tylka</surname> <given-names>AJ</given-names></name> <name><surname>Dietrich</surname> <given-names>WF</given-names></name></person-group>. <article-title>Ionization effect of solar particle GLE events in low and middle atmosphere</article-title>. <source>Atmos Chem Phys.</source> (<year>2011</year>) <volume>11</volume>:<fpage>1979</fpage>&#x02013;<lpage>88</lpage>. <pub-id pub-id-type="doi">10.5194/acp-11-1979-2011</pub-id></citation></ref>
<ref id="B9">
<label>9.</label>
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Velinov</surname> <given-names>PIY</given-names></name> <name><surname>Mishev</surname> <given-names>A</given-names></name> <collab>Iop Comparison of ionization effect in the atmosphere of the Earth due to GLE 65 and GLE 69</collab></person-group>. In: <source>23rd European Cosmic Ray Symposium (and 32nd Russian Cosmic Ray Conference)</source>. <publisher-loc>Moscow</publisher-loc>. (<year>2013</year>). p. <fpage>409</fpage>. <pub-id pub-id-type="doi">10.1088/1742-6596/409/1/012211</pub-id></citation></ref>
<ref id="B10">
<label>10.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Farrell</surname> <given-names>WM</given-names></name> <name><surname>Desch</surname> <given-names>MD</given-names></name></person-group>. <article-title>Solar proton events and the fair weather electric field at ground</article-title>. <source>Geophys Res Lett.</source> (<year>2002</year>) <volume>29</volume>:<fpage>37</fpage>&#x02013;<lpage>34</lpage>. <pub-id pub-id-type="doi">10.1029/2001GL013908</pub-id></citation></ref>
<ref id="B11">
<label>11.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kokorowski</surname> <given-names>M</given-names></name> <name><surname>Sample</surname> <given-names>J</given-names></name> <name><surname>Holzworth</surname> <given-names>R</given-names></name> <name><surname>Bering</surname> <given-names>E</given-names></name> <name><surname>Bale</surname> <given-names>S</given-names></name> <name><surname>Blake</surname> <given-names>J</given-names></name> <etal/></person-group>. <article-title>Rapid fluctuations of stratospheric electric field following a solar energetic particle event</article-title>. <source>Geophys Res Lett.</source> (<year>2006</year>) <volume>33</volume>:<fpage>L20105</fpage>. <pub-id pub-id-type="doi">10.1029/2006GL027718</pub-id></citation></ref>
<ref id="B12">
<label>12.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Smirnov</surname> <given-names>SE</given-names></name> <name><surname>Mikhailova</surname> <given-names>GA</given-names></name> <name><surname>Kapustina</surname> <given-names>OV</given-names></name></person-group>. <article-title>Variations in electric and meteorological parameters in the near-Earth&#x00027;s atmosphere at Kamchatka during the solar events in October 2003</article-title>. <source>Geomagnetism Aeronomy.</source> (<year>2014</year>) <volume>54</volume>:<fpage>240</fpage>&#x02013;<lpage>7</lpage>. <pub-id pub-id-type="doi">10.1134/S0016793214020182</pub-id></citation></ref>
<ref id="B13">
<label>13.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shumilov</surname> <given-names>OI</given-names></name> <name><surname>Kasatkina</surname> <given-names>EA</given-names></name> <name><surname>Frank-Kamenetsky</surname> <given-names>AV</given-names></name></person-group>. <article-title>Effects of extraordinary solar cosmic ray events on variations in the atmospheric electric field at high latitudes</article-title>. <source>Geomagnetism Aeronomy.</source> (<year>2015</year>) <volume>55</volume>:<fpage>650</fpage>&#x02013;<lpage>7</lpage>. <pub-id pub-id-type="doi">10.1134/S0016793215050151</pub-id></citation></ref>
<ref id="B14">
<label>14.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tacza</surname> <given-names>J</given-names></name> <name><surname>Raulin</surname> <given-names>JP</given-names></name> <name><surname>Mendonca</surname> <given-names>RRS</given-names></name> <name><surname>Makhmutov</surname> <given-names>VS</given-names></name> <name><surname>Marun</surname> <given-names>A</given-names></name> <name><surname>Fernandez</surname> <given-names>G</given-names></name></person-group>. <article-title>Solar effects on the atmospheric electric field during 2010&#x02013;2015 at low latitudes</article-title>. <source>J Geophys Res Atmos.</source> (<year>2018</year>) <volume>123</volume>:<fpage>11970</fpage>&#x02013;<lpage>9</lpage>. <pub-id pub-id-type="doi">10.1029/2018JD029121</pub-id></citation></ref>
<ref id="B15">
<label>15.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Golubenko</surname> <given-names>K</given-names></name> <name><surname>Rozanov</surname> <given-names>E</given-names></name> <name><surname>Mironova</surname> <given-names>I</given-names></name> <name><surname>Karagodin</surname> <given-names>A</given-names></name> <name><surname>Usoskin</surname> <given-names>I</given-names></name></person-group>. <article-title>Natural sources of ionization and their impact on atmospheric electricity</article-title>. <source>Geophys Res Lett.</source> (<year>2020</year>) <volume>47</volume>:<fpage>e2020GL088619</fpage>. <pub-id pub-id-type="doi">10.1029/2020GL088619</pub-id></citation></ref>
<ref id="B16">
<label>16.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rycroft</surname> <given-names>MJ</given-names></name> <name><surname>Israelsson</surname> <given-names>S</given-names></name> <name><surname>Price</surname> <given-names>C</given-names></name></person-group>. <article-title>The global atmospheric electric circuit, solar activity and climate change</article-title>. <source>J Atmos Solar Terrestrial Phys.</source> (<year>2000</year>) <volume>62</volume>:<fpage>1563</fpage>&#x02013;<lpage>76</lpage>. <pub-id pub-id-type="doi">10.1016/S1364-6826(00)00112-7</pub-id></citation></ref>
<ref id="B17">
<label>17.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tinsley</surname> <given-names>BA</given-names></name> <name><surname>Burns</surname> <given-names>GB</given-names></name> <name><surname>Zhou</surname> <given-names>L</given-names></name></person-group>. <article-title>The role of the global electric circuit in solar and internal forcing of clouds and climate</article-title>. <source>Adv Space Res.</source> (<year>2007</year>) <volume>40</volume>:<fpage>1126</fpage>&#x02013;<lpage>39</lpage>. <pub-id pub-id-type="doi">10.1016/j.asr.2007.01.071</pub-id></citation></ref>
<ref id="B18">
<label>18.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rycroft</surname> <given-names>MJ</given-names></name> <name><surname>Harrison</surname> <given-names>RG</given-names></name> <name><surname>Nicoll</surname> <given-names>KA</given-names></name> <name><surname>Mareev</surname> <given-names>EA</given-names></name></person-group>. <article-title>An overview of earth&#x00027;s global electric circuit and atmospheric conductivity</article-title>. <source>Space Sci Rev.</source> (<year>2008</year>) <volume>137</volume>:<fpage>83</fpage>&#x02013;<lpage>105</lpage>. <pub-id pub-id-type="doi">10.1007/s11214-008-9368-6</pub-id></citation></ref>
<ref id="B19">
<label>19.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rycroft</surname> <given-names>MJ</given-names></name> <name><surname>Nicoll</surname> <given-names>KA</given-names></name> <name><surname>Aplin</surname> <given-names>KL</given-names></name> <name><surname>Giles Harrison</surname> <given-names>R</given-names></name></person-group>. <article-title>Recent advances in global electric circuit coupling between the space environment and the troposphere</article-title>. <source>J Atmos Solar Terrestrial Phys.</source> (<year>2012</year>) <volume>90&#x02013;1</volume>:<fpage>198</fpage>&#x02013;<lpage>211</lpage>. <pub-id pub-id-type="doi">10.1016/j.jastp.2012.03.015</pub-id></citation></ref>
<ref id="B20">
<label>20.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nicoll</surname> <given-names>KA</given-names></name></person-group>. <article-title>Space weather influences on atmospheric electricity</article-title>. <source>Weather.</source> (<year>2014</year>) <volume>69</volume>:<fpage>238</fpage>&#x02013;<lpage>41</lpage>. <pub-id pub-id-type="doi">10.1002/wea.2323</pub-id></citation></ref>
<ref id="B21">
<label>21.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lethbridge</surname> <given-names>MD</given-names></name></person-group>. <article-title>Cosmic rays and thunderstorm frequency</article-title>. <source>Geophys Res Lett.</source> (<year>1981</year>) <volume>8</volume>:<fpage>521</fpage>&#x02013;<lpage>2</lpage>. <pub-id pub-id-type="doi">10.1029/GL008i005p00521</pub-id></citation></ref>
<ref id="B22">
<label>22.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gurevich</surname> <given-names>AV</given-names></name> <name><surname>Milikh</surname> <given-names>GM</given-names></name> <name><surname>Roussel-Dupre</surname> <given-names>R</given-names></name></person-group>. <article-title>Runaway electron mechanism of air breakdown and preconditioning during a thunderstorm</article-title>. <source>Phys Lett A.</source> (<year>1992</year>) <volume>165</volume>:<fpage>463</fpage>&#x02013;<lpage>8</lpage>. <pub-id pub-id-type="doi">10.1016/0375-9601(92)90348-P</pub-id></citation></ref>
<ref id="B23">
<label>23.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gurevich</surname> <given-names>AV</given-names></name> <name><surname>Zybin</surname> <given-names>KP</given-names></name> <name><surname>Roussel-Dupre</surname> <given-names>RA</given-names></name></person-group>. <article-title>Lightning initiation by simultaneous effect of runaway breakdown and cosmic ray showers</article-title>. <source>Phys Lett A.</source> (<year>1999</year>) <volume>254</volume>:<fpage>79</fpage>&#x02013;<lpage>87</lpage>. <pub-id pub-id-type="doi">10.1016/S0375-9601(99)00091-2</pub-id></citation></ref>
<ref id="B24">
<label>24.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chronis</surname> <given-names>T</given-names></name></person-group>. <article-title>Investigating possible links between incoming cosmic ray fluxes and lightning activity over the United States</article-title>. <source>J Climate.</source> (<year>2009</year>) <volume>22</volume>:<fpage>1</fpage>. <pub-id pub-id-type="doi">10.1175/2009JCLI2912.1</pub-id></citation></ref>
<ref id="B25">
<label>25.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kristj&#x000E1;nsson</surname> <given-names>J</given-names></name> <name><surname>Stjern</surname> <given-names>C</given-names></name> <name><surname>Stordal</surname> <given-names>F</given-names></name> <name><surname>Fjaeraa</surname> <given-names>A</given-names></name> <name><surname>Myhre</surname> <given-names>G</given-names></name> <name><surname>J&#x000F3;nasson</surname> <given-names>K</given-names></name> <etal/></person-group>. <article-title>Cosmic rays, cloud condensation nuclei and clouds&#x02014;A reassessment using MODIS data</article-title>. <source>Atmos Chem Phys.</source> (<year>2008</year>) <volume>8</volume>:<fpage>7373</fpage>&#x02013;<lpage>87</lpage>. <pub-id pub-id-type="doi">10.5194/acp-8-7373-2008</pub-id></citation></ref>
<ref id="B26">
<label>26.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Okike</surname> <given-names>O</given-names></name> <name><surname>Umahi</surname> <given-names>AE</given-names></name></person-group>. <article-title>Cosmic ray&#x02014;global lightning causality</article-title>. <source>J Atmos Solar Terrestrial Phys.</source> (<year>2019</year>) <volume>189</volume>:<fpage>35</fpage>&#x02013;<lpage>43</lpage>. <pub-id pub-id-type="doi">10.1016/j.jastp.2019.04.002</pub-id></citation></ref>
<ref id="B27">
<label>27.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Holzworth</surname> <given-names>RH</given-names></name> <name><surname>McCarthy</surname> <given-names>MP</given-names></name> <name><surname>Brundell</surname> <given-names>JB</given-names></name> <name><surname>Jacobson</surname> <given-names>AR</given-names></name> <name><surname>Rodger</surname> <given-names>CJ</given-names></name></person-group>. <article-title>Global distribution of superbolts</article-title>. <source>J Geophys Res Atmos.</source> (<year>2019</year>) <volume>124</volume>:<fpage>9996</fpage>&#x02013;<lpage>10005</lpage>. <pub-id pub-id-type="doi">10.1029/2019JD030975</pub-id></citation></ref>
<ref id="B28">
<label>28.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dowden</surname> <given-names>RL</given-names></name> <name><surname>Brundell</surname> <given-names>JB</given-names></name> <name><surname>Rodger</surname> <given-names>CJ</given-names></name></person-group>. <article-title>VLF lightning location by time of group arrival (TOGA) at multiple sites</article-title>. <source>J Atmos Solar Terrestrial Phys.</source> (<year>2002</year>) <volume>64</volume>:<fpage>817</fpage>&#x02013;<lpage>30</lpage>. <pub-id pub-id-type="doi">10.1016/S1364-6826(02)00085-8</pub-id></citation></ref>
<ref id="B29">
<label>29.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mialdea-Flor</surname> <given-names>I</given-names></name> <name><surname>Segura-Garcia</surname> <given-names>J</given-names></name> <name><surname>Felici-Castell</surname> <given-names>S</given-names></name> <name><surname>Garcia-Pineda</surname> <given-names>M</given-names></name> <name><surname>Alcaraz-Calero</surname> <given-names>MJ</given-names></name> <name><surname>Navarro-Camba</surname> <given-names>E</given-names></name></person-group>. <article-title>Development of a low-cost iot system for lightning strike detection and location</article-title>. <source>Electronics.</source> (<year>2019</year>) <volume>8</volume>:<fpage>1512</fpage>. <pub-id pub-id-type="doi">10.3390/electronics8121512</pub-id></citation></ref>
<ref id="B30">
<label>30.</label>
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rodger</surname> <given-names>C</given-names></name> <name><surname>Werner</surname> <given-names>S</given-names></name> <name><surname>Brundell</surname> <given-names>J</given-names></name> <name><surname>Lay</surname> <given-names>E</given-names></name> <name><surname>Thomson</surname> <given-names>N</given-names></name> <name><surname>Holzworth</surname> <given-names>R</given-names></name> <etal/></person-group>. <article-title>Detection efficiency of the VLF World-Wide Lightning Location Network (WWLLN): initial case study</article-title>. <source>Annal Geophys.</source> (<year>2006</year>) <volume>24</volume>:<fpage>3197</fpage>&#x02013;<lpage>214</lpage>. <pub-id pub-id-type="doi">10.5194/angeo-24-3197-2006</pub-id></citation></ref>
<ref id="B31">
<label>31.</label>
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Usoskin</surname> <given-names>I</given-names></name> <name><surname>Ibragimov</surname> <given-names>A</given-names></name> <name><surname>Shea</surname> <given-names>M</given-names></name> <name><surname>Smart</surname> <given-names>D</given-names></name></person-group>. <article-title>Database of Ground Level Enhancements (GLE) of high energy solar proton events</article-title>. In: <source>34th International Cosmic Ray Conference (ICRC2015)</source>. <publisher-loc>Hague</publisher-loc>. (<year>2016</year>). p. <fpage>34</fpage>. <pub-id pub-id-type="doi">10.22323/1.236.0054</pub-id></citation></ref>
</ref-list>
<fn-group>
<fn id="fn0001"><p><sup>1</sup>WWLLN: <ext-link ext-link-type="uri" xlink:href="http://wwlln.net/">http://wwlln.net/</ext-link></p></fn>
<fn id="fn0002"><p><sup>2</sup>International GLE database: <ext-link ext-link-type="uri" xlink:href="https://gle.oulu.fi/&#x00023;/">https://gle.oulu.fi/&#x00023;/</ext-link></p></fn>
</fn-group>
<fn-group>
<fn fn-type="financial-disclosure"><p><bold>Funding.</bold> This work was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (CAS, China), grant no. XDA17010301; the National Natural Science Foundation of China, grant nos. 41874203, 41574169, 41574159, and 41731070; the Young Elite Scientists Sponsorship Program of the China Association for Science and Technology, grant no. 2016QNRC001; the Strategic Pioneer Program on Space Science of CAS, China, grant no. XDA15052500; and in part by the Specialized Research Fund for State Key Laboratories of China. HL is also supported by the Youth Innovation Promotion Association of CAS, China and the research fund for key development directions of the National Space Science Center of CAS, China. CW is also supported by the Key Research Program of Frontier Sciences of CAS, China, grant no. QYZDJ-SSW-JSC028.</p></fn>
</fn-group>
</back>
</article>