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 649 records · Page 36

On the Correspondence between Mean Forecast Errors and Climate Errors in CMIP5 Models

The present study examines the correspondence between short- and long-term systematic errors in five atmospheric models by comparing the 16 five-day hindcast ensembles from the Transpose Atmospheric Model Intercomparison Project II (Transpose-AMIP II) for July–August 2009 (short term) to the climate simulations from phase 5 of the Coupled Model Intercomparison Project (CMIP5) and AMIP for the June–August mean conditions of the years of 1979–2008 (long term). Because the short-term hindcasts were conducted with identical climate models used in the CMIP5/AMIP simulations, one can diagnose over what time scale systematic errors in these climate simulations develop, thus yielding insights into their origin through a seamless modeling approach. The analysis suggests that most systematic errors of precipitation, clouds, and radiation processes in the long-term climate runs are present by day 5 in ensemble average hindcasts in all models. Errors typically saturate after few days of hindcasts with amplitudes comparable to the climate errors, and the impacts of initial conditions on the simulated ensemble mean errors are relatively small. This robust bias correspondence suggests that these systematic errors across different models likely are initiated by model parameterizations since the atmospheric large-scale states remain close to observations in the first 2–3 days. However, biases associated with model physics can have impacts on the large-scale states by day 5, such as zonal winds, 2-m temperature, and sea level pressure, and the analysis further indicates a good correspondence between short- and long-term biases for these large-scale states. Therefore, improving individual model parameterizations in the hindcast mode could lead to the improvement of most climate models in simulating their climate mean state and potentially their future projections.

54 ENVIRONMENTAL SCIENCES↗

Shallow-to-Deep Transition of Continental Moist Convection: Cold Pools, Surface Fluxes, and Mesoscale Organization

Large-eddy simulation is used to investigate the effects of cold pools driven by rain evaporation on the shallow-to-deep convection transition over land. The physically consistent methodologies are developed to obtain a time-dependent reference ensemble without cold pools and to apply interactive surface heat fluxes without modeling of surface energy and water budgets. Three different simulation ensembles are contrasted. The reference ensemble, in the spirit of one-dimensional single-column models, eliminates cold pools by horizontally homogenizing negative buoyancy production due to rain evaporation. The additional ensembles complement the reference cold-pool-free ensemble by including cold pools and by applying either interactive or prescribed surface fluxes. Contrasting these ensembles suggests possible improvements of convection parameterization in large-scale models of weather and climate. Without cold pools, the reference ensemble preserves key features of buoyancy-driven cellular convection associated with a field of convective plumes, as assumed in a typical convection parameterization. With cold pools, a significant enhancement of surface heat and moisture fluxes and about an hour delay of their daily maximum is simulated. Cold pools enhance near-surface temperature and moisture standard deviations as well as maxima of the near-surface updraft velocity. They also lead to the reduction of cloud lateral entrainment, deeper vertical development of the cloud layer, and a few-times-larger accumulated surface precipitation. Interactive surface fluxes provide a damping mechanism that noticeably suppresses all these effects. Perhaps surprisingly, cold pools do not significantly change the cloud-base convective mass flux that approximately follows the evolution of surface heat fluxes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accelerating Monte Carlo molecular simulations by reweighting and reconstructing Markov chains: Extrapolation of canonical ensemble averages and second derivatives to different temperature and density conditions

Accurate determination of thermodynamic properties of petroleum reservoir fluids is of great interest to many applications, especially in petroleum engineering and chemical engineering. Molecular simulation has many appealing features, especially its requirement of fewer tuned parameters but yet better predicting capability; however it is well known that molecular simulation is very CPU expensive, as compared to equation of state approaches. We have recently introduced an efficient thermodynamically consistent technique to regenerate rapidly Monte Carlo Markov Chains (MCMCs) at different thermodynamic conditions from the existing data points that have been pre-computed with expensive classical simulation. This technique can speed up the simulation more than a million times, making the regenerated molecular simulation almost as fast as equation of state approaches. In this paper, this technique is first briefly reviewed and then numerically investigated in its capability of predicting ensemble averages of primary quantities at different neighboring thermodynamic conditions to the original simulated MCMCs. Moreover, this extrapolation technique is extended to predict second derivative properties (e.g. heat capacity and fluid compressibility). The method works by reweighting and reconstructing generated MCMCs in canonical ensemble for Lennard-Jones particles. In this paper, system's potential energy, pressure, isochoric heat capacity and isothermal compressibility along isochors, isotherms and paths of changing temperature and density from the original simulated points were extrapolated. Finally, an optimized set of Lennard-Jones parameters (ε, σ) for single site models were proposed for methane, nitrogen and carbon monoxide.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Role of Coupled Feedbacks in the Decadal Variability of the Southern Hemisphere Eddy-Driven Jet

