SEARCH · Engineering Papers

Results for “Ensemble simulations”

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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Regional Earthquake Ground Motion Simulations for Southern California With EQSIM: Insights From the 2008 Chino Hills, 2024 Highland Park, and 2021 Carson Earthquakes

This study presents physics-based, 3D simulations using the EQSIM framework for several earthquakes in the Los Angeles region. The primary objective was to assess the ability of deterministic physics-based ground motion simulations to reproduce the observed motions from historical events. The selected events included the mathematical equation M w 5.4 2008 Chino Hills, the mathematical equation M w 4.4 2024 Highland Park, and the mathematical equation M w 4.3 2021 Carson events. The simulated motions were evaluated by comparing the recorded and simulated seismograms, as well as the Fourier amplitude spectra, across multiple seismic stations. The SCEC 3D velocity model, CVM-S4.26.M01, was used to represent the regional geology, and ground motion simulations were carried out with a resolution of up to 5 Hz. The results indicate that the simulated motions captured the recorded motions up to approximately 4 Hz. While careful iterations regarding source parameters and corner frequencies were required, and, for the case of the Highland Park event, some of the near-source stations had relatively low accuracy, the present study established a positive step toward the utilization of physics-based simulations in practical applications. The computational efficiencies exhibited by EQSIM, especially on GPU clusters, further supported this assertion, as wall-clock times of simulations involving more than 10 billion grid points were as low as mathematical equation minutes. This permits ensemble simulations for a considered scenario event so that modeling uncertainties (e.g., source and geology) can be bracketed.

EQSIM↗

DOE OSTI2020

FATES size- and age-dependent mortality simulation output analysis scripts

This dataset contains python scripts to analyse outputs from the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) and accompanies the paper “Needham, J.F., Chambers, J., Fisher, R., Knox, R., and Koven, C. D., Forest responses to simulated elevated CO2 under alternate hypotheses of size- and age-dependent mortality, 2020, Global Change Biology”. These scripts process FATES outputs from size- and age-dependent mortality simulations which were run to test the effect of different mechanisms of mortality on forest response to elevated CO2 (eCO2). Specifically, these scripts will process single plant functional type (PFT) simulations in which mortality is either a constant background rate, size-dependent or age-dependent. In each case, a simulation with constant woody NPP is compared to a simulation in which woody NPP increases by 25% to simulate the growth response of forests to eCO2. In addition, the data package contains scripts to analyse ensemble simulations with size- and age-dependent mortality and two PFTs, that were run to test the impact of different demographic rates on coexistence and the forest response to eCO2. Finally, the dataset includes scripts for processing simulations testing the sensitivity of results to allometry and to the recruitment scheme.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Robust projection of East Asian summer monsoon rainfall based on dynamical modes of variability

Abstract The Asian monsoon provides the freshwater that a large population in Asia depends on, but how anthropogenic climate warming may alter this key water source remains unclear. This is partly due to the prevailing point-wise assessment of climate projections, even though climate change patterns are inherently organized by dynamics intrinsic to the climate system. Here, we assess the future changes in the East Asian summer monsoon precipitation by projecting the precipitation from several large ensemble simulations and CMIP6 simulations onto the two leading dynamical modes of internal variability. The result shows a remarkable agreement among the ensembles on the increasing trends and the increasing daily variability in both dynamical modes, with the projection pattern emerging as early as the late 2030 s. The increase of the daily variability of the modes heralds more monsoon-related hydrological extremes over some identifiable East Asian regions in the coming decades.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2011

MoSST DAS: The First Working Geomagnetic Data Assimilation System

