We are grateful for the support of VORTEX2 by the National Science Foundation (NSF) and National Oceanic and Atmospheric Administration. The authors were supported by awards AGS-0801035, AGS-0801041, and AGS-1036237 made to PSU, CSWR, and Rasmussen Systems, respectively. We also thank the countless number of VORTEX2 PIs, students, and other participants, without which the project would not have been possible. Radar data were edited using the SOLOII software provided by the Earth Observing Laboratory at the National Center for Atmospheric Research (NCAR). Jim LaDue provided the photograph in Fig. 9b. The video frames in Figs. 9c,d are courtesy of Nolan Atkins, Roger Wakimoto, and the Lyndon State University–NCAR photogrammetry group. The DOW radars are NSF Lower Atmospheric Observing Facilities supported by AGS-0734001. Lastly, we appreciate the selfless donation of time and constructive comments provided by the reviewers (Chris Weiss, Morris Weisman, and one anonymous reviewer) and editor (George Bryan).
Brandes, E. A., 1993: Tornadic thunderstorm characteristics determined with Doppler radar. The Tornado: Its Structure, Dynamics, Prediction, and Hazards, Geophys. Monogr., Vol. 79, Amer. Geophys. Union, 143–159.
Byko, Z., , P. Markowski, , Y. Richardson, , J. Wurman, , and E. Adlerman, 2009: Descending reflectivity cores in supercell thunderstorms observed by mobile radars and in a high-resolution numerical simulation. Wea. Forecasting, 24, 155–186.
Davies-Jones, R. P., , R. J. Trapp, , and H. B. Bluestein, 2001: Tornadoes and tornadic storms. Severe Convective Storms, Meteor. Monogr., No. 50, Amer. Meteor. Soc., 167–221.
Doviak, R. J., , P. S. Ray, , R. G. Strauch, , and L. J. Miller, 1976: Error estimation in wind fields derived from dual-Doppler radar measurement. J. Appl. Meteor., 15, 868–878.
Dowell, D. C., , and A. Shapiro, 2003: Stability of an iterative dual-Doppler wind synthesis in Cartesian coordinates. J. Atmos. Oceanic Technol., 20, 1552–1559.
Dowell, D. C., , Y. Richardson, , and J. Wurman, 2002: Observations of the formation of low-level rotation: The 5 June 2001 Sumner County, Kansas tornado. Preprints, 21st Conf. on Severe Local Storms, San Antonio, TX, Amer. Meteor. Soc., 12.3. [Available online at http://ams.confex.com/ams/pdfpapers/47335.pdf.]
Emanuel, K. A., 1994: Atmospheric Convection. Oxford University Press, 580 pp.
Gal-Chen, T., 1978: A method for the initialization of the anelastic equations: Implications for matching models with observations. Mon. Wea. Rev., 106, 587–606.
Grzych, M. L., , B. D. Lee, , and C. A. Finley, 2007: Thermodynamic analysis of supercell rear-flank downdrafts from Project ANSWERS. Mon. Wea. Rev., 135, 240–246.
Kennedy, A. D., , E. N. Rasmussen, , and J. M. Straka, 2007a: A visual observation of the 6 June 2005 descending reflectivity core. Electron. J. Severe Storms Meteor., 2 (6).
Kennedy, A. D., , J. M. Straka, , and E. N. Rasmussen, 2007b: A statistical study of the association of DRCs with supercells and tornadoes. Wea. Forecasting, 22, 1192–1199.
Koch, S. E., , M. DesJardins, , and P. J. Kocin, 1983: An interactive Barnes objective map analysis scheme for use with satellite and conventional data. J. Climate Appl. Meteor., 22, 1487–1503.
Leise, J. A., 1982: A multidimensional scale-telescoped filter and data extension package. NOAA Tech. Memo. ERL WPL-82, 19 pp. [Available from NOAA Office of Oceanic and Atmospheric Research, Silver Spring Metro Center, Bldg 3, Room 11627, Silver Spring, MD 20910.]
Lemon, L. R., , and C. A. Doswell III, 1979: Severe thunderstorm evolution and mesocyclone structure as related to tornadogenesis. Mon. Wea. Rev., 107, 1184–1197.
Majcen, M., , P. Markowski, , Y. Richardson, , D. Dowell, , and J. Wurman, 2008: Multipass objective analyses of Doppler radar data. J. Atmos. Oceanic Technol., 25, 1845–1858.
Markowski, P. M., , and Y. P. Richardson, 2008: The structure and evolution of vortex lines in supercell thunderstorms. Preprints, 24th Conf. on Severe Local Storms, Savannah, GA, Amer. Meteor. Soc., 4.1. [Available online at http://ams.confex.com/ams/24SLS/techprogram/paper_141560.htm.]
Markowski, P. M., , J. M. Straka, , and E. N. Rasmussen, 2002: Direct surface thermodynamic observations within the rear-flank downdrafts of nontornadic and tornadic supercells. Mon. Wea. Rev., 130, 1692–1721.
Markowski, P. M., , E. N. Rasmussen, , J. M. Straka, , R. P. Davies-Jones, , Y. Richardson, , and J. Trapp, 2008: Vortex lines within low-level mesocyclones obtained from pseudo-dual-Doppler radar observations. Mon. Wea. Rev., 136, 3513–3535.
Markowski, P. M., , M. Majcen, , and Y. Richardson, 2010: Near-surface vortex genesis in idealized three-dimensional numerical simulations involving a heat source and a heat sink in a vertically sheared environment. Preprints, 25th Conf. on Severe Local Storms, Denver, CO, Amer. Meteor. Soc., 15.1. [Available online at http://ams.confex.com/ams/25SLS/techprogram/paper_176112.htm.]
Markowski, P. M., , M. Majcen, , Y. Richardson, , J. Marquis, , and J. Wurman, 2011: Characteristics of the wind field in three nontornadic low-level mesocyclones observed by the Doppler on Wheels radars. Electron. J. Severe Storms Meteor., 6 (3).
Markowski, P. M., and Coauthors, 2012: The pretornadic phase of the Goshen County, Wyoming, supercell of 5 June 2009 intercepted by VORTEX2. Part II: Intensification of low-level rotation. Mon. Wea. Rev.,140, 2916–2938.
Marquis, J., , Y. Richardson, , J. Wurman, , P. Markowski, , and D. Dowell, 2008: Mobile radar observations of tornadic supercells with multiple rear-flank gust fronts. Preprints, 24th Conf. on Severe Local Storms, Savannah, GA, Amer. Meteor. Soc., 19.3. [Available online at http://ams.confex.com/ams/24SLS/techprogram/paper_142112.htm.]
Marquis, J., , Y. Richardson, , P. Markowski, , D. Dowell, , and J. Wurman, 2012: Tornado maintenance investigated with high-resolution dual-Doppler and EnKF analysis. Mon. Wea. Rev., 140, 3–27.
Okubo, A., 1970: Horizontal dispersion of floatable particles in the vicinity of velocity singularities such as convergences. Deep-Sea Res., 17, 445–454.
Pauley, P. M., , and X. Wu, 1990: The theoretical, discrete, and actual response of the Barnes objective analysis scheme for one- and two-dimensional fields. Mon. Wea. Rev., 118, 1145–1163.
Rasmussen, E. N., , J. M. Straka, , M. S. Gilmore, , and R. Davies-Jones, 2006: A preliminary survey of rear-flank descending reflectivity cores in supercell storms. Wea. Forecasting, 21, 923–938.
Ray, P. S., , C. L. Ziegler, , W. Bumgarner, , and R. J. Serafin, 1980: Single- and multiple-Doppler radar observations of tornadic storms. Mon. Wea. Rev., 108, 1607–1625.
Shabbott, C. J., , and P. M. Markowski, 2006: Surface in situ observations within the outflow of forward-flank downdrafts of supercell thunderstorms. Mon. Wea. Rev., 134, 1422–1441.
Straka, J. M., , E. N. Rasmussen, , and S. E. Fredrickson, 1996: A mobile mesonet for fine-scale meteorological observations. J. Atmos. Oceanic Technol., 13, 921–936.
Straka, J. M., , E. N. Rasmussen, , R. P. Davies-Jones, , and P. M. Markowski, 2007: An observational and idealized numerical examination of low-level counter-rotating vortices toward the rear flank of supercells. Electron. J. Severe Storms Meteor., 2 (8).
Wakimoto, R. M., , N. T. Atkins, , and J. Wurman, 2011: The LaGrange tornado during VORTEX2. Part I: Photogrammetric analysis of the tornado combined with single-Doppler radar data. Mon. Wea. Rev., 139, 2233–2258.
Waugh, S., , and S. E. Fredrickson, 2010: An improved aspirated temperature system for mobile meteorological observations, especially in severe weather. Preprints, 25th Conf. on Severe Local Storms, Denver, CO, Amer. Meteor. Soc., P5.2. [Available online at http://ams.confex.com/ams/25SLS/techprogram/paper_176205.htm.]
Weiss, C. C., , and J. L. Schroeder, 2008: StickNet—A new portable, rapidly deployable, surface observing system. Preprints, 24th Conf. on Interactive Information and Processing Systems, New Orleans, LA, Amer. Meteor. Soc., 4A.1. [Available online at http://ams.confex.com/ams/88Annual/techprogram/paper_134047.htm.]
Wurman, J., , J. Straka, , E. Rasmussen, , M. Randall, , and A. Zahrai, 1997: Design and deployment of a portable, pencil-beam, pulsed, 3-cm Doppler radar. J. Atmos. Oceanic Technol., 14, 1502–1512.
Wurman, J., , Y. Richardson, , C. Alexander, , S. Weygandt, , and P. F. Zhang, 2007: Dual-Doppler and single-Doppler analysis of a tornadic storm undergoing mergers and repeated tornadogenesis. Mon. Wea. Rev., 135, 736–758.
Wurman, J., , K. Kosiba, , P. Markowski, , Y. Richardson, , and D. Dowell, 2010: Finescale single- and dual-Doppler analysis of tornado intensification, maintenance, and dissipation in the Orleans, Nebraska, tornadic supercell. Mon. Wea. Rev., 138, 4439–4455.
Wurman, J., , D. Dowell, , Y. Richardson, , P. Markowski, , E. Rasmussen, , D. Burgess, , L. Wicker, , and H. Bluestein, 2012: The Second Verification of the Origins of Rotation in Tornadoes Experiment: VORTEX2. Bull. Amer. Meteor. Soc., 93, 1147–1170.
The Goshen County storm also has been referred to as the “LaGrange storm” (e.g., Wakimoto et al. 2011).
Circulation centers were identified at the location of the minimum in the field of Okubo–Weiss number (W; Okubo 1970; Weiss 1991) following Markowski et al. (2011). The W field was further smoothed in order to filter submesocyclone-scale features. Scales smaller than 2 km were filtered using a three-step Leise (1982) filter, which suppresses wavelengths ≤8Δ, where Δ is the grid spacing.
Although the mean motion was faster and less rightward in the 2100–2130 UTC period, the grid to which raw radar observations from this period were interpolated was translated at the same constant velocity as in the 2130–2148 UTC period. Given that the main reason for using a translating grid was to reduce time interpolation errors in the trajectory calculations, we were not concerned with a grid-relative drift of the storm during the single-Doppler observation period.
The azimuthal wind shear is defined as r−1∂υR/∂ϕ, where r is the range from the radar, υR is the radial velocity, and ϕ is the radar azimuth (not to be confused with the mathematical coordinate).
The time series of maximum azimuthal shear at low levels also is much noisier than at midlevels. One possibility is that the smaller scale of the low-level mesocyclone, relative to the midlevel mesocyclone, makes its intensity more dependent on the location of the radar beams relative to the axis of rotation. The relatively noisy time series of maximum low-level azimuthal shear also could be attributable to the influence of the gust front. Azimuthal shear can be associated with deformation or rotation. Given that gust fronts can be associated with significant deformation and rotation, whereas mesocyclones are dominated by rotation, one might expect low-level azimuthal shear to be intrinsically more variable in time than midlevel shear.
Vortex lines were computed using a fourth-order Runge–Kutta algorithm. Velocity derivatives were computed using fourth-order, centered differences, except near data boundaries where second-order, uncentered differences were used. The qualitative characteristics of the vortex lines are robust for reasonable ranges of smoothing [e.g., see the appendix in Markowski et al. (2008)]. Vortex-line calculations originating within the midlevel mesocyclone (e.g., the blue lines in Fig. 10) are qualitatively insensitive to the observed variations in the magnitude of the horizontal gradient of vertical velocity, |∇hw|, within the midlevel mesocyclone [vortex lines were computed from additional locations within the midlevel mesocyclone (not shown) spanning a wide range of |∇hw|], which suggests that the vortex lines computed from the fixed lattices of points in the midlevel mesocyclone also are qualitatively insensitive to the |∇hw| errors (errors in the retrieved w, which can be large at midlevels, potentially could be accompanied by significant errors in |∇hw|).
Even at times when the data horizon is low, the errors in the vertical velocity field accumulate with altitude (e.g., Doviak et al. 1976); thus, all retrieved midlevel vertical velocities in this study should be interpreted only qualitatively (e.g., the maximum vertical velocities evident in Fig. 13 are not credible).
Even if the mobile mesonet could have reached the storm earlier, it is unlikely that they would have been able to sample this region given the paucity of roads in the area; what few roads existed were generally unpaved and unfriendly to data collection in heavy precipitation.
The trajectories were integrated using a fourth-order Runge–Kutta algorithm (as was done for the vortex lines). A time step of 15 s was used, and the 3D wind fields were assumed to vary linearly in time between dual-Doppler analyses. The spatial interpolation was trilinear.
Bolton’s (1980) formula for pseudoequivalent potential temperature was used. When referring to θe throughout the article, we technically are referring to pseudoequivalent potential temperature, which is conserved for dry adiabatic and pseudoadiabatic processes.
This would require a hydrometeor mixing ratio of roughly 7 g kg−1.
To put this gradient into perspective, a parcel of air residing in a 1 K km−1 horizontal θρ gradient would acquire a horizontal vorticity of 0.01 s−1 in 5 min.