Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Land 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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Beyond Population and Environment: Household Life Cycle Demography and Land Use Allocation among Small Farm Colonists in the Amazon

Most research featuring demographic factors in environmental change has focused on processes operating at the level of national or global populations. This paper focuses on household-level demographic life cycles among colonists in the Amazon, and evaluates the impacts on land use allocation. The analysis goes beyond prior research by including a broader suite of demographic variables, and by simultaneously assessing their impacts on multiple land uses with different economic and ecological implications. We estimate a system of structural equations that accounts for endogeneity among land uses, and the findings indicate stronger demographic effects than previous work. These findings bear implications for modeling land use, and the place of demography in environmental research.

Perz, Stephen G.↗
DOE OSTI2024

Advancing earth system model calibration: a diffusion-based method

Understanding of climate impact on ecosystems globally requires site-specific model calibration. Here we introduce a novel diffusion-based uncertainty quantification (DBUQ) method for efficient model calibration. DBUQ is a score-based diffusion model that leverages Monte Carlo simulation to estimate the score function and evaluates a simple neural network to quickly generate samples for approximating parameter posterior distributions. DBUQ is stable, efficient, and can effectively calibrate the model given diverse observations, thereby enabling rapid and site-specific model calibration on a global scale. This capability significantly advances Earth system modeling and our understanding of climate impacts on Earth systems. We demonstrate DBUQ's capability in E3SM land model calibration at the Missouri Ozark AmeriFlux forest site. Both synthetic and real-data applications indicate that DBUQ produces accurate parameter posterior distributions similar to those generated by Markov Chain Monte Carlo sampling but with 30X less computing time. This efficiency marks a significant stride in model calibration, paving the way for more effective and timely climate impact analyses.

Liu, Yanfang↗
DOE OSTI2020

Simulating coupled surface–subsurface flows with ParFlow v3.5.0: capabilities, applications, and ongoing development of an open-source, massively parallel, integrated hydrologic model

Surface flow and subsurface flow constitute a naturally linked hydrologic continuum that has not traditionally been simulated in an integrated fashion. Recognizing the interactions between these systems has encouraged the development of integrated hydrologic models (IHMs) capable of treating surface and subsurface systems as a single integrated resource. IHMs are dynamically evolving with improvements in technology, and the extent of their current capabilities are often only known to the developers and not general users. This article provides an overview of the core functionality, capability, applications, and ongoing development of one open-source IHM, ParFlow. ParFlow is a parallel, integrated, hydrologic model that simulates surface and subsurface flows. ParFlow solves the Richards equation for three-dimensional variably saturated groundwater flow and the two-dimensional kinematic wave approximation of the shallow water equations for overland flow. The model employs a conservative centered finite-difference scheme and a conservative finite-volume method for subsurface flow and transport, respectively. ParFlow uses multigrid-preconditioned Krylov and Newton–Krylov methods to solve the linear and nonlinear systems within each time step of the flow simulations. The code has demonstrated very efficient parallel solution capabilities. ParFlow has been coupled to geochemical reaction, land surface (e.g., the Common Land Model), and atmospheric models to study the interactions among the subsurface, land surface, and atmosphere systems across different spatial scales. This overview focuses on the current capabilities of the code, the core simulation engine, and the primary couplings of the subsurface model to other codes, taking a high-level perspective.

58 GEOSCIENCES↗
DOE OSTI2021

Modeling Profiles of Micrometeorological Variables in a Tropical Premontane Rainforest Using Multi-Layered CLM (CLM-ML)