The Southern Hemisphere summertime eddy-driven jet and storm tracks have shifted poleward over the recent few decades. In previous studies, explanations have mainly stressed the influence of external forcing in driving this trend. Here we examine the role of internal tropical SST variability in controlling the austral summer jet’s poleward migration, with a focus on interdecadal time scales. The role of external forcing and internal variability are isolated by using a hierarchy of Community Earth System Model version 1 (CESM1) simulations, including the pre-industrial control, large ensemble, and pacemaker runs. Model simulations suggest that in the early twenty-first century, both external forcing and internal tropical Pacific SST variability are important in driving a positive southern annular mode (SAM) phase and a poleward migration of the eddy-driven jet. Additionally, Tropical Pacific SST variability, associated with the negative phase of the interdecadal Pacific oscillation (IPO), acts to shift the jet poleward over the southern Indian and southwestern Pacific Oceans and intensify the jet in the southeastern Pacific basin, while external forcing drives a significant poleward jet shift in the South Atlantic basin. In response to both external forcing and decadal Pacific SST variability, the transient eddy momentum flux convergence belt in the middle latitudes experiences a poleward migration due to the enhanced meridional temperature gradient, leading to a zonally symmetric southward migration of the eddy-driven jet. This mechanism distinguishes the influence of the IPO on the midlatitude circulation from the dynamical impact of ENSO, with the latter mainly promoting the subtropical wave-breaking critical latitude poleward and pushing the midlatitude jet to higher latitudes.

54 ENVIRONMENTAL SCIENCES↗

Optimization of a Nucleic Acids united-RESidue 2-Point model (NARES-2P) with a maximum-likelihood approach

Coarse-grained models are useful tools to investigate the structural and thermodynamic properties of biomolecules. They are obtained by merging several atoms into one interaction site. Such simplified models try to capture as much as possible information of the original biomolecular system in all-atom representation but the resulting parameters of these coarse-grained force fields still need further optimization. In this paper, a force field optimization method, which is based on maximum-likelihood fitting of the simulated to the experimental conformational ensembles and least-squares fitting of the simulated to the experimental heat-capacity curves, is applied to optimize the Nucleic Acid united-RESidue 2-point (NARES-2P) model for coarse-grained simulations of nucleic acids recently developed in our laboratory. The optimized NARES-2P force field reproduces the structural and thermodynamic data of small DNA molecules much better than the original force field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Role of Tropical Variability in Driving Decadal Shifts in the Southern Hemisphere Summertime Eddy-Driven Jet

The Southern Hemisphere summertime eddy-driven jet and storm tracks have shifted poleward over the recent few decades. In previous studies, explanations have mainly stressed the influence of external forcing in driving this trend. Here we examine the role of internal tropical SST variability in controlling the austral summer jet’s poleward migration, with a focus on interdecadal time scales. The role of external forcing and internal variability are isolated by using a hierarchy of Community Earth System Model version 1 (CESM1) simulations, including the pre-industrial control, large ensemble, and pacemaker runs. Model simulations suggest that in the early twenty-first century, both external forcing and internal tropical Pacific SST variability are important in driving a positive southern annular mode (SAM) phase and a poleward migration of the eddy-driven jet. Tropical Pacific SST variability, associated with the negative phase of the interdecadal Pacific oscillation (IPO), acts to shift the jet poleward over the southern Indian and southwestern Pacific Oceans and intensify the jet in the southeastern Pacific basin, while external forcing drives a significant poleward jet shift in the South Atlantic basin. In response to both external forcing and decadal Pacific SST variability, the transient eddy momentum flux convergence belt in the middle latitudes experiences a poleward migration due to the enhanced meridional temperature gradient, leading to a zonally symmetric southward migration of the eddy-driven jet. This mechanism distinguishes the influence of the IPO on the midlatitude circulation from the dynamical impact of ENSO, with the latter mainly promoting the subtropical wave-breaking critical latitude poleward and pushing the midlatitude jet to higher latitudes.

54 ENVIRONMENTAL SCIENCES↗

