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 523 records · Page 29

Quantum Simulations of Radiation Damage in a Molecular Polyethylene Analog

Abstract An atomic‐level understanding of radiation‐induced damage in simple polymers like polyethylene is essential for determining how these chemical changes can alter the physical and mechanical properties of important technological materials such as plastics. Ensembles of quantum simulations of radiation damage in a polyethylene analog are performed using the Density Functional Tight Binding method to help bind its radiolysis and subsequent degradation as a function of radiation dose. Chemical degradation products are categorized with a graph theory approach, and occurrence rates of unsaturated carbon bond formation, crosslinking, cycle formation, chain scission reactions, and out‐gassing products are computed. Statistical correlations between product pairs show significant correlations between chain scission reactions, unsaturated carbon bond formation, and out‐gassing products, though these correlations decrease with increasing atom recoil energy. The results present relatively simple chemical descriptors as possible indications of network rearrangements in the middle range of excitation energies. Ultimately, the work provides a computational framework for determining the coupling between nonequilibrium chemistry in polymers and potential changes to macro‐scale properties that can aid in the interpretation of future radiation damage experiments on plastic materials.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Multimodel Ensembles of Wheat Growth: Many Models are Better than One

Crop models of crop growth are increasingly used to quantify the impact of global changes due to climate or crop management. Therefore, accuracy of simulation results is a major concern. Studies with ensembles of crop model scan give valuable information about model accuracy and uncertainty, but such studies are difficult to organize and have only recently begun. We report on the largest ensemble study to date, of 27 wheat models tested in four contrasting locations for their accuracy in simulating multiple crop growth and yield variables. The relative error averaged over models was 2438 for the different end-of-season variables including grain yield (GY) and grain protein concentration (GPC). There was little relation between error of a model for GY or GPC and error for in-season variables. Thus, most models did not arrive at accurate simulations of GY and GPC by accurately simulating preceding growth dynamics. Ensemble simulations, taking either the mean (e-mean) or median (e-median) of simulated values, gave better estimates than any individual model when all variables were considered. Compared to individual models, e-median ranked first in simulating measured GY and third in GPC. The error of e-mean and e-median declined with an increasing number of ensemble members, with little decrease beyond 10 models. We conclude that multimodel ensembles can be used to create new estimators with improved accuracy and consistency in simulating growth dynamics. We argue that these results are applicable to other crop species, and hypothesize that they apply more generally to ecological system models.

model intercomparison↗

Changes in Precipitation From North Atlantic Tropical Cyclones Under RCP Scenarios in the Variable-Resolution Community Atmosphere Model

Decreasing climate models' grid spacing improves the representation of tropical cyclones at decadal time scales. In this study, a variable-resolution (VR) version of the Community Atmosphere Model 5 (CAM5-VR) is utilized to study North Atlantic tropical cyclone climatology in ensemble historical climate simulations and under two Representative Concentration Pathway (RCP) projections (RCP4.5 and RCP8.5). Basin-wide tropical cyclone counts decrease in the RCP simulations, although landfalling storm counts do not show as straightforward of a pattern, especially when focusing on regional changes. Lifetime maximum intensity metrics suggest that tropical cyclones increase in strength in the RCP ensembles. However, despite increases in tropical cyclone-related precipitation rates and the amount of precipitation produced per storm with warming, the annual average Rx5day from tropical cyclones over the eastern United States decreases due to less landfalling storms. This work is part of a continued effort to quantify how tropical cyclone-induced hazards may change in future climates.

54 ENVIRONMENTAL SCIENCES↗

A teleconnection between Atlantic sea surface temperature and eastern and central North Pacific tropical cyclones

The El Niño–Southern Oscillation (ENSO) is a major source of seasonal tropical cyclone (TC) predictability in both local and remote ocean basins. Unusually warm eastern-central equatorial Pacific sea surface temperature (SST) during El Niño tends to enhance eastern and central North Pacific (ECNP) TCs and suppress Atlantic TCs. Here in this paper, we demonstrate that Atlantic SST variability likewise influences remote TC activity in the eastern-central Pacific through a Walker Circulation-type response analogous to the ENSO-Atlantic TC teleconnection, using observations and 27 km resolution tropical channel model (TCM) simulations. Observed and simulated ECNP TC activity is reduced during the positive Atlantic Meridional Mode (AMM), which is characterized by warm northern and cool southern tropical Atlantic SST anomalies, and vice versa during the negative AMM. Large ensembles of TCM simulations indicate that SST variability, rather than internal atmospheric variability, drives extreme ECNP hurricane seasons.

