A Review of Atmospheric Electricity Research in China from 2011 to 2018

In the study of atmospheric electricity, various electrical processes taking place in the atmosphere are considered, usually within a framework of a global electric circuit, in which thunderstorm regions and fair-weather regions, together with the earth and the lower ionosphere, are combined into a distributed electric circuit. A central issue concerning research on atmospheric electricity in China is lightning and its effects in a wide range of aspects. (Qie, 2012) reviewed the research progress in atmospheric electricity in China during the period 2006-10. In the present paper, the period 2011-18 is reviewed.

Most achievements in the study of atmospheric electricity in the last few years have come from the coordinated research efforts of the "Program on Dynamic, Microphysical and Electrical Processes in Severe Thunderstorms and Lightning Hazards" (Strom973), which was funded as a National Key Basic (973) Research Program by the Ministry of Science and Technology, China. The Storm973 program carried out six major research topics: (1) Coordinated Observation of Lightning and Thunderstorms based on New Detection Techniques; (2) Evolution of Dynamical Processes and Structures of Strong Thunderstorms; (3) Microphysical Characteristics and their Impacts on Electrical Charge Distribution Inside Thunderstorms; (4) Charge Structure and Lightning Development in Severe Thunderstorms; (5) Physics of Lightning Discharge and Hazard Causes; and (6) Lightning Data Assimilation and Forecasting of Thunderstorms and Lightning.

This review paper covers the following five aspects: (1) lightning detection and location techniques; (2) discharge processes and parameters associated with rocket-triggered lightning; (3) physical processes in natural lightning and attachment to ground; (4) lightning activities and charge structure in different thunderstorms; and (5) effects of thunderstorms on the upper atmosphere. In addition, some outstanding questions for future research are outlined.

New lightning observational capabilities have been critical in allowing new understanding on the nature of lightning. Besides steady improvements in electronics and signal processing, a major advance for lightning detection has been the advent of the GPS constellation, which provides extremely accurate timing worldwide. These new capabilities have enabled several groups in China to develop the 3D positioning systems for lightning radiation pulses based on the time-of-arrival technique or interferometry techniques.

In addition to operational cloud-to-ground (CG) lightning location networks covering all of China, several 3D lightning location research systems have been developed and put forward in lightning detection and research. Figure 1 shows the distribution of major research lightning networks in China's mainland. (Zhang et al., 2010) developed the first 3D very-high-frequency (VHF) lightning radiation source mapping technique in China, which was first installed in Zhanhua, Shandong Province. It was then moved to Qinghai Province to investigate the lightning characteristics and charge structure associated with plateau thunderstorms. The system is operated at 270 MHz with a 3-dB bandwidth of 6 MHz and processes peak events in a consecutive time window of 50 μs (Li et al., 2017b). (Zhang et al., 2015a) estimated the location error for this system using a balloon-borne VHF transmitter. The system showed a horizontal error of 12-48 m and a vertical uncertainty of 20-78 m for radiation sources over the network. A broadband electric field (E-field) location system working at a bandwidth from 1.5 kHz to 10 MHz is installed synchronously (Li et al., 2013).

Figure 1.  Distribution of major research lightning networks (red stars) and sites of rocket-triggering lightning experiments (blue triangles) in China's mainland.

A multiband 3D lightning location network has been installed in Beijing, referred to as the Beijing Lighting Network (BLNet). BLNet is a regional lightning location network with one data center and sixteen substations. Each station is equipped with fast antenna, used for real-time location, while two systems from slow antenna, magnetic antenna and VHF antenna are selected to be equipped at each station (Wang et al., 2015c, 2016d), covering a bandwidth from very low frequency (VLF) to VHF. According to (Srivastava et al., 2017), the average detection efficiency (DE) of BLNet for total flashes is 93.2%, and 73.9% for CG flashes. The location error in the horizontal direction ranges from 52.9 m to 250.0 m, based on high-tower lightning flashes.

A low frequency (LF) E-field Detection Array (LFEDA) is installed in Conghua, Guangdong, aimed at locating the 3D positions of lightning pulse discharge events (Shi et al., 2017). After checking the LFEDA locations for triggered lightning, it was found that LFEDA has an average horizontal location error of 102 m for return strokes (RSs), with the DE for flashes and RSs being 100% and 95%, respectively. (Fan et al., 2018c) developed a new 3D location method for the LFEDA, which significantly increased the capability of the system. For a fatal "bolt from the blue", the mean location error for five strokes was 27 m.

Another lightning location network was set up in Chongqing. The lightning sensors in the network detect VLF/ LF and VHF sources radiated by lightning (Liu et al., 2018b). The network is composed of 14 lightning sensors. The minimum and maximum distance between two sensors is about 10 km and 100 km. Each sensor records the VLF/LF and VHF signals continuously at sample rates of 5 MHz and 20 MHz, respectively.

In addition, there are several lightning location networks operating in China and it is important to evaluate the detection performance of these systems. (Chen et al., 2012) investigated the performance of three lightning location systems (LLSs) in Guangdong Province based on tower lightning and rocket-triggered lightning. The location error, flash DE, stroke DE, and peak current absolute percentage errors are around 710 m, 94%, 60% and 16.3%, respectively, for the Guangdong Power Grid LLS. For the Guangdong-Hong Kong-Macao LLS (GHMLLS), the flash DE, stroke DE, location error are 74%, 96% and 532 m, respectively. However, the GHMLLS-inferred peak currents were found to be about 37% lower than the actual value during 2012-15. The corresponding values for the Earth Networks Total Lightning Network are 77%, 76%, 685 m and 39% (Zhang et al., 2016b).

The performance of the World Wide Lightning Location Network (WWLLN) was evaluated over the Tibetan Plateau and the ocean. (Fan et al., 2018a) found that the location accuracy of the WWLLN was around 10 km and the DE was 9.37% for CG and 2.58% for total flashes over the Tibetan Plateau. The DE of the WWLLN increased year by year from 4.3% to 19.1% over the Northwest Pacific from 2005 to 2016 (Pan et al., 2013; Wang et al., 2018b; Zhang et al., 2018b). It was roughly 16.8% in Beijing in 2015 and 2016 (Srivastava et al., 2017).

