Introduction

Recent trends in global and United States temperature have provoked questions about the impact on frequency, intensity, timing, and location of tornadoes. When removing many non-meteorological factors, it is shown that the annual frequency of United States tornadoes through the most reliable portions of the historical record has remained relatively constant.1,2,3,4 The most notable trends in tornado frequency are associated with increasing annual variability and a recent tendency for more tornadoes on any given tornado day.5,6,7,8,9 There has also been a measurable shift in the annual timing of tornado season, but it is currently unknown if this is due to rising global temperatures or natural variability.10,11 Two recent studies show increased tornado frequency in the eastern United States and decreased activity generally in the central United States, but analyses were limited to epochs.7,12 In addition, decreasing tornado frequency trends in the Great Plains and increasing trends in the Southeast are present in annual counts of E(F)1 + tornadoes.13 Despite these changes, detecting spatial shifts in tornado frequency has proved challenging given the relatively small spatial and temporal scales associated with tornado-producing thunderstorms and the well-documented deficiencies accompanying the tornado reporting database.1 The work herein seeks to examine spatial trends in tornado environments and tornado frequency at finer-scale resolutions in space and time.

These deficiencies have motivated researchers to focus on developing and analyzing atmospheric derived indices relating environmental conditions to tornado occurrence. Such indices are often used in an ingredients-based approach commonly utilized in operational weather forecasting.14 The significant tornado parameter (STP) is one such index, originally designed as a statistical discriminator between non-tornado and significant ([E]F2 or greater using the Fujita damage scale) tornado environments.15 STP is a commonly used operational weather forecasting index, but no research has used STP to examine the covariant relationship with aspects of the United States tornado report climatology.

A derived covariate, such as STP, is complementary to tornado reports for climatological studies. For example, it has been shown that variables like storm relative helicity and convective precipitation adequately represent United States monthly tornado frequency.16 Reports are subject to a human reporting process, whereas environmental covariates allow for an objective model-derived climatology. One major limitation of the covariate methodology for tornadoes is the lack of accounting for a thunderstorm initiation metric. It is well documented that tornadoes occur in environments characteristic of adequate vertical wind shear, relatively high surface water vapor mixing ratios, low static stability, and directional change in the lowest 1-km of tropospheric wind.17,18 Yet, even with a supportive atmosphere, thunderstorms may not form in the absence of an initial lifting mechanism needed to commence the convective process. The space and time scales of various lifting sources are often below the spatial scale of current reanalysis observation capabilities, and thus, not accounted for by most covariate environmental metrics. This often leads to an overestimation of favorable tornado environments and a lack of explanatory capability regarding observed tornado reports, especially when convective inhibition is present.19

Recent increases in computational efficiency have led several studies to explicitly depict severe convective storms in convective permitting regional climate model configurations to evade the previously discussed environmental caveats.20,21,22,23,24,25 These dynamical downscaling approaches (sometimes referred to as model telescoping) have proven reliable for the replication of spatial frequency and temporal timing of the annual peak of boreal severe weather activity. Potential changes in characteristics of late 21st century severe weather activity have also been posited by these methods. However, due to the resource intensive nature of these approaches, long model-derived climatologies are often unfeasible and it is not yet clear if the explanatory capability is significantly better than environments to justify the resources necessary to conduct such simulations.

This research focuses on an environmental covariate approach to examine potential changes in United States tornado frequency. In particular, we begin by analyzing an aggregated covariate methodology in an attempt to explain the variance associated with monthly and annual United States tornado reports. We further examine the spatio-temporal changes in this covariate over the past ~40 years. This allows us to draw conclusions about the regional trends in frequency of tornado environments. We conclude by comparing the covariate trend analyses to actual tornado reporting trends to strengthen the presented results. Ultimately, our research highlights that trends in United States tornado frequency are regional in nature and this change in physical risk has implications for tornado exposure and vulnerability.

Results

Despite the aforementioned environmental caveats, annual accumulated diurnal max STP is a statistically significant covariate to annual United States tornado reports (Fig. 1). This degree of explanatory power (44% of the variance in annual tornado counts is explained by accumulated diurnal maximum STP) suggests that STP may be used as a climatological proxy for annual tornado report counts. At the monthly scale, accumulated STP exhibits a seasonal trend in explanatory power. Monthly standardized anomalies of accumulated diurnal maximum STP for January, February, March, May, and December all explain roughly half of the variance associated with the standardized anomaly of monthly United States tornado counts (Fig. 2). August displays a minimum in variance explanation (9%), and it is speculated this is a function of the mesoscale nature of late boreal summer tornado environments and larger contributions to the monthly count climatologies from tropical storms and hurricanes. Standardized anomaly regression slopes reveal that monthly accumulated STP is a conservative estimator of monthly tornado counts, and this bias is greatest in the boreal summer (Fig. 2). The null hypothesis of, “Monthly accumulated STP does not explain statistically significant amounts of variance in monthly tornado counts” was rejected at the 0.01 significance level for all months except August. Thus, the STP covariate can be applied at the monthly scale, but caution should be used due to the seasonal changes in the coefficient of determination and level of statistical significance.