54 ENVIRONMENTAL SCIENCES↗

Computational chemistry

With the advent of supercomputers, modern computational chemistry algorithms and codes, a powerful tool was created to help fill NASA's continuing need for information on the properties of matter in hostile or unusual environments. Computational resources provided under the National Aerodynamics Simulator (NAS) program were a cornerstone for recent advancements in this field. Properties of gases, materials, and their interactions can be determined from solutions of the governing equations. In the case of gases, for example, radiative transition probabilites per particle, bond-dissociation energies, and rates of simple chemical reactions can be determined computationally as reliably as from experiment. The data are proving to be quite valuable in providing inputs to real-gas flow simulation codes used to compute aerothermodynamic loads on NASA's aeroassist orbital transfer vehicles and a host of problems related to the National Aerospace Plane Program. Although more approximate, similar solutions can be obtained for ensembles of atoms simulating small particles of materials with and without the presence of gases. Computational chemistry has application in studying catalysis, properties of polymers, all of interest to various NASA missions, including those previously mentioned. In addition to discussing these applications of computational chemistry within NASA, the governing equations and the need for supercomputers for their solution is outlined.

Arnold, J. O.↗

A Regional Model Study of Synoptic Features Over West Africa

Synoptic weather features over West Africa were studied in simulations by the regional simulation model (RM) at the NASA/Goddard Institute for Space Studies. These pioneering simulations represent the beginning of an effort to adapt regional models for weather and climate prediction over West Africa. The RM uses a cartesian grid with 50 km horizontal resolution and fifteen vertical levels. An ensemble of four simulations was forced with lateral boundary conditions from ECMWF global analyses for the period 8-22 August 1988. The simulated mid-tropospheric circulation includes the skillful development and movement of several African wave disturbances. Wavelet analysis of mid-tropospheric winds detected a dominant periodicity of about 4 days and a secondary periodicity of 5-8 days. Spatial distributions of RM precipitation and precipitation time series were validated against daily rain gauge measurements and ISCCP satellite infrared cloud imagery. The time-space distribution of simulated precipitation was made more realistic by combining the ECMWR initial conditions with a 24-hr spin-up of the moisture field and also by damping high frequency gravity waves by dynamic initialization. Model precipitation "forecasts" over the Central Sahel were correlated with observations for about three days, but reinitializing with observed data on day 5 resulted in a dramatic improvement in the precipitation validation over the remaining 9 days. Results imply that information via the lateral boundary conditions is not always sufficient to minimize departures between simulated and actual precipitation patterns for more than several days. In addition, there was some evidence that the new initialization may increase the simulations' sensitivity to the quality of lateral boundary conditions.

Druyan, Leonard M.↗

Historical and future trends in South Asian monsoon low pressure systems in a high-resolution model ensemble

Historical trends in monsoon low pressure systems (LPS), the dominant rain-bearing weather system of South Asia, have been difficult to assess due to changes in the observing network. Future projections have also remained uncertain because prior studies concluded that many coarse-resolution climate models do not accurately simulate LPS. Here, we examine changes in South Asian monsoon LPS simulated by an ensemble of global models, including some with high spatial resolution, that we show skillfully represent LPS. In the ensemble mean, the number of strong LPS (monsoon depressions) decreased over the last 65 years (1950–2014) by about 15% while no trend was detected for weaker LPS (monsoon lows). The reduction in depression counts then moderated, yielding no trend in the periods 1980–2050 or 2015–2050. The ensemble mean projects a shift in genesis from ocean to land and an increase in LPS precipitation of at least 7% K −1 , which together contribute to a projected increase in seasonal mean and extreme precipitation over central India.

54 ENVIRONMENTAL SCIENCES↗

Robust Spring Drying in the Southwestern U.S. and Seasonal Migration of Wet/Dry Patterns in a Warmer Climate