Data and Scripts Associated with "Modeling Ecohydrological Responses of Vegetation to Urban Microclimates Using the E3SM Land Model"

This dataset supports the study of vegetation ecohydrological responses to urban microclimates using the land component of the Energy Exascale Earth System Model (ELM) at four urban sites in Knoxville, Tennessee, USA. It includes the model inputs, simulation outputs, and associated scripts for running ELM simulations and analyzing the resulting data. The Model_Inputs folder includes static surface data, satellite-derived phenology (i.e., leaf area index), and atmospheric forcing data used to drive ELM simulations. Detailed descriptions of these datasets are provided in Section 2.3.2 of the associated manuscript. The Model_Outputs folder contains simulation results for the baseline, treatment, and ensemble experiments. Outputs from the baseline and treatment simulations are provided as raw ELM NetCDF files. Because the raw outputs from the 4,000-member ensemble are prohibitively large, the ensemble results are provided as summarized CSV files, which also serve as the source data for Figure 5 of the associated manuscript. The Scripts folder contains three components: E3SM, the core codebase of the Energy Exascale Earth System Model (E3SM); elm-olmt, the Offline Land Model Testbed (OLMT) used to perform the simulations; and knoxville_elm, which contains the analysis scripts used to process model outputs and generate the figures and results presented in the associated manuscript. Additional information is provided in Scripts_readme.txt within the Scripts directory.

Lu, Xiaoman [ORNL] (ORCID:0000000306698780)↗

Multi-Conformer Ensemble Docking to Difficult Protein Targets

We investigate large-scale ensemble docking using five proteins from the Directory of Useful Decoys (DUD, dud.docking.org) for which docking to crystal structures has proven difficult. Molecular dynamics trajectories are produced for each protein and an ensemble of representative conformational structures extracted from the trajectories. Docking calculations are performed on these selected simulation structures and ensemble-based enrichment factors compared with those obtained using docking in crystal structures of the same protein targets or random selection of compounds. We also found simulation-derived snapshots with improved enrichment factors that increased the chemical diversity of docking hits for four of the five selected proteins. A combination of all the docking results obtained from molecular dynamics simulation followed by selection of top-ranking compounds appears to be an effective strategy for increasing the number and diversity of hits when using docking to screen large libraries of chemicals against difficult protein targets.

59 BASIC BIOLOGICAL SCIENCES↗

Learning-Based Inversion-Free Model-Data Integration to Advance Ecosystem Model Prediction

Ecosystem model prediction is important for understanding ecosystem responses to climate change and for management support. Model prediction and quantification of predictive uncertainty of ecosystems have long been investigated. The traditional workflow, which calibrates models to match observations and then uses the calibrated models for predictions, relies heavily on inverse modeling to constrain uncertain parameters in complex forward models. This inversion-based prediction approach is infeasible for complex models with heterogeneous parameter uncertainties and incapable of rapid integration of streaming and multiple sources of data because of the difficulty and computational cost in the model inversion, which is typically ill-posed and can require hundreds of thousands of expensive forward simulations to be performed iteratively. We propose to circumvent inverse modeling by precomputing an ensemble of unconstrained forward simulations and then using machine learning (ML) methods to learn the statistical relationship between simulated observation and prediction quantities. Once the ML model has learned the relationship, it can be used to make predictions of future system behavior with uncertainty quantification based on observations. The proposed learningbased inversion-free model prediction (LIMP) framework is computationally efficient which only requires a few thousands of fully parallelizable forward simulations. Additionally, LIMP can continually update predictions based on streaming observations from multiple locations and sources without necessarily requiring extra model simulations. In this study, we apply LIMP to a regional terrestrial ecosystem model with 47 parameters for testing, refining, and evaluating the approach. We demonstrate that LIMP can be used for efficient model prediction, rapid data assimilation, and cost-effective experimental design for improving robust predictive understanding of ecosystems.

Lu, Dan↗

Linear Time-Invariant Models of a Large Cumulus Ensemble

Abstract Methods in system identification are used to obtain linear time-invariant state-space models that describe how horizontal averages of temperature and humidity of a large cumulus ensemble evolve with time under small forcing. The cumulus ensemble studied here is simulated with cloud-system-resolving models in radiative–convective equilibrium. The identified models extend steady-state linear response functions used in past studies and provide accurate descriptions of the transfer function, the noise model, and the behavior of cumulus convection when coupled with two-dimensional gravity waves. A novel procedure is developed to convert the state-space models into an interpretable form, which is used to elucidate and quantify memory in cumulus convection. The linear problem studied here serves as a useful reference point for more general efforts to obtain data-driven and interpretable parameterizations of cumulus convection.

