Search NASA⌕ Search

SEARCH · Search NASA

Results for “spatially explicit modeling”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Evaluating ecosystem water use efficiency under drought stress: a case study of the Helan Mountain region, northwest China

Context Water use efficiency (WUE) is a fundamental ecological indicator links carbon assimilation and water loss in terrestrial ecosystems. Understanding its responses to drought stress is essential for adaptive ecosystem management, particularly in climate-sensitive mountain landscapes. Objectives This study aimed to investigate drought-driven variations in WUE across major vegetation types in the Helan Mountain region of Northwest China. Specifically, we sought to identify dominant ecological drivers of WUE variability and to disentangle their relative importance and causal pathways. Methods We quantified WUE using the Moderate Resolution Imaging Spectroradiometer (MODIS) products and the Drought Severity Index (DSI) data from 2001 to 2020. To examine WUE – drought relationships across contrasting vegetation types, we employed a spatially explicit analytical framework integrating Random Forest (RF) modeling, partial correlation analysis, and structural equation modeling (SEM). Results Regional WUE exhibited relatively stable interannual dynamics, yet pronounced spatial heterogeneity that was strongly modulated by drought conditions. Vegetation properties, particularly Leaf Area Index (LAI) and Normalized Difference Vegetation Index (NDVI), emerged as the dominant determinants of WUE, with NDVI alone explaining over 20% of its spatial variance in forest and grassland during non-drought periods. SEM analyses revealed that climate forcing influenced WUE mainly through indirect pathways mediated by soil moisture availability and vegetation structural dynamics, rather than through direct climatic controls. Among all regulating factors, LAI acted as the central control node governing ecosystem carbon–water coupling. In contrast, short-term climatic stress, especially atmospheric demand and drought duration, exerted weak or negative direct effects on WUE. Ecosystem-specific responses were observed, with croplands mainly regulated by soil water availability, whereas forests and grasslands showed more sensitive to atmospheric drought stress. Together, these results reveal a hierarchical control framework where soil–vegetation interactions mediate climate impacts on WUE, driving strong spatial heterogeneity in drought responses across mountain landscapes. Conclusions Our findings highlight the pivotal role of indirect drought effects mediated by vegetation and soil processes in shaping ecosystem WUE. The identified soil–vegetation–climate regulatory hierarchy provides mechanistic insight into landscape–scale drought sensitivity and supports integrated modeling approaches for evaluating ecosystem resilience and sustainable management in arid mountain regions.

China↗
DOE OSTI2021

Tropical Cirrus in Global Storm-Resolving Models. Part II: Cirrus Life Cycle and Top-of-Atmosphere Radiative Fluxes

Cirrus clouds of various thicknesses and radiative characteristics extend over much of the tropics, especially around deep convection. They are difficult to observe due to their high altitude and sometimes small optical depths. They are also difficult to simulate in conventional global climate models, which have coarse grid spacings and simplified parameterizations of deep convection and cirrus formation. We investigate the representation of tropical cirrus in global storm-resolving models (GSRMs), which have higher spatial resolution and explicit convection and could more accurately represent cirrus cloud processes. This study uses GSRMs from the DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains (DYAMOND) project. The aggregate life cycle of tropical cirrus is analyzed using joint albedo and outgoing longwave radiation (OLR) histograms to assess the fidelity of models in capturing the observed cirrus cloud populations over representative tropical ocean and land regions. The proportions of optically-thick deep convection, anvils, and cirrus vary across models and are portrayed in the vertical distribution of cloud cover and top-of-atmosphere radiative fluxes. Model differences in cirrus populations, likely driven by subgrid processes such as ice microphysics, dominate over regional differences between convectively-active tropical land and ocean locations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Application of a Distributed Element Roughness Model to Additively Manufactured Internal Cooling Channels