Fig. 1
figure 1

Scatter plot of United States annual tornado counts vs. United States annual accumulated significant tornado parameter (1979–2017) from NARR

Fig. 2
figure 2

1979–2017 Unites States monthly tornado reports and accumulated diurnal maximum STP (expressed here as a standardized anomaly by month) from NARR. Regression r2 and p values are listed next to their respective months

Examining annual tornado report counts and accumulated annual STP through time confirm this explanatory relationship and give further confidence in the use of STP as a climatological tornado report proxy (Fig. 3). Both tornado counts and STP indicate a slight positive trend through time over the analysis period, which starts at the beginning of the modern Geostationary Operational Environmental Satellite (GOES) era. Consistent tornado reporting data exists back to the early 1950s, but there are currently no datasets that allow for the space and time scales necessary to adequately retroactively compute STP back to 1950, and hence, we are limited here to the analysis period 1979–2017. Most of the increase in tornado reporting is observed from the late 1980s to the early 2000s, with a steady or even declining trend thereafter which is consistent with previous research.3

Fig. 3
figure 3

Interannual variability of United States tornado reports (blue; axis left) and the United States annual sum of daily max significant tornado parameter (red; axis right) from NARR

Of greatest interest here is the potential for shifts in the spatial location of tornado frequency. Research has identified evidence of a “Dixie Alley,” which represents an eastward extension of the traditional "Tornado Alley" in the central Great Plains.26 However, this could be due to the techniques used to smooth reports in space and time.27 An eastward shift in tornado frequency was also suggested, but results were limited to two arbitrary regions, coarse spatial environmental covariates, and lacked a robust environmental tornado report proxy such as STP.12,13 Using Theil-Sen slope analysis—a robust linear trend estimator which is insensitive to outliers—and Kendall’s τ test for statistical significance, we show a robust downward trend in the annual accumulation of STP across the central and southern Great Plains. Meanwhile, a robust upward trend is found in portions of the Southeast, Midwest, and Northeast (Fig. 4). No significant increase (decrease) in tornado environments is observed west (east) of the 95th meridian. We believe these trends in tornado environments are significant and have not been documented with this level of detail by previous research.

Fig. 4
figure 4

Theil-Sen slope analysis of 1979–2017 annual grid-point sum of daily max STP from NARR. p values are hatched at values ≤ 0.05 significance using Kendall’s τ statistic. Slope units are sum of daily max STP per year

Given these environmental tendencies, one would also suspect to find such evidence of trends in tornado reports given the explanatory power of the STP covariate. Indeed, significant trends in the spatial frequency of annual tornado reports are found over the analysis period (Fig. 5). Both tornado reports and environments indicate significant decreasing trends in frequency over portions of Texas, Oklahoma, and northeast Colorado. Agreement in the sign of the Theil-Sen slope is also noted for significant increasing trends in portions of Mississippi, Alabama, Arkansas, Missouri, Illinois, Indiana, Tennessee, and Kentucky. We are most confident about the trends in these areas, where both tornado reports and environments share the same significant sign of Theil-Sen’s slope. Confidence in the trends of tornado frequency are reduced in areas where the sign of Theil-Sen’s slope is opposite or do not share statistical significance. For example, significant increases in tornado reporting trends are found across Kansas, but environments indicate a neutral or downward trend. Additionally, tornado environments indicate a robust increasing trend in portions of Arkansas and Louisiana that are not supported by tornado reporting trends. We cannot offer a physical reason for these discrepancies at the current time, but they may be a function of non-continuous regional/seasonal efficiencies of the STP covariate (e.g., biases shown in Fig. 2) and biases in the tornado report database. We speculate the non-continuous regional/seasonal efficiencies of STP are a function of convective initiation and forcing for ascent. In addition, STP was designed with the intention of evaluating environments favorable for right-moving supercells in the Great Plains of the U.S. and not necessarily non-mesocyclone tornadoes.15,28 However, results associated with reports herein are generally consistent with an increasing trend found in the Southeast United States due to a greater propensity for tornado outbreaks east of the Mississippi River.29,30

Fig. 5
figure 5

Theil-Sen slope analysis of 19792017 annual gridded tornado reports. p values are hatched at values ≤ 0.05 significance using Kendall’s τ statistic. Slope units are reports per year * 10−1

Further examination reveals that spatial trends in tornado environment frequency are modulated by seasons (Fig. 6). Positive trends in tornado environment frequency are shown to progress northward through the annual cycle. Most of the decreasing trends across Texas are due to the months of March, April, and May (Fig. 6b), and most of the Great Plains has shown a decreasing trend in tornado environments during June, July, and August (Fig. 6c). It should be noted that our results do not indicate that maximum United States tornado environment and report frequency is no longer in the Great Plains. Rather, we demonstrate there has been a robust downward trend of tornado environments and reports in portions of the Great Plains concurrent with a simultaneous robust increasing trend in areas surrounding the Mississippi River Valley and the Midwest.

