Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic Vegetation”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Microscale spatiotemporal patterns of water, soil organic carbon, and enzymes in plant litter detritusphere

Understanding biophysical and biochemical processes in detritusphere is critical to quantifying and modeling plant residue decomposition dynamics and subsequent soil organic carbon (C) accrual. The objectives of the study were to explore (i) micro-environmental conditions within the detritusphere formed around decomposing corn and soybean leaves, and (ii) the relationships between moisture distribution, enzyme activity, and C dynamics during decomposition. We assessed spatial and temporal dynamics of moisture distribution using X-ray and neutron computed tomography, and activities of β-glucosidase and chitinase, two enzymes involved in soil C and N processing, using zymography. We used 13 C labeled residue to track residue contribution to atmospheric CO 2 and soil organic C. Moisture redistribution pattern varied depending on the residue type. While the water was immediately absorbed by the corn leaves and maintained afterward, switchgrass leaves absorbed water more slowly and created water-deficient zones within ~1mm from the residue. This initial moisture depletion led to lower chitinase activity and residue-derived CO 2 emissions in switchgrass. In contrast to chitinase, β-glucosidase activity was influenced by a combination of vegetation history and residue type, and it was higher when the origin of the residue matched the vegetation history of the soil. Pore size had an opposite impact on the studied enzymes, supporting the notion that contrasting soil pore architecture can stimulate activities of different enzymes through a selection of dominant enzyme producers. We concluded that the decomposition dynamics of plant residues is not only a simple function of residue chemistry, but rather a combined effect of the vegetation history, in part through its effect on microbial community composition, the plant residue chemical and likely physical characteristics, and the soil pore structure in the detritusphere. Together, they create temporally dynamic micro-environmental conditions influencing decomposition. Specifically, our study demonstrated that the initial micro-environment formulated in detritusphere can play an important role in enzyme activities and consequent C dynamics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A machine learning approach targeting parameter estimation for plant functional type coexistence modeling using ELM-FATES (v2.0)

Tropical forest dynamics play a crucial role in the global carbon, water, and energy cycles. However, realistically simulating the dynamics of competition and coexistence between different plant functional types (PFTs) in tropical forests remains a significant challenge. This study aims to improve the modeling of PFT coexistence in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), a vegetation demography model implemented in the Energy Exascale Earth System Model (E3SM) land model (ELM), ELM-FATES. Specifically, we explore (1) whether plant trait relationships established from field measurements can constrain ELM-FATES simulations and (2) whether machine learning (ML)-based surrogate models can emulate the complex ELM-FATES model and optimize parameter selections to improve PFT coexistence modeling. We conducted three ensembles of ELM-FATES experiments at a tropical forest site near Manaus, Brazil. By comparing the ensemble experiments without (Exp-CTR) and with (Exp-OBS) consideration of observed trait relationships, we found that accounting for these relationships slightly improves the simulations of water, energy, and carbon variables when compared to observations but degrades the simulation of PFT coexistence. Using ML-based surrogate models trained on Exp-CTR, we optimized the trait parameters in ELM-FATES and conducted another ensemble of experiments (Exp-ML) with these optimized parameters. The proportion of PFT coexistence experiments significantly increased from 21 % in Exp-CTR to 73 % in Exp-ML. After filtering the experiments that allow for PFT coexistence to agree with observations (within 15 % tolerance), 33 % of the Exp-ML experiments were retained, which is a significant improvement compared to the 1.4 % in Exp-CTR. Exp-ML also accurately reproduces the annual means and seasonal variations in water, energy, and carbon fluxes and the field inventory of aboveground biomass. This study represents a reproducible method that utilizes machine learning to identify parameter values that improve model fidelity against observations and PFT coexistence in vegetation demography models for diverse ecosystems. Our study also suggests the need for new mechanisms to enhance the robust simulation of coexisting plants in ELM-FATES and has significant implications for modeling the response and feedbacks of ecosystem dynamics to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Dynamic ecosystem assembly and escaping the “fire trap” in the tropics: insights from FATES_15.0.0

