Search NASA⌕ Search

SEARCH · Search NASA

Results for “Separation of variables”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

47 Tuc in Rubin Data Preview 1. Exploring Early LSST Data and Science Potential

We present analyses of the early data from Rubin Observatory’s Data Preview 1 (DP1) for the field of the globular cluster 47 Tuc. The DP1 data set for 47 Tuc includes four nights of observations from the Rubin Commissioning Camera (LSSTComCam), covering multiple bands (ugriy). We address challenges of crowding in the inner region of the cluster and toward the SMC in DP1, and demonstrate improved star–galaxy separation by fitting fifth-degree polynomials to the stellar loci in color–color diagrams and applying multidimensional sigma clipping. We compile a catalog of 3576 probable 47 Tuc member stars selected via a combination of isochrone, Gaia proper-motion, and color–color space matched filtering. We explore the sources of photometric scatter in the 47 Tuc color–color sequence, evaluating contributions from various potential sources, including differential extinction within the cluster. Finally, of the 72 well-characterized variables in the field, we recover three known variable stars, including two RR Lyrae and one eclipsing binary, in the coadd-based object catalog, and identify 62 in the difference image-based object catalog. Although the DP1 lightcurves have sparse temporal sampling, they appear to follow the patterns of densely sampled literature lightcurves well. Despite some data limitations for crowded-field stellar analysis, DP1 demonstrates the promising scientific potential for future LSST data releases.

Choi, Yumi [NSF National Optical-Infrared Astronom↗
DOE OSTI2020

Markovian transport processes in a heterogeneous, variably saturated watershed: A multi-domain spatial Markov model

Transport processes in watersheds remain challenging to describe accurately with upscaled models because of their intricate complexities at multiple scales that can lead to multi-modal, non-Fickian breakthrough curves. One of the recent advances in solute transport modeling has been the Lagrangian spatial Markov model (SMM), which describes transport using velocity increments and models the correlated transitions between them. The method has been applied previously to saturated porous media and fractured media, but these successful applications suggest it may be useful for describing watershed scale processes as well. An existing, 3-D, heterogeneous variably saturated model of a headwaters catchment and classical random walk particle tracking were used to generate a forward simulation and the resulting transport was analyzed with a spatial Markov framework. Existing SMM methods were used to analyze the transport paths through the complex, unconfined domain and we found that the vadose zone plays an important role in the correlation structure. As such, a new kind of SMM is proposed, termed a multi-domain SMM (MD-SMM), that uses separate correlation models for auto-transitions in the vadose zone and the saturated zone, and for the cross-transitions between them. The MD-SMM improves representation of transport relative to the single domain model and can naturally delineate times when a particle is in the vadose zone or the saturated zone, adding an additional degree of realism to the upscaled model. Lastly, we show that the transition probability matrix for transport in the saturated portion of this unconfined system is approximately Markovian, which is an important validation of the proposed framework.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Confined and deconfined chaos in classical spin systems

Weakly perturbed integrable many-body systems are typically chaotic, and thermal at late times. However, there are distinct relationships between the timescales for thermalization and chaos. The typical relationship is confined chaos: when trajectories are still confined to regions in phase space with constant conserved quantities (actions), the conjugate angle variables are already unstable. Chaotic instabilities thus far precede thermalization. In a different relationship, which we term deconfined chaos, chaotic instabilities and thermalization occur on the same timescale. We investigate these two qualitatively distinct scenarios through numerical and analytical studies of two perturbed integrable classical spin models: the Ishimori spin chain (confined chaos), and the central spin model with XX interactions (deconfined chaos). We analytically establish (super)-integrability in the latter model in a microcanonical shell. Deconfined chaos emerges through the separation of phase space into large quasi-integrable regions and a thin chaotic manifold. The latter leads to chaos and thermalization on the fastest possible timescale, which is proportional to the inverse perturbation strength. This behavior is reminiscent of the quantum SYK models and strange metals.

Chaotic Dynamics (nlin.CD)↗
DOE OSTI2021

Economic and Environmental Assessment of Biological Conversions of Agile BioFoundry (ABF) Bio-Derived Chemicals

Bio-derived chemicals are an essential part of the growing bioeconomy. They possess the potential to boost a new domestic bioproduct industry, improve the sustainability of integrated biomanufacturing, and reduce U.S. dependence on fossil energy. The Agile BioFoundry (ABF) consortium, a Department of Energy (DOE)-sponsored collaboration, is investigating biobased pathways to produce advantaged products including advanced biofuels, fuel intermediates, and bioproducts. ABF is integrating advanced computational tools for biological engineering, process design and data analysis, and economic/sustainability modeling tools into a comprehensive and dynamic platform for biomanufacturing of microbes and using them to produce key metabolic intermediates or beachhead molecules which may be derivatized into several distinct bioproducts of industrial interest. In this presentation, we first discuss a methodology to select a single exemplar product molecule to represent each beachhead pathway based on similarities with other end-molecule options (yields/titers/rates, fermentation operation mode, oxygen requirements, general separations challenges, etc.), in order to maintain a reasonable number of cases for rigorous process simulation. We then use techno-economic analysis (TEA) and life-cycle analysis (LCA) to highlight sensitivities and trends around economic performance and environmental footprint, reflecting examples for two selected ABF technology pathways to bio-derived chemicals: 1) adipic acid production via muconic acid fermentation from mixed sugars with Pseudomonas putida and 2) cineole via geranyl diphosphate with Rhodosporidium toruloides. We present multi-variable scan plots highlighting key drivers on costs and greenhouse gas emissions in order to identify the parameter space in which each pathway could ultimately achieve economic and environmental sustainability meeting or exceeding commodity product benchmarks, thus prioritizing future R&D focus areas.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Oxygen Storage Incorporated Into Net Power and the Allam–Fetvedt Oxy-Fuel sCO2 Power Cycle—Techno-Economic Analysis

