Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “global 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 325 records · Page 18

LANL Institutional Computing Highlight: Local Time-Stepping in Global to Regional Ocean Modeling (Annual Report) [Slides]

Multi-resolution global ocean simulations are possible with the Model for Prediction Across Scales-Ocean (MPAS-Ocean), although the model is currently constrained to a uniform time step. The time step is determined by the size of the smallest cell in the grid, according to the CFL condition. To overcome this issue, a local time stepping (LTS) scheme has been designed for MPAS-Ocean. We have made progress on implementing the LTS scheme, as well as quantitative results and performance tests.

58 GEOSCIENCES↗

The effect of geographic sampling on evaluation of extreme precipitation in high-resolution climate models

Traditional approaches for comparing global climate models and observational data products typically fail to account for the geographic location of the underlying weather station data. For modern global high-resolution models with a horizontal resolution of tens of kilometers, this is an oversight since there are likely grid cells where the physical output of a climate model is compared with a statistically interpolated quantity instead of actual measurements of the climate system. In this paper, we quantify the impact of geographic sampling on the relative performance of high-resolution climate model representations of precipitation extremes in boreal winter (December–January–February) over the contiguous United States (CONUS), comparing model output from five early submissions to the HighResMIP subproject of the CMIP6 experiment. We find that properly accounting for the geographic sampling of weather stations can significantly change the assessment of model performance. Across the models considered, failing to account for sampling impacts the different metrics (extreme bias, spatial pattern correlation, and spatial variability) in different ways (both increasing and decreasing). We argue that the geographic sampling of weather stations should be accounted for in order to yield a more straightforward and appropriate comparison between models and observational data sets, particularly for high-resolution models with a horizontal resolution of tens of kilometers. While we focus on the CONUS in this paper, our results have important implications for other global land regions where the sampling problem is more severe.

54 ENVIRONMENTAL SCIENCES↗

Modeling noise in global Mølmer-Sørensen interactions applied to quantum approximate optimization

Many-qubit Mølmer-Sørensen (MS) interactions applied to trapped ions offer unique capabilities for quantum information processing, with applications including quantum simulation and the quantum approximate optimization algorithm (QAOA). Here, we develop a physical model to describe many-qubit MS interactions under four sources of experimental noise: vibrational mode frequency fluctuations, laser power fluctuations, thermal initial vibrational states, and state preparation and measurement errors. The model parametrizes these errors from simple experimental measurements, without free parameters. Herein, we validate the model in comparison with experiments that implement sequences of MS interactions on two 171 Yb + ions. The model shows reasonable agreement after several MS interactions as quantified by the reduced chi-squared statistic $χ^{2}_{red}$≈2. As an application we examine MaxCut QAOA experiments on three and six ions. The experimental performance is quantified by approximation ratios that are 91% and 83% of the optimal theoretical values. Our model predicts $0.93^{+0.03}_{-0.02}$ and $0.95^{+0.04}_{-0.03}$, respectively, with disagreement in the latter value attributable to secondary noise sources beyond those considered in our analysis. With realistic experimental improvements to reduce measurement error and radial trap frequency variations, the model achieves approximation ratios that are 99% of the optimal. Incorporating these improvements into future experiments is expected to reveal new aspects of noise for future modeling and experimental improvements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Increasing Functional Diversity in a Global Land Surface Model Illustrates Uncertainties Related to Parameter Simplification

Abstract Simulations of the land surface carbon cycle typically compress functional diversity into a small set of plant functional types (PFT), with parameters defined by the average value of measurements of functional traits. In most earth system models, all wild plant life is represented by between five and 14 PFTs and a typical grid cell (≈100 × 100 km) may contain a single PFT. Model logic applied to this coarse representation of ecological functional diversity provides a reasonable proxy for the carbon cycle, but does not capture the non‐linear influence of functional traits on productivity. Here we show through simulations using the Energy Exascale Land Surface Model in 15 diverse terrestrial landscapes, that better accounting for functional diversity markedly alters predicted total carbon uptake. The shift in carbon uptake is as great as 30% and 10% in boreal and tropical regions, respectively, when compared to a single PFT parameterized with the trait means. The traits that best predict gross primary production vary based on vegetation phenology, which broadly determines where traits fall within the global distribution. Carbon uptake is more closely associated with specific leaf area for evergreen PFTs and the leaf carbon to nitrogen ratio in deciduous PFTs.

58 GEOSCIENCES↗

Aboveground biomass density models for NASA’s Global Ecosystem Dynamics Investigation (GEDI) lidar mission