Abstract. Fire is a fundamental part of the Earth system, with impacts on vegetation structure, biomass, and community composition, the latter mediated in part via key fire-tolerance traits, such as bark thickness. Due to anthropogenic climate change and land use pressure, fire regimes are changing across the world, and fire risk has already increased across much of the tropics. Projecting the impacts of these changes at global scales requires that we capture the selective force of fire on vegetation distribution through vegetation functional traits and size structure. We have adapted the fire behavior and effects module, SPITFIRE (SPread and InTensity of FIRE), for use with the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), a size-structured vegetation demographic model. We test how climate, fire regime, and fire-tolerance plant traits interact to determine the biogeography of tropical forests and grasslands. We assign different fire-tolerance strategies based on crown, leaf, and bark characteristics, which are key observed fire-tolerance traits across woody plants. For these simulations, three types of vegetation compete for resources: a fire-vulnerable tree with thin bark, a vulnerable deep crown, and fire-intolerant foliage; a fire-tolerant tree with thick bark, a thin crown, and fire-tolerant foliage; and a fire-promoting C4 grass. We explore the model sensitivity to a critical parameter governing fuel moisture and show that drier fuels promote increased burning, an expansion of area for grass and fire-tolerant trees, and a reduction of area for fire-vulnerable trees. This conversion to lower biomass or grass areas with increased fuel drying results in increased fire-burned area and its effects, which could feed back to local climate variables. Simulated size-based fire mortality for trees less than 20 cm in diameter and those with fire-vulnerable traits is higher than that for larger and/or fire-tolerant trees, in agreement with observations. Fire-disturbed forests demonstrate reasonable productivity and capture observed patterns of aboveground biomass in areas dominated by natural vegetation for the recent historical period but have a large bias in less disturbed areas. Though the model predicts a greater extent of burned fraction than observed in areas with grass dominance, the resulting biogeography of fire-tolerant, thick-bark trees and fire-vulnerable, thin-bark trees corresponds to observations across the tropics. In areas with more than 2500 mm of precipitation, simulated fire frequency and burned area are low, with fire intensities below 150 kW m−1, consistent with observed understory fire behavior across the Amazon. Areas drier than this demonstrate fire intensities consistent with those measured in savannas and grasslands, with high values up to 4000 kW m−1. The results support a positive grass–fire feedback across the region and suggest that forests which have existed without frequent burning may be vulnerable at higher fire intensities, which is of greater concern under intensifying climate and land use pressures. The ability of FATES to capture the connection between fire disturbance and plant fire-tolerance strategies in determining biogeography provides a useful tool for assessing the vulnerability and resilience of these critical carbon storage areas under changing conditions across the tropics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Soil water percolation and nutrient fluxes as a function of topographical, seasonal and soil texture variation in Central Amazonia, Brazil

Abstract Understanding soil water dynamics and transport of nutrients is challenging in tropical rainforests due to the uniqueness of several properties related to soils, vegetation and seasonality that make relating patterns found in temperate environments to tropical sites difficult. We address the need for better edaphic characterization in tropical environments by investigating soil water percolation rates and chemistry across topographic, soil texture and seasonal gradients in a mature tropical rainforest in Central Amazonia, Brazil. We utilized a passive wick flux meter (e.g., drainage lysimeter) to directly measure real‐time percolation fluxes at 60‐cm depth, and to sample a suite of chemical species across plateau, slope and valley topographic positions. We found percolation flux volume and chemical exports generally increase with decreasing elevation and clay content, which was lowest in the valley. Daily percolation flux was observed to be 2.39 ± 0.44 in plateau, 3.01 ± 0.50 in slope and 6.16 ± 0.83 mm in valley. Most solutes were present in small amounts of <1 mg L −1 , such as PO₄ 3− , Fe 2+ /Fe 3+ and Mn 2+ ; however, NO 3 − concentrations were >20 mg L −1 , even exceeding 100 mg L −1 in the valley. Based on additional isotopic analysis, we speculate high NO 3 − concentrations are partially an artefact of root decomposition following installation of the flux meters. The empirical relationships we show among percolation volume and nutrient exports under varying topographies and soil textures can improve Earth System Model performance by better constraining ecohydrological relationships to nutrient fluxes, which can in‐turn better illuminate the important factors that govern their behaviour.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Monitoring agroecosystem productivity and phenology at a national scale: A metric assessment framework

