Proper prediction of the inflow layer of deep convective storms is critical for understanding their potential updraft properties and likelihood of producing severe weather. In this study, an existing forecast metric known as the effective inflow layer (EIL) is evaluated with an emphasis on its performance for supercell thunderstorms, where both buoyancy and dynamic pressure accelerations are common. A total of 15 idealized simulations with a range of realistic base states are performed. Using an array of passive fluid tracers initialized at various vertical levels, the proportion of simulated updraft core air originating from the EIL is determined. Results suggest that the EIL metric performs well in forecasting peak updraft origin height, particularly for supercell updrafts. Moreover, the EIL metric displays consistent skill across a range of updraft core definitions. The EIL has a tendency to perform better as convective available potential energy, deep-layer shear, and EIL depth are increased in the near-storm environment. Modifications to further constrain the EIL based on the most-unstable parcel height or storm-relative flow may lead to marginal improvements for the most stringent updraft core definitions. Finally, effects of the near-storm environment on low-level and peak updraft forcing and intensity are discussed.
The properties of the air within the inflow to supercell updrafts strongly influence vertical accelerations, vertical velocities, and consequently the ability of an updraft to produce severe weather and flash flooding. Sounding-derived parameters from the near-storm environment have long been used to predict various attributes of the inflow to a storm (e.g., Showalter and Fulks 1943; Fawbush and Miller 1954; Beebe 1958; Maddox 1976; Kerr and Darkow 1996; Rasmussen and Blanchard 1998; Thompson et al. 2003; Davies 2004). Most of these parameters use the thermodynamic characteristics of a level or layer from a sounding that is assumed to be representative of the inflow layer to calculate convective available potential energy (CAPE) and convective inhibition (CIN), and relate these CAPE and CIN calculations to the potential for a convective updraft and the peak strength of an updraft if it were to form. The simplest method is a surface-based parcel, where the temperature and dewpoint are taken from a surface observation (or lowest model grid level). Slightly more complicated mixed-layer or most-unstable parcel methods recognize that deep convective updrafts ingest air originating from a deeper layer, and that the characteristics at the surface are often unrepresentative of that layer. Such methods emphasize the effects of buoyancy on updraft intensity, neglecting the potential contributions of vertical perturbation pressure gradient forces (VPPGFs) to vertical velocity, particularly in supercells (e.g., Rotunno and Klemp 1982, 1985; Weisman and Rotunno 2000). In the context of mesoscale convective systems (MCSs), where the presence of a strong cold pool is integral to the self-organization of convection through forced lifting of inflow air, recent studies have developed and evaluated layer-lifting indices that incorporate the interaction between cold pools and ambient lower-tropospheric shear in determining the potential buoyancy of inflow in deep convection (e.g., Alfaro and Khairoutdinov 2015; Alfaro 2017; Alfaro and Coniglio 2018).
In response to the limitations of arbitrary, fixed estimations of the inflow layer, Thompson et al. (2007, hereafter T07) developed a method for more precisely determining what they refer to as the effective inflow layer (EIL; herein EIL refers specifically to the T07 parameter or our alterations to it) to a storm. This parameter was formulated based on their stated assumption that “only lifted parcels associated with CAPE will sustain a deep thunderstorm updraft, whereas parcels associated with very large CIN and/or negligible CAPE will ultimately result in storm demise.” Using a database of 835 Rapid Update Cycle (RUC) proximity soundings, T07 tested several combinations of CAPE and CIN thresholds to define the EIL, concluding that the definition that resulted in a nonzero EIL depth in the largest number of cases (i.e., 95%) was most desirable. They settled on the following definition that is widely used today: the contiguous layer wherein lifted parcels would have at least 100 J kg−1 of CAPE and CIN < −250 J kg−1. More stringent thresholds of CAPE and CIN led to a nonzero EIL depth in a fewer number of cases, and these alternate criteria were thus discarded. Population EIL depths were relatively consistent (median values between 1400 and 1800 m) between supercells and nonsupercell environments, but T07 identified a slight preference toward deeper EILs for significantly tornadic supercells over nonsupercellular convection.
Bunkers et al. (2002) highlight the importance of choosing a representative level or layer in computing both CAPE and shear for convective storms forecasting. Thus, it is not surprising that parameters using the EIL or integrated over the depth of the EIL have been developed to provide insight into the characteristics or potential severity of deep convection. Sounding-derived measures of vertical wind shear or storm-relative helicity (SRH) are generally relied on to predict convective organization and rotation, with the implicit assumption that greater organization and rotation may increase updraft strength through contributions from VPPGFs. T07 improved on fixed layer (e.g., 0–1 or 0–3 km) SRH by integrating SRH over the depth of the EIL, showing that the effective storm relative helicity (ESRH) more clearly discriminated between significantly tornadic and nontornadic supercells. T07 also showed that using the effective bulk wind difference as a measure of deep-layer shear better discriminated between tornadic and nontornadic supercells than the 0–6-km bulk wind difference. In recognition of these improvements, EIL-derived parameters have been included in updated versions of the supercell composite parameter and significant tornado parameter as well as numerous other studies, sounding-derived parameters, and prediction systems (e.g., Smith et al. 2012; Cintineo et al. 2014; Bunkers et al. 2014; Sherburn and Parker 2014; Davenport and Parker 2015; Van Den Broeke 2016; Hart and Cohen 2016; Gallo et al. 2018). Indeed, commonly used sounding analysis software such as the Sounding and Hodograph Analysis and Research Program in Python [SHARPPy; Blumberg et al. (2017)] now include numerous EIL-based parameters, including CAPE based on an EIL-averaged parcel. Despite the increasing ubiquity of EIL-based parameters, uncertainty remains as to its best applications. For instance, Coffer et al. (2019) show shallow-layer (500 m) SRH better discriminates between tornadic and nontornadic supercells than ESRH, but ESRH or other EIL-derived parameters may still be useful in convective hazards more directly related to peak updraft strength (e.g., hail, flash flooding, and lightning).
The assumption underlying the T07 approach to defining the inflow layer recognizes the primary importance of buoyancy in dictating updraft strength in deep convection. However, this approach implicitly allows for the potential for VPPGFs in low-level updrafts, as parcels with −250–0 J kg−1 of CIN require a force other than buoyancy to lift them to their levels of free convection (LFC). Indeed Nowotarski et al. (2011), found that in simulations of supercells with EIL bases above the surface, midlevel supercell updrafts contained parcels originating from within shallow near-surface stable layers. Within the T07 study only 55 of the 835 total supercells featured an EIL base above the surface. Because of this small sample size of potentially elevated storms relative to surface-based storms, it is comparatively less clear how relevant the EIL definition is for elevated storms. Given the exclusive dependence of the EIL parameter on buoyancy, it also seems reasonable to contemplate how well the EIL approximates the true inflow layer of supercell thunderstorms that often feature particularly strong low-level VPPGFs—particularly in situations where storms may be elevated over a surface-based layer of stable air. For instance, there is the possibility that layers not meeting the EIL criteria actively participate in updrafts given strong dynamic lifting. Alternatively, some parcels that meet the EIL criteria, but do not experience adequate upward VPPGFs may not overcome their CIN and participate in the updraft above their LFC. In the latter possibility, vertical wind shear and SRH metrics that depend on the EIL may be overestimated such that storm severity is over predicted. Finally, air entrained from above the EIL may obtain adequate buoyancy through mixing such that it participates in the updraft core.
The uncertainties described above motivate a verification study of the EIL, particularly for situations where elevated supercells are likely. This could better inform the use of EIL-derived parameters as forecasting tools. While an observation-based verification would be ideal, the lack of direct observations and limitations of dual-Doppler datasets in determining updraft parcel trajectories with adequate detail prevent a robust verification in a range of environments. Thus, we turn to numerical simulations of supercell and nonsupercell convection in a variety of real kinematic and thermodynamic environments to test the performance of the EIL parameter. Specifically, we address the following research questions:
How accurate is the EIL formulation of T07 in predicting the origin heights of air in supercell and nonsupercell updrafts cores?
How does the EIL performance vary as the definition of updraft core is changed?
Which environmental parameters best predict EIL performance, and more broadly, low-level and peak updraft forcing and intensity in real environments?
and the corresponding hypotheses:
The EIL definition is more accurate in the case of supercells than in the case of nonsupercells because of the comparatively large midtropospheric entrainment rates in the case of nonsupercells (e.g., Peters et al. 2019, 2020b).
The probability of detection of the updraft core composition by the EIL increases with the magnitude of the vertical velocity threshold used to define the updraft core.
The EIL probability of detection improves as environmental parameters related to updraft intensity, such as CAPE and vertical wind shear, increase.
a. Simulation base-state selection
To answer our research questions, a series of 15 idealized, horizontally homogeneous numerical simulations were performed using CM1 version 19.8 (Bryan and Fritsch 2002). To assess the EIL over a range of supercell and nonsupercell storm environments, we used the T07 database of RUC proximity soundings. Our goal was to generate an archetypal base-state temperature and wind profile for the significantly tornadic supercell (ST), weakly tornadic supercell (WT), nontornadic supercell (NT), and nonsupercell (NS) storm types included in the T07 database. Rather than randomly selecting soundings of each storm type, we first classified the database using the four-variable (temperature, dewpoint temperature, u-wind component, υ-wind component) self-organizing map (SOM) method described by Nowotarski and Jones (2018). A 16-node SOM, considering vertical profiles of all four variables derived from all RUC proximity soundings in the database, including nonsupercell storms was constructed, and we calculated the percent of each storm type matching each resulting node.
Our initial plan was to use the average vertical profile of all storms (regardless of type) matching the node with the highest percentage of each storm type (relative to other storm types within that node) to represent that storm type in our range of simulation base states. However, the node-averaged soundings generally did not result in sustained deep convection in our simulations. We attribute this failure to considerable low- to midlevel tropospheric dryness in the node averages (translucent profiles in Fig. 1) as compared with the Weisman and Klemp (1982) profile and RUC soundings that ultimately did result in sustained deep convection. Thus, we ultimately selected one RUC profile associated with each of the four storm types that matched the node where that storm type was most common. For instance, for our ST sounding, we selected an ST RUC profile that matched the node with the highest percentage of ST cases. This selection process was semirandom; however, care was taken to ensure that the characteristics of the selected profile closely matched the pattern identified by that node and that the selected profile could sustain deep convection in a simulation, even with the modifications discussed below. In practice, this led to the selection of RUC profiles that had considerably more low- to midlevel tropospheric moisture than the node averages, but similar hodograph shapes. This resulted in one base state for each storm type, which is relatively representative of the node average and other profiles matching that node except for the aforementioned moisture differences (cf. node average profiles to real-case RUC profiles in Fig. 1). In addition, we also used the analytic Weisman and Klemp (1982) thermodynamic and wind profiles (WK) as a base state to compare with the realistic profiles.
This process resulted in five base-state thermodynamic and kinematic soundings (Fig. 1). Despite a range of EIL depths, most-unstable parcel heights, and resulting CAPE, CIN, and shear parameters consistent with the resulting storm types, each of the unmodified base states had large surface-based CAPE (Table 1) and a surface-based EIL. This, in part, was due to the large percentage (93%) of surface-based storms included in the T07 database. Thus, to fully test the EIL for elevated supercells (and nonsupercells), each base-state thermodynamic profile was modified in order to stabilize the low-level environment and raise the EIL base. First, a shallow (S) surface-based inversion with a temperature lapse rate of −20°C km−1 was applied over a depth of 300 m to each base state (Fig. 2b). Then, in a deep (D) inversion set of simulations, the same lapse rate was applied over a depth of 600 m (Fig. 2c). To avoid supersaturation in the cases where the surface temperature modification resulted in a modified temperature below the original dewpoint temperature at that layer, the dewpoint temperature was set to 1°C less than the modified temperature.
Table 1 shows relevant thermodynamic and kinematic sounding-derived parameters for all 15 simulation base states.1 In general, the WK base-state characteristics fall somewhere between those of the corresponding ST and WT base states. For all storm types, the addition of the near-surface stable layers reduces or eliminates surface-based CAPE such that the EIL base raises from the surface, with the highest EIL base at 400 m. The most-unstable (MU) parcel height ranges from the surface to 1304 m, with varying positions within the EIL across simulations. Deep-layer shear (represented by the 0–6 km bulk wind difference) is relatively strong for all but the NS type, and SRH decreases progressively from ST to NS types. As a consequence of the reduction in the EIL thickness, ESRH decreases as the depth of the stable layer increases. Finally, both the supercell composite parameter (SCP) and the significant tornado parameter (STP) generally decrease when transitioning from ST to NS storm types and as the depth of the stable layer is increased. In summary, a reasonably wide range of base-state thermodynamic and kinematic environments supportive of supercell and nonsupercellular convection resulted from the 15 base states, including a relatively wide range of EIL bases, tops, depths and MU parcel levels relative to the EIL. Thus, the performance of the EIL can be tested over a realistic range of environments.
b. Numerical model configuration
Each of the storm environments described above serves as the horizontally homogeneous base state for a 2-h simulation. The simulations were performed on a 120 km × 120 km × 20 km domain with horizontal grid spacing of 250 m and a vertical grid spacing of 50 m below 4 km stretched to a grid spacing of 250 m above 13 km. The upper and lower boundary conditions were rigid and free slip with Rayleigh damping applied above 15 km and open-radiative horizontal boundary conditions. Simulations were performed with CM1 in the large-eddy simulation configuration with the Deardorff (1980) subgrid-scale TKE scheme, fifth-order advection, and a time-splitting pressure solver with acoustic wave terms solved on smaller time steps (Klemp and Wilhelmson 1978). The Coriolis force, radiation, and surface fluxes of heat, moisture, and momentum are omitted. Precipitation microphysics were parameterized using the NSSL two-moment scheme including both graupel and hail (Mansell et al. 2010). In each simulation, deep convection was initialized at the center of the domain using an updraft nudging technique following Naylor and Gilmore (2012). The only changes in the updraft nudging from Naylor and Gilmore (2012) were a modified vertical position (z = 1000 m) and vertical radius (1000 m) of the spheroid over which the updraft forcing is applied as well as reduction of the period of forcing to the first 1200 s of the model integration. Grid translation was applied to each simulation based on the unique storm motion resulting from each base-state wind profile (Fig. 1) in order to keep the resulting deep convection approximately centered in the domain.
c. Passive tracer method for updraft parcel origin analysis
To track the vertical motion of air originating at various layers, each simulation was initialized with an array of passive tracers. For each 100-m-thick layer below 5 km and each 500-m-thick layer above 5 km, a unique passive tracer was defined with the initial concentration set to 100% within that layer and 0% at all heights outside that layer. Thus, the initial condition was stratified with 80 unique passive fluid tracers, which were advected by the simulation as it integrated. Consequently, at any time in the simulation, the concentration of each passive tracer at any grid point represents the amount of air at that grid point originating from the corresponding layer in the base-state. By comparing the tracer origin heights to the EIL, the proportion of air at any location in the domain originating within, below, or above the EIL can then be calculated. More specifically, if the midpoint height of the passive tracer layer was greater than or equal to the bottom height of the EIL or less than or equal to the EIL top height, that tracer was considered part of the EIL.
Next, the updraft core of each simulated supercell was identified by locating the maximum vertical velocity associated with the dominant right-moving supercell at each simulation time after 60 min (for the nonsupercell simulations, the strongest vertical velocity within the domain was used). A 30 km × 30 km subdomain centered on each maximum updraft coordinate was defined, then for each vertical level, any grid point with vertical velocity exceeding a specified percentile of vertical velocity over 0.1 m s−1 was classified as part of the updraft core at that level. By comparing concentrations of tracers from each level at these grid points for various updraft percentile thresholds used to define the updraft core and any height within the updraft, we can determine the accuracy of the T07 EIL definition in predicting the composition of the updraft core at that level. For most subsequent analysis, passive tracer (i.e., updraft origin level) concentrations were averaged over the 60–90-min period of each simulation.2
d. Calculation of vertical perturbation pressure gradient forcing terms
In addition to buoyancy, vertical accelerations due to vertically oriented dynamic perturbation pressure gradients (DPA) have been shown to enhance supercell updrafts (e.g., Weisman and Klemp 1984; Weisman and Rotunno 2000), whereas vertical accelerations due to a vertically oriented buoyancy perturbation pressure gradient (BPA) are often downward oriented and hinder updraft intensity. These acceleration terms may be understood in terms of the anelastic vertical momentum equation:
where B ≡ −g(ρ′/ρ0) − gqi is buoyancy (where qi is the mixing ratio of the ith hydrometeor species) and DPA is the sum of the vertical acceleration terms due to linear dynamic pressure perturbations (LDPA) and nonlinear dynamic pressure perturbations (NLDPA). The decomposed pressure terms in Eq. (1) are defined as follows:
where ρ0(z) is the density from the initial model sounding, V0 is the horizontal wind from the initial model sounding, V′ ≡ V − V0 is the storm modified wind field and V is the total wind field, is referred to as “buoyancy pressure perturbation,” is referred to as “nonlinear dynamic pressure perturbation,” and is referred to as “linear dynamic pressure perturbation.” All pressure perturbations and their associated forcing terms were calculated using the in-line code and direct output variables now available in CM1. With open radiative boundary conditions, as set here, CM1 uses the method of images in a horizontal Fourier transform to enforce zero-gradient lateral boundary conditions.
Both DPA and BPA are also thought to be influenced by properties of the EIL. For instance, BPA is primarily dependent on updraft buoyancy, which itself results from the release of CAPE. DPA, on the other hand, is connected to complex flow interactions between the environment (including the inflow) characteristics and the updraft, which include (but are not limited to) the tilting of ambient inflow horizontal vorticity into the vertical direction and the associated dynamic pressure perturbations. In the context of supercell updrafts, DPA is often associated with the vertical distribution of vertical vorticity and attendant nonlinear dynamic negative pressure perturbation, wherein a mesocyclone may enhance updrafts below the level of maximum vertical vorticity (e.g., Rotunno and Klemp 1982, 1985; Weisman and Rotunno 2000). Because of these connections between the inflow layer, vertical accelerations, and updraft speed, it is important that we understand what background environmental and cloud-scale processes determine the extent of a storm’s inflow layer.
3. Results and discussion
All 15 simulations employing the base states described above resulted in sustained deep convection. The unmodified supercell base states (WK, ST, WT, and NT) all resulted in a relatively steady, right-moving cyclonic supercell with generally one large updraft and at least transient, if not permanent, reflectivity features characteristic of supercells (e.g., hook echoes, inflow notches, forward flank reflectivity cores; Figs. 3a,d,g,j). In general, as the stable layers were added to the base states, convection tended to become weaker and less obviously supercellular (cf. Figs. 3d,e,f), though the WK suite of simulations remained quite strong with well-defined supercellular characteristics. Moreover, the peak updraft cores in the WK simulations tended to be more coherent in the sense that the 95th percentile of vertical velocity was confined to one contiguous region with relatively smooth edges, particularly when compared with the WT suite of simulations (cf. Figs. 3a,g). The NS base states resulted in multicellular convection with clusters of more transient, smaller updraft cores (Fig. 3m).
a. Comparison of EIL performance for varied storm types
To address our first research question, we examine the distribution of the origins of air within the updraft core at each altitude for each simulation as compared with the sounding-derived EIL for each simulation (Fig. 4). For an updraft core threshold definition of the 95th percentile of vertical velocity exceeding 0.1 m s−1 at each height, it is clear that the vast majority of air within the updraft core above ~2 km originates within the T07 EIL for the WK, ST, and NT sets of simulations, even when the EIL is elevated (Figs. 4a–f,j,k,l). The notable exceptions are the WT and NS classes of simulations, which generally have weaker updrafts and much lower levels of peak vertical velocity (indeed, many of these updrafts do not reach the tropopause; Figs. 4h,i,m,n,o). While much of the updraft core air in these simulations originates within the EIL, a considerable fraction originates above the EIL, which is representative of entrainment of midtropospheric air into the updraft core. It also appears that the T07 EIL may under predict depth of the true inflow layer for the NS storm types given the extension of the simulated inflow layer above the EIL top to 2 km (Figs. 4m,n,o).
While, at least qualitatively, the EIL appears to capture the majority of air participating in many of the updrafts, it appears to “overforecast” the inflow layer in many supercell cases, where significant portions of the T07 EIL do not participate in the updraft core (e.g., Fig. 4a). It seems that in many simulations, the true inflow layer is focused on a relatively narrow layer centered on or near the MU parcel layer (Figs. 4a–i). There is also evidence of low-level updrafts “recycling” outflow air in some simulations where low-level (below 2 km) updrafts have origins extending above 4 km (Figs. 4b,c,j–l).
To clarify the physical features leading to differences in updraft origin distributions, Fig. 5 illustrates a stark contrast between the updraft characteristics and airflow between the nonsupercell (left column, NSD) and typical3 supercell (right column, WKD) with a surface inversion. First, the NSD updraft is comprised of a series of discrete rising thermals with relatively weak (w < 20 m s−1) maximum vertical velocity, whereas the WKD supercell updraft is a wide, continuous, plume of upward motion with a relatively strong (w > 50 m s−1) peak updraft (cf. Figs. 5a,b). These results are consistent with the structural differences found between supercells and nonsupercells in Peters et al. (2020b).
Peters et al. (2019) found that supercell peak updrafts are greater than those of nonsupercells largely because of their enhanced width in regimes with greater low-level, storm-relative flow, which protects the updraft core (composed of positively buoyant air originating in the EIL) from entrainment of free tropospheric air, resulting in greater core buoyancy. Despite variations in CAPE between simulations, our results support this finding. The lower thermals of the NSD updraft are comprised of air that is 60%–80% from within the EIL (Fig. 5c), but this fraction falls below 50% for the upper, weaker thermals that have a significant fraction of air entrained from the 4–6-km layer above the EIL (Fig. 5a). As might be expected, there is relatively little DPA in the NSD updraft, with only small areas of DPA mostly associated with the dynamic low pressure perturbations that form in toroidal updraft circulations in the middle troposphere consistent with Morrison and Peters (2018). Thus, virtually no air from below the EIL participates in this updraft (Fig. 5e). The supercell (WKD) updraft, by contrast, is composed of nearly pure (>90%) EIL air as far aloft as 11 km (Fig. 5d). On the WKD updraft periphery at mid- to upper levels, particularly on the eastern flank, some air is ingested from above the EIL, but this air does not participate in the updraft core (Fig. 5b). Finally, the lower portion of the WKD updraft is collocated with widespread areas of positive DPA, which is associated with lofting of air from below the EIL as far as 3–4 km aloft, particularly in the outflow region of the supercell, but this air does not participate in the peak updraft farther aloft (Fig. 5f). Below the EIL, BPA is also generally positive in both simulations, but particularly at the leading edge of the WKD outflow (Fig. 5d), providing additional upward forcing for air from below the EIL and within the cold pool, consistent with results of Brown and Nowotarski (2019). BPA is generally negative within the WKD updraft core above 1.5 km.
While Fig. 4 provides a qualitative analysis of the T07 EIL accuracy, we perform a more quantitative analysis of updraft origins relative to the predicted EIL regardless of if the EIL is surface based or elevated for each storm type. Figure 6 shows the percentage of the 95th percentile updraft core originating within, below, and above the EIL at each height averaged over the 60–90-min period of the three simulations for each storm type. First, the WK simulations (Fig. 6a) show the strongest updraft core of any storm type (cf. dashed black lines with other panels), consistent with having the largest magnitudes of 0–6-km shear (Table 1). The WK simulations also have the smallest proportion of air originating above the EIL, with nearly pure EIL air composing the updraft core above 2 km (cf. the red line values to black line values for each panel). The ST simulations (Fig. 6b) show significantly more entrainment of air from above the EIL into the updraft core than the WK simulations (cf. red lines in Figs. 6a,b), and a somewhat weaker core updraft, despite having comparable base-state CAPE and deep-layer shear with significantly greater SRH. Indeed, all the supercell simulations based on real environments (ST, WT, NT) have significantly greater fractions of above-EIL air within the updraft core than the WK simulations, suggesting that the WK analytic sounding may result in unrealistically small entrainment (also evident in Figs. 4a–c). The WT simulations have significantly weaker updraft cores and greater above-EIL entrainment than the other supercells (Fig. 6c), which is perhaps unsurprising given that these base states have the weakest CAPE and deep-layer shear of all the supercell simulations, including the NT supercells (Fig. 6d). As might be expected from the analysis of Fig. 5, the NS simulations (Fig. 6e) exhibit considerable entrainment of above-EIL air; in fact, more air within the core updraft above 4 km comes from above the EIL than within it (e.g., the black and red lines have similar values between 4 and 8 km).
Regardless of storm type, air from below the EIL does not form a significant fraction of the core updraft above 2–3 km. This is not surprising, given this air is excluded from the T07 EIL metric precisely because it does not have appreciable CAPE. Nowotarski et al. (2011) showed that some of the air entering the 4-km updraft in their simulations originated within low-level stable layers, but their analysis did not examine whether this air had appreciable CAPE, and they adopted a weak vertical velocity threshold for their definition of “updraft.” This analysis suggests that below-EIL air does not meaningfully contribute to the updraft core in the mid- to upper troposphere. However, the fraction of the low-level (1 km) updraft composed of below-EIL air varies between storm types. The WK and tornadic supercells loft below-EIL air higher than the NT and NS simulations, which would be consistent with stronger low-level dynamic lifting of this negatively buoyant air, hypothesized to be critical for tornadogenesis.
We can consider the T07 EIL metric as a forecast of the air that will participate in the updraft core. Within this forecast framework, the percentage of updraft core air that originates from within the sounding-derived EIL is equivalent to the probability of detection (POD; i.e., nonzero tracer concentration within the red lines on Fig. 4 and black lines in Fig. 6), and the percentage of the sounding-derived EIL that does not meaningfully participate in the updraft core4 is equivalent to the false alarm ratio (FAR; i.e., annotated white space in Fig. 4a and gray lines in Fig. 6). These results suggest that, at least for a definition of updraft core using the 95th percentile of vertical velocity as the threshold, the POD of the EIL is limited by entrainment of air from above the EIL, and not air from below it (Fig. 7d). It should be noted that this forecast framework (e.g., using POD and FAR concepts) is intended to provide an objective analysis of EIL performance rather than an explicit measure of practical forecast skill. For instance, we recognize that the T07 EIL parameter was never intended to account for entrainment and forecasters do not use it to do so; however, entrained air that enters the updraft core still serves to dilute the updraft core and alter its properties from those that might be expected if the EIL were the only origins of updraft core air. For this updraft-core threshold (Fig. 8d), the POD for the peak (i.e., maximum vertical velocity at any level) updraft core ranges from 40% to 99%, with the poorest performance (i.e., greatest entrainment) in the nonsupercell simulations. Alternatively, the FAR ranges from 0% to 45% for all simulations, with the lowest FAR for the nonsupercell and deep-stable layer simulations, wherein the EIL depth (i.e., EIL top minus EIL base) is generally the smallest (Table 1).
Using the above definitions of POD and FAR, the critical success index (CSI) or threat score can be defined as (Gerapetritis and Pelissier 2004)
Thus, the CSI is maximized when the POD is large and the FAR is small, such that forecast skill is best with a CSI equal to 1. The distribution of the peak updraft core EIL CSI for the 95th percentile updraft core definition shows a median CSI for these simulations of 0.64. The maximum EIL CSI is for the NTD simulation (CSI = 0.78) and the minimum is for the NS simulation (CSI = 0.44). It should be noted that the CSI equally weights false alarms (higher FAR) and missed events (lower POD) such that the index alone gives an incomplete understanding of forecast skill. Following Roebber (2009), we plot the POD versus the success ratio (1 − FAR) in performance diagrams (Fig. 8). Thus, while forecast skill (i.e., CSI) increases to the upper right of the diagram, changes in CSI can be placed in the context of changes in both POD and FAR, commonly referred to as forecast bias. There is no clear trend as to which storm types or stable layer depths exhibit the best EIL forecast skill, but it does seem clear that the EIL generally has the poorest forecast skill for the nonsupercell simulations, consistent with their large entrainment (Fig. 6e) and low POD (Fig. 8d) with this updraft core definition.
b. Comparison of EIL performance for varied updraft core definitions
So far, our discussion of EIL performance has been limited to forecasts of the origins of air comprising the updraft core, where updraft core is defined based on the 95th percentile of vertical velocity exceeding 0.1 m s−1 occurring within a 30 km × 30 km subdomain centered on the right-moving supercell (or maximum vertical velocity for nonsupercell simulations) at each height. This definition was chosen in part because visual analysis of this threshold suggests it identifies the main updraft in each supercell simulation with little ancillary convection. Perhaps more relevant to severe weather forecasting, the average vertical profile of this threshold averaged over all simulations reaches a peak value of approximately 20 m s−1 (Fig. 7d), roughly corresponding to the terminal velocity of a severe (1 in.) hailstone (Markowski and Richardson 2010). Thus, for the supercell simulations here, the 95th percentile corresponds to updraft speeds capable of producing severe weather. However, other forecasting applications (e.g., lightning, flash flooding, severe winds) may require forecasts of updraft origins over a different definition of updraft core. Therefore, to answer our second research question, we examine EIL performance over several updraft core threshold definitions: the 50th, 75th, 90th, 95th, and 99th percentile of vertical velocity at each height.
In general, as the updraft core threshold is reduced the fraction of entrained air into the core updraft increases. For low updraft core thresholds (50th and 75th percentiles), the corresponding vertical velocity is quite weak (<5 m s−1) and about half of the air comprising the mid- to upper-tropospheric updraft core is entrained from above the EIL (Figs. 7a,b). Thresholds at and above the 90th percentile (Figs. 7c,d,e), exhibit a reduction in the entrained air above the EIL, though even for the most strict definition of updraft core (99th percentile), approximately 20% of the peak updraft air is entrained from above the EIL, on average. Similar analysis as that of Fig. 6 for other thresholds (not shown) suggests, however, that most of this entrainment occurs in the NS and WT simulations. Regardless of updraft core threshold, the amount of EIL air that does not participate in the updraft tends to increase with height. Again, this trend is likely driven by the relatively weak and shallow updrafts in the NS and WT simulations. In general, for simulations with relatively high percentages of “unused” EIL, air that does not participate in the peak updraft core tends to come from the upper few hundred meters of the EIL (not shown), where CAPE (CIN) is substantially decreased (increased) from its maximum (minimum) value, but still meets the T07 EIL criteria.
Returning to the EIL-as-forecast framework described in the previous section, we can compare the T07 EIL metric forecast skill for each updraft core threshold. Figure 8 shows performance diagrams where the FAR is equivalent to the percentage of EIL unused by the updraft core at the level of maximum vertical velocity in each simulation averaged over the 60–90 min simulation period (i.e., the “peak” updraft) and the POD is equivalent to the percentage of the peak updraft core composed of air from the T07 EIL for each threshold.5 In general, there is a transition from lower FAR (i.e., higher success ratio) and lower POD toward higher POD and higher FAR (i.e., lower success ratio) as the updraft core threshold is increased, suggesting a trade-off that limits overall forecast skill. In other words, the bias increases as the updraft core threshold becomes more stringent. Consistent with earlier analysis showing notably low entrainment, the WK simulations (hollow boxes in Fig. 8) tend to have the highest POD for all thresholds, but their FAR also tends to be highest, consistent with the greatest EIL depth. A notable exception to the trend of increasing FAR (i.e., decreasing success ratio) at greater updraft core thresholds are the NS simulations, which show a reduction in the percentage of the EIL unused by the updraft core (i.e., FAR) as the threshold increases (asterisks in Fig. 8). For the supercell simulations initialized in real environments (ST, WT, NT), there is no notable trend toward greater or less skill in any one storm type across all updraft thresholds. Finally, there is also no notable pattern between the depth of the surface inversion (i.e., red, blue, and black symbols) and either the POD or FAR.
Consistent with the aforementioned trade-off between POD and FAR, the median CSI is similar regardless of the updraft core definition (0.6–0.7; Figs. 8 and 9). The spread in CSI tends to decrease as the updraft core threshold increases, particularly above the 50th percentile definition. While the median CSI is slightly higher for the 90th percentile threshold, the minimum, maximum, and quartile values are all higher for the 95th percentile threshold. Apart from the generally poor performance for the NS, NSS, and NSD simulations, there is no clear trend toward greater CSI in any particular storm type (Fig. 8). As might be expected from the previous discussion of decreasing FAR in the NS simulations, improvements in CSI (particularly in the lowest quartile) are driven by increasing performance for the NS simulations. Thus, these results suggest that the T07 EIL metric exhibits robust skill in predicting the origins of peak updraft across a range of updraft core definitions, EIL depths, and EIL base heights, particularly for supercell storm types. It should be noted here that we have not evaluated the EIL performance for mesoscale convective systems (MCSs), where the large influence of the cold pool on updraft organization and propagation may alter our results.
c. Alternative EIL formulations
While the POD of the T07 EIL metric generally improved as the updraft core threshold increased, the concurrent increase in FAR results in an upper limit to the forecast skill of the standard T07 EIL metric. Thus, as potential improvements to the sounding-derived EIL forecast metric, we examined two ways to further constrain the EIL depth to reduce the FAR while maintaining the high POD of larger updraft core thresholds.
Our first new constraint was motivated by analysis of the distribution of passive tracer origin heights within the updraft core (Fig. 4), which suggest that this distribution is generally centered on or near the MU parcel (MUP) level. To incorporate this finding, we added an additional thermodynamic criterion to limit the EIL depth—requiring that, in addition to meeting the original T07 CAPE and CIN criteria, parcels originate within a fixed vertical distance of the MUP within a sounding. Figure 9b shows the distribution of CSI for EIL with the additional 1-km from MUP criterion, which resulted in the best skill of any distance tested (200 m, 500 m, 1 km). Compared to the standard T07 EIL metric (Fig. 9a), the spread in skill across updraft core thresholds increases. For lower updraft core thresholds (50th, 75th, 90th percentiles), the new EIL median CSI decreases substantially. At the most stringent updraft core threshold (99th percentile), the median CSI for the new EIL metric increases relative to the T07 EIL, and the interquartile range decreases substantially. The steady increase in skill as the updraft threshold increases indicates that this method is successful in limiting the additional FAR at larger updraft thresholds seen with the T07 EIL.
Our second attempt to further constrain the EIL metric was based on recent research that suggests updraft purity (i.e., the fraction of air from within the EIL) increases as the low-level, storm-relative (SR) flow increases (e.g., Peters et al. 2019). As such, we added a kinematic criterion to the T07 EIL metric—requiring that, in addition to meeting the original T07 CAPE and CIN criteria, parcels must also originate within the deepest contiguous layer with SR flow [defined based on the Bunkers et al. (2000) sounding-derived storm motion] exceeding a specified threshold within a sounding. We tested three SR-flow thresholds: 5, 10, and 15 m s−1. In cases where no contiguous layer of air met both the T07 thermodynamic criteria and our SR-flow threshold, we reverted to the standard T07 EIL metric (this was typical in the NS base state and occasional in the NT base state). The greatest skill came with the 10 m s−1 SR-flow threshold (Fig. 9c), perhaps reflective of earlier work showing this threshold is generally exceeded for supercell thunderstorms (Droegemeier et al. 1993). Like our MUP-based modification, the SR-flow criterion increases the range in forecast CSI and substantially decreases the median skill at the lower updraft core thresholds relative to the standard T07 EIL metric. Again, steadily increasing CSI as the updraft threshold is increased indicates that POD increases while holding FAR relatively constant. Finally, like the MUP-based modification, the best performance is at the 99th updraft core threshold, where the median is higher than the T07 EIL and the interquartile range is reduced.
Overall, both modifications suggest that EIL FAR can be reduced from the standard T07 EIL metric, though it often comes at the expense of decreased POD for some simulations. Indeed, simulations with the least entrainment, highest POD, and deepest EILs (the WK simulations) see considerable CSI improvements (CSI > 0.8 in some cases, not shown). The greatest potential improvements by constraining the EIL based on either proximity to the MUP or SR-flow come at the highest updraft thresholds, suggesting that any such modifications would primarily benefit forecasting supercellular convection with severe hazards. However, these are not substantial improvements from the T07 EIL definition, which has the added benefit of consistent CSI across a range of updraft definitions. Applying both the SR-flow and MUP criteria simultaneously does not result in significant changes in the EIL depth or skill for our sample of simulation base states, since the majority of layers within 1 km of the MUP level have SR flow exceeding 10 m s−1.
d. Dependence of EIL performance and updraft forcing on environmental parameters
Our final research question is concerned with the aspects of the near-storm environment that might dictate EIL performance and updraft strength overall. Again using the 95th percentile threshold for our definition of updraft core, we compared the EIL CSI for predicting the peak updraft origin as a function of a wide range of sounding-derived parameters (many of which are shown in Table 1) that might potentially influence the peak updraft strength and EIL. The majority of relationships have statistically insignificant correlations (using a p < 0.05 threshold for statistical significance), but all figures displaying a correlation coefficient have a statistically significant correlation.6 EIL CSI is moderately correlated with both MUCAPE (Fig. 10a) and 0–6-km bulk wind difference (Fig. 10b). These variables are also the strongest predictors of peak updraft strength (Fig. 11), suggesting along with earlier results, that EIL performance may increase with peak updraft strength. This is consistent with the result that stronger peak updrafts in the WK, ST, and NT simulations generally have the smallest fraction of entrained air from above the EIL (entrainment would reduce updraft speed via dilution of core buoyancy) and highest EIL POD (Fig. 8).
Peters et al. (2019) showed that low-level SR flow magnitude also increases peak updraft strength. Not surprisingly, average SR flow over the 0–3-km layer is also moderately correlated (r = 0.55) with EIL CSI in our simulations. There is also a strong correlation between 0–6-km shear and 0–3-km SR flow (Table 1; r = 0.95), consistent with Peters et al. (2019) finding that the link between deep-layer shear and supercell updraft intensity is primarily due to its association with stronger low-level SR flow and wider updrafts. Somewhat surprisingly there is not a statistically significant correlation between peak updraft speed and EIL CSI (not shown). This is because EIL CSI is also moderately correlated with EIL depth (Fig. 10c), and the simulations with the deepest EILs do not necessarily exhibit the strongest peak updrafts. Regardless of the connection to peak updraft strength, these results indicate that the T07 EIL metric is likely to exhibit higher skill in environments with stronger CAPE and deep-layer shear.
The results discussed above suggest that the link between EIL performance and the storm environment may be associated with factors that dictate updraft forcing and intensity. Therefore, to better understand the connection between the near-storm environment and low-level and peak updraft intensity, we compared sounding-derived parameters with the 1-km updraft and peak updraft strength (Fig. 11) and updraft forcing terms (Figs. 12–14). The 1-km updraft strength is strongly correlated with deep-layer shear (Fig. 11a), strongly correlated with 0–3 km and effective layer average SR flow (r = 0.81, r = 0.78, respectively), and moderately correlated with the ESRH (Fig. 11c). We also examined the relationship between 0–1 and 0–3 km fixed-layer SRH with the 1-km updraft, but found no statistically significant correlation with either—thus highlighting the utility of using SRH integrated over the EIL, particularly in environments with surface-based stable layers. Interestingly, the 1-km updraft strength has a positive, statistically significant correlation with MUCAPE (not shown) and MLCAPE (Fig. 11e). This is surprising because updraft parcels are typically negatively buoyant at 1 km (Fig. 12), which is below the LFC for most simulations. Indeed there is a negative (though not statistically significant) relationship between the 1-km updraft strength and the buoyancy averaged over the 1-km core updraft (Fig. 13a). There are two possible interpretations of this result. One is that CAPE and shear may be inadvertently related in our small sample size. This is possible considering the simulations with the highest MUCAPE (NT, NTS, NTD) have significantly larger MUCAPE than the simulations with the weakest deep-layer shear (NS, NSS, NSD) as shown in Table 1. Indeed, this may be the most likely explanation as there is a weak positive, albeit insignificant correlation between MUCAPE and deep-layer shear in the 15 base states (not shown).
If the potential relationship between CAPE and shear in our sample size does not fully explain the positive correlation between MUCAPE and 1-km updraft strength, it is possible that there is a physical explanation. One such explanation is that the BPA, which tends to be positive at low levels, especially in the WK and ST simulations (Figs. 12a,b), increases at low levels as MUCAPE increases. However, there is only a very weak, statistically insignificant correlation between the 1-km updraft core BPA and CAPE (not shown) and a weak statistically insignificant correlation between 1-km updraft strength and BPA averaged over the 1-km updraft core (not shown), limiting the likelihood of this scenario. In any case, there is a moderate, statistically significant correlation (r = 0.64) between the percentage of air from below the EIL that participates in the 1-km updraft core and MLCIN, which may be partly explained by the positive (but not statistically significant) correlation between MLCIN and BPA averaged over the 1-km updraft core (not shown). Thus the potential role of BPA in low-level updraft forcing remains uncertain, and should be explored in more detail in future studies with a larger sample of cases.
Regardless of the potential role of BPA, it seems clear that dynamic updraft forcing (DPA) plays a large role in the low-level updraft forcing. Every storm type displays positive nonlinear DPA in the lowest few hundred meters of the updraft core (Fig. 12), potentially caused by either a dynamic pressure maximum at the surface associated with convergence along the gust front and/or dynamic pressure minimum aloft driven by elevated rotation. Indeed, nonlinear DPA tends to be the largest updraft forcing over the lowest 2 km of each supercell storm type, consistent with findings of Weisman and Rotunno (2000) for idealized hodograph shapes. Interestingly, the two tornadic storm types (ST, WT; it is not clear if the WK base-state would produce a tornado in nature) display slightly deeper layers of enhanced nonlinear DPA (Figs. 12b,c). The linear DPA term is close to zero at all levels of each simulation updraft core. This is not surprising, as linear dynamic pressure perturbations exist somewhat symmetrically around the updraft core in areas of strong horizontal gradients in vertical velocity, thus their contribution to the DPA averages out over the updraft maximum. The 1-km DPA shown to be largely responsible for low-level updrafts in these simulations is strongly correlated with environmental ESRH (Fig. 14a) and moderately correlated with deep-layer shear (Fig. 14b), likely explaining the relationship between these parameters and low-level updraft strength. DPA at 1-km is also positively correlated with 0–1 km SRH (r = 0.56; not shown) and 0–3 km SRH (r = 0.52; not shown), but these are notably weaker relationships with than those with ESRH (r = 0.75; Fig. 14a), perhaps explaining the lack of a significant correlation between fixed-layer SRH values and low-level updraft strength in these simulations. The overall relationship between SRH (regardless of the layer of integration) and low-level DPA is consistent with existing studies of supercell low-level updraft forcing in more realistic environments (e.g., Coffer et al. 2017; Parker 2017).
The peak updraft strength in these simulations is strongly correlated with deep-layer shear (Fig. 11b) and moderately correlated with MUCAPE (Fig. 11f), but not significantly correlated with ESRH (Fig. 11d). The strong correlation with deep-layer shear might suggest a prominent role for DPA in the peak updraft forcing; however, recent work has suggested this relationship is largely a function of stronger low-level SR flow leading to wider updrafts with cores that are protected from entrainment-driven dilution and ultimately are able to realize a greater amount of buoyancy (Peters et al. 2019). In our simulations, DPA is considerable throughout the lowest half of the troposphere in the WK, ST, and NT simulations (Fig. 12), but may be offset by negative BPA as suggested by Peters et al. (2019) in at least the WK and NT simulations (Figs. 12a,d) and downward DPA above the midlevel rotation maximum. Buoyancy tends to dominate the updraft forcing within the updraft core above 2 km for all of the supercell types (Figs. 12a–d). Consistent with the structural differences of the NS updraft (Fig. 5), the NS updraft core displays pockets of weak buoyancy at varying altitudes characteristic of a chain of rising thermals (Fig. 12e). Ultimately, the peak updraft strength in these simulations is strongly correlated with both the maximum updraft core buoyancy and DPA (Figs. 13b,d) at any height below the peak updraft, but values of maximum DPA are often much less than those of buoyancy forcing.
As might be expected, the maximum updraft core buoyancy is positively correlated with both MLCAPE and deep-layer shear (Figs. 14e,f), with the latter relationship potentially a consequence of wider updrafts in high-shear environments consistent with a moderate correlation between maximum updraft core buoyancy and both 0–3 km and effective layer SR flow (r = 0.59, r = 0.63, respectively). There are positive (though not statistically significant) relationships between deep-layer shear and updraft core area, updraft area and percentage of updraft originating in the EIL, and updraft area and buoyancy (not shown) that are consistent with the results of Peters et al. (2019). Maximum DPA is also well-correlated with deep-layer shear as might be expected (Fig. 14d), but it is also strongly correlated with MLCAPE (Fig. 14c), which is somewhat unexpected. A potential explanation for this relationship is the dependence of DPA on vertical velocity—linear dynamic pressure perturbations are directly dependent on (gradients in) vertical velocity, and nonlinear perturbations may be indirectly dependent on vertical velocity through the stretching term in the vertical vorticity equation. Thus, greater updraft buoyancy forcing resulting from more CAPE and/or wider updrafts in strongly sheared environments may also increase DPA, as suggested by Peters et al. (2020a).
4. Summary and conclusions
In this paper, we seek to verify the performance of the Thompson et al. (2007) formulation for the effective inflow layer (EIL) of deep convective storms using a suite of idealized numerical simulations. Analytic soundings and near-storm environments obtained from model analyses of real events were chosen for a range of supercell storm types (significantly tornadic, weakly tornadic, nontornadic) and a nonsupercell severe storm. To modify the EIL base, surface-based temperature inversions of varying depth were added to each base state, resulting in 15 simulations with varying EIL depths, MUP heights, low-level and deep-layer shear, and CAPE. Using an array of passive tracers initialized at varying vertical levels in each simulation, the proportion of air within the resulting updraft cores originating below, within, and above the forecast EIL were determined. Considering the EIL as a forecast of updraft parcel origin, the skill of the EIL was determined for varying updraft core definitions and as a function of near-storm environments. Finally, the role of the near-storm environment on low-level and peak updraft forcing was discussed.
The primary conclusions of our study, in regards to the specific science questions we sought to answer are as follows:
How accurate is the EIL as formulated by T07 in predicting the true inflow layer (i.e., the origin height of air participating in the “updraft core”) of supercell and nonsupercell updrafts in a variety of realistic environments?
The EIL tends to perform reasonably well for most storm types. The worst performance was, as hypothesized, for nonsupercells that are characterized by weaker, intermittent, thermal-like updrafts, with significant entrainment of air from above the EIL.
How does the EIL performance vary as the definition of “updraft core” is changed?
As the vertical velocity threshold used to define the updraft core is increased, the probability of the EIL detecting the origin of peak updraft parcels increases as hypothesized, but the chance of the EIL falsely identifying updraft parcels also increases such that overall EIL forecast skill is steady across a range of updraft core definitions.
Alternative definitions of the EIL wherein the depth of the EIL is constrained by proximity to the MUP level or SR flow may increase performance for the most stringent definitions of updraft core, but generally result in worse performance for weaker updraft core thresholds.
Which environmental parameters best predict EIL performance, and more broadly, low-level and peak updraft forcing and intensity in real environments?
As hypothesized, EIL performance improves as CAPE, deep-layer shear, and EIL depth increase.
Low-level updraft strength is most strongly modulated by dynamic forcing, which is strongly correlated with ESRH and, to a lesser extent, deep-layer shear.
Peak updraft strength is strongly modulated by both buoyancy and dynamic forcing. Maximum updraft core buoyancy and dynamic updraft forcing are strongly related to both CAPE and deep-layer shear, but not ESRH.
The importance of DPA in low-level updraft forcing is largely consistent with prior studies (e.g., Weisman and Klemp 1984; Weisman and Rotunno 2000; Markowski and Richardson 2014; Coffer et al. 2017; Parker 2017). However, there is some indication that larger CAPE may be associated with increased upward buoyancy pressure accelerations below the LFC, thereby also increasing the low-level updraft despite its negative buoyancy. This finding awaits further testing over a larger sample of simulations, including in rotating updrafts embedded in MCSs, where Nielsen and Schumacher (2018) found similar updraft forcing distributions. The role of deep-layer shear in intensifying maximum buoyancy and thereby controlling the peak updraft intensity remains uncertain, though our results are consistent with earlier work by Peters et al. (2019) suggesting that supercell updraft cores are comprised of a much higher fraction of EIL air and less entrainment of air from above the EIL. Moreover, the physical explanation behind the finding that dynamic forcing is related to environmental CAPE is uncertain, but we hypothesize that this is because both nonlinear and linear pressure perturbations are directly or indirectly dependent on updraft strength, such that DPA may increase if updrafts are intensified through buoyancy. We note that these results apply to the mature, relatively steady phase of simulated supercells, but the roles of entrainment and DPA may vary at earlier, developing stages of supercells.
To our knowledge, this is the first study conducting a comprehensive analysis of simulated supercell updraft core origins at all levels in a variety of environments. The results of this paper confirm that the commonly used EIL metric as formulated by Thompson et al. (2007) is largely an accurate reflection of the air ingested by isolated deep convective updrafts over a relatively wide range of realistic thermodynamic and kinematic environments, particularly those associated with elevated supercells. This finding is somewhat at odds with those of Nowotarski et al. (2011) who showed that supercell updrafts could ingest air from near-ground stable layers; however, their study employed a less stringent definition of “updraft” than those tested here. Our results suggest (particularly at weaker updraft core thresholds) that the vast majority of non-EIL air participating in updraft cores aloft is via entrainment, not stable air from below the EIL. Entrainment is particularly notable in nonsupercellular convection. We also note that the commonly used Weisman and Klemp (1982) thermodynamic and kinematic profile results in supercell updrafts with unrealistically small entrainment (at least compared with the real supercell environments simulated herein).
Future work should explore how the EIL performs for multicellular/MCS storm modes, where updraft organization is strongly influenced by outflow and previous studies have shown a layer-lifted approach that integrates outflow interactions with low-level shear over the evolution of individual storms may be more appropriate (Alfaro and Coniglio 2018). Indeed, such an approach may improve inflow approximations of supercellular convection, but testing such methods is beyond the scope of this study. Finally, despite the apparent success of the T07 EIL parameter, our results suggest that the strong dependence of supercell updraft properties on both CAPE and shear warrant continued exploration and efforts to incorporate kinematic variables into parameters that forecast supercell updraft characteristics.
We are indebted to George Bryan for his development and continued support of CM1. We thank Rich Thompson and Roger Edwards for providing their proximity sounding dataset used to initialize some of the simulations in this study. J. Peters’s and C. Nowotarski’s efforts were supported by the National Science Foundation (NSF) Grants AGS-1928666 and AGS-1928319, respectively. Additionally, J. Peters and J. Mulholland were partially supported by NSF Grants AGS-1841674 and Department of Energy Atmospheric System Research Grant DE-SC0000246356. Simulations were performed on Texas A&M High Performance Research Computing resources. Many sounding-derived parameters and sounding figures were calculated/plotted using SHARPpy (Blumberg et al. 2017). We also thank the editor and reviewers whose constructive comments helped to clarify the manuscript.
Data availability statement. All simulation source code, configuration files, and analysis code are archived locally and available upon request to the corresponding author.
Most sounding-derived parameters were computed using the Sounding/Hodograph Analysis and Research Program in Python [SHARPpy; Blumberg et al. (2017)], which uses the modified Bunkers right-moving storm motion (Bunkers et al. 2014). This storm motion uses a pressure-weighted mean wind with the effective inflow layer as a base. Thus, some fixed-layer, storm-relative quantities vary slightly for simulations with identical hodographs, but varied near-surface stability.
The results at any particular time during the mature phase of deep convection are qualitatively similar to the 60–90-min average results.
The features of the WKD updraft and airflow shown here are generally representative of most of the supercell simulations.
Herein we consider any passive tracer level with a concentration of less than 2% as “not meaningfully participating” in the updraft core at any height. These parcels may still be lifted or even become positively buoyant, but we consider them “unused” in this evaluation if they do not participate in the updraft core, as defined by a given vertical velocity percentile threshold.
By convention, a performance diagram plots the success ratio on the x axis, but the success ratio = 1 − FAR, such that an increasing success ratio represents a decreasing FAR.
Though it is often difficult to ascribe qualitative importance to correlation values beyond their statistical significance, we classify r < 0.50 as a weak correlation, 0.50 < r < 0.75 as moderate, and r > 0.75 as strong.