Search NASA⌕ Search

Engineering topics

Peters, John M.

Publications and source records attributed to Peters, John M..

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Cumulonimbus Clouds Convert a Smaller Fraction of CAPE into Kinetic Energy in a Warmer Atmosphere

Abstract This study investigates how entrainment’s diluting effect on cumulonimbus updraft buoyancy is affected by the temperature of the troposphere, which is expected to increase by the end of the century. A parcel model framework is constructed that allows for independent variations in the temperature ( T ), the entrainment rate ε , the free-tropospheric relative humidity (RH), and the convective available potential energy (CAPE). Using this framework, dilution of buoyancy is evaluated with T and RH independently varied and with CAPE either held constant or increased with temperature. When CAPE is held constant, buoyancy decreases as T increases, with parcels in warmer environments realizing substantially smaller fractions of their CAPE as kinetic energy (KE). This occurs because the increased moisture difference between an updraft and its surroundings at warmer temperatures drives greater updraft dilution. Similar results are found in midlatitude and tropical conditions when CAPE is increased with temperature. With the expected 6%–7% increase in CAPE per kelvin of warming, KE only increases at 2%–4% K −1 in narrow updrafts but tracks more closely with CAPE at 4%–6% in wider updrafts. Interestingly, the rate of increase in the KE with T becomes larger than that of CAPE when the later quantity increases at more than 10% K −1 . These findings emphasize the importance of considering entrainment in studies of moist convection’s response to climate change, as the entrainment-driven dilution of buoyancy may partially counteract the influence of increases in CAPE on updraft intensity. Significance Statement Cumulonimbus clouds mix air with their surrounding environment through a process called entrainment, which controls how efficiently environmental energy is converted into upward speed in thunderstorm updrafts. Our research shows that warmer temperatures will exacerbate the moisture difference between cumulonimbus updrafts and their surroundings, leading to greater mixing and less efficient conversion of environmental energy into updraft speeds. This effect should be considered in future research that investigates how climate change will affect cumulonimbus clouds.

Peters, John M.↗
DOE OSTI2024

Should We Conserve Entropy or Energy when Computing CAPE with Mixed-Phase Precipitation Physics?

Abstract The rapidly increasing resolution of global atmospheric reanalysis and climate model datasets necessitates finding methods for computing convective available potential energy (CAPE) both efficiently and accurately. To this end, this article compares two common methods for computing CAPE which conserve either energy or entropy. Inaccuracies in these computations arise from both physical and numerical errors. For instance, computing CAPE with entropy conserved results in physical errors from nonequilibrium phase transitions but minimizes numerical errors because solutions are analytic at each height. In contrast, computing CAPE with energy conserved avoids these physical errors, but accumulates numerical errors that are grid-resolution-dependent because the numerical integration of a differential equation is required. Analysis of CAPE computed with large databases of soundings from the tropical Amazon and midlatitude storm environments shows that physical errors from the entropy method are typically 1%–3% as large as CAPE, which is comparable to the numerical errors from conserving energy with grid spacing of 25 and 250 m using explicit first-order and second-order integration schemes, respectively. Errors in entropy-based CAPE calculations are also insensitive to vertical grid spacing, in contrast to energy-based calculations whose error strongly scales with the grid spacing. It is shown that entropy-based methods are advantageous when intercomparing datasets with differing vertical resolution because they produce accurate and reasonably fast results that are insensitive to grid resolution, whereas a second-order energy-based method is advantageous when analyzing data with a consistent vertical resolution because of its superior computational efficiency. Significance Statement Convective available potential energy (CAPE) is a measure of instability in the atmosphere that helps forecasters and researchers understand when and where thunderstorms will form. The purpose of this article is to identify the most efficient and accurate methods for computing CAPE. Two methods are considered here, one that relates to the entropy (a measure of thermodynamic disorder) of an air parcel and one that relates to the energy of an air parcel. Results indicate that the entropy method is most accurate and insensitive to the resolution of the data used for the calculation (which can vary considerably), whereas the energy method uses the least computation time.

Peters, John M.↗
DOE OSTI2024

How Does Vertical Wind Shear Influence Updraft Characteristics and Hydrometeor Distributions in Supercell Thunderstorms?

