Cloud-to-ground Research Articles

Initiation of Downward Positive Leader Beneath the Negative Leader Channel

AbstractDue to the weak radiation generated by positive leaders, our understanding of how positive leaders are initiated and reach the ground remains limited. This study investigated positive cloud‐to‐ground (CG) lightning induced after the long intracloud lightning based on high‐speed video and fast antenna mapping results. Below the cloud base in the field of view, re‐breakdowns of decayed negative branches occurred consecutively beneath the horizontal negative channel in the form of bidirectional leaders. These bidirectional leaders advanced along the same channel, and eventually, the third reached the ground, generating a positive return stroke (RS) with a peak current of 157 kA. Following the RS, new negative discharges emerge adjacent to the vertical return stroke channel, persisting and propagating to form a continuing current (CC) lasting over 200 milliseconds.

Virtual Reflection Height of Nighttime Equatorial Ionosphere Estimated With Low‐Frequency Magnetic Sferics Measured in Malacca

AbstractThe return stroke of cloud‐to‐ground (CG) lightning is an impulsive radiator of very low‐frequency/low‐frequency (VLF/LF) electromagnetic signals allowing for the remote sensing of lower ionosphere over large spatial coverage. In this study, we examined the LF magnetic fields measured in Malacca, Malaysia, to probe reflection heights of the lower ionosphere near the equator on three different nights in 2021. The results show that the virtual ionospheric height at nighttime typically ranged from 82.0 to 90.0 km, with a mean value of 85.3 km. Our measurements also revealed significant variations in the virtual ionospheric height across different regions over a spatial scale of about 800 km. The maximum height difference was about 5.0 km. Moreover, the fluctuation characteristics are observed in both estimated ionospheric height and calculated peak reflection ratio during similar periods. This fluctuation may be related to atmospheric gravity waves in the nighttime ionosphere. In addition, we compared the virtual ionospheric height estimated from CG strokes of different polarities, and the results showed that the virtual reflection height for positive CG strokes is lower than that for negative ones.

Total lightning signatures in a tornadic thunderstorm over the Pearl River Delta of Southern China

Tornado rises rotating air columns extending from the bottom of cumulonimbus clouds to the ground, often accompanied by strong convection with thunderstorms, hail, and short-term heavy precipitation, which seriously threaten people's lives and property safety. Studying the characteristics of thunderstorm activity during a tornado is crucial for comprehending the atmospheric electrical mechanisms involved in its generation process and developing cooperative methods for tornado warning and forecasting. Based on the VLF/LF total lightning location system and radar echo data, this paper analyzes the spatial-temporal evolution, frequency, polarity, and height of total lightning during a strong tornado (EF3) in the Pearl River Delta region on June 3, 2014. The lightning activity lasted for about 3.5 h, with a total of 41,117 total lightning flashes, of which intra-cloud (IC) flashes accounted for 74.7%, cloud-to-ground (CG) flashes accounted for 21.6%, and narrow bipolar events (NBEs) accounted for 3.7%. Connected-Component Labeling (CCL) algorithm was used to divide the thunderstorm into four stages: initiation, development, maturation, and dissipation, and it was observed that the occurrence of tornadoes was closely related to total lightning activities. The observed characteristics of lightning activity in the tornado are summarized through comparison with other studies.

Evaluation of Detection Efficiency of World Wide Lightning Location Network in Southeast Asian Region

AbstractThe detection efficiency (DE) of the World Wide Lightning Location Network (WWLLN) is evaluated in Southeast Asia by comparing WWLLN data with the Earth Network Total Lightning Network (ENTLN) data taking into account time, distance, and peak‐current parameters. The performance of WWLLN is evaluated in the months of November and December in two different years (2020–2021). These periods are selected to assess the change (if any) in DE overtime and the inclusion of new stations. The strokes between the two networks were considered matched if they fell within a time criterion of 100 µs and a location difference of 25 km. Using this criterion, 5.2 × 106 WWLLN strokes were matched with ENTLN cloud‐to‐ground (CG) lightning strokes in November‐December 2020, resulting in a DE of 32.9%. Similarly, 4.6 × 106 WWLLN strokes were found to match in November‐December 2021, yielding a DE of 36.5%. Analysis of the peak‐currents reveals that DE is lowest (<10%) for a peak‐current below ±10 kA. However, for peak‐current exceeding ±50 kA, the DE increases to ∼60%. During November‐December 2021, WWLLN reported 38.95 × 106 lightning strokes globally; amongst them, Dhaka station detected 0.5 × 106 strokes, contributing to a 1.3% increase in the global DE. Dhaka station detects most lightning strokes within 8 × 103 km, which diminishes to zero after 10 × 103 km. The Dhaka station recorded a larger number of strokes at longer detection distances during midnight (22:00–02:00) than during noon (10:00–14:00). The results signify a positive impact of the Dhaka station on WWLLN's DE during the mentioned period.

