Search NASA⌕ Search

SEARCH · Search NASA

Results for “parcels”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

Fallowed agricultural lands dominate anthropogenic dust sources in California

Air pollution remains a major problem in many parts of California, significantly impacting public health and regional climate. However, the contribution of anthropogenic dust from agricultural sources, among major pollutants in California’s semi-arid Central Valley, remains largely unclear. Here, we used the Cropland Data Layer from the U.S. Department of Agriculture to identify crop types and land use/cover and leveraged satellite-derived estimates of major dust events between 2008 and 2022 over California. We identified fallowed land—an unplanted agricultural land parcel—as a key anthropogenic dust source in California. Specifically, we find that the Central Valley accounts for about 77% of total fallowed land areas in California, where they are associated with about 88% of major anthropogenic dust events. We also find that the geographic coverage of these fallowed lands expanded between 2008 and 2022 with associated increasing anthropogenic dust activities. Additionally, these anthropogenic dust activities are sensitive to the drought severity over the fallowed lands, with potential cumulative effects on downstream dust burden during prolonged multi-year drought conditions. Overall, our results have important implications for public health, including increased risk for Valley fever and for regional climates, such as increases in extreme precipitation and snowmelt over the Sierra Nevada.

Climate sciences↗

The role of turbulent fluctuations in aerosol activation and cloud formation

Aerosol indirect effects are one of the leading contributors to cloud radiative properties relevant to climate. Aerosol particles become cloud droplets when the ambient relative humidity (saturation ratio) exceeds a critical value, which depends on the particle size and chemical composition. In the traditional formulation of this problem, only average, uniform saturation ratios are considered. Using experiments and theory, we examine the effects of fluctuations, produced by turbulence. Our measurements, from a multiphase, turbulent cloud chamber, show a clear transition from a regime in which the mean saturation ratio dominates to one in which the fluctuations determine cloud properties. The laboratory measurements demonstrate cloud formation in mean-subsaturated conditions (i.e., relative humidity <100%) in the fluctuation-dominant activation regime. The theoretical framework developed to interpret these measurements predicts a transition from a mean- to a fluctuation-dominated regime, based on the relative values of the mean and standard deviation of the environmental saturation ratio and the critical saturation ratio at which aerosol particles activate or become droplets. The theory is similar to the concept of stochastic condensation and can be used in the context of the atmosphere to explore the conditions under which droplet activation is driven by fluctuations as opposed to mean supersaturation. Overall, it provides a basis for future development of cloud droplet activation parameterizations that go beyond the internally homogeneous parcel calculations that have been used in the past.

54 ENVIRONMENTAL SCIENCES↗

Transport and chemistry of isoprene and its oxidation products in deep convective clouds

Deep convective clouds can transport trace gases from the planetary boundary layer into the upper troposphere where subsequent chemistry may impact aerosol particle formation and growth. In this modelling study, we investigate processes that affect isoprene and its oxidation products injected into the upper troposphere by an isolated deep convective cloud in the Amazon. We run a photochemical box model with coupled cloud microphysics along hundreds of individual air parcel trajectories sampled from a cloud-resolving model simulation of a convective event. The box model simulates gas-phase chemical reactions, gas scavenging by liquid and ice hydrometeors, and turbulent dilution inside a deep convective cloud. The results illustrate the potential importance of gas uptake to anvil ice in regulating the intensity of the isoprene oxidation and associated low volatility organic vapour concentrations in the outflow. Isoprene transport and fate also depends on the abundance of lightning-generated nitrogen oxide radicals (NO x = NO + NO 2 ). If gas uptake on ice is efficient and lightning activity is low, around 30% of the boundary layer isoprene will survive to the cloud outflow after approximately one hour of transport, while all the low volatile oxidation products will be scavenged by the cloud hydrometeors. If lightning NO x is abundant and gas uptake by ice is inefficient, then all isoprene will be oxidised during transport or in the immediate outflow region, while several low volatility isoprene oxidation products will have elevated concentrations in the cloud outflow. Reducing uncertainties associated with the uptake of vapours on ice hydrometeors, especially HO 2 and oxygenated organics, is essential to improve predictions of isoprene and its oxidation products in deep convective outflows and their potential contribution to new particle formation and growth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

AutoAtlas: Neural Network for 3D Unsupervised Partitioning and Representation Learning

