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 91 records · Page 5

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↗

The DOE E3SM coupled model version 1: Overview and evaluation at standard resolution

This work documents the first version of the U.S. Department of Energy (DOE) new Energy Exascale Earth System Model (E3SMv1). We focus on the standard resolution of the fully--coupled physical model designed to address DOE mission-relevant water cycle questions. Its components include atmosphere and land (110km grid spacing), ocean and sea ice (60km in the mid-latitudes and 30km at the equator and poles), and river transport (55km) models. This base configuration will also serve as a foundation for additional configurations exploring higher horizontal resolution as well as augmented capabilities in the form of biogeochemistry and cryosphere configurations. The performance of E3SMv1 is evaluated by means of a standard set of Coupled Model Intercomparison Project Phase 6 Diagnosis, Evaluation, and Characterization of Klima (CMIP6 DECK) simulations consisting of a long pre-industrial control, historical simulations (ensembles of fully coupled and prescribed SSTs) as well as idealized CO 2 forcing simulations. The model performs well overall with biases typical of other CMIP-class models, although the simulated Atlantic Meridional Overturning Circulation is weaker than many CMIP-class models. While the E3SMv1 historical ensemble captures the bulk of the observed warming between pre-industrial (1850) and present-day, the trajectory of the warming diverges from observations in the second half of the 20th century with a period of delayed warming followed by an excessive warming trend. Using a two-layer energy balance model, we attribute this divergence to the model's strong aerosol-related effective radiative forcing (ERF ari+aci = -1.65 W m -2 ) and high equilibrium climate sensitivity (ECS = 5.3 K).

58 GEOSCIENCES↗

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↗

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↗

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↗

Model Resolution Sensitivity of the Simulation of North Atlantic Oscillation Teleconnections to Precipitation Extremes

We evaluate a high-resolution (0.25°), four-member ensemble simulation of the global climate (1979–2005) with the U.S. Department of Energy's Energy Exascale Earth System Model v0.3—forced with observed ocean surface temperatures and sea ice extent—for its ability to represent the North Atlantic Oscillation (NAO) teleconnections to winter precipitation extremes over western Europe. As compared to the low-resolution model (1°), it simulates a stronger impact of NAO on daily precipitation extremes over the western slopes of mountain ranges over southwestern Norway, northwestern United Kingdom, and the Western Balkan states. Precipitation extremes and their linear relationship with NAO are quantified using the generalized extreme value distribution. NAO-dependent large-scale (stratiform) precipitation intensity strengthens in the high-resolution model on seasonal time scales and plays a dominant role during simulated daily precipitation extremes. Improvements in the high-resolution model over these varied-topography regions largely appear to be due to finer resolved scales of motion that amplify NAO-dependent seasonal vertical moisture fluxes and enhance stable condensation. However, the high-resolution model simulates a weaker than observed impact of NAO on extratropical cyclone activity and total precipitable water, generally underperforming the low-resolution model These effects possibly offset the impact of enhanced vertical moisture fluxes on NAO-dependent precipitation extremes in the high-resolution model in these regions. Over the southwestern Iberian peninsula, the high-resolution model underperforms the low-resolution model simulating weaker than observed impact of NAO on precipitation extremes. Finally, this appears to be due to the reduction in total precipitable water despite an increase in NAO-dependent extratropical activity there.

54 ENVIRONMENTAL SCIENCES↗

How Would the Twenty-First-Century Warming Influence Pacific Decadal Variability and Its Connection to North American Rainfall: Assessment Based on a Revised Procedure for the IPO/PDO

Decadal climate variability of sea surface temperature (SST) over the Pacific Ocean can be characterized by interdecadal Pacific oscillation (IPO) or Pacific decadal oscillation (PDO) based on empirical orthogonal function (EOF) analysis. Although the procedures to derive the IPO and PDO indices differ in their regional focuses and filtering methods to remove interannual variability, the IPO and PDO are highly correlated in time and are often used interchangeably. Studies have shown that the IPO and PDO conjointly (IPO/PDO for conciseness) play a vital role in modulating the pace of global warming. It is less clear, however, how externally forced global warming may, in turn, affect the IPO/PDO. One obstacle to revealing this effect is that the conventional definitions of the IPO/PDO fail to account for the spatial heterogeneity of the background warming trend, which causes the IPO/PDO to be conflated with the warming trend, especially for the twenty-first-century simulation when the forced change is likely to be more dominant. Using a large-ensemble simulation in the Community Earth System Model, version 1 (CESM1), it is shown here that a better practice of detrending prior to EOF analysis is to remove the local and nonlinear trend, defined as the ensemble-mean time series at each grid box (or simply as the quadratic fit of the local time series if such an ensemble is not available). The revised IPO/PDO index is purely indicative of internal decadal variability. In the twenty-first-century warmer climate, the IPO/PDO has a weaker amplitude in space, a higher frequency in time, and a muted impact on global and North American temperature and rainfall.

54 ENVIRONMENTAL SCIENCES↗

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↗

Surface and top-of-atmosphere radiative feedback kernels for CESM-CAM5

Radiative kernels at the top of the atmosphere are useful for decomposing changes in atmospheric radiative fluxes due to feedbacks from atmosphere and surface temperature, water vapor, and surface albedo. Here we describe and validate radiative kernels calculated with the large-ensemble version of CAM5, CESM1.1.2, at the top of the atmosphere and the surface. Estimates of the radiative forcing from greenhouse gases and aerosols in RCP8.5 in the CESM large-ensemble simulations are also diagnosed. As an application, feedbacks are calculated for the CESM large ensemble. The kernels are freely available at https://doi.org/10.5065/D6F47MT6 , and accompanying software can be downloaded from https://github.com/apendergrass/cam5-kernels .

54 ENVIRONMENTAL SCIENCES↗

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↗

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↗

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↗

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↗

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↗

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

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↗

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↗

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↗