The Earth possesses an internal magnetic field (geomagnetic field) generated by convection in the outer core (geodynamo). Previous efforts have been focused along two distinct paths: (1) numerical geodynamo modeling to understand the origin of the geomagnetic field, and the mechanisms of geomagnetic secular variations (SV); and (2) geomagnetic field modeling to map the spatial/temporal variations of the field from geomagnetic data, and to derive core properties, e.g. inversion of core flow near the core-mantle boundary (CMB). Geomagnetic data assimilation is a new approach emerged over the past 5 years: surface observations are assimilated with geodynamo models for better understanding of the core dynamical state, and accurately prediction of SV. In collaboration with several geomagnetic research groups, we have developed the first working geomagnetic data assimilation system, Modular, Scalable, Self-consistent, and Three-dimensional (MoSST) DAS, that includes the MoSST numerical dynamo model; 7000 years of geomagnetic field maps from several field models utilizing satellite and ground observatory data, historical magnetic records and archeo/paleo magnetic data; and an ensemble based optimal interpolation (01) assimilation algorithm. With this system, we have demonstrated clearly that the assimilated core dynamical state is substantially different from those of pure geodynamo simulations. Ensemble assimilation runs also show the convergence of the assimilated solutions inside the core, suggesting that the simulation state is pulled closer to the truth via data assimilation. The forecasts from this system are also very accurate: the 5-year forecast of the geomagnetic field agrees very well with the observations; and the 5-year secular variation forecast is more accurate than the IGRF SV forecast models in the past. Using geomagnetic records up to 2009, we have made an SV forecast for the period from 2010-2015, and is a candidate SV model for IGRF-11.

Kuang, Weijia↗

NASA NTRS2017

El Nino-Induced Tropical Ocean/Land Energy Exchange in MERRA-2 and M2AMIP

Studies have shown the correlation and connection of surface temperatures across the globe, ocean and land, related to Tropical SSTs especially El Nino. This climate variability greatly influences regional weather and hydroclimate extremes (e.g. drought and flood). In this paper, we evaluate the relationship of temperatures across the tropical oceans and continents in MERRA-2, and also in a newly developed MERRA-2 AMIP ensemble simulation (M2AMIP). M2AMIP uses the same model and spatial resolution as MERRA-2, producing the same output diagnostics over 10 ensemble members. Composite El Nino temperature data are compared with observations to evaluate the land/sea contrast, variations and phase relationship. The temperature variations are related to surface heat fluxes and the atmospheric temperatures and transport, to identify the processes that lead to the lagged redistribution of heat in the tropics and beyond. Discernable cloud, radiation and data assimilation changes accompany the onset of El Nino affecting continental regions through the progression to and following the peak values. While the model represents these variations in general, regional strengths and weaknesses can be identified.

Bosilovich, Michael G.↗

DOE OSTI2024

Evaluation of a Stochastic Mixing Scheme in the Deep Convective Gray Zone Using a Tropical Oceanic Deep Convection Case Study

Abstract A stochastic horizontal subgrid‐scale mixing scheme is evaluated in ensemble simulations of a tropical oceanic deep convection case using a horizontal grid spacing (Δ h ) of 3 km. The stochastic scheme, which perturbs the horizontal mixing coefficient according to a prescribed spatiotemporal autocorrelation scale, is found to generally increase mesoscale organization and convective intensity relative to a non‐stochastic control simulation. Perturbations applied at relatively short autocorrelation scales induce differences relative to the control that are more systematic than those from perturbations applied at relatively long scales that yield more variable outcomes. A simulation with mixing enhanced by a constant factor of 4 significantly increases mesoscale organization and convective intensity, while turning off horizontal subgrid‐scale mixing decreases both. Total rainfall is modulated by a combination of mesoscale organization, areal coverage of convection, and convective intensity. The stochastic simulations tend to behave more similarly to the constant enhanced mixing simulation owing to greater impacts from enhanced mixing as compared to reduced mixing. The impacts of stochastic mixing are robust, ascertained by comparing the stochastic mixing ensembles with a non‐stochastic mixing ensemble that has grid‐scale noise added to the initial thermodynamic field. Compared to radar observations and a higher resolution Δ h = 1 km simulation, stochastic mixing seemingly degrades the simulation performance. These results imply that stochastic mixing produces non‐negligible impacts on convective system properties and evolution but does not lead to an improved representation of convective cloud characteristics in the case studied here.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Roles of External Forcing and Internal Variability in Precipitation Changes of a Sub-Region of the U.S. Mid-Atlantic During 1979–2019

