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 217 records · Page 12

The MJO-QBO Relationship in a GCM with Stratospheric Nudging

Observational studies show a strong connection between the intraseasonal Madden-Julian oscillation (MJO) and the stratospheric quasi-biennial oscillation (QBO): the boreal winter MJO is stronger, more predictable, and has different teleconnections when the QBO in the lower stratosphere is easterly versus westerly. Despite the strength of the observed connection, global climate models do not produce an MJO-QBO link. Here the authors use a current-generation ocean-atmosphere coupled NASA Goddard Institute for Space Studies global climate model (Model E2.1) to examine the MJO-QBO link. To represent the QBO with minimal bias, the model zonal mean stratospheric zonal and meridional winds are relaxed to reanalysis fields from 1980-2017. The model troposphere, including the MJO, is allowed to freely evolve. The model with stratospheric nudging captures QBO signals well, including QBO temperature anomalies. However, an ensemble of nudged simulations still lacks an MJO-QBO connection.

Madden-Julian oscillation (MJO)↗

Role of forest regrowth in global carbon sink dynamics

Although the existence of a large carbon sink in terrestrial ecosystems is well-established, the drivers of this sink remain uncertain. It has been suggested that perturbations to forest demography caused by past land-use change, management, and natural disturbances may be causing a large component of current carbon uptake. Here we use a global compilation of forest age observations, combined with a terrestrial biosphere model with explicit modeling of forest regrowth, to partition the global forest carbon sink between old-growth and regrowth stands over the period 1981–2010. For 2001–2010 we find a carbon sink of 0.85 (0.66–0.96) Pg/year located in intact old-growth forest, primarily in the moist tropics and boreal Siberia, and 1.30 (1.03–1.96) Pg/year located in stands regrowing after past disturbance. Approaching half of the sink in regrowth stands would have occurred from demographic changes alone, in the absence of other environmental changes. These age-constrained results show consistency with those simulated using an ensemble of demographically-enabled terrestrial biosphere models following an independent reconstruction of historical land use and management. We estimate that forests will accumulate an additional 69 (44–131) Pg C in live biomass from changes in demography alone if natural disturbances, wood harvest, and reforestation continue at rates comparable to those during 1981–2010. Our results confirm that it is not possible to understand the current global terrestrial carbon sink without accounting for the sizeable sink due to forest demography. They also imply that a large portion of the current terrestrial carbon sink is strictly transient in nature.

forest↗

Comparison of Dust Optical Depth From Multi-Sensor Products and MONARCH (Multiscale Online Non-hydrostatic AtmospheRe CHemistry) Dust Reanalysis Over North Africa, the Middle East, and Europe