Temporal characteristics of positive cloud-to-ground flashes from tropical thunderstorms

The occurrence of positive cloud-to-ground (CG) lightning flashes proves to be a challenging task when compared to the negative CG flash which common happen in tropical region. This due to the positively charged clouds to reach higher altitudes to enable the transfer of electrical charges to the Earth's surface. With several sensors that deployed in the Malacca region (2.314077°N, 102.318282°E) such as the wide band fast and slow antenna, magnetic field sensors, and the electric field mill (EFM-100). This research study focuses on the temporal characteristics of positive CG lightning activity in the tropical region during September and October 2021. A total of 203 positive cloud-to-ground flashes were meticulously analyzed in the context of 16 thunderstorm events. Key findings highlight the consistent initiation of positive cloud-to-ground flashes through initial breakdown (IB) pulses, followed by the development of stepped leaders and subsequent return strokes. Notably, most positive cloud-to-ground events (72%) were observed as single strokes and the maximum multiplicity recorded being four subsequent return strokes. Furthermore, the analysis focuses into the temporal aspects of these lightning events, revealing average durations for various parameters, including rise-time (6.26 μs), zero-crossing (26.43 μs), pulse duration of the first RS (139.76 μs), and the interval from the first IB to the first return stroke (134.19 ms). Additionally, the recorded data from the Electric Field Mill (EFM) instruments together with Constant Altitude Plan Position Indicator (CAPPI) radar demonstrate a notable correlation. Specifically, the electric field (E-field) values exhibit a discernible increase when positive cloud-to-ground lightning events are detected near the lightning sensor.

Investigating Temporal Characteristics of Polarimetric and Electrical Signatures in Three Severe Storms: Insights from the VORTEX-Southeast Field Campaign

Abstract The dual-polarization radar characteristics of severe storms are commonly used as indicators to estimate the size and intensity of deep convective updrafts. In this study, we track rapid fluctuations in updraft intensity and size by objectively identifying polarimetric fingerprints such as ZDR and KDP columns, which serve as proxies for mixed-phase updraft strength. We quantify the volume of ZDR and KDP columns to evaluate their utility in diagnosing temporal variability in lightning flash characteristics. Specifically, we analyze three severe storms that developed in environments with low-to-moderate instability and strong 0–6-km wind shear in northern Alabama during the 2016–17 VORTEX-Southeast field campaign. In these three cases (a tornadic supercell embedded in stratiform precipitation, a nontornadic supercell, and a supercell embedded within a quasi-linear convective system), we find that the volume of the KDP columns exhibits a stronger correlation with the total flash rate. The higher covariability of the KDP column volume with the total flash rate suggests that the overall electrification and precipitation microphysics were dominated by cold cloud processes. The lower covariability with the ZDR column volume indicates the presence of nonsteady updrafts or a less prominent role of warm rain processes in graupel growth and subsequent electrification. Furthermore, we observe that the majority of cloud-to-ground (CG) lightning strikes a carried negative charge to the ground. In contrast to findings from a tornadic supercell over the Great Plains, lightning flash initiations in the Alabama storms primarily occurred outside the footprint of the ZDR and KDP column objects. Significance Statement This study quantifies the correlation between mixed-phase updraft intensity and total lightning flash rate in three severe storms in northern Alabama. In the absence of direct updraft velocity measurements, we use polarimetric signatures, such as ZDR and KDP columns, as proxies for updraft strength. Our analysis of polarimetric radar and lightning mapping array data reveals that the lightning flash rate is more highly correlated with the KDP column volume than with the ZDR column volume in all three storms examined. This contrasts with previous findings in storms over the central Great Plains, where the ZDR column volume showed higher covariability with flash rate. Interestingly, lightning initiation in the Alabama storms mainly occurred outside the ZDR and KDP column areas, contrary to previous findings.

Observation of tropical positive cloud-to-ground flashes accompanied by chaotic and regular pulse trains

