Abstract

Winter and spring tornado activity tends to be heightened during the La Niña phase of the El Niño/Southern Oscillation and suppressed during the El Niño phase. Despite these tendencies, some La Niña seasons have fewer tornadoes than expected and some El Niño seasons have more than expected. To gain insight into such anomalous seasons, the two La Niña winters and springs with the fewest tornadoes and the two El Niño winters and springs with the most tornadoes between 1979 and 2016 are identified and analyzed in this study. The relationships between daily tornado count and the Global Wind Oscillation and Madden-Julian Oscillation in these anomalous seasons are also explored. Lastly, seasonal and daily composites of upper-level flow, low-level flow and humidity, and atmospheric instability are generated to describe the environmental conditions in the anomalous seasons. The results of this study highlight the potential for large numbers of tornadoes to occur in a season if favorable conditions emerge in association with individual synoptic-scale events, even during phases of the El Niño/Southern Oscillation, Global Wind Oscillation, and Madden-Julian Oscillation that seem to be unfavorable for tornadoes. They also highlight the potential for anomalously few tornadoes in a season even when the oscillations are in favorable phases.

1. Introduction

More tornadoes occur in the United States (US) per year than in any other country [1], and these tornadoes are capable of producing incredible economic and human loss [25]. Despite their relatively common occurrence, there is notable intra- and interannual variability in the number of US tornadoes [610]. Due to the danger and variability of tornado activity, there is a need to improve seasonal outlook capabilities. One approach is to identify relationships between seasonal tornado count and climate oscillations, such as the El Niño/Southern Oscillation (ENSO) (e.g., [1119]).

Numerous studies over the past couple of decades have analyzed the relationship between tornadoes and ENSO. An early study by Monfredo [15] reported that strong and violent tornadoes were more common during the La Niña (LN) phase of ENSO and less common during the El Niño (EN) phase. Cook and Schaefer [16] later reported that winter tornado outbreaks were stronger and more frequent during the neutral phase (N) of ENSO, followed by the LN then EN phases. The most recent studies report that tornado activity is heightened during the LN phase [1214, 1719]. Recent studies also illustrate that the seasons with the most tornadoes tend to be classified as LN. Lee et al. [17], for example, analyzed the number of tornadoes rated 3 or higher on the Fujita or Enhanced Fujita damage scales (hereafter referred to as E(F)) and reported that the five most extreme tornado outbreak years were characterized as persistent LN events or LN events that were transitioning to a different phase. Allen et al. [12] similarly reported that most of the seasons with a tornado count >100% of climatology occurred during the LN phase, whereas most that were <75% of climatology occurred during the EN phase.

The relatively consistent reporting of heightened tornado activity with the LN phase of ENSO elicits the possibility of predicting the likelihood of below- or above-normal tornado activity based on ENSO conditions. Elsner et al. [18] developed a statistical model to evaluate variability in tornado activity in which ENSO was the most important predictor. Their model also indicated that tornado activity was heightened in the Midwest and Southeast regions of the US during LN and in the Great Plains region during EN. Allen et al. [12] and Lepore et al. [13] developed extended logistic regression models using ENSO to predict the likelihood of below-normal, normal, or above-normal tornado activity in spring. Their models similarly showed that tornado activity is increased during the LN phase, particularly in the south central US.

Anomalous seasons occur. Some LN seasons have fewer than expected tornadoes and some EN seasons have more than expected. Lee et al. [19], for example, noted that an anomalously large number of tornadoes occurred in the 2015-2016 EN winter. Elsner et al. [18] also noted the occurrence of anomalous seasons and attributed them to additional unknown factors. While previous studies have noted these anomalous seasons, none have provided detailed descriptions and analyses of them. More research, including work focusing on the anomalous seasons, is needed to improve our understanding of the tornado-ENSO link and to refine statistical models.

A case study approach was taken here to describe and analyze the two EN winters and springs with the most tornadoes, and the two LN winters and springs with the fewest tornadoes. The states of the Madden-Julian Oscillation (MJO) and Global Wind Oscillation (GWO) are of particular interest here. Both of these oscillations vary over subseasonal timescales (i.e., shorter than the timescale over which ENSO varies) and have been linked to variations in tornado activity [2023]. Seasonal and daily composites are also generated to document the synoptic patterns that accompanied the anomalous seasons.