Abstract Design for cooling effectiveness in turbine blades relies on accurate models for dynamic losses and heat transfer of internal cooling passages. Metal additive manufacturing (AM) has expanded the design space for these configurations, but can give rise to large-scale roughness features. The range of roughness length scales in these systems makes morphology resolved computational fluid dynamics (CFD) impractical. However, volumetric roughness models can be leveraged, as they have computational costs orders of magnitude lower. In this work, a discrete element roughness model (DERM), based on the double-averaged Navier–Stokes equations, is presented and applied to additively manufactured rough channels, representative of gas turbine blade cooling passages. Unique to this formulation of DERM is a generalized sheltering-based treatment of drag, a two-layer model for spatially averaged Reynolds stresses, and explicit treatment of dispersion. Six different AM rough surface channel configurations are studied, with roughness trough to peak sizes ranging from 15% to 60% nominal channel passage half-width, and the roughness Reynolds number ranges from Rek = 60 to 300. DERM predictions for spatially and temporally averaged mean flow quantities are compared to previously reported direct numerical simulation results. Good agreement in the mean velocity profiles, stress balances, and drag partitions are observed. While DERM models are typically calibrated to specific deterministic roughness morphologies at comparatively small roughness Reynolds numbers, the present more generalized DERM formulation has wider applicability. Here, it is demonstrated that the model can accommodate random roughness of large scale, typical of AM.

Engineering↗
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 OSTI2024

Data for "Divergent urban land trajectories under alternative population projections within the shared socioeconomic pathways"

Datasets supporting findings and visualizations behind McManamay et al. (2024) Divergent urban land trajectories under alternative population projections within the shared socioeconomic pathways. Environmental Research Letters, DOI: 10.1088/1748-9326/ad2eec Using a spatial modeling experiment at high resolution (1-km), this study compared how two alternative US population projections, varying in the spatially explicit nature of demographic patterns and migration, affect urban land dynamics simulated by the Spatially Explicit, Long-term, Empirical City development (SELECT) model for SSP2, SSP3, and SSP5. The numerical experiment, inputs and outputs are fully described in McManamay et al. (2024). The datasets summarize SELECT model simulations within urban areas and rural areas of the conterminous United States. For code to reproduce the results, please refer to https://github.com/IMMM-SFA/mcmanamay_etal_2024_erl Please refer to the README file provided in Files for more details. Descriptions of the datasets are provided below. Dataset(s) Descriptions: Urban_Land_delta_UA_County.csv: Urban land area for original (default) population, updated population, and urban land delta (difference between the two) according to Shared Socioeconomic Pathways (SSPs) and year. Urban land is summarized for U.S. counties (FIPS code) and urban areas (GEOID) as unique areas within various counties. JO_Gao_pop_ufdelta.csv: Comparison of population differences and urban land differences between the origial (default) and updated population projections. City_Case_studies.csv: Differences in urban land areas arising from different population projections summarized for selected cities (Atlanta, Los Angeles, Houston, and New York) and their surrounding rural areas UA_county_FID.zip: Zipped folder of .tiff file depicting unique combination of urban areas and US counties used for summarizing urban land differences in urban and rural areas based on different population projections. ULD_clusters.zip: Zipped folder of .csv files depicting normalized changes in urban land delta (difference in urban land arising from population projections) and the respective Ward's hierarchical clusters, which group urban and rural areas based on similarities in temporal trends. Urban_Percent_change_2100.zip: Zipped folder of .csv files depicting percent changes in urban land area for year 2100 based on differences in the original (default) population and updated population. Each file is for SSPs and urban or rural areas.

Land↗
DOE OSTI2020

Spatially and Temporally Explicit Life Cycle Environmental Impacts of Soybean Production in the U.S. Midwest

