Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vegetation dynamics modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

Changes in high-latitude surface energy balance driven by snowpack and vegetation dynamics under warmer climate

With rapid climate warming, expected changes in snowpack and vegetation will alter the seasonal surface albedo of high-latitude ecosystems. The extent to which these albedo changes may affect surface energy balances and thus soil temperatures is uncertain, but represents a potentially important component of ecosystem feedbacks to climate change. Here, we apply a well-tested process-rich ecosystem model, ecosys , to examine changes in seasonal surface albedo and soil temperature driven by climate-induced snowpack and vegetation changes across Alaska under a warmer twenty-first century climate. Under the Representative Concentration Pathway 8.5 climate change scenario, the modeled changes in surface albedo exhibited large seasonal and spatial variations. We found spring albedo decreases driven by increases in snow-free periods (>20 d) and an extended growing season length that resulted in greater gains in leaf area index (LAI) in most parts of Alaska. In contrast, we modeled increases in summer and winter albedo (despite modeled increases in LAI) across much of the boreal forest due to an increased proportion of aspen, which has a higher leaf albedo than the currently dominant black spruce. Modeled latent heat fluxes generally increase across the twenty-first century, particularly during the spring and summer. Overall, climate warming and changes in surface energy fluxes resulted in a 3.5 ± 0.50 °C increase in spatial- and annual-averaged top 10 cm surface soil temperatures across Alaskan ecosystems by the year 2100, with larger increases in tundra than boreal forest regions. We conclude that under warmer climates, seasonal variations in albedo and surface energy fluxes are particularly pronounced during the spring and summer, driven by changes in snowpack and vegetation dynamics.

54 ENVIRONMENTAL SCIENCES↗

Expansion of high-latitude deciduous forests driven by interactions between climate warming and fire: Modeling Archive

This Modeling Archive is in support of an NGEE Arctic publication. More information in the related references. The dataset is a model (ecosys) output used in a study that examined the potential future shifts in high-latitude vegetation types (PFTs) driven by changes in 21st century climate and fire. The study highlights the importance of mechanistically representing differences in plant functional traits on current and possible 21st century trajectories of Alaskan PFTs. This dataset package includes 13 files in pdf, csv, and nc formats.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Understanding water and energy fluxes in the Amazonia: Lessons from an observation‐model intercomparison

Abstract Tropical forests are an important part of global water and energy cycles, but the mechanisms that drive seasonality of their land‐atmosphere exchanges have proven challenging to capture in models. Here, we (1) report the seasonality of fluxes of latent heat (LE), sensible heat ( H ), and outgoing short and longwave radiation at four diverse tropical forest sites across Amazonia—along the equator from the Caxiuanã and Tapajós National Forests in the eastern Amazon to a forest near Manaus, and from the equatorial zone to the southern forest in Reserva Jaru; (2) investigate how vegetation and climate influence these fluxes; and (3) evaluate land surface model performance by comparing simulations to observations. We found that previously identified failure of models to capture observed dry‐season increases in evapotranspiration (ET) was associated with model overestimations of (1) magnitude and seasonality of Bowen ratios (relative to aseasonal observations in which sensible was only 20%–30% of the latent heat flux) indicating model exaggerated water limitation, (2) canopy emissivity and reflectance (albedo was only 10%–15% of incoming solar radiation, compared to 0.15%–0.22% simulated), and (3) vegetation temperatures (due to underestimation of dry‐season ET and associated cooling). These partially compensating model‐observation discrepancies (e.g., higher temperatures expected from excess Bowen ratios were partially ameliorated by brighter leaves and more interception/evaporation) significantly biased seasonal model estimates of net radiation ( R n ), the key driver of water and energy fluxes (LE ~ 0.6 R n and H ~ 0.15 R n ), though these biases varied among sites and models. A better representation of energy‐related parameters associated with dynamic phenology (e.g., leaf optical properties, canopy interception, and skin temperature) could improve simulations and benchmarking of current vegetation–atmosphere exchange and reduce uncertainty of regional and global biogeochemical models.

Restrepo‐Coupe, Natalia↗

Multiple RGB ortho-mosaics and digital surface models in 2017 and 2018 across the Lower Montane site in the East River Watershed, Colorado