This study investigates the occurrence of chaotic pulse trains (CPTs) and regular pulse trains (RPTs) in tropical positive cloud-to-ground (CG) lightning flashes. These flashes are categorized into four types based on the initial polarity of the initial breakdown (IB) pulses and their relationship to the first return stroke (RS). A total of 71 positive CG flashes from five thunderstorm events were analyzed. The analysis reveals instances of CPTs and RPTs both before and after the first positive RS, along with the occurrence of mixed polarities in RPTs. Variations in IB pulse polarities and the presence of CPTs and RPTs before the first positive RS were observed, contrasting with previous findings in negative CG flashes. All positive CG flashes have been detected when cloud top height occurrences were between 12 and 18 km. In contrast, for negative CG flashes with CPTs and RPTs the cloud top height occurrences were between 5 and 12 km. It is interesting that CPTs and RPTs can be detected during IB process of positive CG flashes at relatively high altitude in the thundercloud. Perhaps due to low pressure at higher altitudes in the cloud, electrical process associated with CPTs and RPTs are easily discharged before the occurrence first positive return stroke. The altitudes of cloud top heights for the inverse polarity of IB pulses were located between 16 and 18 km. This research enhances the understanding of positive CG lightning initiation process and their relationship with CPTs and RPTs, as well as the occurrence of recoil leaders.

Inter-comparison of World Wide Lightning Location Network (WWLLN) and Lightning Detection Sensor Network (LDSN) data over India

Atmospheric lightning is a dominant weather related disaster which requires continuous monitoring efforts. While the world wide lightning location network (WWLLN) measures the very low frequency (VLF) signals emanating from the lightning, other regional networks, such as Lightning Detection Sensor Network (LDSN) use the broad band radio frequency signals. Present investigation compares the cloud-to-ground (CG) lightning occurrences recorded by the WWLLN and LDSN during 02 October – October 09, 2022 over India. Results show that though the WWLLN monitored the CG lightning occurrences about 20% lesser than the LDSN, spatial distribution of the occurrences and trends were in very good agreement.

A 3D Numerical Model to Estimate Lightning Types for PyroCb Thundercloud

Pyrocumulonimbus (pyroCb) thunderclouds, produced from extreme bushfires, can initiate frequent cloud-to-ground (CG) lightning strikes containing extended continuing currents. This, in turn, can ignite new spot fires and inflict massive harm on the environment and infrastructures. This study presents a 3D numerical thundercloud model for estimating the lightning of different types and its striking zone for the conceptual tripole thundercloud structure which is theorized to produce the lightning phenomenon in pyroCb storms. More emphasis is given to the lower positive charge layer, and the impacts of strong wind shear are also explored to thoroughly examine various electrical parameters including the longitudinal electric field, electric potential, and surface charge density. The simulation outcomes on pyroCb thunderclouds with a tripole structure confirm the presence of negative longitudinal electric field initiation at the cloud’s lower region. This initiation is accompanied by enhancing the lower positive charge region, resulting in an overall positive electric potential increase. Consequently, negative surface charge density appears underneath the pyroCb thundercloud which has the potential to induce positive (+CG) lightning flashes. With wind shear extension of upper charge layers in pyroCb, the lightning initiation potential becomes negative to reduce the absolute field value and would generate negative (−CG) lightning flashes. A subsequent parametric study is carried out considering a positive correlation between aerosol concentration and charge density to investigate the sensitivity of pyroCb electrification under the influence of high aerosol conditions. The suggested model would establish the basis for identifying the potential area impacted by lightning and could also be expanded to analyze the dangerous conditions that may arise in wind energy farms or power substations in times of severe pyroCb events.

Observations of very unusual behavior of lightning discharges

Based on recent observations at the Lightning Observatory in Gainesville (LOG), Florida, we present three examples of very unusual behavior of cloud-to-ground (CG) lightning flashes: (a) One +CG and one tower-terminated -CG separated by 11 km, were briefly coupled, so that positive charge was in effect drawn from the ground at the position of –CG and injected into the ground at the position of +CG. That unusual interaction was accomplished by the initiation of M-component incident wave in the +CG channel by the leader of the first stroke of the –CG. (b) A branched positive leader intermittently extended via bidirectional transients separated by inactive intervals of about 5 ms, on average. One of the bidirectional transients, moving negative charge in the backward direction, triggered a downward negative leader, resulting in a three-stroke –CG. (c) A subsequent-stroke negative leader, moving through virgin air, entered the lower part of residual (non-luminous) first-stroke channel, in addition to creating a new ground termination about 950 m away. Studying unusual lightning behavior helps to improve our understanding of the variability of basic lightning processes and to identify new potential lightning hazards to various objects and systems.

