Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “land surface model”

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

Coupling Noah-Multiparameterization land-surface Model with Energy Research and Forecasting Model

The Energy Research and Forecasting (ERF) model is a high-performance atmospheric model built on the AMReX adaptive mesh refinement (AMR) framework, enabling efficient simulations on heterogeneous computing platforms that combine multicore processors with hardware accelerators. To support land–atmosphere interactions within ERF’s AMR-based environment, a land-surface model must be capable of operating directly on hierarchically refined meshes. In this work, we present a methodology for coupling the Fortran-based Noah-Multiparameterization (Noah-MP) land-surface model with ERF’s C++ codebase. Rather than rewriting Noah-MP, we construct a Fortran–C interoperability layer using CodeScribe, a tool that leverages large language models (LLMs) to automate the generation of interface code. CodeScribe applies structured prompting techniques to generate bindings that support efficient data exchange and function calls between ERF and Noah-MP. The coupling framework also incorporates AMR-aware data handling strategies, allowing NoahMP to operate seamlessly within ERF’s hierarchical mesh structure. This work provides a structured approach for integrating legacy Fortran models into modern C++-based modeling systems using LLM-assisted code generation.

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 OSTI2021

Land Surface Modeling 2.0 for agricultural climate change impact assessments

Focus Area: This white paper addresses DOE AI4ESP focus area #2 by providing a sketch blueprint for a next generation, hybrid AI/process-based global Land Surface Modeling (LSM) framework to improve projections of climate change impacts on the land surface system including agriculture and the hydrological cycle. Science Challenge: Climate change impacts on agriculture are highly uncertain: -50% to +150% global production changes for major grains under high-end climate change, for example. The core challenge in land surface system predictability is the large gap between the scale at which the relevant biological processes act and the scale at which the risks need to be assessed. Conventional empirical and process-based modeling approaches are insufficient. A new multi-scale modeling paradigm which employs AI/ML to learn from new streams of remote sensing data and targeted ’gene-to-global’ simulations is needed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Sensitivity of Near-Surface Variables in the RUC Land Surface Model in the Weather Research and Forecasting Model

In this study, we investigate the parametric sensitivity of near-surface variables, such as sensible heat flux, latent heat flux, ground heat flux, hub-height wind speed and land surface temperature, to the parameters used in the Rapid Update Cycle (RUC) land surface model (LSM) during a wintertime and summertime period. The model simulations are compared with observations collected from the second Wind Forecast Improvement Project (WFIP2) field campaign. The results suggest that parameters related to snow/ice and thermal processes can have significant impact on the simulated near-surface variables. Out of the 11 examined parameters, only 6 of them have considerable influences on the model behaviors and explain about 60 ~ 80 % of the estimated total variance of the simulated variables. In addition, the magnitude of the parametric sensitivity varies with season. For instance, parameters associated with snow/ice processes are dominant during the wintertime whereas those associated with thermal processes are more important during the summertime. Furthermore, the impact of the identified parameters on the simulated variables is highly related to the topography. There is a high degree of sensitivity to the parameter values over the slope region. This points out the importance of collecting field observations over steep areas to better quantity the appropriate values of key parameters. Overall, our findings provide a better understanding of the RUC LSM behavior associated with parameter uncertainties and can be used to improve the forecasting skill of land surface processes via calibration of the most uncertain model parameters.

17 WIND ENERGY↗
DOE OSTI2022

Evaluation and Optimization of Snow Albedo Scheme in Noah-MP Land Surface Model Using In Situ Spectral Observations in the Colorado Rockies

The Biosphere-Atmosphere Transfer Scheme (BATS) ground snow albedo algorithm is commonly used in land-surface models (LSM), weather forecasting and research applications. This study addresses key uncertainties in BATS simulated ground snow albedo within the Noah-MP LSM framework through evaluation and optimization of the Noah-MP BATS ground snow albedo formulation using 2-band (visible and near-infrared (NIR)) in situ albedo observations at Rocky Mountain field stations. The Noah-MP BATS ground snow albedo scheme is extremely sensitive to its input parameters. Namely, an ensemble generated by varying BATS input parameters within potentially plausible ranges provides an average daily range (maximum ensemble member minus minimum ensemble member) of ground snow albedo exceeding 0.45 in visible and NIR bands. Parameter optimization improves agreement between simulated and in situ observed ground snow albedo in visible, NIR and broadband spectrums. Importantly, optimized parameters result in reduced biases relative to observed fresh-snow albedo and better agreement with observed albedo decay. Our analysis across different sites supports that the optimized BATS ground snow albedo parameters are appropriate to transfer in space and time, at least within the region studied (the central-southern Rocky Mountains). The primary error source remaining after parameter optimization is that observed fresh-snow albedo is highly variable, particularly in the NIR spectrum, whereas BATS fresh-snow albedo is constant, an issue which requires further investigation. This study shows significant correlations between observed fresh-snow albedo and surface meteorological conditions (e.g., downward shortwave radiation and temperature) which can support future model development that attempts to include a time-varying formulation for fresh-snow albedo.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Extending a land-surface model with <i>Sphagnum</i> moss to simulate responses of a northern temperate bog to whole ecosystem warming and elevated CO<sub>2</sub>