Interferometry location technology, which measures the phase difference or time difference of the lightning signals from different antennas with very short baselines, has developed quickly in China. This system is usually composed of a four orthogonal broadband antenna array with baselines from around 10 m to several tens of meters (Dong et al., 2002). (Cao et al., 2012) developed the short-baseline 2D LLS with a bandwidth from 125 MHz to 200 MHz to locate lightning discharge processes. (Sun et al., 2013) proposed an algorithm using a time delay estimation based on generalized cross correlation together with wavelet transform to reduce the noise interfaces. The location accuracy for the discharge channel was highly improved. The system has been used to map rocket-triggered lightning and natural lightning discharges (Sun et al., 2014a, b, 2015). Figure 2 shows a VHF mapping example for rocket-triggered lightning with comparison to an optical image. It can be seen that the VHF-mapped discharge channel is very well matched to the high-speed video image in the view of the camera, with a bias of less than one Charge Coupled Device (CCD) element view, which is less than 10 m at a distance of 980 m (Sun et al., 2014b). (Liu et al., 2018b) conducted an observation of two broadband interferometers that were about 8.15 km apart and obtained 3D lightning location discharge results. The results indicated that the system effectively mapped the lightning radiation sources and were in good agreement with optical images.

Figure 2.  A triggered lightning discharge channel mapped by a VHF interferometer with comparison to an optical image; modified from (Sun et al., 2014b).

Rocket-triggered lightning experiments using the rocket-wire technique were conducted continuously in Zhanhua, Shandong (SHATLE) and Conghua, Guangdong (GCOELD). The locations of the two sites are shown on the map in Fig. 1. Some new insights into the lightning process and the corresponding electromagnetic (EM) effects have been obtained, having profited from the ability to detect the discharge current and the simultaneous close EM field and channel luminosity evolution. Figure 3 shows site pictures in Zhanhua and Conghua of the triggered lightning experiments using the rocket trailing wire technique. The triggering facilities and main observation devices are also marked. A schematic view of both the triggering site and main observation site in the Zhanhua experiment is also shown in Fig. 3. The primary installation of current and EM field measuring systems are similar at both sites. A schematic view of the Conghua experiment is available in the literature (e.g., Zhang et al., 2016b; Zheng et al., 2017).

Figure 3.  Installation of rocket-triggering lightning experiment: upper left panel——in SHATLE; upper right panel——in GCOELD; lower panel——schematic view in the Zhanhua experiment.

Rocket-triggered lightning is usually negative in polarity and initiated with an upward positive leader (UPL). The propagating features of UPLs in the initial stage (IS) of triggered lightning were investigated both for SHATLE (Qie et al., 2011, Qie et al., 2012, Qie et al., 2014a, Qie et al., 2017; Wang et al., 2012a; Jiang et al., 2013a) and GCOELD (Zhang et al., 2014c, Zhang et al., 2016b; Zheng et al., 2017). The average 2D speed of a UPL was found to be roughly 1.0× 10⁵ m s^-1 and showed a tendency to accelerate with height (Jiang et al., 2013a). The peak current for IS pulses ranged from tens of amperes to about 150 A. After a UPL develops to several hundreds of meters high, some impulses can be observed that are superimposed on the initial continuous current (ICC) and the associated E-field (Lu et al., 2014; Fan et al., 2017).

(Zhang et al., 2017d) found with a VHF interferometer that precursor pulses appear first during a period of hundreds of milliseconds before the UPL's sustained development in the IS of a triggered lightning discharge. The individual precursor pulses were initiated by weak upward positive breakdown, meters in scale, followed by fast, energetic downward negative breakdown, tens of meters in scale. Their average speed was about 5×10⁶ m s^-1 and 3× 10⁷ m s^-1, respectively.

Lu et al. (2014, 2018) and Fan et al. (2017, 2018d) introduced sensitive LF magnetic sensors into SHATLE in 2013. The UPL propagation was found to be associated with a burst of EM impulses when it entered the in-cloud negative charge center, possibly suggesting a stepped pattern of positive leader propagation even in the cloud (Lu et al., 2014; Zheng et al., 2018). The UPL propagation is active in LF radiation in this period, while it is relatively quiet in the VHF band, implying a possible fundamental change in the progression mechanism of the positive leader. (Lu et al., 2016a) classified these magnetic field signals into impulsive and ripple pulses according to their waveform features. At the stage of the impulsive type, the leader channel had a negligible length, so the associated current pulses mainly propagated downward along the high conductivity steel wire into the ground. Based on the above fact, (Fan et al., 2018d) simulated the pulse of the impulsive type by using a transmission-line model, suggesting that the pulses are generated by leader current pulses propagating downward along the steel wire. The waveform of impulsive current pulses is changed to ripple pulses after traveling through the high impedance of the prolonging leader channel.

Two cases of rare positive lightning flashes were triggered in SHATLE with an upward negative leader (UNL) initiation, but just followed by the ICC stage and without an RS. Figure 4 shows composite UNL images from a high-speed video camera, together with the UPL images in rocket-triggered lightning. A branched channel structure is apparent for the UNL, while there are no obvious branches for the UPL. (Pu et al., 2017) found that the UNL branched shortly after its initiation, and consequently the stepping-caused current pulses became multi-peaked and wider in their waveforms. The negative charge of an individual leader step was nearly an order of magnitude larger than that of the positive counterpart of the triggered lightning. The main negative leader tip propagated at an average speed of 2.1× 10⁵ m s^-1. Two cases of positive lightning flashes without an RS were also triggered in GCOELD (Zheng et al., 2017). The geometric mean (GM) value of the 2D speed of the UNL was 1.79× 10⁵ m s^-1, and the ISs showed much larger currents, charge transfers and action integrals relative to their negative counterpart in GCOELD.

Figure 4.  Composite images of a UNL (left panel) and UPL in rocket-triggered lightning (right panel).

Dart leaders usually propagate fast, and a high-speed video camera cannot resolve them well. The VHF/UHF lightning interferometer was upgraded to operate at a continuous recording mode or sequential mode, and is capable of capturing weak radiation signals from lightning, like UPLs, and the discharge channels for downward negative dart leaders can be mapped clearly too (Sun et al., 2014b, 2016). For a triggered negative lightning flash with 16 leader and RS sequences, the dart leader speed was mostly on the order of 10⁶-10⁷ m s^-1. For dart-stepped leaders, it was of 10⁵ m s^-1.

Bidirectional leader propagation was detected for the first time in triggered lightning discharges by (Qie et al., 2017). The bileader initiated almost immediately below a decaying downward leader, which stopped midway before reaching the ground. The average speed of the UPL end was 1.3× 10⁶ m s^-1 and 2.2× 10⁶ m s^-1 for the two cases, respectively——roughly two times larger than the negative downward leader end. The bidirectional leader might be considered as a polarity-reversed recoil leader with its positive leader end retrogressing along the existing negative channel, and the polarity of the recoil leader here was completely reversed to the traditional one. (Zheng et al., 2012) reported an abnormal triggered lightning flash containing two upward propagation processes, with the second upward leader being triggered by a downward aborted leader. They propagated along the same channel but with different polarities and opposite directions.