The impact of aerosols as ice nucleating particle on microphysics and electrification in the cumulus model

Here, a heterogeneous ice nucleation parameterization associated with aerosol acting as ice nucleating particle (INP) was implanted into a two-dimensional numerical cumulus model. To explore the impact of INP on microphysical and electrical processes, a comparison was conducted with the original approach, which employs an empirical formula. Simulation results indicate that INP greatly impacted microphysical evolution in the heterogeneous ice nucleation process, reducing total liquid precipitation amounts and causing a slight precipitation delay, as well as increasing the diameter and mixing ratio of ice crystals and the expansion of the vertical distribution of ice crystals. This led to a notable change in electrification in thunderstorms. The increase of ice crystal diameter was the dominant contributor to the enhancement of electrification using the new parameterization. Additionally, unlike the structure of thunderclouds in a mature stage, which always retains a normal dipole structure adopting the empirical formula, a tripole structure developed a lower positive charge, and the polarity inversion with upper negative and lower positive charges occurred when the new parameterization was adopted. This predominately was the result of abundant ice crystals present below the reversal temperature. The electrification characteristics of thunderstorms may have a close connection with lightning activity. It has been found that charge structure changed significantly in the two cases, with the tripolar charge structure facilitating the production of inverted intra-cloud (IC) flashes and negative Cloud-to-ground (CG) flashes simulated by the new parameterization; additionally, clouds simulated using the empirical formula may be able to develop normal IC lightning and positive CG flashes. Therefore, it will be very meaningful to obtain greater insight into the characteristics of thunderstorms electrification in cumulus model with aerosols.

Study of dynamics-microphysical-lightning activity characteristics in a tropical hailstorm

Based on three-dimensional lightning locations and S-band polarimetric Doppler radar data, a hailstorm that occurred on Hainan Island on June 4, 2021 was analyzed in detail. The analysis revealed that the flash activity of hailstorms is mainly intracloud (IC) flashes, the cloud-to-ground (CG) flashes frequency were less, positive CG flashes (+CG) were domination in the CG activity before hailfall. The flash frequency appeared “lightning jumps” when hailstorms merging with convective cells, ice hydrometeors also rapidly increased in the clouds of the mixed-phase area in this period, which was beneficial for efficient electrification. The Zdr column and evolution of echo top indicated presence of an intense updraft during the mature stage of the hailstorm. Zdr column depth about 4 km close to 20 °C during the mature stage. Near the time of the hailfall, Zdr column was contraction, due to sinking hail generated downdraft within the cloud. Zdr column after hailfall may be produced by the convergence of downdraft and sea breeze. In combination with the spatiotemporal distribution of the hydrometeors in hailstorms, flash activity was most closely related to graupel and it plays a important role for flashes electrification. The hail plays a minimum role in charge separation. Radiation sources were concentrated between 50 dBZ and 55 dBZ. The spatial distribution of radiation sources indicated possible complex charge structure within hail storm clouds.

Identification of cloud-to-ground lightning and intra-cloud lightning based on their radiated electric field signatures using different types of neural networks and machine learning classifiers

The scope of this paper is to test the efficiency of three different neural networks in the area of cloud-to-ground/intra-cloud lightning classification using a not very large data set of electric field data. Effective neural network learning verified for a small lightning event dataset but with whole 2-seconds particular record length is an important feature of fast and efficient procedures implemented into lightning location systems among which discrimination between cloud-to-ground and intra-cloud lightning is one of the most important task. Measurement data obtained at the Lightning Observatory in Rzeszow were used for the analysis. About 100 extremely low/medium frequency bandwidth electric field of cloud-to ground lightning registrations and the same number of intra-cloud lightning events were taken into account. During the study 81% of those registrations were dedicated for learning of neural networks while the remaining 19% were put into the testing set. All data were preselected manually with consideration of the electric field variation and information from the lightning location system database. Lightning events were selected uniformly with the well-pronounced presence of their basic components and reported distance to the registration station which ensured a wide variation of expected lightning electric field signatures. The multilayer perceptron neural network, the radial basis function neural network and the convolutional neural network were tested and compared in different configurations while distinguishing lightning. The results were compared with other conventional classification approaches, such as traditional machine learning methods as well as one a little more contemporary architecture, the long short-term memory neural network. The multilayer perceptron neural network achieved the best detection accuracy overall. Presented research is an extension of lightning identification studies available so far mainly based on the convolutional neural network by the results achieved using the multilayer perceptron neural network, the radial basis function neural network, the long short-term memory neural network and several machine learning-based classification methods.

