1. Introduction
Presently, much attention is paid to the study of the regional features of mesoscale convective systems and trends of lightning activity. This is due to an increase in the frequency of atmospheric hazards in a changing climate [
1,
2]. Severe convective storms are destructive weather events that cause substantial damage and fatalities, which are commonly observed in different parts of the world [
3,
4].
Thunderstorms themselves influence the climate system via the redistribution of heat, moisture, and momentum in the atmosphere [
5]. The formation of thunderclouds is accompanied by electrification processes, which give rise to lightning discharges. Lightning is a major hazard to the public, nature, aviation, and power companies that causes great damage and death.
The prediction of thunderstorms and their characteristics remains one of the most difficult problems in modern weather forecasting. The need to develop short-term forecasting systems in the territory of Russia is explained by an increase in convective precipitation, which has been observed in Russian regions over the last 50 years [
6]. Moreover, Chernokulsky et al. [
7] predict that the risk of severe convective events will continue to rise in the 21st century.
Various instability indices are used to forecast convection-driven hazardous weather events and thunderstorm activity [
8,
9,
10,
11]. There are also methods for identifying thunderclouds using meteorological radar products, which are based on counting cells with certain parameters (primarily radar reflectivity) as indicators of thunderstorms (see, for example, [
12,
13,
14,
15,
16]). Such indicators, as a rule, are quite sensitive to regional climatic and geographical features, as well as to the equipment used; therefore, their use requires adjustment for a specific region and equipment. One of the most common methods for predicting local atmospheric phenomena is numerical mesoscale modeling (see, for example, [
17,
18,
19]). Analyzing the microphysical and electric parameters of thunderclouds obtained through numerical simulations and field observations is a key factor in improving the prediction capabilities of a system with strong convection.
The presence of aerosol particles affects the microphysical and electrical parameters of convective systems. As cloud condensation nuclei lift up to the mixed-phase region, more liquid water is formed and, hence, more ice-phase particles are produced, which impacts charge separation and lightning activity (see, for example, [
17,
20,
21]. Based on the Weather Research and Forecasting (WRF) model simulations, Khain et al. [
22] demonstrated that accumulated aerosols acting as cloud condensation nuclei lead to an increase in the number of cloud droplets and the release of latent heat. The study of aerosol effects, both in individual convective events and climate trends, has recently received increased attention (see, for example, [
23,
24,
25,
26]). An increase in the growth rate of snow and graupel and, as a result, an intensification of electric processes due to an increased aerosol load were shown in [
19]. In addition, Zhao et al. [
27] pointed to the correlation between the concentration of sulfates and lightning activity in the plateau and basin regions of Sichuan (China). Some observational studies demonstrated the impact of increased cloud condensation nuclei on lightning activity (see, for example, [
28,
29]).
Sun et al. [
30] noted a need to carry out numerical experiments aimed at studying the effect of aerosols acting as ice nuclei on lightning activity. Zheng et al. [
31] pointed out that the characteristics of lightning activity in severe thunderstorms are different from those in regular thunderstorms. For this reason, a detailed study of the characteristics of individual powerful thunderstorm events is of great interest.
In this paper, we analyze the number and characteristics of the thunderstorms during the convective season of 2020 in the Upper Volga region, which encompasses the majority of the Nizhny Novgorod region, using data from the World-Wide Lightning Location Network (WWLLN) and electric field mill (EFM) measurements; identify the most powerful convective event; and study its characteristics using numerical modeling, meteorological radar observations, WWLLN data, and electric field mill measurements. Thus, the purpose of this work is to find general trends in the lightning activity and particular features of thunderstorm formation in the Upper Volga region.
2. Materials and Methods
The primary focus of this work is the study of thunderstorms observed during the 2020 convective season over the Upper Volga region (with the center in Nizhny Novgorod) and the most powerful convective event of this season. The analysis of the convective season as a whole focuses on identifying the powerful thunderstorm events and determining their characteristics and general trends. The analysis of the most powerful thunderstorm is based on the use of a wide range of experimental and numerical tools to study the characteristics and features of the convective system’s development.
2.1. Weather Archive Data
The general meteorological situation in the area under consideration was analyzed by the open-access data from the <<rp5>> weather archive [
32] at the meteorological station <<Nizhny Novgorod>> (56.32
N, 44.00
E) from April to September 2020. The following parameters were analyzed (all of them were measured every 3 h): air temperature and relative humidity at a height of 2 m above the ground; wind direction and speed at a height of 10–12 m, averaged over a 10-min period; current meteorological events and weather in the previous 3 h; total cloudiness; height of the cloud base; and precipitation rate. The performed analysis made it possible to identify thunderstorm events during the convective season and also to determine their approximate duration and the accompanying meteorological conditions.
2.2. Meteorological Radar
We used data from the Doppler meteorological radar <<DMRL-C>> located on the roof of the <<Oka>> hotel 56.29
N, 43.98
E (see
Figure 1, right panel). The radar’s area of view is 11 km in height and 250 km in radius from the observation point. The DMRL-C scans the atmosphere with a time resolution of 10 min and provides primary measurements such as the radar reflectivity of the vertical and horizontal polarizations, the linear depolarization ratio, the radial velocity, the width of the Doppler spectrum of radial velocity, and other measurements typical of radars. At a peak radiated power of 15 kW, the maximum dynamic range of each polarization was not less than 100 dB [
33].
In this study, we used only the radar reflectivity Z measured in decibels (dBZ). The radar reflectivity was determined on a 3D grid with a vertical step of 1 km and a horizontal step of 4 km, i.e., each grid cell has a size of km.
2.3. Electric Field Measurements
Multipoint long-term measurements of the atmospheric electric potential gradient have been conducted in Nizhny Novgorod since 2012.
Figure 1 (right panel) shows a layout of the measurement points in Nizhny Novgorod: the electric field mill Boltek EFM-100 [
34] installed on the roof of the IAP RAS 56
19
25
N, 44
01
21
E, and the meteorological radar installed on the roof of the <<Oka>> hotel for electromagnetic field reception and recording in the territory of Nizhny Novgorod. The Boltek EFM-100 is used for recording the electric field and has an exact reference to UTC time, which enables joint analysis with other time series, for example, themeasurements of global lightning detection networks. The rated discretization frequency of measurements and the sensitivity of the EFM amount to 20 Hz and ±80 kV/m, respectively.
The efficiency of the tool is illustrated in
Figure 2 using the lightning activity in the Nizhny Novgorod region on 13–14 July 2020 as an example. All lightning events detected by the WWLLN within a radius of 100 km from the observation point are shown in the left panel by the blue color. The discharges that coincided in time with the discharges determined according to the EFM measurements are marked by the red color. The ratio of the number of lightning events registered by both the EFM and WWLLN to the number of all lightning events detected by the WWLLN is demonstrated in the right panel.
Figure 2 confirms the EFM documentation [
34], which states that the maximum radius for detecting lightning is 38 km.
The continuous atmospheric electric potential gradient records are interesting from the viewpoint of thunderstorm climatology. The comparative analysis of the thunderstorm electricity observation results and the meteorological data are also of great significance. Electrostatic fluxmeter data records allowed us to evaluate the number of discharges by a fairly simple calculation method using the specified threshold of the leap of the filtered electric field strength derivative [
35]. First, the signal was filtered to extract the high-frequency component, and then the discharge itself was separated. A set of the following conditions was taken as the discharge criterion: (1) more than 0.8 s elapsed from the last detected discharge, (2) the time-dependence sharpness of the field intensity did not exceed 5 kV/(m·s) for 0.1 s before the assumed discharge time, and (3) the field intensity varied by more than 1.5 kV/m for 0.25 s from the assumed discharge time. The criteria were selected taking into account local characteristics using the long-term observations potential gradient database and verified using DMRL-C records and lightning-direction-finding data.
2.4. Lightning Data
We use observational lightning data from the World-Wide Lightning Location Network [
36]. This network provides information about thunderstorm activity over the entire globe by registering the electromagnetic signals from thunderstorms in the VLF range (3–30 kHz). Currently, it consists of more than 60 detectors, including those in Russia. According to various studies, the network provides an opportunity to localize lightning discharges on a global scale with a high spatial resolution and makes it possible to obtain estimates of the spatial and temporal variability of lightning activity for any region of the planet. The lightning data contain the location, time of discharge, time uncertainty, energy of the strike, energy uncertainty, and the number of stations that have recorded the discharge.
For our research, the WWLLN data for the internal domain of the simulations are analyzed. The number of lightning strikes is summed up over km squares every 3 h. To analyze lightning activity during the convective season, the distribution of the number of discharges per day in the area under consideration is plotted vs. the time. To analyze the lightning activity of an individual convective event, we consider the time dependence of the number of discharges summed over 10 min.
2.5. Numerical Simulation
The numerical modeling is based on the Weather Research and Forecasting (WRF) model [
37,
38]. It is a collaborative project designed for the development of a mesoscale weather forecasting model to gain better insight into mesoscale weather-related processes. The WRF model includes many parametrizations, including the parametrizations of cloud physics, radiation, microphysical processes, the planetary boundary layer, and the underlying surface. Thus, the WRF model is suitable for numerical weather prediction and atmospheric modeling in various research and practical tasks from small-scale to global modeling.
The numerical simulation was performed on two nested domains: the 1890 km × 1890 km external domain with a horizontal grid spacing of 3 km and the 210 km × 210 km internal domain with a horizontal grid spacing of 1 km. The center of the computational area is chosen to coincide with the location of the Nizhny Novgorod aerological station (56.29 N, 43.98 E), equipped with a meteorological radar and involved in monitoring the state of the atmosphere. The internal domain in this work is used as the area under consideration unless otherwise noted.
The computational domains used to simulate the thunderstorms in the Nizhny Novgorod region are shown in
Figure 1 (left panel). The area outlined in green corresponds to the external domain, and the area outlined in blue to the internal domain. The vertical grid contains 40 levels from 0 to 20 km. The calculations were made for 36 h, with the start time such that the thunderstorms did not fall into the model’s adjustment period to the initial conditions. For each convective event, the start time of the calculation was set as a standard meteorological hour that is at least 6 h before the beginning of the thunderstorm.
Two types of Thompson’s microphysical parametrizations were used in the numerical simulations by the WRF-ARW model for studying the impact of aerosols on thundercloud development. The first of them [
39] described only solid hydrometeors, which include graupel, snow, and ice. The second parametrization [
24] also took into account the presence of aerosol particles in the atmosphere. In this parametrization, the initial distributions of the aerosol mass mixing ratios are taken from multi-year global model simulations by the Goddard Chemistry Aerosol Radiation and Transport (GOCART) model [
40]. We use background mass mixing ratios from the GOCART model as the initial distribution of aerosols in our research.
Since a rather weak instantaneous correlation between the simulation and the observational data is obvious, in this paper, we propose visualizing the characteristics of a thunderstorm using a function that best resembles the trace of a cloud (such as the trace that a tornado leaves behind). This function corresponds to the 2D horizontal distribution of the maximum value of the parameter in the vertical column for the selected time interval. In this paper, we use 3 h observation intervals with the initial times at 0, 3, 6, 9, 12, 15, 18, and 21 UTC.
2.6. Electric Parameters Calculation
The electric potential and electric field are the parameters of thunderclouds that directly affect the occurrence of electric discharges and lightning activity. For a detailed analysis of the features of a convective event’s development, the electric parameters of the clouds are calculated. The parametrization of the electric processes described in [
41] is used in this work to calculate the electric parameters of thunderstorms.
This parametrization is based on the assumption of non-inductive charging of solid hydrometeors. The main charge-transfer particles are small ice particles and graupel pellets [
42]. The electric charge accumulated on snow particles is much less than on ice crystals or graupel pellets [
43]; therefore, it is not considered as a charge carrier in this work. The reversal temperature is assumed to be equal to −15
C. The ice particles have a negative charge and the graupel pellets have a positive charge up to the reverse point and vice versa above it. The spatial distribution of the electric charge at each time step is calculated from the hydrometeors mixing ratios. It is supposed that the charge density of each type of hydrometeor is proportional to the mass per unit volume of these hydrometeors:
where
is the mixing ratio of the graupel and ice particles (in g·g
),
is their maximum value,
is the maximum charge that can accumulate on a graupel pellet (in C), and
is the maximum graupel concentration (in m
). It is also assumed that charge separation occurs mainly due to collisions of solid hydrometeors in vertical air flows; therefore,
is the coefficient (in C·m
) that is found from the condition of general electric neutrality in a vertical column where charge separation occurs. The maximum charge on the graupel is taken to be equal to
C, while the maximum graupel concentration is assumed to be 200 m
[
44].
The temporal resolution of the WRF output data is 10 min, which allows us to consider the pattern of the electric processes at each time step to be a quasi-static one. To calculate the electric potential
, the three-dimensional Poisson equation is solved:
where
is the electric charge density, which is the sum of the charge densities on the graupel and ice particles, and
is the vacuum permittivity.
The electric field at each grid point is defined as the electric potential gradient with an opposite sign:
The discharge processes are not directly included in the description of the electric processes in thunderclouds. However, the growth of the electric field strength is limited by the maximum electric charge density on the particles. It should also be noted that due to the coarse spatiotemporal resolution (compared to the field measurements), the obtained distributions of the electric parameters are averaged.
4. Discussion
In this paper, we analyze the lightning activity and characteristics of thunderstorms in the Upper Volga region. It was revealed before that about 34 convective events in Nizhny Novgorod (the center of the considered region) and its surroundings occur annually. According to the data presented in
Figure 4, about 15 thunderstorms occurred in the 2020 convective season. The difference in the number is due to the presence of weak local thunderstorms with a small number of discharges, which are not displayed on the graph against the background of more powerful events and were not recorded by meteorological observation stations.
Every year one or two severe thunderstorms similar to the thunderstorm of 13–14 July 2020 studied in this article occur in the territory of the Upper Volga region [
45]. For example, a similar thunderstorm that occurred on 1–2 June 2015 was studied in [
52]. The forecast for that thunderstorm in the territory of the Nizhny Novgorod region was based on the numerical simulations of the WRF model with the same non-aerosol parametrization of the microphysics as in the present work. The same algorithm for calculation of the electric parameters was used. The data from continuous EFM measurements and weather radar data were used for verification of the simulation results. The powerful continuous thunderstorm of 1–2 June 2015 passed through the Nizhny Novgorod region from south to north. According to the EFM measurements, the thunderstorm began at about 20:00 UTC and lasted almost 8 h; the total area of the thundercloud, according to radar observations and modeling results, reached 20,000 km
. During the first stage of the thunderstorm, which was about three hours of the developed thunderstorm activity, there was a maximum correspondence between the simulation results and the real parameters of the thunderstorm. The agreement between the simulated and observed parameters worsened towards the end of the thunderstorm activity.
Thus, in both cases, we dealt with thunderstorms that were the most powerful intra-mass events in the convective season usually starting in the afternoon local time, moving into the next day, and continuing through the night. In addition, these thunderstorms usually come from the southern direction. Interestingly, when simulating both thunderstorms (1–2 June 2015 and 13–14 July 2020) using the WRF model, there was no final stage of the thunderstorm; thereby, the modeled thunderstorms ended too early. Thus, it can be assumed that such long-term thunderstorms can be classified as a separate type and similar parameters can be expected to be obtained both as a result of observations and numerical simulation.
In this work, to visualize the spatial distributions of the characteristics of a thundercloud, a 3 h trace is used, which is the distribution on the plane of the maximum value of the parameter in a vertical column for a selected time interval (3 h). The use of 3 h averages (see, for example, [
10]) greatly smooths the real picture of a thunderstorm and, from our point of view, is not correct when studying the parameters of a thundercloud, which in fact determine the risk of lightning activity in the area under consideration. The 3 h trace of thundercloud proposed in this paper seems to be a more appropriate function for comparison with the WWLLN data summed over the same period.
The analysis of the simulated distributions of the electric potential and electric field for the thunderstorm of 13–14 July allows us to make several remarks. The aerosol simulation reflects an earlier development of a thundercloud with high values for the electric field and potential and a longer duration of the thunderstorm. The results of the modeling using non-aerosol microphysical parametrization, in turn, show a later start of a thunderstorm compared to aerosol modeling and approximately the same end time of the thunderstorm.
The comparison of the simulation results with the radar and WWLLN data shows that the features of the initial stage of thunderstorm development are more correctly described by non-aerosol modeling; the end time of the thunderstorm is underestimated in simulations with both parametrizations of microphysics; the duration of the thunderstorm (disregarding the time offset) is correctly predicted by aerosol simulations. It is worth noting that on 14 July 2020, there was also a powerful thunderstorm in Moscow (located about 400 km from Nizhny Novgorod). If we consider these events as parts of the same cyclone, then we can assume that the error in the aerosol modeling is due to the fact that the main part of the thunderstorm shifted toward a more powerful aerosol island.
In this work, an attempt was made to count cells with high reflectivity in the radar data, which can probably be associated with the electroactive part of the cloud (that is, the part in which lightning activity is most likely). The applied thresholds have not been validated yet, so we will only briefly discuss such an observational tool. The volume of thunderclouds (the sum of volume cells) observed by a meteorological radar with a radar reflectivity of more than 45 dBZ is shown in
Figure 11 (upper panel). Temperature and altitude were not taken into account when calculating the volume of thunderclouds with high radar reflectivity. There is a fairly good agreement between the observed time dependence and the WWLLN data on the number of lightning discharges when such a threshold is used (see the lower panel of
Figure 11). In addition, it should be noted that the characteristic volume of a single thunderstorm cell is about 300 km
with characteristic dimensions of 5–7 km in radius and 2–3 km in height. For example, Mason [
53] gives an estimate of the horizontal scales of the cell of 5–14 km when observing the mature stage of the cell. A similar cell with a size of 260 km
was observed during the thunderstorm at 18:00 UTC on 14 July 2020.
When comparing the radar reflectivity obtained from observations and simulations, we found that the model overestimates the reflectivity values. The threshold value was selected according to the maximum volume of thunderstorm cells. The volume of simulated thunderclouds with a radar reflectivity of more than 60 dBZ is shown in
Figure 11 (central panel). The blue curve is used for the non-aerosol simulation, and the green curve is used for the aerosol one. When counting cell volume with a reflectivity above 45 dBZ (as for the radar data), we obtained values that are an order of magnitude larger than for radar observations. To obtain comparable maximum values, we had to raise the reflectivity threshold to 60 dBZ (see the central panel of
Figure 11). Previously, the results of radar reflectivity simulations using the WRF model and the WWLLN data for several convective seasons in the Nizhny Novgorod region were compared in [
15]. It was shown that, for the area under consideration, the best correlation was observed at a threshold value of 55 dBZ. Thus, some differences in the maximum values of the radar reflectivity obtained in the simulations and according to the observational data (see
Figure 6) can be associated with a systematic overestimation of this parameter in the simulation of convective events in this area.
As seen from the simulation results, both considered parametrizations of the microphysics have their advantages and disadvantages. For an accurate representation of all stages of thunderstorm development, the advantages of both parametrizations should be combined. In the aerosol parametrization of the microphysics considered in this paper, the initial distributions of aerosol mass mixing ratios are taken from multi-year global model simulations by the GOCART model. However, the same aerosol parametrization lets us multiply the initial distribution of the aerosol mass mixing ratio by any factor (see, for example, [
54]). Thus, as a further development of this study, it is possible to perform a set of simulations with different configurations of the initial distribution of the aerosol mass mixing ratio. If we assume that the initial distributions of the aerosol mass mixing ratios taken from the GOCART model are too rough (the spatial resolution of the GOCART model is 0.5 degrees in longitude and 1.25 degrees in latitude) or do not correspond sufficiently well to the observed ones in the territory under consideration, then we can affect the aerosol load by changing the initial distribution of the aerosols (for example, by increasing their number in the immediate neighborhood of the metropolis and decreasing on the periphery). Since the thunderstorm under consideration comes from the south of the region, this approach will probably make it possible to bring the results of modeling the beginning of a thunderstorm closer to the results of non-aerosol modeling. At the same time, a high aerosol load will remain over Nizhny Novgorod, which will probably maintain the lifetime of the convective system for a sufficiently long time.
5. Conclusions
In this paper, the annual variability in the number of lightning discharges in the territory of the Nizhny Novgorod region has been analyzed using the WWLLN data. It has been shown that the total number of discharges for the convective seasons of 2016–2021 increased. For the convective season of 2020, the daily number of discharges in the neighborhood of Nizhny Novgorod was calculated according to the WWLLN data and the EFM measurements. The most powerful thunderstorm was observed on 13–14 July 2020. The total number of discharges in the Nizhny Novgorod area with a radius of 38 km on 13–14 July exceeded 800 lightning events, which was several times higher than the number of lightning events during the other thunderstorms of the considered season. A few years earlier, a similar thunderstorm was studied [
52]. These thunderstorms have a number of similar features, which is why it can be assumed that such convective events can be classified as a separate type of thunderstorms, and similar parameters can be expected to be obtained both as a result of observations and numerical simulation.
A numerical simulation of the thunderstorm was performed using the WRF model coupled with a calculation of the intracloud electric parameters to study the parameters of the convective system. The simulation was carried out using two parametrizations of the microphysical processes to reveal the features associated with the allowance for aerosol particles. The data from the DMRL-C were used to validate the radar reflectivity simulation results of the WRF model. The results of modeling radar reflectivity are in good agreement with the observational data in space. According to the simulation using the aerosol parametrization of the microphysics, the thunderstorm begins earlier and lasts longer compared to the non-aerosol simulation. It is noted that the simulated thundercloud also has a wider thunderstorm front than seen on the radar.
The electric potential and the intracloud electric field were calculated to study the electric characteristics of the thunderstorm. In the aerosol simulation, the thunderstorm started developing faster, and the duration of the thunderstorm was longer. The WWLLN data on the detected lightning events were used to test the obtained simulation results. A comparison of the aerosol simulations with the WWLLN data did not confirm such an early development of lightning-active regions. In terms of the start time of the lightning activity, the results of the modeling using non-aerosol microphysics were in better agreement with the WWLLN data, but the duration of the thunderstorm obtained using non-aerosol parametrization was too short. The maximum value of the electric potential for the non-aerosol modeling exceeded the same parameter for the aerosol simulations. Thus, we can conclude that, when modeling severe thunderstorms in the Upper Volga region, the non-aerosol microphysical parametrization gives a more correct description of the thunderstorm development (with the exception of the final stage of the convective system development). For a convective event to have the required duration (more than 6 h), it is necessary to take aerosol particles into account.