Abstract With the planned future reliance on variable renewable energy, the ability to store energy for prolonged time periods will be required to reduce the disruption of market fluctuations. This paper presents a method to analyze a hybrid liquid-oxygen (LOx) storage/direct-fired supercritical carbon dioxide (sCO2) power cycle and optimize the economic performance over a diverse range of scenarios. The system utilizes a modified version of the NET Power process to produce energy when energy demand exceeds the supply while displacing much of the cost of the air separation unit (ASU) energy requirements through cryogenic storage of oxygen. The model uses marginal cost of energy data to determine the optimal times to charge and discharge the system over a given scenario. The model then applies ramp rates and other time-dependent factors to generate an economic model for the system without storage considerations. The size of the storage system is then applied to create a realistic model of the plant operation. From the real plant operation model, the amount of energy charged and discharged, the capital expenditures (CAPEX) of each system, energy costs and revenue and other parameters can be calculated. The economic parameters are then combined to calculate the net present value (NPV) of the system for the given scenario. The model was then run through the SMPSO genetic algorithm in Python for a variety of geographic regions and large-scale scenarios (high solar penetration) to maximize the NPV based on multiple parameters for each subsystem. The LOx storage requirements will also be discussed.

Engineering↗
DOE OSTI2026

WINDPROF: Merged Best-Estimate Wind Profile Data – Cape Cod (WFIP3 Campaign)

WINDPROF provides 10-minute wind and turbulence profiles that integrate a scanning Doppler lidar, a profiling Doppler lidar, and two sonic anemometers across Northeast U.S. coastal and offshore sites during the WFIP3 campaign. Variables include wind speed, wind direction, vertical velocity, turbulence intensity, and turbulent kinetic energy, each with per-instrument quality control and inter-instrument agreement validation. Profiles are mapped to a standardized height grid – dedicated near-surface levels at the sonic measurement heights (4 m and 10 m AGL), 20 m spacing to 100 m, and 30 m spacing above – and carry component and derived uncertainty estimates. Heights are reported above ground level; the site's ground elevation is stored separately. The Cape Cod dataset covers 1 February 2024 – 6 September 2025.

17 WIND ENERGY↗
DOE OSTI2026

Fine-Grained Power and Energy Attribution on AMD GPU/APU-Based Exascale Nodes