Quantifying the spatially and temporally explicit life cycle environmental impacts of crop production is critical for designing sustainable supply chains for biofuel and animal sectors. This study assessed life cycle environmental impacts of soybean production in around 1000 Midwest counties over 9 years. A combination of Environmental Policy Integrated Climate model and process-based life cycle assessment model was used to estimate the spatially and temporally explicit life cycle global warming (GW), eutrophication (EU) and acidification (AD) impacts. Sequentially, a machine learning approach was applied to identify the top influential factors among soil, climate and farming practices for spatially and temporally explicit life cycle environmental impacts. The results indicated that significant variations existed in life cycle environmental impacts among counties and across years. Life cycle GW impacts ranged from -11.4 to 22.0 kg CO2-eq. kg soybean-1, whereas life cycle EU and AD impacts varied by a factor of 302 and 44, respectively. Nitrogen application rates, temperature in March and soil texture were the top influencing factors for the spatial-temporal variations in life cycle GW impacts. In contrast, soil organic content and nitrogen application rate were the top influencing factors for the spatial-temporal variations in life cycle EU and AD impacts.

Soybean, life cycle analysis, sustainability, spat↗
DOE OSTI2024

Fertilizer management for global ammonia emission reduction

Crop production is a large source of atmospheric ammonia (NH 3 ), which poses risks to air quality, human health and ecosystems. However, estimating global NH 3 emissions from croplands is subject to uncertainties because of data limitations, thereby limiting the accurate identification of mitigation options and efficacy. In this report we develop a machine learning model for generating crop-specific and spatially explicit NH 3 emission factors globally (5-arcmin resolution) based on a compiled dataset of field observations. We show that global NH 3 emissions from rice, wheat and maize fields in 2018 were 4.3 ± 1.0 Tg N yr -1 , lower than previous estimates that did not fully consider fertilizer management practices. Furthermore, spatially optimizing fertilizer management, as guided by the machine learning model, has the potential to reduce the NH3 emissions by about 38% (1.6 ± 0.4 Tg N yr -1 ) without altering total fertilizer nitrogen inputs. Specifically, we estimate potential NH3 emissions reductions of 47% (44–56%) for rice, 27% (24–28%) for maize and 26% (20–28%) for wheat cultivation, respectively. Under future climate change scenarios, we estimate that NH3 emissions could increase by 4.0 ± 2.7% under SSP1–2.6 and 5.5 ± 5.7% under SSP5–8.5 by 2030–2060. However, targeted fertilizer management has the potential to mitigate these increases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Spatial Microbial Respiration Variations in the Hyporheic Zones Within the Columbia River Basin

Abstract While the hyporheic zone (HZ) accounts for a significant portion of whole stream CO 2 concentrations, HZ respiration modeling studies are lacking in quantifying their contributions to the total CO 2 at large watershed/basin scales. Quantifying the contribution of anaerobic respiration is also underappreciated. This study used a carbon‐nitrogen‐coupled river corridor model to quantify HZ aerobic and anaerobic respiration and determined the key factors controlling their spatial variability within the Columbia River Basin (CRB). The modeled respiration patterns showed high spatial variability. Among the nine sub‐basins composing the CRB, the Lower Columbia and the Willamette, which receive higher precipitation, had higher respiration. Medium‐sized rivers (fourth to sixth orders) produced the highest aerobic and anaerobic respiration among reaches of different sizes. At the basin scale, aerobic respiration is dominant, representing approximately 98.7% of the total respiration across the CRB. While most of the reaches were dominant with aerobic respiration, reaches in agricultural land showed a relatively higher anaerobic respiration (18%) ratio. A variable importance analysis showed that hyporheic exchange flux controlled most of the spatial variability of HZ respiration, dominating over other physical variables such as residence time, stream dissolved organic carbon (DOC), nitrate, and dissolved oxygen (DO). The influence of substrate concentration (DOC and DO) is larger in modeling anaerobic respiration than aerobic respiration. Future efforts will focus on improving the estimation of the HZ exchange flux and the implementation of spatially explicit parameterizations for the reactions of interest to reduce model uncertainty.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

IM3 Phase 2 Official Simulations: GCAM-Demeter-SELECT Annualized Land Use and Land Cover, Wood Harvest and Fertilization Data with Dynamic Urbanization Harmonized to CLM Land Definitions at 0.125 Degrees