Figure 5 shows an example of the whole current waveform and partial enlarged waveforms of a classical triggered lightning event in Conghua, Guangdong. The flash lasted about 800 ms. The large pulses were produced by 11 RSs, with the peak current from 4 kA to 21 kA. Based on the current and simultaneous close EM field, the discharge processes and corresponding waveform parameters of current pulses were analyzed statistically (Jiang et al., 2011; Wang et al., 2012a; Zhao et al., 2011; Qie et al., 2014a; Zhang et al., 2014c; Zheng et al., 2017). Four categories of current pulses have been found to be associated with triggered flashes: RSs, ICC pulses (ICCPs), M-components superimposed on the continuing current (CC) following the RS, and mixed return stroke- M component (RM) pulses featuring both RS and M-event characteristics. Figure 6 shows the current and corresponding E-field at 60 m produced by these impulsive processes. Differing behavior of the current and E-field can be found between these pulses.

Figure 5.  Example of the current waveform of a classical triggered lightning flash in Conghua, Guangdong: (a) whole current waveform; (b) IS current; (c, d) enlarged current waveforms of the third and seventh RS and the following CCs and M-component.

Figure 6.  Waveform example of current and close electric field at 60 m from the discharge channel for the RS, M-event, ICCP and RM event documented in Zhanhua, Shandong.

Table 1 compares the RS parameters of rocket-triggered lightning from the two sites. The GM values of the peak current, half peak width, risetime from 10% to 90% peak, and charge transfer for the RSs were found to be 12.1 kA, 14.8 μs, 1.0 μs and 0.86 C, respectively, while they were 0.28 kA, 242 μs, 251 μs and 0.10 C for the M-components, and 0.09 kA, 712 μs, 437 μs and 0.10 C for the ICC pulses in SHATLE (Qie et al., 2014a). Some of the M-components might exhibit a peak current with several kilo-amperes (Jiang et al., 2011), comparable to some of the weak RSs. The charge transferred to the ground by an individual lightning discharge ranged from 6.3 C to 68.1 C.

(Zheng et al., 2017) presented the current characteristics of the IS and RS of triggered negative lightning in GCOELD. The IS had GMs of 347.9 ms for duration, a 132.5 A average current, 45.1 C charge transfer and 10.0× 10³ A² s action integral, with larger values than those reported elsewhere. The RS featured a greater peak current, charge transfer and action integral within 1 ms, but a shorter 10% to 90% rise time than elsewhere. A triggered lightning flash containing RSs tended to last longer and neutralize fewer charges during its IS than that without RSs. The RSs seem stronger in Conghua than in Zhanhua in view of the peak current and charge transfer.

(Zheng et al., 2013) investigated the currents of IS RSs (ISRSs) in two altitude-triggered lightning flashes with gaps of 17 m and 50 m from the wire lower end to the ground, and the peak current was 10.19 kA and 9.03 kA, respectively. They concluded that the current waves of ISRSs were similar to those of typical RSs in classical triggered flashes in terms of multiple parameters, except for the smaller charge transfer, half-peak width, and action integral of ISRSs.

(Zhang et al., 2016d) found that M-components with small amplitude (<0.5 kA) were necessary for a long CC, and named the phenomenon the "restricted zone". Some M-events with current exceeding 1 kA might result from fast positive streamers with a speed of about 10⁷ m s^-1, or by a dart leader from other branches of the lightning channel while the channel of the CC is still existing (Zhang et al., 2018c). (Zhou et al., 2013) calculated the correlation between current and luminosity of the lightning channel during the period of the ICC process and CC process following RSs, and the same variation trend between the channel CC current and the luminosity was found.

(Jiang et al., 2013b) modified the M-components model based on Rakov's "two-wave" theory and confirmed that the M-component evolution through the lightning CC channel involved a downward process transferring negative charge to the lower channel from the upper part, and an upward process that neutralized the charge deposited by the downward process. (Wang et al., 2012b) simulated the subsidiary peak in the current waveform of triggered lightning strokes, and four possibilities of subsidiary peaks were suggested: channel branching, the corona current, and flashover along the triggering wire in a previous unsuccessful launch, or reflection of the current.

Observation of lightning flashes has been conducted in several regions with different climatic or orographic features in China, including the Tibetan Plateau (Cao et al., 2011b; Wang et al., 2011b), Daxing'Anling (Wang et al., 2011a; Lu et al., 2013; Wu et al., 2013a), Beijing (Wu et al., 2016c, d; Li et al., 2017a), Qinghai (Li et al., 2012; Qu et al., 2012a, b; Zhang et al., 2015a), Guangdong (Lan et al., 2011; Zhang et al., 2016d) and the Yangtze-Huaihe river basin (Zhu et al., 2014, Zhu et al., 2015). Some of the main results are outlined in the following subsections.

The preliminary breakdown process (PBP) and initial radiation pulses in CG lightning have been studied in terms of initiation height and pulse characteristics (Cao et al., 2011a; Zhang et al., 2013b, Zhang et al., 2015c; Wang et al., 2014a; Wu et al., 2016a). (Wu et al., 2016a) divided the PBP into an initial PBP and a subsequent PBP. They found that when multiple pulse clusters were included in the initial PBPs, the initial streamer exhibited a discontinuous channel with a stepped-manner propagation traveling downward or upward. Each step corresponded to a pulse cluster. (Wang et al., 2016d), using the BLNet 3D location results, found a clear branched 3D structure of the PBP in a negative CG (-CG) flash. The PBP started at an altitude of about 6 km and developed downward to about 3 km with a speed of 5.9× 10⁵ m s^-1. An intercloud (IC) discharge initiated from about 7.2 km and propagated upward to 10 km with a velocity of 4.8× 10⁵ m s^-1.

Three types of initial breakdown pulse trains in positive CG (+CG) lightning were recognized by (Zhang et al., 2013b). The category of the pulse trains is the same in Beijing and Guangzhou, but with different percentages of each category. (Zhang et al., 2015c) analyzed the characteristics of bipolar pulse trains involved in E-field waveforms of lightning discharges. They indicated that the ratio of the largest peak of the PB pulse train to that of the first RS was obviously larger in Beijing, perhaps because of the different scales of the lower positive charge region (LPCR) in thunderstorms in Beijing and Guangzhou.