Modern exascale GPU- and APU-based systems provide multiple power and energy sensors, but differences in scope, update rate, timing, and filtering complicate the attribution of short-lived accelerator activity. This paper presents a methodology to characterize and correct these effects on Cray EX systems with AMD Instinct MI250X GPUs (Frontier) and MI300A APUs (Portage). Using controlled square-wave workloads, we quantify update intervals, delay, aliasing, and variability across up to 512 GPUs and 480 APUs with on-chip (rocm-smi/amd-smi) and off-chip Cray Power Management sensors. We reconstruct power from cumulative energy counters to achieve faster response times, validate it against on-chip, off-chip, and node-level sensors, and integrate the resulting streams into a Score-P/PAPI-based tool for time-aligned, phase-level attribution. Applied to rocHPL, rocHPL-MxP, and HPG-MxP, the method separates energy savings due to reduced runtime from changes in power. Mixed precision reduces node energy on Frontier by 79% for rocHPL-MxP and 31% for HPG-MxP, with similar trends on Portage. These results provide portable guidance for sensor validation and power-aware optimization on current and future exascale systems.

Mcdaniel, Adam [ORNL] (ORCID:000000016926028X)↗
DOE OSTI2023

Swelling and permeability effects during propellant cookoff

Large rocket motors may violently explode when exposed to accidental fires. Even hot metal fragments from a nearby accident may penetrate the propellant and ultimately cause thermal ignition. A mechanistic understanding of heated propellants leading to thermal runaway is a major unsolved problem. Here we show that thermal ignition in propellants can be predicted using a universal cookoff model coupled to a micromechanics pressurization model. Our model predicts the time to thermal ignition in cookoff experiments with variable headspace volumes. Furthermore, we found that experiments with headspace volumes are more prone to deformation which distorts pores and causes increased permeability when the propellant expands into this headspace. Delayed ignition with larger headspace volume correlates with lower headspace pressures during decomposition. We found that our predictions matched experimental measurements best when the initial propellant was impermeable to gas flow rather than being permeable. Similar behavior is expected with other energetic materials with rubbery binders. Our model is validated using data from a separate laboratory. We also present an uncertainty analysis using Latin Hypercube Sampling (LHS) of thermal ignition caused by a steel fragment embedded in the propellant.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Climatology of the Vertical Profiles of Polarimetric Radar Variables and Retrieved Microphysical Parameters in Continental/Tropical MCSs and Landfalling Hurricanes

Most existing cloud models tend to overestimate the size of cloud ice particles and underestimate their concentration. This emphasizes the need to provide a reliable observational reference to optimize cloud model performance, particularly in areas of high concentration of ice at high altitudes. The dual-polarization radars give the community a unique opportunity to quantify cloud ice with a good accuracy using polarimetric radar retrievals. Here, in this study, we utilize the network of operational WSR-88D radars to build a climatology of the vertical profiles of radar variables, such as radar reflectivity Z, differential reflectivity Z DR , and specific differential phase K DP as well as the radar-retrieved vertical profiles of ice water content (IWC) above the melting layer and liquid water content below it, mean volume diameter D m , and total number concentration N t of ice and liquid particles. Such climatology was created for continental/marine mesoscale convective systems (MCSs) and tropical cyclones including hurricanes. The dataset includes 13 continental MCSs, 10 marine MCSs, and 11 tropical cyclones. Separate statistics of the “background” vertical profiles and the ones associated with high IWC aloft have been obtained in the course of this study. It is shown that continental MCSs exhibit larger size of ice in lower concentration aloft compared to the marine MCSs and especially tropical cyclones/hurricanes. A combination of high KDP and low Z aloft signifies lower D m , higher N t , and often substantial IWC.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Temporal Study 2021-2022: Sensor-Based Time Series of Surface Water Temperature, Specific Conductance, Total Dissolved Solids, Turbidity, pH, and Dissolved Oxygen from across Multiple Watersheds in the Yakima River Basin in Washington, USA (v2)

