Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ensemble simulations”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

Projecting exposure to extreme climate impact events across six event categories and three spatial scales

The extent and impact of climate‐related extreme events depend on the underlying meteorological, hydrological, or climatological drivers as well as on human factors such as land use or population density. Here we quantify the pure effect of historical and future climate change on the exposure of land and population to extreme climate impact events using an unprecedentedly large ensemble of harmonized climate impact simulations from the Inter‐Sectoral Impact Model Intercomparison Project phase 2b. Our results indicate that global warming has already more than doubled both the global land area and the global population annually exposed to all six categories of extreme events considered: river floods, tropical cyclones, crop failure, wildfires, droughts, and heatwaves. Global warming of 2°C relative to preindustrial conditions is projected to lead to a more than five‐fold increase in cross‐category aggregate exposure globally. Changes in exposure are unevenly distributed, with tropical and subtropical regions facing larger increases than higher latitudes. The largest increases in overall exposure are projected for the population of South Asia.

climate‐related extreme events↗
NASA NTRS2021

Contribution of the World's Main Dust Source Regions to the Global Cycle of Desert Dust

Even though desert dust is the most abundant aerosol by mass in Earth's atmosphere, the relative contributions of the world's major source regions to the global dust cycle remain poorly constrained. This problem hinders accounting for the potentially large impact of regional differences in dust properties on clouds, the Earth's energy balance, and terrestrial and marine biogeochemical cycles. Here, we constrain the contribution of each of the world's main dust source regions to the global dust cycle. We use an analytical framework that integrates an ensemble of global aerosol model simulations with observationally informed constraints on the dust size distribution, extinction efficiency, and regional dust aerosol optical depth (DAOD). We obtain a dataset that constrains the relative contribution of nine major source regions to size-resolved dust emission, atmospheric loading, DAOD, concentration, and deposition flux. We find that the 22–29 Tg (1 standard error range) global loading of dust with a geometric diameter up to 20 µm is partitioned as follows: North African source regions contribute ∼ 50 % (11–15 Tg), Asian source regions contribute ∼ 40 % (8–13 Tg), and North American and Southern Hemisphere regions contribute ∼ 10 % (1.8–3.2 Tg). These results suggest that current models on average overestimate the contribution of North African sources to atmospheric dust loading at ∼ 65 %, while underestimating the contribution of Asian dust at ∼ 30 %. Our results further show that each source region's dust loading peaks in local spring and summer, which is partially driven by increased dust lifetime in those seasons. We also quantify the dust deposition flux to the Amazon rainforest to be ∼ 10 Tg yr−1, which is a factor of 2–3 less than inferred from satellite data by previous work that likely overestimated dust deposition by underestimating the dust mass extinction efficiency. The data obtained in this paper can be used to obtain improved constraints on dust impacts on clouds, climate, biogeochemical cycles, and other parts of the Earth system.

Dust↗
NASA NTRS2022

Uncertainty Estimation for SMAP Level-1 Brightness Temperature Assimilation at Different Timescales

Soil Moisture Active Passive (SMAP) mission brightness temperature (T(b) ) observations are assimilated into NASA’s Catchment Land Surface Model using an Ensemble Kalman filter to update simulations of surface and root-zone soil moisture. Different time series components of the T(b) observations are assimilated including anomalies, inter-annual variations, and high frequency variations. To optimize the weights that the data assimilation (DA) puts on the observations, the ratio between the uncertainties of modeled and observed T(b) is approximated using modeled and observed soil moisture uncertainties estimated using triple collocation analysis. In a benchmark experiment, T(b) observations are assimilated using a spatially constant 4 Kelvin (K) observation uncertainty, as in the operational SMAP Level-4 algorithm. All DA experiments exhibit notable skill improvements in most regions. Improvements are largest for the inter-annual variations in the simulations of both surface and root-zone soil moisture (mean improvements in terms of Pearson correlation (-) are 0.08 and 0.06, respectively). Anomaly simulations improve similarly (0.07), and improvements in the high-frequency variations are only observed for surface soil moisture simulations (0.06). No notable difference in skill - neither improvement nor deterioration - is observed between the experiments that use optimized observation uncertainty parameters and the 4 K benchmark experiment. This may be explained by the presence of large observation operator errors, which are analytically shown to have the potential to render post-update uncertainty insensitive to inaccuracies in estimates of the Kalman gain. These results have important implications for the design of soil moisture DA systems, in particular for parameterizing model and observation uncertainties.