NASA's Global Ecosystem Dynamics Investigation (GEDI) is collecting spaceborne full waveform lidar data with a primary science goal of producing accurate estimates of forest aboveground biomass density (AGBD). This paper presents the development of the models used to create GEDI's footprint-level (~25 m) AGBD (GEDI04_A) product, including a description of the datasets used and the procedure for final model selection. The data used to fit our models are from a compilation of globally distributed spatially and temporally coincident field and airborne lidar datasets, whereby we simulated GEDI-like waveforms from airborne lidar to build a calibration database. We used this database to expand the geographic extent of past waveform lidar studies, and divided the globe into four broad strata by Plant Functional Type (PFT) and six geographic regions. GEDI's waveform-to-biomass models take the form of parametric Ordinary Least Squares (OLS) models with simulated Relative Height (RH) metrics as predictor variables. From an exhaustive set of candidate models, we selected the best input predictor variables, and data transformations for each geographic stratum in the GEDI domain to produce a set of comprehensive predictive footprint-level models. We found that model selection frequently favored combinations of RH metrics at the 98th, 90th, 50th, and 10th height above ground-level percentiles (RH98, RH90, RH50, and RH10, respectively), but that inclusion of lower RH metrics (e.g. RH10) did not markedly improve model performance. Second, forced inclusion of RH98 in all models was important and did not degrade model performance, and the best performing models were parsimonious, typically having only 1-3 predictors. Third, stratification by geographic domain (PFT, geographic region) improved model performance in comparison to global models without stratification. Fourth, for the vast majority of strata, the best performing models were fit using square root transformation of field AGBD and/or height metrics. There was considerable variability in model performance across geographic strata, and areas with sparse training data and/or high AGBD values had the poorest performance. These models are used to produce global predictions of AGBD, but will be improved in the future as more and better training data become available.

54 ENVIRONMENTAL SCIENCES↗

Experiments and modeling to evaluate global reaction kinetics of three-way catalyst light off for net-zero carbon fuels and selected fuel blends

Internal combustion engine exhaust composition can have a significant impact on the light-off characteristics of CO, NO, and organics in the three-way-catalyst during cold-start. The U.S. Department of Energy’s Co-Optimization of Fuels and Engines initiative was aimed at developing new high-performance fuels, which raises the question of how current three-way-catalysts will perform if the exhaust composition is significantly changed. This work is aimed at quantifying the performance of the catalyst under a variety of different fuel components and mixtures. To this end, a reduced-order model is developed wherein the set of organics examined are placed into molecular groups based on their reactivities and to model each organic within that group with a common kinetic expression. In conclusion, this strategy allows for a first approximation to the light-off characteristics of all organics that fall within that group and to make predictions for fuel blends.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Decarbonization Scenarios in the United States: Comparing Biofuels Growth in Two Models - GCAM and BSM