(Jiang et al., 2017) found that branching in downward negative leaders resulted from the multiple connections of clustered space leaders with the same parent channel, and connections occurred successively or almost simultaneously as some of the space leaders resurrected after their termination. The GM value of the step length was 4.4 m. The distance from the space leader center to the main channel tip was about 3.6 m. More than 50% of steps occurred within ±30^° from the direction of the advancing leader. (Qi et al., 2016) detected a total of 23 space stems/leaders from a natural downward negative lightning discharge about 350 m away. The mean 2D length was 5 m, varying from 1 to 13 m. The mean distance between the bright segments and the main leader channel tip was about 4 m, ranging from 1 to 8 m.

(Lan et al., 2011) reported direct measurements of the broadband EM field radiated from chaotic pulse trains (CPTs) associated with the leader-subsequent strokes in -CG lightning, and found that the corresponding channel extension speed was 2× 10⁷ m s^-1. (Zhang et al., 2016c) found that 44% of subsequent strokes were preceded by CPTs and almost all CPTs were characterized by strong luminosity. (Qiu et al., 2015) found that the chaotic leader propagated at a speed of (1.0-2.9) ×10⁷ m s^-1 and the duration of the CPT was 300-700 μs. The chaotic fluctuations were believed to be related to the continuous VHF radiation.

High spatial and temporal resolution lightning data have enabled many details involved in upward lightning from high structures in China to be revealed (Lu et al., 2012, Lu et al., 2013, Lu et al., 2015, Lu et al., 2016b; Gao et al., 2014; Jiang et al., 2014a, 2014b, 2014c; Qi et al., 2016; Zhang et al., 2017a; Wang et al., 2016e; 2018c).

(Lu et al., 2013) found a negative lightning discharge exhibiting unexpected attachment behaviour with the tip of the downward leader connecting to the lateral surface of an upward connecting leader (UCL). The connecting behaviour of the upward leaders and the downward leader involved in the attachment process preceding the first RS could be divided into three types: a downward leader tip to the UCL tip, a downward leader tip to the UCL lateral surface, and a hybrid type that is a combination of the above two. The three types accounted for 42%, 50% and 8% of all the events, respectively. (Lu et al., 2012) found that the unconnecting upward leader (UUL), initiated by the downward leader before a new striking point, propagated 1 km from the inception point. The 3D length for the six UCLs ranged from 180 m to 818 m, while the 3D speed ranged from 0.8× 10⁵ m s^-1 to 14.3× 10⁵ m s^-1 (Gao et al., 2014). The average 3D speed can be 1.3 times larger than the corresponding 2D speed. The 3D speed for both the UCL and UUL exhibited a trend of increasing after their inception, while the downward leader did not show a clear variation trend (Lu et al., 2015).

(Jiang et al., 2014b) and (Yuan et al., 2017) conducted observations of lightning flashes striking at a 325 m high tower. They found that 95% of the 20 cases were initiated from the tower, and most of them were induced by nearby +CG flashes. The approach of in-cloud horizontal negative leader during the +CG is a vital condition for the inception of the upward leader from the tower. The other-triggered lightning-was more likely to occur in the thunderstorm dissipation stage, with a relatively lower cloud top and weaker radar echo, whereas self-initiating lightning tended to occur in the mature stage of the thundercloud and the stratiform region of the thundercloud trailed by the convective line was above the tower.

(Wang et al., 2016e) presented a detailed picture of the stepping process in a UPL initiated from the 325 m high tower. It was found that the stepping process of the UPL occurred as an abrupt extension of the leader channel accompanied by a burst of corona streamer from the fresh leader top. The burst of corona streamer at the end of the stepping process was significantly brighter than the enhancing streamer zone when the stepping process started. The mean value of the step length was 4.9 m, and the interstep interval was 62 μs. The average speed of the leader propagation was 8.1× 10⁴ m s^-1, while the mean speed of channel extension for an individual step should be higher than 7.3× 10⁵ m s^-1, which is almost an order of magnitude larger than the average leader propagation speed. (Jiang et al., 2014a) documented the bidirectional leader propagation through a pre-existing channel created by a UPL from a high structure, similar to the bidirectional leader in the rocket-triggered lightning mentioned above (Qie et al., 2017).

The existence of high structures affects the lightning activity around them. (Zhang et al., 2017a) found that the flash or stroke density showed a considerable increase around the Canton Tower within a radius of 1 km and a clear radial decrease from 1 km to 4 km. The LLS-inferred stroke peak current occurred around the Canton Tower showed an obvious increase within a radius of 1 km.

Substantial progress has been made in the past few years in understanding the processes of lightning attachment. The attachment process and associated leader behavior for natural lightning were captured by high-speed camera and VHF interferometer measurements. For stepped leader-first RSs, (Jiang et al., 2015) found that the UCLs were more likely to be initiated by those brighter downward branches with lower tips, and there was a greater possibility that they might achieve attachment. The UCLs in two leader-subsequent stroke processes exhibited relatively long lengths, at 340 m and 105 m for the two cases, respectively. Evidence concerning the mechanism for multiple groundings has been provided by sensitive broadband VHF interferometers. It was revealed that multiple-ground terminations of a negative lightning event resulted from either the multiple termination strokes or forked stroke with the same main leader channel, or from the different branches of the PBP inside the cloud (Sun et al., 2016).

Although just around 10% of all natural CG lightning flashes are of positive polarity, they tend to be more intense and possess greater potential damage than negative flashes. (Qie et al., 2013) investigated +CG lightning characteristics in the Daxing'Anling forest region. They documented 185 +CG flashes containing 196 RSs and showed that 71.9% of all +CG discharges contained a CC process, and 94.6% of all the +CG flashes were single-stroke flashes. About 15.1% of the +CG flashes were a kind of byproduct of IC flashes.

(Kong et al., 2015) found that the speed of five positive downward leaders increased as they approached the ground, with an average 2D speed ranging from 0.3× 10⁵ m s^-1 to 2.0× 10⁵ m s^-1. Approximately 67.4% of the 89 observed +CG flashes accompanied a prior long-lasting intense IC discharge ranging from about 100 ms to 973 ms, suggesting that a previous in-cloud discharge is conducive to +CG flashes.