Effective measurement of seasonal variations in the timing and amount of production is critical to managing spatially heterogeneous agroecosystems in a changing climate. Although numerous technologies for such measurements are available, their relationships to one another at a continental extent are unknown. Using data collected from across the Long-Term Agroecosystem Research (LTAR) network and other networks, we investigated correlations among key metrics representing primary production, phenology, and carbon fluxes in croplands, grazing lands, and crop-grazing integrated systems across the continental U.S. Metrics we examined included gross primary productivity (GPP) estimated from eddy covariance (EC) towers and modelled from the Landsat satellite, Landsat NDVI, and vegetation greenness (Green Chromatic Coordinate, GCC) from tower-mounted PhenoCams for 2017 and 2018. Overall, our analysis compared production dynamics estimated from three independent ground and remote platforms using data for 34 agricultural sites constituting 51 site-years of co-located time series. Pairwise sensor comparisons across all four metrics revealed stronger correlation and lower root mean square error (RMSE) between end of season (EOS) dates (Pearson R ranged from 0.6 to 0.7 and RMSE from 32.5 to 67.8) than start of season (SOS) dates (0.46 to 0.69 and 40.4 to 66.2). Overall, moderate to high correlations between SOS and EOS metrics complemented one another except at some lower productivity grazing land sites where estimating SOS can be challenging. Growing season length estimates derived from 16-day satellite GPP (179.1 days) were significantly longer than those from PhenoCam G CC (70.4 days, p adj < 0.0001) and EC GPP (79.6 days, p adj < 0.0001). Landscape heterogeneity did not explain differences in SOS and EOS estimates. Annual integrated estimates of productivity from EC GPP and PhenoCam G CC diverged from those estimated by Landsat GPP and NDVI at sites where annual production exceeds 1000 gC/m –2 yr –1 . Based on our results, we developed a “metric assessment framework” that articulates where and how metrics from satellite, eddy covariance and PhenoCams complement, diverge from, or are redundant with one another. The framework was designed to optimize instrumentation selection for monitoring, modeling, and forecasting ecosystem functioning with the ultimate goal of informing decision-making by land managers, policy-makers, and industry leaders working at multiple scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Modeling Spatial Heterogeneity in Surface Turbulent Heat Flux in the U.S. Southern Great Plains

Advances in numerical modeling of cloud dynamics are driving a need for improved land model prediction at convective storm scales. Here, satellite and ground-based vegetation remote sensing data were combined with land model experiments to more accurately characterize land surface spatial heterogeneity in the Community Land Model (CLM4.0). The new subgrid classification of plant functional types (PFT) and leaf area index (LAI) enables consistent comparison between models and ground-based flux measurements in the U.S. southern Great Plains. Errors in vegetation data sets (inferred from comparison between 250 m satellite and ground-based LAI), while large, had less impact on the simulated characteristics of spatial heterogeneity compared to errors in model representation of surface energy partitioning (between latent and sensible heat flux) and its relationship to LAI. Predicted spatial heterogeneity in surface energy partitioning was enhanced after replacing soil and stomatal resistance parameters with a new set that better predicts the observed relationship to LAI. These modifications increase the number of smaller (mesoscale) dry land patches having higher sensible heat flux. The parameter experiments suggest that vegetation state and processes (transpiration) act to broaden the size spectrum of surface heat flux heterogeneity, which can influence clouds and convective initiation. Moreover, improvements in vegetation input data and model parameters had partially compensating effects on surface flux heterogeneity, indicating the importance of evaluating input data and parameterizations together to improve prediction at higher spatial resolutions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Floods due to Atmospheric Rivers along the U.S. West Coast: The Role of Antecedent Soil Moisture in a Warming Climate

Precipitation extremes are projected to become more frequent along the U.S. West Coast due to increased atmospheric river (AR) activity, but the frequency of less intense precipitation events may decrease. Antecedent soil moisture (ASM) conditions can have a large impact on flood responses, especially if prestorm precipitation decreases. Taken together with increased antecedent evaporative demand due to warming, this would result in reduced soil moisture at the onset of extreme precipitation events. In this work, we examine the impact of ASM on AR-related floods in a warming climate in three basins that form a transect along the U.S. Pacific Coast: the Chehalis River basin in Washington, the Russian River basin in Northern California, and the Santa Margarita River basin in Southern California. Additionally, we ran the Distributed Hydrology Soil Vegetation Model (DHSVM) over the three river basins using forcings downscaled from 10 global climate models (GCMs). We examined the dynamic role of ASM by comparing the changes in the largest 50, 100, and 150 extreme events in two periods, 1951–2000 and 2050–99. In the Chehalis basin, the projected fraction of AR-related extreme discharge events slightly decreases. In the Russian basin, this fraction increases, however, and more substantially so in the Santa Margarita basin. This is due to increases in AR-related extreme precipitation events, as well as the fact that the relationship of extreme precipitation to extreme discharge is strengthened by projected increases in year-to-year volatility of annual precipitation in California, which increases the likelihood of concurrent occurrence of large storms and wet ASM conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Data used in “Enguehard et al. 2022, Machine-Learning Functional Zonation Approach for Characterizing Terrestrial–Aquatic Interfaces: Application to Lake Erie”