Assessing Flash Characteristics in Lightning-Initiated Wildfire Events between 1995 and 2020 within the Contiguous United States

Abstract Twenty-six years of lightning data were paired with over 68 000 lightning-initiated wildfire (LIW) reports to understand lightning flash characteristics responsible for ignition in between 1995 and 2020. Results indicate that 92% of LIW were started by negative cloud-to-ground (CG) lightning flashes and 57% were single stroke flashes. Moreover, 62% of LIW reports did not have a positive CG within 10 km of the start location, contrary to the science literature’s suggestion that positive CG flashes are a dominant fire-starting mechanism. Nearly ⅓ of wildfire events were holdovers, meaning 1 or more days elapsed between lightning occurrence and fire report. However, fires that were reported less than a day after lightning occurrence statistically burned more acreage. Peak current was not found to be a statistically significant delineator between fire starters and non–fire starters for negative CGs but was for positive CGs. Results highlighted the need for reassessing the role of positive CG lightning and subsequently long-continuing current in wildfire ignition started by lightning. One potential outcome of this study’s results is the development of real-time tools to identify ignition potential during lightning events to aid in fire mitigation efforts.

A 13-year long strokes statistical analysis over the Central Mediterranean area

This paper presents the first detailed analysis of cloud-to-ground (CG) and intra-cloud (IC) strokes characteristics from the Lightning Detection Network (LINET) over Italy and the Central Mediterranean area, a lightning active area in south Europe. We study the strokes over a 13-year period from 2010 to 2022, aiming to understand how it varies with different temporal scales (hourly, monthly, seasonally, and yearly), surface types (sea and land), and ground levels (0–100 m; 100–200 m; 200–400 m; 400–800 m; 800–1200 m; 1200–2000 m and above 2000 m).We found that the stroke's maximum activity was observed in August; specifically, July has the maximum activity over the land with a maximum diurnal peak in the afternoon, while in September, the convection shifts over the sea with a secondary daily maximum in the morning. The largest current intensities are observed in January, over sea and during nighttime. Moreover we found that stroke current intensities, polarity and IC height emissions are influenced by ground altitude level. Our paper provides new insights into the spatio-temporal patterns and characteristics of lightning over Italy and the Central Mediterranean area, which can be useful for improving weather forecasting, climate modeling, risk assessment, and damage mitigation strategies in this area.

How Much Lightning Actually Strikes the United States?

Abstract The number of cloud-to-ground (CG) flashes over the contiguous United States (CONUS) has been estimated to be from as small as 25 million per year to as many as 40 million. In addition, many CG flashes contact the ground in more than one place. To clarify these values, recent data from the National Lightning Detection Network (NLDN) have been examined since the network is performing well enough to make precise updates to the number of CG flashes and their associated ground contact points. The average number of CG flashes is calculated to be about 23.4 million per year over the CONUS, and the average number of ground contact points is calculated as 36.8 million per year. Knowledge of these two parameters is critical to lightning protection standards, as well as better understanding of the effects of lightning on forest fire initiation, geophysical interactions, human safety, and applications that benefit from knowing that a single flash may transfer charge to the ground in multiple, widely spaced locations. Sensitivity tests to assess the effects of misclassification of CG and in-cloud (IC) lightning are also made to place bounds on these estimates, and the likely uncertainty is a few percent.

Cloud-to-Ground Lightning in Cities: Seasonal Variability and Influential Factors