The specific objectives of this study are to(1)identify and describe the two EN winters and springs with the most tornadoes and the two LN winters and springs with the fewest tornadoes;(2)analyze GWO and MJO activity during these anomalous seasons;(3)document the synoptic patterns that accompany these anomalous seasons.

2. Data and Methods

Tornado data were taken from the Storm Prediction Center’s Severe Weather Database (specifically, the 1950–2016_actual_tornadoes.csv file) [24]. This dataset was a subset to the 18,251 E(F)1+ tornadoes that occurred over the period December 1978–November 2016 in the contiguous US. The time series of E(F)1+ tornadoes is less influenced by detection and reporting changes over time than is the E(F)0+ series, which includes a notable upward trend [25, 26]. The beginning of the period of record is consistent with recent studies [12, 13] and corresponds to the first year of the North American Regional Reanalysis (NARR) product [27], which is used to generate seasonal and daily atmospheric composites. Seasonal tornado counts were generated for winter (December of previous year, January, and February (DJF)) and spring (March, April, and May (MAM)). DJF and MAM are the focus of this study because previous studies illustrate that tornado counts and ENSO are related in these seasons [1214, 1719].

Multiple indices represent ENSO conditions. Some represent mostly the oceanic component of ENSO (e.g., the Niño 1.2, 3, 3.4, and 4 regions and the Oceanic Niño Index (ONI)), some represent mostly the atmospheric component (e.g., the Southern Oscillation Index (SOI)), and some combine both (e.g., the Multivariate ENSO Index (MEI) and Bivariate ENSO Index (BEST)). The classification of seasons as LN or EN will likely change in some cases, depending on the index. In this study, ONI data were obtained from the Climate Prediction Center [28] to represent ENSO conditions, which is consistent with recent efforts to predict the likelihood of an active or inactive season based on the state of ENSO [12, 13]. ONI is also one of the commonly used indices by the Climate Prediction Center to diagnose ENSO conditions and to place current events into historical perspective. It is given as anomalies of 3-month running means of sea surface temperature in the Niño 3.4 region (5°N–5°S, 120°W–170°W). Series of ONI for DJF and MAM were extracted and merged with the seasonal tornado counts. Once merged, seasons were classified as LN if ONI ≤−0.5 and EN if ONI ≥0.5. Seasons with −0.5 < ONI < 0.5 were classified as neutral (N). The two most active EN DJFs and MAMs (i.e., with the most tornadoes) and two least active LN DJFs and MAMs (i.e., with the fewest tornadoes) were extracted, yielding eight seasons as the focus of this study.

Climate oscillations other than ENSO have been linked to tornado activity. Those that vary over subseasonal timescales may provide insight into the intraseasonal distribution of tornadoes in the anomalous seasons. Two such oscillations are the MJO and GWO. The MJO is a tropical wave of enhanced convection that propagates the globe from west to east at approximately 5 m·s−1 [29]. It is linked to variability in tornado activity in the midlatitudes through changes in global relative atmospheric angular momentum (AAM), which arise from tropical convective forcing and teleconnected Rossby wave propagation in the midlatitudes [20, 23, 3032]. Other factors, such as friction and mountain torque, affect AAM in addition to tropical convection [3032]. The GWO combines these factors and, therefore, provides a comprehensive representation of AAM and midlatitude circulation patterns [31, 32]. Daily GWO data were obtained from the Earth System Research Laboratory’s GWO dataset [33]. The attributes of the GWO used in this study are AAM anomalies, AAM time tendency, and GWO phase. The relationship between AAM anomalies and AAM time tendency determines the phase of GWO, which run from 1 to 8. Phases 2 and 3 represent anomalously low AAM, whereas phases 6 and 7 represent anomalously high AAM. Phases 1, 4, 5, and 8 represent transition states. Daily MJO data were obtained from the Australian Bureau of Meteorology’s Real-Time Multivariate (RMM) MJO dataset [34, 35]. The attributes used here are the first two principal components of empirical orthogonal functions that consist of a meridional average of the 200 mb zonal wind, the 850 mb zonal wind, and outgoing longwave radiation between 15°N and 15°S (referred to as RMM1 and RMM2) and the MJO phase. Similar to GWO, the MJO index has eight phases, but they represent the location of enhanced tropical convection. Phase 1 of the MJO indicates that the enhanced convection is located near eastern Africa and the subsequent phases represent its eastward progression.

