Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grid cells”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

TROPOMI observations allow for robust exploration of the relationship between solar-induced chlorophyll fluorescence and terrestrial gross primary production

Solar-induced chlorophyll fluorescence (SIF) observed by satellites has advanced the monitoring of terrestrial photosynthesis regionally and globally. The relationship between SIF and gross primary production (GPP) at leaf, canopy, and ecosystem scales has received tremendous attention in recent years. It remains controversial whether the SIF-GPP relationship at the ecosystem scale is universal or dependent upon vegetation type. New SIF observations from the TROPOspheric Monitoring Instrument (TROPOMI) with unprecedented high spatial and temporal resolution provide a new opportunity to elucidate the SIF-GPP relationship. Here, we examine the SIF-GPP relationship for seven major vegetation types across the U.S. with TROPOMI SIF and in-situ GPP data for 83 eddy covariance flux sites. We find that TROPOMI SIF shows a strong and consistent relationship with tower-based GPP at both satellite footprint and grid-cell levels. The slope of the SIF-GPP relationship is similar among all the vegetation types except croplands, demonstrating a nearly universal (converging to ~13.5 g C m -2 d -1 /W m -2 μm -1 sr -1) rather than vegetation type-specific SIF-GPP relationship. This confirms that TROPOMI SIF can be used as a proxy for GPP across a wide variety of vegetation types and can also be used to quantify GPP by avoiding uncertainty associated with land cover maps. The C 4 crops have a much higher slope than the C 3 crops, and therefore croplands tend to have a higher slope than C 3 -dominated vegetation types (e.g., forests, shrublands, savannas). We also find that the TROPOMI SIF is well correlated with GPP under normal or wetter conditions, while their relationship becomes weaker under water stress. In conclusion, our TROPOMI-based study could improve our understanding of the SIF-GPP relationship at the ecosystem scale and advance the mapping of GPP globally with SIF observations from space.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The impact of urban configuration types on urban heat islands, air pollution, CO 2 emissions, and mortality in Europe: a data science approach

The world is becoming increasingly urbanized. As cities around the world continue to grow, it is important for urban planners and policymakers to understand how different urban configuration patterns affect the environment and human health. We aimed at identifying European urban configuration types, based on the Local Climate Zones categories and street design variables from Open Street Map, and evaluating their association with motorized traffic flows, Surface Urban Heat Island (SUHI) intensities, tropospheric nitrogen dioxide (NO 2 ), CO 2 per capita emissions and age-standardized mortality. We considered 946 European cities from 31 countries for the analysis defined in the 2018 Urban Audit database, of which 919 European cities were analysed. Data were collected at a 250 m × 250 m grid cell resolution. We divided all cities into five concentric rings based on the Burgess concentric urban planning model and calculated the mean values of all variables for each ring. First, to identify distinct urban configuration types, we applied the Uniform Manifold Approximation and Projection for Dimension Reduction method, followed by the k-means clustering algorithm. Next, statistical differences in exposures (including SUHI) and mortality between the resulting urban configuration types were evaluated using a Kruskal–Wallis test followed by a post-hoc Dunn's test. We identified four distinct urban configuration types characterising European cities: compact high density (n=246), open low-rise medium density (n=245), open low-rise low density (n=261), and green low density (n=167). Compact high density cities were a small size, had high population densities, and a low availability of natural areas. In contrast, green low-density cities were a large size, had low population densities, and a high availability of natural areas and cycleways. The open low-rise medium and low-density cities were a small to medium size with medium to low population densities and low to moderate availability of green areas. Motorised traffic flows and NO 2 exposure were significantly higher in compact high density and open low rise medium density cities when compared with green low density and open low-rise low density cities. Additionally, green low-density cities had a significantly lower SUHI effect compared with all other urban configuration types. Per person CO 2 emissions were significantly lower in compact high density cities compared with green low density cities. Lastly, green low density cities had significantly lower mortality rates when compared with all other urban configuration types. Our findings indicate that, although the compact city model is more sustainable, European compact cities still face challenges related to poor environmental quality and health. Our results have notable implications for urban and transport planning policies in Europe and contribute to the ongoing discussion on which city models can bring the greatest benefits for the environment, climate, and health.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2019

Evaluation of Implicit‐Explicit Additive Runge‐Kutta Integrators for the HOMME‐NH Dynamical Core

