1. Introduction
The role of vertical wind shear (hereafter “shear”) in squall-line evolution has been extensively addressed over the last half century. For instance, the strength and direction of low-level (LL) shear that occurs over the depth of the cold pool regulates updraft tilt and speed (e.g., Hane 1973; Thorpe et al. 1982; Rotunno et al. 1988; Fovell and Dailey 1995; Grady and Verlinde 1997; Davis et al. 2004; Coniglio et al. 2006). This connection exists because the ratio of the LL shear magnitude to the integrated cold pool buoyancy (hereafter simply “cold pool intensity”) strongly influences the character of the lifting along the outflow boundary (OFB) that underlies the buoyant updrafts within the squall line (e.g., Rotunno et al. 1988; Weisman and Rotunno 2000; Bryan and Rotunno 2014). When the cold pool intensity overwhelms the LL shear or the LL shear overwhelms the cold pool intensity, ascent occurs in a slantwise manner along the OFB within the squall line’s buoyant updrafts. Alternatively, when LL shear and cold pool intensity approximately balance each other, OFB ascent and buoyant updrafts are comparatively upright. These ideas are encapsulated by “RKW” theory (Rotunno et al. 1988). Updrafts that are more slanted have larger downward-oriented pressure gradient accelerations (PGA) that weaken their ascent rates relative to less slanted updrafts, all else being equal (e.g., Parker 2010; Peters 2016; Peters et al. 2019a). Furthermore, the magnitude of dynamically driven lifting along an OFB positively correlates with the LL shear magnitude, so long as the LL shear magnitude does not “overpower” the cold pool. These studies suggest that squall lines in environments with stronger LL shear should have stronger updrafts than those in environments with weaker LL shear, all else being equal.
Oftentimes, substantial upper-level (UL) shear is also found within squall-line environments (e.g., see Fig. 6 in Davis et al. 2004; see Fig. 21 in Coniglio et al. 2006). Fovell and Dailey (1995) and Coniglio et al. (2006) showed that both LL and UL shear work together to result in a deeper overturning layer, and stronger maximum vertical velocities (wMAX), than would occur if UL shear were absent and only LL shear were present. In addition, Moncrieff and Liu (1999) noted that certain orientations of UL shear result in deep convective updrafts that “anchor” vertically above an OFB, which is an optimal situation for squall-line maintenance. Other UL shear orientations prevent the anchoring of storms along an OFB, which may result in less efficient triggering and maintenance of deep convection. Strong mid- to upper-level winds also influence the cold pool speed via downward momentum transport (Mahoney et al. 2009), which in turn alters the balance between the cold pool and the LL shear. Furthermore, strong UL shear may also influence squall-line motion via changes in steering flows aloft, and thus, LL OFB-relative flow and LL horizontal mass flux. Finally, studies of model-derived soundings taken in proximity to squall lines suggest that their probability of maintenance increases as UL shear increases (e.g., Coniglio et al. 2007). Most of the aforementioned numerical modeling studies varied LL (UL) shear while holding UL (LL) shear fixed, so the importance of variations in both LL and UL shear collectively on squall-line updraft properties is largely unknown, as stated in Coniglio et al. (2006).
One updraft process that has seldom been discussed in the context of squall lines is entrainment. Entrainment involves the mixing of updraft air with surrounding environmental air. Entrainment typically results in dilution of updraft buoyancy (B), leading to a reduction in vertical accelerations and vertical velocities (w). Understanding entrainment within deep convective updrafts is vital, as this process modulates updraft w, in turn influencing how updrafts vertically redistribute heat, moisture, momentum, and chemical constituents. Updraft w also plays a key role in the generation of large hail (e.g., Browning 1963; Wakimoto et al. 2004), damaging winds (e.g., Marion and Trapp 2019), tornadoes (e.g., Markowski and Richardson 2014; Coffer and Parker 2015), lightning (e.g., Lang and Rutledge 2002; Tessendorf et al. 2005; Wiens et al. 2005; Deierling and Petersen 2008), and heavy rainfall (e.g., Doswell et al. 1996; Nielsen and Schumacher 2018).
Canonical studies of squall lines have typically depicted a structure with deep convective updrafts along an OFB that are driven by a “tongue” of high equivalent potential temperature and conditionally unstable air that slopes up and over the cold pool—colloquially known as front-to-rear flow (e.g., Zipser 1977; Houze 1989; Kingsmill and Houze 1999; Bryan and Fritsch 2000). This type of layer lifting is commonly referred to as “slab-like” ascent (e.g., Zipser 1977; Chong et al. 1987; Redelsperger and Lafore 1988; Houze 1989; Jorgensen et al. 1997; Trier et al. 1997; Kingsmill and Houze 1999; Bryan and Fritsch 2000; Mechem et al. 2002; James et al. 2005; Alfaro and Khairoutdinov 2015; Alfaro 2017). The role of entrainment in these regions of slab-like ascent is presently unclear, and it is possible that entrainment minimally affects updraft air because of the laminar flow within these layer-lifting regions. However, this slab-like airstream often breaks down into discrete, quasi-spherical updraft pulses—commonly known as thermals—in the middle to upper troposphere (e.g., Redelsperger and Lafore 1988; Weisman et al. 1988; Fovell and Tan 1998; Bryan and Fritsch 2000; Bryan et al. 2003, 2007; Bryan and Morrison 2012; Lebo and Morrison 2015). The diameters of these thermals are typically on the order of 1–2 km, which is more than an order of magnitude smaller than the horizontal scale of the slab-like airstream; however, it is conceivable that the horizontal scale of the ascending slab-like airstream “sets the stage” for the width of thermal-like updrafts aloft in squall lines, as hypothesized by James et al. (2005). The structure of thermals generally plays an essential role in the entrainment behavior of deep convection, and there is evidence that thermal-like updrafts are more susceptible to the deleterious effects of entrainment-driven dilution of B than other updraft organizational structures, such as plumes (e.g., Morton et al. 1956; Scorer 1957; Morrison 2017; Morrison et al. 2020; Peters et al. 2020a,c). Consequently, the thermal-like structure of middle- to upper-tropospheric squall-line updrafts suggest that they are susceptible to the deleterious effects of entrainment.
It is possible that shear influences squall-line entrainment characteristics, in addition to the dynamical connections that have been described by previous authors and are outlined above, though this potential connection has received little attention in previous literature. For instance, the influence of shear on supercellular updraft width and intensity has received growing attention in recent literature. Studies such as Dennis and Kumjian (2017), Trapp et al. (2017), Warren et al. (2017), Marion and Trapp (2019), Peters et al. (2019b, 2020a,b), and Mulholland et al. (2020) have all shown supercell updraft width to positively scale with shear magnitude (i.e., stronger shear is correlated with wider supercell updrafts). Specifically, Peters et al. (2019b) provided theoretical and idealized numerical model simulation results that showed increasing shear led to increasing LL storm-relative inflow (i.e., greater LL horizontal mass flux) and, through a mass continuity argument (see section 4 therein), resulted in wider supercell updrafts. Wider supercell updrafts were shown to be more resistant to entrainment-driven dilution of B, and therefore had larger vertical accelerations and w (e.g., Peters et al. 2019b, 2020b). It is natural to wonder if this concept also applies to squall lines. An idealized numerical modeling study of squall lines by French and Parker (2010) showed that increasing LL shear in squall line environments equated to larger vertical mass flux in the deep convective line. An idealized numerical modeling study of squall lines by Alfaro (2017) found that increasing LL shear resulted in a greater fraction of inflowing high-CAPE air to the total inflowing air, resulting in greater B aloft in simulated squall lines. Alfaro (2017) argued that the greater fraction of inflowing high-CAPE air to the total inflowing air resulted in less dilution of updrafts aloft in the simulated squall lines, although that study did not explicitly show variations in entrainment-driven dilution and/or changes in updraft width with changes in LL shear. We thus build upon the Alfaro (2017) results and pose the following central research questions: 1) Do environments with stronger shear also foster wider squall-line updrafts than environments with weaker shear, and are squall-line updrafts in stronger-shear environments therefore less susceptible to entrainment-driven dilution? 2) If such a relationship exists, what are the relative influences of LL and UL shear on this relationship? Some of these research questions were further motivated by the pathways for future work outlined in Marinescu et al. (2017) and results shown in Phinney et al. (2017). Specifically, Marinescu et al. (2017) suggested that future squall-line simulation studies use large-eddy-resolving scales to further evaluate how environmental factors, such as wind shear, modulate entrainment within squall lines.
In response to these research questions, we evaluate the following hypothesis: stronger LL shear in squall-line environments equate to stronger LL OFB-relative flow, greater LL horizontal mass flux into the edge of the cold pool, and thus wider squall-line updrafts with less entrainment-driven dilution than in environments where LL shear is comparatively weak. The wider squall-line updrafts in environments with stronger LL shear equate to larger updraft B, and consequently w, when compared to squall-line updrafts amid environments with weaker LL shear. Furthermore, we expect a similar relationship between UL shear, which may alter squall-line steering currents and OFB-relative horizontal mass flux, updraft width, entrainment, and intensity. Here we show that, in addition to the dynamics discussed by previous authors, the specific connection between LL shear, updraft width, and entrainment is an essential missing piece explaining why LL shear tends to modulate the intensity (i.e., wMAX) of deep convective updrafts in squall lines. Finally, although we expect LL and UL shear to have similar impacts on squall-line updraft properties, we hypothesize that LL shear should have a greater influence on these aforementioned processes than UL shear owing to the more obvious connection between LL shear and both LL OFB-relative flow and LL horizontal mass flux (e.g., Parker and Johnson 2004b; French and Parker 2010; Alfaro 2017). To address this hypothesis, we analyzed a suite of idealized numerical model simulations of squall lines with variations solely in the background LL and UL shear profiles for a singular thermodynamic environment (outlined in section 2). Results from the suite of idealized numerical model simulations are detailed in section 3 and a summary, discussion, and conclusions are given in section 4.
2. Experimental design
a. Numerical modeling setup
We conducted idealized numerical model simulations using Cloud Model 1 (CM1; Bryan and Fritsch 2002) release 19.7. All simulations used isotropic grid spacings of 250 m in the horizontal and vertical directions. The domain top was set to 25 km above ground level (AGL) with a Rayleigh dampening layer (coefficient = 3.33 × 10−3 s−1) applied above 20 km AGL. The north–south domain length was 99 km with periodic lateral boundary conditions while the east–west domain length was 420 km with open radiative lateral boundary conditions. Coarser horizontal grid spacing (Δx, Δy = 1 km) simulations (with and without east–west horizontal grid stretching) were conducted with various domain sizes; the results of those simulations were qualitatively similar to those presented in section 3 for the higher-resolution simulations. The lower and upper boundaries were free slip with a rigid domain top. In all simulations, we neglected radiation, surface fluxes, friction, terrain, and Coriolis acceleration. Microphysical processes were parameterized with the Morrison two-moment scheme (Morrison et al. 2009) with hail set as the prognostic rimed ice hydrometeor species. The simulations used a time step of 3.5 s with integration lengths of 8 h and output written every 5 min. A summary of all CM1 settings is given in Table 1.
Summary of the CM1 configuration. LBC—Lateral boundary condition.
b. Description of base-state environments
The thermodynamic base-state environments for all simulations were identical and used the Weisman and Klemp (1982, hereafter WK82) analytic equations (Fig. 1). Modifications to the original WK82 sounding included increasing the surface potential temperature θ from 300 to 302 K and reducing the relative humidity (RH) above 2.5 km AGL to a constant value of 45%. The slight increase in surface θ (from 300 to 302 K) was made to reduce the amount of errant convection that occurred ahead of the squall-line OFB in the simulations, while the reduction in RH above 2.5 km AGL to 45% was to create an environment more in-line with observed squall-line environments (e.g., Fig. 7 in Coniglio et al. 2012). The modified planetary boundary layer (PBL) mixing ratio was set to 14 g kg−1, which resulted in surface-based, 0–1 km mean layer, and most-unstable convective available potential energy (CAPE) [convective inhibition (CIN)] of 2043 (−69), 2191 (−44), and 2615 (−3) J kg−1, respectively—all of which fall within the range of observed values for squall-line environments (e.g., see Fig. 9 in Evans and Doswell 2001; see Fig. 8 in Cohen et al. 2007). The approximate effective inflow layer (EIL; Thompson et al. 2007) depth following the CAPE and CIN criteria of Thompson et al. (2007) (i.e., largest contiguous vertical layer with CAPE > 100 J kg−1 and CIN < −250 J kg−1) was ~2.5 km (cf. Figs. 2c,d).
Deep convection was initiated in the simulations with a cold pool covering the domain from south to north and extending from the center of the domain westward. The minimum surface θ perturbation of the cold pool was set to −5 K and decreased in magnitude toward 0 K with increasing height through the 2.5 km depth of the cold pool (Fig. 2a). To facilitate the development of turbulence, ±0.25 K random θ perturbations were seeded throughout the depth and horizontal extent of the cold pool in the initial conditions (e.g., as in Parker and Johnson 2004b). Two sensitivity simulations (not shown) with initial ±0.25 K random θ perturbations placed throughout the entire domain below 2.5 km AGL, instead of just confined within the cold pool, yielded nearly identical results as the configuration we used here.
A set of various initial vertical wind profiles were used in the simulations to address our hypothesis. All wind profiles featured unidirectional west-to-east shear. This unidirectional “shear” profile was chosen for its simplicity and extensive use in previous similar squall-line simulations (e.g., Rotunno et al. 1988; Weisman et al. 1988; Fovell and Ogura 1989; Skamarock et al. 1994; Fovell and Dailey 1995; Weisman et al. 1997; Fovell and Tan 1998; Parker and Johnson 2004a,b,c; Weisman and Rotunno 2004; James et al. 2005; Bryan et al. 2006; Coniglio et al. 2006; James et al. 2006; Lane and Moncrieff 2010; Alfaro and Khairoutdinov 2015; Lebo and Morrison 2015; Peters 2016; Alfaro 2017). Three different LL (0–2.5 km AGL) and UL (2.5–10 km AGL) vertical wind shear1 magnitudes constituted our simulation matrix, all of which generally fall within the range of shear magnitudes used in prior squall-line simulation studies (e.g., Parker and Johnson 2004a,b,c; Coniglio et al. 2006; Alfaro and Khairoutdinov 2015; Alfaro 2017) and also in observations (e.g., Bluestein and Jain 1985; Parker and Johnson 2000; Evans and Doswell 2001; Coniglio et al. 2004; Parker and Johnson 2004a; Cohen et al. 2007; Coniglio et al. 2007; Bryan and Parker 2010; Lombardo and Colle 2012; Hitchcock et al. 2019) (Fig. 2b; Table 2). The naming convention for the various simulations uses numbers from 1 to 3 to indicate the relative LL and UL shear magnitudes. For example, the LL1_UL3 simulation is characterized by the weakest LL shear (7 m s−1) and the strongest UL shear (10 m s−1) shown in Table 2. It should be noted that our LL3 shear profile likely represents the far high end of the LL shear spectrum for observed squall lines in environments with substantial CAPE (e.g., see Fig. 13 in Evans and Doswell 2001; see Figs. 1 and 2 in Cohen et al. 2007; see Fig. 2 in Coniglio et al. 2007). Furthermore, our strongest UL shear profile (UL3) is likely at the low to moderate end of the UL shear spectrum for observed squall lines in environments with substantial CAPE (e.g., see Figs. 1 and 2 in Cohen et al. 2007; see Fig. 2 in Coniglio et al. 2007). Finally, to approximately center the squall-line OFBs within the domains throughout the length of the simulations, a unique uniform zonal wind speed was subtracted from the entire domain for each simulation, which included regions inside and outside of the cold pool.
Magnitude of LL (0–2.5 km AGL) and UL (2.5–10 km AGL) shear (calculated as a simple bulk wind difference and expressed in units of m s−1) for all nine squall-line simulations. The naming conventions follow that outlined in section 2.
c. Quantifying entrainment-driven dilution
To quantify entrainment-driven dilution of updraft cores in the squall lines, a passive tracer (PT) was initialized in the base state to the east of the cold pool between 0 and 1.5 km AGL (light red shading in Figs. 2c and 2d), or approximately throughout the lower half of the EIL in which the largest magnitudes of CAPE were located (e.g., Fig. 2c; similar PT setup as in Phinney et al. 2017). In subsequent analyses, the degree to which the PT concentration in updrafts falls below 100% is used to quantify the degree of entrainment-driven updraft dilution.
d. In-line forward parcel trajectories
As a final test of our hypothesis, we computed in-line forward trajectories2 in the simulations with the strongest and weakest LL and UL shear combinations (i.e., LL1_UL1, LL1_UL3, LL3_UL1, and LL3_UL3). At four hours into the simulations, 1000 trajectories were initialized ahead of the squall-line OFBs. Trajectories were placed in a 10 km south-to-north line at 100 km east of the domain center and spanned between 25 and 835 m AGL, covering the lower EIL. Trajectories were calculated with a time step consistent with the model simulations (3.5 s) with output written every time step. To test if alterations in shear affected squall-line updraft B, and ultimately wMAX, we compared the contribution of B versus the contribution of vertical perturbation pressure gradient accelerations (VPG) to wMAX along the trajectories among the four different simulations.
3. Results
a. General attributes of simulations
The 4–7-h3 averaged lowest model level (125 m AGL) reflectivity (computed relative to the OFB4 each 5-min output time) revealed that all nine simulations produced squall lines with many characteristic features of observed squall lines, such as a stratiform region that trailed the leading deep convective line (Fig. 3; also see radar animations in Fig. S1 in the online supplemental material). Simulations with the weakest LL shear displayed notably broken and cellular reflectivity features displaced rearward of their OFBs (e.g., far left column of Fig. 3 and Fig. S1). Increasing LL shear led to more widespread and coherent reflectivity structures with a nearly continuous south-to-north deep convective region that was located closer to the OFB (black dashed lines in Fig. 3 and Fig. S1), in agreement with Hane (1973) and numerous other studies. Conversely, alterations in UL shear had relatively muted impacts on the squall lines’ organizational structures, especially in LL2 simulations. Interestingly, increasing UL shear in the LL1 simulations actually resulted in slightly weaker reflectivity magnitudes closer in proximity to the OFBs (e.g., compare LL1_UL2 and LL1_UL3 in Fig. 3 and Fig. S1). Alternatively, increasing UL shear in the LL3 simulations resulted in slightly stronger reflectivity magnitudes and more three-dimensional structures (e.g., bowing segments in LL3_UL3)—in agreement with the squall-line simulations of Coniglio et al. (2006).
West-to-east oriented, OFB-relative vertical cross sections of PT concentration, w, and reflectivity averaged −5 to +5 km around the center of the domain (y = 49.5 km) and over the 4–7-h period revealed stark differences in the convective features of the squall lines (Fig. 4; also see animations in Fig. S2). Simulations with weaker LL shear had correspondingly weaker and shallower updrafts whereas simulations with stronger LL shear had stronger and deeper updrafts. Furthermore, as noted in the plan view reflectivity composites (Fig. 3), the low-level OFBs were closer in proximity to the leading deep convective line in the stronger LL shear simulations, whereas the OFBs were effectively detached from the middle- to upper-tropospheric updrafts in LL1 simulations. Last, weaker LL shear simulations displayed greater dilution of the EIL PT aloft compared to the stronger LL shear simulations.
To provide context for comparisons with previous studies of squall lines, we assessed the relative “balance” between the cold pool-induced circulation and LL shear-induced circulation for the simulated squall lines. This was accomplished by calculating the C/Δu ratio (Rotunno et al. 1988), where C is the integrated cold pool strength equal to
b. Influence of shear on updraft width and entrainment-driven dilution
A first step toward evaluating our hypothesis was to determine whether the simulations with stronger shear displayed wider updrafts than the simulations with weaker shear. Updraft area was quantified by first calculating contiguous horizontal regions that met a certain w threshold at each model vertical level. We then took an average of all such updraft area “objects” at each model vertical level and at each time, excluding areas of small-scale vertical motion (i.e., contiguous areas ≤ 1 km2). Finally, we time averaged updraft area “objects” at each model vertical level over the 4–7-h period. Updraft area “objects” were characterized using two different w thresholds: 1) “broad” (w = 5 m s−1) and 2) “core” (w = 20 m s−1) updrafts. Our results were relatively insensitive to small variations (~few meters per second) in these w thresholds. Broad updraft area increased with increasing LL shear, especially below 5 km AGL (Fig. 6a). Differences in broad updraft area due to variations in UL shear were less pronounced than for variations in LL shear, with a slight increase in broad updraft area with stronger UL shear in LL2 and LL3 simulations and a slight decrease in broad updraft area with stronger UL shear in LL1 simulations. Similar to broad updraft area, core updraft area also generally increased with increasing LL shear, with lesser differences due to variations in UL shear (Fig. 6b). The large “spikes” in broad updraft area below ~2.5 km AGL were associated with the front-to-rear slab-like ascending airstreams that characterized these squall lines. The larger regions of slab-like ascent in the stronger LL shear simulations—which likely owe to stronger LL shear and stronger, deeper cold pools (see Fig. 5; e.g., James et al. 2005)—appear to “set the stage” for wider thermals aloft. Statistically significant correlations (linear correlation coefficient CC = 0.81; p value < 0.05) between average 0–5 km AGL broad updraft area and average 5–10 km AGL core updraft area across all simulations over the 4–7-h period show strong support for this idea (Fig. 7). A heat map5 of the median 4–7-h average 0–5 km AGL broad updraft area shows that as LL shear increased, broad updraft area increased (Fig. 8a). A similar trend was noted for median 4–7-h average 5–10 km AGL core updraft area (Fig. 8b). Interestingly, stronger UL shear with LL1 simulations resulted in slightly smaller broad and core updraft area than weaker UL shear with LL1 simulations (far left columns in Figs. 8a and 8b), whereas stronger UL shear with LL2 and LL3 simulations resulted in comparable broad and core updraft area relative to weaker UL shear with LL2 and LL3 simulations. Thus, updraft width appears to be more strongly influenced by LL shear than by UL shear.
To quantify the potential influence of updraft width on entrainment-driven dilution, we compared PT concentrations aloft across the simulations. A heat map constructed using the median 4–7-h 5–10 km AGL layer maximum PT concentrations show that stronger LL shear equated to larger PT concentrations aloft, relative to weaker LL shear (Fig. 9). In contrast, but in a consistent manner with the trends in updraft width, stronger UL shear in the LL1 simulations actually resulted in smaller PT concentrations aloft when compared to weaker UL shear in the LL1 simulations. However, there was no statistically significant connection between UL shear and PT concentration in the LL2 and LL3 simulations, which is again consistent with the trends in updraft width. Correlations between 4–7-h 5–10 km AGL layer maximum PT concentrations and average 5–10 km AGL core updraft area were statistically significant (i.e., p value <0.05), with a log-scaled CC = 0.83 (Fig. 10). These results reveal that squall-line updrafts in environments with stronger LL shear displayed less entrainment-driven dilution compared to squall-line updrafts in environments with weaker LL shear. The connection between LL shear and updraft width is evident in the statistically significant correlations (CC = 0.79; p value <0.05) between 4–7-h 0–1 km AGL horizontal mass flux and average 5–10 km AGL core updraft area across all simulations (Fig. 11). The LL shear is directly connected to the LL horizontal mass flux (e.g., French and Parker 2010) and, in turn, connected to the vertical mass flux via mass continuity. In an analogous manner to what occurs in supercells (e.g., Peters et al. 2019b), the increase in vertical mass flux equates to wider updrafts.
c. Theoretical support for the relationship between updraft width and entrainment-driven dilution
Axisymmetric updrafts were assumed in deriving the theoretical expression for PTTH concentration. Because our simulations are fully 3D, we assume that R is proportional to the effective radius Reff of a circular region with area equal to the simulated 5–10 km AGL core updraft area A. In other words, R = cReff and Reff = (A/π)1/2. The constant of proportionality c is found by approximately matching the solution for PTTH concentration from Eq. (1) with the PT concentration from the simulation with the largest core A and Reff, and then applying this value of c across the full range of A and Reff from the simulations. Note that modifying c simply shifts the PTTH concentration curve up or down, but does not modify the shape of the relationship between PTTH concentration with A. For simplicity, we set height z in Eq. (1) to 5 km AGL; while there is variability in height of maximum PT concentrations at different times, it generally occurs at ~5 km AGL (e.g., Fig. S3).
The closeness of the PTTH concentration curve (green line in Fig. 10) and the 4–7-h 5–10 km AGL layer maximum PT concentrations from the simulations supports a 1/R2 scaling for fractional entrainment rate and dilution of the simulated updrafts. The theoretical 1/R2 scaling comes directly from the assumption of constant L in Eq. (2). If instead L is assumed to be proportional to R, the classical 1/R scaling of fractional entrainment rate from dimensional analysis of dry plumes and thermals (e.g., Morton et al. 1956; Scorer 1957) is obtained from Eq. (2). However, this scaling gives a much poorer comparison of PTTH concentration to the simulated PT concentration (gold line in Fig. 10). Specific reasons for this behavior are unclear and are beyond the scope of the study. Regardless, this analysis strongly supports an inverse relationship between updraft area and entrainment-driven dilution of the simulated squall-line updrafts, which is in line with our primary hypothesis.
d. Impact of differences in updraft width and entrainment-driven dilution on maximum vertical velocities
Next, we establish a relationship between updraft dilution and wMAX among the simulations, which is a key to our hypothesis. A heat map displaying a matrix of median 4–7-h wMAX magnitudes for the various LL and UL shear combinations was analyzed first (Fig. 12). Visual inspection of the heat map in Fig. 12 reveals that as LL shear increased, wMAX increased. Consistent with earlier analysis, stronger UL with LL1 simulations actually resulted in slightly weaker wMAX relative to weaker UL shear with LL1 simulations (far left column of Fig. 12), whereas the trend between UL shear and wMAX in the LL2 and LL3 simulations was ill defined. To provide a baseline for comparison with parcel theory, the “thermodynamic speed limit” (i.e.,
Vertical profiles of wMAX averaged over the 4–7-h period depicted an increase in the magnitude and height of wMAX with stronger LL shear (Fig. 13). Stronger LL shear equates to stronger wMAX along the LL OFBs (along and below ~2.5 km AGL) across all simulations when compared to weaker LL shear, as previously shown in numerous studies. Note, however, that the height of wMAX, and the largest variations in wMAX between the various LL and UL shear simulations, corresponded with thermals aloft (>5 km AGL) and not with the dynamically driven low-level ascent along the OFB. Variations in UL shear had the largest impact in LL3 simulations, with a monotonic increase in both the height and magnitude of wMAX with increasing UL shear (red shaded lines in Fig. 13). Furthermore, we examined the maximum forward parcel trajectory displacements, as it is a combination of both w and the time it takes parcels to ascend that determines parcel displacements. Figure 14 displays these results and shows that, in line with the aforementioned wMAX analyses, maximum parcel displacements were substantially larger in the stronger LL shear simulations. Thus, updrafts experiencing stronger LL shear were faster and deeper than those experiencing weaker LL shear. Variations in UL shear had lesser impact on maximum parcel displacements, although there still was a notable shift in the maximum parcel displacements upward between the LL3_UL1 and LL3_UL3 simulations (Fig. 14).
The 4–7-h wMAX magnitudes were statistically significantly correlated (CC = 0.71; p value <0.05) with 4–7-h average 5–10 km AGL core updraft areas (Fig. 15), consistent with less entrainment-driven dilution associated with wider updrafts in the stronger LL shear simulations. This suggests that the wider updrafts in the simulations with stronger LL shear were more buoyant because of their decreased dilution, relative to the simulations with weaker LL shear.
e. Impact of upper-level shear variations on squall-line updraft characteristics
As discussed above, variations in LL shear led to the most drastic changes in updraft width, entrainment characteristics, B, and ultimately wMAX, whereas variations in UL shear had comparatively lesser impacts. LL3 simulations with increasing UL shear generally showed minor increases in updraft area (e.g., Figs. 6–8), PT concentrations aloft (e.g., Figs. 9 and 10), LL horizontal mass flux (e.g., Fig. 11), wMAX (e.g., Figs. 12, 13, 15), and B (e.g., Fig. 5). LL2 simulations with increasing UL shear generally had small changes in squall-line updraft characteristics. Interestingly, LL1 simulations with increasing UL shear had opposite changes in squall-line updraft characteristics compared to the LL3 simulations, including minor decreases in updraft area (e.g., Figs. 6–8), PT concentrations aloft (e.g., Figs. 9 and 10), LL horizontal mass flux (e.g., Fig. 11), wMAX (e.g., Figs. 12, 13, 15), and B (e.g., Fig. 5). Changes in wMAX from variations in UL shear for a given value of strong LL shear (i.e., our LL3 simulations) were generally comparable to the results of Coniglio et al. (2006). For example, the difference in median 3–6-h wMAX between the weak UL shear and moderate-to-strong UL shear simulations in Coniglio et al. (2006) was roughly 40% (see their Fig. 10a). Similarly, the differences in median 4–7-h wMAX between our LL3_UL1 and LL3_UL3 simulations was roughly 27%.
4. Summary, discussion, and conclusions
This study investigated the influence of LL and UL shear variations on squall-line updraft characteristics for a singular thermodynamic environment (i.e., fixed CAPE and CIN). Our simulations revealed that as LL shear increased, LL horizontal mass flux increased, leading to wider squall-line updrafts that were stronger and deeper. Changes in squall-line updraft width, entrainment, and intensity characteristics owing to variations in UL shear were relatively minor. The wider updrafts in stronger LL shear environments displayed the least amount of entrainment-driven dilution of B, consistent with larger concentrations aloft of a PT that originated from the EIL. Statistically significant correlations were found between LL horizontal mass flux (and by extension, LL shear), updraft width, PT transport aloft, and wMAX. Forward trajectory analyses revealed that the larger B in wider, less dilute updrafts resulted in larger vertical accelerations and wMAX when compared to narrower and more dilute updrafts. Thus, our results highlight a previously unexplored mechanism for how shear (especially LL shear) influences squall-line updraft intensity. In an interesting (albeit unsurprising) contrast with previous work that has established the central role of deep-layer shear in regulating supercell updraft width, our results suggest that LL shear is more important than UL shear to determining squall-line updraft width.
Alfaro’s (2017) results are particularly pertinent to our study. He showed that as the ratio of high-CAPE inflowing air to the total inflowing air increased—whether through variations in LL shear and/or LL thermodynamic variations—squall-line updrafts became less dilute, had greater B, and thus were stronger and deeper. Our work adds to this conceptual model by explaining how differences in LL shear affect updraft properties, such as updraft width, within squall lines for a given thermodynamic environment. Finally, we studied possible impacts of varying UL shear on squall line updraft properties, similar to Coniglio et al. (2006). Our results arguably agree with Coniglio et al. (2006); increasing UL shear (in our case, to 10 m s−1 over the 2.5–10 km AGL layer) led to slightly stronger squall-line updrafts in the strongest LL shear simulations. Interestingly, the weakest LL shear simulations in combination with increasing UL shear led to weaker squall-line updrafts. Further exploration into the reasons for this behavior of squall-line updrafts within various LL shear environments in combination with a much more expansive set of UL shear magnitudes (i.e., stronger) is warranted for future work.
Clearly, CAPE is strongly tied to wMAX and future simulations across a range of thermodynamic environments (i.e., various CAPE and CIN magnitudes) are necessary to evaluate the sensitivity of the aforementioned relationship to CAPE. It is also possible that differences in the cold pool propagation speed and the resultant magnitude of cold pool relative flow outside of the cold pool, in addition to vertical wind shear, influence the tilt, entrainment characteristics, and intensity of squall-line deep convective updrafts. This is yet another potential mechanism that may modulate squall-line intensity and structure that future studies should investigate. Furthermore, future simulated squall-line entrainment studies should include stable layers, drier middle- to upper-troposphere environments, stronger UL shear magnitudes, different three-dimensional wind profiles, and modified initial and lateral boundary conditions to see if these trends extend beyond the environments considered herein. We did not examine how the relationship between environmental shear, entrainment, and updraft characteristics might impact severe weather, which is left to future work.
Acknowledgments
The authors thank Chris Nowotarski, Geoff Marion, Bowen Pan, Peter Marinescu, Susan van den Heever, Rachel Phinney, Johana Lambert, and Janice Mulholland for their insightful comments, coding assistance, and stimulating conversations related to the material presented in this article. The authors would also like to thank the exceptionally helpful comments and suggestions from the three anonymous reviewers. We greatly appreciate George Bryan’s tireless efforts in maintaining CM1. We would like to acknowledge high-performance computing support from Cheyenne (doi:10.5065/D6RX99HX) provided by NCAR’s Computational and Information Systems Laboratory, sponsored by the National Science Foundation. Funding for this research was supported by National Science Foundation Grants AGS-1928666 and AGS-1841674 and Department of Energy Atmospheric System Research Grants DE-SC0020104 and DE-SC0000246356.
Data availability statement
All code, model configuration, and model output are available upon request to the corresponding author.
REFERENCES
Alfaro, D. A., 2017: Low-tropospheric shear in the structure of squall lines: Impacts on latent heating under layer-lifting ascent. J. Atmos. Sci., 74, 229–248, https://doi.org/10.1175/JAS-D-16-0168.1.
Alfaro, D. A., and M. Khairoutdinov, 2015: Thermodynamic constraints on the morphology of simulated midlatitude squall lines. J. Atmos. Sci., 72, 3116–3137, https://doi.org/10.1175/JAS-D-14-0295.1.
Bluestein, H. B., and M. H. Jain, 1985: Formation of mesoscale lines of precipitation: Severe squall lines in Oklahoma during the spring. J. Atmos. Sci., 42, 1711–1732, https://doi.org/10.1175/1520-0469(1985)042<1711:FOMLOP>2.0.CO;2.
Browning, K. A., 1963: The growth of large hail within a steady updraught. Quart. J. Roy. Meteor. Soc., 89, 490–506, https://doi.org/10.1002/qj.49708938206.
Bryan, G. H., and M. J. Fritsch, 2000: Moist absolute instability: The sixth static stability state. Bull. Amer. Meteor. Soc., 81, 1207–1230, https://doi.org/10.1175/1520-0477(2000)081<1287:MAITSS>2.3.CO;2.
Bryan, G. H., and M. J. Fritsch, 2002: A benchmark simulation for moist nonhydrostatic numerical models. Mon. Wea. Rev., 130, 2917–2928, https://doi.org/10.1175/1520-0493(2002)130<2917:ABSFMN>2.0.CO;2.
Bryan, G. H., and M. D. Parker, 2010: Observations of a squall line and its near environment using high-frequency rawinsonde launches during VORTEX2. Mon. Wea. Rev., 138, 4076–4097, https://doi.org/10.1175/2010MWR3359.1.
Bryan, G. H., and H. Morrison, 2012: Sensitivity of a simulated squall line to horizontal resolution and parameterization of microphysics. Mon. Wea. Rev., 140, 202–225, https://doi.org/10.1175/MWR-D-11-00046.1.
Bryan, G. H., and R. Rotunno, 2014: The optimal state for gravity currents in shear. J. Atmos. Sci., 71, 448–468, https://doi.org/10.1175/JAS-D-13-0156.1.
Bryan, G. H., J. C. Wyngaard, and J. M. Fritsch, 2003: Resolution requirements for the simulation of deep moist updrafts. Mon. Wea. Rev., 131, 2394–2416, https://doi.org/10.1175/1520-0493(2003)131<2394:RRFTSO>2.0.CO;2.
Bryan, G. H., J. C. Knievel, and M. D. Parker, 2006: A multimodel assessment of RKW theory’s relevance to squall-line characteristics. Mon. Wea. Rev., 134, 2772–2792, https://doi.org/10.1175/MWR3226.1.
Bryan, G. H., R. Rotunno, and J. M. Fritsch, 2007: Roll circulations in the convective region of a simulated squall line. J. Atmos. Sci., 64, 1249–1266, https://doi.org/10.1175/JAS3899.1.
Chen, Q., J. Fan, S. Hagos, W. I. Gustafson Jr., and L. K. Berg, 2015: Roles of wind shear at different vertical levels: Cloud system organization and properties. J. Geophys. Res. Atmos., 120, 6551–6574, https://doi.org/10.1002/2015JD023253.
Chong, M., P. Amayenc, G. Scialom, and J. Testud, 1987: A tropical squall line observed during the COPT 81 experiment in West Africa. Part I: Kinematic structure inferred from dual-Doppler radar data. Mon. Wea. Rev., 115, 670–694, https://doi.org/10.1175/1520-0493(1987)115<0670:ATSLOD>2.0.CO;2.
Coffer, B. E., and M. D. Parker, 2015: Impacts of increasing low-level shear on supercells during the early evening transition. Mon. Wea. Rev., 143, 1945–1969, https://doi.org/10.1175/MWR-D-14-00328.1.
Cohen, A. E., M. C. Coniglio, S. F. Corfidi, and S. J. Corfidi, 2007: Discrimination of mesoscale convective system environments using sounding observations. Wea. Forecasting, 22, 1045–1062, https://doi.org/10.1175/WAF1040.1.
Coniglio, M. C., D. J. Stensrud, and M. B. Richman, 2004: An observational study of derecho-producing convective systems. Wea. Forecasting, 19, 320–337, https://doi.org/10.1175/1520-0434(2004)019<0320:AOSODC>2.0.CO;2.
Coniglio, M. C., D. J. Stensrud, and L. J. Wicker, 2006: Effects of upper-level shear on the structure and maintenance of strong quasi-linear mesoscale convective systems. J. Atmos. Sci., 63, 1231–1252, https://doi.org/10.1175/JAS3681.1.
Coniglio, M. C., H. E. Brooks, S. J. Weiss, and S. F. Corfidi, 2007: Forecasting the maintenance of quasi-linear mesoscale convective systems. Wea. Forecasting, 22, 556–570, https://doi.org/10.1175/WAF1006.1.
Coniglio, M. C., S. F. Corfidi, and J. S. Kain, 2012: Views on applying RKW theory: An illustration using the 8 May 2009 derecho-producing convective system. Mon. Wea. Rev., 140, 1023–1043, https://doi.org/10.1175/MWR-D-11-00026.1.
Davis, C., and Coauthors, 2004: The bow echo and MCV experiment: Observations and opportunities. Bull. Amer. Meteor. Soc., 85, 1075–1094, https://doi.org/10.1175/BAMS-85-8-1075.
Deierling, W., and W. A. Petersen, 2008: Total lightning activity as an indicator of updraft characteristics. J. Geophys. Res., 113, D16210, https://doi.org/10.1029/2007JD009598.
Dennis, E. J., and M. R. Kumjian, 2017: The impact of vertical wind shear on hail growth in simulated supercells. J. Atmos. Sci., 74, 641–663, https://doi.org/10.1175/JAS-D-16-0066.1.
Doswell, C. A., H. E. Brooks, and R. A. Maddox, 1996: Flash flood forecasting: An ingredients-based methodology. Wea. Forecasting, 11, 560–581, https://doi.org/10.1175/1520-0434(1996)011<0560:FFFAIB>2.0.CO;2.
Durran, D. R., and J. B. Klemp, 1983: A compressible model for the simulation of moist mountain waves. Mon. Wea. Rev., 111, 2341–2361, https://doi.org/10.1175/1520-0493(1983)111<2341:ACMFTS>2.0.CO;2.
Evans, J. S., and C. A. Doswell III, 2001: Examination of derecho environments using proximity soundings. Wea. Forecasting, 16, 329–342, https://doi.org/10.1175/1520-0434(2001)016<0329:EODEUP>2.0.CO;2.
Fovell, R. G., and Y. Ogura, 1989: Effect of vertical wind shear on numerically simulated multicell storm structure. J. Atmos. Sci., 46, 3144–3176, https://doi.org/10.1175/1520-0469(1989)046<3144:EOVWSO>2.0.CO;2.
Fovell, R. G., and P. S. Dailey, 1995: The temporal behavior of numerically simulated multicell-type storms. Part I: Modes of behavior. J. Atmos. Sci., 52, 2073–2095, https://doi.org/10.1175/1520-0469(1995)052<2073:TTBONS>2.0.CO;2.
Fovell, R. G., and P.-H. Tan, 1998: The temporal behavior of numerically simulated multicell-type storms. Part II: The convective cell life cycle and cell regeneration. Mon. Wea. Rev., 126, 551–577, https://doi.org/10.1175/1520-0493(1998)126<0551:TTBONS>2.0.CO;2.
French, A. J., and M. D. Parker, 2010: The response of simulated nocturnal convective systems to a developing low-level jet. J. Atmos. Sci., 67, 3384–3408, https://doi.org/10.1175/2010JAS3329.1.
Grady, R. L., and J. Verlinde, 1997: Triple-Doppler analysis of a discretely propagating, long-lived, High Plains squall line. J. Atmos. Sci., 54, 2729–2748, https://doi.org/10.1175/1520-0469(1997)054<2729:TDAOAD>2.0.CO;2.
Hane, C. E., 1973: The squall line thunderstorm: Numerical experimentation. J. Atmos. Sci., 30, 1672–1690, https://doi.org/10.1175/1520-0469(1973)030<1672:TSLTNE>2.0.CO;2.
Hernandez-Deckers, D., and S. C. Sherwood, 2016: A numerical investigation of cumulus thermals. J. Atmos. Sci., 73, 4117–4136, https://doi.org/10.1175/JAS-D-15-0385.1.
Hernandez-Deckers, D., and S. C. Sherwood, 2018: On the role of entrainment in the fate of cumulus thermals. J. Atmos. Sci., 75, 3911–3924, https://doi.org/10.1175/JAS-D-18-0077.1.
Hitchcock, S. M., R. S. Schumacher, G. R. Herman, M. C. Coniglio, M. D. Parker, and C. L. Ziegler, 2019: Evolution of pre- and postconvective environmental profiles from mesoscale convective systems during PECAN. Mon. Wea. Rev., 147, 2329–2354, https://doi.org/10.1175/MWR-D-18-0231.1.
Houze, R. A., Jr., 1989: Observed structure of mesoscale convective systems and implications for large-scale heating. Quart. J. Roy. Meteor. Soc., 115, 425–461, https://doi.org/10.1002/qj.49711548702.
James, R. P., J. M. Fritsch, and P. M. Markowski, 2005: Environmental distinctions between cellular and slabular convective lines. Mon. Wea. Rev., 133, 2669–2691, https://doi.org/10.1175/MWR3002.1.
James, R. P., P. M. Markowski, and J. M. Fritsch, 2006: Bow echo sensitivity to ambient moisture and cold pool strength. Mon. Wea. Rev., 134, 950–964, https://doi.org/10.1175/MWR3109.1.
Jorgensen, D. P., M. A. LeMone, and S. B. Trier, 1997: Structure and evolution of the 22 February 1993 TOGA COARE squall line: Aircraft observations of precipitation, circulation, and surface energy fluxes. J. Atmos. Sci., 54, 1961–1985, https://doi.org/10.1175/1520-0469(1997)054<1961:SAEOTF>2.0.CO;2.
Kingsmill, D. E., and R. A. Houze Jr., 1999: Kinematic characteristics of air flowing into and out of precipitating convection over the west Pacific warm pool: An airborne Doppler radar survey. Quart. J. Roy. Meteor. Soc., 125, 1165–1207, https://doi.org/10.1002/qj.1999.49712555605.
Lane, T. P., and M. W. Moncrieff, 2010: Characterization of momentum transport associated with organized moist convection and gravity waves. J. Atmos. Sci., 67, 3208–3225, https://doi.org/10.1175/2010JAS3418.1.
Lang, T. J., and S. A. Rutledge, 2002: Relationships between convective storm kinematics, precipitation, and lightning. Mon. Wea. Rev., 130, 2492–2506, https://doi.org/10.1175/1520-0493(2002)130<2492:RBCSKP>2.0.CO;2.
Lebo, Z. J., and H. Morrison, 2015: Effects of horizontal and vertical grid spacing on mixing in simulated squall lines and implications for convective strength and structure. Mon. Wea. Rev., 143, 4355–4375, https://doi.org/10.1175/MWR-D-15-0154.1.
Lombardo, K. A., and B. A. Colle, 2012: Ambient conditions associated with the maintenance and decay of quasi-linear convective systems crossing the northeastern U.S. coast. Mon. Wea. Rev., 140, 3805–3819, https://doi.org/10.1175/MWR-D-12-00050.1.
Mahoney, K. M., G. M. Lackmann, and M. D. Parker, 2009: The role of momentum transport in the motion of a quasi-idealized mesoscale convective system. Mon. Wea. Rev., 137, 3316–3338, https://doi.org/10.1175/2009MWR2895.1.
Marinescu, P. J., S. C. van den Heever, S. M. Saleeby, S. M. Kreidenweis, and P. J. De Mott, 2017: The microphysical roles of lower-tropospheric versus midtropospheric aerosol particles in mature-stage MCS precipitation. J. Atmos. Sci., 74, 3657–3678, https://doi.org/10.1175/JAS-D-16-0361.1.
Marion, G. R., and R. J. Trapp, 2019: The dynamical coupling of convective updrafts, downdrafts, and cold pools in simulated supercell thunderstorms. J. Geophys. Res. Atmos., 124, 664–683, https://doi.org/10.1029/2018JD029055.
Marion, G. R., and R. J. Trapp, 2021: Controls of quasi-linear convective system tornado intensity. J. Atmos. Sci., 78, 1189–1205, https://doi.org/10.1175/JAS-D-20-0164.1.
Markowski, P. M., and Y. P. Richardson, 2014: The influence of environmental low-level shear and cold pools on tornadogenesis: Insights from idealized simulations. J. Atmos. Sci., 71, 243–275, https://doi.org/10.1175/JAS-D-13-0159.1.
Mechem, D. B., R. A. Houze Jr., and S. S. Chen, 2002: Layer inflow into precipitating convection over the western tropical Pacific. Quart. J. Roy. Meteor. Soc., 128, 1997–2030, https://doi.org/10.1256/003590002320603502.
Moncrieff, M. W., and C. Liu, 1999: Convection initiation by density currents: Role of convergence, shear, and dynamical organization. Mon. Wea. Rev., 127, 2455–2464, https://doi.org/10.1175/1520-0493(1999)127<2455:CIBDCR>2.0.CO;2.
Morrison, H., 2017: An analytic description of the structure and evolution of growing deep cumulus updrafts. J. Atmos. Sci., 74, 809–834, https://doi.org/10.1175/JAS-D-16-0234.1.
Morrison, H., and J. M. Peters, 2018: Theoretical expressions for the ascent rate of moist convective thermals. J. Atmos. Sci., 75, 1699–1719, https://doi.org/10.1175/JAS-D-17-0295.1.
Morrison, H., G. Thompson, and V. Tatarskii, 2009: Impact of cloud microphysics on the development of trailing stratiform precipitation in a simulated squall line: Comparison of one- and two-moment schemes. Mon. Wea. Rev., 137, 991–1007, https://doi.org/10.1175/2008MWR2556.1.
Morrison, H., J. M. Peters, A. C. Varble, W. M. Hannah, and S. E. Giangrande, 2020: Thermal chains and entrainment in cumulus updrafts. Part I: Theoretical description. J. Atmos. Sci., 77, 3637–3660, https://doi.org/10.1175/JAS-D-19-0243.1.
Morrison, H., J. M. Peters, and S. C. Sherwood, 2021: Comparing growth rates of simulated moist and dry convective thermals. J. Atmos. Sci., 78, 797–816, https://doi.org/10.1175/JAS-D-20-0166.1.
Morton, B. R., G. I. Taylor, and J. S. Turner, 1956: Turbulent gravitational convection from maintained and instantaneous sources. Proc. Roy. Soc. London, 234A, 1–23, https://doi.org/10.1098/rspa.1956.0011.
Mulholland, J. P., S. W. Nesbitt, R. J. Trapp, and J. M. Peters, 2020: The influence of terrain on the convective environment and associated convective morphology from an idealized modeling prospective. J. Atmos. Sci., 77, 3929–3949, https://doi.org/10.1175/JAS-D-19-0190.1.
Nielsen, E. R., and R. S. Schumacher, 2018: Dynamical insights into extreme short-term precipitation associated with supercells and mesovortices. J. Atmos. Sci., 75, 2983–3009, https://doi.org/10.1175/JAS-D-17-0385.1.
Parker, M. D., 2010: Relationship between system slope and updraft intensity in squall lines. Mon. Wea. Rev., 138, 3572–3578, https://doi.org/10.1175/2010MWR3441.1.
Parker, M. D., and R. H. Johnson, 2000: Organizational modes of midlatitude mesoscale convective systems. Mon. Wea. Rev., 128, 3413–3436, https://doi.org/10.1175/1520-0493(2001)129<3413:OMOMMC>2.0.CO;2.
Parker, M. D., and R. H. Johnson, 2004a: Structures and dynamics of quasi-2D mesoscale convective systems. J. Atmos. Sci., 61, 545–567, https://doi.org/10.1175/1520-0469(2004)061<0545:SADOQM>2.0.CO;2.
Parker, M. D., and R. H. Johnson, 2004b: Simulated convective lines with leading precipitation. Part I: Governing dynamics. J. Atmos. Sci., 61, 1637–1655, https://doi.org/10.1175/1520-0469(2004)061<1637:SCLWLP>2.0.CO;2.
Parker, M. D., and R. H. Johnson, 2004c: Simulated convective lines with leading precipitation. Part II: Evolution and maintenance. J. Atmos. Sci., 61, 1656–1673, https://doi.org/10.1175/1520-0469(2004)061<1656:SCLWLP>2.0.CO;2.
Peters, J. M., 2016: The impact of effective buoyancy and dynamic pressure forcing on vertical velocities within two-dimensional updrafts. J. Atmos. Sci., 73, 4531–4551, https://doi.org/10.1175/JAS-D-16-0016.1.
Peters, J. M., W. Hannah, and H. Morrison, 2019a: The influence of vertical wind shear on moist thermals. J. Atmos. Sci., 76, 1645–1659, https://doi.org/10.1175/JAS-D-18-0296.1.
Peters, J. M., C. J. Nowotarski, and H. Morrison, 2019b: The role of vertical wind shear in modulating maximum supercell updraft velocities. J. Atmos. Sci., 76, 3169–3189, https://doi.org/10.1175/JAS-D-19-0096.1.
Peters, J. M., C. J. Nowotarski, and G. L. Mullendore, 2020a: Are supercells resistant to entrainment because of their rotation? J. Atmos. Sci., 77, 1475–1495, https://doi.org/10.1175/JAS-D-19-0316.1.