Annualized land use land cover data, including wood harvest and fertilizer use data from the Global Change Analysis Model (GCAM) downscaled to 0.125 degrees for couping with the Community Land Model (CLM). GCAM here refers to GCAM-USA v5.3.im3 which has an enhanced electricity sector and an updated data system needed to represent regional to local scale dynamics. Data is also harmonized with future urbanization projections from the Spatially-Explicit, Long-term, Empirical City developmenT (SELECT) model. Projections/Data are generated using the demeter land use and land cover downscaling model. Original projections were generated at 0.05 degrees before being aggregared to 0.125 degrees. Projections are available for 8 alternative scenarios. Two versions of final data are included- one with managed forests or harvested forest area per pixel broken out and one with the same aggregated into total forests. Following folders are included: demeter_78_PFT_output:This is the final output of dynamic land use land cover change for 78 PFTs as required by CLM raw_outputs_incl_managed_forest: This is the final output but with managed forests broken out as a different PFT. Essentially a 79th PFT is added. wood_harvest_outputs: Wood harvest output per pixel in gC/m2 fertilization_outputs: Fertilizer use per pixel in gN/m2 Each NetCDF file in each folder represents a projection for a separate year, scenario. Land use outputs are organized as PFT level data saved as subdata. Link to GCAM version used- https://data.msdlive.org/records/yb23g-44274 Link to SELECT documentation -https://www.sciencedirect.com/science/article/pii/S1364815219301707 Link to CLM documentation- - https://www.cesm.ucar.edu/models/clm In case of questions contact- kanishka.narayan@pnnl.gov

GCAM-USA↗
DOE OSTI2026

Spatial Optimization of Multiscale Biorefinery Deployment for a Diversified Bioeconomy in the United States

Strategic biorefinery siting is critical for a diversified bioeconomy, yet industry, policy, and research often focus on either large-scale biofuel plants or smaller-scale specialty bioproduct facilities, with limited coordination across scales. We address this gap by modeling biorefinery deployment spanning a 28-fold difference in capacity. We developed an open-source, spatially explicit framework integrating techno-economic analysis with logistics and refinery cost surrogate models to evaluate multiscale miscanthus-derived biorefineries across the rainfed U.S. for the production of ethanol, succinic acid, lactic acid, potassium sorbate, and acrylic acid. Overall costs change little as feedstock density increases, while transport distances decrease by ∼30 to 67% (∼100 km) and siting flexibility improves. Specifically, a 5-fold feedstock density increase (2% to 10% of suitable land) reduces minimum selling prices by <10% (e.g., 0.27 USD·gal –1 for ethanol). This limited economic sensitivity suggests dense planting is not required for competitive deployment, particularly for smaller-scale facilities. Representing collection areas as irregular rather than circular expands the feasible space under low-density scenarios. While large-scale refineries anchor regional supply chains, smaller facilities retain spatial flexibility even when large refineries are established. These findings highlight the importance of spatial representation and multiscale coordination for robust, regionally tailored biomanufacturing networks to advance renewable carbon integration without extensive land conversion.

biorefinery siting↗
DOE OSTI2020

Willow Biomass Crops Are a Carbon Negative or Low-Carbon Feedstock Depending on Prior Land Use and Transportation Distances to End Users

Few life cycle assessments (LCAs) on willow biomass production have investigated the effects of key geographically specific parameters. This study uses a spatial LCA model for willow biomass production to determine spatially explicit greenhouse gas (GHG) emissions and energy return on investment (EROI), including land use conversion from pasture and cropland or grassland. There were negative GHG emissions on 92% of the land identified as suitable for willow biomass production, indicating this system’s potential for climate change mitigation. For willow planted on cropland or pasture, life cycle GHG emissions ranged from −53.2 to −176.9 kg CO2eq Mg-1. When willow was grown on grassland the projected decrease in soil organic carbon resulted in a slightly positive GHG balance. Changes in soil organic carbon (SOC) associated with land use change, transportation distance, and willow yield had the greatest impacts on GHG emissions. Results from the uncertainty analysis exhibited large variations in GHG emissions between counties arising from differences in these parameters. The average EROI across the entire region was 19.2. Willow biomass can be a carbon negative or low-carbon energy source with a high EROI in regions with similar infrastructure, transportation distances, and growing conditions such as soil characteristics, land cover types, and climate.