AAM anomalies, AAM time tendency, RMM1, and RMM2 were used to construct GWO and MJO phase space diagrams for each of the anomalous seasons. Daily tornado counts were generated (where a day is defined as 0000–2359 UTC) and plotted in the phase space diagrams along with the progression of GWO and MJO. Daily tornado counts were also assessed across the phases of MJO and GWO, and Kruskal–Wallis tests were used to determine if the mean rank of the tornado counts significantly varied across the phases of the oscillations. This nonparametric test was chosen because of the skewed nature of daily tornado counts.

The synoptic patterns associated with the anomalous seasons were characterized using gridded composites from NARR [27]. The following variables were chosen to be comparable to composites in previous research (e.g., [12, 20, 21]): upper-level flow was characterized using 300 mb geopotential heights; low-level flow was characterized using 850 mb geopotential heights; low-level moisture was characterized using 850 mb specific humidity; and atmospheric stability was characterized using surface-based convective available potential energy (CAPE). Surface-based CAPE was used because Gensini et al. [36] illustrate that surface-based parcels are more accurate than 100 mb mixed layer parcels, largely due to errors in low-level moisture fields. The seasonal mean composites are averages of the variables for the dates under consideration (e.g., mean 300 mb geopotential height in DJF 2016 is generated by averaging the 300 mb geopotential height for the days in DJF 2016). The seasonal anomaly composites represent the difference between the mean composites and the 1981–2010 climatology. Daily composites are averages of the 3 hr NARR data. All composites were generated with NOAA’s Earth System Research Laboratory’s online plotting tool [37, 38]. Once generated, the Network Common Data Form (NetCDF) of the composites was imported into ArcMAP [39] for plotting and display.

3. Results and Discussion

3.1. Identification and Description of the Anomalous Seasons

The relationships between the number of E(F)1+ tornadoes and ONI in DJF and MAM are depicted in Figure 1. The negative relationship reported by others, whereby tornado frequency tends to be greater during the LN phase and lesser during the EN phase [1214, 1719], is visually apparent, but the correlation is not statistically significant in this subset. The central tendencies in Table 1 also illustrate the negative association by showing that the mean and median tornado counts of the LN seasons exceed those of the EN seasons.

The two most inactive LN DJF and MAM seasons are marked with blue diamonds in Figure 1, and the two most active ENs are marked with red diamonds. The two most active EN DJFs, which occurred in 2016 and 1983, are especially noticeable. There were 133 E(F)1+ tornadoes in both seasons. These counts are 183% greater than the mean number of tornadoes in EN DJFs and 359% greater than the median. The two most active EN MAM seasons were 1982 and 1983, when 340 and 251 E(F)1+ tornadoes occurred, respectively. MAM 1982 had approximately 79% more E(F)1+ tornadoes than the mean EN MAMs (76% greater than the median). MAM 1983 exceeded the mean and median counts of EN MAM seasons by approximately 30%.

The two DJF seasons with the fewest E(F)1+ tornadoes were 1985 and 1986, with 9 and 21 tornadoes, respectively. DJF 1985 had approximately 80% fewer tornadoes than the mean and median LN DJF, and DJF 1986 had 64% and 54% fewer than the mean and median counts, respectively. The least active LN MAM seasons were 1985 when 179 E(F)1+ tornadoes occurred and 2000 when 178 E(F)1+ occurred. These seasons had approximately 43% fewer tornadoes than the mean LN MAM tornado count and approximately 24% fewer than the median count.