Aerial imagery was collected at the Lower Montane site (Pumphouse) in the East River Watershed, Colorado during the spring, summer, and fall seasons of 2017 and 2018 to improve the understanding of seasonal vegetation dynamics and their drivers. The datasets include Red-Green-Blue (RGB) ortho-mosaics and digital surface models (DSMs) inferred from the Unoccupied Aerial System (UAS) acquired aerial RGB imagery for June 3, June 19, July 7, and August 14, 2017, and for March 14, April 26, June 1, June 18, July 6, and August 7, 2018. Real-Time Kinematic Global Positioning System (RTK-GPS) surveyed Ground control points (GCPs) were used to increase the reconstruction accuracy. The reconstructed RGB mosaics and DSMs have been trimmed to cover a similar spatial domain. The accuracy of the RGB mosaics is considered high (~10 cm). DSM accuracy is highest (~10 cm) where sufficient GCPS are available, and more difficult to assess elsewhere (see reconstruction reports for uncertainty estimates). The dataset includes a total of 20 GeoTIFF (.tif) files, 10 PDF (.pdf) files, 3 data CSV (.csv) files, and 2 metadata CSV (.csv) files. Feel free to contact the authors with any questions or collaboration interests.This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗

An experimental and kinetic modeling study of the pyrolysis of isoprene, a significant biogenic hydrocarbon in naturally occurring vegetation fires

Isoprene dominates the carbon flux emitted by vegetation and constitutes 40% of non-methane biogenic emissions worldwide. Despite pyrolysis experiments at temperatures above 1000 K showing a link between isoprene combustion and aromatic species formation, comprehensive mechanistic research on isoprene is scarce in the literature. Here, in this work, we carry out an experimental and theoretical study to build, for the first time, a chemical kinetic model describing isoprene pyrolysis. The formation of polycyclic aromatic hydrocarbon (PAH) precursor species, often observed in vegetation fire plumes, is partially explained by isoprene pyrolysis experiments and theoretical modeling. Molecular dynamics (MD) simulations unveil reaction pathways from allylic isoprenyl radicals to allene and cyclopentadiene (CPD) intermediates, two relevant species detected in the experiments. Rate constants for these identified pathways are calculated using variational transition state theory to update the kinetic model, which is validated against single-pulse shock tube (SPST), and jet-stirred reactor (JSR) experimental data in the temperature range of 850–1690 K. The kinetic model presents satisfactory agreement with the SPST experimental data, and a reaction pathway analysis shows that association of propargyl radicals results in benzene formation. The JSR pathway analysis also identifies the prominent reactions for CPD, benzene, styrene, and toluene formation. Our model does not reproduce the CPD experimental profiles, indicating that additional studies are necessary. Overall, our findings advance the understanding of isoprene pyrolysis and its related atmospheric pollutants in naturally occurring vegetation fires where smoldering and oxygen-deficient combustion processes are present.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Response of Ecosystem Productivity to High Vapor Pressure Deficit and Low Soil Moisture: Lessons Learned From the Global Eddy-Covariance Observations

Although there is mounting concern about how high vapor pressure deficit (VPD) and low soil moisture (SM) affect ecosystem productivity, their relative importance is still under debate. Here, we comprehensively quantified the relative impacts of these two factors on ecosystem gross primary production (GPP) using observations from a global network of eddy-covariance towers and two approaches (sensitivity analysis and linear regression model). Both approaches agree that a higher percentage of sites experience GPP reduction from high VPD than from low SM over the growing season. However, the constraint of high VPD and low SM on GPP reduction is tightly linked with climates and plant functional types. Humid and mesic ecosystems including forests and grasslands are dominated by VPD, while the semi-arid and arid ecosystems including shrublands and savannas are dominated by SM. The varying dominant role of these two factors on GPP is closely related to plant stomatal behavior, as predicted by a stomatal conductance model. Additionally, we highlight the non-linear impact of SM on GPP during droughts and the possible underestimation of the SM effects for deep-rooted plants when only using surface-layer SM. Our results shed light on a better understanding of the impacts of VPD and SM on vegetation productivity, with important implications for modeling the response and feedback of ecosystem dynamics to current and future climates.