Meteorology & Atmospheric Sciences↗

FATES Parameters and Output for Parameter Sensitivity at the Panama Barro Colorado Island Testbed

This data package contains a set of Functionally Assembled Terrestrial Ecosystem Simulator (FATES) history files and parameter values for the set of model simulations used in Koven et al. Manuscript "Benchmarking and Parameter Sensitivity of Physiological and Vegetation Dynamics using the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) at Barro Colorado Island, Panama." It includes several ensembles of FATES runs, where each ensemble is comprised of a set of 576 runs that have different values of PFT traits. The ensembles include 1 PFT, 2 PFT, 3 PFT, and 10 PFT simulations. Different ensembles also have different values of disturbance and height sorting parameters, and structural differences including running FATES within different LSMs.

54 ENVIRONMENTAL SCIENCES↗

Mt. Pinatubo's Impacts on the GEOS Forecasting System

The eruption of Mount Pinatubo in June 1991 had dramatic effects on the global climate, introducing a sudden and extreme forcing on the radiative budget. In this presentation we analyze the effects of the Mt. Pinatubo eruption on the seasonal forecast performed with Goddard Earth Observing System (GEOS), an Earth System Model that includes a bulk aerosol model coupled to radiation and an interactive ocean.We performed three ten-member ensembles of 12-month simulations (June 1991-May 1992). These simulations were performed with a 0.5 longitude by 0.5 latitude horizontal resolution. Out of the three ensembles, one excludes the eruption, representing the control experiment. The other two ensembles include the eruption of Mt. Pinatubo, varying in the effective radius of the volcanic sulfate: one assumed the effective radius to be equal to 0.35 micron (as for tropospheric aerosol) and the other has the effective radius set to 0.6 micron (closer to natural observation). Through the analysis of the aerosol forcing results derived from the two assumptions, we can show how this forcing acts on the seasonal forecast system. In particular, we will focus on the impacts to the surface and ocean temperatures and precipitation.

Mukherjee, Nikita↗

An evaluation of experimental decadal predictions using CCSM4

This study assesses retrospective decadal prediction skill of sea surface temperature (SST) variability in initialized climate prediction experiments with the Community Climate System Model version 4 (CCSM4). Ensemble forecasts initialized with two different historical ocean and sea-ice states are evaluated and compared to an ensemble of uninitialized coupled simulations. Both experiments are subject to identical twentieth century historical radiative forcings. Each forecast consists of a 10-member ensemble integrated over a 10-year period. One set of historical ocean and sea-ice conditions used for initialization comes from a forced ocean-ice simulation driven by the Coordinated Ocean-ice Reference Experiments interannually varying atmospheric dataset. Following the Coordinated Model Intercomparison Project version 5 (CMIP5) protocol, these forecasts are initialized every 5 years from 1961 to 1996, and every year from 2000 to 2006. A second set of initial conditions comes from historical ocean state estimates obtained through the assimilation of in-situ temperature and salinity data into the CCSM4 ocean model. These forecasts are only available over a limited subset of the CMIP5 recommended start dates. Both methods result in retrospective SST prediction skill over broad regions of the Indian Ocean, western Pacific Ocean and North Atlantic Ocean that are significantly better than reference skill levels from a spatio-temporal auto-regressive statistical model of SST. However the subpolar gyre region of the North Atlantic stands out as the only region where the CCSM4 initialized predictions outperform uninitialized simulations. Some features of the ocean state estimates used for initialization and their impact on the forecasts are discussed.

Meteorology & Atmospheric Sciences↗

Emergent anthropogenic trends in California Current upwelling

Upwelling in the California Current System (CCS) sustains a productive ecosystem and is mediated by alongshore, equatorward wind stress. A decades-old hypothesis proposes that global warming will accelerate these upwelling favorable winds. Recent analyses provide empirical support for upwelling intensification in the poleward portion of the CCS. However, these studies rely on proxies for upwelling and are limited in their ability to distinguish anthropogenic forcing from internal climate variability. Here in this paper we estimate simulated changes in CCS upwelling from 1920 to 2100 using monthly output from a single climate model ensemble, where divergences among simulations can be attributed entirely to internal climate variability. Our projections suggest that CCS upwelling will become more intense in the spring and less intense in the summer as a result of anthropogenic climate change. Anthropogenic changes in upwelling will emerge primarily in the second half of the century.