As with most seasons, the tornadoes in the anomalous seasons were not uniformly distributed (Figure 2). In DJF 2016, for example, 50% of the tornadoes occurred on 4 days, each with 10+ tornadoes. There were 6 days in DJF 1983 with 10+ tornadoes that account for 68% of the tornadoes in that season. The tornadoes in the anomalous MAM seasons were spread over a larger number of days, but there were still clusters of activity. Twelve days in MAM 1982 and eight in MAM 1983 had 10+ tornadoes. The tornadoes on these days account for 60% and 52% of the seasonal count, respectively.

3.2. Role of GWO and MJO in the Anomalous Seasons

Previous studies illustrate that tornado activity is heightened during certain phases of the GWO and MJO. Gensini and Marinaro [21] reported that daily tornado anomalies in spring (March–June) are greatest on GWO phase 1 and 2 days when AAM is negative. Moore [22] similarly reported a tendency for tornado frequency to be greater in MAM seasons when GWO phase 2, 3, and 4 days are more common. Moore [22] also showed this to be true in DJF. Barrett and Gensini [23] reported that tornado days in April are most common on phase 6 and 8 days of the MJO and less common on phase 3, 4, and 7 days. They also reported that tornado days are most common with phases 5 and 8 and less common with phases 2 and 3 in May. Thompson and Roundy [20] reported that violent tornado outbreaks in MAM are most common on MJO phase 2 days and least common on phase 8 days. GWO and MJO vary on subseasonal timescales. They are, therefore, capable of modulating tornado activity within a given season and may provide insight into some of the subseasonal periods of suppressed and heightened tornado activity during these anomalous seasons (as seen in Figure 2).

AAM anomalies were positive throughout both of the anomalously active EN DJFs (Figures 3(a) and 3(c)). The mean AAM anomaly was 2.6 kg·m2·s−2 in DJF 1983 and 1.9 kg·m2·s−2 in DJF 2016 (Table 1). AAM anomalies were also positive throughout the active EN MAM 1983 season, when the mean was 1.8 kg·m2·s−2 (Figure 3(g), Table 1). Tornadoes in these seasons, therefore, occurred on GWO phase 5–8 days when AAM was anomalously high and the amplitude was most often >1. Positive anomalies were expected in these seasons because AAM tends to be heightened during the EN phase of ENSO [40], but it is unexpected, based on previous studies linking enhanced tornado activity to anomalously low AAM [21, 22], that all of the days with tornadoes in DJF 1983, DJF 2016, and MAM 1983 had anomalously high AAM. Despite the positive AAM anomalies throughout these seasons, the tendency of AAM was volatile, and tornadoes occurred on days when AAM tendency was increasing and decreasing.

In MAM 1982, there were periods when AAM was anomalously low, which resulted in a lower seasonal mean of 0.1 kg·m2·s−2 (Figure 3(e), Table 1). In this season, numerous tornadoes occurred on GWO phase 1–4 in addition to 5–8 days. However, there is not a significant difference in daily tornado counts across the phases of GWO in this season or any of the others (Table 2).

AAM is often relatively low during the LN phase of ENSO [40]. It is, therefore, not surprising that AAM was anomalously low throughout most of the inactive LN seasons (Figure 4). The mean AAM anomaly was −1.5·kg·m2·s−2 in MAM 2000 and −0.9 kg·m2·s−2 in DJF 1985 (Table 1). AAM fluctuated between negative and positive anomalies in MAM 1985 and DJF 1986, which led to higher mean values of 0.1 kg·m2·s−2 in each season (Figures 4(c) and 4(e); Table 2). Similar to the active EN seasons, there is not a significant difference in daily tornado counts across the GWO phases in any of the inactive LN seasons. Also, similar to the tornadoes in the active EN seasons, those in the inactive LN seasons occurred during periods of increasing and decreasing AAM (Figure 4).

