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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Meridional Propagation of the MJO/ISO and Asian Monsoon Variability

In this study we examine the links between tropical heating, the Madden Julian Oscillation (MJO)/Intraseasonal Oscillation (ISO), and the Asian monsoon. We are particularly interested in isolating the nature of the poleward propagation of the ISO/MJO in the monsoon region. We examine both observations and idealized "MJO heating" experiments employing the NASA Seasonal-Interannual Prediction Project (NSIPP) atmospheric general circulation model (AGCM). In the idealized 10-member ensemble simulations, the model is forced by climatological SST and an idealized eastward propagating heating profile that is meant to mimic the canonical heating associated with the MJO in the Indian Ocean and western Pacific. In order to understand the impact of SST on the off equatorial convection (or Rossby-wave response), a second set of 10-member ensemble simulations is carried out with the climatological SSTs shifted in time by 6-months. The observational analysis highlights the strong link between the Indian summer monsoon and the tropical ISO/MJO activity and heating. This includes the well-known meridional propagation that affects the summer monsoons of both hemispheres. The AGCM experiments with the idealized eastward propagating MJO-like heating reproduce the observed meridional propagation including the observed seasonal differences. The impact of the SSTs are to enhance the magnitude of the propagation into the summer hemispheres. The results suggest that the winter/summer differences associated with the MJO/ISO are auxiliary features that depend on the MJO's environment (basic state and boundary conditions) and are not the result of fundamental differences in the MJO itself.

Wu, Man Li C.↗
NASA NTRS2018

Quantifying the Uncertainties in an Ensemble of Decadal Climate Predictions

Meaningful climate predictions should be accompanied by the corresponding uncertainty range. Common methods for estimating the uncertainty range are based on the spread of ensemble predictions. However, a simulation ensemble is not necessarily a proper sample of the real distribution of the climate, and therefore, the ensemble spread cannot be interpreted as the actual uncertainty. We propose a new method that links between the ensemble spread and the uncertainty without relying on any assumptions regarding the distribution of the ensemble predictions. The method is tested using CMIP5 1981-2010 decadal predictions and is shown to outperform other common methods.

Strobach, Ehud↗
NASA NTRS2001

Statistical Ensemble of Large Eddy Simulations

A statistical ensemble of large eddy simulations (LES) is run simultaneously for the same flow. The information provided by the different large scale velocity fields is used to propose an ensemble averaged version of the dynamic model. This produces local model parameters that only depend on the statistical properties of the flow. An important property of the ensemble averaged dynamic procedure is that it does not require any spatial averaging and can thus be used in fully inhomogeneous flows. Also, the ensemble of LES's provides statistics of the large scale velocity that can be used for building new models for the subgrid-scale stress tensor. The ensemble averaged dynamic procedure has been implemented with various models for three flows: decaying isotropic turbulence, forced isotropic turbulence, and the time developing plane wake. It is found that the results are almost independent of the number of LES's in the statistical ensemble provided that the ensemble contains at least 16 realizations.

Carati, Daniele↗
NASA NTRS2010

Effects of Irrigation on Global Climate During the 20th Century

Various studies have documented the effects of modern ]day irrigation on regional and global climate, but none, to date, have considered the time ]varying impact of steadily increasing irrigation rates on climate during the 20th century. We investigate the impacts of observed irrigation changes over this century with two ensemble simulations using an atmosphere general circulation model. Both ensembles are forced with transient climate forcings and observed sea surface temperatures from 1902 to 2000; one ensemble includes irrigation specified by a time ]varying data set of irrigation water withdrawals. Early in the century, irrigation is primarily localized over southern and eastern Asia, leading to significant cooling in boreal summer (June.August) over these regions. This cooling spreads and intensifies by century fs end, following the rapid expansion of irrigation over North America, Europe, and Asia. Irrigation also leads to boreal winter (December.February) warming over parts of North America and Asia in the latter part of the century, due to enhanced downward longwave fluxes from increased near ]surface humidity. Precipitation increases occur primarily downwind of the major irrigation areas, although precipitation in parts of India decreases due to a weaker summer monsoon. Irrigation begins to significantly reduce temperatures and temperature trends during boreal summer over the Northern Hemisphere midlatitudes and tropics beginning around 1950; significant increases in precipitation occur in these same latitude bands. These trends reveal the varying importance of irrigation ]climate interactions and suggest that future climate studies should account for irrigation, especially in regions with unsustainable irrigation resources.

