Numerical weather prediction ensembles are routinely used for operational weather forecasting. The members of these ensembles are individual simulations with either slightly perturbed initial conditions or different model parameterizations, or occasionally both. Multi-member ensemble output is usually large, multivariate, and challenging to interpret interactively. Forecast meteorologists are interested in understanding the uncertainties associated with numerical weather prediction; specifically variability between the ensemble members. Currently, visualization of ensemble members is mostly accomplished through spaghetti plots of a single midtroposphere pressure surface height contour. In order to explore new uncertainty visualization methods, the Weather Research and Forecasting (WRF) model was used to create a 48-hour, 18 member parameterization ensemble of the 13 March 1993 "Superstorm". A tool was designed to interactively explore the ensemble uncertainty of three important weather variables: water-vapor mixing ratio, perturbation potential temperature, and perturbation pressure. Uncertainty was quantified using individual ensemble member standard deviation, inter-quartile range, and the width of the 95% confidence interval. Bootstrapping was employed to overcome the dependence on normality in the uncertainty metrics. A coordinated view of ribbon and glyph-based uncertainty visualization, spaghetti plots, iso-pressure colormaps, and data transect plots was provided to two meteorologists for expert evaluation. They found it useful in assessing uncertainty in the data, especially in finding outliers in the ensemble run and therefore avoiding the WRF parameterizations that lead to these outliers. Additionally, the meteorologists could identify spatial regions where the uncertainty was significantly high, allowing for identification of poorly simulated storm environments and physical interpretation of these model issues.

The term “Tornado Alley” is a gross approximation of the most tornado-prone region in the United States. Depending on calculation methods, Tornado Alley can vary dramatically across the area between the Rocky and Appalachian Mountains. There is some evidence that multiple alleys of peak tornado activity exist around the country, including “Dixie Alley” in the Southeast. Therefore, we assess the spatial tornado risk and seek any regions of elevated tornado risk that are distinctly separate from the traditional Tornado Alley of the Great Plains. Results show there are no tornado risk areas statistically separate from Tornado Alley, but there are large portions of the Southeast that experience more tornadoes than the rest of the country. It appears that Tornado Alley and Dixie Alley are part of a single large region of high tornado risk with a relative minimum near the middle due to the Ozark and Ouachita Mountains. Placement of the maximum tornado density in Mississippi, along with other regions of relative maxima across the Southeast, may warrant modification of the traditional tornado risk map that focuses only on the Great Plains. Understanding such patterns is important for preparing the public and mitigating tornado hazards.

This paper presents an efficient visualization and exploration approach for modeling and characterizing the relationships and uncertainties in the context of a multidimensional ensemble dataset. Its core is a novel dissimilarity-preserving projection technique that characterizes not only the relationships among the mean values of the ensemble data objects but also the relationships among the distributions of ensemble members. This uncertainty-aware projection scheme leads to an improved understanding of the intrinsic structure in an ensemble dataset. The analysis of the ensemble dataset is further augmented by a suite of visual encoding and exploration tools. Experimental results on both artificial and real-world datasets demonstrate the effectiveness of our approach.

Tornadic and nontornadic outbreaks occur within the United States and elsewhere around the world each year with devastating effect. However, few studies have considered the physical differences between these two outbreak types. To address this issue, synoptic-scale pattern composites of tornadic and nontornadic outbreaks are formulated over North America using a rotated principal component analysis (RPCA). A cluster analysis of the RPC loadings group similar outbreak events, and the resulting map types represent an idealized composite of the constituent cases in each cluster. These composites are used to initialize a Weather Research and Forecasting Model (WRF) simulation of each hypothetical composite outbreak type in an effort to determine the WRF’s capability to distinguish the outbreak type each composite represents. Synoptic-scale pattern analyses of the composites reveal strikingly different characteristics within each outbreak type, particularly in the wind fields. The tornado outbreak composites reveal a strong low- and midlevel cyclone over the eastern Rockies, which is likely responsible for the observed surface low pressure system in the plains. Composite soundings from the hypothetical outbreak centroids reveal significantly greater bulk shear and storm-relative environmental helicity values in the tornado outbreak environment, whereas instability fields are similar between the two outbreak types. The WRF simulations of the map types confirm results observed in the composite soundings.

Tornadoes often strike as isolated events, but many occur as part of a major outbreak of tornadoes. Nontornadic outbreaks of severe convective storms are more common across the United States but pose different threats than do those associated with a tornado outbreak. The main goal of this work is to distinguish between significant instances of these outbreak types objectively by using statistical modeling techniques on numerical weather prediction output initialized with synoptic-scale data. The synoptic-scale structure contains information that can be utilized to discriminate between the two types of severe weather outbreaks through statistical methods. The Weather Research and Forecast model (WRF) is initialized with synoptic-scale input data (the NCEP–NCAR reanalysis dataset) on a set of 50 significant tornado outbreaks and 50 nontornadic severe weather outbreaks. Output from the WRF at 18-km grid spacing is used in the objective classification. Individual severe weather parameters forecast by the model near the time of the outbreak are analyzed from simulations initialized at 24, 48, and 72 h prior to the outbreak. An initial candidate set of 15 variables expected to be related to severe storms is reduced to a set of 6 or 7, depending on lead time, that possess the greatest classification capability through permutation testing. These variables serve as inputs into two statistical methods, support vector machines and logistic regression, to classify outbreak type. Each technique is assessed based on bootstrap confidence limits of contingency statistics. An additional backward selection of the reduced variable set is conducted to determine which variable combination provides the optimal contingency statistics. Results for the contingency statistics regarding the verification of discrimination capability are best at 24 h; at 48 h, modest degradation is present. By 72 h, the contingency statistics decline by up to 15%. Overall, results are encouraging, with probability of detection values often exceeding 0.8 and Heidke skill scores in excess of 0.7 at 24-h lead time.