1. Introduction
Tornadoes rated F3 or higher on the Fujita (1981) scale are rare in eastern New York and western New England. Since 1950, there have been only six occurrences of F3 tornadoes or stronger (Table 1), two of which were killer tornadoes. On 28 August 1973, an F4 tornado struck Columbia County, New York, and Berkshire County, Massachusetts, killing four and injuring 36. The same counties were struck again on 29 May 1995 by an F3 tornado producing three deaths and 29 injuries. On 31 May 1998, an F3 tornado formed just west of Mechanicville, Saratoga County, New York, and traveled over 50 km before reaching southwest Vermont as an F2 tornado (Fig. 1). Since significant tornadoes have the potential to produce severe damage and great loss of life, it is incumbent on the National Weather Service (NWS) to provide as much warning as possible for these storms. Although 68 people were injured in the Mechanicville tornado, accurate and timely severe weather outlooks, tornado watches, and tornado warnings, and rapid communication of hazardous weather information likely contributed to no lives being lost. Across the entire Northeast on 31 May, there were 5 fatalities and 127 injuries, associated with 32 tornadoes, 264 reports of wind damage, and 84 reports of 3/4-in. (1.9 cm) diameter or larger hail (Fig. 2a).
Although severe weather in the Northeast is infrequent compared to that in the Great Plains, it does occur as is evident from climatological studies (e.g., Brooks et al. 2003) and a review of research on severe weather in the Northeast (e.g., Johns and Dorr 1996). The unique terrain of the Northeast can also influence the observed distribution of severe weather (e.g., Riley and Bosart 1987; Wasula et al. 2002). A perspective on operational severe weather forecasting in the Northeast and an operational severe weather forecasting algorithm for the NWS Albany, New York, county warning area is presented in Maglaras and LaPenta (1997).
Since tornado outbreaks of this magnitude are relatively rare in the Northeast, it is important to examine the evolution of this outbreak, and in particular the Mechanicville tornado. Carefully reviewing past storms is a key step in improving a forecaster’s ability to diagnose severe weather potential and accurately forecast and warn for tornadoes and severe thunderstorms. This paper will examine the synoptic and mesoscale conditions that resulted in this widespread severe weather outbreak. Although the Mechanicville tornado developed within a synoptic-scale environment already favorable for tornadoes and severe thunderstorms, we will present some limited evidence that terrain-channeled low-level flow in north–south river valleys and the interaction of the isolated supercell responsible for the Mechanicville tornado with an overtaking squall line may have enhanced the likelihood of tornadogenesis. The paper is organized as follows. The synoptic overview appears in section 2. The mesoscale environment is described in section 3, while the details of the Mechanicville tornado appear in section 4. Speculations about the formation of the Mechanicville tornado are presented in section 5 with a summary in section 6.
2. Synoptic overview
Severe weather was widespread from the northern plains across the Great Lakes to New England during late May and early June 1998 (Fig. 2). The 200-hPa height pattern over North America during this period was dominated by above normal heights from north of Hawaii northward across the Gulf of Alaska, eastward toward Greenland, and over much of the southern half of the United States (Fig. 3). Below normal heights were found west of California and over Canada to the south of Hudson Bay. Given the strong trough (ridge) south of Hudson Bay (over the southern United States), the pattern shown in Fig. 3 was conducive to a strong westerly jet from the northern Rockies to New England. The severe weather threat was maximized where the ascent zones associated with disturbances embedded in the jet intercepted low-level baroclinic regions and warm, moist unstable air.
A 500-hPa closed low, destined to play a significant role in the severe weather shown in Fig. 2, moved southward off the west coast of North America, turned northeastward, caused weak cyclogenesis over Nevada at 0000 UTC 29 May, and reached central Idaho by 1200 UTC 30 May (Fig. 4). A weak 500-hPa trough to the southeast of this landfalling closed low formed over western Colorado at 1200 UTC 30 May, moved to central Nebraska 12 h later, and subsequently initiated weak cyclogenesis over northwest Kansas at 1200 UTC 30 May (a weak lee cyclone was present over the Texas and Oklahoma panhandles during the previous 24-h period).
The key player in the eventual 31 May severe weather outbreak over the Northeast was a 500-hPa vorticity center that was located over northern Canada at 1200 UTC 29 May (Fig. 4). This disturbance moved southward and then southeastward around the western periphery of a closed low situated over Hudson Bay (not shown). During the late afternoon and evening of 30 May, it helped produce severe weather over South Dakota, where the town of Spencer was devastated by an F4 tornado (NCDC 1998). A derecho evolved from this convection and moved eastward across Minnesota, Iowa, and Wisconsin at more than 30 m s−1, reaching lower Michigan by 0800 UTC 31 May. This derecho produced widespread wind damage (Fig. 2) in association with wind gusts of 30–40 m s−1 (NCDC 1998).
By 1200 UTC 31 May, the above disturbance was located over the upper peninsula of Michigan, where it began to phase with the weak disturbance moving northeastward across Colorado. In response, the surface cyclone over the northern Rockies accelerated eastward and merged with the weaker cyclone moving northeastward out of Kansas. Subsequently, the 500-hPa disturbance turned northeastward, stalled in the vicinity of James Bay, and was absorbed by the Hudson Bay closed low. As the associated surface cyclone moved northeastward to southeast of James Bay, it deepened 13 hPa in the 24-h period ending 0000 1 June. The severe weather over the northeastern United States on 31 May occurred in the warm sector of this deepening cyclone. Additional disturbances that rotated around the closed low also interacted with the strong upper-level jet (Fig. 3) and helped to produce two additional major severe weather outbreaks (29 May and 2 June) over the northeastern United States in addition to the 31 May storms (Figs. 2b and 2c).
A series of surface and upper-air maps were crafted from the Global Reanalysis Data produced by the National Oceanic and Aeronautics Administration’s (NOAA) National Centers for Environmental Prediction (NCEP) and the National Center for Atmospheric Research (NCAR) (Kalnay et al. 1996; Kistler et al. 2001) to provide a synoptic overview of this event. A surface cyclone located over South Dakota at 1200 UTC 30 May (Fig. 4) intensified as it moved east-northeast to a position north of New York and New England by 0000 UTC 1 June (Figs. 4, 5a–c). Model-diagnosed upward vertical motion at 700 hPa, from the Great Lakes northward into Canada at 1200 UTC 31 May, moved east across New York, New England, and Quebec during the next 12 h (Figs. 6d–f).
The observed 850-hPa winds at Detroit, Michigan (DTX), and Wilmington, Ohio (ILN), at 1200 UTC 31 May were 33 and 23 m s−1, respectively. This strong low-level flow produced significant warm-air advection across the Northeast (Figs. 5d–f). In addition, it helped to advect moisture rapidly into eastern New York during the morning. A tongue of high equivalent potential temperature (θe) was drawn northeastward toward the lower Great Lakes and then New York and New England, increasing the convective available potential energy (CAPE; Figs. 5d–f). The strong 850-hPa flow also contributed to strong shear and storm-relative helicity (SRH) in the lower troposphere (section 3). The severe storms that were initiated across the Northeast in the 6 h ending 0000 UTC 1 June occurred in a warm-sector environment of 700-hPa ascent and relatively dry air, a situation conducive to destabilization, with drier air aloft over high-θe air below (Figs. 6e,f). By 0000 UTC 1 June, a deepening 500-hPa trough and associated cyclonic vorticity axis had acquired a negative tilt (Figs. 7a–c). The surface cyclone had deepened to 985 hPa and had turned northward as the trailing cold front began to cross New York (Figs. 4, 5c).
At 250 hPa, a coupled jet pattern was evident across the Northeast (Figs. 7d–f). Divergence was concentrated in the left exit region of the trailing jet and the right entrance region of the leading jet, especially by 0000 UTC 1 June (Fig. 7f). The aforementioned derecho that affected the Great Lakes region originated hours earlier as isolated storms over South Dakota within the left exit region of the trailing 60 m s−1 250-hPa jet that was crossing extreme northern Colorado and southern Wyoming. This trailing jet shifted rapidly eastward to Iowa by 1200 UTC 31 May (Fig. 7a), at which time the remnant convection from the derecho over lower Michigan at 0800 UTC was situated near a region of 250-hPa divergence located between two coupled jets (Fig. 7d). By 0000 UTC 1 June, the 250-hPa divergence core was spread in an arc from eastern Pennsylvania to southern Quebec between the coupled jets and above the southern part of the axis of 700-hPa ascent (Figs. 6f, 7f). The strongest warm advection and strongest ascent were well north of the severe weather areas, but nonetheless, upward motions were sufficient to enhance severe storm potential, and allow parcels to reach the level of free convection and release the instability.
3. Mesoscale environment
a. Morning destabilization
The 1200 UTC 31 May sounding from Detroit, Michigan (DTX; Fig. 8a), was taken just to the west of the derecho that crossed the Great Lakes overnight and, thus, was representative of the immediate post-derecho mesoscale environment. A strong corridor of westerly winds (peak speeds > 35 m s−1) was found in conjunction with the deep mixed layer from 825 to 700 hPa that was possibly the signature of a descending rear-inflow jet behind the derecho (Weisman 1992). The DTX sounding also indicated that the air mass below the mixed layer had been stabilized, likely in response to rain-cooled air. Heating and moistening the low levels of the DTX sounding (Fig. 8a), to make it representative of the warmer and moister conditions observed in New York during the afternoon, would produce a very unstable air mass. Unfortunately, no representative sounding was available to confirm this hypothesis because there was no Buffalo, New York radiosonde launch for 1200 UTC 31 May, owing to equipment failure. Instead, Fig. 8b shows the 1200 and 1800 UTC 31 May soundings from Albany, New York (ALY). A stable pre–warm front environment at 1200 UTC warmed and moistened considerably by 1800 UTC. During this 6-h period, the surface-based CAPE (SBCAPE) increased from 0 to more than 750 J kg−1 as higher-θe air to the southwest began to arrive in conjunction with strengthening low-level southerly winds.
A comparison of the 1200 and 1800 UTC 31 May ALY hodographs (Fig. 9) shows that the 0–6-km above ground level (AGL) shear1 and SRH increased (shear from 24 to 29 m s−1 and SRH from about 250 to over 500 m2 s−2) in response to a strengthening 25 m s−1 west-southwesterly jet between 750 and 500 hPa and a strengthening of the low-level southerly flow below 850 hPa (note that the wind direction near 800 hPa in the 1800 UTC ALY sounding was likely erroneous; Fig. 8b). The end result was a hodograph indicative of favorable shear conditions for severe weather, by 1800 UTC 31 May.
Between 1500 and 2000 UTC, the air mass over New York and Pennsylvania rapidly destabilized (Figs. 10 –12). At 1500 UTC, most of the Hudson Valley region of New York and all of New England were under a south and southeasterly flow of relatively cool marine air, with surface dewpoints generally under 15°C (59°F; Fig. 10). A southwesterly flow across Pennsylvania and western-central New York marked the leading edge of the advancing warmer, more moist, and unstable air mass (Fig. 10). The eastward advance of the cloud shield across the upper Hudson Valley into western New England between 1200 and 1500 UTC and the clear area over eastern Pennsylvania, New Jersey, and the southwestern Catskills of New York at 1500 UTC show the eastward movement of the warm, humid, and unstable air mass (Figs. 10, 13a–c). Surface θe values exceeded 340 K by 1500 UTC across western New York and Pennsylvania, with an area in excess of 348 K over southeastern Pennsylvania, as solar heating combined with the transport of increasingly warm and moist air to destabilize the air mass [Fig. 10; surface-based lifted indexes (LIs) had decreased to −5 to −8 in this region].
Satellite visible imagery showed thick cloudiness with a large area of showers, and some embedded thunderstorms, from northern New York to northeastern Ohio by 1700 UTC (Fig. 13c). This area of clouds and precipitation was associated with the remnants of the derecho that crossed portions of the Great Lakes the previous night. As this area of convection moved eastward across New York, thunderstorms re-ignited between 1700 and 2000 UTC (Figs. 13c,d). The surface map for 1800 UTC (Fig. 11) showed the strengthened southerly flow up the Hudson Valley as higher-dewpoint air [∼18°C (64°F) ] reached the region. Immediately to the west in the clear region (Fig. 13b), θe values exceeded 340 K over a wide area. Surface-based LIs were between −6 and −10 in the clear air in response to continued heating and airmass destabilization.
b. Environment for supercell growth over complex terrain
As the afternoon of 31 May 1998 progressed, the atmosphere across eastern New York became increasingly favorable for severe weather. By 2000 UTC the warm and moist, unstable air had overspread the Hudson Valley and much of interior southern New England with θe values > 348 K being common (Fig. 12). At Glens Falls (KGFL) and Albany (KALB), New York, sustained winds in the southerly flow in the Hudson Valley were exceeding 10 m s−1. Surface temperatures in the clear air southeast of the developing thunderstorms were in the 27°–32°C (81°–90°F) range, with dewpoints of 18°–20°C (64°–68°F). The 1800 UTC radiosonde observation at ALY (Fig. 8b), located about 6 km south of KALB, was modified for the observed KALB surface temperature and dewpoint temperature at 2000 UTC. The SBCAPE [mean layer CAPE (MLCAPE)] increased to about 2000 J kg−1 (1600 J kg−1). The mean layer lifted condensation level (MLLCL, lifted condensation level based on lifting a parcel with conditions averaged over the lowest 100 hPa) was about 1200 m. Thompson et al. (2003) showed median MLLCLs for significant tornado, weak tornado, and non–tornadic supercell cases, of about 1000, 1200, and 1350 m, respectively.
The shear magnitude was 29 m s−1 in the 0–6-km AGL level. Observations and modeling studies (e.g., Weisman and Rotunno 2000) indicate that vertical wind shear greater than 20–25 m s−1 is favorable for the development of long-lived convective structures such as supercells. Based on observed storm motion of 260° at 21 m s−1, the SRH was greater than 500 m2 s−2. Davies-Jones et al. (1990) indicated this SRH value would be large enough to support significant (F2 and greater) tornadoes. Recent studies using Rapid Update Cycle model generated proximity soundings (Thompson et al. 2003) and observed proximity soundings (Rasmussen 2003) suggest that shear and SRH in the 0–1-km AGL layer are better at discriminating supercells that produce significant tornadoes from nontornadic supercells than deep-layer (0–6-km shear and 0–3-km SRH) measures of these parameters. The SRH (shear) in the 0–1-km AGL layer was about 140 m2 s−2 (14 m s−2). In Thompson et al. (2003) more than 75% of the significantly tornadic supercells had 0–1-km SRH > 75 m2 s−2 while 70% of the nontornadic supercells occurred with lesser values. The 0–1-km SRH was well above the median value (89 m2 s−2) for tornadic storms in Rasmussen (2003). Finally, based on Thompson et al. (2003) and Rasmussen (2003), the combination of MLCAPE and 0–1-km SRH pointed to the possibility of significant tornadoes.
Although conditions in the warm sector of this storm were favorable for the development of severe storms and tornadoes, terrain may have played a role in enhancing the tornadic potential in the immediate Hudson Valley. A number of significant physiographic features are present in eastern New York and western New England. To the west of the north–south-oriented Hudson Valley are the Catskill and Adirondack Mountains (Fig. 1). These two mountain ranges are separated by the east–west-oriented Mohawk Valley. To the east of the Hudson Valley in western New England are the Green (Vermont) and Berkshire (Massachusetts) Mountains. On this day, terrain channeling of the southerly low-level flow within the Hudson Valley may have enhanced the transport of warm, moist air northward. At KALB, the dewpoint increased steadily through the day, and the wind remained out of the south until convection arrived shortly after 2000 UTC (Fig. 14a). This southerly flow transported high-θe air northward. The 1900 UTC θe analysis2 shows an axis of higher-θe air extending northward up the Hudson Valley (Fig. 15). The radar reflectivity pattern has been included to help delineate areas of rain-cooled air. This pattern of higher-θe (instability; Fig. 16) air in the Hudson Valley with a south or southwest low-level wind flow is common, and is also seen in higher-resolution operational numerical models.
The near surface flow west of the Hudson Valley was generally out of the south-southwest. This south-southwest flow was channeled into a more southerly direction at low levels within the immediate Hudson Valley. Binghamton, New York (KBGM), consistently had a more westerly component to the surface wind than Hudson Valley locations [KALB and Schenectady, New York (SCH); Fig. 14b]. At Utica, New York (KUCA), the surface wind was from the southeast (up the Mohawk Valley; see Fig. 1) prior to 1500 UTC. Once in the warm sector (after 1600 UTC), however, the wind at KUCA showed a greater westerly component than at Hudson Valley locations (Fig. 14b). The channeling of the surface wind in the Hudson Valley to a more southerly direction resulted in an increase in the shear and SRH. The channeling of a southwest flow to a southerly flow within the Hudson Valley is a common occurrence, as indicated by the wind rose for KALB in Fig. 14c taken from Wasula et al. (2002). The wind rose for KALB at 850 hPa in Wasula et al. (2002); not shown) does not show a southerly directional bias, which illustrates the decreasing influence of the terrain with height on wind direction.
Figure 17 shows estimated hodographs based on the 1800 UTC ALY sounding for a location in the Hudson Valley, and for a location west of the Hudson Valley. The lower 2 km of the observed hodograph were modified using the 1900 UTC observed surface winds at KALB and KBGM, and applying a nonlinear Ekman spiral approach to interpolate the winds between the surface and 2000 m. The channeling of the wind in the Hudson Valley hodograph resulted in an increase in 0–3-km (0–1 km) SRH from about 400 (300) to about 530 (420) m2 s−2. There was also an increase in shear in the lowest kilometer. There is no way to know the actual hodographs at this time, and a different methodology for estimating the hodograph could produce a different result. To test the validity of the estimated hodographs, the 1200 UTC 31 May Eta Model soundings for 1900 UTC for KALB (in the Hudson Valley) and KBGM (west of the Hudson Valley) were compared to the estimated hodographs (Fig. 17b). At 1900 UTC, both KBGM and KALB were in the warm sector of the low pressure area, and ahead of the thunderstorms approaching from the west (Figs. 14b and 15). The winds in the two model sounding hodographs were similar at 2 and 3 km AGL, but were markedly different below. The KALB model sounding was similar to the Hudson Valley hodograph, while the KBGM sounding showed somewhat less curvature (Fig. 17a). The model sounding winds showed a more westerly component at 2 km AGL than in the observed winds. The hodographs indicate there might have been enhancement of SRH and shear due to terrain effects.
A 6-h NCEP Eta Model forecast (the model was initialized at 1500 UTC with 32-km grid spacing) of 975-hPa winds valid at 2100 UTC (Fig. 18) showed evidence of terrain channeling. The winds in the Hudson Valley, especially near and just south of Mechanicville, were more southerly than in surrounding areas to the west and to the east. This helped maximize the 0–6-km shear (29 m s−1) just south of Mechanicville (Fig. 18). While mesoscale conditions supported the possibility of tornadoes (e.g., Rasmussen 2003; Johns and Doswell 1992; Davies-Jones et al. 1990), a storm moving from the northern Catskills into the Hudson Valley would encounter a somewhat more unstable air mass with greater SRH and shear in the lower troposphere, due to terrain channeling.
4. Evolution of the Mechanicville tornado
a. Radar analysis
The remnants of the derecho that crossed Michigan during the early morning hours of 31 May reached western New York about 1600 UTC. The convection then intensified as it moved east, producing locally severe weather as it reached central New York by 1824 UTC (Figs. 2 and 19a), about 2 h prior to the Mechanicville tornado (2022 UTC). As the line of thunderstorms passed KUCA, it produced a 7°C (13°F) drop in temperature between 1900 and 2000 UTC. At 1824 UTC, a small cluster of isolated storms developed about 55 km to the east of the line. One of the storms in this cluster intensified by 1906 UTC as it moved east-northeast (Fig. 19b). By 1942 UTC, the advancing and intensifying line began to bow down the Mohawk Valley as it was catching up with the leading supercell storm (Fig. 19c). This leading storm, only about 30 km ahead of the line racing down the Mohawk Valley, had become severe (downed trees were reported) as it approached the Hudson Valley north of ALY (Fig. 19c).
North of the apex of the bowing squall line, there was evidence of a line echo wave pattern (LEWP; Fig. 20a). A rotational velocity couplet (13 m s−1 at 1946 UTC) along the Fulton–Montgomery County line (upper circle in Fig. 20b) was found near the northern end of the advancing squall line, in a location where a northern bookend vortex would be found. A separate velocity couplet (13 m s−1 at 1946 UTC) was also apparent just south of the Fulton–Montgomery County line (lower circle in Fig. 20b). Two volume scans later, at 1957 UTC, the rotational velocity associated with the mesoscale vortex had decreased (Figs. 20c,d). Just to the east, the leading supercell intensified between 1957 and 2012 UTC (see discussion below) as it approached the confluence of the Mohawk and Hudson Valleys (Figs. 19d, 20a,c). This intensification may be related to the presence of a slightly more unstable and more highly sheared environment in the Hudson Valley and/or the ongoing interaction of the overtaking squall line with the leading supercell (Figs. 14a, 15 –18, 19d).
A tornado developed at 2022 UTC, about the time the outflow from the line of storms to the west reached the leading severe storm (Fig. 19d). This outflow was likely channeled east-southeastward down the Mohawk Valley. KENX base velocity data at 1957 UTC (25 min before tornadogenesis) indicated an area of higher winds in this storm-initiated outflow moving down the Mohawk Valley (Fig. 21). At 2002 UTC, 20 min prior to the tornado, weak convergence was indicated in the SRM data about 1 km above mean sea level (MSL) in the KENX 0.5° beam angle volume scan (not shown). Above 1 km, cyclonic rotation increased with height, reaching 21 m s−1 at 3–4 km MSL. By 2012 UTC, strong rotation (21–23 m s−1) had descended toward the surface, but at the lowest elevation scan (about 1 km MSL), convergent rotation was present. At that time, the outflow from the line of severe storms approaching from the west had reached the storm. At 2017 UTC, strong rotation reached the lowest elevation scan (0.5°), but inbound and outbound velocity maxima were still separated by about 2 km. The first tornado damage occurred at 2022 UTC and, at that time, the radar showed a maximum gate-to-gate rotational velocity of 28 m s−1 at 0.5° (Fig. 22a). The Weather Surveillance Radar-1988 Doppler (WSR-88D) first identified a tornado vortex signature (TVS) based on the volume scan beginning at 2022 UTC, although it should be noted that the TVS information was not available to forecasters until several minutes after the completion of the volume scan, or at about 2030 UTC.
At 2022 UTC there was strong gate-to-gate rotation at the lowest elevation scan (about 1 km MSL), which was likely indicative of the low-level tornadic mesocyclone (Fig. 22b). Also apparent is the thunderstorm mesocyclone with inbound and outbound velocity maxima separated by about 8 km. The tornado reached its greatest intensity shortly after formation in Saratoga County, New York. It moved east across northern Rensselaer County, New York, before dissipating in Bennington County, Vermont (F2 intensity), with a total pathlength of 48 km (Fig. 1).
b. Lightning analysis
Cloud-to-ground (CG) lightning flash rates and distributions were also examined to document the interaction of the overtaking squall line with the isolated supercell that eventually produced the Mechanicville tornado. Our analysis was motivated by previous CG lightning studies (e.g., Kane 1991; Seimon 1993; Doswell and Brooks 1993, MacGorman and Burgess 1994; Galarneau et al. 2000; Gilmore and Wicker 2002) that suggested that some tornadic supercell thunderstorms can have distinct lightning flash rates and polarity characteristics. We used CG flash data obtained from the National Lightning Detection Network (NLDN) to construct a time series of CG flashes in the advancing squall line and the isolated supercell beginning approximately 2 h prior to tornado formation at 2022 UTC 31 May. We also mapped the storm-relative distribution of CG flashes in 10-min periods, beginning at 2002 UTC, in an effort to identify possible spatial variations in CG flash activity during the interaction of the squall line and the isolated supercell.
In order to construct Fig. 23a, CG lightning flash data obtained from the NLDN were plotted on a topographic map. These CG flashes were assigned to either the squall line or the isolated supercell. We estimate that the subjective tabulation procedure had a small margin of error (±1–2 flashes per count) prior to squall line and supercell merger, and a larger margin of error (±5–10 flashes per count) subsequently. The separation distance between the squall line and supercell was determined using the KENX Doppler radar reflectivity images and was measured from the core of the highest reflectivities in the squall line to the supercell mesocyclone. The 10-min distributions of CG lightning flashes shown in Fig. 23b were plotted onto a horizontal map relative to the position of the KENX Doppler radar-indicated mesocyclone prior to tornadogenesis. The analysis domain was 70 km × 80 km (5600 km2) and was centered on the mesocyclone (or tornado). The number of CG lightning flashes, denoted in the lower-right corner of each panel, was counted objectively. The CG flash numbers were not normalized because the relative length of the squall line within the analysis domain did not change during the 40-min period shown in Fig. 23b (i.e., the overall area of the horizontal map covered by electrically active convection changed little with time).
The results of our CG flash analysis are shown in Fig. 23. Until about 1855 UTC, the CG flash rate in the isolated supercell remained under 10 flashes per 5 min, during which time the squall line remained 40–45 km to the west of the supercell (Fig. 23a). Subsequent to 1855 UTC, the squall line began to accelerate eastward relative to the isolated supercell, closing to within 10 km by 2005 UTC. Over this 80-min period, the supercell flash rate first increased slowly to 10–15 CG flashes per 5 min period, and then more rapidly to more than 20 CG flashes per 5 min after 1950 UTC. Beginning at 2010 UTC, the CG flash rate increased to 20 flashes per 2 min (note time interval reduction beginning at 2010 UTC) through 2014 UTC, during which time the advancing squall line reached to within 6–7 km of the supercell. A brief decrease in the CG flash rate from 2014 to 2018 UTC was followed by an increase to almost 25 flashes per 2 min by 2022 UTC. Subsequent to tornadogenesis at 2022 UTC, the CG flash rate continued to increase to more than 35 flashes per 2 min. The percentage of positive flashes generally remained under 10% (except 12% near 1905 UTC), but was briefly near zero at the time of tornadogenesis.
A clear separation can be seen between the CG flashes in the overtaking squall line and the leading isolated supercell between 2002 and 2012 UTC (Fig. 23b). The few positive CG flashes were concentrated mostly along the northern half of the squall line. By 2022 UTC (tornado formation), the overall number of CG flashes in 10 min had increased from 253 to 461. This increase in CG flash rate occurred as the squall line overtook the supercell (it was still possible to distinguish the CG activity with each feature). Between 2022 and 2032 UTC, the overall number of CG flashes increased further to 664 but it was more difficult to distinguish between squall line- and supercell-related CG flashes, as merger of the two systems had occurred. Likewise, the number and percentage of positive CG flashes increased through 2042 UTC, particularly in close proximity (within 10–15 km) of the tornado. The CG flashes were concentrated mostly to the left of the storm track, primarily in the forward-flank downdraft region of the supercell prior to tornadogenesis (Lemon and Doswell 1979) and merger with the squall line. Subsequent to tornadogenesis and merger, the CG flashes were concentrated in the left-rear and right-rear quadrants relative to the storm track. Overall, the results from Fig. 23 suggest that as the squall line bowed down the Mohawk Valley and overtook the leading supercell, there was a significant increase in the CG flash rate, which continued beyond 2022 UTC.
5. Speculations about the genesis of the Mechanicville tornado
While most supercells produce severe weather, less than half produce tornadoes (e.g., Thompson et al. 2003; Brooks et al. 1994). Several studies (e.g., Thompson et al. 2003; Rasmussen 2003; Brooks et al. 1994) have attempted to identify meteorological conditions that differentiate environments that produce tornadic supercells and nontornadic supercells. Although atmospheric conditions on 31 May 1998 favored the development of supercells and tornadoes, the critical question is what physical mechanisms may have contributed to the actual tornadogenesis process. We have attempted to address this question through an examination of the available observations. Our results indicate terrain-channeled southerly flow existed in the Hudson Valley that was backed relative to the more southwesterly flow to the west of the valley. This contributed to increased SRH and shear in the Hudson Valley. The highest peaks of the Catskill Mountains rise to a little above 1200 m MSL, indicating the impact of the terrain-channeled flow will be most significant in the critical 0–1-km layer. Increased SRH and shear in the lowest 1 km above the ground have been shown to be favorable for tornadogenesis in idealized modeling studies (e.g., Wicker 2000; Wilhelmson and Wicker 2001), and in studies using model-generated (Thompson et al 2003) and observed soundings (Rasmussen 2003) near supercells. As the supercell that produced the Mechanicville tornado reached the juxtaposition of the western edge of the Hudson Valley with the eastern end of the Mohawk Valley, it also may have been influenced by the higher-θe air being advected northward up the valley.
Another factor that may have played a role in tornadogenesis was the interaction of an overtaking squall line with a leading isolated supercell. The storm that produced the Mechanicville tornado developed over 50 km ahead of a line of severe storms to the west (1824 UTC). Just prior to the tornado (2022 UTC), the outflow from the line of storms to the west reached the supercell. The timing of the merger suggests the possibility that it could have contributed to the formation of the Mechanicville tornado, or at least the intensification of the supercell that produced it.
Przybylinski (1995) examined seven cases where isolated convective elements developed 50–80 km ahead of bowing lines. The lead cell intensified as the line approached and, in over half of the cases, tornadoes were observed when the merger occurred. One case produced several strong tornadoes. Goodman and Knupp (1993), Bullas and Wallace (1988), Hamilton (1969), and Cook (1961) have documented cases where tornadoes formed after this type of interaction. Goodman and Knupp (1993) examined the Huntsville, Alabama, tornado (F4) of 15 November 1989. The Edmonton, Alberta, Canada, tornado (F4) on 31 July 1987 developed in a similar fashion (Bullas and Wallace 1988). Hamilton (1969) summarized the work of Cook (1961) concerning cases where cells developed ahead of, and eventually merged with, a line of thunderstorms. Tornadoes occurred in 8 of 11 cases where the advance cell was intense. Other studies have also shown the importance of boundary interaction in tornado formation, although in these studies squall line–supercell interactions were not present in all cases. In a study of 86 tornadoes in the northeastern United States, LaPenta et al. (2000) indicated radar-detectable boundaries were present in about a fifth of the cases. Nearly 70% of significant tornadoes in the VORTEX-95 experiment (Markowski et al. 1998) occurred near preexisting mesoscale boundaries. It should be noted that unlike the Mechanicville event, the VORTEX cases involved old outflow boundaries, and the air on the cool side of these boundaries had been substantially modified. All of the studies cited point to the complexities involved in the tornado formation process and indicate that boundaries can be an important source of environmental heterogeneity.
A different perspective on the possible interaction of a squall line and an isolated supercell comes from an idealized simulation of a quasi-linear convective system by Weisman and Davis (1998). In an experiment with Coriolis effects included, they showed that well-defined backed flow was prevalent at the 2-km level ahead of a northern bookend vortex 4.5 and 6 h into the simulation (their Fig. 4). This result, repeated in Weisman (2001, his Fig. 22), prompts us to ask if such a region of low-level backed flow ahead of a northern bookend vortex were to begin overtaking an isolated supercell downstream, that the immediate supercell environment might be especially favorable for tornadogenesis? A future idealized simulation of a squall line with a northern bookend vortex overtaking a leading isolated supercell, modeled after the Mechanicville tornado case, might prove to be enlightening. A schematic synopsis of how the overtaking northern bookend vortex and squall line might have interacted with the leading isolated supercell and low-level terrain-channeled southerly flow in the Hudson Valley is offered in Fig. 24.
There are some similarities between the Mechanicville tornado event and the 29 May 1995 F3 tornado in Great Barrington, Massachusetts (GBR), as described by Bosart et al. (2004). Both events occurred on days when the environment over the northeastern United States was favorable for the development of supercells and tornadoes. In addition, both events were associated with terrain-channeled low-level flow in the north–south-oriented Hudson Valley and external forcing may have played a role in the evolution of the tornadoes. The possible importance of external forcing on tornadogenesis in both the Mechanicville and GBR events should be tested by idealized simulations as noted above.
6. Summary
During late May and early June 1998, a large-scale circulation pattern evolved into a configuration that was unusually favorable for widespread severe weather from the northern plains eastward across the Great Lakes into New England. A series of disturbances, embedded in a strong upper-level flow, helped to initiate episodic severe weather. On 31 May 1998, one episode resulted in 32 tornadoes over the Northeast, one of which was the F3 Mechanicville, New York, tornado and was the focus of this paper. This tornado was examined in detail because it was one of four F3 events in the Northeast, it caused significant damage and many injuries, and storms of its severity are rare in the Northeast. The 31 May disturbance was associated with a coupled jet structure at 250 hPa, with wind speeds exceeding 65 m s−1. A strong low-level jet (33 m s−1) advected warm, moist air into the Northeast; helped to destabilize the atmosphere (CAPE ∼2000 J kg−1); and contributed to high 0–6-km shear (>25 m s−1) and SRH (>450 m2 s−2). The Mechanicville event also bore some similarity to the westerly flow tornado cases defined in the Johns and Dorr (1996) study of strong and violent tornadoes in New England and eastern New York.
Forecast issues were also raised by the Mechanicville tornado event. Recognizing that synoptic-scale conditions were favorable for severe thunderstorms and tornadoes on 31 May 1998 was the first step in alerting the public to the potential for dangerous weather. The Storm Prediction Center had indicated a high risk of severe thunderstorms and tornadoes for parts of the Northeast. A local study (Maglaras and LaPenta 1997) supported the tornadic threat, and strongly worded statements were issued alerting the public to the dangerous situation. Previous research on the GBR tornado event (Bosart et al. 2004), a storm similar to the Mechanicville tornado, provided meteorologists an understanding of the potential role of terrain and mesoscale processes in supercell intensification. It also provided experience in radar-based detection of tornado formation. Being aware of the mesoscale environment, and the role terrain can play in modifying it, can be important in maximizing warning lead time. This is especially true in locations where tornadoes are infrequent. A tornado watch was issued nearly 2 h prior to the Mechanicville tornado and a tornado warning was issued with 39-min lead time. The large lead time can be attributed to a number of factors. The city of Mechanicville lies in the extreme eastern portion of Saratoga County, New York. As the storm moved east and neared the western portion of the county, about 40 km to the west of Mechanicville, a warning decision had to be made. Radar and spotter reports indicated that at least a severe thunderstorm warning was necessary. While the radar did not indicate a tornado at that time, a tornado warning was issued based on a number of factors. A tornado watch was in effect and the SPC had indicated a high risk of severe storms. The storm was already severe and a spotter reported a tornado about 9 km west of the Saratoga County line. The KENX radar is located at an elevation of about 550 m MSL and the radar beam can overshoot shallow features in valleys, so it was difficult to rule out this report. Also, experience and previous research (Bosart et al. 2004) indicated the storm was entering an area of somewhat increased instability and shear. The storm continued to intensify as it crossed the county line and as radar indicated the high probability of a tornado, a severe weather statement was issued refining the tornado track and indicating that a “dangerous tornado was moving across Saratoga County. . . .” Tornadoes rated F3 or higher in upstate New York are relatively rare. Although 68 people were injured in the Mechanicville tornado, accurate and timely outlooks, tornado watches and tornado warnings, and rapid communication of hazardous weather information likely contributed to no lives being lost.
Acknowledgments
The authors are grateful to COMET, which provided support for research (UCAR COMET Partners Grant S9773888) that was extremely valuable in supporting the NWS’s timely warnings for the Mechanicville tornado, and was the impetus for the continued research on tornadoes in the Northeast. Research on this storm has continued as parts of COMET Cooperative Grant UCAR-09915806 and NOAA Grant 1007941-1-012365 awarded to The University at Albany as part of the CSTAR program. The authors would also like to thank the National Weather Service for supporting the cost of publication. Anton Seimon is thanked for making critical comments on the manuscript. The authors would also like to thank Kristen Corbosiero for helping with the figures, Celeste Iovinella for helping to put the manuscript in final form, and Heather Hauser for her initial review of the manuscript. Finally, the authors would like to express their gratitude to the anonymous reviewers for many constructive suggestions for improving the manuscript.
REFERENCES
Bosart, L F., LaPenta K D. , Seimon A. , and Dickinson M. , 2004: Terrain-influenced tornadogenesis in the northeastern United States: An examination of the 29 May 1995 Great Barrington, Massachusetts, tornado. Preprints, 22d Conf. on Severe Local Storms, Hyannis, MA, Amer. Meteor. Soc., CD-ROM, 14.5.
Brooks, H E., Doswell C A. III, and Cooper J. , 1994: On the environment of tornadic and nontornadic mesocyclones. Wea. Forecasting, 9 , 606–618.
Brooks, H E., Doswell C A. III, and Kay M P. , 2003: Climatological estimates of local daily tornado probability. Wea. Forecasting, 18 , 626–640.
Bullas, J M., and Wallace A F. , 1988: The Edmonton tornado, July 31, 1987. Preprints, 15th Conf. on Severe Local Storms, Baltimore, MD, Amer. Meteor. Soc., 438–443.
Cook, B J., 1961: Some radar LEWP observations and associated severe weather. Proc. of the Ninth Radar Conf., Boston, MA, Amer. Meteor. Soc., 181–185.
Davies-Jones, R P., Burgess D. , and Foster M. , 1990: Test of helicity as a tornado forecast parameter. Preprints, 16th Conf. on Severe Local Storms, Kananaskis Park, AB, Canada, Amer. Meteor. Soc., 588–592.
Doswell C A. III, , and Brooks H E. , 1993: Comments on “Anomalous cloud-to-ground lightning in an F5-tornado-producing supercell thunderstorm on 28 August 1990.”. Bull. Amer. Meteor. Soc., 74 , 2213–2218.
Fujita, T T., 1981: Tornadoes and downbursts in the context of generalized planetary scales. J. Atmos. Sci., 38 , 1511–1534.
Galarneau T J. Jr, , Honikman S F. , Cacciola A C. , Bosart L F. , LaPenta K D. , and Quinlan J S. , 2000: Lightning in tornadic thunderstorms in the northeastern United States. Preprints, 20th Conf. on Severe Local Storms, Orlando, FL, Amer. Meteor. Soc., 108–109.
Gilmore, M S., and Wicker L J. , 2002: Influences of the local environment on supercell cloud-to-ground lightning, radar characteristics, and severe weather on 2 June 1995. Mon. Wea. Rev., 130 , 2349–2372.
Goodman, S J., and Knupp K R. , 1993: Tornadogenesis via squall line and supercell interaction: The November 15, 1989, Huntsville, Alabama, tornado. The Tornado: Its Structure, Dynamics, Prediction, and Hazards, Geophys. Monogr., No. 79, Amer. Geophys. Union, 183–200.
Hamilton, R E., 1969: A review of use of radar in detection of tornadoes and hail. Tech. Memo. WBTM-ER-34, Weather Bureau Eastern Region Headquarters, Garden City, NY, 64 pp.
Hart, J A., and Korotky W. , 1991: The SHARP workstation v1.50 users guide. National Weather Service, 30 pp. [Available from Scientific Services Division, NWS Eastern Region Headquarters, 630 Johnson Ave., Bohemia, NY 11716.].
Hart, J A., and Janish P. , cited. 1999: SeverePlot v2.0, historical severe weather report database. Version 2.0. NWS Storm Prediction Center, Norman, OK. [Available online at http://www.spc.noaa.gov/software/svrplot2/.].
Johns, R H., and Doswell C A. III, 1992: Severe local storms forecasting. Wea. Forecasting, 7 , 588–612.
Johns, R H., and Dorr R A. Jr., 1996: Some meteorological aspects of strong and violent tornado episodes in New England and eastern New York. Natl. Wea. Dig., 20 , 4,. 2–12.
Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-Year Reanalysis Project. Bull. Amer. Meteor. Soc., 77 , 437–471.
Kane, R J., 1991: Correlating lightning to severe local storms in the northeastern United States. Wea. Forecasting, 6 , 3–12.
Kistler, R., and Coauthors, 2001: The NCEP–NCAR 50-Year Reanalysis: Monthly means CD-ROM and documentation. Bull. Amer. Meteor. Soc., 82 , 247–268.
LaPenta, K D., Maglaras G M. , Quinlan J S. , Johnson H J. , Bosart L F. , and Galarneau T J. Jr., 2000: Radar observations of northeastern United States tornadoes. Preprints, 20th Conf. on Severe Local Storms, Orlando, FL, Amer. Meteor. Soc., 356–359.
Lemon, L R., and Doswell C A. III, 1979: Severe thunderstorm evolution and mesocyclone structure as related to tornadogenesis. Mon. Wea. Rev., 107 , 1184–1197.
MacGorman, D R., and Burgess D W. , 1994: Positive cloud-to-ground lightning in tornadic storms and hailstorms. Mon. Wea. Rev., 122 , 1671–1697.
Maglaras, G J., and LaPenta K D. , 1997: Development of a forecast equation to predict the severity of thunderstorm events in New York State. Natl. Wea. Dig., 21 , 3,. 3–9.
Markowski, P M., Rasmussen E N. , and Straka J M. , 1998: The occurrence of tornadoes in supercells interacting with boundaries during VORTEX-95. Wea. Forecasting, 13 , 852–859.
NCDC, 1998: Storm Data. Vol. 40, No. 5, 408 pp. [Available from National Climatic Data Center, 151 Patton Ave., Asheville, NC 28801.].
Przybylinski, R W., 1995: The bow echo: Observations, numerical simulations, and severe weather detection methods. Wea. Forecasting, 10 , 203–218.
Rasmussen, E N., 2003: Refined supercell and tornado forecast parameters. Wea. Forecasting, 18 , 530–535.
Riley, G T., and Bosart L F. , 1987: The Windsor Locks, Connecticut tornado of 3 October 1979: An analysis of an intermittent severe weather event. Mon. Wea. Rev., 115 , 1655–1677.
Seimon, A., 1993: Anomalous cloud-to-ground lightning in an F5-tornado-producing supercell thunderstorm on 28 August 1990. Bull. Amer. Meteor. Soc., 74 , 189–203.
Thompson, R L., Edwards R. , Hart J A. , Elmore K L. , and Markowski P M. , 2003: Close proximity soundings within supercell environments obtained from the Rapid Update Cycle. Wea. Forecasting, 18 , 1243–1261.
Wasula, A C., Bosart L F. , and LaPenta K D. , 2002: The influence of terrain on the severe weather distribution across interior eastern New York and western New England. Wea. Forecasting, 17 , 1277–1289.
Weisman, M L., 1992: The role of convectively generated rear-inflow jets in the evolution of long-lived mesoconvective systems. J. Atmos. Sci., 49 , 1826–1847.
Weisman, M L., 2001: Bow echoes: A tribute to T. T. Fujita. Bull. Amer. Meteor. Soc., 82 , 97–116.
Weisman, M L., and Davis C A. , 1998: Mechanisms for the generation of mesoscale vortices within quasi-linear convective systems. J. Atmos. Sci., 55 , 2603–2622.
Weisman, M L., and Rotunno R. , 2000: The use of vertical wind shear versus helicity in interpreting supercell dynamics. J. Atmos. Sci., 57 , 1452–1472.
Wicker, L J., 2000: The role of low-level shear, mid-level shear, and buoyancy in the intensity of modeled low-level mesocyclones. Preprints, 20th Conf. on Severe Local Storms, Orlando, FL, Amer. Meteor. Soc., 587–590.
Wilhelmson, R B., and Wicker L J. , 2001: Numerical modeling of severe local storms. Severe Convective Storms, Meteor. Monogr., No. 50, Amer. Meteor. Soc., 123–166.
Tornadoes of F3 or greater intensity since 1950 in east-central New York and adjacent western New England.
Shear values are based on output from the NWS SHARP workstation (Hart and Korotky 1991) and represent a summation of 500-m hodograph shear segments, from 0 to 6 km, showing clockwise curvature.
Analysis taken from the NWS Advanced Weather Interactive Processing System’s (AWIPS) Local Analysis and Prediction System (LAPS) analysis.