Enhancement of Cloud-to-Ground Lightning Activity Caused by the Urban Effect: A Case Study in the Beijing Metropolitan Area

To investigate the possible impact of urban development on lightning activity, an eight-year (2010–2017) cloud-to-ground (CG) lightning dataset provided by the National-Wide Lightning Detection Network in China was analyzed to characterize the CG lightning activity in the metropolitan area of Beijing. There is a high CG flash density area over the downtown of Beijing, but different from previous studies, the downwind area of Beijing is not significantly enhanced. Compared with the upwind area, the CG flash density in the downtown area was enhanced by about 50%. Negative CG flashes mainly occurred in the downtown and industrial area, while positive CG flashes were distributed evenly. The percentage of positive CG flashes with Ipeak ≥ 75 kA is more than six times that of the corresponding negative CG flashes in the Beijing area. The enhancement of lightning activity varies with season and time. About 98% of CG flashes occurred from May to September, and the peak of CG diurnal variation is from 1900 to 2100 local time. Based on the analysis of thunderstorm types in Beijing, it is considered that the abnormal lightning activity is mainly responsible for an enhancement of the discharge number in frontal systems rather than the increase of the number of local thunderstorms. In addition, there is a non-linear relationship between pollutant concentrations and CG flash number, which indicates that there are other critical factors affecting the production of lightning.