This study updates the multi-layered Community Land Model (CLM-ml) for hillslopes and compares predictions from against observations collected in tropical montane rainforest, Costa Rica. Modifications are made in order to capture a wider array of vertical leaf area distributions, predict CO 2 profiles, account for soil respiration, and adjust wind forcings for difficult topographic settings. Test results indicate that the modified multi-layer CLM model can successfully replicate the shape of various micrometeorological profiles (humidity, CO 2 , temperature, and wind speed) under the canopy. In the single-layer models (CLM4.5 and CLM5), excessive day-to-night differences in leaf temperature and leaf wetness were originally noted, but CLM-ml significantly improved these issues, decreasing the amplitudes of diurnal cycles by 67% and 47%. Sub-canopy considerations, such as canopy shapes and turbulent transfer parameters, also played a significant role in model performance. More importantly, unlike single layer models, the results that CLM-ml produces can be compared to variables measured within the canopy to provide far more detailed diagnostic information. Further observations and model developments, aimed at reflecting surface heterogeneity, will be necessary to adequately capture the complexity and the features of the tropical montane rainforest.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Hourly land-use regression models based on low-cost PM monitor data

Highlights: • A network of 23 low-cost PM monitors was deployed across an urban area. • Land use regression models predict the hourly intra-urban PM variations. • Coefficients of determination for hourly PM predictions ranged 0.66–0.76. • Low-cost PM monitor data need to be carefully evaluated before use. • The combination of low-cost monitors and modeling provides useful PM estimates. Land-use regression (LUR) models provide location and time specific estimates of exposure to air pollution and thereby improve the sensitivity of health effects models. However, they require pollutant concentrations at multiple locations along with land-use variables. Often, monitoring is performed over short durations using mobile monitoring with research-grade instruments. Low-cost PM monitors provide an alternative approach that increases the spatial and temporal resolution of the air quality data. LUR models were developed to predict hourly PM concentrations across a metropolitan area using PM concentrations measured simultaneously at multiple locations with low-cost monitors. Monitors were placed at 23 sites during the 2015/16 heating season. Monitors were externally calibrated using co-located measurements including a reference instrument (GRIMM particle spectrometer). LUR models for each hour of the day and weekdays/weekend days were developed using the deletion/substitution/addition algorithm. Coefficients of determination for hourly PM predictions ranged from 0.66 and 0.76 (average 0.7). The hourly-resolved LUR model results will be used in epidemiological studies to examine if and how quickly, increases in ambient PM concentrations trigger adverse health events by reducing the exposure misclassification that arises from using less time resolved exposure estimates.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2021

Validation Assessment for the Soil Moisture Active Passive (SMAP) Level 4 Carbon (L4_C) Data Product Version 5

