Donald J. Emanuel - NWSFO Great Falls, MT
The 500 mb chart from 1200 UTC 11 June showed an upper-level ridge axis over eastern Montana. By 0000 UTC 12 June, the axis had shifted to the Montana-Dakotas border (Figs. 2-3). This placed north-central Montana under southwest flow aloft, a favorable regime for summer convective activity. The convective indices from the 1200 UTC Great Falls sounding on the morning of 11 June indicated marginal stability over the area. Personal Computer Gridded Interactive Display and Diagnostic System (PCGRIDDS) analysis of the Eta model, as well as Skew-T/Hodograph Analysis and Research Program (SHARP) output, showed a Lifted Index (LI) value of zeroøC and a Convective Available Potential Energy (CAPE) value of zero J kg-1. Storm relative helicity in the lowest three kilometers (km) was a low 15 m2s-2. Dew point temperatures were in the 40 to 45 degree Fahrenheit range. The airmass west of Great Falls was destabilizing throughout the afternoon hours as indicated by the Great Falls sounding of 0000 UTC 12 June. PCGRIDDS and SHARP output indicated LI values of -2 to -3øC as well as a CAPE value of 754 J kg-1 and a helicity value of 105 m2s-2 . Dew point temperatures were now in the lower 50s (0F).
During the mid-afternoon of 11 June, WSR-88D radar at Great Falls showed a number of thunderstorms north and east of the city. At 2214 UTC, an outflow boundary became detectable on both base reflectivity and base velocity images. This outflow boundary originated from thunderstorms in eastern Choteau County and was not visible on radar until it reached Teton County (Fig. 4), northwest of Great Falls.
McCollum et al., has found that drainage winds can help initiate convection by providing a mechanism for lift. In this case, a convergence zone intensified as southwest drainage winds collided with northeast surface winds, originating from the outflow boundary (Fig. 4). Furthermore, research by Maddox has found that the presence of a dry adiabatic lapse rate within the lower atmosphere can help generate a greater response in the vertical motion field for parcels that are being acted upon and lifted by drainage winds (personal communication, Maddox). The 0000 UTC 12 June Great Falls sounding did show a fairly dry adiabatic lapse rate throughout the lower atmosphere (Fig. 5). As stated previously, the airmass through which the outflow boundary propagated continued to destabilize. PCGRIDDS analysis of the 0000 UTC 12 June Eta model output indicated an axis of weak convergence through the surface to 700 mb layer. Additionally, weak divergence was occurring through the 400-300 mb layer (not shown).
Also indicated by the Eta model was a large area of isentropic potential vorticity. The concept of isentropic potential vorticity anomalies has been widely researched and documented (Bluestein, 1993). The mathematical foundation for their existence is beyond the scope of this TA. In general terms, the theory behind potential vorticity anomalies holds that they exhibit a specific circulation (much like an electric charge). Positive anomalies are associated with cyclones and negative anomalies with anticyclones. The values are expressed in terms of Potential Vorticity Units (1 PVU=10-6m-2s-1K kg-1).
Analysis of the 300-200 mb layer from 1200 UTC 11 June showed values of +3 PVUs across northern Montana with values of +2 PVU stretching from west-central Montana to the southwest corner of the state (Figs. 6-7). Researchers agree that potential vorticity values greater than 1.5 PVU are usually associated with air of stratospheric origin (Bluestein, 1993). By 0000 UTC 12 June, the area of +3 PVU had moved further south and east (Fig. 8). This would indicate that stratospheric intrusion of potential vorticity had occurred in this area. Analysis of the 1200 UTC Great Falls sounding from 11 June showed the tropopause at a height of 10.4 km MSL. By 0000 UTC 12 June, the tropopause was at 10.9 km MSL, indicating the intrusion had lessened somewhat, as corroborated in (Fig. 9).
The depth of these stratospheric intrusions must also be considered. The depth of penetration is somewhat dependent upon the static stability of the layer. Low static stability , i.e., an unstable lapse rate, would allow a larger penetration than if the atmosphere were more stable. PCGRIDDS analysis of the 850-400 mb layer from 0000 UTC 12 June showed a static stability value of 7.00 C km-1 over the area west of Great Falls (Fig. 10). By 0600 UTC, this rather unstable lapse rate remained in this same area.
Additional analysis of 0000 UTC 12 June data offered a few more clues. Wind data from 250 mb showed the core of the jet streak/stream to be moving across southern Montana (Fig. 11). Isohypsic data from the same level indicated divergence was occurring, as quantified by PC GRIDDS data. At 200 mb, the data suggested positive vorticity advection (PVA) was occurring in the presence of a weak short wave (not shown).
The southern end of the outflow boundary was less conspicuous on radar imagery and difficult to discern. It seemed likely that convection was initiated along the southern edge of the boundary due to combined effects of evening drainage flows and the instability of the airmass due, in part, to the stratospheric intrusion of potential vorticity. Once convection initiated at the boundary's southern end, it developed rapidly. New thunderstorms developed west of Great Falls around 0500 UTC with a maximum reflectivity of 47 dBz and vertically integrated liquid water content (VIL) values of 27 k gm-2. Echo tops were fairly constant at 38000 feet above ground level (AGL). By 0607 UTC, a second thunderstorm developed five miles west of Great Falls. This storm closely followed the track of the first while rapidly intensifying. Echo tops for this second storm reached 40000 feet AGL, a rather impressive height given the time of day. Hail measuring 3/4" was reported at the NWS office at 0622 UTC (Fig. 13). Reflectivities remained in the 60 to 64 dBz range as the storm moved rapidly northeastward, dissipating several hours later over eastern Montana. No additional severe weather was reported after these storms moved beyond the Great Falls area.
Barnes, S. L. and B. R. Colman, 1994: Diagnosing an operational numerical model using
Bluestein, H. B., 1993: Synoptic-Dynamic Meteorology in Midlatitudes, Vol. II., Oxford
McCollum, D. M., R. A. Maddox, and K. W. Howard, 1995: Case study of a severe
Thaler, E., 1995: Isentropic analysis and IPV, COMAP Notes class 95-1.