Observations of the seasonal mean precipitation averaged over a sub-region of the U.S. mid-Atlantic show significant increases in summer but no changes in winter during 1979–2019. Clarifying the drivers of these observed changes is important for understanding and interpreting future projections of regional precipitation by climate models. Herein, this study examines the relative effects of external forcing and internal variability in the observed multi-decadal changes using the historical and RCP8.5/SSP585 simulations from 37 models in the Coupled Model Intercomparison Project phase 5 (CMIP5), 35 models in CMIP6, as well as large ensemble simulations from 6 models, each including more than 20 members. During 1979–2019, enhanced warming over land relative to ocean due to external forcing induces low pressure anomalies over land that contribute to the wetting summer. While similar wetting also occurs during winter, such effect is offset by the drying associated with the positive phase of the North Atlantic Oscillation (NAO) in recent decades. Internal variability produces large uncertainty in the multi-decadal precipitation changes simulated by climate models through changes in large-scale circulation associated with the NAO in winter and the North Atlantic Subtropical High in summer. These multi-decadal relationships between precipitation and large-scale circulation under internal variability also manifest on inter-annual timescale. These relationships highlight enhanced moisture transport from the Gulf of Mexico and North Atlantic in winter and near the mid-Atlantic coast in summer that modulates precipitation change in the study region.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Assessing mean climate change signals in the global CORDEX-CORE ensemble

The new Coordinated Output for Regional Evaluations (CORDEX-CORE) ensemble provides high-resolution, consistent regional climate change projections for the major inhabited areas of the world. It serves as a solid scientific basis for further research related to vulnerability, impact, adaptation and climate services in addition to existing CORDEX simulations. The aim of this study is to investigate and document the climate change information provided by the CORDEX-CORE simulation ensemble, as a part of the World Climate Research Programme (WCRP) CORDEX community. An overview of the annual and monthly mean climate change information in selected regions in different CORDEX domains is presented for temperature and precipitation, providing the foundation for detailed follow-up studies and applications. Initially, two regional climate models (RCMs), REMO and RegCM were used to downscale global climate model output. The driving simulations by AR5 global climate models (AR5-GCMs) were selected to cover the spread of high, medium, and low equilibrium climate sensitivity at a global scale. The CORDEX-CORE ensemble has doubled the spatial resolution compared to the previously existing CORDEX simulations in most of the regions (25 $\,\mathrm {km}$ (0.22 $^{\circ }$ /jats:inline-formula>) versus 50 $\,\mathrm {km}$ (0.44 $^{\circ }$ )) leading to a potentially improved representation of, e.g., physical processes in the RCMs. The analysis focuses on changes in the IPCC physical climate reference regions. The results show a general reasonable representation of the spread of the temperature and precipitation climate change signals of the AR5-GCMs by the CORDEX-CORE simulations in the investigated regions in all CORDEX domains by mostly covering the AR5 interquartile range of climate change signals. The simulated CORDEX-CORE monthly climate change signals mostly follow the AR5-GCMs, although for specific regions they show a different change in the course of the year compared to the AR5-GCMs, especially for RCP8.5, which needs to be investigated further in region specific process studies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Estimating the Likelihood of GHG Concentration Scenarios From Probabilistic Integrated Assessment Model Simulations