54 ENVIRONMENTAL SCIENCES↗

Biogeochemistry simulations for the Salt Marsh Accretion Response to Temperature eXperiment (SMARTX)

Coastal ecosystems have been largely ignored in Earth system models but are important zones for carbon and nutrient processing. Interactions between water, microbes, soil, sediments, and vegetation are important for mechanistic representation of coastal processes and ecosystem function. To investigate the role of these feedbacks, we used a reactive transport model (PFLOTRAN) that has the capability to be connected to the Energy Exascale Earth System Model (E3SM). PFLOTRAN was used to incorporate redox reactions and track chemical species important for coastal ecosystems as well as define simple representations of vegetation dynamics. Our goal was to incorporate oxygen flux, salinity, pH, sulfur cycling, and methane production along with plant-mediated transport of gases and tidal flux. Using porewater profile and incubation data for model calibration and evaluation, we were able to create depth-resolved biogeochemical soil profiles for saltmarsh habitat and use this updated representation to simulate direct and indirect effects of elevated CO2 and temperature on subsurface biogeochemical cycling. We found that simply changing the partial pressure of CO2 or increasing temperature in the model did not fully reproduce observed changes in the porewater profile, but the inclusion of plant or microbial responses to CO2 and temperature manipulations was more accurate in representing porewater concentrations. This indicates the importance of characterizing tightly coupled vegetation-subsurface processes for developing predictive understanding and the need for measurement of plant-soil interactions on the same time scale to understand how hotspots or moments are generated.Included in this data package are PFLOTRAN input (PFLOTRAN input files and chemical database) files for simulating single column biogeochemistry, root, and tide interactions at the Global Change Research Wetland (Kirkpatrick Marsh; Edgewater, MD). The biogeochemical network includes soil organic matter decomposition, nitrogen, iron, and sulfur cycling, and methanogenesis. Reduced species can be oxidized and plant processes include oxygen and nutrient priming, methane release, and nutrient uptake.Inputs:TAI_database.dat - geochemical database for reactions, more information on database structure and variables can be found here https://www.pflotran.org/documentation/user_guide/cards/pages/geochemical_database.htmlswamp.in - input file for biogeochemical network in PFLOTRANswamp_eCO2.in - input file for biogeochemical network in PFLOTRAN with input gas partial pressures/concentrations adjusted for elevated CO2 treatmentsOutputs:swamp_obs_0.tec - hourly porewater concentrations from from multiple depths in the soil columnswamp_eCO2_obs_0.tec - hourly porewater concentrations from multiple depths in the soil column for elevated CO2 treatmentsPFLOTRAN code access: https://github.com/fmyuan/pflotran-elm-interface.git

54 ENVIRONMENTAL SCIENCES↗

Developing a Redox Network for Coastal Saltmarsh Systems in the PFLOTRAN Reaction Model

Coastal ecosystems have been largely ignored in Earth system models but are important zones for carbon and nutrient processing. Interactions between water, microbes, soil, sediments, and vegetation are important for mechanistic representation of coastal processes and ecosystem function. To investigate the role of these feedbacks, we used a reactive transport model (PFLOTRAN) that has the capability to be connected to the Energy Exascale Earth System Model (E3SM). PFLOTRAN was used to incorporate redox reactions and track chemical species important for coastal ecosystems as well as define simple representations of vegetation dynamics. Here our goal was to incorporate oxygen flux, salinity, pH, sulfur cycling, and methane production along with plant-mediated transport of gases and tidal flux. Using porewater profile and incubation data for model calibration and evaluation, we were able to create depth-resolved biogeochemical soil profiles for saltmarsh habitat and use this updated representation to simulate direct and indirect effects of elevated CO 2 and temperature on subsurface biogeochemical cycling. We found that simply changing the partial pressure of CO 2 or increasing temperature in the model did not fully reproduce observed changes in the porewater profile, but the inclusion of plant or microbial responses to CO 2 and temperature manipulations was more accurate in representing porewater concentrations. This indicates the importance of characterizing tightly coupled vegetation-subsurface processes for developing predictive understanding and the need for measurement of plant-soil interactions on the same time scale to understand how hotspots or moments are generated.