Abstract Vertical wind shear is known to affect supercell thunderstorms by displacing updraft hydrometeor mass downshear, thereby facilitating the storms’ longevity. Shear also impacts the size of supercell updrafts, with stronger shear leading to wider, less dilute, and stronger updrafts with likely greater hydrometeor production. To more clearly define the role of shear across different vertical layers on hydrometeor concentrations and displacements relative to supercell updrafts, a suite of idealized numerical model simulations of supercells was conducted. Shear magnitudes were systematically varied across the 0–1, 1–6, and 6–12 km AGL layers, while the thermodynamic environment was held fixed. Simulations show that as shear magnitude increases, especially from 1 to 6 km, updrafts become wider and less dilute with an increase in hydrometeor loading, along with an increase in the low-level precipitation area/rate and total precipitation accumulation. Even with greater updraft hydrometeor loading amid stronger shear, updrafts are more intense in stronger shear simulations due to larger thermal buoyancy owing to wider, less dilute updraft cores. Furthermore, downshear hydrometeor displacements are larger in environments with stronger 1–6-km shear. In contrast, there is relatively less sensitivity of hydrometeor concentrations and displacements to variations in either 0–1- or 6–12-km shear. Results are consistent across free tropospheric relative humidity sensitivity simulations, which show an increase in updraft size and hydrometeor mass with increasing free tropospheric relative humidity owing to a reduction in entrainment-driven dilution for wider updrafts in moister environments. Significance Statement Rotating thunderstorms, known as supercells, are able to persist for multiple hours. One common explanation is that large changes in wind speed and/or direction with height, or shear, transport rain/hail away from supercell updrafts, supporting their maintenance. The strong shear within supercell environments, however, may also lead to greater rail/hail amounts, thereby leading to weaker storms due to this extra mass of water/ice within updrafts. Furthermore, the impact of shear across different height layers on supercell rain/hail characteristics has not been thoroughly investigated. In this study, computer simulations of supercells were conducted to determine that shear occurring between 1 and 6 km above ground level has a large impact on rain/hail distribution in supercells and that stronger shear in this layer leads to wider/stronger supercells with greater rain/hail accumulations at the surface. Additionally, some of the extra mass of water/ice is transported farther away from updrafts due to the stronger environmental storm-relative winds.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Simplified Approximations of Direct Cumulus Entrainment and Detrainment

Abstract In recent years, direct calculations of simulated cumulus entrainment and detrainment have facilitated new physical insights into these highly elusive but critically important processes. However, these calculations require substantial computational resources that may limit their widespread usage. To facilitate such calculations, two simplified approximations of direct cumulus entrainment and detrainment are examined herein. The first approximation, termed the “semidirect” method, follows a standard bulk approach but makes more realistic assumptions about the sources of entrained and detrained air near the cloud edges. In contrast, the second approximation (the “projection” method) uses the governing equations of motion to project whether grid points near the cloud edge will entrain or detrain as the mean cloud ascends by one grid point. Verification exercises using large-eddy simulations reveal that both methods generally agree better with corresponding direct entrainment/detrainment estimates than the traditional bulk formulation, with the projection method outperforming the semidirect method. The two methods can be used in a synergistic fashion, with the semidirect method helping to optimize the projection method, to suit a wide range of applications. Because the latter incorporates the essential dynamics of entrainment and detrainment at the local scale, it can be used to gain physical insight into the causal mechanisms regulating these complex processes.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Modulation of Tropical Convection-Circulation Interaction by Aerosol Indirect Effects in Convective Self-Aggregation Simulations of a Gray Zone Global Model

Disentangling the response of tropical convective updrafts to enhanced aerosol concentrations has been challenging. Leading theories for explaining the influence of aerosol concentrations on tropical convection are based on the dynamical response of convection to changes in cloud microphysics, neglecting possible changes in the environment. In recent years, global convection-permitting models (GCPM) have been developed to circumvent problems arising from imposing artificial scale separation on physical processes associated with deep convection. Here, we use a global model in the convective gray zone that partially simulates deep convection to investigate how enhanced concentrations of aerosols that act as cloud condensate nuclei (CCN) impact tropical convection features by modulating the convection-circulation interaction. Results from a pair of idealized non-rotating radiative-convective equilibrium simulations show that the enhanced CCN concentration leads to weaker large-scale circulation, the closeness of deep convective systems to the moist cluster edges, and more mid-level cloud water at an equilibrium state in which convective self-aggregation occurred. Correspondingly, the enhanced CCN concentration modulates how the physical processes that support or oppose convective aggregation maintain the aggregated state at equilibrium. Overall, the enhanced CCN concentration facilitates the development of deep convection in a drier environment but reduces mean precipitation. Our results emphasize the importance of allowing atmospheric phenomena to evolve continuously across spatial and temporal scales in simulations when investigating the response of tropical convection to changes in cloud microphysics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

What Controls the Entrainment Rate of Dry Buoyant Thermals with Varying Initial Aspect Ratio?