The post-launch Cal/Val phase of the SMAP mission is guided by two primary objectives for each science product team: 1) to calibrate, verify, and improve the performance of the science algorithms, and 2) validate accuracies of the science data products as specified in the SMAP Level-1 mission science requirements. Algorithm science and product maintenance activities during the SMAP extended mission phase have also involved periodic algorithm calibration and product refinements to maintain or enhance product consistency and performance as well as science utility. This report provides an assessment of the latest (Version 5) SMAP Level 4 Carbon (L4_C) product. The L4_C Version 5 (v5) global record now spans more than six years (March 2015 – present) of SMAP operations and has benefited from five major reprocessing updates to the operational product. These reprocessing events and L4_C product release updates have incorporated various algorithm refinements and calibration adjustments to account for similar refinements to the upstream GEOS land model assimilation system, SMAP brightness temperatures, and MODIS vegetation inputs used for L4_C processing. The SMAP L4_C algorithms utilize a terrestrial carbon flux model informed by daily surface and root zone soil moisture information contributed from the SMAP Level 4 Soil Moisture (L4_SM) product along with optical remote sensing-based (e.g. MODIS-based) land cover and canopy fractional photosynthetic active radiation (fPAR), and other ancillary biophysical data. The carbon flux model estimates global daily net ecosystem CO2 exchange (NEE) and the component carbon fluxes, namely, vegetation gross primary production (GPP) and soil heterotrophic respiration (Rh). Other L4_C product elements include surface (~0-5 cm depth) soil organic carbon (SOC) stocks and associated environmental constraints to these processes, including soil moisture-related controls on GPP and ecosystem respiration (Kimball et al. 2014, Jones et al. 2017). The L4_C product addresses SMAP carbon cycle science objectives by: 1) providing a direct link between terrestrial carbon fluxes and underlying freeze/thaw and soil moisture-related constraints to these processes, 2) documenting primary connections between terrestrial water, energy and carbon cycles, and 3) improving understanding of terrestrial carbon sink activity. The SMAP L4_C algorithms and operational product are mature and at a CEOS Validation Stage 4 level (Jackson et al. 2012) based on extensive validation of the multi-year record against a diverse array of independent benchmarks, well characterized global performance, and systematic refinements gained from five major reprocessing events. There are no Level-1 mission science requirements for the L4_C product; however, self-imposed requirements have been established focusing on NEE as the primary product field for validation, and on demonstrating L4_C accuracy and success in meeting product science requirements (Jackson et al. 2012). The other L4_C product fields also have strong utility for carbon science applications (e.g., Liu et al. 2019, Endsley et al. 2020); however, analysis of these other fields is considered secondary relative to primary validation activities focusing on NEE. The L4_C targeted accuracy requirements are to meet or exceed a mean unbiased root-mean-square error (ubRMSE, or standard deviation of the error) for NEE of 1.6 g C m-2 d-1 and 30 g C m-2 yr-1, emphasizing northern (≥45°N) boreal and arctic ecosystems; this accuracy is similar to that of tower eddy covariance measurement-based observations (Baldocchi 2008). Methods used for the latest v5 L4_C product performance and validation assessment have been established from the SMAP Cal/Val plan and previous studies (Jackson et al. 2012, Jones et al. 2017) and include: 1) consistency evaluations of the product fields against earlier product releases (version 4 or earlier); 2) comparisons of daily carbon flux estimates with independent tower eddy 4covariance measurement-based daily carbon (CO2) flux observations from core tower validation sites (CVS); and 3) consistency checks against other global carbon products, including soil carbon inventory records, global GPP records derived from tower observation upscaling methods, and satellite-based observations of canopy solar induced chlorophyll fluorescence (SIF) as a surrogate for GPP. Metrics used to evaluate relative agreement between L4_C product fields and observational benchmarks include correlation (r-value), RMSE differences, bias and model sensitivity diagnostics. Following these validation criteria, the present report provides a validation assessment of the latest L4_C product release (v5). Detailed descriptions of the L4_C algorithm and additional global product accuracy and performance results are given elsewhere (Jones et al. 2017, Endsley et al. 2020). The v5 L4_C product replaces earlier product versions and continues to show: (i) accuracy and performance levels meeting or exceeding SMAP L4_C science requirements; (ii) improvement over the previous product version (version 4); and (iii) suitability for a diversity of science applications. Example L4_C applications from the recent literature include clarifying environmental trends and controls on the northern terrestrial carbon sink (Liu et al. 2019), diagnosing drought-related impacts on ecosystem productivity (Li et al. 2020), and regional monitoring of cropland conditions for projecting annual yields (Wurster et al. 2020). 2EXPECTED L4_C ALGORITHM AND PRODUCTPERFORMANCE The L4_C algorithm performance, including variance and uncertainty estimates of model outputs, was determined during the mission pre-launch phase through spatially explicit model sensitivity studies using available model inputs similar to those currently being used for operational production and evaluating the resulting model simulations over the observed range of northern (≥45 °N) and global conditions (Kimball et al. 2012, Entekhabi et al. 2014). The L4_C algorithm options were also evaluated during the mission prelaunch phase, including deriving canopy fPAR from lower order NDVI (Normalized Difference Vegetation Index) inputs in lieu of using MODIS (MOD15) fPAR; and including an explicit model representation of boreal fire disturbance recovery impacts. These results indicated that the L4_C accuracy requirements (i.e., NEE ubRMSE ≤ 30 g C m-2 yr-1or ≤ 1.6 g C m-2 d-1) could be met from the baseline algorithms over more than 82% and 89% of global and northern vegetated land areas, respectively (Yi et al. 2013, Kimball et al. 2014). The global L4_C algorithm error budget for NEE derived during the mission prelaunch phase indicated that the estimated NEE ubRMSE uncertainty is proportional to GPP and is therefore larger in higher biomass productivity areas, including forests and croplands (Kimball et al. 2014). Likewise, NEE ubRMSE uncertainty is expected to be lower in less-productive areas, including grasslands and shrublands. Expected model NEE ubRMSE levels were also generally within targeted accuracy levels for characteristically less-productive boreal and Arctic biomes, even though relative model error as a proportion of total productivity (NEE RMSE / GPP) may be large in these areas. The estimated NEE uncertainty was lower than expected in some warmer tropical high biomass productivity areas (e.g. Amazon rainforest) because of reduced low temperature and moisture constraints to the L4_C respiration calculations so that the bulk of model uncertainty is contributed by GPP in these areas. Model NEE uncertainty in the African Congo was estimated to