The climate scenarios that form the basis for current climate risk assessments have no assigned probabilities, and this impedes the analysis of future climate risks. This paper proposes an approach to estimate the probability of carbon dioxide (CO 2 ) concentration scenarios used in key climate change modeling experiments. It computes the CO 2 emissions compatible with the concentrations prescribed by Coupled Model Intercomparison Project Phase 5 (CMIP5) and CMIP6 experiments. The distribution of these compatible cumulative emissions is interpreted as the likelihood of future emissions given a concentration pathway. Using Bayesian analysis, the probability of each pathway can be estimated from a probabilistic sample of future emissions. The approach is demonstrated with five probabilistic CO 2 emission simulation ensembles from four Integrated Assessment Models (IAM), leading to independent estimates of the likelihood of the CO 2 concentration of Representative Concentration Pathways (RCP) and Shared Socioeconomic Pathways (SSP). Results suggest that SSP5-8.5 is unlikely for the second half of the 21st century, but offer no clear consensus on which of the remaining scenarios is most likely. Estimates of likelihoods of CO 2 concentrations associated with RCP and SSP scenarios are affected by sampling errors, differences in emission sources simulated by the IAMs, and a lack of a common experimental framework for IAM simulations. These shortcomings, along with a small IAM ensemble size, limit the applicability of the results presented here. Novel joint IAM and the Earth System Model experiments are needed to deliver actionable probabilistic climate risk assessments.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

GrainNN: A neighbor-aware long short-term memory network for predicting microstructure evolution during polycrystalline grain formation

High fidelity simulations of grain formation in alloys are an indispensable tool for process-to-mechanical-properties characterization. Such simulations, however, can be computationally expensive as they require fine spatial and temporal discretizations. Their cost becomes an obstacle to parametric studies and ensemble runs and ultimately makes downstream tasks like optimal control and uncertainty quantification challenging. To enable such downstream tasks, we introduce GrainNN, an efficient and accurate reduced-order model for epitaxial grain growth in additive manufacturing conditions. GrainNN is a sequence-to-sequence long-short-term-memory (LSTM) deep neural network that evolves the dynamics of manually crafted features. Its innovations are (1) an attention mechanism with grain-microstructure-specific transformer architecture; and (2) an overlapping combination of several clones of the network to generalize to grain configurations that are different from those used for training. This design enables GrainNN to predict grain formation for unseen physical parameters, grain number, domain size and geometry. Furthermore, GrainNN not only reconstructs the quantities of interest but also can be pointwise accurate. In our numerical experiments, we use a polycrystalline phase field method to both generate the training data and assess GrainNN. For multiparametric, ensemble simulations with many grains, GrainNN can be orders of magnitude faster than phase field simulations, while delivering 5%–15% pointwise error. Additionally, this speedup includes the cost of the phase field simulations for generating training data.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Accelerated Probabilistic Marching Cubes by Deep Learning for Time-Varying Scalar Ensembles

Visualizing the uncertainty of ensemble simulations is challenging due to the large size and multivariate and temporal features of en-semble data sets. One popular approach to studying the uncertainty of ensembles is analyzing the positional uncertainty of the level sets. Probabilistic marching cubes is a technique that performs Monte Carlo sampling of multivariate Gaussian noise distributions for positional uncertainty visualization of level sets. However, the technique suffers from high computational time, making interactive visualization and analysis impossible to achieve. This paper introduces a deep-learning-based approach to learning the level-set uncertainty for two-dimensional ensemble data with a multivariate Gaussian noise assumption. We train the model using the first few time steps from time-varying ensemble data in our workflow. We demonstrate that our trained model accurately infers uncertainty in level sets for new time steps and is up to 170X faster than that of the original probabilistic model with serial computation and 10X faster than that of the original parallel computation.

Han, Mengjiao↗

NASA NTRSDate not supplied

The Role of Internal Variability and Feedbacks Controlling AMOC Stability

A bi-stable mode of the Atlantic Meridional Overturning Circulation (AMOC) is found in a 10-member ensemble simulation of the SSP2-4.5 scenario using the NASA GISS-E2-1-G climate model. Local feedbacks in the subpolar North Atlantic region in conjunction with internal variability in sea-ice transport and melt play a critical role in causing the divergent behavior of the AMOC in the ensemble members. While other fully coupled models have demonstrated the important role of surface freshening in leading to AMOC shutdown, either through hosing experiments or increased precipitation and greenhouse gas warming at high latitudes, in the GISS simulations, there are no external freshwater perturbations. This is the first time that a CMIP-class model has shown such a bifurcation across an initial condition ensemble.