54 ENVIRONMENTAL SCIENCES↗

Building a machine learning surrogate model for wildfire activities within a global Earth system model

Abstract. Wildfire is an important ecosystem process, influencing land biogeophysical and biogeochemical dynamics and atmospheric composition. Fire-driven loss of vegetation cover, for example, directly modifies the surface energy budget as a consequence of changing albedo, surface roughness, and partitioning of sensible and latent heat fluxes. Carbon dioxide and methane emitted by fires contribute to a positive atmospheric forcing, whereas emissions of carbonaceous aerosols may contribute to surface cooling. Process-based modeling of wildfires in Earth system land models is challenging due to limited understanding of human, climate, and ecosystem controls on fire counts, fire size, and burned area. Integration of mechanistic wildfire models within Earth system models requires careful parameter calibration, which is computationally expensive and subject to equifinality. To explore alternative approaches, we present a deep neural network (DNN) scheme that surrogates the process-based wildfire model with the Energy Exascale Earth System Model (E3SM) interface. The DNN wildfire model accurately simulates observed burned area with over 90 % higher accuracy with a large reduction in parameterization time compared with the current process-based wildfire model. The surrogate wildfire model successfully captured the observed monthly regional burned area during validation period 2011 to 2015 (coefficient of determination, R2=0.93). Since the DNN wildfire model has the same input and output requirements as the E3SM process-based wildfire model, our results demonstrate the applicability of machine learning for high accuracy and efficient large-scale land model development and predictions.

58 GEOSCIENCES↗

Whole plant traits of coastal dune vegetation and implications for interactions with dune dynamics

Abstract Coastal dunes are important protective features against sea level rise and coastal storms. Interactions between dune plant aboveground structures and sediment trapping that allow for dune building and maintenance are well established. More recently, studies documenting belowground biomass for promoting erosion resistance in dominant dune species have been conducted, yet a knowledge gap remains regarding species‐specific characterization of whole plants, specifically with respect to roots, rhizomes, and belowground stems. Our objective was to quantify above‐ and belowground traits of four dominant dune grasses to document the potential for species‐specific effects on dune growth, maintenance, and erosion resistance. We examined above‐ and belowground traits among four prominent dune grasses of the Atlantic and Gulf Coasts of North America: Ammophila breviligulata , Panicum amarum , Spartina patens , and Uniola paniculata . Whole plant samples of each species were collected from the foredune at the US Army Engineer Research and Development Center's Field Research Facility in Duck, North Carolina, USA, and quantified for several above‐ and belowground traits (e.g., stem height, rhizome number and length, root surface area by diameter class, root tensile strength, and mycorrhizal percent infection). Belowground factors known to impact important dune processes, such as rhizome length, mycorrhizal percent infection, and root traits, differed substantially among species. When visualized in multivariate space, all species significantly differed in suites of above‐ and belowground traits. When considering belowground only, Ammophila and Spartina were similar, despite differences in biomass allocation. Species separated along axes related to mycorrhizal association, biomass allocation, and root construction. The four co‐occurring dune grass species were dissimilar in suites of plant traits. Belowground trait differences were driven by those describing root construction, biomass allocation, and mycorrhizal infection. Dissimilarity in above‐ and belowground suites of traits may demonstrate different approaches for surviving the dune environment. Incorporating belowground traits into modeling will enhance predictions of dune response to climate change through interactions between vegetation and dune dynamics that facilitate coastal resistance and resilience.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of obstacle modelling approaches for resource assessment and small wind turbine siting: case study in the northern Netherlands