SMAP↗
DOE OSTI2013

Implementation and Initial Evaluation of the Glimmer Community Ice Sheet Model in the Community Earth System Model

The Glimmer Community Ice Sheet Model (Glimmer-CISM) has been implemented in the Community Earth System Model (CESM). Glimmer-CISM is forced by a surface mass balance (SMB) computed in multiple elevation classes in the CESM land model and downscaled to the ice sheet grid. Ice sheet evolution is governed by the shallow-ice approximation with thermomechanical coupling and basal sliding. This paper describes and evaluates the initial model implementation for the Greenland Ice Sheet (GIS). The ice sheet model was spun up using the SMB from a coupled CESM simulation with preindustrial forcing. The model's sensitivity to three key ice sheet parameters was explored by running an ensemble of 100 GIS simulations to quasi equilibrium and ranking each simulation based on multiple diagnostics. With reasonable parameter choices, the steady-state GIS geometry is broadly consistent with observations. The simulated ice sheet is too thick and extensive, however, in some marginal regions where the SMB is anomalously positive. The top-ranking simulations were continued using surface forcing from CESM simulations of the twentieth century (1850–2005) and twenty-first century (2005–2100, with RCP8.5 forcing). In these simulations the GIS loses mass, with a resulting global-mean sea level rise of 16 mm during 1850–2005 and 60 mm during 2005–2100. This mass loss is caused mainly by increased ablation near the ice sheet margin, offset by reduced ice discharge to the ocean. Projected sea level rise is sensitive to the initial geometry, showing the importance of realistic geometry in the spun-up ice sheet.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

A modeling study of irrigation effects on global surface water and groundwater resources under a changing climate

Abstract This study investigates the effects of irrigation on global water resources by performing and analyzing Community Land Model 4.0 (CLM4) simulations driven by downscaled/bias‐corrected historical simulations and future projections from five General Circulation Models (GCMs). For each climate scenario, three sets of numerical experiments were performed: (1) a CTRL experiment in which all crops are assumed to be rainfed; (2) an IRRIG experiment in which the irrigation module is activated using surface water (SW) to feed irrigation; and (3) a PUMP experiment in which a groundwater pumping scheme coupled with the irrigation module is activated for conjunctive use of surface water and groundwater (GW) for irrigation. The parameters associated with irrigation and groundwater pumping are calibrated based on a global inventory of census‐based water use compiled by the Food and Agricultural Organization (FAO). Our results suggest that irrigation could lead to two major effects: SW (GW) depletion in regions with irrigation primarily fed by SW (GW), respectively. Furthermore, irrigation depending primarily on SW tends to have larger impacts on low‐flow than high‐flow conditions, suggesting increased vulnerability to drought. By the end of the 21st century, combined effect of increased irrigation water demand and amplified temporal‐spatial variability of water supply may lead to severe local water scarcity for irrigation. Regionally, irrigation has the potential to aggravate/alleviate climate‐induced changes of SW/GW although such effects are negligible when averaged globally. Our study highlights the need to account for irrigation effects and sources in assessing regional climate change impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Investigating the impact of aerosol deposition on snowmelt over the Greenland Ice Sheet using a large-ensemble kernel