This study compares climate simulations over the United States produced by a regional climate model with the driving global climate simulations as well as a large multi-model ensemble of global climate simulations to investigate robust changes in water availability (precipitation (P) – evapotranspiration (E)). A robust spring dry signal across multiple models is identified in the Southwest that results from a decrease in P and an increase in E in the future. In the boreal winter and summer, the prominent changes in P – E are associated with a north – south dipole pattern, while in spring, the prominent changes in P – E appear as an east – west dipole pattern. The progression of the north – south and east – west dipole patterns through the seasons manifests clearly as a seasonal “clockwise” migration of wet/dry patterns, which is shown to be a robust feature of water availability changes in the US consistent across regional and global climate simulations.

CMIP5↗

PROTON KINETIC EFFECTS IN VLASOV AND SOLAR WIND TURBULENCE

Kinetic plasma processes are investigated in the framework of solar wind turbulence, employing hybrid Vlasov-Maxwell (HVM) simulations. Statistical analysis of spacecraft observation data relates proton temperature anisotropy T /T {sub ∥} and parallel plasma beta β{sub ∥}, where subscripts refer to the ambient magnetic field direction. Here, this relationship is recovered using an ensemble of HVM simulations. By varying plasma parameters, such as plasma beta and fluctuation level, the simulations explore distinct regions of the parameter space given by T /T {sub ∥} and β{sub ∥}, similar to solar wind sub-datasets. Moreover, both simulation and solar wind data suggest that temperature anisotropy is not only associated with magnetic intermittent events, but also with gradient-type structures in the flow and in the density. This connection between non-Maxwellian kinetic effects and various types of intermittency may be a key point for understanding the complex nature of plasma turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗

Resolving Mesoscale Convective Systems: Grid Spacing Sensitivity in the Tropics and Midlatitudes

Abstract Mesoscale convective systems (MCSs) are a critical global water cycle component and drive extreme precipitation events in tropical and midlatitude regions. However, simulating deep convection remains challenging for modern numerical weather and climate models due to the complex interactions of processes from microscales to synoptic scales. Recent models with kilometer‐scale horizontal grid spacings offer notable improvements in simulating deep convection compared to coarser‐resolution models. Still, deficiencies in representing key physical processes, such as entrainment, lead to systematic biases. Additionally, evaluating model outputs using process‐oriented observational data remain difficult. This study presents an ensemble of MCS simulations with spanning the deep convective gray zone ( from 12 km to 125 m) in the Southern Great Plains of the U.S. and the Amazon Basin. Comparing these simulations with Atmospheric Radiation Measurement (ARM) wind profiler observations, we find greater sensitivity in the Amazon Basin compared to the Great Plains. Convective drafts converge structurally at sub‐kilometer scales, but some deficiencies remain. In both regions, simulated up and downdrafts are too deep and extreme downdrafts are not strong enough. Furthermore, Amazonian updrafts are too strong. Overall, we observe higher sensitivity in the tropics, including an artificial buildup in vertical kinetic energy at scales of , suggesting a need for 250 m in this region. Nevertheless, bulk convergence—agreement of storm‐average statistics—is achievable with kilometer‐scale simulations within a 10% error margin with 1 km providing a good balance between accuracy and computational cost.

54 ENVIRONMENTAL SCIENCES↗

Timescales of Soil Moisture Anomalies: Results from Two GCMs

Soil moisture anomalies dissipate over timescales that may span weeks to months. Characterizing the geographical and seasonal variations in these timescales can have important practical benefit; significant soil moisture "memory" allows long-lead forecasts of soil moisture, which have been found in recent studies to be essential for useful Ion--lead forecasts of precipitation in many regions. In this talk, we will present and compare the soil moisture timescales derived in two separate general circulation model (GCM) studies. Both studies employ multiple ensembles of short-term climate simulations. Timescales at a given point are effectively estimated by determining how quickly the soil moisture distribution generated in one ensemble of simulations (characterized by a unique set of initial soil moisture conditions) approaches that produced by another ensemble (characterized by a different set of initial soil moisture conditions). The talk will include a discussion of why the timescales produced by the two GCMs differ in some regions, and it will also describe the impact of soil moisture memory on simulated precipitation.

