Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ensemble models”

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 757 records · Page 42

A Process-Model Perspective on Recent Changes in the Carbon Cycle of North America

Continental North America has been found to be a carbon (C) sink over recent decades by multiple studies employing a variety of estimation approaches. However, several key questions and uncertainties remain with these assessments. Here we used results from an ensemble of 19 state-of-the-art dynamic global vegetation models from the TRENDYv9 project to improve these estimates and study the drivers of its interannual variability. Our results show that North America has been a C sink with a magnitude of 0.37 ± 0.38 (mean and one standard deviation) PgC year−1 for the period 2000–2019 (0.31 and 0.44 PgC year−1 in each decade); split into 0.18 ± 0.12 PgC year−1 in Canada (0.15 and 0.20), 0.16 ± 0.17 in the United States (0.14 and 0.17), 0.02 ± 0.05 PgC year−1 in Mexico (0.02 and 0.02) and 0.01 ± 0.02 in Central America and the Caribbean (0.01 and 0.01). About 57% of the new C assimilated by terrestrial ecosystems is allocated into vegetation, 30% into soils, and 13% into litter. Losses of C due to fire account for 41% of the interannual variability of the mean net biome productivity for all North America in the model ensemble. Finally, we show that drought years (e.g., 2002) have the potential to shift the region to a small net C source in the simulations (−0.02 ± 0.46 PgC year−1). Our results highlight the importance of identifying the major drivers of the interannual variability of the continental-scale land C cycle along with the spatial distribution of local sink-source dynamics.

Guillermo Murray-Tortarolo↗

Leveraging AI and Spatial Data to Unlock Pipeline Integrity Insights: NETL’s Advanced Infrastructure Integrity Model (AIIM)

Maintaining the integrity of natural gas infrastructure plays a critical role in ensuring energy security. Robust, data-driven foundational AI models for pipeline integrity can help address risk management and mitigation issues. Trusted foundational models can help with industry adoption and accelerate innovation by enhancing integrity predictions, reduce costs, and informing infrastructure build-out. The AIIM dashboard was released in 2022 and utilizes multi-ML models for ensemble-type insights. It was expanded to include analytics on reported incidents. It was developed as an ESRI Dashboard to support data visualization & interrogation and contains pipeline data and model results.

Advanced Infrastructure Integrity Model (AIIM)↗

Turbofan noise generation. Volume 1: Analysis

Computer programs were developed which calculate the in-duct acoustic modes excited by a fan/stator stae operating at subsonic tip speed. Three noise source mechanisms are included: (1) sound generated by the rotor blades interacting with turbulence ingested into, or generated within, the inlet duct; (2) sound generated by the stator vanes interacting with the turbulent wakes of the rotors blades; and (3) sound generated by the stator vanes interacting with the mean velocity deficit wakes of the rotor blades. The fan/stator stage is modeled as an ensemble of blades and vanes of zero camber and thickness enclosed within an infinite hard-walled annular duct. Turbulence drawn into or generated within the inlet duct is modeled as nonhomogeneous and anisotropic random fluid motion, superimposed upon a uniform axial mean flow, and convected with that flow. Equations for the duct mode amplitudes, or expected values of the amplitudes, are derived.

Ventres, C. S.↗

Parameterizing Subgrid Variations of Land Surface Heat Fluxes to the Atmosphere Improves Boreal Summer Land Precipitation Simulation With the NCAR CESM1.2

Subgrid horizontal variations of land surface heat fluxes to the atmosphere resulting from subgrid land cover heterogeneity are important in land-atmosphere interaction in global climate models (GCMs). To incorporate it, a parameterization using stochastic sampling based on truncated normal distributions diagnosed from the land model and internal ensemble mean of multiple calls to the planetary boundary layer and deep convection schemes is developed. After its implementation in the National Center for Atmospheric Research (NCAR) Community Earth System Model version 1.2 (CESM1.2), large changes of simulated land precipitation are found to occur in regions where large subgrid variations of surface heat fluxes exist. The simulated precipitation is improved in boreal summer, especially over eastern China and the coastal areas of the Bay of Bengal. The improved precipitation is mainly a result of improved large-scale moisture convergence and advection by altered vertical diffusion and convection.