09 BIOMASS FUELS↗
DOE OSTI2026

Forest aboveground biomass estimation through integration of sentinel-2 and PALSAR-2 time series: assessing models trained on GEDI and field inventory benchmarks

Accurate and spatially explicit forest Aboveground Biomass (AGB) mapping through remote sensing is critical for quantifying terrestrial carbon stocks and informing effective forest management strategies. However, AGB estimation in dense forests with complex terrain remains challenging due to satellite sensor signal saturation problem (saturation issue occurs in high biomass forests), structural complexity, and limited ground truth for calibration. This study presents a novel framework that integrates multi-temporal Sentinel-2 optical imagery, ALOS PALSAR-2 Synthetic Aperture Radar (SAR) data, and topographic variables with explainable Machine Learning to map AGB across mountainous forests within subtropical and temperate oceanic climate zones of Mexico. We evaluate the effects of temporal granularity and sensor synergy by comparing multiple temporal inputs and sensor configurations (Sentinel-2, PALSAR-2, and their fusion), and assess model performance using two reference datasets: NASA GEDI LiDAR-derived biomass and Mexico’s National Forest and Soil Inventory (INFyS). Our results showed that models trained on INFyS consistently outperformed those trained on GEDI, highlighting limitations in GEDI’s reliability in biomass estimates within this study region. Furthermore, the integration of Sentinel-2 and PALSAR-2 provided improved predictions compared to single-sensor models, particularly when combined with temporally explicit yearly statistics. The best-performing model, which was trained on INFyS data, and considered both Sentinel-2 and PALSAR-2 yearly statistics, as well as topographic variables, achieved an R2 of 0.64, RMSE of 51.10 Mg/ha, and relative RMSE (rRMSE) of 58.69%. Explainable ML analysis identified Sentinel-2 spectral indices and topographic features as key predictors, while PALSAR-2 metrics provided complementary information, partially mitigating saturation effects in high-biomass areas. Specifically, integrating both sensors substantially improved AGB estimation in high biomass forest (≥200 Mg/ha), yielding 98% gains over optical-only model, with resulting estimates exceeding GEDI L4B by 29% and ESA-CCI-BIOMASS by 174%. Terrain-stratified analysis indicated close agreement with GEDI in low-slope areas, with increasing divergence as slope steepness increased, while estimates remained consistently higher than ESA-CCI-BIOMASS across all slope classes. The proposed approach advances multi-sensor fusion and temporal feature engineering for AGB mapping using open-access satellite datasets, providing a scalable and reproducible framework for annual biomass monitoring in topographically complex mountainous forests. The resulting 25 m resolution biomass product has the potential to provide spatially detailed information for forest monitoring and may support applications in carbon accounting and forest management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Doubling protected land area may be inefficient at preserving the extent of undeveloped land and could cause substantial regional shifts in land use

Abstract Projection of land use and land‐cover change is highly uncertain yet drives critical estimates of carbon emissions, climate change, and food and bioenergy production. We use new, spatially explicit land availability data in conjunction with a model sensitivity analysis to estimate the effects of additional land protection on land use and land cover. The land availability data include protected land and agricultural suitability and is incorporated into the Moirai land data system for initializing the Global Change Analysis Model. Overall, decreasing land availability is relatively inefficient at preserving undeveloped land while having considerable regional land‐use impacts. Current amounts of protected area have little effect on land and crop production estimates, but including the spatial distribution of unsuitable (i.e., unavailable) land dramatically shifts bioenergy production from high northern latitudes to the rest of the world, compared with uniform availability. This highlights the importance of spatial heterogeneity in understanding and managing land change. Approximately doubling the current protected area to emulate a 30% protected area target may avoid land conversion by 2050 of less than half the newly protected extent while reducing bioenergy feedstock land by 10.4% and cropland and grazed pasture by over 3%. Regional bioenergy land may be reduced (increased) by up to 46% (36%), cropland reduced by up to 61%, pasture reduced by up to 100%, and harvested forest reduced by up to 35%. Only a few regions show notable gains in some undeveloped land types of up to 36%. Half of the regions can reach the target using only unsuitable land, which would minimize impacts on agriculture but may not meet conservation goals. Rather than focusing on an area target, a more robust approach may be to carefully select newly protected land to meet well‐defined conservation goals while minimizing impacts to agriculture.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