The package contains the data layers used in “Enguehard et al. 2022, Machine-Learning Functional Zonation Approach for Characterizing Terrestrial–Aquatic Interfaces: Application to Lake Erie”. Spatial data layers include: topography, wetland vegetation cover, time series of Landsat’s enhanced vegetation index (EVI) between 1990 and 2020. The study aims to characterize coastal wetlands with particular focus on the co-variability between plant dynamics, topography, soil, and other environmental factors. We proposed a functional zonation approach based on machine learning clustering to identify the spatial regions, i.e., zones that capture these co-varied properties. This approach was applied to publicly available datasets along Lake Erie, in the Great Lakes Region

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Effects of spatial variability in vegetation phenology, climate, landcover, biodiversity, topography, and soil property on soil respiration across a coastal ecosystem

Coastal terrestrial-aquatic interfaces (TAIs) are crucial contributors to global biogeochemical cycles and carbon exchange. A systematic evaluation of the interaction between coastal catchment properties and carbon dioxide (CO2) emission by soil respiration is significant for assessing carbon dynamics and predicting the future trajectory of atmospheric CO2 concentrations in coastal TAIs. The soil CO2 efflux in these transition zones is however poorly understood due to the high spatiotemporal dynamics of TAIs, as various sub-ecosystems in this region are compressed and expanded by complex influences of tides, changes in river levels, climate, and land use. We focus on the Chesapeake Bay region to (i) investigate the spatial heterogeneity of the coastal ecosystem and identify spatial zones with similar environmental characteristics based on the spatial data layers, including vegetation index (kNDVI), climate, landcover, diversity, topography, soil property, and relative tidal elevation; (ii) understand the primary driving factors affecting soil respiration within sub-ecosystems of the coastal ecosystem. Specifically, we employed hierarchical clustering analysis to identify spatial regions with distinct environmental characteristics, followed by the determination of main driving factors using Random Forest regression and SHapley Additive exPlanations. Maximum and minimum temperature are the main drivers common to all sub-ecosystems, while each region also has additional unique major drivers that differentiate them from one another. Precipitation exerts an influence on vegetated lands, while soil pH value holds importance specifically in forested lands. In croplands characterized by high clay content and low sand content, the significant role is attributed to bulk density. Wetlands demonstrate the importance of both elevation and sand content, with clay content being more relevant in non-inundated wetlands than in inundated wetlands. The topographic wetness index significantly contributes to the mixed vegetation areas, including shrub, grass, pasture, and forest. Additionally, our research reveals that dense vegetation land covers and urban/developed areas exhibit distinct soil property drivers. Overall, there is no one-size-fits-all approach to modeling carbon fluxes in coastal TAIs, and our study highlights the importance of further research and monitoring practices to improve our understanding of carbon dynamics and promote the sustainable management of coastal TAIs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A new framework to map fine resolution cropping intensity across the globe: Algorithm, validation, and implication