Hydrology↗
NASA NTRS2023

Impacts of Stratospheric Ozone Recovery on Southern Ocean Temperature and Heat Budget

The impacts of stratospheric ozone recovery on Southern Ocean surface and interior temperature, heat content, heat uptake, and heat transport are investigated by contrasting two ensemble chemistry-climate model simulations in 2005–2099: one with fixed ozone depleting substances (ODSs) and another with decreasing ODSs. In our simulations ozone recovery significantly affects Southern Ocean temperature, with large latitudinal and vertical variations. Ozone recovery causes a dipole change of the full-depth ocean heat content (OHC) with an increase south of 60°S and a decrease between 45°S and 60°S. Integrated over latitudes south of 40°S, OHC decreases in response to ozone recovery. This ocean heat loss is shown to be driven by weakened poleward ocean heat transport (OHT) across 40°S, which is partly canceled by enhanced heat uptake. The weakening of poleward OHT into the Southern Ocean is caused by the ozone-induced equatorward shift of the meridional overturning circulation.

stratospheric ozone↗
DOE OSTI2025

An Image-Plane Approach to Gravitational Lens Modeling of Interferometric Data

Strong gravitational lensing acts as a cosmic telescope, enabling the study of the high-redshift universe. Astronomical interferometers, such as the Atacama Large Millimeter/submillimeter Array (ALMA), have provided high-resolution images of strongly lensed sources at millimeter and submillimeter wavelengths. To model the mass and light distributions of lensing and source galaxies from strongly lensed images, strong lens modeling for interferometric observations is conventionally performed in the visibility space, which is computationally expensive. In this paper, we implement an image-plane lens modeling methodology for interferometric dirty images by accounting for noise correlations. We show that the image-plane likelihood function produces accurate model values when tested on simulated ALMA observations with an ensemble of noise realizations. We also apply our technique to ALMA observations of two sources selected from the South Pole Telescope survey, comparing our results with previous visibility-based models. Our model results are consistent with previous models for both parametric and pixelated source-plane reconstructions. We implement this methodology for interferometric lens modeling in the open-source software package lenstronomy.