Abstract. Growth in adoption of distributed wind turbines for energy generation is significantly impacted by challenges associated with siting and accurate estimation of the wind resource. Small turbines, at hub heights of 40 m or less, are greatly impacted by terrestrial obstacles such as built structures and vegetation that can cause complex wake effects. While some progress in high-fidelity complex fluid dynamics (CFD) models has increased the potential accuracy for modelling the impacts of obstacles on turbulent wind flow, these models are too computationally expensive for practical siting and resource assessment applications. To understand the efficacy of available models in situ, this study evaluates classic and commonly used methods alongside new state-of-the-art lower-order models derived from CFD simulations and machine learning approaches. This evaluation is conducted using a subset of an extensive original dataset of measurements from more than 300 operational wind turbines in the northern Netherlands. The results show that data-driven methods (e.g. machine learning and statistical modelling) are most effective at predicting production at real sites with an average error in annual energy production of 2.5 %. When sufficient data may not be available de novo to support these data-driven approaches, models derived from high-fidelity simulations show promise and reliably outperform classic methods. On average these models have 6.3 %–11.5 % error compared with 26 % for classic methods and 27 % baseline error for reanalysis data without obstacle correction. While more performant on average, these methods are also sensitive to the quality of obstacle descriptions and reanalysis inputs.

17 WIND ENERGY↗

The persistence and conversion of coastal foredune and swale vegetation community distributions 63 years later

Vegetation shifts can directly alter habitat dynamics and indirectly impact habitat stability relative to disturbance response. Barrier island dune habitats exhibit spatiotemporally dynamic topography that is affected by vegetation. However, vegetation distribution data can be rare and vegetation persistence is largely unknown such that species turnover over in communities can occur unnoticed. This is true despite concerns and documented cases of woody encroachment related to climate change in these ecogeomorphic habitats where physical stability to resist storm erosion varies with vegetation distribution and density. In 1956 and 1957, the vegetation of Island Beach State Park, NJ, was mapped as a permanent record for subsequent ecological study. In June 2020, we remapped the vegetation of 5.4 of 17 km north to south, seaward of the thicket community boundary. We maintained the same classification system as the historic record, physically mapping vegetation patches via GPS. We quantified changes in thicket, heather, and grass community distribution between the two time periods. Habitat persistence and conversion varied in the 63 years. Heath communities saw 80%–97% habitat loss in conversion to woody thicket. Conversely, thicket community distribution drastically increased, replacing heath where it was previously prevalent. This represents the first known documented instances of woody encroachment for Morella pensylvanica. When they did not expand, thicket communities receded landward or underwent turnover to an invasive species. Woody species of interest for dune stabilization occupied areas of similar habitat characteristics. Dune vegetation distribution persistence from 1957 to 2020 was relatively consistent and stable with the exception of heath habitat conversion. Barrier island stability is directly related to vegetation stability such that understanding where community shifts might occur over time can aid in managing and modeling efforts surrounding these dynamic ecogeomorphic habitats.

54 ENVIRONMENTAL SCIENCES↗

DUET - Distribution of Understory using Elliptical Transport: A mechanistic model of leaf litter and herbaceous spatial distribution based on tree canopy structure

Heterogeneity in surface fuels produced by overstory trees and understory vegetation is a major driver of fire behavior and ecosystem dynamics. Previous attempts at predicting tree leaf and needle litter accumulation over time have been constrained in scope to probabilistic models that consider a limited number of key factors influencing tree litter dispersal patterns and decomposition processes. Here we present a mechanistic model for estimating variation in surface fuels called the Distribution of Understory using Elliptical Transport (DUET). DUET uses a pre-generated voxelated canopy array and distributes the leaf litter from the trees in an elliptical shape based on tree species characteristics, wind data, and location-specific features. DUET then calculates grass growth between trees and around the litter patterns and sets moisture levels based on relative humidity. In this study, we use DUET to simulate spatial distributions of leaf litter and illustrate model sensitivity to several parameters, providing inputs for testing fuel spatial heterogeneity impacts on fire behavior using process-based models.

54 ENVIRONMENTAL SCIENCES↗

Integrating New Arctic Plant Functional Types in a Land Surface Model Using Above- and Belowground Field Observations: Modeling Archive