(Chen et al., 2015) observed a downward bipolar lightning discharge hit at a 90 m tall object. The six negative strokes developed along the same discharge channel with the first positive RS. The velocity of the RS ranged from 1.0 to 1.3× 10⁸ m s^-1. The leader before the positive RS developed downward with a 2D velocity of 2.5× 10⁶ m s^-1 without obvious branches. (Tian et al., 2016) found that the onset of the bipolar flash was followed by extension of several positive leaders near the cloud bottom, and one of them developed downward culminating with a positive stroke with a long CC. Another positive leader propagated horizontally toward a possible far negative charge center and initiated several recoil leaders that sporadically retrograded in the discharge channel. Finally, three recoil leaders developed along the positive stroke path and induced negative strokes respectively, which ultimately resulted in the polarity reversal of the following RSs.

The lightning-radiated EM signal propagates along the ground surface, and its far EM field is an important parameter for accurately locating lightning. The effects of lightning EM field propagation along a ground surface path with different characteristics (Zhang et al., 2012a, c), and a rough surface (Zhang et al., 2012b; Li et al., 2014), have been widely studied. (Zhang et al., 2015b) studied the propagation effect of frequency-dependent soil on the far vertical E-fields radiated by subsequent lightning strikes to tall objects. They found that the field propagation attenuation along frequency-dependent soil is obviously less than the case where the parameters are assumed to be constant. Hou et al. (2018a, b) presented a new approximating method for lightning-radiated extremely low frequency (ELF) and VLF ground wave propagation over intermediate ranges within 1500 km from lightning strike points, and it was found that the curvature of the earth had a much greater effect on the field propagation at intermediate ranges than the earth's finite conductivity, and the lightning-radiated ELF/VLF E-field peak value (V m^-1) at intermediate ranges yielded a propagation distance d (km) dependence of d^-1.32. The effect of a tall tower on the conversion factors from the EM field to the channel current has also been studied, by Zhang et al. (2014a, 2015b).

Triggered lightning has become an important technique in testing the direct and indirect effects of lightning. (Zhang et al., 2013a) investigated the induced overvoltage on a vertical signal line of an aerovane from an automatic weather station when lightning was triggered nearby. (Chen et al., 2016) analyzed the features of the residual voltage in the Surge Protective Devices (SPD) connected to an overhead distribution line and ground during triggered lightning. The residual voltage of the SPD linked to the overhead line was affected by the ground potential rise (GPR) at the SPD grounding point. The SPD residual voltage is usually determined by the induced voltage on the overhead distribution line by the RS and big M-component processes (Chen et al., 2018b). A long-duration induction current might come to the SPD from the ground because of the GPR effect. The SPD's residual voltage might last up to an order of milliseconds and damage the SPD. The average gradient of the triggering lightning current is a very important parameter in determining the peak current and surge energy through SPDs (Chen et al., 2018a). (Liu et al., 2015) proposed a tangential E-field technique based on the EM field to estimate the critical E-field of soil ionization, which is significant for the performance of a grounding system of lightning protection. The instant moment of the soil ionization can be found by checking the tangential component of the horizontal E-field.

Ground-based lightning location data and space-based lightning data have been widely used to study statistically the characteristics of lightning in different regions (Wu et al., 2013b, Wu et al., 2016b; Xia et al., 2015; Li et al., 2016; Zheng et al., 2016a), or lightning activities in different kinds of thunderstorms in China (Wu et al., 2013a; Xie et al., 2015; Zhang et al., 2017c).

Based on spaced-based Lightning Imaging Sensor and Optical Transient Detector (LIS/OTD) data, (Zhu et al., 2013) found that the global flash density was 46.2 fl s^-1 and nearly 78.1% of them occurred in the region of 30^°S-30^°N. The ratio of lightning density over land to that over the ocean was 9.6:1. Variation in lightning density with altitude was represented by two peaks and three valleys, with the former appearing at 100-2400 m and 3300-4600 m, and the latter appearing below 100 m, at 2400-3300 m, and above 4600 m. Zheng et al. (2016a, b) found that CG flashes over land showed only one peak in the afternoon, while those over offshore waters showed peaks in the morning and afternoon. The diurnal variations of small-current CG lightning (peak current ≤ 50 kA) and large-current CG lightning (peak current >50 kA, >75 kA and >100 kA) were much more different over the ocean, but similar over the continent.

The Northwest Pacific region is affected by diversified climatic systems and exhibits unique weather patterns. (Zhang et al., 2018b) analyzed the climatology of lightning activity over the Northwest Pacific using WWLLN data for the period 2005-15. It was found that the maximum lightning densities were along the coast of Southeast Asia and the tropical islands. The strongest lightning activity appeared in the monsoon season in the central and southern South China Sea. (Yuan et al., 2016) found that the lightning anomaly showed significant positive correlation with Oceanic NiÑo Index over both East China and Indonesia during El NiÑo episodes, while no obvious correlation was found during La NiÑa episodes. Significantly increased lightning activity was also found for larger convective available potential energy (CAPE).

The relationship between lightning and precipitation for short-duration rainfall events in the Beijing region was analyzed during the summer seasons from 2006 to 2007 by Wu et al. (2017, 2018). They developed a nowcasting method for short-duration rainfall (i.e., <6 h) events in terms of rainfall and lightning jumps. It was found that the approach provided early warnings for the associated short-duration rainfall events from the regional scale to meso-γ scale. (Xia et al., 2018) classified mesoscale convective systems (MCSs) into four categories based on their high/low convective rainfall rates (HR/LR) and high/low CG lightning frequencies (HL/LL). They found that the HRHL, HRLL, LRHL, and LRLL categories exhibited orders of the highest-to-smallest CAPE and perceptible water, but the smallest-to-largest convective inhibition and lifted indices.

Many unique features of lightning and charge structure in the Chinese Inland Plateau or Tibetan Plateau region have been discovered since the 1980s. Recently, in-situ E-field sounding and 3D lightning mapping results showed some promising new results on these issues.

(Li et al., 2017b) analyzed the evolution of the charge structure inferred from 3D VHF lightning mapping data for a small thunderstorm in Qinghai region. During the initial developing stage and the mature stage of the thundercloud, the charge structure exhibited a steady negative dipole, i.e., a negative above a positive charge region. During the dissipation stage of the thundercloud, the charge structure changed to a more complicated one, resulting from the merger of two convective cells. A positive-charge dipole, negative-charge dipole, and tripole charge structure coexisted in different regions of the thunderstorm during its late stage. (Fan et al., 2018b) found a positive tripole charge structure during most stages of a hailstorm with a larger LPCR in the thunderstorm, resulting in a low +CG rate. (Wang et al., 2013) found some kind of interaction between nearby lightning discharges, i.e., the following lightning discharges would be suppressed by the lightning discharges before. However, the prior lightning discharge might provide a preconditioned channel to facilitate the following lightning discharges. (Fan et al., 2014) developed an overlap and progressive method using either a one-point charge model or point dipole model, and found that the charge neutralized in the negative charge layer by the CC of lightning flashes usually ranged from 2.5 km to 4.7 km above the ground in Qinghai region.

