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.
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.
Modeling land surface processes over a mountainous rainforest in Costa Rica using CLM4.5 and CLM5
This study compares the performance of the Community Land Models (CLM4.5 and CLM5) against tower and ground measurements from a tropical montane rainforest in Costa Rica. The study site receives over 4000 mm of mean annual precipitation and has high daily levels of relative humidity. The measurement tower is equipped with eddy-covariance and vertical profile systems able to measure various micrometeorological variables, particularly in wet and complex terrain. In this work, results from point-scale simulations for both CLM4.5 and its updated version (CLM5) are compared to observed canopy flux and micrometeorological data. Both models failed to capture the effects of frequent rainfall events and mountainous topography on the variables of interest (temperatures, leaf wetness, and fluxes). Overall, CLM5 alleviates some errors in CLM4.5, but CLM5 still cannot precisely simulate a number of canopy processes for this forest. Soil, air, and canopy temperatures, as well as leaf wetness, remain too sensitive to incoming solar radiation rates despite updates to the model. As a result, daytime vapor flux and carbon flux are overestimated, and modeled temperature differences between day and night are higher than those observed. Slope effects appear in the measured average diurnal variations of surface albedo and carbon flux, but CLM5 cannot simulate these features. This study suggests that both CLMs still require further improvements concerning energy partitioning processes, such as leaf wetness process, photosynthesis model, and aerodynamic resistance model for wet and mountainous regions.
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.
Streamlining land surface model Initialization: Automated data retrieval for VELMA using HMS REST API and GDAL
Not Available
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.
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.
Improving the ESA CCI daily soil moisture time series with physically-based land surface model datasets using a Fourier time-filtering method
Explore the source record for details and available documents.
ELM–Wet: Inclusion of a Wet–Landunit With Sub–Grid Representation of Eco–Hydrological Patches and Hydrological Forcing Improves Methane Emission Estimations in the E3SM Land Model (ELM)
Wetlands are the largest emitters of biogenic methane (CH 4 ) and represent the highest source of uncertainty in global CH 4 budgets. Here, we aim to improve the realism of wetland representation in the U.S. Department of Energy's Exascale Earth System Model land surface model, ELM, thereby reducing uncertainty of CH 4 flux predictions. We develop an updated version, ELM-Wet, where we activate a separate landunit for wetlands that handles multiple wetland-specific eco-hydrological patch functional types. We introduce more realistic hydrological forcing through prescribing site-level constraints on surface water elevation, which allows resolving different sustained inundation depth for different patches, and if data exists, prescribing inundation depth. We modified the calculation of aerenchyma transport diffusivity based on observed conductance per leaf area for different vegetation types. We use Bayesian Optimization to parameterize CO 2 and CH 4 fluxes in the developed wet-landunit. Site-level simulations of a coastal non-tidal freshwater wetland in Louisiana were performed with the updated model. Eddy covariance observations of CO 2 and CH 4 fluxes from 2012 to 2013 were used to train the model and data from 2021 were used for validation. Patch-specific chamber flux observations and observations of CH 4 concentration profiles in the soil porewater from 2021 were used for evaluation of the model performance. Our results show that ELM-Wet reduces the model's CH 4 emission root mean squared error by up to 33% and is able to represent inter-daily CO 2 and CH 4 flux variability across the wetland's eco-hydrological patches, including during periods of extreme dry or wet conditions.
Formulation of a Consistent Multi‐Species Canopy Description for Hydrodynamic Models Embedded in Large‐Scale Land‐Surface Representations of Mixed‐Forests
Abstract The plant hydrodynamic approach represents a recent advancement to land surface modeling, in which stomatal conductance responds to water availability in the xylem rather than in the soil. To provide a realistic representation of tree hydrodynamics, hydrodynamic models must resolve processes at the level of a single modeled tree, and then scale the resulting fluxes to the canopy and land surface. While this tree‐to‐canopy scaling is trivial in a homogeneous canopy, mixed‐species canopies require careful representation of the species properties and a scaling approach that results in a realistic description of both the canopy and individual‐tree hydrodynamics, as well as leaf‐level fluxes from the canopy and their forcing. Here, we outline advantages and pitfalls of three commonly used approaches for representing mixed forests in land surface models and present a new framework for scaling vegetation characteristics and fluxes in mixed forests. The new formulation scales fluxes from the tree to canopy level in an energy‐ and mass‐conservative way and allows for a consistent multi‐species/multi‐type canopy description by hydrodynamic models.
WRF-ELM v1.0: a regional climate model to study land–atmosphere interactions over heterogeneous land use regions
Abstract. The Energy Exascale Earth System Model (E3SM) Land Model (ELM) is a state-of-the-art land surface model that simulates the intricate interactions between the terrestrial land surface and other components of the Earth system. Originating from the Community Land Model (CLM) version 4.5, ELM has been under active development, with added new features and functionality, including plant hydraulics, radiation–topography interaction, subsurface multiphase flow, and more explicit land use and management practices. This study integrates ELM v2.1 with the Weather Research and Forecasting (WRF; WRF-ELM) model through a modified Lightweight Infrastructure for Land Atmosphere Coupling (LILAC) framework, enabling affordable high-resolution regional modeling by leveraging ELM's innovative features alongside WRF's diverse atmospheric parameterization options. This framework includes a top-level driver for variable communication between WRF and ELM and Earth System Modeling Framework (ESMF) caps for the WRF atmospheric component and ELM workflow control, encompassing initialization, execution, and finalization. Importantly, this LILAC–ESMF framework demonstrates a more modular approach compared to previous coupling efforts between WRF and land surface models. It maintains the integrity of ELM's source code structure and facilitates the transfer of future developments in ELM to WRF-ELM. To test the ability of the coupled model to capture land–atmosphere interactions over regions with a variety of land uses and land covers, we conducted high-resolution (4 km) WRF-ELM ensemble simulations over the Great Lakes region (GLR) in the summer of 2018 and systematically compared the results against observations, reanalysis data, and WRF-CTSM (WRF coupled with the Community Terrestrial Systems Model). In general, the coupled WRF-ELM model has reasonably captured the spatial distribution of surface state variables and fluxes across the GLR, particularly over the natural vegetation areas. The evaluation results provide a baseline reference for further improvements in ELM in the regional application of high-resolution weather and climate predictions. Our work serves as an example to the model development community for expanding an advanced land surface model's capability to represent fully-coupled land–atmosphere interactions at fine spatial scales. The development and release of WRF-ELM marks a significant advancement for the ELM user community, providing opportunities for fine-scale regional representation, parameter calibration in coupled mode, and examination of new schemes with atmospheric feedback.
Topographic hydro-conditioning to resolve surface depression storage and ponding in a fully distributed hydrologic model
Land surface depressions play a central role in the transformation of rainfall to ponding, infiltration and runoff, yet digital elevation models (DEMs) used by spatially distributed hydrologic models that resolve land surface processes rarely capture land surface depressions at spatial scales relevant to this transformation. Methods to generate DEMs through processing of remote sensing data, such as optical and light detection and ranging (LiDAR) have favored surfaces without depressions to avoid adverse slopes that are problematic for many hydrologic routing methods. Here, in this study, we present a new topographic conditioning workflow, Depression-Preserved DEM Processing (D2P) algorithm, which is designed to preserve physically meaningful surface depressions for depression-integrated and efficient hydrologic modeling. D2P includes several features: (1) an adaptive screening interval for delineation of depressions, (2) the ability to filter out anthropogenic land surface features (e.g., bridges), (3) the ability to blend river smoothing (e.g., a general downslope profile) and depression resolving functionality. From a case study in the Goodwin Creek Experimental Watershed, D2P successfully resolved 86% of the ponds at a DEM resolution of 10 m. Topographic conditioning was achieved with minimum impact as D2P reduced the number of modified cells from the original DEM by 51% compared to a conventional algorithm. Furthermore, hydrologic simulation using a D2P processed DEM resulted in a more robust characterization on surface water dynamics based on higher surface water storage as well as an attenuated and delayed peak streamflow.
More Realistic Intermediate Depth Dry Firn Densification in the Energy Exascale Earth System Model (E3SM)
Abstract Earth system models account for seasonal snow cover, but many do not accommodate the deeper snowpack on ice sheets (aka firn) that slowly transforms to ice under accumulating snowfall. To accommodate and resolve firn depths of up to 60 m in the Energy Exascale Earth System Model's land surface model (ELM), we add 11 layers to its snowpack and evaluate three dry snow compaction equations in multi‐century simulations. After comparing results from ELM simulations (forced with atmospheric reanalysis) with empirical data, we find that implementing into ELM a two‐stage firn densification model produces more accurate dry firn densities at intermediate depths of 20–60 m. Compared to modeling firn using the equations in the (12 layer) Community Land Model (version 5), switching to the two‐stage firn densification model (with 16 layers) significantly decreases root‐mean‐square errors in upper 60 m dry firn densities by an average of 41 kg m −3 (31%). Simulations with three different firn density parameterizations show that the two‐stage firn densification model should be used for applications that prioritize accurate upper 60 m firn air content (FAC) in regions where the mean annual surface temperature is greater than roughly −31°C. Because snow metamorphism, firn density, and FAC are major components in modeling ice sheet surface albedo, melt water retention, and climatic mass balance, these developments advance broader efforts to simulate the response of land ice to atmospheric forcing in Earth system models.
Gridded Sub-daily Climate Forcings for North America Based on Daymet and ERA5 (Daymet-ERA5)
To support high spatial and temporal resolution land surface modeling, this dataset provides hourly time step historic weather forcing at 1-km spatial resolution for the entire North America. The latest Daymet V4 data provides gridded historic daily weather observations at 1-km spatial resolution from 1980 to 2024. Using sub-daily temporal information from the ECMWF ERA5 reanalysis, Daymet was further temporally downscaled to hourly time steps and provided in the format required for land surface model simulations. The process of temporal downscaling preserves the relative magnitude in each hourly time step from ERA5 while maintaining the total and average values from Daymet for each day. This results in a blended 1980-2024 Daymet-ERA5 dataset. Available variables include surface air temperature, precipitation, specific humidity, shortwave and longwave radiation, wind speed, and pressure. These data can be used as a high-resolution meteorological forcing dataset to support high-resolution land surface modeling where accurate meteorological forcing datasets built from historic observations and/or reanalysis datasets are desirable.
Improved Representations of Land‐Atmosphere Interactions Over the Continental U.S. Through Dynamic Root Modeling
Previous studies have identified the oversimplified root system representation as a key factor leading to inaccuracies in vegetation‐atmosphere feedbacks. In this study, a dynamic root water uptake scheme in the Noah‐MP land surface model has been coupled to the Weather Research and Forecasting (WRF) model to investigate its impact on the surface hydroclimate variables and land‐atmosphere interactions. To evaluate the impact of the dynamic root, two coupled simulations were conducted, one with the dynamic root water uptake scheme (DynRt) and one with the static root water uptake scheme (StcRt), which is based on the default root representation in Noah‐MP, with slight modifications, primarily in vegetation‐related parameters. Both DynRt and StcRt simulations were conducted with a small ensemble of three members to account for variations in physical parameterizations, initial and boundary forcing and model setup. When compared with reference data sets, the DynRt simulations show improved results than the StcRt simulations, reducing biases in the simulated leaf area index, surface energy fluxes, soil moisture and precipitation. Two different mechanisms through which roots affect land‐atmosphere coupling have been identified. Over the transitional climate zone between the dry and wet climate, the dynamic root scheme affects surface climate and land‐atmosphere coupling mainly through changes in soil moisture through hydraulic redistribution by plant root system. Over the energy‐limited mesic zone, the dynamic root affects regional land‐atmosphere coupling mainly through changes in carbon allocation. This work highlights the importance of dynamic root representation in improving vegetation‐atmosphere simulations by enhancing predictions of water, energy, and carbon fluxes.
Maps of growing season gross primary production and net ecosystem exchange for Council Road Mile Marker 71, Seward Peninsula, Alaska, [2017-2023]
This data archive is in support of the Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) publication "Integrating Characteristic Arctic Vegetation in a Land Surface Model Improves Representation of Carbon Dynamics Across a Tundra Landscape", by Murphy et al. (2025a). Murphy et al. (2025a) evaluated whether incorporating observed Arctic vegetation heterogeneity into ELM, the land model of the Department of Energy’s Energy Exascale Earth System Model (E3SM), improved simulations of tundra carbon cycling. The associated model archive can be found at Murphy et al. (2025b). The study focused on the spatial patterns and net landscape-level growing season productivity and carbon uptake. As part of this evaluation, observationally derived maps of average growing season (June–August) net ecosystem exchange (NEE) and gross primary production (GPP) were developed for the same domain. These maps, which form the dataset described here, integrate eddy covariance flux tower, remote sensing, and vegetation community data to provide spatially explicit benchmarks for model evaluation. The maps provide spatially explicit estimates of average growing season NEE and GPP across 13 tundra vegetation communities within the study domain. By combining flux tower observations with Airborne Visible-Infrared Imaging Spectrometer-Next Generation (AVIRIS-NG) hyperspectral imagery and drone-based normalized difference vegetation index (NDVI), these maps capture the heterogeneity of carbon fluxes associated with different Arctic vegetation types. While they represent average seasonal conditions rather than interannual variability, the maps provide a unique dataset for evaluating model performance, comparing vegetation community contributions to landscape-scale carbon cycling, and supporting regional analyses of Arctic carbon dynamics. This data archive contains 5 m resolution maps of vegetation communities, vegetation community average growing season GPP, and vegetation community average growing season NEE (three *.tif files), a User’s Guide (*pdf file), and Table 1 of the User’s Guide displaying vegetation community coverage and average growing season NEE and GPP values (*.csv file).
Large Ensemble Diagnostic Evaluation of Hydrologic Parameter Uncertainty in the Community Land Model Version 5 (CLM5)
Abstract Land surface models such as the Community Land Model version 5 (CLM5) seek to enhance understanding of terrestrial hydrology and aid in the evaluation of anthropogenic and climate change impacts. However, the effects of parametric uncertainty on CLM5 hydrologic predictions across regions, timescales, and flow regimes have yet to be explored in detail. The common use of the default hydrologic model parameters in CLM5 risks generating streamflow predictions that may lead to incorrect inferences for important dynamics and/or extremes. In this study, we benchmark CLM5 streamflow predictions relative to the commonly employed default hydrologic parameters for 464 headwater basins over the conterminous United States (CONUS). We evaluate baseline CLM5 default parameter performance relative to a large (1,307) Latin Hypercube Sampling‐based diagnostic comparison of streamflow prediction skill using over 20 error measures. We provide a global sensitivity analysis that clarifies the significant spatial variations in parametric controls for CLM5 streamflow predictions across regions, temporal scales, and error metrics of interest. The baseline CLM5 shows relatively moderate to poor streamflow prediction skill in several CONUS regions, especially the arid Southwest and Central U.S. Hydrologic parameter uncertainty strongly affects CLM5 streamflow predictions, but its impacts vary in complex ways across U.S. regions, timescales, and flow regimes. Overall, CLM5's surface runoff and soil water parameters have the largest effects on simulated high flows, while canopy water and evaporation parameters have the most significant effects on the water balance.
Influence of sun zenith angle on canopy clumping and the resulting impacts on photosynthesis
Addressing the impact of vegetation architecture on shortwave radiation transfer in land surface models is important for accurate weather forecasting, carbon budget estimates, and climate predictions. This paper investigates to what extent it is possible to retrieve structural parameters of two different parameterization schemes from direct transmittance derived from digital hemispherical photography and 3D radiative transfer modeling for two study sites with different vegetation canopy architectures. Neglecting the representation of 3D canopy structure in radiative transfer schemes leads to significant errors in shortwave radiation partitioning (up to 3.5 times more direct transmittance in the 3D model). Structural parameters, referred to as whole canopy ‘clumping indices’, were obtained in order to evaluate the impact of angular variation in clumping on shortwave radiation transfer. Impacts on photosynthesis were evaluated at site level with the UKESM land surface model, JULES. A comparison between flux tower derived and modeled photosynthesis indicates that considering zenith angular variations of structural parameters in the radiative transfer scheme of the UKESM land surface model significantly improves photosynthesis prediction in light limited ecosystems (from RMSE = 2.91 μmol CO 2 .m -2 .s -1 to RMSE = 1.51 μmol CO 2 .m -2 .s -1 , 48% smaller), typically with enhanced photosynthesis from bottom layers.