Koster, Randal D.↗

Grain boundary structure search by using an evolutionary algorithm with effective mutation methods

Grain boundaries (GBs) accommodate the misorientation between adjacent grains in a polycrystalline material. GBs are geometrically described by the macroscopic and microscopic degrees of freedom. Besides, at the atomistic level, GBs exhibit complicated behaviors under varying thermodynamic conditions. The complexity of atomistic GB structures demands stochastic searching for possible states. The effectiveness of stochastic search methods relies on techniques to recreate and select atomistic structures. In this work, we developed a new mutation operator that can induce direct and collective atomistic structure changes to boost the search efficiency of exploring GB structures with evolutionary algorithms (EA). We implemented the mutation methods along with innovative selection, crossover, boundary condition preprocessing methods to form an EA-based package to explore GB structures in grand canonical ensembles with atomistic simulations. We used this package to study the [001] symmetric tilt grain boundaries (STGBs) in FCC copper (Cu), the [110] STGBs in BCC tungsten (W), and the $[1\bar{2}10]$ STGBs in HCP magnesium (Mg). The results show that our design and implementation based on new mutation procedures, selection, and boundary conditions provide a high-quality search of atomistic GB structures in the grand canonical ensemble for different crystal lattices.

36 MATERIALS SCIENCE↗

Iterative sampling of expensive simulations for faster deep surrogate training

Deep neural network (DNN) surrogates of expensive physics simulations are enabling a rapid change in the way that common experimental design and analysis tasks are approached. Surrogate models allow simulations to be performed in parallel and separately from downstream tasks, thereby enabling analyses that would be impossible with the simulation in-the-loop; surrogates based on DNNs can effectively emulate diverse non-scalar data of the types collected in fusion and laboratory-astrophysics experiments. The challenge is in training the surrogate model, for which large ensembles of physics simulations must be run, preferably without wasting computational effort on uninteresting simulations. Here, in this paper, we present an iterative sampling scheme that can preferentially propose simulations in interesting regions of parameter space without neglecting unexplored regions, allowing high-quality and wide-ranging surrogate models to be trained using 2–3 times fewer simulations compare to space-filling designs. Our approach uses an explicit importance function defined on the simulation output space, balanced against a measure of simulation density which serves as a proxy for surrogate accuracy. It is easy to implement and can be tuned to find interesting simulations early in the study, allowing surrogates to be trained quickly and refined as new simulations become available; this represents an important step towards the routine generation of deep surrogate models quickly enough to be truly relevant to experimental work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulating the Transport and Rupture of Pollen in the Atmosphere

Pollen, one type of primary biological aerosol particle (PBAP), is emitted from the terrestrial biosphere and can undergo physical changes in the atmosphere via particle rupture. To examine the fate of pollen and its atmospheric processing, a pollen emission and transport scheme is coupled to the Weather Research and Forecasting Model with Chemistry (WRF-Chem). We simulate the emission of pollen and its impacts on the cloud properties and precipitation in the Southern Great Plains from 12 to 19 April 2013, a period with both high pollen emissions and convective activity. We conduct a suite of ensemble runs that simulate primary pollen and three different pollen rupture mechanisms that generate subpollen particles, including (a) high humidity-induced surface rupture, (b) high humidity-induced in-atmosphere plus surface rupture, and (c) lightning-induced rupture, where in-cloud and cloud-to-ground lightning strikes trigger pollen rupture events. When relative humidity is high (>80%), coarse primary pollen (~1 μg m –3 ) is converted into fine subpollen particles (~1.2e –4 μg m –3 ), which produces 80% more subpollen particles than lightning-induced rupture. The in-atmosphere humidity-driven rupture predominantly produces subpollen particles, which is further enhanced during a frontal thunderstorm. During strong convection, vertical updrafts lift primary pollen and subpollen particles (~0.5e –4 μg m –3 ) to the upper troposphere (~12 km) and laterally transports the ruptured pollen in the anvil top outflow. In regions of high pollen and strong convection, ruptured pollen can influence warm cloud formation by decreasing low cloud (<4 km) cloud water mixing ratios and increasing ice phase hydrometeors aloft (>10 km)