58 GEOSCIENCES↗

Stormtime ring current and radiation belt ion transport: Simulations and interpretations

We use a dynamical guiding-center model to investigate the stormtime transport of ring current and radiation-belt ions. We trace the motion of representative ions' guiding centers in response to model substorm-associated impulses in the convection electric field for a range of ion energies. Our simple magnetospheric model allows us to compare our numerical results quantitatively with analytical descriptions of particle transport, (e.g., with the quasilinear theory of radial diffusion). We find that 10-145-keV ions gain access to L approximately 3, where they can form the stormtime ring current, mainly from outside the (trapping) region in which particles execute closed drift paths. Conversely, the transport of higher-energy ions (approximately greater than 145 keV at L approximately 3) turns out to resemble radial diffusion. The quasilinear diffusion coefficient calculated for our model storm does not vary smoothly with particle energy, since our impulses occur at specific (although randomly determined) times. Despite the spectral irregularity, quasilinear theory provides a surprisingly accurate description of the transport process for approximately greater than 145-keV ions, even for the case of an individual storm. For 4 different realizations of our model storm, the geometric mean discrepancies between diffusion coefficients D(sup sim, sub LL) obtained from the simulations and the quasilinear diffusion coefficient D(sup ql, sub LL) amount to factors of 2.3, 2.3, 1.5, and 3.0, respectively. We have found that these discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) can be reduced slightly by invoking drift-resonance broadening to smooth out the sharp minima and maxima in D(sup ql, sub LL). The mean of the remaining discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) for the 4 different storms then amount to factors of 1.9, 2.1, 1.5, and 2.7, respectively. We find even better agreement when we reduce the impulse amplitudes systematically in a given model storm (e.g., reduction of all the impulse amplitudes by half reduces the discrepancy factor by at least its square root) and also when we average our results over an ensemble of 20 model storms (agreement is within a factor of 1.2 without impulse-amplitude reduction). We use our simulation results also to map phase-space densities f in accordance with Liouville's theorem. We find that the stormtime transport of approximately greater than 145-keV ions produces little change in f-bar the drift-averaged phase-space density on any drift shell of interest. However, the stormtime transport produces a major enhancement from the pre-storm phase-space density at energies approximately 30-145 keV, which are representative of the stormtime ring current.

Lyons, Larry R.↗

Future Global Mortality from Changes in Air Pollution Attributable to Climate Change

Ground-level ozone and fine particulate matter (PM (sub 2.5)) are associated with premature human mortality; their future concentrations depend on changes in emissions, which dominate the near-term, and on climate change. Previous global studies of the air-quality-related health effects of future climate change used single atmospheric models. However, in related studies, mortality results differ among models. Here we use an ensemble of global chemistry-climate models to show that premature mortality from changes in air pollution attributable to climate change, under the high greenhouse gas scenario RCP (Representative Concentration Pathway) 8.5, is probably positive. We estimate 3,340 (30,300 to 47,100) ozone-related deaths in 2030, relative to 2000 climate, and 43,600 (195,000 to 237,000) in 2100 (14 percent of the increase in global ozone-related mortality). For PM (sub 2.5), we estimate 55,600 (34,300 to 164,000) deaths in 2030 and 215,000 (76,100 to 595,000) in 2100 (countering by 16 percent the global decrease in PM (sub 2.5)-related mortality). Premature mortality attributable to climate change is estimated to be positive in all regions except Africa, and is greatest in India and East Asia. Most individual models yield increased mortality from climate change, but some yield decreases, suggesting caution in interpreting results from a single model. Climate change mitigation is likely to reduce air-pollution-related mortality.

Atmohspheric chemistry↗

Coordinating AgMIP Data and Models Across Global and Regional Scales for 1.5°C and 2.0°C Assessments