Abstract The nonhydrostatic High‐Order Method Modeling Environment (HOMME‐NH) atmospheric dynamical core supports acoustic waves that propagate significantly faster than the advective wind speed, thus greatly limiting the time step size that can be used with standard explicit time integration methods. Resolving acoustic waves is unnecessary for accurate climate and weather prediction. This numerical stiffness is addressed herein by considering implicit‐explicit additive Runge‐Kutta (ARK IMEX) methods that can treat the acoustic waves in a stable manner without requiring implicit treatment of nonstiff modes. Various ARK IMEX methods are evaluated for their efficiency in producing accurate solutions, ability to take large time step sizes, and sensitivity to grid cell length ratio. Both the gravity wave test and baroclinic instability test from the 2012 Dynamical Core Model Intercomparison Project are used to recommend 5 of the 27 ARK IMEX methods tested for use in HOMME‐NH.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A Global LIS/OTD Climatology of Lightning Flash Extent Density

Previous lightning climatologies derived from Lightning Imaging Sensor (LIS) and Optical Transient Detector (OTD) total lightning measurements have quantified lightning frequency as a Flash Rate Density (FRD). This approach assumes that lightning flashes can be represented as points and quantifies the frequency of lightning centered in each grid cell. However, lightning has a finite extent that can reach hundreds of kilometers. A new climatology based on Flash Extent Density (FED) is constructed for LIS (including ISS-LIS) and OTD that accounts for the horizontal dimension of lightning. The FED climatology documents the frequency that an observer can expect lightning to be visible overhead—regardless of where the flash began or ended. This new FED climatology confirms and elaborates on the previous global LIS/OTD FRD and Americas-only Geostationary Lightning Mapper (GLM) findings. The FED climatology maintains Lake Maracaibo as the global lightning hotspot with an average of 389 flashes/day but designates Karabre in the Democratic Republic of the Congo as the global thunderstorm duty (percent of the total viewtime where lightning is observed) hotspot at 7.29%. Meanwhile, Kuala Lumpur is the national capital city with the most lightning, and its airport (KUL) is the top major airport affected by lightning (in terms of duty). The FED seasonal cycle and month-to-month changes in the “center of lightning” for the three continental chimney regions are also discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Topographical Controls on Hillslope-Scale Hydrology Drive Shrub Distributions on the Seward Peninsula, Alaska

Observations indicate shrubs are expanding across the Arctic tundra, mainly on hillslopes and primarily in response to climate warming. However, the impact topography exerts on hydrology, nutrient dynamics, and plant growth can make untangling the mechanisms behind shrub expansion difficult. We examined the role topography plays in determining shrub expansion by applying a coupled transect version of a mechanistic ecosystem model (ecosys) in a tundra hillslope site in the Seward Peninsula, Alaska. Modeled biomass of the dominant plant functional types agreed well with field measurements (R 2 = 0.89) and accurately represented shrub expansion over the past 30 years inferred from satellite observations. In the well-drained crest position, canopy water potential and plant nitrogen (N) uptake was modeled to be low from plant and microbial water stress. Intermediate soil water content in the mid-slope position enhanced mineralization and plant N uptake, increasing shrub biomass. The deciduous shrub growth in the mid-slope position was further enhanced by symbiotic N 2 fixation primed by increased root carbon allocation. The gentle slope in the poorly drained lower-slope position resulted in saturated soil conditions that reduced soil O 2 concentrations, leading to lower root O 2 uptake and lower nutrient uptake and plant biomass. Additionally, a simulation that removed topographical interconnectivity between grid cells resulted in (1) a 28% underestimate of mean shrub biomass and (2) over or underestimated shrub productivity at the various hillslope positions. Our results indicate that land models need to account for hillslope-scale coupled surface and subsurface hydrology to accurately predict current plant distributions and future trajectories in Arctic ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

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↗
DOE OSTI2025

Atmospheric River Frequency-Category Characteristics Shape U.S. West Coast Runoff