Accelerating surface melt on the Greenland Ice Sheet (GrIS) has led to a doubling of Greenland's contribution to global sea level rise during recent decades. Black carbon (BC), dust, and other light-absorbing impurities (LAIs) darken the surface and enhance snowmelt by boosting the absorption of solar energy. It is therefore important for coupled aerosol–climate and ice sheet models to include snow darkening effects from LAI, and yet most do not. In this study, we conduct several thousand simulations with the Community Land Model (CLM) component of the Community Earth System Model (CESM) to characterize changes in melt runoff due to variations in the amount, timing, and nature (wet or dry) of BC deposition on the GrIS. From this large matrix of simulations, we develop a kernel relating runoff to the location, month, year (from 2006 to 2015), and magnitudes of BC concentration within precipitation and dry deposition flux. BC deposition during June–August causes the largest increase in annually integrated runoff, but winter deposition events also exert large (roughly half as great) runoff perturbations due to reexposure of impurities at the snow surface during summer melt. Current BC deposition fluxes simulated with the atmosphere component of CESM induce a climatological-mean increase in GrIS-wide runoff of ∼ 8 Gt yr −1 , or +6.8 % relative to a paired simulation without BC deposition. We also provide linear equations that relate the increase in total runoff to GrIS-wide wet and dry BC deposition fluxes. It is our hope that the runoff kernel and simple equations provided here can be used to extend the utility of state-of-the-art aerosol models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Regionalization of land use impact models for life cycle assessment: Recommendations for their use on the global scale and their applicability to Brazil

Life Cycle Assessment (LCA) is the main technique for evaluate the environmental impacts of product life cycles. A major challenge in the field of LCA is spatial and temporal differentiation in Life Cycle Impact Assessment (LCIA) methods, especially impacts resulting from land occupation and land transformation. Land use characterization modeling has advanced considerably over the last two decades and many approaches have recently included crucial aspects such as geographic differentiation. Nevertheless, characterization models have so far not been systematically reviewed and evaluated to determine their applicability to South America. Given that Brazil is the largest country in South America, this paper analyzes the main international characterization models currently available in the literature, with a view to recommending regionalized models applicable on a global scale for land use life cycle impact assessments, and discusses their feasibility for regionalized assessment in Brazil. The analytical methodology involves classification based on the following criteria: midpoint/endpoint approach, scope of application, area of data collection, biogeographical differentiation, definition of recovery time and reference situation; followed by an evaluation of thirteen scientific robustness and environmental relevance subcriteria. The results of the scope of application are distributed among 25% of the models developed for the European context, and 50% have a global scope. There is no consensus in the literature about the definition of parameters such biogeographical differentiation and reference situation, and our review indicates that 35% of the models use ecoregion division while 40% use the concept of potential natural vegetation. Four characterization models show high scores in terms of scientific robustness and environmental relevance. These models are recommended for application in land use life cycle impact assessments, and also to serve as references for the development or adaptation of regional methodological procedures for Brazil. - Highlights: • A discussion is made on performing regionalized impact assessments using spatial differentiation in LCA. • A review is made of 20 characterization models for land use impacts in Life Cycle Impact Assessment. • Four characterization models are recommended according to different land use impact pathways for application in Brazil.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2010

Comparison of Satellite-Derived TOA Shortwave Clear-Sky Fluxes to Estimates from GCM Simulations Constrained by Satellite Observations of Land Surface Characteristics