Scenarios for deep decarbonization rely on biomass for biofuels, biopower, and bioproducts, often including negative emissions via carbon capture and storage or utilization. Despite the prominence of biomass in many deep decarbonization pathways, critical questions remain about biomass allocation, effects of transportation electrification, the pace of growth, and implications for agriculture and land use. We address these questions through a unique comparison of carbon pricing effects on the growth of biomass utilization and its effects on land use in the United States by comparing results from a multisectoral integrated assessment model, the Global Change Analysis Model [GCAM], with results from a biomass-to-biofuels system dynamics model, the Biomass Scenario Model [BSM]. We contribute to model comparison efforts by analyzing the biomass deployment needed for a scenario consistent with a "Middle of the Road" Shared Socioeconomic Pathway [SSP2] and a representative concentration pathway of 2.6 W/m2. The GCAM scenarios solve for global equilibrium conditions that are consistent with this pathway, including demands for biomass across all economic sectors and representing bioenergy with carbon capture and storage as a technology option. The BSM scenarios assess those biomass and biofuel results for the United States and identify challenges associated with that pace and amount of expansion. In the scenario analysis, we harmonize key factors such as carbon price trajectory, domestic ethanol fuel demand, ethanol blending, and arable land availability, and vary them in both models. In GCAM, we vary the carbon price, transportation electrification, ethanol blending constraints, and arable land availability inputs and the value of the carbon in land; in BSM, in addition to directly inputting certain GCAM results, we vary the maximum rate of biorefinery construction, flexibility of feedstock types across conversion processes, and policy incentives such as tax credits and renewable identification number payments. The selected carbon price trajectory results in a rapid increase in biofuel production in the United States, reaching about 9.4 EJ/year in 2060 in the highest scenario analyzed in GCAM. Results differ between the two models in timing and ultimate quantity of biomass and biofuel production. GCAM biofuel quantities generally exceed BSM amounts because CCS is applied to biofuel pathways in GCAM, and because of differences in capacity expansion and related dynamics of land allocation, biomass production, and price dynamics. These dynamics include rapid biorefinery capacity expansion in high demand cases. To satisfy this biomass demand, GCAM rapidly equilibrates land allocation, but the BSM limits the rate at which this re-allocation can occur. A further contrast with the equilibrium approach in GCAM is that the BSM represents a delay between planting and harvesting woody biomass resources. As a result of these model contrasts, feedstock costs in BSM increase more than in GCAM, and the absence of CCS in the BSM also reduces the relative economic attractiveness of biofuels production. The bottlenecks, lags, and price increases also lead to potential for volatility in feedstock price and land allocation to biomass in the BSM. GCAM has more biomass production than BSM in all scenarios, partly because of the broader, economy-wide coverage of GCAM, in contrast to BSM's exclusive focus on biofuels. In both models, trends like those of biofuels production were observed for biomass production: minimal growth without a carbon price and policy incentives, and increases with a carbon price, particularly with carbon capture and storage, because the inputs assume that biopower and biofuels decrease greenhouse gas emissions. In high policy scenarios, biomass demand is high, and the consequent high biomass prices due to the land re-allocation bottleneck in the BSM limit biofuel production even if the biorefinery capacity is expanded. However, because biomass prices do not increase as much in the low policy scenario, growth is slower and the land-reallocation bottleneck no longer dominates, such that the effect of increased capacity can be seen. Across both the models, a change in assumptions from less to more land availability increases biofuel production in both GCAM and BSM, as the upward pressure on feedstock price and volatility are both reduced.

biofuels↗

Improving Stratocumulus Cloud Amounts in a 200‐m Resolution Multi‐Scale Modeling Framework Through Tuning of Its Interior Physics

Abstract High‐Resolution Multi‐scale Modeling Frameworks (HR)—global climate models that embed separate, convection‐resolving models with high enough resolution to resolve boundary layer eddies—have exciting potential for investigating low cloud feedback dynamics due to reduced parameterization and ability for multidecadal throughput on modern computing hardware. However low clouds in past HR have suffered a stubborn problem of over‐entrainment due to an uncontrolled source of mixing across the marine subtropical inversion manifesting as stratocumulus dim biases in present‐day climate, limiting their scientific utility. We report new results showing that this over‐entrainment can be partly offset by using hyperviscosity and cloud droplet sedimentation. Hyperviscosity damps small‐scale momentum fluctuations associated with the formulation of the momentum solver of the embedded large eddy simulation. By considering the sedimentation process adjacent to default one‐moment microphysics in HR, condensed phase particles can be removed from the entrainment zone, which further reduces entrainment efficiency. The result is an HR that can produce more low clouds with a higher liquid water path and a reduced stratocumulus dim bias. Associated improvements in the explicitly simulated sub‐cloud eddy spectrum are observed. We report these sensitivities in multi‐week tests and then explore their operational potential alongside microphysical retuning in decadal simulations at operational 1.5° exterior resolution. The result is a new HR having desired improvements in the baseline present‐day low cloud climatology, and a reduced global mean bias and root mean squared error of absorbed shortwave radiation. We suggest it should be promising for examining low cloud feedbacks with minimal approximation.

54 ENVIRONMENTAL SCIENCES↗

Evaluating and Improving Convective Parameterization for GCMs Using ARM Observations. Final report

This technical report summarizes the achievements during the project period funded by the Atmospheric System Research (ASR) program. The general goal of the project is to improve the representation of atmospheric deep convection in global climate models. By using DOE Atmospheric Radiation Measurement (ARM) program observations and a hierarchy of models from cloud-resolving to global climate models, we evaluated and improved many aspects of convection parameterization schemes, such as the trigger function for convection and closure condition to determine the amount of convection. We also incorporated the stochasticity of convection into a convection scheme. These improvements, when implemented into the National Center for Atmospheric Research (NCAR) Community Earth System Model (CESM) and the DOE Exascale Energy Earth System Model (E3SM), led to improved simulation of precipitation characteristics, including precipitation intensity, Intertropical Convergence Zone, and Madden-Julian oscillation. The project also improved the representation of microphysical processes in convective clouds. It enables the interaction of cloud microphysics with aerosols in weather and climate models.