Here, this study uses theory and numerical simulations to analyze the nondimensional spreading rate a (change in radius with height) of buoyant thermals as they rise and entrain surrounding environmental fluid. A focus is on how a varies with initial thermal aspect ratio A r , defined as height divided by width of the initial buoyancy perturbation. An analytic equation for thermal ascent rate w t that depends on a is derived from the thermal-volume-averaged momentum bud get equation. The thermal top height when w t is maximum, defining a critical height z c , is inversely proportional to a. The height z c also corresponds to the thermal top height when buoyant fluid along the thermal’s vertical axis is fully replaced by entrained nonbuoyant environmental fluid rising from below the thermal. The time scale for this process is controlled by the vertical velocity of parcels rising upward through the thermal’s core. This parcel vertical velocity is approximated from Hill’s analytic spherical vortex, yielding an analytic inverse relation between a and A r . Physically, this a–A r relation is connected to changes in circulation as A r is modified. Numerical simulations of thermals with A r varied from 0.5 to 2 give a values close to the analytic theoretical relation, with a factor of ~3 decrease in a as A r is increased from 0.5 to 2. The theory also explains why a of initially spherical thermals from past laboratory and modeling studies is about 0.15. Overall, this study provides a theoretical underpinning for understanding the entrainment behavior of thermals, relevant to buoyantly driven atmospheric flows.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

An Analytic Formula for Entraining CAPE in Midlatitude Storm Environments

Abstract This article introduces an analytic formula for entraining convective available potential energy (ECAPE) with an entrainment rate that is determined directly from an environmental sounding, rather than prescribed by the formula user. Entrainment is connected to the background environment using an eddy diffusivity approximation for lateral mixing, updraft geometry assumptions, and mass continuity. These approximations result in a direct correspondence between the storm-relative flow and the updraft radius and an inverse scaling between the updraft radius squared and entrainment rate. The aforementioned concepts, combined with the assumption of adiabatic conservation of moist static energy, yield an explicit analytic equation for ECAPE that depends entirely on state variables in an atmospheric profile and a few constant parameters with values that are established in past literature. Using a simplified Bernoulli-like equation, the ECAPE formula is modified to account for updraft enhancement via kinetic energy extracted from the cloud’s background environment. CAPE and ECAPE can be viewed as predictors of the maximum vertical velocity w max in an updraft. Hence, these formulas are evaluated using w max from past numerical modeling studies. Both of the new formulas improve predictions of w max substantially over commonly used diagnostic parameters, including undiluted CAPE and ECAPE with a constant prescribed entrainment rate. The formula that incorporates environmental kinetic energy contribution to the updraft correctly predicts instances of exceedance of by w max , and provides a conceptual explanation for why such exceedance is rare among past simulations. These formulas are potentially useful in nowcasting and forecasting thunderstorms and as thunderstorm proxies in climate change studies. Significance Statement Substantial mixing occurs between the upward-moving air currents in thunderstorms (updrafts) and the surrounding comparatively dry environmental air, through a process called entrainment. Entrainment controls thunderstorm intensity via its diluting effect on the buoyancy of air within updrafts. A challenge to representing entrainment in forecasting and predictions of the intensity of updrafts in future climates is to determine how much entrainment will occur in a given thunderstorm environment without a computationally expensive high-resolution simulation. To address this gap, this article derives a new formula that computes entrainment from the properties of a single environmental profile. This formula is shown to predict updraft vertical velocity more accurately than past diagnostics, and can be used in forecasting and climate prediction to improve predictions of thunderstorm behavior and impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Supercell Low-Level Mesocyclones: Origins of Inflow and Vorticity

Abstract The development and intensification of low-level mesocyclones in supercell thunderstorms have often been attributed, at least in part, to augmented streamwise vorticity generated baroclinically in the forward flank of supercells. However, the ambient streamwise vorticity of the environment (often quantified via storm-relative helicity), especially near the ground, is particularly skillful at discriminating between nontornadic and tornadic supercells. This study investigates whether the origins of the inflow air into supercell low-level mesocyclones, both horizontally and vertically, can help explain the dynamical role of environmental versus storm-generated vorticity in the development of low-level mesocyclone rotation. Simulations of supercells, initialized with wind profiles common to supercell environments observed in nature, show that the air bound for the low-level mesocyclone primarily originates from the ambient environment (rather than from along the forward flank) and from very close to the ground, often in the lowest 200–400 m of the atmosphere. Given that the near-ground environmental air comprises the bulk of the inflow into low-level mesocyclones, this likely explains the forecast skill of environmental streamwise vorticity in the lowest few hundred meters of the atmosphere. The low-level mesocyclone does not appear to require much augmentation from the development of additional horizontal vorticity in the forward flank. Instead, the dominant contributor to vertical vorticity within the low-level mesocyclone is from the environmental horizontal vorticity. This study provides further context to the ongoing discussion regarding the development of rotation within supercell low-level mesocyclones. Significance Statement Supercell thunderstorms produce the majority of tornadoes, and a defining characteristic of supercells is their rotating updraft, known as the “mesocyclone.” When the mesocyclone is stronger at lower altitudes, the likelihood of tornadoes increases. The purpose of this study is to understand if the rotation of the mesocyclone in supercells is due to horizontal spin present in the ambient environment or whether additional horizontal spin generated by the storm itself primarily drives this rotation. Our results suggest that inflow air into supercells and low-level mesocyclone rotation are mainly due to the properties of the environmental inflow air, especially near the ground. This hopefully provides further context to how our community views the development of low-level mesocyclones in supercells.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Entrainment Makes Pollution More Likely to Weaken Deep Convective Updrafts Than Invigorate Them

