Impacts of negative PNA and warm GoM on winter US tornado outbreaks
We first carried out a composite analysis to examine the roles of negative PNA and warm GoM in winter US tornado activity. As shown in Fig. 2a, the negative PNA is associated with the development of an anomalous atmospheric ridge around the eastern US. The associated anticyclonic circulation anomalies enhance the low-level southwesterly wind anomalies, increasing the flux of warm and moist low-level air from the GoM to the central and eastern US and LLWS in these regions (Fig. 2b). As a result, near-surface temperature 27–28 and atmospheric instability (e.g., LI) increase over the central and eastern US (Fig. 2b), leading to an increase in winter tornadogenesis in and around the Ohio Valley (Fig. 2c).
Interestingly, the warm GoM is also linked to a negative PNA-like pattern (Fig. 2d). An atmospheric ridge prevails along the eastern US, similar in location albeit weaker than the negative PNA composite (Fig. 2a), with associated anticyclonic circulation anomalies that enhance the low-level southwesterly wind anomalies. In this case, however, the low-level southwesterly anomalies strengthen predominantly over the southern US; thus, LLWS, the flux of warm and moist low-level air from the GoM, and negative LI are present along the Gulf states (i.e., Texas, Louisiana, Mississippi, and Alabama; Fig. 2e). This suggests that warm GoM is linked to enhanced tornado activity over the southern US and a portion of the southern Ohio Valley (Fig. 2f). The composites of SSTAs and geopotential height at 500 hPa derived from Community Earth System Model Large Ensemble Simulation (CESM-LENS) are largely consistent with the results from observations and reanalysis datasets (Fig. S1).
Relationship between PNA and GoM SSTAs
The negative PNA-like pattern in the warm GoM composite (Fig. 2d) suggests that GoM SSTAs can be forced by the PNA phase. Indeed, it is noted that the GoM is also warm in the negative PNA composite (Fig. 2a). To further examine the relationship between the PNA and GoM SSTAs, we show a histogram of GoM SSTAs during negative, positive, and neutral PNA phases in Fig. 3a. It clearly shows an inverse relationship between the PNA and GoM SSTAs. Specifically, during the negative PNA, the occurrence of warm GoM is more frequent than that during the positive PNA phase. For instance, the probability of extremely warm GoM (> 1.2K) increases by ~ 300% during the negative PNA compared to that during the positive PNA. During the positive PNA, on the other hand, the occurrence of cold GoM is more frequent than that during the negative PNA. This result supports the hypothesis that GoM SSTAs are physically linked to the PNA phases. The histogram derived from CESM-LENS also supports this result from observational and reanalysis datasets (Fig. S2a).
To further examine the physical relationship between the PNA and GoM SSTAs, we carried out a lead-lag correlation analysis between the PNA and GoM SSTAs. As shown in Fig. 3b, the negative PNA leads to warm GoM with a time-lag of approximately six days. A consistent result is found from CESM-LENS (Fig. S2b). The lagged relationship between the PNA and GoM SSTAs is largely due to PNA-induced surface heat flux anomalies over the GoM (Fig. 3b). Further analysis indicates that surface turbulent heat flux (i.e., latent and sensible) anomalies modulated by PNA-induced surface wind speed are the major contributors to the surface heat flux anomalies (Fig. S3).
Impact of long- versus short-lived PNA on winter US tornado outbreaks
The lead-lag correlation analysis (Fig. 3b) suggests that if the negative PNA and associated atmospheric ridge over the eastern US persist for six days or longer, GoM warms up significantly. In contrast, if the negative PNA dissipates within six days, GoM SSTAs are not significantly changed. This indicates that long-lived PNA helps warm up the GoM and thus provides a set of ideal environmental conditions for winter US tornado activity. To test this hypothesis, we carried out a composite analysis for the long-lived negative PNA (i.e., exceeding half a standard deviation and persisting for at least six days; Table S2) versus the short-lived negative PNA (i.e., exceeding half a standard deviation and dissipating in less than six days; Table S2).
As shown in Fig. 4a and 4e, 500 hPa geopotential anomalies exhibit a similar spatial pattern during the short-lived and long-lived negative PNA. However, the associated anomalous ridge over the US eastern seaboard is much more robust during the long-lived PNA. During the short-lived PNA, the atmospheric ridge extends meridionally across the eastern US, producing weak low-level southerly wind anomalies over the central US. This, in turn, increases LLWS, the flux of warm and moist air from the GoM, and decreases LI over the central US, producing favorable conditions for tornadogenesis in this region (Fig. 4b and c). In contrast, during the long-lived PNA, the anomalous ridge is tilted in such a way to produce strong low-level southwesterly anomalies across both the central and eastern US. This increases warm and moist air flux from the GoM and the eastern US coast into the US while increasing LLWS over the eastern US (Fig. 4f). This, in turn, increases atmospheric instability (i.e., LI) over the central and eastern US (Fig. 4g). During both the short- and long-lived PNA, surface heat flux is enhanced over the GoM as well as along the eastern US coast (Fig. 4c and g). However, since the reduction of near-surface wind speed is greater during the long-lived negative PNA than during the short-lived negative PNA (Fig. S4), the surface flux anomalies are much stronger and more extensive during the long-lived negative PNA (Fig. 4c and g). This explains why the long-lived PNA is more frequently accompanied by warm GoM. These large-scale environmental conditions associated with the short- and long-lived PNA are consistent with those derived from CESM-LENS (Fig. S5).
Largely consistent with the large-scale environmental patterns, winter US tornadogenesis is increased over a limited region in the central US (i.e., Arkansas, Missouri, and Illinois) during the short-lived negative PNA (Fig. 4d), while the area of increased winter tornadogenesis is much more extensive across the central and eastern US, including Mississippi, Tennessee, Kentucky, and Indiana, during the long-lived negative PNA (Fig. 4h). Indeed, the area of increased winter tornadogenesis during the long-lived negative PNA is surprisingly consistent with the area affected by the quad-state tornado outbreak (Fig. 4h).
The role of La Niña in the development of long- and short-lived PNA
Thus far, we have shown that a long-lived negative PNA pattern provides a set of atmosphere and ocean conditions ideal for fostering winter US tornado outbreaks similar to the quad-state tornado outbreak on 10–11 December 2021. Previous studies have shown that the PNA is a mode of internal variability and thus external forcing (e.g., ENSO) is not required for its development 25–26, 29. However, during La Niña, suppression of deep convection in the tropical Pacific promotes a negative PNA-like pattern through the atmospheric bridge 29–32. Consequently, cold SSTAs tend to appear in the tropical Pacific during both the short- and long-lived PNA events (Fig. 4a and 4e).
A related outstanding question is if the 2021-22 La Niña promoted the unusually prolonged negative PNA event that in turn fostered the quad-state tornado outbreak. While it is very difficult to find a definitive answer to this question, we may instead explore statistical relationships between the histogram of negative PNA days and ENSO conditions. As expected, the number of negative PNA days in a winter season (i.e., 92 days from November 1 to January 31) increases from 26.4 days per one El Niño season to 33.1 days per one La Niña season (a 25% increase). However, as shown in Fig. 5, the persistence of negative PNA is greater during La Niña than during El Niño. Specifically, about 39 out of 100 negative PNA events are long-live events (≥ 6 days) during El Niño conditions, whereas about 52 out of 100 negative PNA events are long-live events during La Niña conditions. Interestingly, unusually prolonged negative PNA events (≥ 20 days), as in the case of the quad-state tornado outbreak, occurred only during La Niña conditions. These results, also supported by the histogram derived from CESM-LENS (Fig. S6), suggest that the exceptionally long-lived negative PNA event during December 2, 2021 - January 4, 2022 (a total of 34 days) was promoted by the 2021-22 La Niña.
Concluding remarks
This study shows that the quad-state tornado outbreak on December 10–11, 2021 is closely tied to the long-lived negative PNA and associated warm GoM that produced a set of highly favorable environmental conditions for a winter tornado outbreak (i.e., increased LLWS and flux of warm and moist air from the GoM, and decreased LI) in and around the Ohio Valley. Further analysis indicates that the long-lived negative PNA produces a robust anomalous atmospheric ridge along the eastern US seaboard, producing anomalous low-level southwesterly wind and flux of warm and moist low-level air from the GoM into the central and eastern US. Through the associated surface turbulent heat flux, the long-lived negative PNA warms the GoM and thus further enhances the flux of warm and moist low-level air to the US. These atmospheric environments lead to a large increase in atmospheric instability over a broad region of the central and eastern US, increasing the likelihood of tornadogenesis in the region including the Ohio Valley where the quad-state tornado outbreak predominantly occurred. It is also shown that the occurrence of long-lived negative PNA events is more frequent during La Niña condition than during El Niño condition. This supports the idea that the 2021-22 La Niña promoted the exceptionally long-lived negative PNA event associated with the quad-state tornado outbreaks.
In contrast, the short-lived negative PNA produces a relatively weak and meridionally elongated atmospheric ridge over the eastern US. The atmospheric ridge produces southerly low-level wind anomalies, increasing LLWS and atmospheric instability over the central US. However, due to its short duration, the GoM SSTAs remain relatively neutral. During the short-lived negative PNA, tornadogenesis is increased over the central US. However, due to the relatively weak amplitude of the atmospheric instability (i.e., LI) and the near-neutral GoM SSTAs, winter tornadogenesis during the short-lived PNA is increased in a limited region in Arkansas, Missouri, and Illinois.
There remain several questions that are unexplored in this study. First, it is unclear through what processes the long-lived negative PNA is generated. Previous studies have shown that the negative PNA can be forced by the suppression of tropical Pacific convection (e.g., during La Nina) and associated extratropical Rossby wave propagation 29–31. The negative PNA can be also triggered by transient eddy vorticity fluxes over the extratropics 23, 28 and by extraction of kinetic energy from the zonally varying basic state 24, 29, 31, 33–35. Further analysis is required to better understand the potential link between these drivers and the persistence of the PNA. Another important question is if the frequency of long-lived negative PNA will change in the future. According to several previous studies 36–38, the amplitude of PNA is projected to increase partly due to the increasing variance of SSTAs in the tropical Pacific in the future, while the frequency of PNA is not expected to change significantly. However, it is still unclear whether the duration of the PNA will change in the future climate. Addressing the above questions could potentially help achieve a useful forecast capability for future winter US tornado outbreaks.