Here we present a novel neural network architecture called AutoAtlas for fully unsupervised partitioning and representation learning of 3D brain Magnetic Resonance Imaging (MRI) volumes. AutoAtlas consists of two neural network components: one neural network to perform multi-label partitioning based on local texture in the volume, and a second neural network to compress the information contained within each partition. We train both of these components simultaneously by optimizing a loss function that is designed to promote accurate reconstruction of each partition, while encouraging spatially smooth and contiguous partitioning, and discouraging relatively small partitions. We show that the partitions adapt to the subject specific structural variations of brain tissue while consistently appearing at similar spatial locations across subjects. AutoAtlas also produces very low dimensional features that represent local texture of each partition. We demonstrate prediction of metadata associated with each subject using the derived feature representations and compare the results to prediction using features derived from FreeSurfer anatomical parcellation. Since our features are intrinsically linked to distinct partitions, we can then map values of interest, such as partition-specific feature importance scores onto the brain for visualization.

42 ENGINEERING↗

Large Eddy Simulation of Gasoline Sprays in a Lagrangian–Eulerian Framework Using the High-Order Spectral Element Method

Predicting the spray evolution using simulations requires accurate modeling of the turbulent gas-phase flow field. Here, in this study, the high-order spectral-element method (SEM), implemented in the code Nek5000, was used to provide highly resolved solutions to the turbulent flow field. Spray modeling capabilities were implemented into the Nek5000 code. The spray is modeled in a Lagrangian–Eulerian (LE) framework, where the liquid is represented by discrete parcels of droplets. The method for coupling liquid and gas in the context of SEM is described, which allows for very fine meshes to be used without affecting the stability of the solution. Large-eddy simulations (LES) of the eight-hole ECN Spray G gasoline injector were conducted. Numerical results are compared against experimental data for liquid penetration, droplet size and gas velocity. The morphology of the multiplume spray is compared against experimental data. The effect of different spray injection inputs is analyzed. It was found that using a plume direction of 33 deg and an injection cone angle of 30 deg produced the best results overall. This work shows the applicability of SEM for spray modeling applications, where use of a high-order flow solver can help us understand the multiplume spray aerodynamics and how it leads to plume collapse under certain conditions. Results also highlight the need for tuning spray input parameters in the LE framework, even when high-fidelity gas flow solutions are possible.

33 ADVANCED PROPULSION SYSTEMS↗

Field To Farm Aggregation For Agricultural Systems

The Fields to Farms methodology illustrates the generation of farm parcels from the Crop Data Layer (CDL), a raster dataset containing 133 categories representing various crop types and land uses. This methodology involves two primary steps: Field Delineation and Farm Aggregation. The code specifically addresses the aggregation of pre-delineated fields within a county to form farms, adhering to predefined criteria for farm size categories. It is assumed that the field delineation process, which involves creating vector polygons from CDL raster, has been completed beforehand, possibly through external tools or methods. Upon initialization, the script processes county-level fields, preparing them for farm aggregation. In the Farm Aggregation phase, the code iteratively combines delineated fields into farms based on specified criteria, continuing until the aggregated farm size meets predefined thresholds derived from data from the 2017 National Agricultural Statistics Service (NASS) census. Throughout this iterative process, the script dynamically adjusts the aggregation to ensure alignment with the desired distribution reported by NASS. The resulting output of the script is a GeoDataFrame containing classified farms, which are subsequently saved as GeoPackage files. These files enable further analysis and visualization, facilitating comprehensive exploration of the farm landscape generated through the methodology.