By conducting balloon-borne E-field sounding in Pingliang, Gansu Province, (Zhang et al., 2018a) found five charge centers in the cloud region colder than 0^°C, and the polarity of charge regions alternated vertically, with the lowest region being positive. Two other charge layers were located near the cloud base, with the negative region below the positive. (Zhang et al., 2012d) suggested that the LPCR intensity of thunderstorms in the Chinese Inland Plateau region was likely determined by the quantity of graupel, based on the vertical distribution of each type of precipitation particle inferred from data from an X-band dual linear polarization Doppler radar.

Figure 7.  Schematic illustrations of the charge structure evolution of a typical plateau thunderstorm.

Figure 7 shows a schematic illustration of the charge structure evolution of a typical plateau thunderstorm based on previous studies. The outstanding characteristics of the plateau thunderstorm charge structure are the development of an inverted dipole (positive charge region below negative charge region) and the upper positive charge region may develop in the mature stage of the storm to form a tripole charge structure with a larger-than-usual lower positive charge region (refer to Figure 7a), or may not develop to form just an inverted dipole charge structure (refer to Figure 7b). In both cases, the inverted polarity IC flashes are observed frequently in the Chinese Inland and Tibetan plateau regions. In the later stage or the dissipating stage, the charge structure of the thunderstorm could be a normal dipole or normal tripole, with descent of the positively charged graupel particles from the lower part of the thundercloud.

Several cases of squall line systems have been analyzed based on SAFIR3000 or BLNet whole lightning data in the Beijing area (Liu et al., 2011, Liu et al., 2013a, b; Xu et al., 2016c; Sun et al., 2018a; Xu et al., 2018). (Liu et al., 2013b) found that most of the lightning radiation sources were located horizontally in the squall line with higher radar reflectivity. The distribution of radiation sources developed vertically from two layers into three layers. In the developing stage of the convective system, the upper-layer center of the radiation sources was located at about 11 km, while the lower layer centered at about 4 km. In the mature stage of the convective system, the entire convective line was characterized by a multilayer charge structure with three positive charge layers located at 5 km, 9.5 km, and 13 km, and two negative charge layers at 7 km and 11 km, respectively. (Xu et al., 2018) discussed the lightning activity in squall lines with a cell merging process, and found that lightning peaked after the cell merging in most of the studied cases. About 85% of lightning flashes were distributed within the 10 km range of the convection line, while quite a few were in the stratiform region. The majority of lightning was concentrated in the region with strong vertical wind shear (VWS), and the positive CG (PCG) was more likely to occur in the transition zone and stratiform region. (Wang et al., 2016b) found that the dominant polarity of the stratiform CG flashes was different from that of the CG flashes in the convective region. The current of the first RS of the stratiform CG flashes was usually greater than that of the convective CG flashes. The location of the stratiform CG flashes was always near or at the edge of a region with the bright band, but never directly below the reflectivity core in the bright band.

(Wang et al., 2017a) found that the grounding location of nearly 79.1% of CG flashes was underneath the cloud region with vertical velocity at the 0^°C layer ranging from -5 m s^-1 to 5 m s^-1, and the majority of the locations were under the weak updraft region. The bright band showed a close relationship with the lightning initiation near the 0^°C isotherm in the stratiform region (Wang et al., 2018a).

Some hailstorms have been found to contain two stages with high-rate lightning activity divided by a low lightning rate (Zheng et al., 2009; Xu et al., 2016c). The first peak is usually before the hailfall and the second peak after the hailfall, with the greater peak in the second stage. In two hailstorms, the dominant polarity of CG lightning around the first lightning peak was positive and changed to negative around the second lightning peak. The change in polarity resulted from the evolution of the charge structure from an inverted one during the first stage to a normal one during the second stage.

Supporting by 3D lightning location data with the capability of describing the channel structure of flashes, (Zhang et al., 2017e) reported that the mean flash size was inversely correlated to the flash rate, with a correlation coefficient of -0.87 following a unary power function. It was more likely that thunderstorms with relatively weak convection and high precipitation corresponded to CG lightning with a large peak current. The charge pockets and horizontally broad charge regions were employed to explain the relationships of flash density with flash size or intensity (Zheng et al., 2016b).

Narrow bipolar events (NBEs), sometimes referred as compact intracloud discharges (CIDs), are a distinct type of intercloud lightning flashes characterized by a narrow bipolar EM pulse with a timescale of 10 μs to 20 μs, and produce powerful VHF radiation.

(Liu et al., 2012) mapped the discharge channel images of CIDs by using VHF interferometers and indicated that the CIDs tended to develop vertically with a scale ranging from 0.40 km to 1.9 km and an average velocity from 0.44× 10⁸ m s^-1 to 1.0× 10⁸ m s^-1. (Wang et al., 2012d) found 32 out of 236 CIDs in isolation, while 204 of them were accompanied by either IC or CG flashes. Among the latter, a total of 130 were more likely to trigger lightning, while 72 were in the middle of lightning discharges, and the other two terminated lightning. The peak power of the CIDs always occurred between an altitude of 7 km and 16 km, and ranged from 12 kW to 781 kW in the frequency band from 267 MHz to 273 MHz.

(Wu et al., 2011) found that the ratio of -CIDs appears to increase with the intensity of the convection. (Zhang et al., 2016a) developed a 3D locating method for NBEs based on a two-axis magnetic sensor and ionospheric reflection pairs. They found that NBEs were mainly produced in the convective region with an echo intensity ≥ 30 dBZ. Positive NBEs were predominately produced between 7 km and 15 km, while -NBEs were above 14 km. The NBEs that occurred in Northeast China were observed to be consistent with the results of those in other regions (Lü et al., 2013a, b)