Aerosol reanalysis datasets are model-based, observationally constrained, continuous 3D aerosol fields with a relatively high temporal frequency that can be used to assess aerosol variations and trends, climate effects, and impacts on socioeconomic sectors, such as health. Here we compare and assess the recently published MONARCH (Multiscale Online Non-hydrostatic AtmospheRe CHemistry) high-resolution regional desert dust reanalysis over northern Africa, the Middle East, and Europe (NAMEE) with a combination of ground-based observations and space-based dust retrievals and products. In particular, we compare the total and coarse dust optical depth (DOD) from the new reanalysis with DOD products derived from MODIS (MODerate resolution Imaging Spectroradiometer), MISR (Multi-angle Imaging SpectroRadiometer), and IASI (Infrared Atmospheric Sounding Interferometer) spaceborne instruments. Despite the larger uncertainties, satellite-based datasets provide a better geographical coverage than ground-based observations, and the use of different retrievals and products allows at least partially overcoming some single-product weaknesses in the comparison. Nevertheless, limitations and uncertainties due to the type of sensor, its operating principle, its sensitivity, its temporal and spatial resolution, and the methodology for retrieving or further deriving dust products are factors that bias the reanalysis assessment. We, therefore, also use ground-based DOD observations provided by 238 stations of the AERONET (AErosol RObotic NETwork) located within the NAMEE region as a reference evaluation dataset. In particular, prior to the reanalysis assessment, the satellite datasets were evaluated against AERONET, showing moderate underestimations in the vicinities of dust sources and downwind regions, whereas small or significant overestimations, depending on the dataset, can be found in the remote regions. Taking these results into consideration, the MONARCH reanalysis assessment shows that total and coarse-DOD simulations are consistent with satellite- and ground-based data, qualitatively capturing the major dust sources in the area in addition to the dust transport patterns. Moreover, the MONARCH reanalysis reproduces the seasonal dust cycle, identifying the increased dust activity that occurred in the NAMEE region during spring and summer. The quantitative comparison between the MONARCH reanalysis DOD and satellite multi-sensor products shows that the reanalysis tends to slightly overestimate the desert dust that is emitted from the source regions and underestimate the transported dust over the outflow regions, implying that the model's removal of dust particles from the atmosphere, through deposition processes, is too effective. More specifically, small positive biases are found over the Sahara desert (0.04) and negative biases over the Atlantic Ocean and the Arabian Sea (−0.04), which constitute the main pathways of the long-range dust transport. Considering the DOD values recorded on average there, such discrepancies can be considered low, as the low relative bias in the Sahara desert (< 50 %) and over the adjacent maritime regions (< 100 %) certifies. Similarly, over areas with intense dust activity, the linear correlation coefficient between the MONARCH reanalysis simulations and the ensemble of the satellite products is significantly high for both total and coarse DOD, reaching 0.8 over the Middle East, the Atlantic Ocean, and the Arabian Sea and exceeding it over the African continent. Moreover, the low relative biases and high correlations are associated with regions for which large numbers of observations are available, thus allowing for robust reanalysis assessment.

Michail Mytilinaios↗

Global Nitrogen Needs to Improve Wheat Yield Under Climate Change

Increasing global food demand will require more food production without further exceeding the planetary boundaries, while at the same time adapting to climate change. We used an ensemble of wheat simulation models, with improved yield traits from the highest-yielding wheat genotypes to quantify potential yield gains and associated N requirements. This was explored for current and climate change scenarios across representative sites of major world wheat producing regions. The improved yield traits increased yield by 16% with current N fertilizer applications under both current climate and mid-century climate change scenarios. To achieve the full yield potential, a 52% increase in global average yield under a mid-century high warming climate scenario (RCP8.5), fertilizer use would need to increase fourfold over current use, which would unavoidably lead to higher environmental impacts from wheat production. Our results show the need to improve soil N availability and N use efficiency, along with yield potential.

Nitrogen↗

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

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↗

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

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↗

Effects of Forcing Differences and Initial Conditions on Inter-Model Agreement in the VolMIP Volc-Pinatubo-Full Experiment

This paper provides initial results from a multi-model ensemble analysis based on the volc-pinatubo-full experiment performed within the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP) as part of the sixth phase of the Coupled Model Intercomparison Project (CMIP6). The volc-pinatubo-full experiment is based on an ensemble of volcanic forcing-only climate simulations with the same volcanic aerosol dataset across the participating models (the 1991–1993 Pinatubo period from the CMIP6-GloSSAC dataset). The simulations are conducted within an idealized experimental design where initial states are sampled consistently across models from the CMIP6-piControl simulation providing unperturbed preindustrial background conditions. The multi-model ensemble includes output from an initial set of six participating Earth system models (CanESM5, GISS-E2.1-G, IPSL-CM6A-LR, MIROC-E2SL, MPI-ESM1.2-LR and UKESM1). The results show overall good agreement between the different models on the global and hemispheric scales concerning the surface climate responses, thus demonstrating the overall effectiveness of VolMIP's experimental design. However, small yet significant inter-model discrepancies are found in radiative fluxes, especially in the tropics, that preliminary analyses link with minor differences in forcing implementation; model physics, notably aerosol–radiation interactions; the simulation and sampling of El Niño–Southern Oscillation (ENSO); and, possibly, the simulation of climate feedbacks operating in the tropics. We discuss the volc-pinatubo-full protocol and highlight the advantages of volcanic forcing experiments defined within a carefully designed protocol with respect to emerging modelling approaches based on large ensemble transient simulations. We identify how the VolMIP strategy could be improved in future phases of the initiative to ensure a cleaner sampling protocol with greater focus on the evolving state of ENSO in the pre-eruption period.

