What is the intrinsic predictability of tornadic supercell thunderstorms?

A 25-member ensemble of relatively high-resolution (75-m horizontal grid spacing) numerical simulations of tornadic supercell storms is used to obtain insight on their intrinsic predictability. The storm environments contain large and directionally varying wind shear, particularly in the boundary layer, large convective available potential energy, and a low lifting condensation level. Thus, the environments are extremely favorable for tornadic supercells. Small random temperature perturbations present in the initial conditions trigger turbulence within the boundary layers. The turbulent boundary layers are given 12 h to evolve to a quasi-steady state before storms are initiated via the introduction of a warm bubble. The spatially averaged environments are identical within the ensemble; only the random number seed and/or warm bubble location is varied. All of the simulated storms are long-lived supercells with intense updrafts and strong mesocyclones extending to the lowest model level. Even the storms with the weakest near-surface rotation probably can be regarded as weakly tornadic. However, despite the statistically identical environments, there is considerable divergence in the finescale details of the simulated storms. The intensities of the tornado-like vortices that develop in the simulations range from EF0 to EF3, with large differences in formation time and duration also being exhibited. The simulation differences only can be explained by differences in how the initial warm bubbles and/or storms interact with turbulent boundary layer structures. The results suggest very limited intrinsic predictability with respect to predicting the formation time, duration, and intensity of tornadoes.