The Agricultural Model Intercomparison and Improvement Project (AgMIP) has developed novel methods for Coordinated Global and Regional Assessments (CGRA) of agriculture and food security in a changing world. The present study aims to perform a proof of concept of the CGRA to demonstrate advantages and challenges of the proposed framework. This effort responds to the request by the UN Framework Convention on Climate Change (UNFCCC) for the implications of limiting global temperature increases to 1.5°C and 2.0°C above pre-industrial conditions. The protocols for the 1.5°C/2.0°C assessment establish explicit and testable linkages across disciplines and scales, connecting outputs and inputs from the Shared Socio-economic Pathways (SSPs), Representative Agricultural Pathways (RAPs), Half a degree Additional warming, Prognosis and Projected Impacts (HAPPI) and Coupled Model Intercomparison Project Phase 5 (CMIP5) ensemble scenarios, global gridded crop models, global agricultural economics models, site-based crop models and within-country regional economics models. The CGRA consistently links disciplines, models and scales in order to track the complex chain of climate impacts and identify key vulnerabilities, feedbacks and uncertainties in managing future risk. CGRA proof-of-concept results show that, at the global scale, there are mixed areas of positive and negative simulated wheat and maize yield changes, with declines in some bread basket regions, at both 1.5°C and 2.0°C. Declines are especially evident in simulations that do not take into account direct CO2 effects on crops. These projected global yield changes mostly resulted in increases in prices and areas of wheat and maize in two global economics models. Regional simulations for 1.5°C and 2.0°C using site-based crop models had mixed results depending on the region and the crop. In conjunction with price changes from the global economics models, productivity declines in the Punjab, Pakistan, resulted in an increase in vulnerable households and the poverty rate. This article is part of the theme issue ‘The Paris Agreement: understanding the physical and social challenges for a warming world of 1.5°C above pre-industrial levels’.

interdisciplinary↗

Excessive Subsidence of the Southeast Atlantic Biomass Burning Aerosol Plume in the NASA Goddard Earth Observing System (GEOS) Model

Global models struggle to accurately represent the meteorological environment over the Southeast Atlantic Ocean, as evidenced by the poor representation of persistent marine stratocumulus decks off the west coast of southern Africa. This complex regime also experiences significant seasonal biomass burning aerosol plumes that are transported from southern Africa over the Atlantic Ocean, which gradually descend via large scale subsidence and interact radiatively and microphysically with the cloudy environment of the region. In comparisons to space borne lidar and aircraft observations we find that the NASA Goddard Earth Observing System (GEOS) model’s MERRA-2 aerosol and meteorology reanalysis shows that the aerosol plume descends too rapidly due to exaggerated large-scale downward vertical motion. This complicates the evaluation of aerosol properties utilizing aircraft observations collected throughout the ORACLES campaign as well as the model’s representation of the heating rate profiles due to the aerosol itself and its radiative interaction with the semi-permanent cloud deck typically located lower in altitude than the smoke. Two types of modeling experiments using GEOS will be presented to diagnose deficiencies in the model related to the excessive subsidence. Ensembles of free-running model simulations with and without radiatively active biomass burning aerosol will be used to quantify the impact of smoke on the atmospheric circulation. It will be demonstrated that heating within the aerosol layer results in self lofting through a reduction of the subsidence but did not alter the vertical placement or speed of the African Easterly jet. The amount of aerosol absorption required to eliminate the bias in the height of the aerosol plume through self-lofting will then be quantified. Finally, a second type of modeling experiment will be discussed in which the meteorology is constrained to observations. This will allow for biases in the height of the aerosol plume due to the general circulation model to be distinguished from biases that may be induced through the data assimilation process.

Allison Collow↗

Impact of the Earth’s Density Profile on Atmospheric Neutrino Oscillations

As atmospheric neutrinos traverse the Earth, the matter potential influences their oscillation probabilities in intricate ways, enhancing and suppressing the conversion of one neutrino’s flavor to another along the propagation. Comprehending their behavior within Earth’s complex density profile, primarily described by the Preliminary Reference Earth Model (PREM) [1], is essential for accurately describing oscillations and determining expected atmospheric neutrino event rates in the DUNE Far Detector. In this study, we address this challenge by considering an ensemble of Earth models constrained by astronomical measurements of the planet’s mass and moment of inertia. We evaluate how variations in densities and layer boundaries can affect oscillation probabilities and event rates, and what implications this could have on various physics analyses.