metadata↗

Global Modeling Activities and NAME

In this talk I will review global modeling activities in the United States that could contribute to and benefit from NAME activities. I will present some preliminary results from several global atmospheric general circulation model simulation experiments for the initial NAME model intercomparison project period of May-Oct 1990. These include an ensemble of medium resolution simulations, and a high resolution (one half degree) simulation. I will also discuss possible high resolution global data assimilation experiments that could be used to help validate the model simulations and assimilate planned NAME observations.

Schubert, Siegfried↗

Global Modeling Activities and NAME

In this talk I will review global modeling activities in the United States that could contribute to and benefit from NAME activities. I will present some preliminary results from several global atmospheric general circulation model simulation experiments for the initial NAME model intercomparison project period of May-Oct 1990. These include an ensemble of medium resolution simulations, and a high resolution (one half degree) simulation. I will also discuss possible high resolution global data assimilation experiments that could be used to help validate the model simulations and assimilate planned NAME observations.

Schubert, Siegfried↗

Dataset for manuscript "Equipartition and the temperature of maximum density of TIP4P/2005 water"

We simulate TIP4P/2005 water in the temperature range of 257 K to 318 K with time-steps 0.25, 0.50, 1.00, 2.00, and 4.00 fs. The density-temperature behavior obtained using 0.25 or 0.50 fs are in excellent agreement with each other but differ from those obtained using time-steps that have been shown earlier to lead to a breakdown of equipartition. The temperature of maximum density (TMD) is 277.15 K with time-step 0.25 or 0.50 fs, but is shifted to progressively lower values for longer time-steps, a trend that holds for different thermostat/barostat combinations. Enhancing the water-water dispersion interaction, as has been recommended for simulating disordered proteins in TIP4P/2005, degrades the description of the liquid-vapor phase envelope. We present a simple physically transparent reasoning to highlight the separation of the time-scales between translational and rotational motion. We also develop a metric, Chi, that we term the equipartition anomaly, to detect equipartition violations in simulations that include molecules that are treated as rigid objects. Calculating Chi is shown to be straightforward and sensitive to equipartition violations. A key takeaway from this study is that using sufficiently short time-steps (less than or equal to 0.5 fs) to preserve equipartition is essential for obtaining meaningful liquid water properties and for producing reliable simulation data, as correct-ensemble sampling is fundamental to ensure reproducibility across codes and simulation alogrithms. The included dataset provides the raw data used in the preparation of the graphs noted in the manuscript.

36 MATERIALS SCIENCE↗

CERESMIP: A Climate Modeling Protocol to Investigate Recent Trends in the Earth's Energy Imbalance

The Clouds and the Earth's Radiant Energy System (CERES) project has now produced over two decades of observed data on the Earth's Energy Imbalance (EEI) and has revealed substantive trends in both the reflected shortwave and outgoing longwave top-of-atmosphere radiation components. Available climate model simulations suggest that these trends are incompatible with purely internal variability, but that the full magnitude and breakdown of the trends are outside of the model ranges. Unfortunately, the Coupled Model Intercomparison Project (Phase 6) (CMIP6) protocol only uses observed forcings to 2014 (and Shared Socioeconomic Pathways (SSP) projections thereafter), and furthermore, many of the ‘observed' drivers have been updated substantially since the CMIP6 inputs were defined. Most notably, the sea surface temperature (SST) estimates have been revised and now show up to 50% greater trends since 1979, particularly in the southern hemisphere. Additionally, estimates of short-lived aerosol and gas-phase emissions have been substantially updated. These revisions will likely have material impacts on the model-simulated EEI. We therefore propose a new, relatively low-cost, model intercomparison, CERESMIP, that would target the CERES period (2000-present), with updated forcings to at least the end of 2021. The focus will be on atmosphere-only simulations, using updated SST, forcings and emissions from 1990 to 2021. The key metrics of interest will be the EEI and atmospheric feedbacks, and so the analysis will benefit from output from satellite cloud observation simulators. The Tier 1 request would consist only of an ensemble of AMIP-style simulations, while the Tier 2 request would encompass uncertainties in the applied forcing, atmospheric composition, single and all-but-one forcing responses. We present some preliminary results and invite participation from a wide group of models.