Tropical cyclones (TCs) are a major weather process in the Northwest Pacific. (Pan et al., 2010) found a peak value of lightning flashes in the eyewall region, a minimum in the inner rain band, and a strong maximum in the outer rain band for seven super typhoons over the Northwest Pacific. This bimodal pattern of lightning distribution has been verified by recent research (Pan et al., 2013, Pan et al., 2014; Wang et al., 2016a, Wang et al., 2017b, Wang et al., 2018b). (Zhang et al., 2015d) analyzed the lightning activities in 116 TCs in the years from 2005 to 2009 in the Northwest Pacific using WWLLN lightning data. They indicated that weak storms (tropical depressions and tropical storms) were more likely to produce lightning than very strong storms (typhoons and super typhoons). TCs that strengthened in the next 24 h exhibited higher lightning density than those that weakened. The lightning density inside the inner core was larger for rapid intensification cases than for rapid weakening cases, suggesting a predictive effect of inner-core lightning on TC rapid intensity change. (Pan et al., 2014) found that the correlation coefficient between the total flashes within a radius of 600 km and the maximum wind speed for weak and super typhoons from 2005 to 2009 was 0.81 and 0.74, respectively. For about 78% of the super typhoons, the peak lightning rate was prior to the maximum wind speed, and the most common lead time was about 30 h, while it was 60 h for 56% of weak typhoons, indicating that the lightning rate could be used as a measure of typhoon intensification.

(Wang et al., 2018b) found that the variation in flash rate with sea surface temperature showed a positive correlation for 280 TCs from 2005 to 2016 over the Northwest Pacific. The VWS dominated the downshear left (right) asymmetric lightning distribution of the inner core (outer rain bands), and this asymmetry became more significant as the wind shear increased. The asymmetry of lightning distribution in the inner core was related to the joint effects of the VWS and TC motion vectors; however, the asymmetry in the outer rainband was more closely related to the VWS. The inner core asymmetry weakened when the TC was moving fast in the direction that was opposite or to the right of the VWS.

Numerical simulation of lightning and electrification inside thunderstorms has been conducted in several groups (Tan et al., 2012, Tan et al., 2014c, Tan et al., 2016; Wang et al., 2016c; Guo et al., 2017). (Liu et al., 2014) adopted the noninductive electrification mechanism, i.e., the Takahashi78 scheme and Saunders98 scheme, and lightning parameterization into the RAMSV6.0 model. The simulation showed that the thunderstorm exhibited a tripole charge structure with the Takahashi78 scheme, while the charge structure changed from dipole to tripole with the Saunders98 scheme. The simulated lightning rate was in good agreement with the observation.

Xu et al. (2012, 2014) and (Zhao et al., 2015) introduced electrification schemes into the WRF model, and developed WRF-Electric models. The numerical experiments on the charge structure evolution of a hailstorm suggested that the inverted charge structure could be produced through a so-called dynamical-derived mechanism (Xu et al., 2016a). Simulation of Typhoon Molave (2009) showed that the charge structure in the eyewall exhibited a positive tripole prior to its landfall by using WRF-Electric (Xu et al., 2016b). It became a negative dipole with negative charge in the middle and positive charge below after Molave reached its maximum intensity. The charge structure in the eyewall was closely related to the typhoon intensity, but not directly correlated to the landfall. The convective cells in the outer rain band exhibited a positive tripole or positive dipole charge structure in different stages of Typhoon Molave.

Wang et al. (2015a, b) found that ice particles with vertical velocity ranging from 1 m s^-1 to 5 m s^-1 accumulated the most charge during all stages of a thunderstorm. An updraft with speeds between -1 m s^-1 and 1 m s^-1 was the most favorable for charge separation. The vertical velocity at the initiation locations of lightning flashes was correlated to the maximum updraft speed with a cubic polynomial. The graupel mixing ratios at the initiation sites correlated linearly with the mixing ratios of the graupel concentration center. This linear correlation was significant in the active and the later stage of the thunderstorm (Wang et al., 2017c).

Aerosol effects on thunderstorm electrification and lightning frequency have been studied using numerical simulation (Shi et al., 2015; Zhao et al., 2015). (Zhao et al., 2015) found that increased aerosol loading resulted in a large increasing rate of snow and graupel, as well as a larger ice particle density, and the electrification process in the thunderstorm enhanced consequently.

(Tan et al., 2014a) studied impacts of the LPCR on lightning type using a lightning-resolving model. They found that the LPCR was crucial in producing -CG and inverted -IC flashes. Increasing the LPCR charge density or extension, lightning changed from +IC-dominated to -CG lightning-dominated, and then to inverted IC flashes-dominated. The charge density of the LPCR played a major role in the determination of lightning type compared with extension of the LPCR. It was only when the LPCR maximum charge density was within a certain range that -CG flashes occurred, and the probability of -CG flashes was almost constant. The potential at the lightning initiation point was a crucial parameter to determine whether the leader propagated to the ground (Tan et al., 2014b). The initiation potential for CG lightning was >30 MV, while that for IC lightning was typically <30 MV. Because the space charge determines the environmental potential, the lightning type also depends on the relative position and magnitude of charge regions near the lightning initiation position. (Sun et al., 2018b) discussed the feedback effect of E-field force on electrification and charge structure in thunderstorms based on the NSSL WRF-Electric model; the overall influences of the E-field force on electrification tended to be positive, and its feedback effect on the thunderstorm structure should not be neglected.

Lightning data assimilation methods have been investigated and tested in cloud-resolving numerical models to improve the forecasting of convection and precipitation. (Qie et al., 2014b) constructed an empirical formula between lightning rate and the mixing ratio of the cloud ice particle, including graupel, ice and snow. The established functions were nudged into the WRF model. The intensity and location of the strong convection were significantly improved 1 h after the lightning assimilation. The precipitation center, amount and coverage were improved too, and they were all much closer to the observation than that without lightning assimilation. (Chen et al., 2017) developed a new lightning data assimilation scheme with comprehensive nudging of water content. Both the low-level water vapor and graupel mass in the mixing phase region were nudged according to the detected total lightning flash rate and model conditions. The bulk Richardson number was adopted to measure the dynamic and thermodynamic conditions and as an index to nudge the lightning. The simulation showed that the new scheme better matched the observation in terms of convection and precipitation. (Zhang et al., 2017b) assimilated lightning data using the 3D variational data assimilation method in the WRF-3DVar system in cycling mode with an interval of 10 min based on empirical functions between lightning frequency and the mixing ratio of water vapor. The results for a squall line system indicated that 60 min was the appropriate time-window for lightning assimilation. Forecasting of accumulated precipitation in 1 h during the assimilation and the accumulated precipitation in the following 3 h were greatly improved.

(Wang et al., 2014b) used the CG lightning rate to adjust the vertical velocity, specific humidity and specific cloud water content in a model using a physical initialization method, based on a relationship between lightning flash density and reflectivity. The forecasted convection was improved with the lightning-proxy reflectivity assimilation scheme. The reflectivity prediction was significantly improved and maintained for about 3 h.