Abstract. Mosses need to be incorporated into Earth system models to better simulate peatland functional dynamics under the changing environment. Sphagnum mosses are strong determinants of nutrient, carbon, and water cycling in peatland ecosystems. However, most land-surface models do not include Sphagnum or other mosses as represented plant functional types (PFTs), thereby limiting predictive assessment of peatland responses to environmental change. In this study, we introduce a moss PFT into the land model component (ELM) of the Energy Exascale Earth System Model (E3SM) by developing water content dynamics and nonvascular photosynthetic processes for moss. The model was parameterized and independently evaluated against observations from an ombrotrophic forested bog as part of the Spruce and Peatland Responses Under Changing Environments (SPRUCE) project. The inclusion of a Sphagnum PFT with some Sphagnum-specific processes in ELM allows it to capture the observed seasonal dynamics of Sphagnum gross primary production (GPP) albeit with an underestimate of peak GPP. The model simulated a reasonable annual net primary production (NPP) for moss but with less interannual variation than observed, and it reproduced aboveground biomass for tree PFTs and stem biomass for shrubs. Different species showed highly variable warming responses under both ambient and elevated atmospheric CO2 concentrations, and elevated CO2 altered the warming response direction for the peatland ecosystem. Microtopography is critical: Sphagnum mosses on hummocks and hollows were simulated to show opposite warming responses (NPP decreasing with warming on hummocks but increasing in hollows), and hummock Sphagnum was modeled to have a strong dependence on water table height. The inclusion of this new moss PFT in global ELM simulations may provide a useful foundation for the investigation of northern peatland carbon exchange, enhancing the predictive capacity of carbon dynamics across the regional and global scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Incorporating Microtopography in a Land Surface Model and Quantifying the Effect on the Carbon Cycle

Abstract Northern peatlands are a terrestrial carbon store, with an annual sink of 0.1 Pg C yr −1 and a total storage estimate of 547 Pg C. Northern peatlands are also major contributors of atmospheric methane. Most land surface models do not accurately represent peatland carbon emissions, partly because they do not represent the hydrologic cycle and/or microtopography adequately. Interactions between water table depth and microtopography in peatlands influence decomposition and modulate CO 2 and CH 4 fluxes. A modified version of the land surface component of the Energy Exascale Earth System Model, was recently created to represent the microtopography and hydrology of a raised dome bog in northern Minnesota, USA. In this study, three microtopographic parameters are analyzed in the modified version: hummock height, hummock‐hollow spacing, and percent hollow. Terrestrial laser scanning observations are used to set uncertainty bounds for these parameters. Our model experiment results suggest that carbon‐related quantities of interest (QOI) were typically the most sensitive to hummock height, and those QOI (especially net ecosystem exchange, NEE) were sensitive to interactions between parameters. Furthermore, NEE was most relatively influenced by microtopographic parameters in the model, varying by 35%. We found that increasing hummock height resulted in more C being stored in plant tissue and less in soil organic matter. This coincided with decreases in Sphagnum and increases in Picea and shrub net primary production. These results suggest that future studies may consider extending prognostic capabilities of carbon cycling by incorporating hummock hollow microtopography into earth system models.

58 GEOSCIENCES↗
DOE OSTI2025

A 1 km soil moisture dataset over eastern CONUS generated by assimilating SMAP data into the Noah-MP land surface model