An Improved Energy Deposition Model in MPACT and Explicit Heat Generation Coupling with CTF

The default energy deposition model in the CASL neutronics code MPACT assumes all fission energy is deposited locally in fuel rods. Furthermore, equilibrium delayed energy release is assumed for both steady-state and transient calculations. These approximations limit the accurate representation of the heat generation distribution in space and its variations over time, which are essential for power distribution and thermal-hydraulic coupling calculations. In this paper, an improved energy deposition model is presented in both the spatial and time domains. Spatially, the energy deposition through fission, neutron capture, and slowing-down reactions are explicitly modeled to account for the heat generation from all regions of a reactor core, and a gamma smearing scheme is developed that utilizes the gamma sources from neutron fission and capture. In the time domain, the delayed energy release is modeled by solving an additional equation of delayed heat emitters, similar to the equation of delayed neutron precursors.To allow the explicit heat generation coupling, the interfaces between MPACT and CTF were updated to transfer separate heat sources for different material regions (fuel, clad, moderator and guide tube). The results show that the distributions of the energy deposition between MPACT and MCNP agree very well for various 2-D assembly and quarter-core problems without TH feedback. The MPACT/CTF coupled calculation for the hot full power quarter-core case exhibited a reduced peak pin power by 2.3% and a reduced peak fuel centerline temperature by 17 K when using the explicit energy deposition and heat transfer. The new model also shows a maximum 100 pcm keff effect on assembly depletion problems and an increased overall energy release by 7% in a PWR reactivity-initiated accident (RIA) problem.

Liu, Yuxuan↗
DOE OSTI2020

Interpenetrating Plasma Simulations

This project investigated multifluid models for simulating interpenetrating plasmas. These models solve the Euler equations for each plasma stream along with an appropriate model for electron dynamics. The project developed a novel code, Euclid, that simulates these models using high-order conservative spatial discretization methods and explicit and semi-implicit time integration methods. Euclid was used to solve several multifluid test cases representative of plasma dynamics in high-energy density physics experiments. Kinetic simulations of interpenetrating flows using the LOKI code revealed novel physical effects of the ion-ion streaming instability. A new multi-flow collision operator was developed for the LOKI code. A simple model of the ion-ion streaming instability was also developed. Several presentations were given on the project results; three journal articles have been published, one has been submitted, and two are planned.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Geospatial assessment of the economic opportunity for reforestation in Maryland, USA