CMIP6↗

RASPA3

RASPA3, a molecular simulation code for computing adsorption and diffusion in nanoporous materials and thermodynamic and transport properties of fluids. It implements force field based classical Monte Carlo/molecular dynamics in various ensembles. RASPA3 is rewritten from the ground up in C++23 with speed and code readability in mind. Transition-matrix Monte Carlo is added to compute the density of states and free energies. The Monte Carlo code for rigid molecules is based on quaternions, and the atomic positions needed in the energy evaluation are recreated from the center of mass position and quaternion orientation. The expanded ensemble methodology for fractional molecules, with a scaling parameter λ between 0 and 1, now also keeps track of analytic expressions of dU/dλ, allowing independent verification of the chemical potential using thermodynamic integration. The source code is freely available under the MIT license on GitHub.

Dubbeldam, David↗

Time-Varying Changes in the Simulated Structure of the Brewer-Dobson Circulation

A series of simulations using the NASA Goddard Earth Observing System Chemistry Climate Model are analyzed in order to assess changes in the Brewer-Dobson Circulation (BDC) over the past 55 years. When trends are computed over the past 55 years, the BDC accelerates throughout the stratosphere, consistent with previous modeling results. However, over the second half of the simulations (i.e., since the late 1980s), the model simulates structural changes in the BDC as the temporal evolution of the BDC varies between regions in the stratosphere. In the mid-stratosphere in the midlatitude Northern Hemisphere, the BDC does not accelerate in the ensemble mean of our simulations despite increases in greenhouse gas concentrations and warming sea surface temperatures, and it even decelerates in one ensemble member. This deceleration is reminiscent of changes inferred from satellite instruments and in situ measurements. In contrast, the BDC in the lower stratosphere continues to accelerate. The main forcing agents for the recent slowdown in the mid-stratosphere appear to be declining ozone-depleting substance (ODS) concentrations and the timing of volcanic eruptions. Changes in both mean age of air and the tropical upwelling of the residual circulation indicate a lack of recent acceleration. We therefore clarify that the statement that is often made that climate models simulate a decreasing age throughout the stratosphere only applies over long time periods and is not necessarily the case for the past 25 years, when most tracer measurements were taken.

stratospheric ozone↗

The Pliocene Model Intercomparison Project Phase 2: Large-scale Climate Features and Climate Sensitivity

The Pliocene epoch has great potential to improve our understanding of the long-term climatic and environmental consequences of an atmospheric CO2 concentration near ∼400 parts per million by volume. Here we present the large-scale features of Pliocene climate as simulated by a new ensemble of climate models of varying complexity and spatial resolution based on new reconstructions of boundary conditions (the Pliocene Model Intercomparison Project Phase 2; PlioMIP2). As a global annual average, modelled surface air temperatures increase by between 1.7 and 5.2 ∘C relative to the pre-industrial era with a multi-model mean value of 3.2 ∘C. Annual mean total precipitation rates increase by 7 % (range: 2 %–13 %). On average, surface air temperature (SAT) increases by 4.3 ∘C over land and 2.8 ∘C over the oceans. There is a clear pattern of polar amplification with warming polewards of 60∘ N and 60∘ S exceeding the global mean warming by a factor of 2.3. In the Atlantic and Pacific oceans, meridional temperature gradients are reduced, while tropical zonal gradients remain largely unchanged. There is a statistically significant relationship between a model's climate response associated with a doubling in CO2 (equilibrium climate sensitivity; ECS) and its simulated Pliocene surface temperature response. The mean ensemble Earth system response to a doubling of CO2 (including ice sheet feedbacks) is 67 % greater than ECS; this is larger than the increase of 47 % obtained from the PlioMIP1 ensemble. Proxy-derived estimates of Pliocene sea surface temperatures are used to assess model estimates of ECS and give an ECS range of 2.6–4.8 ∘C. This result is in general accord with the ECS range presented by previous Intergovernmental Panel on Climate Change (IPCC) Assessment Reports.