An improved fine-scale soil moisture (SM) dataset at 1 km grid spacing, covering much of the eastern continental US, was generated by assimilating 9 km Soil Moisture Active Passive (SMAP) SM data into the v4.0.1 Noah-MP land surface model. With 12 ensemble members, the assimilation was carried out using the ensemble Kalman filter algorithm within NASA's Land Information System. The SM analysis for 2016 was fully validated against in situ observations from four different networks and compared with four other existing datasets. Results indicate that this SM analysis surpasses other datasets in top-layer SM distribution, including a machine-learning-based product, despite all SM estimates being less heterogeneous than observed. The analysis of anomalous errors suggests that large similarity in intrinsic errors is likely due to overlapping data sources among the selected SM datasets. More detailed evaluations were performed over two geographic areas. The observations collected by the Atmospheric Radiation Measurement facility in Oklahoma suggest that soil temperature and surface heat fluxes are concurrently simulated with good accuracy. Investigation into the 2016 southeastern US drought response further indicates drier conditions and higher evapotranspiration estimates compared to GLEAMv4.1. Notably, large errors are associated with grids having clay soil textures, underscoring the need for refined model treatments for specific soil types to further improve SM estimates. The dataset is publicly available on Zenodo at https://doi.org/10.5281/zenodo.14370563 (Tai et al., 2024).