This Modeling Archive is in support of an NGEE Arctic publication "Integrating New Arctic Plant Functional Types in a Land Surface Model Using Above- and Belowground Field Observations". We used field observations of biomass and leaf traits (Salmon et al., 2019b) across a gradient of plant communities on the Seward Peninsula in northwest Alaska (Breen et al., 2020; Iversen et al., 2017) to replace the original ELM configuration with nine, arctic-specific PFTs. Original arctic PFTs in the model were (1) broadleaf deciduous boreal shrub and (2) C3 arctic grass. New PFTs were (1) lichen, (2) bryophyte, (3) evergreen dwarf shrub, (4) deciduous dwarf shrub, (5) deciduous low shrub, (6) deciduous low to tall shrub, (7) alder shrub, (8) forb, and (9) graminoid. New PFTs were parameterized and evaluated using site-level trait and biomass measurements (Salmon et al., 2019b), root C:N ratios and rooting depth patterns (Salmon et al., 2019c), and root biomass (Salmon et al., 2020). This archive contains forcing data and model output from ELM simulations conducted at the Kougarok Hillslope site (Kougarok Road Mile Marker 64). Simulations were conducted for the Kougarok Hillslope site using meteorological driving data from the Scenarios Network for Alaska and Arctic Planning (SNAP) downscaled climate projection dataset using NCAR-CCSM forcing (Bieniek et al., 2020 dataset; Walsh et al., 2018). The archive includes three versions of model simulations representing a progression from the original model configuration to an Arctic-specific configuration (see Output section below). Simulation 1 used grid cell data from a global E3SM configuration, including the fractional area of E3SM PFTs assigned to the grid cell containing the Kougarok Hillslope site in global simulations. This simulation did not distinguish between different plant communities on the Kougarok Hillslope, but instead used a single point simulation to represent the entire area. Simulation 2 used default E3SM PFT definitions combined with adjusted depth to bedrock for plant communities with shallow rocky layers to represent the role of abiotic soil factors in driving site differences (based on depth measurements of Iversen et al., 2019a), and adjusted the relative areas of the broadleaf deciduous boreal shrub and C3 arctic grass PFTs to reflect the observed spatial coverage of shrub and non-shrub PFTs across Kougarok Hillslope plant communities. Simulation 3 used the new arctic PFT definitions and parameterizations based on the vegetation types present at the site, along with the adjusted soil depths used in Simulation 2. Included are *.pdf, *.nc (NetCDF), *tar.gz, *.py, and *.bash files. We would like to thank Mary's Igloo Native Corporation for allowing us to perform this research on their land and giving us the opportunity to share this research with their community. Thanks to Shawn Serbin and Alistair Rogers for helpful discussions. Holly Vander Stel, Joanne Childs, and Stan Wullschleger helped with collection of field data. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Nutrient Dynamics in a Coupled Terrestrial Biosphere and Land Model (ELM-FATES-CNP)

We present a representation of nitrogen and phosphorus cycling in the Functionally Assembled Terrestrial Ecosystem Simulator, a demographic vegetation model within the Energy Exascale Earth System land model. This representation is modular, and designed to allow testing of multiple hypothetical approaches for carbon-nutrient coupling in plants. Novel model hypotheses introduced in this work include, (a) the controls on plant acquisition of aqueous mineralized nutrients in the soil and (b) fairly straight forward methods of allocating nutrients to specific plant organs and their losses through live plant turnover as well as litter fluxes generated through plant mortality. This combines the new with pre-existing hypotheses (such as nitrogen fixation and soil decomposition) into a system that can accommodate plant-soil dynamics for a large number of size- and functional-type-resolved plant cohorts within a time-since-disturbance-resolved ecosystem. Root uptake of nutrients is governed by fine root biomass, and plants vary in their fine root biomass allocation in order to balance carbon and nutrient limitations to growth. We test the sensitivity of the model to a wide range of parameter variations and structural representations, and in the context of observations at Barro Colorado Island, Panama. A key model prediction is that plants in the high-light-availability canopy positions allocate more carbon to fine roots than plants in low-light understory environments, given the widely different carbon versus nutrient constraints of these two niches within a given ecosystem. This model provides a basis for exploring carbon-nutrient coupling with vegetation demography within Earth system models.

54 ENVIRONMENTAL SCIENCES↗

Considering coasts: Adapting terrestrial models to characterize coastal wetland ecosystems