Paudel, Rajiv [Idaho National Laboratory (INL), Id↗

Deep Learning Parameterization of Vertical Wind Velocity Variability via Constrained Adversarial Training

Atmospheric models with typical resolution in the tenths of kilometers cannot resolve the dynamics of air parcel ascent, which varies on scales ranging from tens to hundreds of meters. Small-scale wind fluctuations are thus characterized by a subgrid distribution of vertical wind velocity W with standard deviation σ W . The parameterization of σ W is fundamental to the representation of aerosol–cloud interactions, yet it is poorly constrained. Using a novel deep learning technique, this work develops a new parameterization for σ W merging data from global storm-resolving model simulations, high-frequency retrievals of W , and climate reanalysis products. The parameterization reproduces the observed statistics of σ W and leverages learned physical relations from the model simulations to guide extrapolation beyond the observed domain. Incorporating observational data during the training phase was found to be critical for its performance. The parameterization can be applied online within large-scale atmospheric models, or offline using output from weather forecasting and reanalysis products.

54 ENVIRONMENTAL SCIENCES↗

Numerical Simulations of Orographic Convection across Multiple Gray Zones

Idealized simulations are used to determine the sensitivity of moist orographic convection to horizontal grid spacing Δ h . In simulated mechanically (MECH) and thermally (THERM) forced convection over an isolated ridge, Δ h is varied systematically over both the deep-convection (Δ h ~ 10 –1 km) and turbulence (Δ h ~ 1 km–100 m) gray zones. To aid physical interpretation, a new parcel-based bulk entrainment/detrainment diagnosis for horizontally heterogeneous flows is developed. Within the deep-convection gray zone, the Δ h sensitivity is dominated by differences in parameterized versus explicit convection; the former initiates convection too far upstream of the ridge (MECH) and too early in the diurnal heating cycle (THERM). These errors stem in part from a large underprediction of parameterized entrainment and detrainment. Within the turbulence gray zone, sensitivities to Δh arise from the representation of both subcloud- and cloud-layer turbulence. As Δ h is decreased, MECH exhibits stronger cloud-layer entrainment to enhance the convective mass flux M co , while THERM exhibits stronger detrainment to suppress Mco and delay convection initiation. The latter is reinforced by increased subcloud turbulence at smaller Δ h , which leads to drying and diffusion of the central updraft responsible for initiating moist convection. Numerical convergence to a robust solution occurs only in THERM, which develops a fully turbulent flow with a resolved inertial subrange (for Δ h ≤ 250 m). In MECH, by contrast, turbulent transition occurs within the orographic cloud, the details of which depend on both physical location and Δ h .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluating the Lagrangian Evolution of Subtropical Low Clouds in GCMs Using Observations: Mean Evolution, Time Scales, and Responses to Predictors

In this work, a Lagrangian framework is developed to show the daily-scale time evolution of low clouds over the eastern subtropical oceans. An identical framework is applied to two general circulation models (GCMs): the CAM5 and UKMET and a set of satellite observations. This approach follows thousands of parcels as they advect downwind in the subtropical trade winds, comparing cloud evolution in time and space. This study tracks cloud cover, in-cloud liquid water path (CLWP), droplet concentration $N_d$, planetary boundary layer (PBL) depth, and rain rate as clouds transition from regions with predominately stratiform clouds to regions containing mostly trade cumulus. The two models generate fewer clouds with greater $N_d$ relative to observations. Models show stronger Lagrangian cloud cover decline and greater PBL deepening when compared with observations. In comparing frequency distributions of cloud variables over time, it is seen that models generate increasing frequencies of nearly clear conditions at the expense of overcast conditions, whereas observations show transitions from overcast to cloud amounts between 50% and 90%. Lagrangian decorrelation time scales ($\mathcal{e}$-folding time $\tau$) of cloud cover and CLWP are between 11 and 19 h for models and observations, although they are a bit shorter for models. A Lagrangian framework applied here resolves and compares the time evolution of cloud systems as they adjust to environmental perturbations in models and observations. Increasing subsidence in the overlying troposphere leads to declining cloud cover, CLWP, PBL depth, and rain rates in models and observations. Modeled cloud responses to other meteorological variables are less consistent with observations, suggesting a need for continuing mechanical improvements in GCMs.

58 GEOSCIENCES↗

An evaluation of size-resolved cloud microphysics scheme numerics for use with radar observations Part II: Condensation and evaporation

Accurate numerical modeling of clouds and precipitation is essential for weather forecasting and climate change research. While size-resolved (bin) cloud microphysics models predict particle size distributions without imposing shapes, results are subject to artificial size distribution broadening owing to numerical diffusion associated with various processes. Whereas Part I of this study addressed collision–coalescence, here we investigate numerical diffusion that occurs in solving condensation and evaporation. In a parcel model framework, all of the numerical schemes examined converge to one solution of condensation and evaporation as the mass grid is refined, and the advection-based schemes are recommended over the reassigning schemes. Including Eulerian vertical advection in a column limits the convergence to some extent, but that limitation occurs at a sufficiently fine mass grid, and the number of iterations in solving vertical advection should be minimized to reduce numerical diffusion. Insubstantial numerical diffusion in solving condensation can be amplified if collision–coalescence is also active, which in turn can be substantially diminished if turbulence effects on collision are incorporated. Here, large-eddy simulations of a drizzling stratocumulus field reveal that changes in moments of Doppler spectra obtained using different mass grids are consistent with those obtained from the simplified framework, and that spectral moments obtained using a mass grid designed to effectively reduce numerical diffusion are generally closer to observations. Notable differences between the simulations and observations still exist, and our results suggest a need to consider whether factors other than numerical diffusion in the fundamental process schemes employed can cause such differences.

54 ENVIRONMENTAL SCIENCES↗

New Empirical Formulation for the Sublimational Breakup of Graupel and Dendritic Snow

Abstract Ice fragments are generated by sublimation of ice particles in subsaturated conditions in natural clouds. Conceivably, such sublimational breakup would be expected to cause ice multiplication in natural clouds. Any fragment that survives will grow to become ice precipitation that may sublimate and fragment further. As a first step toward assessing this overlooked process, a formulation is proposed for the number of ice fragments from sublimation of ice particles for an atmospheric model. This is done by amalgamating laboratory observations from previously published studies. The concept of a “sublimated mass activity spectrum” for the breakup is applied to the dataset. The number of ice fragments is determined by the relative humidity over ice and the initial size of the parent ice particles. The new formulation applies to dendritic crystals and heavily rimed particles only. Finally, a thought experiment is performed for an idealized scenario of subsaturation with in-cloud descent. Scaling analysis yields an estimate of an ice enhancement ratio of about 5 (10) within a weak deep convective downdraft of about 2 m s −1 , for an initial monodisperse population of dendritic snow (graupel) particles of 3 L −1 and 2 mm. During descent, there is a dynamic equilibrium between continual emission of fragments and their depletion by sublimation. A simplified bin microphysics parcel model exhibits this dynamical quasi equilibrium, consistent with the thought experiment. The fragments have average lifetimes of around 90 and 70 s for dendrites and graupel, respectively. Sublimational breakup is predicted to cause significant secondary ice production.

54 ENVIRONMENTAL SCIENCES↗

Semi-Analytic Functions to Calculate the Deposition Coefficients for Ice Crystal Vapor Growth in Bin and Bulk Microphysical Models

Numerical cloud models require estimates of the vapor growth rate for ice crystals. Current bulk and bin microphysical parameterizations generally assume that vapor growth is diffusion limited, though some parameterizations include the influence of surface attachment kinetics through a constant deposition coefficient. A parameterization for variable deposition coefficients is provided herein. The parameterization is an explicit function of the ambient ice supersaturation and temperature, and an implicit function of crystal dimensions and pressure. The parameterization is valid for variable surface types including growth by dislocations and growth by step nucleation. Deposition coefficients are predicted for the two primary growth directions of crystals, allowing for the evolution of the primary habits. Comparisons with benchmark calculations of instantaneous mass growth indicate that the parameterization is accurate to within a relative error of 1%. Parcel model simulations using Lagrangian microphysics as a benchmark indicate that the bulk parameterization captures the evolution of mass mixing ratio and fall speed with typical relative errors of less than 10%, whereas the average axis lengths can have errors of up to 20%. The bin model produces greater accuracy with relative errors often less that 10%. The deposition coefficient parameterization can be used in any bulk and bin scheme, with low error, if an equivalent volume spherical radius is provided.

54 ENVIRONMENTAL SCIENCES↗

Evaluating the Conservation of Energy Variables in Simulations of Deep Moist Convection

Abstract It is often assumed in parcel theory calculations, numerical models, and cumulus parameterizations that moist static energy (MSE) is adiabatically conserved. However, the adiabatic conservation of MSE is only approximate because of the assumption of hydrostatic balance. Two alternative variables are evaluated here: MSE − IB and MSE + KE, wherein IB is the path integral of buoyancy ( B ) and KE is kinetic energy. Both of these variables relax the hydrostatic assumption and are more precisely conserved than MSE. This article quantifies the errors that result from assuming that the aforementioned variables are conserved in large-eddy simulations (LES) of both disorganized and organized deep convection. Results show that both MSE − IB and MSE + KE better predict quantities along trajectories than MSE alone. MSE − IB is better conserved in isolated deep convection, whereas MSE − IB and MSE + KE perform comparably in squall-line simulations. These results are explained by differences between the pressure perturbation behavior of squall lines and isolated convection. Errors in updraft B diagnoses are universally minimized when MSE − IB is assumed to be adiabatically conserved, but only when moisture dependencies of heat capacity and temperature dependency of latent heating are accounted for. When less accurate latent heat and heat capacity formulae were used, MSE − IB yielded poorer B predictions than MSE due to compensating errors. Our results suggest that various applications would benefit from using either MSE − IB or MSE + KE instead of MSE with properly formulated heat capacities and latent heats.

54 ENVIRONMENTAL SCIENCES↗

Theory of In-Cloud Activation of Aerosols and Microphysical Quasi-Equilibrium in a Deep Updraft

The microphysical quasi-equilibrium in an ascending adiabatic parcel is elucidated by an analytical theory in 0D with drastic simplifications. The theory predicts how in-cloud activation is most likely to be triggered by the onset of precipitation during sufficient ascent, with the ascent only needing to approach almost twice the cloud-base updraft speed aloft. The precipitation and cloud condensate mass fields are coupled in a closed system in a simplified version of the model. The initial state of no precipitation is unstable with respect to a perturbation. In the 2D phase space of both mass fields, there is a neutral line. Unstable growth of precipitation mass drives the system to cross the line into a regime of stability. An attractor is then approached consisting of precipitation mass balanced by accretion of cloud mass and fallout. The cloud-particle number concentration also approaches a stable equilibrium governed by the balance between in-cloud activation aloft from the inexorably increasing supersaturation during vertical acceleration and accretion of cloud droplets by precipitation. Dimensionless numbers characterizing the microphysical equilibria and their stability are derived mathematically, including a condensation–precipitation efficiency and an in-cloud activation efficiency. The theory explains common observations of the orders of magnitude of liquid water content in convective and stratiform clouds. Finally, sensitivity tests of the numerically integrated theoretical equations are documented with respect to variations in cloud condensation nucleus (CCN) aerosols and updraft speed. In conclusion, it is shown that this theory of in-cloud activation, with an increasing supersaturation during ascent, applies to both ice-only and liquid-only cloud.

54 ENVIRONMENTAL SCIENCES↗

Sources of stochasticity in the growth of cloud droplets: supersaturation fluctuations versus turbulent transport

The role played by fluctuations of supersaturation in the growth of cloud droplets is examined in this study. The stochastic condensation framework and the three regimes of activation of cloud droplets – namely, mean-dominant, fluctuation-influenced, and fluctuation-dominant, are used for analyzing the data from high-resolution large-eddy simulations of the Pi convection-cloud chamber. Based on a detailed budget analysis the significance of all the terms in the evolution of the droplet size distribution equation is evaluated in all three regimes. The analysis indicates that the mean-growth rate is a dominant process in shaping the droplet size distribution in all three regimes. Turbulence introduces two sources of stochasticity, turbulent transport and particle lifetime, and supersaturation fluctuations. The transport of cloud droplets plays an important role in all three regimes, whereas the direct effect of supersaturation fluctuations is primarily related to the activation and growth of the small droplets in the fluctuation-influenced and fluctuation-dominant regimes. Here, we compare our results against the previous studies (experimental and theory) of the Pi chamber, and discuss the limitations of the existing models based on the stochastic condensation framework. Furthermore, we extend the discussion of our results to atmospheric clouds, and in particular focus on recent adiabatic turbulent cloud parcel simulations based on the stochastic condensation framework, and emphasize the importance of entrainment/mixing and turbulent transport in shaping the droplet size distribution.

54 ENVIRONMENTAL SCIENCES↗

The Dual Nature of Entrainment-Mixing Signatures Revealed through Large-Eddy Simulations of a Convection-Cloud Chamber

Abstract Entrainment of subsaturated air into a cloud can influence its optical and microphysical properties in various ways, depending on the droplet evaporation and turbulent mixing time scales. Previous experiments in the Pi convection-cloud chamber have revealed that, given a fixed entrained air property, the mixing of entrained subsaturated air results in complete evaporation of some cloud droplets, with the rest remaining unchanged. This is a signature of inhomogeneous mixing. While comparing the results of entrainment with varying air properties, the mixing signature appears as if the subsaturated air is well mixed with the cloud to evenly reduce the droplets’ size. In other words, taken together, the experiments appear to have the signature of homogeneous mixing. To explore these results in a greater depth, we conduct large-eddy simulations combined with a bin microphysics scheme. Our results reproduce the similar signatures of inhomogeneous and homogeneous mixing, implying that LES can resolve the inhomogeneous mixing when the grid spacing is smaller than the entrained air parcel. Additionally, we observe that increasing the aerosol injection rate enhances the signature of inhomogeneous mixing, while coarser grid spacing diminishes it. Finally, the change in wall fluxes in response to various entrained air properties confirms that the homogeneous signature seen in the analysis of an ensemble of simulations is the result of various equilibrium states. This further strengthens the suggestion that the homogeneous mixing signature found in aircraft observations near the cloud top may result from combining entrainment events of different intensities, possibly caused by various-sized eddies. Significance Statement Large-eddy simulation and size-resolved microphysics can resolve time scales for turbulent mixing and evaporation and, therefore, are well suited for reproducing, extending, and interpreting the entrainment experiment in the Pi convection-cloud chamber. Our simulation results confirm (i) the inhomogeneous mixing signature for an individual entrainment event and (ii) the appearance of homogeneous mixing in an ensemble of entrainment episodes. Furthermore, we demonstrate that the inhomogeneous mixing signature is more pronounced in a polluted cloud, but coarser grid spacing in simulations may compromise the accuracy of this signature. Last, the homogeneous mixing signature results from various equilibrium states established for different entrainment intensities and adjusted wall fluxes, which are challenging to measure experimentally but can be easily analyzed in the simulations.

54 ENVIRONMENTAL SCIENCES↗

Entrainment, Detrainment, and Dilution of Dry and Moist Atmospheric Thermals

Here this study examines the entrainment, detrainment, and dilution of dry and moist (cloud) atmospheric thermals in large-eddy simulations. In a neutrally stable environment (with respect to dry dynamics), moist thermals have an increase in radius R with thermal height z t (α ≡ dR/dz t ) about 4 times smaller compared to dry thermals when density stratification is considered and ~2.4 times smaller without density stratification (i.e., applying the Boussinesq approximation). An analytic expression relating α to several dimensionless parameters is derived from the thermal impulse–circulation relation to clarify the factors impacting α. This expression shows that the difference in buoyancy structure between moist and dry thermals, with buoyancy concentrated in the central cores of moist thermals owing to latent heating, explains their smaller spreading rates. Individual contributions of entrainment and detrainment are analyzed using a direct parcel-based approach in the simulations. Moist thermals have similar fractional detrainment but much smaller fractional entrainment rates compared to dry thermals, consistent with the differences in α. Despite having smaller α, moist thermals are similarly dilute (quantified by a passive tracer) as dry thermals because of their greater mixing efficiency with the environment. Thus, moist thermals are substantially dilute but expand much less in size/volume as they rise compared to dry thermals. The α values for moist thermals in (dry) neutral and statically stable environments are similar, but fractional entrainment and especially detrainment rates are greater in the stable environment. Large detrainment rates are associated with a breakdown of the broader thermal vortex ring structure, especially with low environmental relative humidity, attributed in part to evaporation and buoyancy reversal.

54 ENVIRONMENTAL SCIENCES↗

Changes in Convective Available Potential Energy and Convective Inhibition under Global Warming

Atmospheric convective available potential energy (CAPE) is expected to increase under greenhouse gas–induced global warming, but a recent regional study also suggests enhanced convective inhibition (CIN) over land although its cause is not well understood. In this study, a global climate model is first evaluated by comparing its CAPE and CIN with reanalysis data, and then their future changes and the underlying causes are examined. The climate model reasonably captures the present-day CAPE and CIN patterns seen in the reanalysis, and projects increased CAPE almost everywhere and stronger CIN over most land under global warming. Over land, the cases or times with medium to strong CAPE or CIN would increase while cases with weak CAPE or CIN would decrease, leading to an overall strengthening in their mean values. These projected changes are confirmed by convection-permitting 4-km model simulations over the United States. The CAPE increase results mainly from increased low-level specific humidity, which leads to more latent heating and buoyancy for a lifted parcel above the level of free convection (LFC) and also a higher level of neutral buoyancy. The enhanced CIN over most land results mainly from reduced low-level relative humidity (RH), which leads to a higher lifting condensation level and a higher LFC and thus more negative buoyancy. Over tropical oceans, the near-surface RH increases slightly, leading to slight weakening of CIN. Over the subtropical eastern Pacific and Atlantic Ocean, the impact of reduced low-level atmospheric lapse rates overshadows the effect of increased specific humidity, leading to decreased CAPE.

54 ENVIRONMENTAL SCIENCES↗