This dataset supports a broader study examining the drivers of temporal variability in sediment respiration rates in the Yakima River Basin. The dataset provides periodic (weekly or biweekly) in situ hydrological and water chemistry sensor data, handheld sensor water chemistry data, general environmental context photos, and field metadata collected at six sites within multiple rivers in the Yakima River Basin in Washington, USA. In addition to the sensor data, there are plots of continuous in situ sensor data and R scripts used to generate the plots. Related sample-based water chemistry data are published separately at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1898912.The data package was originally published in September 2022. It was updated in June 2025 (v2; modified files). See the change history section in the readme for more details. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. This dataset is comprised of one main data folder containing two sensor-specific subfolders, and one photographs folder. The main data folder includes file-level metadata (flmd), data dictionary (dd), installation methods, field metadata, handheld sensor data, field data collection protocols, international generic sample number (IGSN) mapping file, and a readme file. Each sensor subfolder (BarotrollAtm and MantaRiverData) contains a subfolder containing sensor timeseries data and plots. The BarotrollAtm Data subfolder contains In Situ Rugged BaroTROLL sensor pressure and air temperature data. The MantaRiverData subfolder contains Eureka Manta+ 35B multisonde temperature, specific conductance, and turbidity. The FieldPhotos folder contains environmental context photographs and videos. All files are .csv, .pdf, .R, .jpg, .jpeg, .heic, .mov, or .mp4.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Applying Particle Swarm Optimization and Extended Kalman Filtering to Model Kaplan Generation Dynamics for Hydropower Systems

Variable renewable generation is increasing the need for hydropower plants to provide fast and flexible grid support, which places new demands on plant-level dynamic models used for monitoring, control, and operational decision-making. This need is especially important for hydroelectric systems, where turbine and generator dynamics are strongly coupled, nonlinear, and time-varying, making accurate real-time representation difficult. To address this problem, this paper develops a digital twin (DT) framework for a synchronous generator–Kaplan turbine system using an explicit separation of slow turbine dynamics and fast generator dynamics. The turbine subsystem is represented by a six-coefficient model, whose parameters are identified offline using particle swarm optimization, while the generator subsystem is updated online through an extended Kalman filter for real-time state and parameter estimation. These models are integrated within a closed-loop simulation that includes a proportional–integral–derivative–double-derivative governor and excitation system, allowing the DT to track plant behavior under realistic operating conditions. Unlike prior studies that treat turbine and generator modeling separately or rely mainly on simulated inputs, the proposed framework is validated using real operational data from a hydropower plant. Results show that the DT reproduces terminal voltage, active power, and reactive power with a normalized root mean square error of approximately 5%. This hybrid offline–online formulation constitutes the main contribution of the work, providing an adaptive and practically deployable DT for hydropower systems with direct relevance to control improvement, performance monitoring, and grid-support applications under high renewable penetration.

13 HYDRO ENERGY↗
DOE OSTI2021

Parallel physics-informed neural networks via domain decomposition

Here we develop a distributed framework for the physics-informed neural networks (PINNs) based on two recent extensions, namely conservative PINNs (cPINNs) and extended PINNs (XPINNs), which employ domain decomposition in space and in time-space, respectively. This domain decomposition endows cPINNs and XPINNs with several advantages over the vanilla PINNs, such as parallelization capacity, large representation capacity, efficient hyperparameter tuning, and is particularly effective for multi-scale and multi-physics problems. Here, we present a parallel algorithm for cPINNs and XPINNs constructed with a hybrid programming model described by MPI + X, where X ∈ {CPUs, GPUs}. The main advantage of cPINN and XPINN over the more classical data and model parallel approaches is the flexibility of optimizing all hyperparameters of each neural network separately in each subdomain. We compare the performance of distributed cPINNs and XPINNs for various forward problems, using both weak and strong scalings. Our results indicate that for space domain decomposition, cPINNs are more efficient in terms of communication cost but XPINNs provide greater flexibility as they can also handle time-domain decomposition for any differential equations, and can deal with any arbitrarily shaped complex subdomains. To this end, we also present an application of the parallel XPINN method for solving an inverse diffusion problem with variable conductivity on the United States map, using ten regions as subdomains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

DEPRECATED AI-Batt-OS (Autonomous Identification of Battery Life Models - Open Source) [SWR 21-17]