54 ENVIRONMENTAL SCIENCES↗

Vegetation Representation Influences Projected Streamflow Changes in the Colorado River Basin

Vegetation parameters for the Variable Infiltration Capacity (VIC) hydrologic model were recently updated using observations from the Moderate Resolution Imaging Spectroradiometer (MODIS). Previous work showed that these MODIS-based parameters improved VIC evapotranspiration simulations when compared to eddy covariance observations. Due to the importance of evapotranspiration within the Colorado River basin, this study provided a basin-by-basin calibration of VIC soil parameters with updated MODIS-based vegetation parameters to improve streamflow simulations. Interestingly, while both configurations had similar historic streamflow performance, end-of-century hydrologic projections, driven by 29 downscaled global climate models under the RCP8.5 emissions scenario, differed between the two configurations. The calibrated MODIS-based configuration had an ensemble mean that simulated little change in end-of-century annual streamflow volume (+0.4%) at Lees Ferry, Arizona, relative to the historical period (1960–2005). In contrast, the previous VIC configuration, which is used to inform decisions about future water resources in the Colorado River basin, projected an 11.7% decrease in annual streamflow. Both VIC configurations simulated similar amounts of evapotranspiration in the historical period. However, the MODIS-based VIC configuration did not show as much of an increase in evapotranspiration by the end of the century, primarily within the upper basin’s forested areas. Differences in evapotranspiration projections were the result of the MODIS-based vegetation parameters having lower leaf area index values and less forested area compared to previous vegetation estimates used in recent Colorado River basin hydrologic projections. These results highlight the need to accurately characterize vegetation and better constrain climate sensitivities in hydrologic models.

54 ENVIRONMENTAL SCIENCES↗

Baseline evaluation of the impact of updates to the MIT Earth System Model on its model parameter estimates

Abstract. For over 20 years, the Massachusetts Institute of Technology Earth System Model (MESM) has been used extensively for climate change research. The model is under continuous development with components being added and updated. To provide transparency in the model development, we perform a baseline evaluation by comparing model behavior and properties in the newest version to the previous model version. In particular, changes resulting from updates to the land surface model component and the input forcings used in historical simulations of climate change are investigated. We run an 1800-member ensemble of MESM historical climate simulations where the model parameters that set climate sensitivity, the rate of ocean heat uptake, and the net anthropogenic aerosol forcing are systematically varied. By comparing model output to observed patterns of surface temperature changes and the linear trend in the increase in ocean heat content, we derive probability distributions for the three model parameters. Furthermore, we run a 372-member ensemble of transient climate simulations where all model forcings are fixed and carbon dioxide concentrations are increased at the rate of 1 % year −1 . From these runs, we derive response surfaces for transient climate response and thermosteric sea level rise as a function of climate sensitivity and ocean heat uptake. We show that the probability distributions shift towards higher climate sensitivities and weaker aerosol forcing when using the new model and that the climate response surfaces are relatively unchanged between model versions. Because the response surfaces are independent of the changes to the model forcings and similar between model versions with different land surface models, we suggest that the change in land surface model has limited impact on the temperature evolution in the model. Thus, we attribute the shifts in parameter estimates to the updated model forcings.

58 GEOSCIENCES↗

Biased Estimates of Changes in Climate Extremes From Prescribed SST Simulations

Large climate model ensembles are widely used to quantify changes in climate extremes. Here we demonstrate that model-based estimates of changes in the probability of temperature extremes at 1.5°C global warming regionally differ if quantified using prescribed sea surface temperatures (SSTs) instead of using a fully coupled climate model. Based on the identical climate model used in two experimental setups, we demonstrate that particularly over the tropics and Australia estimates of the changes in the odds of annual temperature extremes can be up to more than a factor of 5 to 10 larger using prescribed SSTs rather than a fully coupled model configuration. The two experimental designs imply a different perspective on framing projections. If experiments conditional on prescribed observed SSTs are interpreted as unconditional real-world projections, they project changes in extremes that are systematically biased high and overconfident. Our results illustrate the importance of carefully considering experimental design when interpreting projections of extremes.

54 ENVIRONMENTAL SCIENCES↗