Tai, Sheng-Lun [Pacific Northwest National Laborat↗
DOE OSTI2025

Incorporating the Acclimation of Photosynthesis and Leaf Respiration in the Noah‐MP Land Surface Model: Model Development and Evaluation

Abstract Realistic simulation of leaf photosynthetic and respiratory processes is needed for accurate prediction of the global carbon cycle. These two processes systematically acclimate to long‐term environmental changes by adjusting photosynthetic and respiratory traits (e.g., the maximum photosynthetic capacity at 25°C ( V cmax,25 ) and the leaf respiration rate at 25°C ( R 25 )) following increasingly well‐understood principles. While some land surface models (LSMs) now account for thermal acclimation, they do so by assigning empirical parameterizations for individual plant functional types (PFTs). Here, we have implemented an Eco‐Evolutionary Optimality (EEO)‐based scheme to represent the universal acclimation of photosynthesis and leaf respiration to multiple environmental effects, and that therefore requires no PFT‐specific parameterizations, in a standard version of the widely used LSM, Noah MP. We evaluated model performance with plant trait data from a 5‐year experiment and extensive global field measurements, and carbon flux measurements from FLUXNET2015. We show that observed R 25 and V cmax,25 vary substantially both temporally and spatially within the same PFT ( C.V. >20%). Our EEO‐based scheme captures 62% of the temporal and 70% of the spatial variations in V cmax,25 (73% and 54% of the variations in R 25 ). The standard scheme underestimates gross primary production by 10% versus 2% for the EEO‐based scheme and generates a larger spread in r (correlation coefficient) across flux sites (0.79 ± 0.16 vs. 0.84 ± 0.1, mean ± S.D.). The standard scheme greatly overestimates canopy respiration (bias: ∼200% vs. 8% for the EEO scheme), resulting in less CO 2 uptake by terrestrial ecosystems. Our approach thus simulates climate‐carbon coupling more realistically, with fewer parameters.

Ren, Yanghang [Department of Earth System Science ↗
DOE OSTI2023

How does humidity data impact land surface modeling of hydrothermal regimes at a permafrost site in Utqia$\mathrm{\dot{g}}$vik, Alaska?

Humidity is a basic and crucial meteorological indicator commonly measured in several forms, including specific humidity, relative humidity, and absolute humidity. These different forms can be inter-derived based on the saturation vapor pressure (SVP). In past decades, dozens of formulae have been developed to calculate the SVP with respect to, and in equilibrium with, liquid water and solid ice surfaces, but many prior studies use a single function for all temperature ranges, without considering the distinction between over the liquid water and ice surfaces. These different approaches can result in humidity estimates that may impact our understanding of surface-subsurface thermal-hydrological dynamics in cold regions. In this study, we compared the relative humidity (RH) downloaded and calculated from four data sources in Alaska based on five commonly used SVP formulas. These RHs, along with other meteorological indicators, were then used to drive physics-rich land surface models at a permafrost-affected site. Here, we found that higher values of RH (up to 40%) were obtained if the SVP was calculated with the over-ice formulation when air temperatures were below freezing, which could lead to a 30% maximum difference in snow depths. The choice of whether to separately calculate the SVP over an ice surface in winter also produced a significant range (up to 0.2 m) in simulated annual maximum thaw depths. The sensitivity of seasonal thaw depth to the formulation of SVP increases with the rainfall rate and the height of above-ground ponded water, while it diminishes with warmer air temperatures. These results show that RH variations based on the calculation of SVP with or without over-ice calculation meaningfully impact physically-based predictions of snow depth, sublimation, soil temperature, and active layer thickness. Under particular conditions, when severe flooding (inundation) and cool air temperatures are present, care should be taken to evaluate how humidity data is estimated for land surface and earth system modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Impacts of Topography-Based Subgrid Scheme and Downscaling of Atmospheric Forcing on Modeling Land Surface Processes in the Conterminous US

The effects of small-scale topography-induced land surface heterogeneity are not well represented in current Earth System Models (ESMs). In this study, a new topography-based subgrid structure referred to as topographic units (TGU) designed to better capture subgrid topographic effects, and methods to downscale atmospheric forcing to the land TGUs have been implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM). Effects of the subgrid scheme and downscaling methods on ELM simulated land surface processes are evaluated over the conterminous United States (CONUS). For this purpose, ELM simulations are performed using two configurations without (NoD ELM) and with (D ELM) downscaling, both using TGUs derived for the 0.5-degree grids and the same land surface parameters. Simulations using the two ELM configurations are compared over the CONUS domain, regional levels, and at observational sites (e.g., SNOTEL). The CONUS-level results suggest that D ELM simulates more snowfall and snow water equivalent (SWE), higher runoff, and less ET during spring and summer. Regional-level results suggest more pronounced impacts of downscaling over regions dominated by higher elevation TGUs and regions with maximum precipitation occurring during cool seasons. Results at the SNOTEL sites suggest that D ELM has superior capability of reproducing the observed SWE at 83% of the sites, with more pronounced performance over topographically heterogeneous TGUs with their maximum precipitation occurring during cool seasons. The results highlight the importance of improving representation of small-scale surface heterogeneity in ESMs and motivate future research to understand their effects on land-atmosphere interactions, streamflow, and water resources management over mountainous regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Updated respiration routines alter spatio-temporal patterns of carbon cycling in a global land surface model

Abstract We updated the routines used to estimate leaf maintenance respiration (MR) in the Energy Land Model (ELM) using a comprehensive global respiration data base. The updated algorithm includes a temperature acclimating base rate, an updated instantaneous temperature response, and new plant functional type specific parameters. The updated MR algorithm resulted in a very large increase in global MR of 16.1 Pg (38%), but the signal was not geographically uniform. The increase was concentrated in the tropics and humid warm-temperate forests. The increase in MR led to large but proportionally smaller decreases in global net primary production (19%) and in average global leaf area index (15%). The effect on global gross primary production (GPP) was a more modest 5.7 Pg (4%). A detailed site level analysis also demonstrated a wide range of effects the updated algorithm can have on the seasonal cycle of GPP. Output from the updated and old models did not differ markedly in how closely they matched a suite of benchmarks. Given the substantial impact on the land surface carbon cycle, a neutral influence on model benchmarks, and better alignment with empirical evidence, an MR algorithm similar to the one presented here should be adopted into ELM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Integrating NDVI-Based Within-Wetland Vegetation Classification in a Land Surface Model Improves Methane Emission Estimations

Earth system models (ESMs) are a common tool for estimating local and global greenhouse gas emissions under current and projected future conditions. Efforts are underway to expand the representation of wetlands in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) by resolving the simultaneous contributions to greenhouse gas fluxes from multiple, different, sub-grid-scale patch-types, representing different eco-hydrological patches within a wetland. However, for this effort to be effective, it should be coupled with the detection and mapping of within-wetland eco-hydrological patches in real-world wetlands, providing models with corresponding information about vegetation cover. In this short communication, we describe the application of a recently developed NDVI-based method for within-wetland vegetation classification on a coastal wetland in Louisiana and the use of the resulting yearly vegetation cover as input for ELM simulations. Processed Harmonized Landsat and Sentinel-2 (HLS) datasets were used to drive the sub-grid composition of simulated wetland vegetation each year, thus tracking the spatial heterogeneity of wetlands at sufficient spatial and temporal resolutions and providing necessary input for improving the estimation of methane emissions from wetlands. Our results show that including NDVI-based classification in an ELM reduced the uncertainty in predicted methane flux by decreasing the model’s RMSE when compared to Eddy Covariance measurements, while a minimal bias was introduced due to the resampling technique involved in processing HLS data. Our study shows promising results in integrating the remote sensing-based classification of within-wetland vegetation cover into earth system models, while improving their performances toward more accurate predictions of important greenhouse gas emissions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Integrating Characteristic Arctic Vegetation in a Land Surface Model Improves Representation of Carbon Dynamics Across a Tundra Landscape