We report accurate estimation of cropping intensity (CI), an indicator of food production, is well aligned with the ongoing efforts to achieve sustainable development goals (SDGs) under diminishing natural resources. The advancement in satellite remote sensing provides unprecedented opportunities for capturing CI information in a spatially continuous manner. However, challenges remain due to the lack of generalizable algorithms for accurately and efficiently mapping global CI with a fine spatial resolution. In this study, we developed a 30-m planetary-scale CI mapping framework with the reconstructed time series of Normalized Difference Vegetation Index (NDVI) from multiple satellite images. Using a binary crop phenophase profile indicating growing and non-growing periods, we estimated pixel-by-pixel CI by enumerating the total number of valid cropping cycles during the study years. Based on the Google Earth Engine cloud computing platform, we implemented the framework to estimate CI during 2016–2018 in eight geographic regions across continents that are representative of global cropping system diversity. Comparison with PhenoCam network data in four cropland sites suggests that the proposed framework is capable of capturing the seasonal dynamics of cropping practices. Spatially, overall accuracies based on validation samples range from 80.0% to 98.9% across different regions worldwide. Regarding the CI classes, single cropping systems are associated with more robust and less biased estimations than multiple cropping systems. Finally, our CI estimates reveal high agreement with two widely used land surface phenology products, including Vegetation Index and Phenology V004 (VIP4) and Moderate Resolution Imaging Spectroradiometer Land Cover Dynamics (MCD12Q2), meanwhile providing much more spatial details. Due to its robustness, the developed CI framework can be potentially generalized to produce global fine resolution CI products for food security and other applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Boreal–Arctic wetland methane emissions modulated by warming and vegetation activity

Abstract Wetland methane (CH 4 ) emissions over the Boreal–Arctic region are vulnerable to climate change and linked to climate feedbacks, yet understanding of their long-term dynamics remains uncertain. Here, we upscaled and analysed two decades (2002–2021) of Boreal–Arctic wetland CH 4 emissions, representing an unprecedented compilation of eddy covariance and chamber observations. We found a robust increasing trend of CH 4 emissions (+8.9%) with strong inter-annual variability. The majority of emission increases occurred in early summer (June and July) and were mainly driven by warming (52.3%) and ecosystem productivity (40.7%). Moreover, a 2 °C temperature anomaly in 2016 led to the highest recorded annual CH 4 emissions (22.3 Tg CH 4 yr −1 ) over this region, driven primarily by high emissions over Western Siberian lowlands. However, current-generation models from the Global Carbon Project failed to capture the emission magnitude and trend, and may bias the estimates in future wetland CH 4 emission driven by amplified Boreal–Arctic warming and greening.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

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

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

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

High-Resolution Spatio-Temporal Estimation of Net Ecosystem Exchange in Ice-Wedge Polygon Tundra Using In Situ Sensors and Remote Sensing Data

Land-atmosphere carbon exchange is known to be extremely heterogeneous in arctic ice-wedge polygonal tundra regions. In this study, a Kalman filter-based method was developed to estimate the spatio-temporal dynamics of daytime average net ecosystem exchange (NEEday) at 0.5-m resolution over a 550 m by 700 m study site. We integrated multi-scale, multi-type datasets, including normalized difference vegetation indices (NDVIs) obtained from a novel automated mobile sensor system (or tram system) and a greenness index map obtained from airborne imagery. We took advantage of the significant correlations between NDVI and NEEday identified based on flux chamber measurements. The weighted average of the estimated NEEday within the flux-tower footprint agreed with the flux tower data in term of its seasonal dynamics. We then evaluated the spatial variability of the growing season average NEEday, as a function of polygon geomorphic classes; i.e., the combination of polygon types—which are known to present different degradation stages associated with permafrost thaw—and microtopographic features (i.e., troughs, centers and rims). Our study suggests the importance of considering microtopographic features and their spatial coverage in computing spatially aggregated carbon exchange.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Model simulations of Plum Island Ecosystems LTER low marsh site using ELM-PFLOTRAN