Atlantic Meridional Overturning Circulation↗

DOE OSTI2023

Ice Supersaturation Variability in Cirrus Clouds: Role of Vertical Wind Speeds and Deposition Coefficients

Abstract Aircraft measurements reveal ice supersaturation statistics in cirrus (ISSs) with broad maxima around ice saturation and pronounced variance. In this study, processes shaping ISSs in midlatitude and tropical upper tropospheric conditions are systematically investigated. Water vapor deposition and sublimation of size‐resolved ice crystal populations are simulated in an air parcel framework. Mesoscale temperature fluctuations (MTFs) due to gravity waves force the temporal evolution of supersaturation. Various levels of background wave forcing and cirrus thickness are distinguished in stochastic ensemble simulations. Kinetic limitations to ice mass growth are brought about by supersaturation‐dependent deposition coefficients that represent efficient and inefficient growth modes as a function of ice crystal size. The simulations identify a wide range of deposition coefficients in cirrus, but most values stay above 0.01 such that kinetic limitations to water uptake remain moderate. Supersaturation quenching times are long, typically 0.5–2 hr. The wave forcing thus causes a remarkably large variability in ISSs and cirrus microphysical properties except in the thickest cirrus, producing ensemble‐mean ISSs in line with in‐situ measurements. ISS variance is controlled by MTFs and increases with decreasing cirrus integral radii. In comparison, the impact of ice crystal growth rates on ISSs is small. These results contribute to efforts directed at identifying and solving issues associated with ice‐supersaturated areas and non‐equilibrium cirrus physics in global models.

Meteorology & Atmospheric Sciences↗

DOE OSTI2024

Probabilistic projections of the Amery Ice Shelf catchment, Antarctica, under conditions of high ice-shelf basal melt

Abstract. Antarctica's Lambert Glacier drains about one-sixth of the ice from the East Antarctic Ice Sheet and is considered stable due to the strong buttressing provided by the Amery Ice Shelf. While previous projections of the sea-level contribution from this sector of the ice sheet have predicted significant mass loss only with near-complete removal of the ice shelf, the ocean warming necessary for this was deemed unlikely. Recent climate projections through 2300 indicate that sufficient ocean warming is a distinct possibility after 2100. This work explores the impact of parametric uncertainty on projections of the response of the Lambert–Amery system (hereafter “the Amery sector”) to abrupt ocean warming through Bayesian calibration of a perturbed-parameter ice-sheet model ensemble. We address the computational cost of uncertainty quantification for ice-sheet model projections via statistical emulation, which employs surrogate models for fast and inexpensive parameter space exploration while retaining critical features of the high-fidelity simulations. To this end, we build Gaussian process (GP) emulators from simulations of the Amery sector at a medium resolution (4–20 km mesh) using the Model for Prediction Across Scales (MPAS)-Albany Land Ice (MALI) model. We consider six input parameters that control basal friction, ice stiffness, calving, and ice-shelf basal melting. From these, we generate 200 perturbed input parameter initializations using space filling Sobol sampling. For our end-to-end probabilistic modeling workflow, we first train emulators on the simulation ensemble and then calibrate the input parameters using observations of the mass balance, grounding line movement, and calving front movement with priors assigned via expert knowledge. Next, we use MALI to project a subset of simulations to 2300 using ocean and atmosphere forcings from a climate model for both low- and high-greenhouse-gas-emission scenarios. From these simulation outputs, we build multivariate emulators by combining GP regression with principal component dimension reduction to emulate multivariate sea-level contribution time series data from the MALI simulations. We then use these emulators to propagate uncertainty from model input parameters to predictions of glacier mass loss through 2300, demonstrating that the calibrated posterior distributions have both greater mass loss and reduced variance compared to the uncalibrated prior distributions. Parametric uncertainty is large enough through about 2130 that the two projections under different emission scenarios are indistinguishable from one another. However, after rapid ocean warming in the first half of the 22nd century, the projections become statistically distinct within decades. Overall, this study demonstrates an efficient Bayesian calibration and uncertainty propagation workflow for ice-sheet model projections and identifies the potential for large sea-level rise contributions from the Amery sector of the Antarctic Ice Sheet after 2100 under high-greenhouse-gas-emission scenarios.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