Zhang, Nan [Illinois U., Urbana (main)] (ORCID:000↗
NASA NTRS2003

Evidence for Limited Indirect Aerosol Forcing in Stratocumulus

Increases in cloud cover and condensed water contribute more than half of the indirect aerosol effect in an ensemble of general circulation model (GCM) simulations estimating the global radiative forcing of anthropogenic aerosols. We use detailed simulations of marine stratocumulus clouds and airborne observations of ship tracks to show that increases in cloud cover and condensed water in reality are far less than represented by the GCM ensemble. Our results offer an explanation for recent simplified inverse climate calculations indicating that indirect aerosol effects are greatly exaggerated in GCMs.

Ackerman, Andrew S.↗
DOE OSTI2021

Extreme metrics from large ensembles: investigating the effects of ensemble size on their estimates

Abstract. We consider the problem of estimating the ensemble sizes required to characterize the forced component and the internal variability of a number of extreme metrics. While we exploit existing large ensembles, our perspective is that of a modeling center wanting to estimate a priori such sizes on the basis of an existing small ensemble (we assume the availability of only five members here). We therefore ask if such a small-size ensemble is sufficient to estimate accurately the population variance (i.e., the ensemble internal variability) and then apply a well-established formula that quantifies the expected error in the estimation of the population mean (i.e., the forced component) as a function of the sample size n, here taken to mean the ensemble size. We find that indeed we can anticipate errors in the estimation of the forced component for temperature and precipitation extremes as a function of n by plugging into the formula an estimate of the population variance derived on the basis of five members. For a range of spatial and temporal scales, forcing levels (we use simulations under Representative Concentration Pathway 8.5) and two models considered here as our proof of concept, it appears that an ensemble size of 20 or 25 members can provide estimates of the forced component for the extreme metrics considered that remain within small absolute and percentage errors. Additional members beyond 20 or 25 add only marginal precision to the estimate, and this remains true when statistical inference through extreme value analysis is used. We then ask about the ensemble size required to estimate the ensemble variance (a measure of internal variability) along the length of the simulation and – importantly – about the ensemble size required to detect significant changes in such variance along the simulation with increased external forcings. Using the F test, we find that estimates on the basis of only 5 or 10 ensemble members accurately represent the full ensemble variance even when the analysis is conducted at the grid-point scale. The detection of changes in the variance when comparing different times along the simulation, especially for the precipitation-based metrics, requires larger sizes but not larger than 15 or 20 members. While we recognize that there will always exist applications and metric definitions requiring larger statistical power and therefore ensemble sizes, our results suggest that for a wide range of analysis targets and scales an effective estimate of both forced component and internal variability can be achieved with sizes below 30 members. This invites consideration of the possibility of exploring additional sources of uncertainty, such as physics parameter settings, when designing ensemble simulations.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1994

Tropical impacts of SST forcing: A case study for 1987 versus 1988

The response of the NASA/Goddard Institute for Space Studies general circulation model (GCM) to large tropical sea surface temperature (SST) anomalies is investigated by evaluating model simulations of the particularly contrasting summer monsoon seasons 1987 and 1988. These years are representative of the warm and cold phases, respectively, of a recent El Nino-Southern Oscillation (ENSO) event. An ensemble averaging the results of three simulations was considered for each season, using monthly mean observed SST anomalies for June-August 1987 and 1988 as lower boundary forcing. Consistent with the European Center for Medium Weather Forecast (ECMWF)-analyzed winds, the simulations based on 1988 as compared to 1987 SST exhibit stronger upper-tropospheric irrotational circulation between the monsoon regions and the Southern Hemispheric sub-tropical anticyclones, a stronger Pacific Walker cell and a weaker subtropical westerly jet over the South Pacific. In the same vein, the modeled precipitation, indicating a more northerly position of the Pacific Inter-Tropical Convergence Zone (ITCZ) in 1988 compared with 1987, is supported by satellite observations of outgoing longwave radiation and highly reflective clouds.

Druyan, Leonard M.↗
NASA NTRS1999

Additional Contributions to the Development of the New Snow-Physics Scheme for SSiB

The Simplified Simple Biosphere Model (SSiB) had a well-documented problem with snowmelt timing and infiltration. A new snow-physics scheme was developed for use in SSIB. In this, the snow layer is separated from the soil, with its own energy budget and temperature. Solar energy reaching the top of the snowpack is divided into three parts: one, reflected by the snow; two, absorbed by the snow; and three, transmitted to the ground following a simple extinction relation. Heat is exchanged between the ground and snow by conduction and by radiation through an arbitrary air-gap between them. In the GSWP exercise using the GEWEX ISLSCP Initiative I forcing data (hereafter "offline"), it was found that the new snow scheme ameliorated a significant fraction of snowmelt time-delay as compared to observations from satellite. It also produced warmer ground temperatures under the snowpack, which allowed realistic meltwater infiltration, resulting in better simulated spring soil moisture recharge and peak runoff amount as compared to observations. An ensemble of six June-July-August (JJA) simulations for 1987 and 1988 were performed with the NASA Goddard GEOS II GCM coupled with the new snow-physics SSIB using new initial soil moisture (ISM) from the offline simulations. The GCM produced more realistic precipitation in northern regions that had large snowmelt and wetter ISM in response to better snow-physics, as compared to simulations with ISM without the new snow scheme. The new SSiB-GCM also increased the interannual precipitation signal in the Indian monsoon region, resulting from changes in ISM in the Himalayas and central Asia.

Mocko, David M.↗
DOE OSTI2020

The Role of Multiscale Interaction in Tropical Cyclogenesis and Its Predictability in Near-Global Aquaplanet Cloud-Resolving Simulations

Tropical cyclogenesis (TCG) is a multiscale process that involves interactions between large-scale circulation and small-scale convection. A near-global aquaplanet cloud-resolving model (NGAqua) with 4-km horizontal grid spacing that produces tropical cyclones (TCs) is used to investigate TCG and its predictability. This study analyzes an ensemble of three 20-day NGAqua simulations, with initial white-noise perturbations of low-level humidity. TCs develop spontaneously from the northern edge of the intertropical convergence zone (ITCZ), where large-scale flows and tropical convection provide necessary conditions for barotropic instability. Zonal bands of positive low-level absolute vorticity organize into cyclonic vortices, some of which develop into TCs. A new algorithm is developed to track the cyclonic vortices. A vortex-following framework analysis of the low-level vorticity budget shows that vertical stretching of absolute vorticity due to convective heating contributes positively to the vorticity spinup of the TCs. A case study and composite analyses suggest that sufficient humidity is key for convective development. TCG in these three NGAqua simulations undergoes the same series of interactions. The locations of cyclonic vortices are broadly predetermined by planetary-scale circulation and humidity patterns associated with ITCZ breakdown, which are predictable up to 10 days. Whether and when the cyclonic vortices become TCs depend on the somewhat more random feedback between convection and vorticity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Non-Kasha fluorescence of pyrene emerges from a dynamic equilibrium between excited states

Pyrene fluorescence after a high-energy electronic excitation exhibits a prominent band shoulder not present after excitation at low energies. The standard assignment of this shoulder as a non-Kasha emission from the second-excited state (S 2 ) has been recently questioned. To elucidate this issue, we simulated the fluorescence of pyrene using two different theoretical approaches based on vertical convolution and nonadiabatic dynamics with nuclear ensembles. To conduct the necessary nonadiabatic dynamics simulations with high-lying electronic states and deal with fluorescence timescales of about 100 ns of this large molecule, we developed new computational protocols. The results from both approaches confirm that the band shoulder is, in fact, due to S 2 emission. We show that the non-Kasha behavior is a dynamic-equilibrium effect not caused by a metastable S 2 minimum. However, it requires considerable vibrational energy, which can only be achieved in collisionless regimes after transitions into highly excited states. This strict condition explains why the S 2 emission was not observed in some experiments.

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

Uncertainty in Twenty-First-Century CMIP5 Sea Level Projections

The representative concentration pathway (RCP) simulations included in phase 5 of the Coupled Model Intercomparison Project (CMIP5) quantify the response of the climate system to different natural and anthropogenic forcing scenarios. These simulations differ because of 1) forcing, 2) the representation of the climate system in atmosphere-ocean general circulation models (AOGCMs), and 3) the presence of unforced (internal) variability. Global and local sea level rise projections derived from these simulations, and the emergence of distinct responses to the four RCPs depend on the relative magnitude of these sources of uncertainty at different lead times. Here, the uncertainty in CMIP5 projections of sea level is partitioned at global and local scales, using a 164-member ensemble of twenty-first-century simulations. Local projections at New York City (NYSL) are highlighted. The partition between model uncertainty, scenario uncertainty, and internal variability in global mean sea level (GMSL) is qualitatively consistent with that of surface air temperature, with model uncertainty dominant for most of the twenty-first century. Locally, model uncertainty is dominant through 2100, with maxima in the North Atlantic and the Arctic Ocean. The model spread is driven largely by 4 of the 16 AOGCMs in the ensemble; these models exhibit outlying behavior in all RCPs and in both GMSL and NYSL. The magnitude of internal variability varies widely by location and across models, leading to differences of several decades in the local emergence of RCPs. The AOGCM spread, and its sensitivity to model exclusion and/or weighting, has important implications for sea level assessments, especially if a local risk management approach is utilized.