DEPRECATED. This repository was archived by the owner on Jun 30, 2026. It is now read-only. Open source implementation of some of the methods utilized by AI-Batt, a battery lifetime modeling and analysis toolkit provided by the National Laboratory of the Rockies (NLR). This software demonstrates the use of bi-level optimization and symbolic regression techniques to semi-autonomously identify algebraic models predicting the capacity fade of lithium-ion batteries during calendar aging. Modeling the degradation of batteries is a complex task, due to the difficulty in separating the time-dependent and time-independent factors impacting cell level degradation, across multiple data series with different numbers of measurements and/or data quality. Bi-level optimization enables model parameters to be optimized to either the entire data set or to individual data series, allowing statistical disambiguation of global behaviors (data series independent) and local behaviors (data series dependent). Symbolic regression is used to automatically search for optimal low-dimesional models predicting the variation of locally optimized parameters versus time-independent experimental variables from millions of possible models, resulting in a more accurate and repeatable model identification process than is possible by a manual search. The provided tools also implement cross-validation and bootstrap resampling schemes, empowering statistical model comparison/selection and quantification of model uncertainties. An example script replicates the results from the manuscript "Challenging Practices of Algebraic Battery Life Models through Statistical Validation and Model Identification via Machine-Learning", submitted to ECS. All code is written in MATLAB. Requires the Statistics and Machine Learning Toolbox. Contact Dr. Paul Gasper at Paul.Gasper@nlr.gov for any questions.

Gasper, Paul↗
DOE OSTI2020

Analysis of Air-Purifying Respirator (APR) and Powered Air-Purifying Respirator (PAPR) Cartridge Performance Testing on a Hanford AX Tank Farm Exhauster Slipstream Volume 1

Washington River Protection Solutions (WRPS) conducted tests of four types of chemical cartridges for air-purifying respirators (APR) and powered air-purifying respirators (PAPR) to determine the period of time the cartridges would provide adequate performance1 for APRs and PAPRs used to protect workers when exposed to a mixture of Chemicals of Potential Concern (COPCs) from vapors exiting the Hanford AX tank farm exhauster slipstream. The Occupational Safety and Health Administration (OSHA) considers cartridge testing to be a valid approach for establishing a cartridge service life. Testing is applied in situations where mixtures of COPCs exist, and where other approaches, such as manufacturer recommendations and modeling, are less reliable. The tests were designed and conducted to assure measurement and/or control of the key variables OSHA identified as important to estimate the cartridge service life, including temperature, humidity, COPC concentration, breathing rate, and cartridge adsorption capacity. Cartridge testing using vapors from a Hanford AX tank farm exhauster slipstream was conducted from August 25–27, 2017. Vapors from the exhauster slipstream were fed to two respirator cartridge test stands developed by WRPS in collaboration with HiLine Engineering (Richland, Washington). Four different cartridges were assessed. Multipurpose APR cartridges—SCOTT 7422-SD1 and SCOTT 7422-SC1 (SCOTT Safety, Monroe, North Carolina)—were assessed on separate days using an APR cartridge test stand. Multipurpose PAPR cartridges—MSA-TL (TL1) (MSA Safety Inc., Pittsburgh, Pennsylvania) and 3M FR57 (TL2) (3M Company, Maplewood, Minnesota)—also were tested over the same two days using a separate PAPR cartridge test stand. Sample media (i.e., sorbent tubes) were used to collect samples of the vapor stream entering and exiting the respirator cartridges and were subsequently analyzed for COPC concentrations. Pacific Northwest National Laboratory was tasked with conducting an independent analysis of the analytical results and making recommendations based on the results for respiratory cartridge performance and service life. The key conclusions from the analysis are described below.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Post-fire soil emissions of nitric oxide (NO) and nitrous oxide (N 2 O) across global ecosystems: a review

Wildfires may increase soil emissions of trace nitrogen (N) gases like nitric oxide (NO) and nitrous oxide (N 2 O) by changing soil physicochemical conditions and altering microbial processes like nitrification and denitrification. When 34 studies were synthesized, we found a significant increase in both NO and N 2 O emissions up to 1 year post-fire across studies spanning ecosystems globally. However, when fluxes were separated by ecosystem type, we found that individual ecosystem types responded uniquely to fire. Forest soils tended to emit more N 2 O after fire, but there was no significant effect on NO. Shrubland soils showed significant increases in both NO and N 2 O emissions after fires; often with extremely large but short-lived NO pulses occurring immediately after fire. Grassland NO emissions increased after fire, but the size of this effect was small relative to shrublands. N 2 O emissions from burned grasslands were highly variable with no significant effect. To better understand the variation in responses to fire across global ecosystems, more consistent measurements of variables recognized as important controls on soil fluxes of NO and N 2 O (e.g., N cycling rates, soil water content, pH, and substrate availability) are needed across studies. We also suggest that fire-specific elements like burn severity, microbial community succession, and the presence of char be considered by future studies. Our synthesis suggests that fires can exacerbate ecosystem N loss long after they burn, increasing soil emissions of NO and N 2 O with implications for ecosystem N loss, climate, and regional air quality as wildfires increase globally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Investigation of a high-temperature combination heat pump for lower-cost electrification in multifamily buildings