Clear-sky, upwelling shortwave flux at the top of the atmosphere (S(sub TOA raised arrow)), simulated using the atmospheric and land model components of the Community Climate System Model 3 (CCSM3), is compared to corresponding observational estimates from the Clouds and Earth's Radiant Energy System (CERES) sensor. Improvements resulting from the use of land surface albedo derived from Moderate Resolution Imaging Spectroradiometer (MODIS) to constrain the simulations are also examined. Compared to CERES observations, CCSM3 overestimates global, annual averaged S(sub TOA raised arrow) over both land and oceans. However, regionally, CCSM3 overestimates S(sub TOA raised arrow) over some land and ocean areas while underestimating it over other sites. CCSM3 underestimates S(sub TOA raised arrow) over the Saharan and Arabian Deserts and substantial differences exist between CERES observations and CCSM3 over agricultural areas. Over selected sites, after using groundbased observations to remove systematic biases that exist in CCSM computation of S(sub TOA raised arrow), it is found that use of MODIS albedo improves the simulation of S(sub TOA raised arrow). Inability of coarse resolution CCSM3 simulation to resolve spatial heterogeneity of snowfall over high altitude sites such as the Tibetan Plateau causes overestimation of S(sub TOA raised arrow) in these areas. Discrepancies also exist in the simulation of S(sub TOA raised arrow) over ocean areas as CCSM3 does not account for the effect of wind speed on ocean surface albedo. This study shows that the radiative energy budget at the TOA is improved through the use of MODIS albedo in Global Climate Models.

Anantharaj, Valentine G.↗
NASA NTRS2017

Orion Multi-Purpose Crew Vehicle Solving and Mitigating the Two Main Cluster Pendulum Problem

The Orion Multi-purpose Crew Vehicle (MPCV) Orion spacecraft will return humans from beyond earth's orbit, including Mars and will be required to land 20,000 pounds of mass safely in the ocean. The parachute system nominally lands under 3 main parachutes, but the system is designed to be fault tolerant and land under 2 main parachutes. During several of the parachute development tests, it was observed that a pendulum, or swinging, motion could develop while the Crew Module (CM) was descending under two parachutes. This pendulum effect had not been previously predicted by modeling. Landing impact analysis showed that the landing loads would double in some places across the spacecraft. The CM structural design limits would be exceeded upon landing if this pendulum motion were to occur. The Orion descent and landing team was faced with potentially millions of dollars in structural modifications and a severe mass increase. A multidisciplinary team was formed to determine root cause, model the pendulum motion, study alternate canopy planforms and assess alternate operational vehicle controls & operations providing mitigation options resulting in a reliability level deemed safe for human spaceflight. The problem and solution is a balance of risk to a known solution versus a chance to improve the landing performance for the next human-rated spacecraft.

Ali, Yasmin↗
DOE OSTI2020

Linear two-pool models are insufficient to infer soil organic matter decomposition temperature sensitivity from incubations

Abstract Terrestrial carbon (C)-climate feedbacks depend strongly on how soil organic matter (SOM) decomposition responds to temperature. This dependency is often represented in land models by the parameter Q 10 , which quantifies the relative increase of microbial soil respiration per 10 °C temperature increase. Many studies have conducted paired laboratory soil incubations and inferred “active” and “slow” pool Q 10 values by fitting linear two-pool models to measured respiration time series. Using a recently published incubation study (Qin et al. in Sci Adv 5(7):eaau1218, 2019) as an example, here we first show that the very high parametric equifinality of the linear two-pool models may render such incubation-based Q 10 estimates unreliable. In particular, we show that, accompanied by the uncertain initial active pool size, the slow pool Q 10 can span a very wide range, including values as high as 100, although all parameter combinations are producing almost equally good model fit with respect to the observations. This result is robust whether or not interactions between the active and slow pools are considered (typically these interactions are not considered when interpreting incubation data, but are part of the predictive soil carbon models). This very large parametric equifinality in the context of interpreting incubation data is consistent with the poor temporal extrapolation capability of linear multi-pool models identified in recent studies. Next, using a microbe-explicit SOM model (RESOM), we show that the inferred two pools and their associated parameters (e.g., Q 10 ) could be artificial constructs and are therefore unreliable concepts for integration into predictive models. We finally discuss uncertainties in applying linear two-pool (or more generally multiple-pool) models to estimate SOM decomposition parameters such as temperature sensitivities from laboratory incubations. We also propose new observations and model structures that could enable better process understanding and more robust predictive capabilities of soil carbon dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Accelerating the spin-up of the coupled carbon and nitrogen cycle model in CLM4