(Wang et al., 2018c) developed a lightning assimilation scheme based on time-lagged ensembles. With the assimilation method, the background error covariance was calculated using time-lagged ensembles. The graupel mixing ratio was retrieved from the observed lightning rate by using empirical vertical profiles (Wang et al., 2017d). The observation errors were estimated by using uncertainties in lightning data and the retrieved vertical profiles. Assimilating lightning data for a severe convective event exhibited many observed convective cells that did not appear in the control run; plus, it restrained the appearance of false convection and decreased the shifting errors of convection.

(Li et al., 2016) simulated the occurrence of lightning by using a modified lightning potential index and observed lightning data. The lightning density was calculated by differentiating the ice mass of precipitation and non-precipitation. The proposed method was examined for a quasi-linear MCS using the WRF model and 3D variation analysis system in the ARPS model. The simulation suggested that most lightning flashes occurred on the right side and at the front of the MCSs, where the surface wind field converged intensely. The lightning flashes tended to occur in the regions with a large gradient of CAPE.

The middle and upper atmosphere including lower ionosphere is dynamically perturbed by the underlying thunderstorms. Observations on transient luminous events (TLEs), including sprites, elves, blue jets, blue starters, halos and gigantic jets (GJs) occurring above thunderstorms have been continuously conducted since 2007. The relationship of TLEs with the parent lightning and thunderstorm has been the main subject in the last few years. In addition, effects of the thunderstorms on the ionosphere have also been investigated.

(Yang and Feng, 2012) reported a GJ produced by a multi-cell thunderstorm in the coast region (Yellow Sea) in eastern China. The relationship between the GJ occurrence time and the strong convection has been confirmed. However, based on data analysis, (Yang et al., 2018a) also found that the GJ was produced during summer thunderstorms in the midlatitude region without strong updrafts, as indicated by the maximum echo top along the GJ azimuth being lower than the regional tropopause. Different from other summer thunderstorms just producing GJs during their evolutions, two sprites were also observed in a time window containing the GJ in the study of (Yang et al., 2018a), suggesting that the meteorological characteristics of the GJ-producing storms has not yet been fully understood. By using data from the lightning imager onboard the FORMOSAT-2 satellite and a ground-based lightning location network, (Liu et al., 2018a) suggested that the negative NBEs were probably the initial part of the blue jets, as evidenced by the fact that the negative NBEs were concurrent with six of the seven blue jets by less than 1 ms.

Sprites are the most easily observed TLEs in ground-based observation. A sprit-producing thunderstorm usually produces several or even many sprites. However, compared with thunderstorms that may induce hundreds of TLEs in the High Plain region in the United States, thunderstorms in China are usually observed to produce several or tens of sprites (Yang et al., 2013a, b,; Wang et al., 2015d; Huang et al., 2018), and some storms even produce just a single sprite during its life cycle (Yang et al., 2017a). No significant differences in the particle (i.e., precipitation ice, cloud ice and cloud water content) distributions have been found among many-, one- and non-sprite-producing storms. Most sprites are caused by positive strokes (Yang et al., 2013a, b,) and the following CC (Yang et al., 2015, 2017a), but they can also be initiated by -CG flashes (Yang et al., 2017b). Compared with many positive sprites, negative sprite events are rare. (Yang et al., 2018b) reported five sprites over an MCS, and one of them was confirmed to be negative by both the local LLS and ELF magnetic field, and another was only confirmed to be negative by the local LLS data. They also found that sprites with different polarity were located in distinct regions of the MCS; the negative sprites were associated with strong convection and large wind shear, and positive ones were above the stratiform region. Negative sprite events could be produced by lightning strokes with a not very large peak current.

The ionosphere's characteristics can be affected by underlying thunderstorms. (Yu et al., 2015) found that enhancement of the ionospheric sporadic E-layer likely resulted from powerful lightning strokes. (Yu et al., 2017) showed a significant increase in the neutral metal Na layer above thunderstorms, and the thunderstorm-associated gravity waves and E-field could be possible reasons. (Xu et al., 2015) studied two concentric gravity wave events launched by underlying thunderstorms. They showed that the horizontal waves of the first event propagated freely within 300 km from the thunderstorm center, but the second event was induced by two strong thunderstorms exhibiting multiscale waves with different horizontal wavelengths.

High-energy particles related to thunderstorms are also a very interesting topic. (Wang et al., 2012c) found that the E-field of thunderstorm could affect the counting rate of the Neutron Monitor at Yangbajing over Tibetan Plateau region. Because of absorption in the atmosphere, ground-based measurement of high-energy particles is difficult; coordinated ground- and space-based observations and theoretical work are needed to further investigate this issue.

Although substantial progress has been made in recent years in understanding lightning discharge processes and their effects, a number of mysteries still remain. As described above, new detection technology holds promise for future advancements in our understanding of lightning and its relationships with other subjects. The following is a summary of some of the outstanding scientific questions that remain to be addressed:

The first and perhaps most interesting question concerns the detailed nature of lightning initiation and its relationship with terrestrial gamma-ray flashes (TGFs), as well as energetic in-cloud pulses (EIPs) with ultra-high currents >200 kA. It is becoming evident that TGFs are likely to be produced during initial breakdown processes, and EIPs could be a highly energetic form of initial breakdown pulses. The relationship between TGFs and lightning initiation and the detailed processes that produce TGFs are still the largest mystery and should be paid more attention in China.

The second outstanding question concerns the nature of fast breakdown associated with high-power NBEs, which also seems to be responsible for initiating lightning. VHF mapping interferometer and 3D VHF location joint observations with continuously improving capability will play a critical role in revealing more features of fast breakdown and answering how fast breakdown is produced.

The third question concerns how kinematic and microphysical properties of storms affect electrification and lightning. In particular, how does ice microphysics affect electrification processes? Polarimetric radars, in-situ particle and E-field sensors, together with improving broadband interferometers and 3D lightning mapping systems hold promise for addressing these outstanding questions. The FY-4 Geosynchronous Lightning Mapping Imager will allow us for the first time to study lightning continuously in geographic regions that have been inaccessible to ground-based lightning location networks. It is expected to allow us to make considerable progress in our understanding of lightning in oceanic thunderstorms, in TCs and in mountainous thunderstorms.

In conclusion, figuring out the various facets of the mysteries of lightning is obviously a primary subject of interest, and needs to be continuously studied by the atmospheric electricity community.