Model simulations using the E3SM Land Model (ELM) coupled to the PFLOTRAN reactive transport model via the Alquimia interface. The simulations were conducted for a tidal salt marsh at the Plum Island Ecosystems LTER near Rowley, Massachusetts, USA. Model simulations were forced using site-specific tidal cycles and salinity, and the simulations used a biogeochemical reaction network including aerobic decomposition, sulfate reduction, iron reduction, and methanogenesis. Model outputs include simulated carbon stocks, carbon dioxide and methane fluxes, and porewater concentrations of key solutes related to sulfur, iron, and carbon cycling. The model simulations included a saline simulation (with tidal sulfate inputs), a fresh simulation (with low salinity and low sulfate inputs), and a saline simulation with lower vegetation productivity to represent the effect of salinity on vegetation. These simulations were conducted to demonstrate that a new model framework incorporating subsurface redox and biogeochemical interactions into a land surface model could reproduce measured surface greenhouse gas fluxes and biogeochemical dynamics in tidal marsh ecosystems, and to test whether including redox interactions in a land surface model would allow the model to resolve contrasts in biogeochemical cycling and greenhouse gas production between saline and freshwater wetlands.The data package includes gzipped tar archives (which can be expanded using standard tar and gzip utilities) of model outputs from three model configurations: saline subsurface and reduced vegetation productivity related to salinity; saline subsurface with vegetation productivity not reduced; and freshwater. Also included are code for the modified E3SM model, Alquimia interface, and PFLOTRAN reactive transport simulator in gzipped tar format; plain text parameter and configuration files; python code files for visualizing model output and defining model configurations; and model output, tide and salinity forcing, and configuration files in netCDF format. See the README.md file in the data package for a detailed description of all files contained in the package. All files are in netCDF (.nc), gzipped tar archive (.tar.gz or .tgz), or text (all other files).Updated: May 13, 2024. Model output, E3SM code, PFLOTRAN input files, and python codes for visualizing results were updated to reflect changes made for the manuscript revision. The updated archive reflects the code and model output from the final accepted manuscript. Changes included updated reaction parameters reflecting improved parameterization and additional comparisons with field measurements. E3SM code changes included better support for multiple grid cells and improved flow and transport parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Water foraging with dynamic roots in E3SM; The role of roots in terrestrial ecosystem memory on intermediate timescales (Final Technical Report)

Terrestrial ecosystems can show sustained responses to stress events that may last for years after the initial perturbation. This phenomenon is referred to as legacy or memory and it emerges from a set of ecosystem processes that are poorly captured by Earth System and Land Surface Models. Consequently, these models struggle to predict the impact that extreme climate events have on surface energy, carbon and hydrological exchange once a stressor such as drought has relaxed or in response to repeated exposure to stress. In this project, we set out to test how the addition of dynamic root profiles in land surface models impact the capacity for Earth System Models, such as the Department of Energy’s E3SM, to capture realistic legacy effects. Observations have shown that vegetation shifts their root profiles during stress events to forage for water and these altered root profiles may take years to relax back to the initial state. We hypothesized that the alteration of the belowground root structures may be a source of terrestrial ecosystem legacy that is missing from models. To test this idea we undertook three core activities. (1) We developed a global-scale analysis of the impact of dynamic roots on ecosystem legacy building from a recently developed dynamic root module. (2) We developed new root dynamics in the DOE’s Energy Land Model (ELM) that allowed not only root profiles to shift their profile but also simultaneously alter carbon allocation to fine root pools. We then tested the impact of these new dynamics with intensive sensitivity analysis at four long-term AmeriFlux sites. (3) We undertook detailed isotopic analysis of tree rings from these AmeriFlux sites to assess ecophysiological and ecohydrological legacy to stress events that can be used to benchmark the sensitivity experiments with ELM. From these core activities, we report the following key findings. Firstly, on a global scale, the addition of dynamic roots led to chronically water-stressed ecosystems recovering faster to climate stress while wetter ecosystem showed enhanced legacy. This is because across ecosystems, stress events were almost universally associated with water shortages that led to the development of deeper root profiles. These deeper root profiles proved beneficial for recovery from drought stress. While the root dynamics did not universally improve the modeled representation of legacy it showed complex transient dynamics that emerge from the addition of dynamic roots. Secondly, the sensitivity analysis illustrated long term shifts in rooting depths away from the prescribed default profile suggesting that initializing of root profiles could benefit from spin-up simulations that converge on locally optimized root profiles. In addition, by enabling dynamic allocation some of the sites predicted unrealistically low allocation to roots (and high allocation to leaves). While these changes did not dramatically alter modeled gross primary productivity, they illustrated that without more sophisticated root processes in models, there is little penalty to dramatically disinvest in roots. Thirdly, the isotopic analysis of tree rings showed highly distinct legacy responses across species and sites. For example, T. canadensis showed reduced transpiration the year after stress events illustrated by sustained elevated $\delta^{18}$O. In contrast, A saccharum displayed elevated $\delta^{13}$C associated with reduced stomatal conductance in response to the previous years’ stress event. These geochemical signatures of legacy were present despite tree growth returning to normal the year after the stress. The results show the importance of species-level dynamics in legacy that are absent in modeling that assumes common traits within plant functional types. In summary, the work here established new avenues to explore root dynamics in models while illustrating how additional root processes are needed before dynamic carbon allocation can be implemented. Lastly, this project provided mentorship to a postdoctoral fellow, training for an early career scientist and multiple undergraduate students recruited from a minority serving institution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Pinus ponderosa physiology data from the 2023 girdling manipulation field experiment at Game Ridge, Missoula County, Montana, USA