Abstractrunoff response to atmospheric rivers (ARs) over the U.S. West Coast. We focused on runoff time series variations impacted by AR characteristics (e.g., category and frequency) and land preconditions during Northern Hemisphere cool seasons in the period of 1940–2023. Results show that high-category ARs significantly increase local runoff with higher hourly precipitation rates leading to a greater incremental rate and peak runoff. Extreme runoff increases greatly with the AR category with an increase rate up to 12.5 times stronger than non-extreme runoff. Besides the AR category, land preconditions such as soil moisture and snowpack also play crucial roles in modulating runoff response. We found that runoff induced by weak-category ARs is more sensitive to land preconditions than high-category ARs, with high peak runoff occurring when soil is nearly saturated. Additionally, more than 50% of high-peak-runoff events in snow-covered grid cells are associated with rain-on-snow events particularly for the events associated with weaker ARs. Regression analysis reveals that AR precipitation and land preconditions jointly influence runoff, emphasizing the importance of including soil moisture and snowpack levels in AR impact assessments. The study also highlights the intensified runoff response to back-to-back ARs with short intervals, which may become more frequent with climate warming, posing increased flood risks via facilitating wet soil conditions. Our findings have significant implications for AR risk predictions and the development of prediction models for AR-induced runoff.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Uncertainties in Simulating Flooding During Hurricane Harvey Using 2D Shallow Water Equations

Abstract Flooding is one of the most impactful weather‐related natural hazards. Numerical models that solve the two dimensional (2D) shallow water equations (SWE) represent the first‐principles approach to simulate all types of spatial flooding, such as pluvial, fluvial, and coastal flooding, and their compound dynamics. High spatial resolution (e.g., () m) is needed in 2D SWE simulations to capture flood dynamics accurately, resulting in formidable computational challenges. Thus, relatively coarser spatial resolutions are used for large‐scale simulations of flooding, which introduce uncertainties in the results. It is unclear how the uncertainty associated with the model resolution compares to the uncertainties in precipitation data sets and assumptions regarding boundary conditions when channelized flows interact with other water bodies. In this study, we compare these three sources of uncertainties in 2D SWE simulations for the 2017 Houston flooding event. Our results show that precipitation uncertainty and mesh resolution have more significant impacts on the simulated streamflow and inundation dynamics than the choice of the downstream boundary condition at the watershed outlet. We point out the viability to confine the uncertainty of coarsening mesh resolution by using the variable resolution mesh (VRM) which refines critical topographic features with far fewer grid cells. Specifically, in simulations with VRM, the simulated inundation depths over the refined region are comparable to that use the finest uniform mesh. This study contributes to understanding the challenges and pathways for applying 2D SWE models to improve the realism of flood simulations over large scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Spatially explicit terrestrial carbon densities for calibrating the carbon cycle in human-Earth system Models

Soil and vegetation carbon stocks play a critical role in human-Earth system models. These stocks (denominated as densities in MgC/ha) affect variables such as land use change emissions and also influence land use change pathways under climate forcing scenarios where terrestrial carbon is assigned a carbon price. Here we present reharmonized soil and vegetation carbon densities both at the 5-arcmin resolution grid cell level and also aggregated to 235 water sheds for 4 land use types (Cropland, Grazed land, Urban land and unmanaged vegetation) and 15 unmanaged land cover types. Moreover, we use the distribution of carbon within and across pixels to define statistical "states" of carbon, once again differentiated by land type. These statistical states are used to define a range of possible carbon values that can be used for defining initial conditions of soil and vegetation carbon in human-Earth system models. We implement these data in a state-of-the-art multi sector dynamics model, namely the Global Change Analysis Model (GCAM), and show that these new data improve several land use responses, especially when terrestrial carbon is assigned a carbon price.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

20% of US electricity from wind will have limited impacts on system efficiency and regional climate

Abstract Impacts from current and future wind turbine (WT) deployments necessary to achieve 20% electricity from wind are analyzed using high resolution numerical simulations over the eastern USA. Theoretical scenarios for future deployments are based on repowering (i.e. replacing with higher capacity WTs) thus avoiding competition for land. Simulations for the contemporary climate and current WT deployments exhibit good agreement with observed electricity generation efficiency (gross capacity factors (CF) from simulations = 45–48%, while net CF for WT installed in 2016 = 42.5%). Under the scenario of quadrupled installed capacity there is a small decrease in system-wide efficiency as indicated by annual mean CF. This difference is approximately equal to that from the two simulation years and may reflect saturation of the wind resource in some areas. WT modify the local near-surface climate in the grid cells where they are deployed. The simulated impact on near-surface climate properties at both the regional and local scales does not increase with increasing WT installed capacity. Climate impacts from WT are modest compared to regional changes induced by historical changes in land cover and to the global temperature perturbation induced by use of coal to generate an equivalent amount of electricity.