Variability↗
DOE OSTI2023

Assessment of the Pan-Arctic Accelerated Rate of Sea Ice Decline in CMIP6 Historical Simulations

Abstract Arctic sea ice loss in response to a warming climate is assessed in 42 models participating in phase 6 of the Coupled Model Intercomparison Project (CMIP6). Sea ice observations show a significant acceleration in the rate of decline commencing near the turn of the twenty-first century. It is our assertion that state-of-the-art climate models should qualitatively reflect this accelerated trend within the limitations of internal variability and observational uncertainty. Our analysis shows that individual CMIP6 simulations of sea ice depict a wide range of model spread on biases and anomaly trends both across models and among their ensemble members. While the CMIP6 multimodel mean captures the observed sea ice area (SIA) decline relatively well, an individual model’s ability to represent the acceleration in sea ice decline remains a challenge. Seventeen (40%) out of 42 CMIP6 models and 37 (13%) out of the total 286 ensemble members reasonably capture the observed trends and acceleration in SIA decline. In addition, a larger ensemble size appears to increase the odds for a model to include at least one ensemble member skillfully representing the accelerated SIA trends. Simulations of sea ice volume (SIV) show much larger spread and uncertainty than SIA; however, due to limited availability of sea ice thickness data, these are not as well constrained by observations. Finally, we find that models with more ocean heat transport simulate larger sea ice declines, which suggests an emergent constraint in CMIP6 ensembles. This relationship points to the need for better understanding and modeling of ice–ocean interactions, especially with respect to frazil ice growth.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The Worldwide C3S CORDEX Grand Ensemble: A Major Contribution to Assess Regional Climate Change in the IPCC AR6 Atlas