The development of space and water heating combination heat pumps capable of generating water temperatures high enough for convective heat emitters will enable more cost-effective and equitable decarbonization solutions for electrifying multifamily buildings. Here, in this paper, multifamily building models and a charge-sensitive mechanistic cycle model of a combination heat pump are developed, and the system performance is predicted based on the models. Unlike other state-of-the-art residential heat pumping equipment, the modeled combination heat pump using an economized, fluid-injected variable-speed compressor can achieve higher temperature lifts of 40° - 85°C, with lower installation costs and complexity. The model predicted heating coefficient of performance (COP h ) is 2.1 at an ambient temperature of -15°C with a high-temperature lift of nearly 85°C, and a seasonal coefficient of performance in heating mode (SCOP h ) ranges from 2 - 4 for different locations. The system shows 30% - 90% lower CO 2 eq emissions over a condensing gas boiler and 9% - 13% lower projected installation costs than two separate space and water heat pumping appliances.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

High-resolution (30-m) urban land cover projections for Los Angeles California Urban Area: 2010 to 2100 under SSP5

These data represent simulations of future land use and land cover (LULCC) for Los Angeles urban area (U.S. Census Bureau defined area) as raster tiff images at a 30-m pixel resolution and at decadal time steps from 2010 to 2100. LULCC classes in this product follow the National Land Cover Dataset (NLCD) classification. NLCD 21-24 correspond to open developed, low developed, medium developed, and high developed urban land classes, respectively. Only urban land cover classes (NLCD class 21, 22, 23, and 24) are dynamic over time; however, all NLCD classes are included in the final product. Therefore, NLCD classes that do not convert to an urban class will be similar to year 2000. The products were developed using a hybridized statistical and cellular automata approach. Linear mixed models (LMMs) were used to estimate future urban land budgets based on 1-km urban land fraction projections from Gao and Pesaresi (2021), whereas separate generalized linear mixed models (GLMMs) were used to estimate shifts in urban land intensities based on retrospective shifts in NLCD urban class intensities over a 20- year period. Based on urban land allocations from the statistical models, a cellular-automata and downscaling routine was used to simulate dynamic urban land expansion at a 30-m resolution based on suitability criteria. Scenarios of future urban landcover change projections include variant solutions for the Shared Socioeconomic Pathway 5 (SSP5) based on different population assumptions, different land use intensification assumptions, variable land zoning constraints, and iterative adjustments to correct for over allocation of urban expansion across decadal time periods from 2010 to 2100. This results in 320 raster products.

Land↗
DOE OSTI2024

Evaluation of Autoconversion Representation in E3SMv2 Using an Ensemble of Large-Eddy Simulations of Low-Level Warm Clouds

In numerical atmospheric models that treat cloud and rain droplet populations as separate condensate categories, precipitation initiation in warm clouds is often represented by an autoconversion rate (Au), which is the rate of formation of new rain droplets through the collisions of cloud droplets. Being a function of the cloud droplet size distribution (DSD), the local Au is commonly parameterized as a function of DSD moments: cloud droplet number (n c ) and mass (q c ) concentrations. When applied in a large-scale model, the grid-mean Au must also include a correction, or enhancement factor, to account for the horizontal variability of the cloud properties across the model grid. In this study, we evaluate the Au representation in the Energy Exascale Earth System Model version 2 (E3SMv2) climate model using large-eddy simulations (LES), which explicitly resolve cloud droplet spectra, and therefore the local Au, as well as its spatial variability. The analysis of an ensemble of warm low-level cloud cases shows that the E3SMv2 formulation represents the Au reasonably well compared to the horizontally averaged explicitly computed rate from LES. The agreement, however, comes from a combination of an underestimated E3SM-tuned local Au rate and an overestimated subgrid cloud variability enhancement factor. The latter bias is traced to neglecting the horizontal variability of n c and its co-variability with q c in parameterizing the grid-mean Au.

54 ENVIRONMENTAL SCIENCES↗