Afforestation and reforestation have the potential to provide effective climate mitigation through forest carbon sequestration. Strategic reforestation activities, which account for both carbon sequestration potential (CSP) and economic opportunity, can provide attractive options for policymakers who must manage competing social and environmental goals. In particular, forest carbon pricing can incentivize reforestation on private land, but this may require landholders to forego other profits. Here, we utilize an ambitious geospatial approach to quantify economic opportunities for reforestation in the state of Maryland (USA) based on high-resolution remoting sensing, ecosystem modeling, and economic analysis. Our results identify spatially-explicit areas of economic opportunity where the potential revenue from forest carbon outcompetes the expected profit of existing cropland at the hectare scale. Specifically, we find that under a baseline economic scenario of 20 dollars per ton of carbon (5% rental rate) and decadal average crop profitability, a transition to forest on agricultural land would be more profitable than 23.2% of cropland in Maryland under a 20-year land-use commitment. Accounting for variations in carbon and crop pricing, 5.5% to 55.4% of cropland would be immediately outcompeted by expected forest carbon revenue, with the potential for an additional 0.5% to 10.6% of outcompeted cropland within 20 years. Under the baseline economic scenario, an annual allocation of $5.8 million towards a carbon rental program could protect 6.93 Tg C (2.2% of the state’s total CSP) on reforested croplands. This moderate yearly cost is equal to 9.7% of Maryland’s average annual auction proceeds from participation in the Regional Greenhouse Gas Initiative (between 2014-2018), and 19.3% of the average annual subsidy payments for corn, soy, and wheat allocated over the same period. This methodological approach may be useful for state governments, not-for-profit organizations, or regional climate initiatives interested in identifying strategic areas for reforestation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Probabilistic Groundwater Modeling of the RDX Plume at Los Alamos National Laboratory to Support Risk Assessment - 20359

A plume of the contaminant hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) with concentrations greater than the New Mexico tap water drinking standard (9.66 ppb) is present in the regional aquifer near the southwestern boundary of Los Alamos National Laboratory (LANL). A risk assessment is performed for exposure to regional aquifer groundwater, with long-term predictions of RDX concentrations provided by a calibrated, probabilistic, numerical fate and transport model. The structure of the model is hierarchical, with the RDX Regional Aquifer (RA) groundwater model acting as the primary tool for analysis of downgradient RDX concentrations. The RA model is deeply informed by the conceptual site model (CSM) and is calibrated using site RDX concentration data and hydraulic head measurements, along with other analyses. The model is calibrated using data through December 2019. Where data are scarce other lines of evidence are used to inform inputs, including the multiphase RDX Vadose Zone (VZ) model and Pipe and Disk analytical screening tool. Model inputs are described with informative prior distributions using a robust approach to development that incorporates all available lines of evidence for every parameter used as an input to the RA model. Calibration is performed using a classical nonlinear optimization routine, which is then used to initialize a Bayesian calibration. The Bayesian calibration constrains the uncertainty in the classical calibration, ultimately providing posterior distributions for all model parameters. The challenges of the calibration include high-dimensional parameter space, including spatially heterogeneous hydraulic conductivities, comparatively sparse data, and low RDX concentrations. Posterior parameter distributions developed in the calibration process are then used for stochastic predictive model runs into the future. The result of the forward runs is spatially and temporally explicit estimates of head and concentration with uncertainty at all points in space and time. The probabilistic modeling approach presented here includes innovative computational and statistical methods that leverage high-performance computing (HPC) resources. It makes use of prior modeling work performed at the chromium plume site in the central LANL area, with extensive updates. The risk assessment will ultimately be used to support decision-making at the site, using multiple appropriately weighted sources of information, as well as uncertainty, sensitivity, and value-of-information analyses. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Global spatially explicit carbon emissions from land-use change over the past six decades (1961–2020)

Land-use change (LUC) is a major source of global anthropogenic carbon emissions that contribute to climate change. However, current estimates of LUC-induced emissions vary widely with differences in data and models, making it challenging to identify and manage emissions hotspots by location and LUC activity. Here, we estimated spatially explicit carbon fluxes associated with global major gross land-use transitions based on a new bookkeeping model (i.e., LUCE). Between 1961 and 2020, LUC induced 215 and 142 Pg C of global carbon emissions and removals, respectively, resulting in average net emissions (E LUC ) of 1.21 Pg C year −1 . Latin America, Southeast Asia, and Sub-Saharan Africa dominated global E LUC with 69% of 60-year cumulative emissions, or about 86% between 2001 and 2020. Forest-related LUC activities (e.g., deforestation, reforestation) contributed largely to both regional and global carbon fluxes. Our findings shed new light on identifying LUC-induced emissions hotspots and managing land for climate mitigation and conservation.

bookkeeping model↗