The collaboration between the Coordinated Regional Climate Downscaling Experiment (CORDEX) and the Earth System Grid Federation (ESGF) provides open access to an unprecedented ensemble of regional climate model (RCM) simulations, across the 14 CORDEX continental-scale domains, with global coverage. These simulations have been used as a new line of evidence to assess regional climate projections in the latest contribution of the Working Group I (WGI) to the IPCC Sixth Assessment Report (AR6), particularly in the regional chapters and the Atlas. Here, we present the work done in the framework of the Copernicus Climate Change Service (C3S) to assemble a consistent worldwide CORDEX grand ensemble, aligned with the deadlines and activities of IPCC AR6. This work addressed the uneven and heterogeneous availability of CORDEX ESGF data by supporting publication in CORDEX domains with few archived simulations and performing quality control. It also addressed the lack of comprehensive documentation by compiling information from all contributing regional models, allowing for an informed use of data. In addition to presenting the worldwide CORDEX dataset, we assess here its consistency for precipitation and temperature by comparing climate change signals in regions with overlapping CORDEX domains, obtaining overall coincident regional climate change signals. The C3S CORDEX dataset has been used for the assessment of regional climate change in the IPCC AR6 (and for the interactive Atlas) and is available through the Copernicus Climate Data Store (CDS).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

An ensemble data assimilation modeling system for operational outdoor microalgae growth forecasting

Microalgae have received increasing attention as a potential feedstock for biofuel or biobased products. Forecasting the microalgae growth is beneficial for managers in planning pond operations and harvesting decisions. This study proposed a biomass forecasting system comprised of the Huesemann Algae Biomass Growth Model (BGM), the Modular Aquatic Simulation System in Two Dimensions (MASS2), ensemble data assimilation (DA), and numerical weather prediction Global Ensemble Forecast System (GEFS) ensemble meteorological forecasts. The novelty of this study is to seek the use of ensemble DA to improve both BGM and MASS2 model initial conditions with the assimilation of biomass and water temperature measurements and consequently improve short-term biomass forecasting skills. This study introduces the theory behind the proposed integrated biomass forecasting system, with an application undertaken in pseudo-real-time in three outdoor ponds cultured with Chlorella sorokiniana in Delhi, California, United States. Results from all three case studies demonstrate that the biomass forecasting system improved the short-term (i.e., 7-day) biomass forecasting skills by about 60% on average, comparing to forecasts without using the ensemble DA method. Given the satisfactory performances achieved in this study, it is probable that the integrated BGM-MASS2-DA forecasting system can be used operationally to inform managers in making pond operation and harvesting planning decisions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Climate hazard indices projections based on CORDEX-CORE, CMIP5 and CMIP6 ensemble