Puma, M. J.↗
DOE OSTI2024

Uncertainty Visualization Challenges in Decision Systems with Ensemble Data & Surrogate Models

Uncertainty visualization is a key component in translating important insights from ensemble simulation data into actionable decision-making by visually conveying various aspects of uncertainty within a system. With the recent advent of fast surrogate models trained on ensemble data, we can substitute computationally expensive simulations, which allows users to interact with more aspects of data spaces than ever before. However, the use of ensemble data with surrogate models in a decision-making tool brings up new challenges for uncertainty visualization, namely how to reconcile and communicate the new and different types of uncertainties brought in by surrogates and how to utilize these new data estimates in actionable ways. In this work, we examine these issues as they relate to high-dimensional data visualization, the integration of discrete datasets and the continuous representations of those datasets, and the unique difficulties associated with systems that allow users to iterate between input and output spaces. We assess the role of uncertainty visualization in facilitating intuitive and actionable interaction with ensemble data and surrogate models, and highlight key challenges in this new frontier of computational simulation.

ensemble data↗
DOE OSTI2026

Storylines for the 1997 New Year’s Flood: The role of watershed antecedent conditions and future warming in shaping discharge in the Truckee River watershed

The 1997 New Year’s flood was among the most devastating floods in the Truckee River watershed located in western Nevada. This event resulted from complex interactions of flood drivers, such as extreme precipitation, wet antecedent watershed conditions, warm temperatures and rapid snowmelt. We leveraged simulated forcings from the regionally refined mesh capabilities of the Energy Exascale Earth System Model (RRM-E3SM) and a process-based hydrological model to recreate the 1997 New Year’s flood for the Truckee River watershed across four climate warming levels ranging from the current temperatures to + 4° C. For each scenario, we conducted ensemble simulations with the same forcing but with 100 different seasonal watershed antecedent conditions, which were randomly sampled from long-term hydrological simulations. The results show that the 1997 New Year’s flood can be reproduced or exceeded consistently only when the antecedent watershed conditions are wet, specifically when streamflows are above the 75th percentile of the climatological value. There is negligible change in ensemble mean peakflows for Truckee River near Reno; however, there are increases of 18% and 14% under the warming levels of + 3° C and + 4° C, respectively. The increases in peakflows under future climate warming are attributed to wetter antecedent watershed conditions and enhanced snowmelt. Furthermore, the largest increases in peakflows occur at small, high-elevation headwater basins along the Sierra Nevada crest. This study highlights that changes in extreme flood events will result from the complex interplay of multiple flood drivers. It also demonstrates the potential of storyline approaches to analyze future realizations of these extreme events under different climate scenarios.

Climate change impact study↗
DOE OSTI2025

Navigating Uncertainty: Challenges in Visualizing Ensemble Data and Surrogate Models for Decision Systems

Uncertainty visualization plays a critical role in transforming ensemble simulation data into actionable insights by effectively communicating various dimensions of uncertainty within a system. The emergence of artificial intelligence-driven surrogate models trained on multirun ensemble data offers a transformative opportunity to replace computationally intensive simulations with fast estimates, enabling users to explore data spaces with unprecedented depth and interactivity. However, integrating ensemble data and surrogate models into decision-making workflows and tools introduces novel challenges for uncertainty visualization. These include reconciling and clearly communicating the unique uncertainties associated with ensembles and their surrogate model estimates, and leveraging these approximations to inform actionable decisions. This work explores these challenges in the context of high-dimensional data visualization, bridging discrete datasets with their continuous representations and addressing the complexities of systems that support iterative navigation between input and output spaces. We evaluate the role of uncertainty visualization in fostering intuitive, actionable interactions and identify critical hurdles in advancing this frontier of computational simulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Using Multiple Isotope-Labeled Infrared Spectra for the Structural Characterization of an Intrinsically Disordered Peptide

Intrinsically disordered proteins (IDPs) rapidly interconvert between conformers, requiring an ensemble description. This complicates their experimental characterization, and force field limitations pose challenges for their simulation. Here, in this work, we use isotope-labeled and unlabeled infrared (IR) spectra to reweight simulated ensembles of the elastin-like peptide GVGVPGVG, a paradigmatic disordered peptide. By comparing the results obtained with different spectra, we explicitly show that the weights are underdetermined by the ensemble averaged data. We identify which labels and frequency regions maximize structural information while minimizing sensitivity to simulation error and show that these regions report on whether the peptide makes specific interactions. Our work shows the importance of incorporating simulations and simulated spectra at the planning stages of isotope-labeled IR experiments and more generally provides a framework for interpreting IR data for IDPs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2023