17 WIND ENERGY↗
DOE OSTI2024

Evaluation of precipitation across the contiguous United States, Alaska, and Puerto Rico in multi-decadal convection-permitting simulations

Abstract This study is an early effort to generate a multi-decadal convection-permitting regional climate dataset that covers nearly the entire North American continent. We assessed a 20 year dynamically downscaled regional climate simulation at a 4 km spatial resolution with explicit convection across the contiguous United States (CONUS), Alaska, and Puerto Rico. Specifically, we evaluated the model’s performance in representing mean, 95th percentile, and extreme precipitation across regions. Our findings indicate that when compared with ERA5 reanalysis, the forcing data, convection-permitting simulation improves representations of seasonal, 95th percentile, and extreme precipitation over a large portion of the CONUS, Alaska, and Puerto Rico, particularly in areas where precipitation is heaviest. The simulation adds value over its forcing data (ERA5) in up to 53% of all grid cells in the CONUS, 68.8% in Alaska, and 84.0% in Puerto Rico. It is important to note that, however, despite improvements, model errors in Puerto Rico remain large. Similar improvements are observed in extreme indices, including consecutive dry days, maximum 5 days precipitation, and extreme precipitation. Analysis of the diurnal cycle of mean hourly precipitation suggests that representations of convective processes—including onset, dissipation, suppression, downstream propagation, and local circulation—improved overall.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within ≲3σ of the mean of all simulations on a given grid (128 3 and 256 3 grid cells), where σ describes the statistical variation due to the flow’s time dependence. They also agree well with a 512 3 reference run. The 128 3 and 256 3 simulations agree within 9% and 4%, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Numerical simulation of a relativistic magnetron using a fluid electron model

An approach to numerically modeling relativistic magnetrons, in which the electrons are represented with a relativistic fluid, is described. A principal effect in the operation of a magnetron is space-charge-limited (SCL) emission of electrons from the cathode. We have developed an approximate SCL emission boundary condition for the fluid electron model. This boundary condition prescribes the flux of electrons as a function of the normal component of the electric field on the boundary. Here we show the results of a benchmarking activity that applies the fluid SCL boundary condition to the one-dimensional Child–Langmuir diode problem and a canonical two-dimensional diode problem. Simulation results for a two-dimensional A6 magnetron are then presented. Computed bunching of the electron cloud occurs and coincides with significant microwave power generation. Numerical convergence of the solution is considered. Sharp gradients in the solution quantities at the diocotron resonance, spanning an interval of three to four grid cells in the most well-resolved case, are present and likely affect convergence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Unsupervised discovery of extreme weather events using universal representations of emergent organization

Spontaneous self-organization is ubiquitous in systems far from thermodynamic equilibrium. While organized structures that emerge dominate transport properties, universal representations that identify and describe these key objects remain elusive. Here, we introduce a theoretically grounded framework for describing emergent organization that, via data-driven algorithms, is constructive in practice. Its building blocks are spacetime lightcones that embody how information propagates across a system through local interactions. We show that predictive equivalence classes of lightcones—local causal states—capture organized behaviors in complex spatiotemporal systems. Employing an unsupervised physics-informed machine learning algorithm and a high-performance computing implementation, we demonstrate automatically discovering organized structures in two real-world domain science problems. We show that local causal states identify vortices and track their power-law decay behavior in two-dimensional fluid turbulence. We then show how to detect and track familiar extreme weather events—hurricanes and atmospheric rivers—and discover other novel structures associated with precipitation extremes in high-resolution climate data at the grid-cell level.

