Abstract Eighteen lightning flash rate parameterization schemes (FRPSs) were investigated in a Weather Research and Forecasting model coupled with chemistry cloud‐resolved simulation of the 29–30 May 2012 supercell storm system observed during the Deep Convective Clouds and Chemistry (DC3) field campaign. Most of the observed storm's meteorological conditions were well represented when the model simulation included both convective damping and lightning data assimilation techniques. Newly‐developed FRPSs based on DC3 radar observations and Lightning Mapping Array data are implemented in the model, along with previously developed schemes from the literature. The schemes are based on relationships between lightning and various kinematic, structural, and microphysical thunderstorm characteristics (e.g., cloud top height, hydrometeors, reflectivity, and vertical velocity) available in the model. The results suggest the model‐simulated graupel and snow/ice hydrometeors require scaling factors to more closely represent proxy observations. The model‐simulated lightning flash trends and total flashes generated by each scheme over the simulation period are compared with observations from the central Oklahoma Lightning Mapping Array. For this supercell system, 13 of the 18 schemes overpredicted flashes by >100% with the group of FRPSs based on storm kinematics and structure (particularly updraft volume) performing slightly better than the hydrometeor‐based schemes. During the storm's first 4 hr, the upward cloud ice flux FRPS, which is based on the combination of vertical velocity and hydrometeors, well represents the observed total flashes and flash rate trend; while, the updraft volume scheme well represents the observed flash rate peak and subsequent sharp decline in flash rate.

Abstract A cloud‐resolved storm and chemistry simulation of a severe convective system in Oklahoma constrained by anvil aircraft observations of NO x was used to estimate the mean production of NO x per flash in this storm. An upward ice flux scheme was used to parameterize flash rates in the model. Model lightning was also constrained by observed lightning flash types and the altitude distribution of flash channel segments. The best estimate of mean NO x production by lightning in this storm was 80–110 mol per flash, which is smaller than many other literature estimates. This result is likely due to the storm having been a high flash rate event in which flash extents were relatively small. Over the evolution of this storm a moderate negative correlation was found between the total flash rate and flash extent and energy per flash. A longer‐term simulation at 36‐km horizontal resolution with parameterized convection was used to simulate the downwind transport and chemistry of the anvil outflow from the same storm. Convective transport of low‐ozone air from the boundary layer decreased ozone in the anvil outflow by up to 20–40 ppbv compared with the initial conditions, which contained stratospheric influence. Photochemical ozone production in the lightning‐NO x enhanced convective plume proceeded at a rate of 10–11 ppbv per day in the 9–11 km outflow layer over the 24‐hr period of downwind transport to the Southern Appalachians. Photochemical production plays a large role in the restoration of upper tropospheric ozone following deep convection.