The commonly adopted biogeochemistry spin-up process in an Earth system model (ESM) is to run the model for hundreds to thousands of years subject to periodic atmospheric forcing to reach dynamic steady state of the carbon–nitrogen (CN) models. A variety of approaches have been proposed to reduce the computation time of the spin-up process. Significant improvement in computational efficiency has been made recently. However, a long simulation time is still required to reach the common convergence criteria of the coupled carbon–nitrogen model. A gradient projection method was proposed and used to further reduce the computation time after examining the trend of the dominant carbon pools. The Community Land Model version 4 (CLM4) with a carbon and nitrogen component was used in this study. From point-scale simulations, we found that the method can reduce the computation time by 20–69% compared to one of the fastest approaches in the literature. We also found that the cyclic stability of total carbon for some cases differs from that of the periodic atmospheric forcing, and some cases even showed instability. Close examination showed that one case has a carbon periodicity much longer than that of the atmospheric forcing due to the annual fire disturbance that is longer than half a year. The rest was caused by the instability of water table calculation in the hydrology model of CLM4. The instability issue is resolved after we replaced the hydrology scheme in CLM4 with a flow model for variably saturated porous media.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Matrix Approach to Accelerate Spin-Up of CLM5

Numerical models have been developed to investigate and understand responses of biogeochemical cycle to global changes. Steady state, when a system is in dynamic equilibrium, is generally required to initialize these model simulations. However, the spin-up process that is used to achieve steady state pose a great burden to computational resources, limiting the efficiency of global modeling analysis on biogeochemical cycles. This study introduces a new Semi-Analytical Spin-Up (SASU) to tackle this grand challenge. We applied SASU to Community Land Model version 5 and examined its computational efficiency and accuracy. At the Brazil site, SASU is computationally 7 times more efficient than (or saved up to 86% computational cost in comparison with) the traditional native dynamics (ND) spin-up to reach the same steady state. Globally, SASU is computationally 8 times more efficient than the accelerated decomposition spin-up and 50 times more efficient than ND. In summary, SASU achieves the highest computational efficiency for spin-up on site and globally in comparison with other spin-up methods. It is generalizable to wide biogeochemical models and thus makes computationally costly studies (e.g., parameter perturbation ensemble analysis and data assimilation) possible for a better understanding of biogeochemical cycle under climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Leveraging High-resolution Molecular Composition of Soil Organic Matter to Enhance Carbon Cycling Modeling

Soils store more carbon than the atmosphere and vegetation combined, yet Earth system models still struggle to predict how this vast reservoir will respond to environmental change. A central limitation is that most soil biogeochemical models represent organic matter using bulk conceptual pools or chemically homogeneous fractions, preventing direct use of rapidly expanding molecular-scale datasets. Here we develop and test a new soil decomposition framework that explicitly integrates high-resolution information on organic matter composition. First, we construct a molecularly informed litter decomposition module in which plant inputs are partitioned into five functional compound classes—carbohydrates, proteins, lignin-like aromatics, lipids, and carbonyls—using a molecular mixing model calibrated to solid-state 13 C Nuclear Magnetic Resonance (NMR) spectra. Class-specific kinetics, lignin-dependent physical protection, and substrate-driven microbial carbon use efficiency allow the module to capture metabolic tradeoffs associated with enzyme production and nutrient limitation. We then embed this litter module within a microbially explicit whole-soil model that tracks the transformation of these compound classes through particulate organic matter, dissolved organic matter, mineral-associated organic matter, and microbial biomass. High-resolution Fourier Transform Ion Cyclotron Resonance mass spectrometry (FTICR-MS) data are used to link internal pools to measurable soil organic matter fractions and to constrain key process parameters. Applications at soil-core and ecosystem scales demonstrate that the new model reproduces observed soil respiration dynamics while providing mechanistic attribution of CO 2 fluxes to specific chemical classes and pools. Compared to existing frameworks such as the Community Land Model soil biogeochemistry module and the Millennial model, our approach maintains competitive predictive skill while substantially improving interpretability and opportunities for data–model integration. This work illustrates a viable pathway for leveraging molecular-scale observations to reduce structural uncertainty in soil carbon–climate feedback projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