On the additivity of climate responses to the volcanic and solar forcing in the early 19th century

The early 19th century was the coldest period over the past 500 years, when strong tropical volcanic events and a solar minimum coincided. The 1809 unidentified eruption and the 1815 Tambora eruption happened consecutively during the Dalton minimum of solar irradiance; however, the relative role of the two forcing (volcano and solar) agents is still unclear. In this study, we examine the responses from a set of early 19th century simulations with combined and separated volcanic and solar forcing agents, as suggested in the protocol for the past1000 experiment of the Paleoclimate Modelling Intercomparison Project – Phase 4 (PMIP4). From 20-member ensemble simulations with the Max Planck Institute Earth system model (MPI-ESM1.2-LR), we find that the volcano- and solar-induced surface cooling is additive in the global mean/large scale, regardless of combining or separating the forcing agents. The two solar reconstructions (SATIRE (Spectral and Total Irradiance REconstruction-Millennia model) and PMOD (Physikalisch-Meteorologisches Observatorium Davos)) contribute to a cooling before and after 1815 of ~0.05 and ~0.15 K monthly average near-surface air cooling, respectively, indicating a limited solar contribution to the early 19th century cold period. The volcanic events provide the main cooling contributions, inducing a surface cooling that peaks at ~0.82 K for the 1809 event and ~1.35 K for Tambora. After the Tambora eruption, the temperature in most regions increases toward climatology largely within 5 years, along with the reduction of volcanic forcing. In the northern extratropical oceans, the temperature increases slowly at a constant rate until 1830, which is related to the reduction of seasonality and the concurrent changes in Arctic sea-ice extent. The albedo feedback of Arctic sea ice is found to be the main contributor to the Arctic amplification of the cooling signal. Several non-additive responses to solar and volcanic forcing happen on regional scales. In the atmosphere, the stratospheric polar vortex tends to strengthen when combining both volcano and solar forcing, even though the two forcing agents separately induce opposite-sign changes in stratospheric temperatures and zonal winds. In the ocean, when combining the two forcings, additional surface cold water propagates to the northern extratropics from the additional solar cooling in the tropics, which results in regional cooling along the propagation. Overall, this study not only quantifies the surface responses from combinations of the volcano and solar forcing, but also highlights the components that cannot be simply added from the responses of the individual forcing agents, indicating that a relatively small forcing agent (such as solar in early 19th century) can impact the response from the large forcing (such as the 1815 Tambora eruption) when considering regional climates.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Transfer learning of neural surrogates on multifidelity groundwater simulations

Multifidelity data used in the paper published in Advances in Water Resources 206 (2025) 105140, https://doi.org/10.1016/j.advwatres.2025.105140 The code used to process the data is openly available on GitHub at https://github.com/Model-Reduction-and-UQ-Group/Transfer_Learning_K_reconstruction Computationally inexpensive surrogates of process-based models, such as deep neural networks, enable ensemble-based computations used in risk assessment, data assimilation, etc. However, generation of large datasets required to train a neural network can be as expensive as the ensemble simulations themselves. We ameliorate this challenge by using data from multifidelity (MF) groundwater simulations and transfer learning (TL) to reduce data generation costs while maintaining model accuracy. As a computational example, we train a deep convolutional neural network (CNN) to reconstruct permeability fields from saturation maps derived from a multiphase flow model. Starting with very low- and low-fidelity data generated on increasingly coarse meshes, we pretrain the CNN, followed by output-layer training and fine-tuning using only a limited number of high-fidelity samples. We demonstrate the surrogate’s robustness when interpreting low-quality inputs—such as interpolated maps or data affected by noise—which has strong implications for the applicability in practical hydrogeological scenarios. This multilevel MF-TL strategy achieves a favorable trade-off between computational efficiency and predictive accuracy, significantly outperforming high-fidelity-only approaches under the same computational budget.