Introduction
Lightning activity is a common cause of power system failures, property damage, and weather-related deaths. Urbanization can influence the long-term behavior of regional cloud-to-ground (CG) lightning activity [1]. Since the 1980s, many studies have been conducted to investigate the impact of urban development on CG lightning activity in different geographical areas of the world [2][3][4][5][6][7][8][9][10][11][12]. Westcott [2] firstly proposed that the urban effect can enhance the lightning activity over and downwind of cities based on data from the National Lightning Detection Network in the United States. Steiger et al. [3] used data from the same lightning detection network and also reported that a significant enhancement of 45% in the CG flash density was observed over Houston, Texas compared to the nearby surrounding areas, and it was shown that the percentage of positive CG flashes exhibited a significant decrease (12%).
The urban effect is primarily a combination of the thermodynamic effect caused by the urban heat island (UHI) [13,14] and the increase of pollution concentration due to the plentiful anthropogenic activity [15,16]. The heat from downtown can modify local convective processes associated with thunderstorms [6], while pollutants in the boundary layer could uplift the concentration of cloud condensation nuclei and thus suppress the mean droplet size [5,[17][18][19]. The enhancement impact of UHI on the lightning activity (thunderstorm lifetime and thunderstorm initiation) over cities has been relatively well understood [12], but there is no consensus on the role of pollutants on this enhancement. Since air pollution is influenced by different meteorological processes, it is not easy to infer the specific role of pollutants on the increase of CG lightning activity.
Based on CG lightning data from the lightning detection network deployed on the Iberian Peninsula that belongs to the Spanish Meteorological Institute, Soriano and de Pablo [4] analyzed the lightning activity in several small Spanish towns using the data of three years, and the results showed that the higher SO 2 (sulfur dioxide) concentration contributes to the increase in the number of CG flashes, while the PM 10 contribution (i.e., particulate matter with aerodynamic diameter <10 µm) seems to be irrelevant. However, Naccarato et al. [6] showed that the PM 10 concentration is positively correlated with the number of CG flashes over three large metropolitan areas in southeastern Brazil. It is likely because the air pollutants saturate the concentration of aerosol particles at a specific level, inhibiting the intensification of lightning activity [10,20]. Farias et al. [11] found that PM 10 increases the lifetime of thunderstorms and the number of flashes per thunderstorm rather than the flash rate of thunderstorms, combined with pollutant data provided by the São Paulo environment agency. Kar et al. [8] and Kar and Liou [9] also indicated that the higher concentrations of PM 10 and SO 2 contribute to the enhancement of CG lightning activity. Although considerable studies have been conducted on the influence of pollutants on the CG lightning activity, the present understanding is incomplete, and extensive studies are still needed.
With the rapid economic and social development, anthropogenic pollutant emissions and urban sprawl have had an increasing impact on lightning activity. With the increase of the world's population, the number of megacities is also growing rapidly. Most existing studies regarding the urban effect on lightning activity are focused on cities with a population of less than 10 million, while the characteristics of lightning activity in megacities are still rarely reported. Beijing, as the capital of China, has a population of about 21 million as well as an urban area of about 16,400 km 2 , which is much larger than that of previous studies [2,4,5,7,9]. In this paper, we attempt to extend the knowledge of the urban effect on lightning activity over megacities by investigating the characteristics of CG flashes in Beijing. Moreover, the comprehensive understanding of lightning activity with multi-year observation data is also the basis of effective prediction and warning of thunderstorm-related disasters in Beijing City.

Data and Methods
At present, the studies on lightning characteristics in Beijing are still limited. This study used an eight-year (2010-2017) CG lightning dataset from the National-Wide Lightning Detection Network operated by the State Grid Electric Power Research Institute to investigate the characteristics of lightning activity in the metropolitan area of Beijing [21]. The dataset contains information such as time, longitude, latitude, peak current, and polarity of CG strokes. The detection network combines the time-of-arrival method and the magnetic detection method to determine the geographical location of CG flashes [21]. Chen et al. [21] showed that the detection efficiency of CG flashes and CG strokes is about 94% and 60%, respectively; the arithmetic mean and the median of location errors are about 710 m and 489 m, respectively; the percentage of evaluation errors of lightning peak current is between 0.4% and 42%; the arithmetic mean and the median is about 16.3% and 19.1%, respectively. Since the detection principle, parameter setting, locating algorithm, and hardware remain unchanged, it could be used as the reference in this paper.
Previous studies indicated that this type of lightning detection network is usually contaminated by intra-cloud (IC) discharges [22][23][24]. Therefore, this work eliminated the misidentification of IC discharges by removing the detection results with peak current <10 kA and assuming that the rest are all CG strokes. Although the contamination level of IC discharges could still exceed this threshold, the error is almost negligible, especially for positive CG strokes that typically have high peak currents. It should be mentioned that the detection network underwent several configuration changes and performance improvements in 2009 and 2018, respectively. In order to exclude the evaluation error caused by the variation in the detection efficiency, this paper only selected the lightning detection data from 2010 to 2017.
In addition to peak current, there is another metric that has been used widely to evaluate the capacity of charge transfer by a particular lightning stroke, namely the impulse charge moment change (iCMC, defined as the product of charge transferred to ground within 2 ms after the return stroke and the original height of deposited charge in the thundercloud) [25]. Based on a three-year measurement in North America, Cummer et al. [26] provided an analytical fit for the mean iCMC as functions of peak current, by differentiating positive and negative CG strokes, as follows: For negative CG strokes with peak current in the range of −10 to −200 kA, the analytical fit to the mean iCMC is For positive CG strokes with peak current in the range of +40 to +200 kA, the analytical fit to the mean iCMC is For positive CG strokes with peak current in the range of +10 to +40 kA, the analytical fit to the mean iCMC is iCMC − (C·km) = 0.53 I pk (kA) − 0.00086I 2 pk (kA).
As shown in Figure 1a, Beijing is a city sited on the northwestern margin of the North China Plain, backed by the remaining Taihang Mountains and the Yanshan Mountains, and about 150 km from the Bohai Sea in the southeast. It is the world's most populous capital city, with over 21 million residents within an administrative area of 16,410 km 2 (http://tjj. beijing.gov.cn/EnglishSite/ accessed on 21 March 2021). Beijing is mostly surrounded by Hebei Province with the exception of neighboring Tianjin to the southeast. The topography of Beijing is generally high in the northwest and low in the southeast. The city's topography consists of two major geomorphic units: the northwestern mountains and the southeastern plains, of which 10,072 km 2 (61.4%) is mountainous and 6338 km 2 (38.6%) is plain. The average altitude of Beijing is 43.5 m above sea level, with the altitude of the plains ranging from 20 to 60 m, while the mountains are generally 1000 to 1500 m above sea level. Since the North China Plain has less topographic relief, we only plot the topographic map of mountain ranges (i.e., at altitudes >500 m above sea level). The city is about 180 km long from north to south and 160 km wide from east to west (see the inset in Figure 1a). Beijing is located in a warm-temperate semi-humid region, and its climate is a monsoon-influenced humid continental climate, with an average annual precipitation of about 600 mm in the plains. Sixty percent of its precipitation occurring in the summer months of July and August, while the air is dry during the rest of the year. The average annual temperature is about 13 • C. The monthly average temperature in January is about −3 • C, while July has an average temperature of about 27 • C. Thunderstorms in the research area are most active in summer, followed by spring and fall, while there are very few thunderstorms reported in winter. Due to the influence of the East Asian monsoon and the vast Siberian

Analysis and Results
The characteristics of lightning activity considerably depend on the spatial scale of thunderstorms, and thus, we classified the lightning-producing thunderstorms as local thunderstorms (that generally stay in the study area) and frontal thunderstorms (that pass through the study area). During the eight-year time period examined in this paper, a  Figure 1b shows the spatial distribution of the land-surface temperature. The results showed that the UHI effect is obvious in the downtown area of Beijing. With the increase of altitude and latitude, the land-surface temperature decreases. As the result of a correlation between lightning activity and deep convection [27][28][29], we also examined precipitation data (spatial resolution: 0.1 • × 0.1 • ) from the Global Precipitation Measurement (GPM) (https://gpm.nasa.gov/GPM/ accessed on 21 March 2021) and balloon sounding data from the Integrated Global Radiosonde Archive (https://www.ncdc.noaa.gov/ accessed on 21 March 2021) to understand the atmospheric background status. The precipitation rate gradually increased from west to east, especially near the coastline (not shown). Precipitation averages around 570 mm annually, with close to three-quarters of that total falling from June to August. As a result of the strong convective activity in thunderstorms, in order to characterize convective weather systems with sufficient energy, we use only the available potential energy values (CAPE) of convection above the marginally unstable (>1000 J/kg). The annual mean CAPE varies from 1860 to 2820 J/kg, allowing updrafts to bring more pollutants into the cloud [28,29]. Moreover, we also analyzed the summertime wind field of the European Centre for Medium-Range Weather Forecasts (https://www.ecmwf.int/ accessed on 21 March 2021), and the results showed that the northeast of Beijing is the downwind area. The air pollution emissions of automotive vehicles and industrial activities mainly include sulfur oxides, nitrogen oxides, and particulate matter [30]. Therefore, we examined the pollutant concentration data from the Beijing Municipal Environmental Monitoring Center (http://www.bjmemc.com.cn/ accessed on 21 March 2021) to determine a possible correlation between air pollutants (including PM 10 , PM 2.5 , SO 2 , and NO 2 ) and lightning activity. With these analyses, it is desired to more comprehensively understand the urban effect on the CG lightning activity in Beijing City.

Analysis and Results
The characteristics of lightning activity considerably depend on the spatial scale of thunderstorms, and thus, we classified the lightning-producing thunderstorms as local thunderstorms (that generally stay in the study area) and frontal thunderstorms (that pass through the study area). During the eight-year time period examined in this paper, a total of 632 thunderstorms were recorded in the research area, including 63 local thunderstorms and 569 frontal thunderstorms. Compared with the frontal thunderstorms, the local thunderstorms are not only less in number but also shorter in duration and smaller in scale. This indicated that the enhancement of the CG lightning density in Beijing may be mainly caused by the frontal thunderstorms associated with the mid-latitude cyclones, while the isolated thunderstorms produced by the local convection had little effect. More detailed information is shown in Table 1. The annual mean CG flash density centered in the metropolitan area of Beijing is presented in Figure 2a. Over 310,000 CG flashes were recorded in this study region from 2010 to 2017. The CG flash density ranges from less than 1 flash km −2 yr −1 in the northwestern of Beijing to over 4 flashes km −2 yr −1 in the east. The abnormally highvalue area of eastern Beijing may be due to pollutants from industrial sources. Compared with the upwind area, the CG flash density of the downtown area increased by about 50%. It is noteworthy that the results of this paper are not in good agreement with those in the scientific literature [4,9]. Previous studies have shown that CG flash density is Remote Sens. 2021, 13, 1228 6 of 13 significantly enhanced both over and downwind of cities compared with other neighboring areas. However, the density-enhanced area of CG lightning was mainly concentrated in the downtown area of Beijing, and there was no obvious enhancement in the downwind area. This may be because the urbanization development of Beijing is far more than those reported cities. The UHI effect caused by the large urban sprawl has slowed down the local convection effects. Thus, it is reasonable to believe that the enhancement of urban effects on lightning activity can be weakened when the urban circulation transports pollutants over longer distances. Moreover, Tianjin is also a municipality, with a total population estimated at 15,600,000. It is one of the nine national central cities in Mainland China. It has also an obvious UHI in the urban area (Figure 1b), but the intensification of lightning activity is not significant (Figure 2a). This indicates that the UHI effect may not be the only factor in the context of high pollution.  Figure 3 shows the geographical distribution of the positive and negative CG flash density in the research region, with significant differences. The negative CG flashes mainly occurred in the downtown and northeastern of the city, whereas the distribution of the positive CG flashes was relatively scattered. The cloud microphysical processes that generate electric charge cannot be analyzed due to the lack of corresponding dual-polarization radar data. Existing results show that the charge reversal temperature will shift to higher temperatures when the cloud droplets contain available air pollutants [3].  Figure 2b shows that the diurnal variation in the CG lightning number has a significant periodic tendency. Lightning activity increased from 0900 local time (LT), reached a maximum at about 2000 LT, and then slowly decreased. Previous studies have shown that the peak of lightning activity is between 1500 and 1700 LT in the continental area [31][32][33]. The significant increase of lightning activity in the evening is mainly because the surface temperature rises with the increase of solar radiation in the daytime. In the evening, the surface temperature and water vapor in the near-earth layer can reach high levels, which leads to convective instability in the atmosphere and further triggers thunderstorms. The reason why the peak value of lightning activity moved back in Beijing may be due to traffic pollution during the evening rush hour. The presence of high aerosol concentrations can delay the onset of precipitation and lightning activity [34]. Figure 2c shows the monthly variation of the CG flash number with a single-peak feature. About 98% of all detected CG lightning flashes occurred between May and September. The maximum and minimum of the CG flash number occurred in July and February, respectively, which also illustrates the importance of solar radiation in forcing the deep convection and forming the unstable layer. The monthly variation is similar to that reported by Holle et al. [35] who indicated that in the United States, the maximum CG flash number occurred in July, but it is different from that reported by Pinto et al. [31] who found two peaks in February and November. Pinto et al. [31] showed that the air temperature of Southeastern Brazil also has a double peak feature typical of the tropics, which confirms the tendency verified by Pinto et al. [36] of the monthly variation of the number of flashes to follow the variation of the air temperature in the tropics. North China and the Midwest of the United States are both in the middle latitudes, where only one temperature peak occurs in the summer [37,38]. Thus, differences of the temporal distribution may be related to the geographical location and climate background of different cities.
Since there is no available extremely-low-frequency (ELF, ≤1 kHz) data in the Beijing area to evaluate the iCMC of CG strokes [39], this paper used empirical formulas presented by Cummer et al. [26] to perform a rough evaluation (see Equations (1)-(3) in Section 2). We qualitatively estimate the charge transfer of CG strokes by adding the corresponding iCMC together. The amount of charge transfer in the thundercloud is better correlated with the lightning number, rather than the peak current, which shows that the total transferred charge is only related to the number of CG flashes. Notably, the iCMC appears to be inversely proportional to the peak current, and the charge transferring ability of lightning is stronger in winter compared to summer (see Figure 2b,c). This is likely because each discharge will "disturb" the charge distribution in the thundercloud, thus inhibiting the discharge process and causing more charges to accumulate in the charge reservoir. Figure 3 shows the geographical distribution of the positive and negative CG flash density in the research region, with significant differences. The negative CG flashes mainly occurred in the downtown and northeastern of the city, whereas the distribution of the positive CG flashes was relatively scattered. The cloud microphysical processes that generate electric charge cannot be analyzed due to the lack of corresponding dualpolarization radar data. Existing results show that the charge reversal temperature will shift to higher temperatures when the cloud droplets contain available air pollutants [3]. Therefore, the main negative charge region can extend to the base of the cloud and cover the positive charge center below, thereby generating more negative CG flashes. For smallsize graupel/ice crystals, they are better able to migrate from urban areas to suburban areas through processes of transport and dispersion of atmospheric flow. Since small-size graupel/ice crystals are positively charged between −10 and −25 • C [40,41], the positive charge region is more widely distributed and generates more positive CG lightning in suburbs. Similar lightning enhancements have been also observed in the surrounding areas of other cities, such as Baoding, Tangshan, and Cangzhou (not shown).
urban areas through processes of transport and dispersion of atmospheric flow. Since small-size graupel/ice crystals are positively charged between −10 and −25 °C [40,41], the positive charge region is more widely distributed and generates more positive CG lightning in suburbs. Similar lightning enhancements have been also observed in the surrounding areas of other cities, such as Baoding, Tangshan, and Cangzhou (not shown).    (about 100,000), and the proportion of negative CG flashes increased in summer. The diurnal and monthly variation tendency of positive and negation CG flashes is similar to that mentioned above, which will not be described here. It should be noted that the number of positive CG flashes tends to increase/decrease more slowly compared with negative CG flashes, and the mean peak current of positive CG flashes is larger than that of negative CG flashes. This paper defined the high peak current threshold as 75 kA according to the method given by Lyons et al. [42]. As illustrated in Figure 4c, the percentage of negative CG flashes with I peak ≥ 75 kA is 4.18%, while the percentage of positive CG flashes is more than 6 times (26.61%) that of the former. This may Remote Sens. 2021, 13, 1228 9 of 13 be related to the different discharge mechanisms of positive and negative CG flashes [43]. In addition, we also examined the spatial distribution of CG flashes with different peak currents. The distributions of strong CG flashes have distinct regional clustering, and they are more concentrated in downtown and industrial areas.
positive and negation CG flashes is similar to that mentioned above, which will not be described here. It should be noted that the number of positive CG flashes tends to increase/decrease more slowly compared with negative CG flashes, and the mean peak current of positive CG flashes is larger than that of negative CG flashes. This paper defined the high peak current threshold as 75 kA according to the method given by Lyons et al. [42]. As illustrated in Figure 4c, the percentage of negative CG flashes with Ipeak ≥ 75 kA is 4.18%, while the percentage of positive CG flashes is more than 6 times (26.61%) that of the former. This may be related to the different discharge mechanisms of positive and negative CG flashes [43]. In addition, we also examined the spatial distribution of CG flashes with different peak currents. The distributions of strong CG flashes have distinct regional clustering, and they are more concentrated in downtown and industrial areas.

Discussion
It can be seen from Figures 1b and 2a that although the UHI is an important factor in increasing lightning activity, it does not seem to be the only factor in the metropolitan area of Beijing. CG lightning activity is not only enhanced in land-surface temperature areas caused by artificial structures and the lack of vegetation, but also significantly enhanced in areas where human activity is more intense (e.g., traffic and industrial pollutant emissions. Most previous studies showed a positive correlation between pollutant concentrations and lightning activity [6,8,9,11], but a negative correlation when pollutant concentrations are high [2,10]. Although there are several different views on the details of lightning initiation [44,45], most scholars agree on the basic principle of the non-inductive charge separation process [46,47]. The augment of pollutants could affect the cloud electrification by altering the concentration, phase, size, and surface property of graupel/ice crystals in the mixed-phase region (from 0 to −40 • C), resulting in the generation of more CG lightning [41,48].
To ascertain the possible effect of anthropogenic pollutants, we examined the correlation between the CG flash number of 222 thunderstorms passing through the downtown and average concentrations of corresponding pollutants (i.e., SO 2 , NO 2 , PM 2.5 , and PM 10 ). However, no significant correlation was found in scatter plots (Figure 5b-e). We also examined the temporal relation of the CG flash number and pollutant concentrations in the vicinity of air pollution monitoring stations (red rectangle in Figure 5a). Remote Sens. 2021, 13, x FOR PEER REVIEW 12 of 15 relevant to the presentation of this paper, we will make an in-depth analysis of the relevant content in the follow-up work.

Conclusions
Even though the local enhancement of lightning activity in the urban areas has been studied extensively, the urban effect of megacities is still lacking. In this work, an eight-year CG lightning dataset was investigated to study the long-term behavior of lightning activity in the metropolitan area of Beijing. The dataset was analyzed from both spatial and temporal perspectives, and the paper also examined the correlation between lightning activity and other parameters such as land-surface temperature and pollutants. The high CG lightning density area was mainly concentrated in the downtown of Beijing, but there was no remarkable enhancement in the downwind area. The CG flash density in the downtown was about 50% higher than that of the upwind area. About 98% of CG flashes mainly occurred from May to September and the enhancement of CG diurnal variation was most pronounced from 1900 to 2100 LT. Negative CG lightning mainly occurred in the downtown and industrial area, while the distribution of the positive CG flashes was relatively scattered. The percentage of positive CG lightning with Ipeak ≥75 kA is more than six times that of the negative CG flashes. By comparing the spatial relationship between the UHI and the CG lightning density in surrounding cities and industrial areas, the enhancement of lightning activity is influenced by multiple factors. Therefore, longer periods and more precise observational studies from different regions of the globe are required.
Author Contributions: Y.W. (Yongping Wang) was responsible for data curation and management, and wrote the whole of this paper. Y.W. (Yu Wang) was a supervisor of this research project and provided funding, data sources, methodology of data curation. G.L. was a principal investigator and a supervisor of this research project, and provided funding. T.S. and C.P. conducted lit- This paper takes the summer season of 2014 as an example. The peak value of the CG flash number usually corresponds to the peak value of pollutant concentrations in a period of time, and it has a slight delay, especially for particulate matter (Figure 5f-i). However, the corresponding relationship is not absolute, and the occurrence of lightning sometimes does not correlate with the concentration of pollutants. This is because the pollutant concentration varies with different meteorological conditions, so it is difficult to quantify the contribution of air pollution to the enhancement of overall lightning activity. Soriano and de Pablo [4] have also indicated that PM 10 is a non-influential parameter for enhancing urban lightning activity. In addition, the aerosol could reduce the electrical activity by absorbing solar radiation [18,49,50], which also suppresses the deep convection and the relative humidity [7,51]. In summary, the lightning activity is influenced by a variety of factors (e.g., urban heat island circulation, frictional lift, aerosol increase, cloud water content, etc.), and further research is needed in the future work to clarify the physical mechanisms that lead to lightning enhancement.
Since the air pollution monitoring stations are mostly deployed in urban areas and only a few are installed in agricultural and industrial areas (see Figure 5a), they can only represent the situation in a limited area. This is very mismatched with the wide spatial distribution of lightning activity. The suburbs of Beijing have dense infrastructures such as refineries, chemical facilities, steelworks, and power plants, but the overall population is low and the urban development is slight. Thus, we also examined the ultraviolet (270-380 nm) aerosol index (UVAI) data (spatial resolution: 13 km × 12 km) from the Ozone Monitoring Instrument (OMI) aboard the Aura satellite (http://projects.knmi.nl/omi/ accessed on 21 March 2021) to verify the aerosol optical depth in non-urban areas [52,53]. However, unfortunately, the scanning range of the satellite orbit is limited, so the relevant data corresponding to the study area cannot be obtained. We examined the correlation between the concentration of pollutants and the lightning activity, and the results showed there was a good spatial correspondence between the high lightning density and the high aerosol concentration in Baoding (not shown). Industrial pollution may play a greater role than domestic pollution. As it is not relevant to the presentation of this paper, we will make an in-depth analysis of the relevant content in the follow-up work.

Conclusions
Even though the local enhancement of lightning activity in the urban areas has been studied extensively, the urban effect of megacities is still lacking. In this work, an eight-year CG lightning dataset was investigated to study the long-term behavior of lightning activity in the metropolitan area of Beijing. The dataset was analyzed from both spatial and temporal perspectives, and the paper also examined the correlation between lightning activity and other parameters such as land-surface temperature and pollutants. The high CG lightning density area was mainly concentrated in the downtown of Beijing, but there was no remarkable enhancement in the downwind area. The CG flash density in the downtown was about 50% higher than that of the upwind area. About 98% of CG flashes mainly occurred from May to September and the enhancement of CG diurnal variation was most pronounced from 1900 to 2100 LT. Negative CG lightning mainly occurred in the downtown and industrial area, while the distribution of the positive CG flashes was relatively scattered. The percentage of positive CG lightning with I peak ≥ 75 kA is more than six times that of the negative CG flashes. By comparing the spatial relationship between the UHI and the CG lightning density in surrounding cities and industrial areas, the enhancement of lightning activity is influenced by multiple factors. Therefore, longer periods and more precise observational studies from different regions of the globe are required.
Author Contributions: Y.W. (Yongping Wang) was responsible for data curation and management, and wrote the whole of this paper. Y.W. (Yu Wang) was a supervisor of this research project and provided funding, data sources, methodology of data curation. G.L. was a principal investigator and a supervisor of this research project, and provided funding. T.S. and C.P. conducted literature research on urban lightning activity. M.M., revised the manuscript. B.Z. and D.L. provided meteorological data. All authors have read and agreed to the published version of the manuscript.