Arctic warming is altering vegetation and carbon dynamics with global implications, yet Earth System Model (ESM) predictions in the Arctic remain highly uncertain, in part due to historically limited data for model parameterization and validation. As such, ESMs typically represent Arctic ecosystems in an oversimplified manner. Recently, nine plant functional types (PFTs) designed to realistically represent tundra vegetation were integrated into the Energy Exascale Earth System Model (E3SM) Land Model (ELM) and parameterized using plot-scale observations from a single site. Additional evaluation was needed to determine their transferability across the Arctic. Here, in this study, we evaluated whether refined representation of tundra vegetation improved model accuracy by conducting spatially explicit 100 × 100 m resolution ELM simulations on Alaska's Seward Peninsula. Simulations with the default two-PFT configuration and with the nine Arctic-specific PFTs were benchmarked against observations of net ecosystem exchange, gross primary production, and aboveground biomass from multiple data streams including an eddy covariance flux tower, flux chambers, and aircraft and unoccupied aerial system hyperspectral remote sensing. Evaluation revealed that Arctic-specific PFT simulations produced more realistic landscape-level carbon exchanges, and better captured observed heterogeneity in biomass and productivity, explaining 60%–70% of spatial variance (R 2 = 0.6–0.7) compared to just 12%–18% (R 2 = 0.12–0.18) with the default configuration. However, the refined model failed to reproduce observed aboveground biomass for highly productive alder-willow communities, requiring further evaluation of carbon allocation parameterizations for tall shrubs that are increasingly expanding across tundra landscapes. Our results demonstrate that enhanced representation of vegetation heterogeneity boosts predictive understanding of tundra carbon dynamics, facilitating regional to pan-Arctic model and remote-sensing scaling.

Murphy, Bailey A. [Oak Ridge National Laboratory (↗
DOE OSTI2022

Improving and testing machine learning methods for benchmarking soil carbon dynamics representation of land surface models

Representation of soil organic carbon (SOC) dynamics in Earth system models (ESMs) is a key source of uncertainty in predicting carbon climate feedbacks. The magnitude of this uncertainty can be reduced by accurate representation of environmental controllers of SOC stocks in ESMs. In this study, we used data of environmental factors, field SOC observations, ESM projections and machine learning approaches to identify dominant environmental controllers of SOC stocks and derive functional relationships between environmental factors and SOC stocks. Our derived functional relationships predicted SOC stocks with similar accuracy as the machine learning approach. We used the derived relationships to benchmark the coupled model intercomparison project phase six ESM representation of SOC stocks. We found divergent environmental control representation in ESMs in comparison to field observations. Representation of SOC in ESMs can be improved by including additional environmental factors and representing their functional relationships with SOC consistent with observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Observational benchmarks inform representation of soil organic carbon dynamics in land surface models

Abstract. Representing soil organic carbon (SOC) dynamics in Earth system models (ESMs) is a key source of uncertainty in predicting carbon–climate feedbacks. Machine learning models can help identify dominant environmental controllers and establish their functional relationships with SOC stocks. The resulting knowledge can be integrated into ESMs to reduce uncertainty and improve predictions of SOC dynamics over space and time. In this study, we used a large number of SOC field observations (n=54 000), geospatial datasets of environmental factors (n=46), and two machine learning approaches (namely random forest, RF, and generalized additive modeling, GAM) to (1) identify dominant environmental controllers of global and biome-specific SOC stocks, (2) derive functional relationships between environmental controllers and SOC stocks, and (3) compare the identified environmental controllers and predictive relationships with those in models used in Phase 6 of the Coupled Model Intercomparison Project (CMIP6). Our results showed that the diurnal temperature, drought index, cation exchange capacity, and precipitation were important observed environmental predictors of global SOC stocks. While the RF model identified 14 environmental factors that describe climatic, vegetation, and edaphic conditions as important predictors of global SOC stocks (R2=0.61, RMSE = 0.46 kg m−2), current ESMs oversimplify the relationships between environmental factors and SOC, with precipitation, temperature, and net primary productivity explaining > 96 % of the variability in ESM-modeled SOC stocks. Further, our study revealed notable disparities among the functional relationships between environmental factors and SOC stocks simulated by ESMs compared with observed relationships. To improve SOC representations in ESMs, it is imperative to incorporate additional environmental controls, such as the cation exchange capacity, and refine the functional relationships to align more closely with observations.

54 ENVIRONMENTAL SCIENCES↗