AgMIP-Wheat Multi-Model Simulations on Climate Change Impact and Adaptation for Global Wheat

The climate change impact and adaptation simulations from the Agricultural Model Intercomparison and Improvement Project (AgMIP) for wheat provide a unique dataset of multi-model ensemble simulations for 60 representative global locations covering all global wheat mega environments. The multi-model ensemble reported here has been thoroughly benchmarked against a large number of experimental data, including different locations, growing season temperatures, atmospheric CO2 concentration, heat stress scenarios, and their interactions. In this paper, we describe the main characteristics of this global simulation dataset. Detailed cultivar, crop management, and soil datasets were compiled for all locations to drive 32 wheat growth models. The dataset consists of 30-year simulated data including 25 output variables for nine climate scenarios, including Baseline (1980-2010) with 360 or 550 ppm CO2, Baseline +2oC or +4oC with 360 or 550 ppm CO2, a mid-century climate change scenario (RCP8.5, 571 ppm CO2), and 1.5°C (423 ppm CO2) and 2.0oC (487 ppm CO2) warming above the pre-industrial period (HAPPI). This global simulation dataset can be used as a benchmark from a well-tested multi-model ensemble in future analyses of global wheat. Also, resource use efficiency (e.g., for radiation, water, and nitrogen use) and uncertainty analyses under different climate scenarios can be explored at different scales. The DOI for the dataset is 10.5281/zenodo.4027033 (AgMIP-Wheat, 2020), and all the data are available on the data repository of Zenodo (http://doi.org/10.5281/zenodo.4027033). Two scientific publications have been published based on some of these data here.

Agricultural Model Intercomparison and Improvement↗
DOE OSTI2025

Plateau to River Model Predictive Simulations for All Ensemble Realizations to Support Modeling Work in Fiscal Year 2025

The purpose of this environmental calculation file (ECF) is to document predictions of flow and hydraulic head on the Central Plateau of the Hanford Site using the Plateau-to-River (P2R) Model (CP-57037, Model Package Report for the Plateau-to-River Model: Version 9.1). This calculation documents the simulation of the groundwater for the parent model domain of the P2R Model as a basis for use in other applications of the P2R Model. This application is unique from the standpoint that it will simulate all ensemble member models of the P2R Model whereas other applications may only utilize specific ensemble members. These simulations provide results that can be used in the process of selecting an appropriate subset of ensemble members for other applications.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2013

Downscaling Reanalysis over Continental Africa with a Regional Model: NCEP Versus ERA Interim Forcing

Five annual climate cycles (1998-2002) are simulated for continental Africa and adjacent oceans by a regional atmospheric model (RM3). RM3 horizontal grid spacing is 0.44deg at 28 vertical levels. Each of 2 simulation ensembles is driven by lateral boundary conditions from each of 2 alternative reanalysis data sets. One simulation downs cales National Center for Environmental Prediction reanalysis 2 (NCPR2) and the other the European Centre for Medium Range Weather Forecasts Interim reanalysis (ERA-I). NCPR2 data are archived at 2.5deg grid spacing, while a recent version of ERA-I provides data at 0.75deg spacing. ERA-I-forced simulations are recomrp. ended by the Coordinated Regional Downscaling Experiment (CORDEX). Comparisons of the 2 sets of simulations with each other and with observational evidence assess the relative performance of each downscaling system. A third simulation also uses ERA-I forcing, but degraded to the same horizontal resolution as NCPR2. RM3-simulated pentad and monthly mean precipitation data are compared to Tropical Rainfall Measuring Mission (TRMM) data, gridded at 0.5deg, and RM3-simulated circulation is compared to both reanalyses. Results suggest that each downscaling system provides advantages and disadvantages relative to the other. The RM3/NCPR2 achieves a more realistic northward advance of summer monsoon rains over West Africa, but RM3/ERA-I creates the more realistic monsoon circulation. Both systems recreate some features of JulySeptember 1999 minus 2002 precipitation differences. Degrading the resolution of ERA-I driving data unrealistically slows the monsoon circulation and considerably diminishes summer rainfall rates over West Africa. The high resolution of ERA-I data, therefore, contributes to the quality of the downscaling, but NCPR2laterai boundary conditions nevertheless produce better simulations of some features.

Reanalysis↗
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↗
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.↗
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 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 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 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↗