The Energy Exascale Earth System Model (E3SM) simulates fully coupled processes and interactions among water, energy, carbon and nutrient cycles. E3SM connects vegetation and soil dynamics through nutrient uptake, plant production, litterfall and decomposition as a function of abiotic parameters (e.g. temperature and moisture). However, E3SM is designed to characterize terrestrial ecosystems and connects land and open ocean systems using a single streamflow transport term, ignoring the complex dynamics of energy, water, carbon, and nutrients in coastal systems. The goals of our project were to: (1) Parameterize a point version of E3SM to capture coastal wetland habitats and (2) Determine marsh community responses to increased temperature and elevated CO 2 . We adapted a version of the E3SM land model, previously configured to represent forested bog hydrology to a coastal ecosystem using datasets from field experiments conducted at the Smithsonian Environmental Research Center's Global Change Research Wetland (GCReW). Tidal forcing in a marsh environment was simulated using a two-column system in which the columns are connected by lateral hydrologic flows. One column simulates interactions between vegetation and soil while a second column simulates variation in water level (both tidal and sea level rise). The updated model captures many aspects of the field experiments, showing that plant community responses to environmental change are non-linear, non-additive and different between plant types. Elevated CO 2 treatments increased C 3 plant biomass more than C 4 (33% vs 17%). Temperature exacerbated CO 2 responses in C 3 plants (0 °C: 26%, 5.1 °C: 56%). We were more successful at characterizing C 3 than C 4 responses and simulating above rather than belowground biomass production. Furthermore, the next steps will include updates to key physiological parameters such as root:shoot carbon allocation and the addition of mechanistic feedbacks between vegetation and biogeochemical processes.

54 ENVIRONMENTAL SCIENCES↗

Snow Depth Datasets for Snodgrass Catchment, Colorado, Water Year 2022-2023

This data package presents snow depths data from distributed temperature probes at 18 locations near Snodgrass catchment, Colorado. These data show that snow melt-out dates are approximately one or two weeks later under evergreen forests compared to other vegetation types even at the same elevation. These data were collected to understand how snowmelt heterogeneity impacts headwater hydrology, including streamflow and groundwater levels. They were also used to compare with process-based model simulations of snow depth to evaluate whether the model accurately represents snowmelt dynamics and their effects on headwater hydrology. Snow_DTPs_locations.csv includes all probes locations and their associated elevation and vegetation types. Snow_Depth_Snodgrass_WY2022_2023.csv includes processed snow depths datasets for Water Year (WY) 2022 and 2023. This dataset also includes a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata and a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type. Several probes have recordings for WY 2021.

54 ENVIRONMENTAL SCIENCES↗

Hydrology controls thermokarst and alters carbon cycling and methane emissions in peatlands near the southern limit of permafrost

Permafrost peatlands store vast amounts of frozen carbon across northern landscapes. When ground ice melts, surface subsidence produces thermokarst landforms that expand wetlands at the edges of permafrost plateaus. Thermokarst represents an accelerating climate feedback, but uncertainties remain about how ground ice, hydrology, and vegetation interact to shape landscape change and carbon fluxes. We extended the process-based model ecosys to simulate thermokarst dynamics in laterally coupled 2D transects at a well-characterized boreal peatland site in Canada’s Northwest Territories. After benchmarking against site observations, we varied ground ice content and hydrologic boundary conditions across ranges typical near the southern permafrost limit. Simulations revealed distinct degradation regimes governed by the elevation difference between the frost table and the external water table. Rates of lateral retreat, the thaw-driven encroachment of wetlands into adjacent plateaus, ranged from 0 to >2 m yr −1 under identical weather forcing, consistent with observations and highlighting the strong role of WT and ground ice. Simulated vegetation dynamics indicate that black spruce mortality cannot be explained by anoxia alone, pointing to additional stressors such as root damage, pathogens, or physical destabilization. Despite large hydrologic shifts, net ecosystem CO 2 exchange remained a slight sink after collapse, while methane (CH 4 ) emissions rose by one to two orders of magnitude. As a result, lateral retreat substantially increases the greenhouse warming potential of permafrost peatlands (1.7 million km 2 in area), with simulated emissions of 0.1–10 Mt CO 2 -eq decade −1 depending on hydrology and retreat rates. These results underscore the need to account for both ground ice and hydrologic dynamics when assessing thermokarst-driven climate feedbacks.

carbon cycling↗