Are the results of aerosol invigoration studies that neglect entrainment valid for diluted deep convective clouds? We address this question by applying an entraining parcel model to soundings from tropical and midlatitude convective environments, wherein pollution is assumed to increase parcel condensate retention. Invigoration of 5%–10% and <2% is possible in undiluted tropical and midlatitude parcels respectively when freezing is rapid. This occurs because the positive buoyancy contribution from freezing is larger than the negative buoyancy contribution from condensate loading, leading to positive net condensate contribution to buoyancy. However, aerosol-induced weakening is more likely when realistic entrainment rates occur because water losses from entrainment more substantially reduce the latent heating relative to the loading contribution. This leads to larger net negative buoyancy contribution from condensates in polluted than in clean entraining parcels. Our results demonstrate that accounting for entrainment is critical in conceptual models of aerosol indirect effects in deep convection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Seasonal controls on isolated convective storm drafts, precipitation intensity, and life cycle as observed during GoAmazon2014/5

Abstract. Isolated deep convective cloud life cycle and seasonal changes in storm properties are observed for daytime events during the US Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Green Ocean Amazon Experiment (GoAmazon2014/5) campaign to understand controls on storm behavior. Storm life cycles are documented using surveillance radar from initiation through maturity and dissipation. Vertical air velocity estimates are obtained from radar wind profiler overpasses, with the storm environment informed by radiosondes. Dry-season storm conditions favored reduced morning shallow cloud coverage and larger low-level convective available potential energy (CAPE) than wet-season counterparts. The typical dry-season storm reached its peak intensity and size earlier in its life cycle compared with wet-season cells. These cells exhibited updrafts in core precipitation regions (Z>35 dBZ) to above the melting level as well as persistent downdrafts aloft within precipitation adjacent to their cores. Moreover, dry-season cells recorded more intense updrafts to earlier life cycle stages as well as a higher incidence of strong updrafts (i.e., >5 m s−1) at low levels. In contrast, wet-season storms were longer-lived and featured a higher incidence of moderate (i.e., 2–5 m s−1) updrafts aloft. These storms also favored a shift in their most intense properties to later life cycle stages. Strong downdrafts were less frequent within wet-season cells aloft, indicating a potential systematic difference in draft behaviors, as linked to graupel loading and other factors between the seasons. Results from a stochastic parcel model suggest that dry-season cells may expect stronger updrafts at low levels because of larger low-level CAPE in the dry season. Wet-season cells anticipate strong updrafts aloft because of larger free-tropospheric relative humidity and reduced entrainment-driven dilution. Enhanced dry-season downdrafts are partially attributed to increased evaporation, dry-air entrainment mixing, and negative buoyancy in regions adjacent to sampled dry-season cores.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Near-Cloud Atmospheric Ingredients for Deep Convection Initiation

A lack of routine environmental observations located near deepening cumulus congestus clouds limits verification of important theorized and simulated updraft–environment interaction processes occurring during deep convection initiation (CI). We analyze radiosonde profiles collected during several hundred CI events near a mountain range in central Argentina during the CACTI field campaign. Statistical analyses illustrate environmental conditions supporting radar-observed CI outcomes that span a spectrum of convective cell depths, widths, and durations, as well as events lacking precipitating convection. Tested environmental factors include a large variety of sounding-derived measurements of CAPE, CIN, moisture, terrain-relative winds, vertical shear, and lifted parcel properties, with supplemental model reanalysis of background larger-scale vertical motion. CAPE and CIN metrics do not consistently differentiate CI success from failure. Only a few environmental factors contain consistent monotonic relationships among the spectrum of cloud depths achieved during CI: (i) the depth and strength of background ascent, and (ii) the component of low-level flow oriented parallel to the ridgeline. These metrics suggest that the ability of the surrounding flow to lift parcels to their LFC and terrain-modified flow are consistently relevant processes for CI. Low- to midlevel relative humidity strongly discriminated between CI and non-CI events, likely reflecting entrainment-driven dilution processes. However, we could not confidently conclude that relative humidity similarly discriminated robust from marginal CI events. Circumstantial evidence was found linking cell width, an important cloud property governing the probability of CI, to LCL height, boundary layer depth, depth and magnitude of the CIN layer, and ambient wind shear.

Storm environments↗