AbstractUrban-induced land use changes have a significant impact on local weather patterns, leading to increased hydro-meteorological hazards in cities. Despite substantial threats posed to humans, understanding atmospheric hazards related to urbanisation, such as thunderstorms, lightning, and convective precipitation, remains unclear. This study aims to analyse seasonal variability of cloud-to-ground (CG) lightning in the five large metropolitans in Bangladesh utilising six years (2015–2020) of Global Lightning Detection Network (popularly known as GLD360) data. It also investigates factors influencing CG strokes. The analysis revealed substantial seasonal fluctuations in CG strokes, with a noticeable increase in lightning activity during the pre-monsoon months from upwind to metropolitan areas across the five cities. Both season and location appear to impact the diurnal variability of CG strokes in these urban centres. Bivariate regression analysis indicated that precipitation and particulate matter (PM) significantly influence lightning generation, whilst population density, urban size, and mean surface temperature have negligible effects. A sensitivity test employing a random forest (RF) model underscored the pivotal role of PM in CG strokes in four of the five cities assessed, highlighting the enduring impact of extreme pollution on lightning activity. Despite low causalities from CG lightning, the risk of property damage remains high in urban environments. This study provides valuable insights for shaping public policies in Bangladesh, a globally recognised climate hotspot.

Tracking Lightning Through 3D Thunder Source Location With Distributed Acoustic Sensing

AbstractInvestigating lightning is of key significance in understanding the lightning mechanism and mitigating lightning hazards. We reported an experiment of investigating lightning through three‐dimensional (3D) thunder locating using a Distributed Acoustic Sensing (DAS) array in Hefei, China. In this experiment, we recast a 7.7 km long urban telecom optical fiber cable as 3,850 sensors using the DAS technique. From dense DAS recording, we manually identified 101 thunder events during six positive cloud‐to‐ground (CG) lightning flashes within 20 min. The DAS recorded thunder signals are dominated by direct acoustic waves rather than air‐ground coupled surface waves. The thunder events were then located using the arrival times of thunder signals. The locations and amplitudes of thunder events are generally consistent with those from the conventional lightning detection data set and broadband magnetic field. There is likely a correlation between the maximum strength of thunder events and the highest peak current for individual CG flashes. Moreover, the comparison with weather radar observations indicates that lightning usually originated from areas of high reflectivity (e.g., ≥50dBz) with diffusely distributed (from ground surface to ∼5 km altitude) thunder events and extended in a narrow altitude range of 3–5 km to areas with low radar reflectivity.

Electrical Properties of Vertical Dominant Charge Structures Observed in Corsican Thunderstorms With a LMA

AbstractLightning characteristics of Corsican storms with different charge structures are investigated in this study. Observations of an LMA network are used to document the total lightning activity. Complementary lightning observations of the lightning detection network Météorage are also used. A clustering algorithm is used to build a database of electrical cells from June to October 2018. A method is also applied to infer the vertical charge structure, as dominant dipoles, per 10‐min period for each electrical cell. As an example, one cell recorded in July 2018 is discussed. The cell database is then presented as well as the main electrical properties according to the dominant charge structures. For instance, the higher in altitude the dominant dipole, the higher the flash rate. Overall, dominant negative dipole are observed for 25% of the 10‐min periods and can be separated into two categories: (a) low altitude negative dipole class dominated by negative cloud‐to‐ground (CG) flashes with a main positive layer located between 2 and 4 km height and (b) high altitude negative dipole class, dominated by negative intracloud (IC) with a main positive layer at 5 km height. Dominant positive dipole can also be separated into two categories with (a) a dominant positive dipole located between 4.5 and 10 km high, −CG dominance, weak flash rate and (b) higher altitude dominant positive dipole, +IC dominance and a larger +CG fraction. The synergistic use of LMA and Météorage observations independently gives a rational type and polarity classification with regard to the vertical charge structure.

Joint distribution of lightning and wind speed and its impact on the transmission system

Distribution of lightning plays a key role in lightning protection of transmission lines. The design of wind deviation under lightning overvoltage of the transmission tower is an empirical parameter in the transmission line design code. Meteorological data from Hong Kong are investigated to analyze the joint distribution of lightning and wind speed. It has been uncovered that daily cloud-to-ground (CG) flashes follow the Burr distribution, which is highly skewed toward a few days with notable lightning. The lightning and wind follow a Gumbel copula joint distribution. According to empirical and theoretical distribution, there are 239 days with more than 1,000 CG lightning flashes per day, and approximately 20% of these days have a maximum wind speed of approximately or over 15 m/s. In 5 days with the number of CG lightning flashes above 30,000, 3 days have a maximum wind speed of over 15 m/s and the other 2 days have a maximum wind speed of over 10 m/s, which suggests that the severe convection with a squall line contributes much to the likelihood of the days with high wind and lightning storms.