54 ENVIRONMENTAL SCIENCES↗

Ventenata dubia projected to expand in the western United States despite future novel conditions

Abstract Distributions of both native and invasive species are expected to shift under future climate. Species distribution models (SDMs) are often used to explore future habitats, but sources of uncertainty including novel climate conditions may reduce the reliability of future projections. We explore the potential spread of the invasive annual grass ventenata ( Ventenata dubia ) in the western United States under both current and future climate scenarios using boosted regression tree models and 30 global climate models (GCMs). We quantify novel climate conditions, prediction variability arising from both the SDMs and GCMs, and the agreement among GCMs. Results demonstrate that currently suitable habitat is concentrated inside the invaded range of the northwest, but substantial habitat exists outside the invaded range in the Southern Rockies and southwestern US mountains. Future suitability projections vary greatly among GCMs, but GCMs commonly projected decreased suitability in the invaded range and increased suitability along higher elevations of interior mountainous areas. Climate novelty did not appear to undermine the prediction reliability in many cases where the climate–species relationship was fully represented by the occurrence data. GCM‐derived variability resulting from variation in future cool season precipitation and temperature seasonality was greatest in the Rocky Mountains. SDM‐derived variability was higher in currently suitable habitat, and few GCMs projections agreed that these areas would contain future suitable habitat. However, while prediction variability was high, many GCM projections agreed that parts of the Rocky, Wasatch, and Uinta Mountains would contain highly suitable habitat in the future. As disturbances in the interior mountains occur in coming decades, reducing some natural barriers to invasion, land managers, and conservationists will need to monitor for ventenata in post‐disturbance environments. Changes to invasion potential may not play out for several decades, but results related to current potential may have applications for early detection and rapid response planning.

Environmental Sciences & Ecology↗

Global Barotropic Tide Modeling Using Inline Self‐Attraction and Loading in MPAS‐Ocean

Abstract We examine ocean tides in the barotropic version of the Model for Prediction Across Scales (MPAS‐Ocean), the ocean component of the Department of Energy Earth system model. We focus on four factors that affect tidal accuracy: self‐attraction and loading (SAL), model resolution, details of the underlying bathymetry, and parameterized topographic wave drag. The SAL term accounts for the tidal loading of Earth's crust and the self‐gravitation of the ocean and the load‐deformed Earth. A common method for calculating SAL is to decompose mass anomalies into their spherical harmonic constituents. Here, we compare a scalar SAL approximation versus an inline SAL using a fast spherical harmonic transform package. Wave drag accounts for energy lost by breaking internal tides that are produced by barotropic tidal flow over topographic features. We compare a series of successively finer quasi‐uniform resolution meshes (62.9, 31.5, 15.7, and 7.87 km) to a variable resolution (45 to 5 km) configuration. We ran MPAS‐Ocean in a single‐layer barotropic mode forced by five tidal constituents. The 45 to 5 km variable resolution mesh obtained the best total root‐mean‐square error (5.4 cm) for the deep ocean ( 1,000 m) tide compared to TPXO8 and ran twice as fast as the quasi‐uniform 8 km mesh, which had an error of 5.8 cm. This error is comparable to those found in other forward (non‐assimilative) ocean tide models. In future work, we plan to use MPAS‐Ocean to study tidal interactions with other Earth system components, and the tidal response to climate change.

54 ENVIRONMENTAL SCIENCES↗

Deep decarbonization and U.S. biofuels production: a coordinated analysis with a detailed structural model and an integrated multisectoral model

Scenarios for deep decarbonization involve biomass for biofuels, biopower, and bioproducts, and they often include negative emissions via carbon capture and storage or utilization. However, critical questions remain about the feasibility of rapid growth to high levels of biomass utilization, given biomass and land availability as well as historical growth rates of the biofuel industry. We address these questions through a unique coordinated analysis and comparison of carbon pricing effects on biomass utilization growth in the United States using a multisectoral integrated assessment model, the Global Change Analysis Model (GCAM), and a biomass-to-biofuels system dynamics model, the Bioenergy Scenario Model (BSM). We harmonized and varied key factors—such as carbon prices, vehicle electrification, and arable land availability—in the two models. We varied the rate of biorefinery construction, the fungibility of feedstock types across conversion processes, and policy incentives in BSM. The rate of growth in biomass deployment under a carbon price in both models is within the range of current literature. However, the reallocation of land to biomass feedstocks would need to overcome bottlenecks to achieve growth consistent with deep decarbonization scenarios. Investments as a result of near-term policy incentives can develop technology and expand capacity—reducing costs, enabling flexibility in feedstock use, and improving stability—but if biomass demand is high, these investments might not overcome land reallocation bottlenecks. Biomass utilization for deep decarbonization relies on extraordinary growth in biomass availability and industrial capacity. In this paper, we quantify and describe the potential challenges of this rapid change.