Ismerio Oliveira, Marcelo [Rio de Janeiro, Pont. U↗

PTA en route noise measurements

A long-range advanced turboprop en route noise database was obtained with weather, tracking, and onboard measurements. In-flight noise directivity measurements were made. Data repeatability within a test day was excellent. Day-to-day variability existed and is not completely understood and therefore not predicted. Comparison of a two-dimensional ray-tracing propagation model with the ensemble average ground-measured data was good; however, as stated above, the day-to-day data variability was not completely predicted. Future research will include looking at alternative propagation models. Three-dimensional ray tracing, fast field program, and the parabolic equation are possibilities. The effect of turbulence needs to be accessed.

Willshire, William L., Jr.↗

Probabilistic Mesomechanical Fatigue Model

A probabilistic mesomechanical fatigue life model is proposed to link the microstructural material heterogeneities to the statistical scatter in the macrostructural response. The macrostructure is modeled as an ensemble of microelements. Cracks nucleation within the microelements and grow from the microelements to final fracture. Variations of the microelement properties are defined using statistical parameters. A micromechanical slip band decohesion model is used to determine the crack nucleation life and size. A crack tip opening displacement model is used to determine the small crack growth life and size. Paris law is used to determine the long crack growth life. The models are combined in a Monte Carlo simulation to determine the statistical distribution of total fatigue life for the macrostructure. The modeled response is compared to trends in experimental observations from the literature.

Tryon, Robert G.↗

How well are hazards associated with derechos reproduced in regional climate simulations?

Abstract. A 15-member ensemble of convection-permitting regional simulations of the fast-moving and destructive derecho of 29–30 June 2012 that impacted the northeastern urban corridor of the USA is presented. This event generated 1100 reports of damaging winds, generated significant wind gusts over an extensive area of up to 500 000 km2, caused several fatalities, and resulted in widespread loss of electrical power. Extreme events such as this are increasingly being used within pseudo-global-warming experiments to examine the sensitivity of historical, societally important events to global climate non-stationarity and how they may evolve as a result of changing thermodynamic and dynamic contexts. As such it is important to examine the fidelity with which such events are described in hindcast experiments. The regional simulations presented herein are performed using the Weather Research and Forecasting (WRF) model. The resulting ensemble is used to explore simulation fidelity relative to observations for wind gust magnitudes, spatial scales of convection (as is manifest in high composite reflectivity, cREF), and both rainfall and hail production as a function of model configuration (microphysics parameterization, lateral boundary conditions (LBCs), start date, use of nudging, compiler choice, damping, and number of vertical levels). We also examine the degree to which each ensemble member differs with respect to key mesoscale drivers of convective systems (e.g., convective available potential energy and vertical wind shear) and critical manifestations of deep convection, e.g., vertical velocities, cold-pool generation, and how those properties relate to the correct characterization of the associated atmospheric hazards (wind gusts and hail). Use of a double-moment, seven-class scheme with number concentrations for all species (including hail and graupel) results in the greatest fidelity of model-simulated wind gusts and convective structure to the observations of this event. All ensemble members, however, fail to capture the intensity of the event in terms of the spatial extent of convection and the production of high near-surface wind gusts. We further show very high sensitivity to the LBCs employed and specifically that simulation fidelity is higher for simulations nested within ERA-Interim compared to ERA5. Excess convective available potential energy (CAPE) in all ensemble members after the derecho passage leads to excess production of convective cells, wind gusts, cREF > 40 dBZ, and precipitation during a frontal passage on the subsequent day. This event proved very challenging to forecast in real time and to reproduce in the 15-member hindcast simulation ensemble presented here. Future work could examine if simulations with other initial and lateral boundary conditions can achieve greater fidelity.

Shepherd, Tristan (ORCID:0000000186276419)↗

The Misrepresentation of the Southern African Easterly Jet in Models and Its Implications for Aerosol, Clouds, and Precipitation Distributions

The southern African easterly jet (AEJ-S) is an important midtropospheric feature critical to understanding the tropical convective system over central Africa and the aerosol–cloud interactions over the southeast Atlantic Ocean. However, it remains unclear how well models represent the AEJ-S and its influence on aerosol transport, clouds, and precipitation distribution. Here, in this work, we use ground- and satellite-based observations and reanalysis datasets to assess the representation of AEJ-S in the Coupled Model Intercomparison Project phase 6 (CMIP6) models between September and October during the peak of midtropospheric winds, aerosol transport, clouds, and precipitation. We find that most CMIP6 models have difficulty accurately simulating the strength, position, and spatial distribution of the AEJ-S. Specifically, the AEJ-S is relatively weaker and at a slightly lower altitude in the ensemble of CMIP6 models than represented by observation and reanalysis datasets. To assess the influence of the misrepresented the AEJ-S on CMIP6-simulated aerosol, clouds, and precipitation distributions, we performed composite analyses using models with low and high biases based on the estimates of their midtropospheric easterly wind speed. We find that the misrepresentation of the AEJ-S in CMIP6 models is associated with the overestimation of clouds and precipitation over central Africa, the underestimation of clouds over the southeast Atlantic Ocean, and the limitation of aerosol transport over the continent or the deviation of its spatial distribution from the typical zonal transport over the Atlantic Ocean. Because aerosols, clouds, and precipitation are important components of the regional climate system, we conclude that accurate representation of the AEJ-S is essential over central Africa and the southeast Atlantic Ocean.

54 ENVIRONMENTAL SCIENCES↗

Studying Turbulence Using Numerical Simulation Databases: Proceedings of the 1996 Summer Program - Part 6

Topics considered include: New approach to turbulence modeling; Second moment closure analysis of the backstep flow database; Prediction of the backflow and recovery regions in the backward facing step at various Reynolds numbers; Turbulent flame propagation in partially premixed flames; Ensemble averaged dynamic modeling. Also included a study of the turbulence structures of wall-bounded shear flows; Simulation and modeling of the elliptic streamline flow.

Source record↗

SMOS Soil Moisture Data Assimilation in the NASA Land Information System: Impact on LSM Initialization and NWP Forecasts

Land surface models are important components of numerical weather prediction (NWP) models, partitioning incoming energy into latent and sensitive heat fluxes that affect boundary layer growth and destabilization. During warm-season months, diurnal heating and convective initiation depend strongly on evapotranspiration and available boundary layer moisture, which are substantially affected by soil moisture content. Therefore, to properly simulate warm-season processes in NWP models, an accurate initialization of the land surface state is important for accurately depicting the exchange of heat and moisture between the surface and boundary layer. In this study, soil moisture retrievals from the Soil Moisture and Ocean Salinity (SMOS) satellite radiometer are assimilated into the Noah Land Surface Model via an Ensemble Kalman Filter embedded within the NASA Land Information System (LIS) software framework. The output from LIS-Noah is subsequently used to initialize runs of the Weather Research and Forecasting (WRF) NWP model. The impact of assimilating SMOS retrievals is assessed by initializing the WRF model with LIS-Noah output obtained with and without SMOS data assimilation. The southeastern United States is used as the domain for a preliminary case study. During the summer months, there is extensive irrigation in the lower Mississippi Valley for rice and other crops. The irrigation is not represented in the meteorological forcing used to drive the LIS-Noah integration, but the irrigated areas show up clearly in the SMOS soil moisture retrievals, resulting in a case with a large difference in initial soil moisture conditions. The impact of SMOS data assimilation on both Noah soil moisture fields and on short-term (0-48 hour) WRF weather forecasts will be presented.

Blankenship, Clay↗

An uncertain future change in aridity over the tropics

Abstract An ensemble of climate models from phase six of the Coupled Model Intercomparison Project shows that temperature and potential evapotranspiration (PET) are projected to increase globally towards the end of the 21st century. However, climate models show a spatially heterogeneous change in precipitation over the tropics. Consequently, future changes in aridity (a measure of water availability) are complex and location-dependent. We assess future changes in aridity using three climate models and several single-forcing experiments. Near-term (2021–2040) changes in aridity are small, and we focus instead on its long-term (2081–2100) changes. We show that the increase in greenhouse gases (GHG) primarily explains the spatial pattern, magnitude and ensemble spread of the long-term future changes in aridity. On this timescale, the effects of changes in emissions of anthropogenic aerosols are moderate compared to the effects of increases in atmospheric GHG concentrations. Model diversity in the responses to GHG concentration is large over northern Africa and North and South America. We suggest the large uncertainty is due to differences between models in simulating the effects of an increase in GHG concentrations on surface air temperature over the North Atlantic Ocean, on the interhemispheric temperature gradient, and on PET over North and South America.

Environmental Sciences & Ecology↗

Ensemble-function relationships to dissect mechanisms of enzyme catalysis

Decades of structure-function studies have established our current extensive understanding of enzymes. However, traditional structural models are snapshots of broader conformational ensembles of interchanging states. We demonstrate the need for conformational ensembles to understand function, using the enzyme ketosteroid isomerase (KSI) as an example. Comparison of prior KSI cryogenic x-ray structures suggested deleterious mutational effects from a misaligned oxyanion hole catalytic residue. However, ensemble information from room-temperature x-ray crystallography, combined with functional studies, excluded this model. Ensemble-function analyses can deconvolute effects from altering the probability of occupying a state (P-effects) and changing the reactivity of each state (k-effects); our ensemble-function analyses revealed functional effects arising from weakened oxyanion hole hydrogen bonding and substrate repositioning within the active site. Ensemble-function studies will have an integral role in understanding enzymes and in meeting the future goals of a predictive understanding of enzyme catalysis and engineering new enzymes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Models, data, and scripts associated with “Prediction of Distributed River Sediment Respiration Rates using Community-Generated Data and Machine Learning”

This data package is associated with the publication “Prediction of Distributed River Sediment Respiration Rates using Community-Generated Data and Machine Learning’’ submitted to the Journal of Geophysical Research: Machine Learning and Computation (Scheibe et al. 2024). River sediment respiration observations are expensive and labor intensive to obtain and there is no physical model for predicting this quantity. The Worldwide Hydrobiogeochemisty Observation Network for Dynamic River Systems (WHONDRS) observational data set (Goldman et al.; 2020) is used to train machine learning (ML) models to predict respiration rates at unsampled sites. This repository archives training data, ML models, predictions, and model evaluation results for the purposes of reproducibility of the results in the associated manuscript and community reuse of the ML models trained in this project. One of the key challenges in this work was to find an optimum configuration for machine learning models to work with this feature-rich (i.e. 100+ possible input variables) data set. Here, we used a two-tiered approach to managing the analysis of this complex data set: 1) a stacked ensemble of ML models that can automatically optimize hyperparameters to accelerate the process of model selection and tuning and 2) feature permutation importance to iteratively select the most important features (i.e. inputs) to the ML models. The major elements of this ML workflow are modular, portable, open, and cloud-based, thus making this implementation a potential template for other applications. This data package is associated with the GitHub repository found at Please see the file level metadata (flmd; “sl-archive-whondrs_flmd.csv”) for a list of all files contained in this data package and descriptions for each. Please see the data dictionary (dd; “sl-archive-whondrs_dd.csv”) for a list of all column headers contained within comma separated value (csv) files in this data package and descriptions for each. The GitHub repository is organized into five top-level directories: (1) “input_data” holds the training data for the ML models; (2) “ml_models” holds machine learning models trained on the data in “input_data”; (3) “scripts” contains data preprocessing and postprocessing scripts and intermediate results specific to this data set that bookend the ML workflow; (4) “examples” contains the visualization of the results in this repository including plotting scripts for the manuscript (e.g., model evaluation, FPI results) and scripts for running predictions with the ML models (i.e., reusing the trained ML models); (5) “output_data” holds the overall results of the ML model on that branch. Each trained ML model resides on its own branch in the repository; this means that inputs and outputs can be different branch-to-branch. Furthermore, depending on the number of features used to train the ML models, the preprocessing and postprocessing scripts, and their intermediate results, can also be different branch-to-branch. The “main-*” branches are meant to be starting points (i.e. trunks) for each model branch (i.e. sprouts). Please see the Branch Navigation section in the top-level README.md in the GitHub repository for more details. There is also one hidden directory “.github/workflows”. This hidden directory contains information for how to run the ML workflow as an end-to-end automated GitHub Action but it is not needed for reusing the ML models archived here. Please the top-level README.md in the GitHub repository for more details on the automation.

13C↗