Search NASA⌕ Search

SEARCH · Search NASA

Results for “ensemble simulation”

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 541 records · Page 30

PNNL-ANL Hydrometeorological Super Ensemble

The current dataset contains data upload links to the following **hydrologic (water balance), river routing (water management), and hydropower simulation** data over CONUS: * Climate Forcing: **Livneh** (https://www.nature.com/articles/sdata201542) * Simulation Scenario: **Historical** * Simulation Period: **1971-2013 (1972-2013 for Hydropower)** * Simulation Models: **VIC (Variable Infiltration Capacity)**, **mosartwmpy (Model for Scale Adaptive River Transport-Water Management in Python)**, and **PNNL B1Hydro** * Output Format: **NetCDF** and **CSV**

Tidwell, Vincent C [Pacific Northwest National Lab↗

PNNL-ANL Hydrometeorological Super Ensemble

The current dataset contains data upload links to the following hydrologic (water balance), river routing (water management), and hydropower simulation data over CONUS: Climate Forcing: ClimRR (https://climrr.anl.gov/climrrdata) Simulation Scenario: Historical, Mid-Century, End-Century Simulation Period: 1995-2004, 2045-2054, 2085-2094 Simulation Models: VIC (Variable Infiltration Capacity), mosartwmpy (Model for Scale Adaptive River Transport-Water Management in Python), and PNNL B1Hydro Output Format: NetCDF and CSV

Tidwell, Vincent C [Pacific Northwest National Lab↗

Broken Symmetry and Coherent Structure in MHD Turbulence

Absolute equilibrium ensemble theory for ideal homogeneous magnetohydrodynamic (MHD) turbulence is fairly well developed. Theory and Simulation indicate that ideal MHD turbulence non-ergodic and contains coherent structure. The question of applicability real (i.e., dissipative) MHD turbulence is examined. Results from several very long time numerical simulations on a 64(exp 3) grid are presented. It is seen that coherent structure begins to form before decay dominates over nonlinearity. The connection with inverse spectral cascades and selective decay will also be discussed.

Shebalin, John V.↗

Strong Regional Influence of Climatic Forcing Datasets on Global Crop Model Ensembles

We present results from the Agricultural Model Intercomparison and Improvement Project (AgMIP) Global Gridded Crop Model Intercomparison (GGCMI) Phase I, which aligned 14 global gridded crop models (GGCMs) and 11 climatic forcing datasets (CFDs) in order to understand how the selection of climate data affects simulated historical crop productivity of maize, wheat, rice and soybean. Results show that CFDs demonstrate mean biases and differences in the probability of extreme events, with larger uncertainty around extreme precipitation and in regions where observational data for climate and crop systems are scarce. Countries where simulations correlate highly with reported FAO national production anomalies tend to have high correlations across most CFDs, whose influence we isolate using multi-GGCM ensembles for each CFD. Correlations compare favorably with the climate signal detected in other studies, although production in many countries is not primarily climate-limited (particularly for rice). Bias-adjusted CFDs most often were among the highest model-observation correlations, although all CFDs produced the highest correlation in at least one top-producing country. Analysis of larger multi-CFD-multi-GGCM ensembles (up to 91 members) shows benefits over the use of smaller subset of models in some regions and farming systems, although bigger is not always better. Our analysis suggests that global assessments should prioritize ensembles based on multiple crop models over multiple CFDs as long as a top-performing CFD is utilized for the focus region.

Agricultural Model Intercomparison and Improvement↗

Performance of Trajectory Models with Wind Uncertainty

Typical aircraft trajectory predictors use wind forecasts but do not account for the forecast uncertainty. A method for generating estimates of wind prediction uncertainty is described and its effect on aircraft trajectory prediction uncertainty is investigated. The procedure for estimating the wind prediction uncertainty relies uses a time-lagged ensemble of weather model forecasts from the hourly updated Rapid Update Cycle (RUC) weather prediction system. Forecast uncertainty is estimated using measures of the spread amongst various RUC time-lagged ensemble forecasts. This proof of concept study illustrates the estimated uncertainty and the actual wind errors, and documents the validity of the assumed ensemble-forecast accuracy relationship. Aircraft trajectory predictions are made using RUC winds with provision for the estimated uncertainty. Results for a set of simulated flights indicate this simple approach effectively translates the wind uncertainty estimate into an aircraft trajectory uncertainty. A key strength of the method is the ability to relate uncertainty to specific weather phenomena (contained in the various ensemble members) allowing identification of regional variations in uncertainty.

Lee, Alan G.↗

SST Variation Due to Interactive Convective-Radiative Processes

The recent linking of Cloud-Resolving Models (CRMs) to Ocean-Mixed Layer (OML) models has provided a powerful new means of quantifying the role of cloud systems in ocean-atmosphere coupling. This is due to the fact that the CRM can better resolve clouds and cloud systems and allow for explicit cloud-radiation interaction. For example, Anderson (1997) applied an atmospheric forcing associated with a CRM simulated squall line to a 3-D OML model (one way or passive interaction). His results suggested that the spatial variability resulting from the squall forcing can last at least 24 hours when forced with otherwise spatially uniform fluxes. In addition, the sea surface salinity (SSS) variability continuously decreased following the forcing, while some of the SST variability remained when a diurnal mixed layer capped off the surface structure. The forcing used in the OML model, however, focused on shorter time (8 h) and smaller spatial scales (100-120 km). In this study, the 3-D Goddard Cumulus Ensemble Model (GCE; 512 x 512 x 23 cu km, 2-km horizontal resolution) is used to simulate convective active episodes occurring in the Western Pacific warm pool and Eastern Atlantic regions. The model is integrated for seven days, and the simulated results are coupled to an OML model to better understand the impact of precipitation and changes in the planetary boundary layer upon SST variation. We will specifically examine and compare the results of linking the OML model with various spatially-averaged outputs from GCE simulations (i.e., 2 km vs. 10-50 km horizontal resolutions), in order to help understand the SST sensitivity to multi-scale influences. This will allow us to assess the importance of explicitly simulated deep and shallow clouds, as well as the subgrid-scale effects (in coarse-model runs) upon SST variation. Results using both 1-D and 2-D OML models will be evaluated to assess the effects of horizontal advection.

Tao, W.-K.↗

Uncertainty-informed selection of CMIP6 Earth System Model subsets for use in multisectoral and impact models

Earth system models (ESMs) and general circulation models (GCMs) are heavily used to provide inputs to sectoral impact and multisector dynamic models, which include representations of energy, water, land, economics, and their interactions. Therefore, representing the full range of model uncertainty, scenario uncertainty, and interannual variability that ensembles of these models capture is critical to the exploration of the future co-evolution of the integrated human–Earth system. The pre-eminent source of these ensembles has been the Coupled Model Intercomparison Project (CMIP). With more modeling centers participating in each new CMIP phase, the size of the model archive is rapidly increasing, which can be intractable for impact modelers to effectively utilize due to computational constraints and the challenges of analyzing large datasets. In this work, we present a method to select a subset of the latest phase, CMIP6, featuring models for use as inputs to a sectoral impact or multisector dynamics models, while prioritizing preservation of the range of model uncertainty, scenario uncertainty, and interannual variability in the full CMIP6 ensemble results. This method is intended to help impact modelers select climate information from the CMIP archive efficiently for use in downstream models that require global coverage of climate information. This is particularly critical for large-ensemble experiments of multisector dynamic models that may be varying additional features beyond climate inputs in a factorial design, thus putting constraints on the number of climate simulations that can be used. We focus on temperature and precipitation outputs of CMIP6 models, as these are two of the most used variables among impact models, and many other key input variables for impacts are at least correlated with one or both of temperature and precipitation (e.g., relative humidity). Besides preserving the multi-model ensemble variance characteristics, we prioritize selecting CMIP6 models in the subset that preserve the very likely distribution of equilibrium climate sensitivity values as assessed by the latest Intergovernmental Panel on Climate Change (IPCC) report. This approach could be applied to other output variables of climate models and, possibly when combined with emulators, offers a flexible framework for designing more efficient experiments on human-relevant climate impacts. It can also provide greater insight into the properties of existing CMIP6 models.

Snyder, Abigail C.↗

One‐at‐a‐Time Parameter Perturbation Ensemble of the Community Land Model, Version 5.1

Comprehensive land models are subject to significant parametric uncertainty, which can be hard to quantify due to the large number of parameters and high model computational costs. We constructed a large parameter perturbation ensemble (PPE) for the Community Land Model version 5.1 with biogeochemistry configuration (CLM5.1-BGC). We performed more than 2,000 simulations perturbing 211 parameters across six forcing scenarios. This provides an expansive data set, which can be used to identify the most influential parameters on a wide range of output variables globally, by biome, or by plant functional type. We found that parameter effects can exceed scenario effects and that a small number of parameters explains a large fraction of variance across our ensemble. The most important parameters can differ regionally and also based on the forcing scenario. The software infrastructure developed for this experiment has greatly reduced the human and computer time needed for CLM PPEs, which can facilitate routine investigation of parameter sensitivity and uncertainty, as well as automated calibration.

Kennedy, Daniel [NSF National Center for Atmospher↗

Deep learning model for fast, science-based forecasting of fluid migration along faults in geologic carbon storage scenarios

Effective long-term geologic storage depends on robust site selection and credible, science-based forecasting of subsurface behavior to ensure storage integrity. For this work, we develop a deep learning–based reduced-order model (ROM) to quantify potential carbon dioxide (CO₂) and brine migration through geological faults. The ROM combines a Transformer model for binary classification and a Stacked Ensemble for regression, trained on a comprehensive dataset generated from 1400 physics-based reservoir simulations. Key geologic and operational parameters—including fault geometry, reservoir structure, and injection conditions—were systematically varied to capture a wide range of fluid migration scenarios. The ROM accurately predicts the onset of migration, cumulative migration volumes of both CO₂ and brine, and associated migration rates, as compared to an independent set of validation simulations, while significantly reducing computational cost compared to traditional simulation methods. Model performance was evaluated across diverse fault configurations, revealing that shallow reservoir geometry and fault angle are among the most influential factors governing migration behavior. Sensitivity analysis using SHapley Additive exPlanations (SHAP) provided interpretability, revealing distinct patterns in how geological and operational features drive transient versus cumulative migration outcomes. The ROM’s ability to rapidly simulate fault migration scenarios enables efficient sensitivity analyses, scenario evaluations, and decision support for site selection and monitoring design. This approach enhances the safety, scalability, and long-term operational performance of geologic carbon storage (GCS) systems by providing a robust, interpretable tool for predicting subsurface fluid migration and assessing fault-related migration potential.

42 ENGINEERING↗

Predicting permeation of compounds across the outer membrane of P. aeruginosa using molecular descriptors

The ability Gram-negative pathogens have at adapting and protecting themselves against antibiotics has increasingly become a public health threat. Data-driven models identifying molecular properties that correlate with outer membrane (OM) permeation and growth inhibition while avoiding efflux could guide the discovery of novel classes of antibiotics. Here we evaluate 174 molecular descriptors in 1260 antimicrobial compounds and study their correlations with antibacterial activity in Gram-negative Pseudomonas aeruginosa. The descriptors are derived from traditional approaches quantifying the compounds’ intrinsic physicochemical properties, together with, bacterium-specific from ensemble docking of compounds targeting specific MexB binding pockets, and all-atom molecular dynamics simulations in different subregions of the OM model. Using these descriptors and the measured inhibitory concentrations, we design a statistical protocol to identify predictors of OM permeation/inhibition. We find consistent rules across most of our data highlighting the role of the interaction between the compounds and the OM. An implementation of the rules uncovered in our study is shown, and it demonstrates the accuracy of our approach in a set of previously unseen compounds. Our analysis sheds new light on the key properties drug candidates need to effectively permeate/inhibit P. aeruginosa, and opens the gate to similar data-driven studies in other Gram-negative pathogens.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Zero-gravity aerosol behavior

The feasibility and scientific benefits of a zero gravity aerosol study in an orbiting laboratory were examined. A macroscopic model was devised to deal with the simultaneous effects of diffusion and coagulation of particles in the confined aerosol. An analytical solution was found by treating the particle coagulation and diffusion constants as ensemble parameters and employing a transformation of variables. The solution was used to carry out simulated zero gravity aerosol decay experiments in a compact cylindrical chamber. The results demonstrate that the limitations of physical space and time imposed by the orbital situation are not prohibitive in terms of observing the history of an aerosol confined under zero gravity conditions. While the absence of convective effects would be a definite benefit for the experiment, the mathematical complexity of the problem is not greatly reduced when the gravitational term drops out of the equation. Since the model does not deal directly with the evolution of the particle size distribution, it may be desirable to develop more detailed models before undertaking an orbital experiment.

Edwards, H. W.↗

A fundamental model and efficient inference for SAR ocean imagery

Employing a synthetic aperture radar (SAR) imaging model based on fundamental models of nonlinear hydrodynamics, electromagnetic scattering from a two-scale surface, and SAR imaging of a time-variant scene, the optimal (minimum mean-square error) estimates of the parameters of a sinusoidal, long gravity wave, and the short gravity wave ensemble are found in an efficient recursive form and their performance evaluated, generally by numerical simulation, in a one-dimensional stationary version. An application is made to Seasat-SAR complex imagery.

Harger, R. O.↗

Triplets of galaxies: Their dynamics, evolution, and the origin of chaos in them

Recently Karachentsev's group at The Smithsonian Astrophysical Observatory (SAO) (6-meter Telescope Observatory) published a list of 84 triple systems of galaxies with their distances, radial (line of sight) velocities, and angular sizes (Karachentseva et al., 1988). This gives a new ground for studies of the dark matter problem which fills the gap between the large cosmic scales (White, 1987; Dekel and Rees, 1987, and Einasto et al., 1977) and the scale of individual galaxies (Erickson et al., 1987). The data on the typical velocity dispersions and linear dimension of the triplets indicate that they contain considerable amounts of dark matter (see also earlier work of Karachentseva et al., 1979). Numerical simulations show that the statistical characteristics of the Karachentsev triplets can be imitated by model ensembles of triple systems with dark matter masses M sub d = (1-3 x 10(exp 12) M sub O, which is almost ten times greater than the typical mass of stellar galaxies estimated by the standard mass-to-luminosity ration (Kiseleva and Chernin, 1988). Here, the authors report that important information can be drawn from the data on the visible configurations of these systems. The statistics of configurations provide an independent evidence for dark matter in the triplets; moreover, it enables one to argue that dark matter seems to be distributed over the whole volume of the typical triplet forming its common corona rather than concentrated within individual coronae (or haloes) of the member galaxies.

Chernin, Arthur D.↗

The Variability and Predictability of Precipitation over Southeast Asian

The large scale atmospheric moisture budget, including precipitation, moisture flux convergence, and evaporation is described over southeast Asian. ne results are from the Goddard Earth Observing System (GEOS-1) data assimilation system for the period March 1981 - February 1995. On average, in the late spring and early summer season, precipitation is largely balanced by the atmospheric moisture flux convergence over southeast Asian; evaporation is also an important contributor. The seasonal variability of precipitation is best related to variations in the moisture flux convergence. On the monthly time scales, both the precipitable water and precipitation anomalies are most closely related to moisture flux convergence anomalies and less closely related to evaporation anomalies. The impact of the low-level jet (LLJ) on the strong low-level moisture flux convergence will be discussed. An ensemble of GEOS-2 atmospheric General Circulation Model (GCM) seasonal forecasts and long-term simulation are analyzed to assess the controlling influences of boundary forcing and memory of initial conditions. We also assess the variability and potential predictability from these experiments.

Chang, Yehui↗

Boundary-induced classical generalized Gibbs ensemble with angular momentum

We investigate how confinement geometry leads to the emergence of a Generalized Gibbs Ensemble (GGE) in classical systems. Unlike the standard Gibbs ensemble, the GGE includes additional conserved quantities, such as angular momentum, that arise from boundary-induced symmetries. Using analytical arguments based on the maximum entropy principle, we show that circular boundaries preserve angular momentum and drive the system toward a chiral, non-ergodic GGE that violates time-reversal symmetry. This ensemble differs fundamentally from the Gibbs case, producing near-boundary condensation and revealing how geometry alone can alter thermal equilibration. To quantify these effects, we introduce an order parameter measuring deviations from Gibbs behavior and demonstrate that conventional Monte Carlo methods must incorporate angular momentum conservation under such conditions. Our study highlights how geometric constraints shape non-equilibrium statistical ensembles and lead to subtle departures from the Bohr-van Leeuwen theorem. These predictions are validated through detailed simulations of confined classical hard-disk gases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An Observation-Driven Framework for Modeling Post-Fire Hydrologic Response: Evaluation for Two Central California Case Studies

In a warming climate, wildfires are becoming increasingly common, especially in semi-arid environments. Wildfires can disrupt forest ecosystems and induce changes to the land surface. Collectively, these impacts can alter the hydrologic response of a catchment following a fire, resulting in increased potential for surface runoff, reduced evapotranspiration, and, ultimately, a higher risk for flash flooding and mass wasting. The timescale for post-fire recovery of hydrological processes to return to pre-fire conditions is not well established due to the lack of ground measurements. Accurate characterization of fire impacts on hydrologic response is also challenging to simulate, given the complex interplay of various processes. Here, we present a generalized framework to quantify wildfire impacts on runoff generation. We consider the vegetation and soil disturbances as the primary contributors to post-fire floods. Using an ensemble modeling structure to account for parameter uncertainty and consider model sensitivity, remotely sensed leaf area index (LAI) is assimilated into a land surface model (LSM) to simulate vegetation disturbance, and the maximum land surface saturation LSM parameter is decreased to parameterize the soil disturbance following observed fires. We evaluated the impacts of fire-induced changes to LAI and soil saturation on hydrologic states like runoff and evapotranspiration. Analysis of two events in central California shows that 1) the LAI assimilation has a greater impact on water balance; and 2) the soil scheme captures a range of outcomes, with substantial impacts on surface runoff for some ensemble members, which could affect flood potential.

Timothy M Lahmers↗

An Observation-Driven Framework for Modeling Post-Fire Hydrologic Response: Evaluation for Two Central California Case Studies

In a warming climate, wildfires are becoming increasingly common, especially in semi-arid environments. Wildfires can disrupt forest ecosystems and induce changes to the land surface. Collectively, these impacts can alter the hydrologic response of a catchment following a fire, resulting in increased potential for surface runoff, reduced evapotranspiration, and, ultimately, a higher risk for flash flooding and mass wasting. The timescale of post-fire recovery of hydrological processes to return to pre-fire conditions is not well established due to the lack of ground measurements. Accurate characterization of the impacts of fire on hydrologic response is also challenging to simulate, given the complex interplay of various processes. Here, we present a generalized framework to quantify the impacts of wildfire on runoff generation. We consider the disturbances in the vegetation and soil as the two main factors contributing to post-fire floods. Using an ensemble modeling structure to account for parameter uncertainty, remotely sensed leaf area index (LAI) is assimilated into a land surface model (LSM) to simulate vegetation disturbance, and the maximum land surface saturation LSM parameter is decreased to parameterize the soil disturbance following observed fires. We consider the impacts of fire-induced changes to LAI and soil saturation on hydrologic states like runoff and evapotranspiration for two case studies. These case studies demonstrate the general applicability of hydrophobicity formulation to serve as a guideline for exploring the range of hydrologic responses post-fire.

Fire Hydrology↗

Development and Evaluation of High-Resolution Climate Simulations Over the Mountainous Northeastern United States

The mountain regions of the northeastern United States are a critical socioeconomic resource for Vermont, New York State, New Hampshire, Maine, and southern Quebec. While global climate models (GCMs) are important tools for climate change risk assessment at regional scales, even the increased spatial resolution of statistically downscaled GCMs (commonly approximately 1/ 8 deg) is not sufficient for hydrologic, ecologic, and land-use modeling of small watersheds within the mountainous Northeast. To address this limitation, an ensemble of topographically downscaled, high-resolution (30"), daily 2-m maximum air temperature; 2-m minimum air temperature; and precipitation simulations are developed for the mountainous Northeast by applying an additional level of downscaling to intermediately downscaled (1/ 8 deg) data using high-resolution topography and station observations. First, observed relationships between 2-m air temperature and elevation and between precipitation and elevation are derived. Then, these relationships are combined with spatial interpolation to enhance the resolution of intermediately downscaled GCM simulations. The resulting topographically downscaled dataset is analyzed for its ability to reproduce station observations. Topographic downscaling adds value to intermediately downscaled maximum and minimum 2-m air temperature at high-elevation stations, as well as moderately improves domain-averaged maximum and minimum 2-m air temperature. Topographic downscaling also improves mean precipitation but not daily probability distributions of precipitation. Overall, the utility of topographic downscaling is dependent on the initial bias of the intermediately downscaled product and the magnitude of the elevation adjustment. As the initial bias or elevation adjustment increases, more value is added to the topographically downscaled product.

climatology↗