Pliocene Model Intercomparison Project↗

"Intelligent Ensemble" Projections of Precipitation and Surface Radiation in Support of Agricultural Climate Change Adaptation

Earth's climate is changing and will continue to change into the foreseeable future. Expected changes in the climatological distribution of precipitation, surface temperature, and surface solar radiation will significantly impact agriculture. Adaptation strategies are, therefore, required to reduce the agricultural impacts of climate change. Climate change projections of precipitation, surface temperature, and surface solar radiation distributions are necessary input for adaption planning studies. These projections are conventionally constructed from an ensemble of climate model simulations (e.g., the Coupled Model Intercomparison Project 5 (CMIP5)) as an equal weighted average, one model one vote. Each climate model, however, represents the array of climate-relevant physical processes with varying degrees of fidelity influencing the projection of individual climate variables differently. Presented here is a new approach, termed the "Intelligent Ensemble, that constructs climate variable projections by weighting each model according to its ability to represent key physical processes, e.g., precipitation probability distribution. This approach provides added value over the equal weighted average method. Physical process metrics applied in the "Intelligent Ensemble" method are created using a combination of NASA and NOAA satellite and surface-based cloud, radiation, temperature, and precipitation data sets. The "Intelligent Ensemble" method is applied to the RCP4.5 and RCP8.5 anthropogenic climate forcing simulations within the CMIP5 archive to develop a set of climate change scenarios for precipitation, temperature, and surface solar radiation in each USDA Farm Resource Region for use in climate change adaptation studies.

Taylor, Patrick C.↗

Huge ensembles – Part 1: Design of ensemble weather forecasts using spherical Fourier neural operators

Abstract. Simulating low-likelihood high-impact extreme weather events in a warming world is a significant and challenging task for current ensemble forecasting systems. While these systems presently use up to 100 members, larger ensembles could enrich the sampling of internal variability. They may capture the long tails associated with climate hazards better than traditional ensemble sizes. Due to computational constraints, it is infeasible to generate huge ensembles (comprised of 1000–10 000 members) with traditional, physics-based numerical models. In this two-part paper, we replace traditional numerical simulations with machine learning (ML) to generate hindcasts of huge ensembles. In Part 1, we construct an ensemble weather forecasting system based on spherical Fourier neural operators (SFNOs), and we discuss important design decisions for constructing such an ensemble. The ensemble represents model uncertainty through perturbed-parameter techniques, and it represents initial condition uncertainty through bred vectors, which sample the fastest-growing modes of the forecast. Using the European Centre for Medium-Range Weather Forecasts Integrated Forecasting System (IFS) as a baseline, we develop an evaluation pipeline composed of mean, spectral, and extreme diagnostics. With large-scale, distributed SFNOs with 1.1 billion learned parameters, we achieve calibrated probabilistic forecasts. As the trajectories of the individual members diverge, the ML ensemble mean spectra degrade with lead time, consistent with physical expectations. However, the individual ensemble members' spectra stay constant with lead time. Therefore, these members simulate realistic weather states during the rollout, and the ML ensemble passes a crucial spectral test in the literature. The IFS and ML ensembles have similar extreme forecast indices, and we show that the ML extreme weather forecasts are reliable and discriminating. These diagnostics ensure that the ensemble can reliably simulate the time evolution of the atmosphere, including low-likelihood high-impact extremes. In Part 2, we generate a huge ensemble initialized each day in summer 2023, and we characterize the simulations of extremes.

Mahesh, Ankur↗