Rupe, Adam [Pacific Northwest National Laboratory ↗
DOE OSTI2024

Changing windstorm characteristics over the US Northeast in a single model large ensemble

Abstract Extreme windstorms pose a significant hazard to infrastructure and public safety, particularly in the highly populated US Northeast (NE). However, the influence climate change and changing land use will have on these events remains unclear. A large ensemble generated using the Max-Planck Institute (MPI) Earth system model is used to generate projections of NE windstorms under different shared socioeconomic pathways (SSPs) and to attribute changes to projected land use land cover (LULC) change, externally forced changes and internal climate variability. To reduce the influence of coarse grid cell resolution and uncertainties in surface roughness lengths, windstorms are identified using simultaneous widespread exceedance of local 99th percentile 10 m wind speeds (U 99 ). Projected declines in forest cover in the NE and the resulting reductions in surface roughness length under SSP3-7.0 lead to projections of large increases in U 99 and derived windstorm intensity and scale. However, these projected changes in regional LULC under SSP3-7.0 are unprecedented in a historical context and may not be realistic. After corrections are applied to remove the influence of LULC on wind speeds, regionally averaged U 99 exhibit declines for most of the single model initial-condition large ensemble (SMILE) members which are broadly proportional to the radiative forcing and global air temperature increase in the SSPs, with a median value of −0.15 ms −1 °C −1 . While weak cyclones are projected to decline in frequency in the NE, intense cyclones and the resulting windstorms and indices of socioeconomic loss do not. Where present, significant trends in these loss indices are positive, and some MPI SMILE members generate future windstorms that are unprecedented in the historical period.

Coburn, Jacob (ORCID:0000000309538117)↗
DOE OSTI2026

Diagnosing the representation of surface and layered soil moisture in Earth system models

Surface soil moisture (mrsos) and vertically integrated soil moisture (mrsol) over the top 10 cm should, by definition, be physically consistent in Earth System Models (ESMs). However, an evaluation of nine CMIP6 models reveals substantial inconsistencies: in some models, mrsos and integrated mrsol agree globally; in others, they align only in specific regions; and in a few, they diverge across all grid cells. These discrepancies arise from a combination of factors, including metadata errors, inconsistent variable definitions, or diagnostic sequencing within the model. We demonstrate how such issues can lead to significant biases, even when both variables are present and seemingly well-defined. As model complexity increases and multi-model comparisons become more common, assumptions about variable equivalence may lead to flawed conclusions. This study highlights the need for routine consistency checks, improved metadata standards, and community-wide practices that ensure reliability of derived variables across ESM outputs, particularly in preparation for CMIP7.

Earth system models↗
DOE OSTI2025

Spatiotemporal predictions of toxic urban plumes using deep learning

Industrial accidents, chemical spills, and structural fires can release large amounts of harmful materials that disperse into urban atmospheres and impact populated areas. Computer models are typically used to predict the transport of toxic plumes by solving fluid dynamical equations. However, these models can be computationally expensive due to the need for many grid cells to simulate turbulent flow and resolve individual buildings and streets. In emergency response situations, alternative methods are needed that can run quickly and adequately capture important spatiotemporal features. Here, we present a novel deep learning model called ST-GasNet inspired by the mathematical equations that govern the behavior of plumes as they disperse through the atmosphere. ST-GasNet learns the spatiotemporal dependencies from a limited set of temporal sequences of ground-level toxic urban plumes generated by a high-resolution large eddy simulation model. On independent sequences, ST-GasNet accurately predicts the late-time spatiotemporal evolution, given the early-time behavior as an input, even when a building splits a large plume into smaller plumes. By incorporating large-scale wind boundary condition information, ST-GasNet achieves a prediction accuracy of at least 90% on test data for the entire prediction period.

Civil and Environmental Engineering↗
DOE OSTI2024

At Risk Population Estimates for Belarus, Poland and Slovakia with Machine Learning

High-resolution gridded population modeling is crucial for various applications, including disaster response planning, infectious disease spread modeling, climate change impact estimation, policy development, and more. Multiple gridded population datasets have been developed, each tailored to meet specific objectives. Among them, LandScan Global dataset is designed to represent ambient and unwarned population distributions. However, this dataset relies on a statistical approach that requires manual adjustments, making it time consuming and labour intensive. Existing machine learning (ML) methods often train and test at different spatial resolutions, potentially leading to inflated results, and they rely on Census population totals for disaggregation. To address these limitations, in this study we developed population estimates using ML models trained and tested at a consistent 30 arc-second resolution (≈1 square kilometer), specifically using Random Forest (RF) and XGBoost. These models were trained on 2020 datum to predict for 2021 for three countries: Belarus, Poland, and Slovakia. Our findings show that both RF (MAE varies from 5.75 to 13.25) and XGBoost (MAE varies from 8.15 to 23.44) model performance is close to LandScan Global estimates. Furthermore, neither of the models performed the best across all grid cells: the RF model was more effective in areas with lower populations, while XGBoost excelled in more densely populated regions. The proposed approach can be used for countries where the Census data is not available.

Lebakula, Viswadeep [ORNL] (ORCID:0000000152935914↗