CLM5 Simulations of Soil Moisture and Grain Carbon for CONUS at 0.125 degrees

This dataset provides 0.125° gridded simulations of soil moisture and crop grain carbon for the Contiguous United States (CONUS), generated using the Community Land Model version 5 (CLM5) with the biogeochemistry module enabled. The data covers a historical baseline (1980–2015) and mid-century future projections (2020–2055). Future projections are organized into two sets of scenarios to distinguish the impacts of different drivers: (1) Atmospheric Only (ATM-only): These scenarios apply future atmospheric forcings while holding land use and land cover change (LULCC) at historical baseline levels. The atmospheric forcings represent moderately versus severely hotter/drier atmospheric conditions (dynamically downscaled perturbed thermodynamics simulations based on CMIP6 SSP245 and SSP585 warming signals), each with cooler versus hotter Earth System Model temperature sensitivity instantiations. These scenarios are identified in the file names as atm45cooler, atm45hotter, atm85cooler, and atm85hotter. (2) Coupled Atmospheric and Land-Use (LAND+ATM): These scenarios apply future atmospheric forcing together with future LULCC by pairing atmospheric pathways with lower versus higher population/economic growth scenarios representing Shared Socioeconomic Pathways 3 and 5 (SSP3 and SSP5). These scenarios are identified in the file names as atm45cooler_ssp3, atm45hotter_ssp3, atm45cooler_ssp5, atm45hotter_ssp5, atm85cooler_ssp3, atm85hotter_ssp3, atm85cooler_ssp5, and atm85hotter_ssp5. Please refer to the README file for detailed information on file structure, variables, units, and data formats.

Yao, Lili [Pacific Northwest National Laboratory] ↗
DOE OSTI2026

Systematic Evaluation of Atmospheric Forcing, Surface Datasets, and Mesh Effects on Kilometer-Scale Land Surface and River Modeling

Earth system models are advancing toward kilometer-scale resolution to capture local climate impacts and extremes. High-resolution land and river modeling depends on multiple factors, including mesh, surface datasets, and atmospheric forcing, but their relative effects at kilometer scales remain unquantified. We evaluated five Energy Exascale Earth System Model land and river configurations over the Mid-Atlantic region using two mesh (1/8° structured versus variable-resolution unstructured mesh), two surface datasets (default versus newly developed), and three atmospheric forcings (NLDAS2, MSWX, GSWP). Evaluation against satellite, reanalysis, and in situ benchmarks across water, energy, and carbon cycles quantifies how these factors affect model performance. Forcing selection produces the largest bias reductions (12-99% across variables), followed by surface datasets (7-75%) and mesh (up to 21%). Forcing effects vary by variable, with MSWX reducing biases for snow water equivalent, evapotranspiration, albedo, temperature, and gross primary productivity, GSWP for snow cover and runoff, and NLDAS for soil moisture and streamflow. The use of newly developed surface datasets improves gross primary productivity (58% bias reduction) and evapotranspiration but increase soil moisture and albedo biases due to current modeling limitations. Variable-resolution unstructured mesh improves the simulation of small-basin streamflow through better capturing drainage networks, though mesh minimally affects other land variables. These findings provide important guidance for high-resolution modeling development and actionable science.

Land and River modeling↗