MJO varied more than GWO throughout the anomalous seasons (Figures 3 and 4). The progression from phase 1 through 8 is apparent, with multiple oscillations in most seasons. Tornadoes concentrate on certain MJO phase days more so than with GWO phases, which led to significant differences in the mean number of tornadoes per day across the phases. In the two active EN DJFs, for example, the mean and mean ranks of the tornado counts were greatest with phases 1 and 6 of the MJO (Table 2). A Kruskal–Wallis test and subsequent post hoc comparisons indicate that the mean rank of phase 6 is significantly greater than the mean ranks of phases 4, 7, and 8; remaining comparisons yielded insignificant differences (see the subscript below Table 2). The percentage of days with tornadoes was also greatest with phase 6—tornadoes occurred on 10 of the 19 (53%) phase 6 days. In the two EN MAMs, the mean and mean ranks of the tornado counts are greatest with phase 4 (Table 2). The statistical tests indicated that the mean rank of phase 4 is significantly greater than that of phase 1 (see the subscript below Table 2). The percentage of days with tornadoes was also greatest with phase 4 (23 of 31 (74%) phase 4 days had tornadoes). The percentage of days with tornadoes was also high with phases 7 and 8 (69% and 68%, resp.). There were not any significant differences in daily tornado counts across the phases of MJO in the inactive LN seasons (Table 2).

3.3. Atmospheric Composites of the Anomalous Seasons

The seasonal composites of 300 mb geopotential height, 850 mb geopotential height and specific humidity, and CAPE are shown in Figures 58. Anomalously high geopotential heights at 300 mb were present across most of the contiguous US in DJF 2016 and across the north central and northeast US in DJF 1983 (Figures 5(a) and 5(b)). The presence of higher than normal heights is more similar to the LN composite reported by Allen et al. [12] rather than their EN composite. Upper-level ridges and troughs were present over the western and central US, respectively, in both of the active DJFs. Anomalously high heights were also present during the MAM 1982 seasons, but only over the southern and eastern US (Figure 6(a)). Low height anomalies were present over the contiguous US during MAM 1983 (Figure 6(b)), which is consistent with the EN composite reported by Allen et al. [12]. The patterns seen in the composites of low-level moisture are inconsistent across the active seasons. For example, anomalously high specific humidity was present in the southeastern US during MAM 1982, but the humidity in this region was anomalously low during MAM 1983 (Figures 6(c) and 6(d)). The humidity was also anomalously low in southeast Texas and Louisiana during both DJF seasons, but was higher to the north and east where most of the tornadoes occurred (Figures 5(c) and 5(d)). CAPE was near normal or anomalously high over most of the eastern US during the DJF 1983, DJF 2016, and MAM 1982 active EN seasons (Figures 5(e), 5(f) , and 6(e)). Similar to some of the patterns of 300 mb geopotential height, these patterns of elevated CAPE are more similar to the CAPE composites shown by Allen et al. [12] in association with the LN phase rather than the EN phase.

The patterns of 300 mb geopotential height varied between the two LN DJFs (Figures 7(a) and 7(b)). A dipole pattern was present in DJF 1985, with anomalously low heights spanning the northern Great Plains southwestward to the Southwest and high heights across the Pacific Northwest and Southeast US. The US was split in DJF 1986, with anomalously high heights across the western region and low heights across the eastern. The low-level humidity anomaly patterns also varied between the DJFs (Figures 7(c) and 7(d)). The upper-level height and low-level humidity patterns were more similar in the two inactive LN MAMs—near normal or anomalously high heights most of the US and anomalously humid conditions across most of its eastern half (Figures 8(a)8(d)). Anomalously low CAPE was present across the Southeast US in DJF 1985 when the fewest tornadoes occurred (Figure 7(e)). Anomalously high CAPE was present across the portions of the US with the most tornadoes in the remaining seasons (Figures 7(f), 8(e), and 8(f)).

Some of the seasonal composites illustrate reasonable patterns. Anomalously high CAPE is shown in both of the active EN DJFs, for example, in Figures 5(e) and 5(f). Most of the eastern US was also anomalously humid in these seasons (Figures 5(c) and 5(d)). Other seasonal composites are unexpected and do not capture the environment on the days when many tornadoes occurred. For example, the negative CAPE anomalies across the eastern US in MAM 1983 are unexpected, given that anomalously many tornadoes occurred in this season (Figure 6(f)). Examination of CAPE on 1-2 May 1983 and 18–20 May 1983, when 39 and 51 tornadoes occurred, respectively, illustrates that elevated CAPE spread from the Gulf of Mexico northward into the eastern US (Figure 9). The negative anomalies of low-level humidity across the south central US in DJF 2016 also did not represent well the environments that were present when many of the tornadoes occurred (Figure 5(d)). Closed shortwave troughs with upstream southerly moisture advection into the US were present on 23 December 2015 and 23-24 February 2016, when 18 and 36 tornadoes occurred, respectively (Figure 10). These cases illustrate that anomalously active seasons can have seemingly unfavorable seasonal composites.