54 ENVIRONMENTAL SCIENCES↗

Ettringite hysteresis under sorption from molecular simulations

Highlights: • Hybrid Grand Canonical Monte Carlo (GCMC) and Molecular Dynamics (MD) simulations are performed. • Anisotropic large strain coefficient of crystal shrinkage of ettringite is provided. • Entropy increases when ettringite converts into metaettringite. • Sorption-induced stress may lead to columns collapse upon drying. • H-bonds involving OH groups in columns play a role in hysteresis. The response of ettringite under sorption is critical for the utilization of this mineral as an energy storage material and to understand the concrete durability problems associated with ettringite formation. We report molecular simulations directly sampling the osmotic ensemble to understand the sorption processes and structural changes in ettringite-metaettringite transition. Desorption branch and sorption-induced volume changes from simulations agree with the experiments. The structural changes at low relative humidity are analyzed, revealing for the first time details of metaettringite structure. The reversibility under sorption observed experimentally is captured by simulation in the ettringite domain. The reversibility in the conversion of metaettringite into ettringite is not captured by the direct simulations in the osmotic ensemble due to the large volume changes associated with ettringite desorption. Finally, we discuss the role of hydrogen bonds on the hysteresis in ettringite.

36 MATERIALS SCIENCE↗

The Surface Energy Budget and Precipitation Efficiency for Convective Systems During TOGA, COARE, GATE, SCSMEX and ARM: Cloud-Resolving Model Simulations

A two-dimensional version of the Goddard Cumulus Ensemble (GCE) Model is used to simulate convective systems that developed in various geographic locations. Observed large-scale advective tendencies for potential temperature, water vapor mixing ratio, and horizontal momentum derived from field campaigns are used as the main forcing. By examining the surface energy budgets, the model results show that the two largest terms are net condensation (heating/drying) and imposed large-scale forcing (cooling/moistening) for tropical oceanic cases. These two terms arc opposite in sign, however. The contributions by net radiation and latent heat flux to the net condensation vary in these tropical cases, however. For cloud systems that developed over the South China Sea and eastern Atlantic, net radiation (cooling) accounts for about 20% or more of the net condensation. However, short-wave heating and long-wave cooling are in balance with each other for cloud systems over the West Pacific region such that the net radiation is very small. This is due to the thick anvil clouds simulated in the cloud systems over the Pacific region. Large-scale cooling exceeds large-scale moistening in the Pacific and Atlantic cases. For cloud systems over the South China Sea, however, there is more large-scale moistening than cooling even though the cloud systems developed in a very moist environment. though For three cloud systems that developed over a mid-latitude continent, the net radiation and sensible and latent heat fluxes play a much more important role. This means the accurate measurement of surface fluxes and radiation is crucial for simulating these mid-latitude cases.

Tao, W.-K.↗

Sampling the isothermal-isobaric ensemble by Langevin dynamics

We present a new method of conducting fully flexible-cell molecular dynamics simulation in isothermal-isobaric ensemble based on Langevin equations of motion. The stochastic coupling to all particle and cell degrees of freedoms is introduced in a correct way, in the sense that the stationary configurational distribution is proved to be consistent with that of the isothermal-isobaric ensemble. In order to apply the proposed method in computer simulations, a second order symmetric numerical integration scheme is developed by Trotter’s splitting of the single-step propagator. Moreover, a practical guide of choosing working parameters is suggested for user specified thermo- and baro-coupling time scales. The method and software implementation are carefully validated by a numerical example.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Strengthening of Ocean Heat Uptake Efficiency Associated with the Recent Climate Hiatus

The rate of increase of global-mean surface air temperature (SAT(sub g)) has apparently slowed during the last decade. We investigated the extent to which state-of-the-art general circulation models (GCMs) can capture this hiatus period by using multimodel ensembles of historical climate simulations. While the SAT(sub g) linear trend for the last decade is not captured by their ensemble means regardless of differences in model generation and external forcing, it is barely represented by an 11-member ensemble of a GCM, suggesting an internal origin of the hiatus associated with active heat uptake by the oceans. Besides, we found opposite changes in ocean heat uptake efficiency (k), weakening in models and strengthening in nature, which explain why the models tend to overestimate the SAT(sub g) trend. The weakening of k commonly found in GCMs seems to be an inevitable response of the climate system to global warming, suggesting the recovery from hiatus in coming decades.

Climate hiatus↗