This dataset contains physiology data on Pinus ponderosa saplings, including: leaf water potential (MPa, megapascals), leaf osmotic potential (MPa, megapascals), leaf pressure potential (MPa, megapascals), leaf relative water content (% saturated mass), as well as soluble sugars (% dry mass),starch (% dry mass), and total nonstructural carbohydrates (% dry mass) found in branch phloem and xylem. Data files and data dictionary(ies) are uploaded as .csv files and .xlsx files. The Users Guide is a .pdf file. Location data can be found in the Google Earth file GameRidge_SitePlotCoordinates.kmz.kml included here. These datasets were collected for Plant Carbohydrate Depletion, Mycorrhizal Networks, and Vulnerability to Drought: An Experimental Test in the Field. This experiment examined the interdependency between plant hydraulics and carbohydrate availability and sought to develop ways to incorporate interactions with below ground symbiotic organisms to better model and quantify forest response to drought.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Plant‐mediated methane transport in emergent and floating‐leaved species of a temperate freshwater mineral‐soil wetland

Abstract Methane flux from freshwater mineral‐soil (FWMS) wetlands and its variability among sites is largely modulated by plant‐mediated transport. However, plant‐mediated transport processes are rarely resolved in land surface models and are poorly parametrized for plants commonly found in FWMS wetlands. Here, relationships between methane flux and CO 2 uptake, as well as plant conductance of methane were evaluated for three plant species and two characteristic functional types: emergent (narrow‐leaved cattail) and floating‐leaved (American lotus and water lily). We found significant but contrasting correlations between methane flux and CO 2 uptake in cattails ( r 2 = 0.51, slope = −0.16, during morning) and water lily ( r 2 = 0.32, slope = 0.064, after midday). This relationship was not significant in American lotus, showing that stomata regulation of methane fluxes is species‐specific and not generalizable across the floating‐leaved plant functional type. Conductance of methane per leaf area showed distinct seasonal dynamics across species. Conductance was similar among the floating‐leaved species (6.2 × 10 −3 m d −1 in lotus and 7.2 × 10 −3 m d −1 in water lily) and higher than conductance in the emergent species (2.7 × 10 −3 m d −1 ). Our results provide direct observations of plant conductance rates and identify the vegetation parameters (leaf area, stomatal conductance) that modify them. Our results further suggest that models of methane emissions from FWMS should parameterize plant‐mediated transport in different plant functional types, scaled by leaf area and with variable seasonal phenological dynamics, and consider possible species‐specific mechanisms that control methane transport through plants.

Villa, Jorge A.↗
DOE OSTI2021

Retreat and Regrowth of the Greenland Ice Sheet During the Last Interglacial as Simulated by the CESM2‐CISM2 Coupled Climate–Ice Sheet Model

Abstract During the Last Interglacial, approximately 129 to 116 ka (thousand years ago), the Arctic summer climate was warmer than the present, and the Greenland Ice Sheet retreated to a smaller extent than its current state. Previous model‐derived and geological reconstruction estimates of the sea‐level contribution of the Greenland Ice Sheet during the Last Interglacial vary widely. Here, we conduct a transient climate simulation from 127 to 119 ka using the Community Earth System Model (CESM2), which includes a dynamic ice sheet component (the Community Ice Sheet Model, CISM2) that is interactively coupled to the atmosphere, land, ocean, and sea ice components. Vegetation distribution is updated every 500 years based on biomes simulated using a monthly climatology to force the BIOME4 equilibrium vegetation model. Results show a substantial retreat of the Greenland Ice Sheet, reaching a minimum extent at 121.9 ka, equivalent to a 3.0 m rise in sea level relative to the present day, followed by gradual regrowth. In contrast, a companion simulation employing static vegetation based on pre‐industrial conditions shows a much smaller ice‐sheet retreat, highlighting the importance of the changes in high‐latitude vegetation distribution for amplifying the ice‐sheet response.

54 ENVIRONMENTAL SCIENCES↗