The CORDEX-CORE initiative was developed with the aim of producing homogeneous regional climate model (RCM) projections over domains world wide. In its first phase, two RCMs were run at 0.22° resolution downscaling 3 global climate models (GCMs) from the CMIP5 program for 9 CORDEX domains and two climate scenarios, the RCP2.6 and RCP8.5. The CORDEX-CORE simulations along with the CMIP5 GCM ensemble and the most recently produced CMIP6 GCM ensemble are analyzed, with focus on several temperature, heat, wet and dry hazard indicators for present day and mid-century and far future time slices. The CORDEX-CORE ensemble shows a better performance than the driving GCMs for several hazard indices due to its higher spatial resolution. For the far future time slice the 3 ensembles project an increase in all temperature and heat indices analyzed under the RCP8.5 scenario. The largest increases are always shown by the CMIP6 ensemble, except for Tx > 35 °C, for which the CORDEX-CORE projects higher warming. Extreme wet and flood prone maxima are projected to increase by the RCM ensemble over the la Plata basin in South America, the Congo basin in Africa, east North America, north east Europe, India and Indochina, regions where a better performance is obtained, whereas the GCM ensembles show small or negligible signals. Furthermore, compound hazard hotspots based on heat, drought and wet indicators are detected in each continent worldwide in region like Central America, the Amazon, the Mediterranean, South Africa and Australia, where a linear relation is shown between the heatwave and drought change signal, and region like Arabian peninsula, the central and south east Africa region (SEAF), the north west America (NWN), south east Asia, India, China and central and northern European regions (WCE, NEU) where the same linear relation is found for extreme precipitation and HW increases. Although still limited, the CORDEX-CORE initiative was able to produce high resolution climate projections with almost global coverage and can provide an important resource for impact assessment and climate service activities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A combined ensemble-volume average homogenization method for lattice structures with defects under dynamic and static loading

In the study of lattices structures, both experiments and numerical simulations are often conducted with small samples. Using combined ensemble and volume averaging, this work introduces a method to extract a macroscopic constitutive response of a lattice material from numerical simulations performed in periodic domains. The domain size needed to obtain statistically accurate results is investigated. Similar to molecular dynamics, the concept of the virial stress is introduced after homogenized equations are derived using the ensemble averaging method. Under static conditions, the virial stress is shown to agree with the volume averaged solid stress. Using the homogenization method, constitutive relations for this stress can be obtained from systems with uniform strains. Application of such obtained constitutive relations to more general cases results in an error proportional to the square of the ratio between the lattice length scale and the macroscopic length scale. Taking advantage of this property, numerical simulations are performed in systems with a uniform gradient of the average velocity. The volume average method is then used to accelerate convergence when studying lattices with defects. To avoid the artificial numerical time scale from the size of a representative volume element divided by the wave speed, a numerical scheme is developed to enforce a spatially uniform velocity gradient within the computational domain while allowing fluctuations of the velocity or displacement to develop naturally. To account for probability distribution of lattice defects, the stress is calculated as the ensemble-volume averaged value. For dynamic systems, energy dissipation properties are also studied.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modes of Variability in E3SM and CESM Large Ensembles

An adequate characterization of internal modes of climate variability (MoV) is a prerequisite for both accurate seasonal predictions and climate change detection and attribution. Assessing the fidelity of climate models in simulating MoV is therefore essential; however, doing so is complicated by the large intrinsic variations in MoV and the limited span of the observational record. Large ensembles (LEs) provide a unique opportunity to assess model fidelity in simulating MoV and quantify intermodel contrasts. Here, in this work, these goals are pursued in four recently produced LEs: the Energy Exascale Earth System Model (E3SM) versions 1 and 2 LEs, and the Community Earth System Model (CESM) versions 1 and 2 LEs. In general, the representation of global coupled modes is found to improve across successive E3SM and CESM versions in conjunction with the fidelity of the base state climate while the patterns of extratropical modes are well simulated across the ensembles. Various persistent shortcomings for all MoV are however identified and discussed. The results both demonstrate the successes of these recent model versions and suggest the potential for continued improvement in the representation of MoV with advances in model physics.

54 ENVIRONMENTAL SCIENCES↗