Chiofalo, Alessia [University of Bologna] (ORCID:0↗

DOE OSTI2021

Supersaturation, buoyancy, and deep convection dynamics

Motivated by recent discussions concerning differences of convective dynamics in polluted and pristine environments, the so-called convective invigoration in particular, this paper provides an analysis of factors affecting convective updraft buoyancy, such as the in-cloud supersaturation, condensate and precipitation loading, and entrainment. We use the deep convective period from simulations of daytime convection development over land discussed in our previous publications. An entraining parcel framework is used in the theoretical analysis. We show that for the specific case considered here, finite (positive) supersaturation noticeably reduces pseudo-adiabatic parcel buoyancy and cumulative convective available potential energy (cCAPE) in the lower troposphere. This comes from keeping a small fraction of the water vapor in a supersaturated state and thus reducing the latent heating. Such a lower-tropospheric impact is comparable to the effects of condensate loading and entrainment in the idealized parcel framework. For the entire tropospheric depth, loading and entrainment have a much more significant impact on the total CAPE. For the cloud model results, we compare ensemble simulations applying either a bulk microphysics scheme with saturation adjustment or a more comprehensive double-moment scheme with supersaturation prediction. We compare deep convective updraft velocities, buoyancies, and supersaturations from all ensembles. In agreement with the parcel analysis, the saturation-adjustment scheme provides noticeably stronger updrafts in the lower troposphere. For the simulations predicting supersaturation, there are small differences between pristine and polluted conditions below the freezing level that are difficult to explain by standard analysis of the in-cloud buoyancy components. By applying the piggybacking technique, we show that the lower-tropospheric buoyancy differences between pristine and polluted simulations come from a combination of temperature (i.e., latent heating) and condensate loading differences that work together to make polluted buoyancies and updraft velocities slightly larger when compared to their pristine analogues. Overall, the effects are rather small and contradict previous claims of a significant invigoration of deep convection in polluted environments.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Convection-Permitting Ensembles of an Isolated Mountain Thunderstorm during RELAMPAGO/CACTI

Abstract The north–south-oriented Sierras de Córdoba (SDC) ridge in central Argentina is noted for initiating thunderstorms that may grow into intense mesoscale convective systems (MCSs). It also initiates more isolated, shorter-lived cells under weaker synoptic forcing. These cells are less impactful than MCSs but may be difficult to predict in convective-scale numerical weather prediction (NWP) due to their strong sensitivities to subgrid and partially resolved processes. To study the mechanisms and predictability of such cells, convection-permitting ensemble simulations were conducted of an isolated, diurnally forced SDC thunderstorm during Cloud, Aerosol, and Complex Terrain Interactions (CACTI)/Remote Sensing of Electrification, Lightning, and Mesoscale/Microscale Processes with Adaptive Ground Observations (RELAMPAGO). The rich observational data facilitated detailed ensemble verification, where dry biases in the surface energy balance and soil moisture were identified. These biases promoted rapid removal of convective inhibition and an early onset of precipitating cells over the SDC that were shallower and weaker than the observed cell. Correction, and then overcorrection, of the soil moisture bias in two successive ensembles was required to rectify the surface energy balance and improve the representation of the SDC cell. Nevertheless, substantial ensemble variability in convective precipitation was found, with some members producing more widespread convection than observed and others producing no deep convection at all. This variability was largely explained by a combination of thermodynamic and dynamic mechanisms, dominated by a positive sensitivity of convective precipitation to preconvective moist instability over the ridge. Secondary sensitivities were found to low-level upward mass flux and midlevel cross-barrier winds, the latter of which caused gravity waves with elevated downdrafts that tended to suppress incipient clouds.

Lopez, Andres [Department of Atmospheric and Ocean↗