09 BIOMASS FUELS↗

Integrating multimodal data through interpretable heterogeneous ensembles

Motivation: Integrating multimodal data represents an effective approach to predicting biomedical characteristics, such as protein functions and disease outcomes. However, existing data integration approaches do not sufficiently address the heterogeneous semantics of multimodal data. In particular, early and intermediate approaches that rely on a uniform integrated representation reinforce the consensus among the modalities but may lose exclusive local information. The alternative late integration approach that can address this challenge has not been systematically studied for biomedical problems. Results: We propose Ensemble Integration (EI) as a novel systematic implementation of the late integration approach. EI infers local predictive models from the individual data modalities using appropriate algorithms and uses heterogeneous ensemble algorithms to integrate these local models into a global predictive model. We also propose a novel interpretation method for EI models. We tested EI on the problems of predicting protein function from multimodal STRING data and mortality due to coronavirus disease 2019 (COVID-19) from multimodal data in electronic health records. We found that EI accomplished its goal of producing significantly more accurate predictions than each individual modality. It also performed better than several established early integration methods for each of these problems. The interpretation of a representative EI model for COVID-19 mortality prediction identified several disease-relevant features, such as laboratory test (blood urea nitrogen and calcium) and vital sign measurements (minimum oxygen saturation) and demographics (age). These results demonstrated the effectiveness of the EI framework for biomedical data integration and predictive modeling.

59 BASIC BIOLOGICAL SCIENCES↗

Challenging a Global Land Surface Model in a Local Socio-Environmental System

Land surface models (LSMs) predict how terrestrial fluxes of carbon, water, and energy change with abiotic drivers to inform the other components of Earth system models. Here, we focus on a single human-dominated watershed in southwestern Michigan, USA. We compare multiple processes in a commonly used LSM, the Community Land Model (CLM), to observational data at the single grid cell scale. For model inputs, we show correlations (Pearson’s R) ranging from 0.46 to 0.81 for annual temperature and precipitation, but a substantial mismatch between land cover distributions and their changes over time, with CLM correctly representing total agricultural area, but assuming large areas of natural grasslands where forests grow in reality. For CLM processes (outputs), seasonal changes in leaf area index (LAI; phenology) do not track satellite estimates well, and peak LAI in CLM is nearly double the satellite record (5.1 versus 2.8). Estimates of greenness and productivity, however, are more similar between CLM and observations. Summer soil moisture tracks in timing but not magnitude. Land surface reflectance (albedo) shows significant positive correlations in the winter, but not in the summer. Looking forward, key areas for model improvement include land cover distribution estimates, phenology algorithms, summertime radiative transfer modelling, and plant stress responses.

54 ENVIRONMENTAL SCIENCES↗

Investigating coastal backwater effects and flooding in the coastal zone using a global river transport model on an unstructured mesh

Abstract. Coastal backwater effects are caused by the downstream water level increase as a result of elevated sea level, high river discharge and their compounding influence. Such effects have crucial impacts on floods in densely populated regions but have not been well represented in large-scale river models used in Earth system models (ESMs), partly due to model mesh deficiency and oversimplifications of river hydrodynamics. Using two mid-Atlantic river basins as a testbed, we perform the first attempt to simulate the backwater effects comprehensively over a coastal region using the MOSART river transport model under an ESM framework, i.e., Energy Exascale Earth System Model (E3SM) configured on a regionally refined unstructured mesh, with a focus on understanding the backwater drivers and their long-term variations. By including sea level variations at the river downstream boundary, the model performance in capturing backwaters is greatly improved. We also propose a new flood event selection scheme to facilitate the decomposition of backwater drivers into different components. Our results show that while storm surge is a key driver, the influence of extreme discharge cannot be neglected, particularly when the river drains to a narrow river-like estuary. Compound flooding, while not necessarily increasing the flood peaks, exacerbates the flood risk by extending the duration of multiple coastal and fluvial processes. Furthermore, our simulations and analysis highlight the increasing strength of backwater effects due to sea level rise and more frequent storm surge during 1990–2019. Thus, backwaters need to be properly represented in ESMs to improve the predictive understanding of coastal flooding.

54 ENVIRONMENTAL SCIENCES↗