The second half of May 2019 was an unusually active period for tornadic thunderstorms across the U.S. Great Plains, Midwest, and lower Great Lakes. While this period typically coincides with the peak climatological frequency of tornadoes, preliminary reports of tornadoes were over triple the expected 30‐year average. Multiple‐day outbreaks of tornadoes are not unprecedented in the United States; however, this event was perhaps the first to be forecast at subseasonal lead times (3–4 weeks) by the Extended Range Tornado Activity Forecast team. This forecast of opportunity was driven, in part, by anomalous convective forcing in portions of the tropical Indian and Pacific Oceans, causing subsequent changes in Northern Hemisphere atmospheric angular momentum. This manuscript analyzes the evolution of hemispheric‐scale circulation features leading up to the event, examines teleconnection processes known to influence U.S. tornadoes, and provides insights into the forecast process at subseasonal lead times.

Nonlinear balance potential vorticity (PV) inversion is used to diagnose the sensitivity of the severe convective parameter space to the amplitude of a subsynoptic-scale PV anomaly on 13 March 1990, a day on which a significant tornado outbreak impacted the Great Plains. PV surgery is used to both amplify and remove the PV anomaly, and the contemporaneous impact on various convective parameters is subsequently quantified by using piecewise PV inversion to compute the changes in those parameters attributable to each PV alteration. It is found that amplifying the anomaly increases the CAPE by amounts typically ranging from 20% to 30% within the atmospheric columns experiencing the maximum PV increase. Ascent is increased slightly downshear of the PV anomaly, consistent with extant conceptual models governing synoptic-scale forcing for vertical motion. Amplifying the PV anomaly increases deep-layer shear over the southern half of the outbreak region and reduces storm-relative helicity over the northern half, primarily through changes in the estimated storm motion vector. Removing the anomaly produces complementary changes of the opposite sign. Thresholds of several commonly used convective parameters are chosen on the basis of prior empirical studies, and the horizontal displacement of these threshold contours produced by the PV alterations reveals that relatively modest subsynoptic-scale PV changes would not likely change the predominant convective mode during the Hesston outbreak.

Large-scale weather patterns favorable for tornado occurrence have been understood for many decades. Yet prediction of tornadoes, especially at extended lead periods of more than a few days, remains an arduous task, partly due to the space and time scales involved. Recent research has shown that tropical convection, sea surface temperatures, and the Earth-relative atmospheric angular momentum can induce jet stream configurations that may increase or decrease the probability of tornado frequency across the United States. Applying this recent theoretical work in practice, on 1 March 2015, the authors began the Extended-Range Tornado Activity Forecast (ERTAF) project, with the following goals: 1) to have a map room–style discussion of the anticipated atmospheric state in the 2–3-week lead window; 2) to predict categorical level of tornado activity in that lead window; and 3) to learn from the forecasts through experience by identifying strengths and weaknesses in the methods, as well as identifying any potential scientific knowledge gaps. Over the last five years, the authors have shown skill in predicting U.S. tornado activity two to three weeks in advance during boreal spring. Unsurprisingly, skill is shown to be greater for forecasts spanning week 2 versus week 3. This manuscript documents these forecasting efforts, provides verification statistics, and shares the challenges and lessons learned from predicting tornado activity on the subseasonal time scale.

A potential vorticity (PV) diagnostic framework is used to explore the sensitivity of the 3 May 1999 Oklahoma City tornado outbreak to the strength of a particular PV anomaly proximate to the geographical region experiencing the tornado outbreak. The results derived from the balanced PV diagnosis agree broadly with those obtained previously in a numerical simulation of the same event, while offering additional insight into the nature of the sensitivity. Similar to the findings of other cases, the balanced diagnosis demonstrates that intensifying (removing) the PV anomaly of interest increases (decreases) the balanced CAPE over the southwestern portion of the outbreak region, reduces (increases) the storm-relative helicity, and increases (reduces) ascent. The latter finding, coupled with the results of the modeling study, demonstrates that intensifying a PV anomaly proximate to an outbreak environment can increase the likelihood that more widespread and possibly less tornadic convection will ensue. The overall results of the balanced diagnosis complement those of other case studies, leading to the formulation of a conceptual model that broadly anticipates how the convective regime will respond to changes in intensity of upper-tropospheric weather features.

Idealized numerical experiments are conducted to understand the effect of upper-tropospheric potential vorticity (PV) anomalies on an environment conducive to severe weather. Anomalies are specified as a single isolated vortex, a string of vortices analogous to a negatively tilted trough, and a pair of string vortices analogous to a position error in a negatively tilted trough. The anomalies are placed adjacent to the tropopause along a strong upper-level jet at a time just prior to a major tornado outbreak and inverted using the nonlinear balance equations.In addition to the expected destabilization beneath and adjacent to a cyclonic PV anomaly, the spatial pattern of the inverted balanced streamfunction and height fields is distorted by the presence of the horizontal PV gradient along the upper-tropospheric jet stream. Streamfunction anomalies are elongated in the cross-jet direction, while height and temperature anomalies are elongated in the along-jet direction. The amplitude of the inverted fields, as well as the changes in CAPE associated with the inverted temperature perturbations, are linearly proportional to the amplitudes of the PV anomalies themselves, and the responses to complex PV perturbation structures are approximated by the sum of the responses to individual simple PV anomalies. This is true for the range of PV amplitudes tested, which was designed to mimic typical 6-h forecast or analysis errors and produced changes in CAPE beneath the trough of well over 100 J kg Ϫ1. Impacts on inverted fields are largest when the PV anomaly is on the anticyclonic shear side of the jet, where background PV is small, compared with the cyclonic shear side of the jet, where background PV is large.