A Year-Long Total Lightning Forecast over Italy with a Dynamic Lightning Scheme and WRF

Lightning is an important threat to life and properties and its forecast is important for several applications. In this paper, we show the performance of the “dynamic lightning scheme” for next-day total strokes forecast. The predictions were compared against strokes recorded by a ground observational network for a forecast period spanning one year. Specifically, a total of 162 case studies were selected between 1 March 2020 and 28 February 2021, characterized by at least 3000 observed strokes over Italy. The events span a broad range of lightning intensity from about 3000 to 600,000 strokes in one day: 69 cases occurred in summer, 46 in fall, 18 in winter, and 29 in spring. The meteorological driver was the Weather Research and Forecasting (WRF) model (version 4.1) and we focused on the next-day forecast. Strokes were simulated by adding three extra variables to WRF, namely, the potential energies for positive and negative cloud to ground flashes and intracloud strokes. Each potential energy is advected by WRF, it is built by the electrification processes occurring into the cloud, and it is dissipated by lightning. Observed strokes were remapped onto the WRF model grid with a 3 km horizontal resolution for comparison with the strokes forecast. Results are discussed for the whole year and for different seasons. Moreover, statistics are presented for the land and the sea. In general, the results of this study show that lightning forecast with the dynamic lightning scheme and WRF model was successful for Italy; nevertheless, a careful inspection of forecast performance is necessary for tuning the scheme. This tuning is dependent on the season. A numerical experiment changing the microphysics scheme used in WRF shows the sensitivity of the results according to the choice of the microphysics scheme.

Evaluation of Different Storm Parameters as the Proxies for Gridded Total Lightning Flash Rates: A Convection-Allowing Model Study

Lightning simulation is important for a variety of applications, including lightning forecast, atmospheric chemical simulation, and lightning data assimilation. In this study, the potential of five storm parameters (graupel volume, precipitation ice mass, radar echo volume, maximum updraft, and updraft volume) to be used as the proxy for the diagnosis of gridded total lightning flash rates has been investigated in a convection-allowing model. A mesoscale convective system occurred in the Guangdong province of China was selected as the test case. Radar data assimilation was used to improve the simulation accuracy of the convective clouds, hence providing strong instantaneous correlations between observed and simulated storm signatures. The areal coverage and magnitude of the simulated lightning flash rates were evaluated by comparing to those of the total lightning observations. Subjective and the Fractions Skill Score (FSS) evaluations suggest that all the five proxies tested in this study are useful to indicate general tendencies for the occurrence, region, and time of lightning at convection-allowing scale (FSS statistics for the threshold of 1 flash per 9 km2 per hour were around 0.7 for each scheme). The FSS values were decreasing as the lightning flash rate thresholds used for FSS computation increased for all the lightning diagnostic schemes with different proxies. For thresholds from 1 to 3 and 16 to 20 flashes per 9 km2 per hour, the graupel contents related schemes achieved higher FSS values compared to the other three schemes. For thresholds from 5 to 15 flashes per 9 km2 per hour, the updraft volume related scheme yielded the largest FSS. When the thresholds of lightning flash rates were greater than 13 flashes per 9 km2 per hour, the FSS values were below 0.5 for all the lightning diagnostic schemes with different proxies.

An Evaluation of Relationships between Radar-Inferred Kinematic and Microphysical Parameters and Lightning Flash Rates in Alabama Storms

Lightning flash rate parameterizations based on polarimetric and multi-Doppler radar inferred microphysical (e.g., graupel volume, graupel mass, 35 dBZ volume) and kinematic (e.g., updraft volume, maximum updraft velocity) parameters have important applications in atmospheric science. Although past studies have established relations between flash rate and storm parameters, their expected performance in a variety of storm and flash rate conditions is uncertain due to sample limitations. Radar network and lightning mapping array observations over Alabama of a large and diverse sample of 33 storms are input to hydrometeor identification, vertical velocity retrieval and flash rate algorithms to develop and test flash rate relations. When applied to this sample, prior flash rate linear relations result in larger errors overall, including often much larger bias (both over- and under-estimation) and root mean square errors compared to the new linear relations. At low flash rates, the new flash rate relations based on kinematic parameters have larger errors compared to those based on microphysical ones. Sensitivity of error to the functional form (e.g., zero or non-zero intercept) is also tested. When considering all factors (e.g., low errors including at low flash rate, consistency with past linear relations, and insensitivity to functional form), the flash rate parameterization based on graupel volume has the best overall performance.