Fig. 6
figure 6

Theil-Sen slope of 1979–2017 annual grid-point sum of daily max STP for a) December, January, February; b) March, April, May; c) June, July, August; and d) September, October, November. p values are hatched at values ≤ 0.05 significance using Kendall’s τ statistic. Slope units are sum of daily max STP per year

Discussion

The presented trends in spatial tornado frequency are important for understanding the tornado as a hazard, which involves examining the comingling of physical risk (i.e., climatological probability of a tornado report) and societal vulnerability (e.g., population density). Such increases in tornado frequency across the Mid-South juxtapose with an already vulnerable population to tornadoes. In fact, this geographic area of the United States already represents a maximum in the occurrence of casualties associated with tornadoes.31 Climate downscaling and societal vulnerability studies also highlight the Mid-South as having the greatest potential for increased tornado disasters by the end of the century.23,25,32,33

Economic losses associated with tornadoes will continue to increase in future years, both due to the spatial changes presented herein and the expanding footprint of the built environment.34 It is noted that the impacts associated with projected 21st century increases in tornado frequency are outweighed by projected growth in the human-built environment. However, the combination of an increase in risk and exposure could lead to a threefold increase in tornado disaster potential.32,33 Thus, it is important to continue monitoring the presented physical trends in tornado frequency and recognize that they have consequences for increases in tornado disaster potential.

At this point, it is unclear whether the observed trends in tornado environment and report frequency are due to natural variability or being altered by anthropogenic forcing on the climate system. Global and regional climate model studies indicate a general increase in severe weather frequency (especially in boreal spring) and variability by the end of the 21st century in portions of the Midwest and Southeast, suggesting that the trends herein may indeed be due to anthropogenic forcing given the accumulating literature examining tornado/severe weather frequency and variability.3 This eastward shift is also consistent with recent work suggesting an eastward shift in the “effective 100th meridian” as the century progresses.35 It is difficult to isolate these natural vs. anthropogenic factors, but perhaps future attribution studies will be able to describe such causal mechanisms.

Methods

Tornado data were retrieved from the Storm Prediction Center’s online database (described online at http://www.spc.noaa.gov/wcm/#data] for the period 1979–2017. All tornadoes were used in this analysis. Data from the North American Regional Reanalysis (NARR; 36) was used to calculate the STP (a unitless derived environmental variable) covariate. 3-hourly data were downloaded from the NCAR/UCAR research data archive portal (https://rda.ucar.edu/datasets/ds608.0/) for the period 1979–2017.

STP was calculated following a fixed-layer15 approach:

$${\mathrm{STP}} = \frac{{{\mathrm{SBCAPE}}}}{{{\mathrm{1500}}}} \times \frac{{\left( {{\mathrm{2000 - SBLCL}}} \right)}}{{{\mathrm{1000}}}} \times \frac{{{\mathrm{SRH01}}}}{{{\mathrm{150}}}} \times \frac{{{\mathrm{BWD06}}}}{{{\mathrm{20}}}}$$

using surface-based convective available potential energy (SBCAPE), surface-based lifting condensation level (SBLCL), 0–1 km storm relative helicity (SRH01), and 0–6 km bulk wind difference (BWD06). Since NARR data use an isobaric vertical coordinate system, data were vertically interpolated from NARR isobaric levels to AGL height coordinates for AGL calculations (e.g., 0–6 km). In addition, STP was masked in the presence of convective inhibition < −50 J kg−1 in an attempt to examine environments most favorable for tornadoes. STP was computed on the native NARR 32-km lambert conformal grid. An accumulated STP metric was derived as a covariate for comparison to the annual tornado report counts by finding the diurnal (1200–1200 UTC) max STP value from NARR 3-hourly data and summing all CONUS land point daily maximums over the year to create annual accumulation values. 1200–1200 UTC was chosen for the daily cycle to capture diurnally driven supercell and tornado climatology.37 Annual diurnal max accumulations were used to explain annual tornado report variability. For monthly regressions, a standardized anomaly was calculated for both monthly tornado counts and monthly accumulation of diurnal maximum STP. Regression slopes of these standardized anomaly relationships were examined to demonstrate monthly explanatory capability, variability, and bias.

For spatial trend analysis of the STP covariate, a Theil-Sen estimator was calculated at each native NARR grid point on the interannual accumulation values. Theil-Sen slope analysis was chosen due its efficient computation and insensitivity to outliers.38 Significance of the Theil-Sen slope was assessed using Kendall’s τ statistic and a p value at the 0.05 significance level. Since tornado reports are point objects in space and time—and not gridded—reports were binned in space on a 1° latitude by 1° longitude grid using the same lambert conformal projection as NARR data and smoothed in space using a 1-σ Gaussian kernel. Theil-Sen slope analysis and Kendall’s τ statistic were then calculated on interannual gridded tornado reports using the same methodology for accumulated STP. The authors note that due to the number of statistical tests, familywise error rate, and inherent autocorrelation in time and space of climatological data, such analyses are potentially impacted by the multiple testing problem.38