4. Conclusions

Previous studies have established a relationship between tornado and ENSO in DJF and MAM, generally with more tornadoes during the LN phase and fewer during the EN phase [1214, 1719]. This study was focused on the seasons that do not fit this relationship—EN seasons with many tornadoes and LN seasons with few tornadoes. Specifically, the two EN DJFs and MAMs with the most tornadoes and the two LN DJFs and MAMs with the fewest tornadoes were described and analyzed. The most anomalous seasons were DJF 1983 and 2016, both of which were active EN seasons with 133 E(F)1+ tornadoes. They were, therefore, 183% (359%) above the mean (median) EN DJF. These seasons illustrate that large numbers of tornadoes are possible even during EN seasons when such large numbers might be unexpected. The other seasons were far less anomalous.

GWO does not explain the anomalous nature of the seasons. Climatological studies show that tornado activity in DJF and MAM tends to be heightened during GWO phases 1–4 when AAM is anomalously low [21, 22], but nearly all of the tornadoes in the active EN seasons occurred on GWO phase 5–8 days when AAM was anomalously high. Furthermore, daily tornado count did not significantly vary across the phases of GWO in any of the seasons. The concentration of tornadoes on high AAM days (GWO phase 5–8 days) during EN seasons, as suggested here, would undoubtedly weaken the statistical relationship between tornadoes and GWO that was reported by others, whereby tornadoes are most common on GWO phase 1–4 days when AAM is anomalously low [21, 22]. Analyzing the tornado-GWO relationship by ENSO phase might amend this relationship and provide additional insight into the interactions between tornado activity, ENSO, and GWO. Another consideration is that the GWO data used in this study are based on globally integrated AAM. As noted by Gensini and Allen [41], this may confound the results. Therefore, it would also be worthwhile to reassess the relationship between tornado activity and GWO with AAM calculated on a hemispheric or latitudinal basis.

Unlike with GWO, daily tornado count did significantly vary across the phases of MJO, but only in the active EN seasons. The mean and mean rank of the tornado counts was greatest with phase 4 in MAM. Clustering of tornadoes on MJO phase 4 days in MAM is unique because climatological studies [20, 23] that aggregated data over a larger number of seasons reported heightened tornado activity with phase 2, 5, and 7 of the MJO in boreal spring months. While only representative of MAM 1982 and 1983, the clustering of tornadoes on MJO phase 4 days in these EN seasons highlights the possibility for a large number of tornadoes when subseasonal climate oscillations are in unfavorable phases. The mean and mean rank of the tornado counts was also greatest with phase 6 of the MJO in the anomalously active EN DJFs. Additional study is needed to determine if the tendency for DJF tornadoes to cluster on phase 6 days is common across a larger number of seasons.

This study highlights EN seasons with anomalously many tornadoes and LN seasons with anomalously few tornadoes. It also highlights days in these seasons with many tornadoes when the subseasonal GWO and MJO were in phases previously unassociated with tornado activity. These results suggest that conditions at the synoptic scale and smaller contributed most to the anomalously active seasons (i.e., favorable conditions emerged despite the seemingly unfavorable states of the climate patterns in the active EN DJFs). Active tornado days in seasons when climate patterns are more favorable are also often driven by synoptic- and subsynoptic-scale conditions, but the tornado-favorable conditions are likely to emerge more often in these seasons. This study only considered the two most anomalously active EN MAMs and DJFs and the two most anomalously inactive LN MAMs and DJFs. Additional study of a larger number of seasons is needed. Study of additional EN seasons is particularly warranted to see if the weakened association between tornado activity and GWO that was seen here is common in other EN seasons. A better understanding of the ways in which ENSO, GWO, and MJO interact across seasons to influence tornado activity will likely improve seasonal and subseasonal outlooks.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.