Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ocean sciences”

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 181 records · Page 10

Evaluating the Trustworthiness of Explainable Artificial Intelligence (XAI) Methods Applied to Regression Predictions of Arctic Sea Ice Motion

Abstract Recent advances in explainable artificial intelligence (XAI) methods show promise for understanding predictions made by machine learning (ML) models. XAI explains how the input features are relevant or important for the model predictions. We train linear regression (LR) and convolutional neural network (CNN) models to make 1-day predictions of sea ice velocity in the Arctic from inputs of present-day wind velocity and previous-day ice velocity and concentration. We apply XAI methods to the CNN and compare explanations to variance explained by LR. We confirm the feasibility of using a novel XAI method [i.e., global layerwise relevance propagation (LRP)] to understand ML model predictions of sea ice motion by comparing it to established techniques. We investigate a suite of linear, perturbation-based, and propagation-based XAI methods in both local and global forms. Outputs from different explainability methods are generally consistent in showing that wind speed is the input feature with the highest contribution to ML predictions of ice motion, and we discuss inconsistencies in the spatial variability of the explanations. Additionally, we show that the CNN relies on both linear and nonlinear relationships between the inputs and uses nonlocal information to make predictions. LRP shows that wind speed over land is highly relevant for predicting ice motion offshore. This provides a framework to show how knowledge of environmental variables (i.e., wind) on land could be useful for predicting other properties (i.e., sea ice velocity) elsewhere. Significance Statement Explainable artificial intelligence (XAI) is useful for understanding predictions made by machine learning models. Our research establishes trustability in a novel implementation of an explainable AI method known as layerwise relevance propagation for Earth science applications. To do this, we provide a comparative evaluation of a suite of explainable AI methods applied to machine learning models that make 1-day predictions of Arctic sea ice velocity. We use explainable AI outputs to understand how the input features are used by the machine learning to predict ice motion. Additionally, we show that a convolutional neural network uses nonlinear and nonlocal information in making its predictions. We take advantage of the nonlocality to investigate the extent to which knowledge of wind on land is useful for predicting sea ice velocity elsewhere.

Hoffman, Lauren [Scripps Institution of Oceanograp↗

Estimating pixel-level uncertainty in ocean color retrievals from MODIS

The spectral distribution of marine remote sensing reflectance, R(rs), is the fundamental measurement of ocean color science, from which a host of bio-optical and biogeochemical properties of the water column can be derived. Estimation of uncertainty in these derived properties is thus dependent on knowledge of the uncertainty in satellite-retrieved R(rs) (u(c)(R(rs))) at each pixel. Uncertainty in R(rs), in turn, is dependent on the propagation of various uncertainty sources through the R(rs) retrieval process, namely the atmospheric correction (AC). A derivative-based method for uncertainty propagation is established here to calculate the pixel-level uncertainty in R(rs), as retrieved using NASA’s multiple-scattering epsilon (MSEPS) AC algorithm and verified using Monte Carlo (MC) analysis. The approach is then applied to measurements from the Moderate Resolution Imaging Spectroradiometer (MODIS) on the Aqua satellite, with uncertainty sources including instrument random noise, instrument systematic uncertainty, and forward model uncertainty. The uc(Rrs) is verified by comparison with statistical analysis of coincident retrievals from MODIS and in situ Rrs measurements, and our approach performs well in most cases. Based on analysis of an example 8-day global products, we also show that relative uncertainty in R(rs) at blue bands has a similar spatial pattern to the derived concentration of the phytoplankton pigment chlorophyll-a (chl-a), and around 7.3%, 17.0%, and 35.2% of all clear water pixels (chl-a ≤ 0.1 mg/cu.m) with valid u(c)(R(rs)) have a relative uncertainty ≤ 5% at bands 412 nm, 443 nm, and 488 nm respectively, which is a common goal of ocean color retrievals for clear waters. While the analysis shows that u(c)(R(rs)) calculated from our derivative-based method is reasonable, some issues need further investigation, including improved knowledge of forward model uncertainty and systematic uncertainty in instrument calibration.

Pixel-level uncertainty↗

A Vertically Resolved Canopy Improves Chemical Transport Model Predictions of Ozone Deposition to North Temperate Forests

Abstract Dry deposition is the second largest tropospheric ozone (O 3 ) sink and occurs through stomatal and nonstomatal pathways. Current O 3 uptake predictions are limited by the simplistic big‐leaf schemes commonly used in chemical transport models (CTMs) to parameterize deposition. Such schemes fail to reproduce observed O 3 fluxes over terrestrial ecosystems, highlighting the need for more realistic treatment of surface‐atmosphere exchange in CTMs. We address this need by linking a resolved canopy model (1D Multi‐Layer Canopy CHemistry and Exchange Model, MLC‐CHEM) to the GEOS‐Chem CTM and use this new framework to simulate O 3 fluxes over three north temperate forests. We compare results with in situ measurements from four field studies and with standalone, observationally constrained MLC‐CHEM runs to test current knowledge of O 3 deposition and its drivers. We show that GEOS‐Chem overpredicts observed O 3 fluxes across all four studies by up to 2×, whereas the resolved‐canopy models capture observed diel profiles of O 3 deposition and in‐canopy concentrations to within 10%. Relative humidity and solar irradiance are strong O 3 flux drivers over these forests, and uncertainties in those fields provide the largest remaining source of model deposition biases. Flux partitioning analysis shows that: (a) nonstomatal loss accounts for 60% of O 3 deposition on average; (b) in‐canopy chemistry makes only a small contribution to total O 3 fluxes; and (c) the CTM big‐leaf treatment overestimates O 3 ‐driven stomatal loss and plant phytotoxicity in these temperate forests by up to 7×. Results motivate the application of fully online vertically explicit canopy schemes in CTMs for improved O 3 predictions.

Vermeuel, Michael P. [Department of Soil, Water, a↗

Global metagenomics reveals plastid diversity and unexplored algal lineages

Photosynthetic organelles in eukaryotes originated through primary endosymbiosis with a cyanobacterium, an event that profoundly shaped the evolutionary landscape of the eukaryotic tree of life. Primary plastids in Archaeplastida, especially in cultivable plants and algae, contribute most to known plastid diversity. Secondary and higher-order endosymbiosis, involving eukaryotic hosts and algal endosymbionts, further spread photosynthesis among protists within the CASH lineages (Cryptophyta, Alveolata, Stramenopila, and Haptophyta). Despite various hypotheses explaining secondary plastid evolution and distribution, empirical support remains limited. Here, we employ cultivation-independent global metagenomics to expand plastid diversity and investigate plastid origins. We capture 1,027 plastid sequences, including 300 novel sequences belonging to previously unsequenced plastids and representing yet-to-be described microeukaryotes. This includes a new lineage that offers insights into plastid evolution in haptophytes and cryptophytes. Our results confirm that Archaeplastida plastids originate from an early branching cyanobacterial lineage closely related to Gloeomargaritales and identify the closest extant relative of Paulinella plastids. Additionally, our findings suggest two independent origins of secondary red-algal plastids, contributing to plastid diversity in CASH lineages and challenging the prevailing model of single secondary plastid origin. Our study highlights the importance of metagenomic data in uncovering biological diversity and advancing understanding of plastid relationships across photosynthetic eukaryotes.

59 BASIC BIOLOGICAL SCIENCES↗

Reduced CO 2 uptake and growing nutrient sequestration from slowing overturning circulation

Current Earth system models (ESMs) project dramatic slowing (28–42% by 2100) of Atlantic Meridional Overturning Circulation and Southern Meridional Overturning Circulation (SMOC) across a range of climate scenarios, with a complete shutdown of SMOC possible by year 2300. Slowing meridional overturning circulation (MOC) differentially impacts the ocean biological and solubility carbon pumps, leaving the net impact on ocean carbon uptake uncertain. Here using a suite of ESMs, we show that slowing MOC reduces anthropogenic carbon uptake by the solubility pump but increases deep-ocean storage of carbon and nutrients by the biological pump. The net effect reduces ocean uptake of anthropogenic CO 2 . The deep-ocean nutrient sequestration will increasingly depress global-scale, marine net primary production over time. As a result, MOC slowdown represents a positive feedback that could extend or intensify peak-warmth climate conditions on multi-century timescales.

58 GEOSCIENCES↗

Knocking out the carboxyltransferase interactor 1 (CTI1) in Chlamydomonas boosted oil content by fivefold without affecting cell growth

Summary The first step in chloroplast de novo fatty acid synthesis is catalysed by acetyl‐CoA carboxylase (ACCase). As the rate‐limiting step for this pathway, ACCase is subject to both positive and negative regulation. In this study, we identify a Chlamydomonas homologue of the plant carboxyltransferase interactor 1 (CrCTI1) and show that this protein interacts with the Chlamydomonas α‐carboxyltransferase (Crα‐CT) subunit of the ACCase by yeast two‐hybrid protein–protein interaction assay. Three independent CRISPR‐Cas9 mediated knockout mutants for CrCTI1 each produced an ‘enhanced oil’ phenotype, accumulating 25% more total fatty acids and storing up to fivefold more triacylglycerols (TAGs) in lipid droplets. The TAG phenotype of the crcti1 mutants was not influenced by light but was affected by trophic growth conditions. By growing cells under heterotrophic conditions, we observed a crucial function of CrCTI1 in balancing lipid accumulation and cell growth. Mutating a previously mapped in vivo phosphorylation site (CrCTI1 Ser108 to either Ala or to Asp), did not affect the interaction with Crα‐CT. However, mutating all six predicted phosphorylation sites within Crα‐CT to create a phosphomimetic mutant reduced this pairwise interaction significantly. Comparative proteomic analyses of the crcti1 mutants and WT suggested a role for CrCTI1 in regulating carbon flux by coordinating carbon metabolism, antioxidant and fatty acid β‐oxidation pathways, to enable cells to adapt to carbon availability. Taken together, this study identifies CrCTI1 as a negative regulator of fatty acid synthesis in algae and provides a new molecular brick for the genetic engineering of microalgae for biotechnology purposes.

Li, Zhongze [Aix‐Marseille Université, CEA, CNRS, ↗

Distinguishing Leptothrix and Sphaerotilus genera by an integrated genomic-phenotypic analysis supported by new Leptothrix genomes

The Sphaerotilus-Leptothrix group of bacteria includes one of the first described microorganisms, Leptothrix ochracea, an uncultured type strain, plus isolates of Leptothrix and Sphaerotilus. This group is unified by the ability to form sheaths and oxidize metals, although L. ochracea exhibits obvious ecological, morphological, and functional differences from the rest of Sphaerotilus-Leptothrix. Recently, there have been calls to combine the group into one genus, Sphaerotilus; however, these studies lacked adequate genomic representation of L. ochracea. Here, we present a comprehensive comparative genomic analysis of the Sphaerotilus-Leptothrix group, including expanded representation of L. ochracea, a closely related novel species, Leptothrix toolikensis, and two new isolates (Leptothrix mechoopdaensis). Analysis of 38 genomes resolves three phylogenetic and functional groups: the ochracea-type Leptothrix (Group 1), the mobilis-type Leptothrix (Group 2), and Sphaerotilus (Group 3). Group 1 genomes form a separate genus based on average nucleotide identity and alignment fraction. The genomes clearly diverge from the rest of Sphaerotilus-Leptothrix in phylogeny, size, and metabolic potential. Group 1 genomes are much smaller (2.59–3.04 Mb) than those of Groups 2 (4.55–6.06 Mb) and 3 (3.94–5.07 Mb), while encoding more metal oxidases and fewer carbohydrate-active enzymes. Group 2 clusters with Group 3 phylogenetically and is similar in organic carbon metabolisms but maintains more metal oxidation genes. Group 2 members lack homogeneity in phenotype and genotype, suggesting that additional isolates and genomes are needed for confident classification. However, Group 1 genomes (L. ochracea and L. toolikensis) show clear divergence, precluding their inclusion in Sphaerotilus and supporting the retention of the genus Leptothrix.

Leptothrix↗

Steric Sea Level Rise and Relationships with Model Drift and Water Mass Representation in GFDL CM4 and ESM4

Abstract Density-driven steric seawater changes are a leading-order contributor to global mean sea level rise. However, intermodel differences in the magnitude and spatial patterns of steric sea level rise exist at regional scales and often emerge during the spinup and preindustrial control integrations of climate models. Steric sea level results from an eddy-permitting climate model, GFDL CM4, are compared with a lower-resolution counterpart, GFDL-ESM4. The results from both models are examined through basin-scale heat budgets and watermass analysis, and we compare the patterns of ocean heat uptake, redistribution, and sea level differ in ocean-only [i.e., Ocean Model Intercomparison Project (OMIP)] and coupled climate configurations. After correcting for model drift, both GFDL CM4 and GFDL-ESM4 simulate nearly equivalent ocean heat content change and global sea level rise during the historical period. However, the GFDL CM4 model exhibits as much as a 40% increase in surface ocean heat uptake in the Southern Ocean and subsequent increases in horizontal export to other ocean basins after bias correction. The results suggest regional differences in the processes governing Southern Ocean heat export, such as the formation of Antarctic Intermediate Water (AAIW), Subpolar Mode Water (SPMW), and gyre transport between the two models, and that sea level changes in these models cannot be fully bias-corrected. Since the process-level differences between the two models are evident in the preindustrial control simulations of both models, these results suggest that the control simulations are important for identifying and correcting sea level–related model biases.

Krasting, John P. [a NOAA/OAR/Geophysical Fluid Dy↗

Lidar-Based Evaluation of HRRR Performance in California’s Diablo Range

The performance of the NOAA High-Resolution Rapid Refresh (HRRR) model for capturing low-level winds near a wind energy production site during summer 2019 is evaluated. This study catalogs the ability of HRRR to predict boundary layer dynamics relevant to wind energy interests over complex terrain, which has presented challenges for weather and energy forecasting. Performance is evaluated by comparing HRRR output to wind-profiling Doppler lidars at Lawrence Livermore National Laboratory Site 300. HRRR captured the diurnal profile of horizontal winds in the observed 150-m layer, despite strong underpredictions (∼4 m s −1 ) during evening and nighttime hours. These underpredictions may be a result of local speedup flows observed by the lidars, which were unresolved in HRRR due to their small spatial extent. HRRR bias magnitude relative to observations was found to be minimal during days with synoptic-scale troughs and strong 850-hPa geopotential gradients, while bias magnitude was maximal during days with synoptic ridging and weak 850-hPa geopotential gradients. To translate wind speed predictions to energy forecasting, generic turbine models were used to estimate power generation for turbines characteristic of the nearby Altamont Pass Wind Resource Area. Results show that HRRR-based energy estimates predicted daytime power generation adequately relative to lidar-based estimates with an 18-h lead time (bias magnitude < 0.4 MW from 0900 to 1400 LT) but overpredicted power during the rest of the diurnal cycle (bias > 1 MW). These results demonstrate conditions under which HRRR performs well for wind energy applications in complex terrain, while highlighting biases that require further investigation to support usage of a high-resolution model for wind energy forecasts.

Boundary layer↗

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↗

Remote Sensing and GIS data at 1km-grid over Chesapeake Bay used in “He et al. 2024, Effects of spatial variability in vegetation phenology, climate, landcover, biodiversity, topography, and soil property on soil respiration across a coastal ecosystem”

The package contains the data layers used in “He et al. 2024, Effects of spatial variability in vegetation phenology, climate, landcover, biodiversity, topography, and soil property on soil respiration across a coastal ecosystem”. The study aims to use multi-source remote sensing and GIS datasets to investigate the spatial heterogeneity and identify spatial zones with similar environmental characteristics and understand the primary driving factors affecting soil respiration within sub-ecosystems of the coastal ecosystem. We employed unsupervised hierarchical clustering analysis to identify spatial regions with distinct environmental characteristics, then determined the main driving factors using Random Forest regression and SHapley Additive exPlanations (SHAP). Spatial data layers include soil respiration, kernel Normalized Difference Vegetation Index (kNDVI) computed from Harmonized Landsat 8 and Sentinel-2 time series, climate variables from the Daymet dataset, land cover, biodiversity, topographical metrics, soil property, and tidal elevation.

54 ENVIRONMENTAL SCIENCES↗

Seedling nitrogen uptake and rhizodeposition between mycorrhizal types

Tree mycorrhizal associations are associated with patterns in N cycling and soil organic matter (SOM) storage, however, we still lack a mechanistic understanding of what tree and fungal traits drive these patterns and how they will respond to global changes in soil N availability. To address this knowledge gap, we investigated how arbuscular mycorrhizal (AM)- and ectomycorrhizal (EcM)-associated seedlings alter rhizodeposition in response to increased inorganic N acquisition. Specifically, we conducted this greenhouse experiment in a sealed labeling chamber with an enriched 13Carbon atmosphere and 15Nitrogen enriched fertilizer over the course of five months from April 2021 - August 2021. To include the variability across tree species, we grew eight species of seedlings belonging to eight families that were either arbuscular (Acer rubrum, Nyssa sylvatica, Thuja occidentalis, and Prunus seritina) or ectomycorrhizal-associated (Quercus rubra, Tilia americana, Pinus strobus, and Betula lenta). We measured rhizodeposition (mg 13C), plant N uptake from fertilizer (mg N), net soil carbon, and the abundance of mycorrhizal fungi (ITS sequencing and qPCR). We also characterized fungal (ITS2) and bacterial (16S) soil communities.The data from this project are ".csv" files that can up downloaded into a folder, and then run in the associated R markdown scripts after changing the source folder location at the top of the script. These data include raw outputs and processed files (using the R markdown files) for seedling growth and biomass, 15N content, soil 13C content, raw reads and processed file versions for fungal and bacterial ASVS, and a final summary file used for modeling. R software is needed to run these data, and the packages needed are listed at the top of the R markdown file.

54 ENVIRONMENTAL SCIENCES↗

COMPASS-FME Synoptic Site Characterization

This dataset contains soil biogeochemical and physicochemical characterization data for the COMPASS-FME synoptic sites.This dataset also contains data for the paper Patel et al. 2025 "Transition zones at the changing coastal terrestrial-aquatic interface", https://doi.org/10.1029/2025JG008978.Coastal soils are a significant but highly uncertain component of global biogeochemical cycles. These systems experience unique spatial and temporal variability in biogeochemical processes, driven by wetland-to-upland gradients and hydrological fluctuations. We studied drivers of coastal soil variability (a) at regional scales and (b) across transects from upland forest to wetland, in two contrasting regions — Lake Erie, a freshwater lacustrine system, and Chesapeake Bay, a saltwater estuarine system. Salinity-related analytes were a key driver of soil variability, not just in the saltwater system, but surprisingly, also in the freshwater system. We had hypothesized linear trends in biogeochemical parameters along the TAI – however, contrary to expectations, transition soils were not consistently intermediate between upland and wetland endmembers; the non-monotonic trends of carbon, phosphorus, iron along our transects suggest that these are key analytes to study in our regions. Rapidly changing soil factors across coastal gradients provide insights into which soil processes may act as precursors to ecosystem shifts. Our comprehensive soil characterization across the coastal transects provides essential data for mechanistic modeling of ecosystem dynamics.The data are provided as processed, csv files. Raw data and processing scripts can be accessed on GitHub (https://github.com/COMPASS-DOE/cmps-soil_characterization).A note on the nomenclature: the experimental design represents three points along the coastal gradient -- upland, transition, and wetland. "wetland" is referred to as "marsh" in the corresponding paper. The two terms can be used interchangeably for the sites in this study.

54 ENVIRONMENTAL SCIENCES↗

COMPASS-FME Terrestrial Ecosystem Manipulation to Probe the Effects of Storm Treatments (TEMPEST) Experiment Tree Inventory

This is the tree inventory (diameter, species, and live/dead status) data from the Terrestrial Ecosystem Manipulation to Probe the Effects of Storm Treatments (TEMPEST) experimental site. This manipulative, ecosystem-scale TEMPEST experiment is part of the COMPASS-FME (Coastal Observations, Mechanisms, and Predictions Across Systems and Scales: Field Measurements and Experiments; see https://compass.pnnl.gov/FME/COMPASSFME) project. It addresses the potential for freshwater and estuarine-water disturbance events to alter tree function, species composition, and ecosystem processes in a deciduous coastal forest in eastern Maryland, USA. The experiment uses a large-unit (2000 m2), un-replicated experimental design, with three 50 m × 40 m plots serving as control, freshwater, and estuarine-water treatments.This dataset includes:- An overall dataset README file.- The tree inventory data in both "wide" and "long" forms. These contain the same information but are structured differently, with the former more useful for human viewers and the latter more amenable for programmatic analyses.- A key to the species/genus codes used, which follow the U.S. Department of Agriculture's PLANTS schema (https://plants.usda.gov/).- A copy of the R code used to generate the wide- and long-form data files.All files are comma-separated value (CSV) and no special software is required to read them.

54 ENVIRONMENTAL SCIENCES↗

COMPASS-FME Synoptic Site Tree Greenhouse Gas Concentrations

These data are tree stem greenhouse gas concentrations collected from tree gas wells at some of the COMPASS-FME (Coastal Observations, Mechanisms, and Predictions Across Systems and Scales; see https://compass.pnnl.gov/) 'synoptic' sites in the Chesapeake Bay region: Moneystump (MSM), Goodwin Islands (GWI), and GCReW (GCW). The sap flow monitoring trees at these sites in the Upland (UP) and Transition (TR) zones were cored and had gas wells installed at breast height. There were also some dead standing trees cored, gas well installed, and sampled at MSM and GWI. The GCW UP samples overlap with the TEMPEST experiment control plot, so the GCW UP data was pulled from the TEMPEST page and included here. These data provide crucial information about possible pathways for the greenhouse gas (carbon dioxide and methane, CO2 and CH4 respectively) production and emission (or in the case of CH4, perhaps taken up from) the atmosphere.All data are plain text CSV (comma separated value) files and require no special software to read.Updated 2025-10-09 to fix two missing dates (lines 77 and 78 in the data file).

54 ENVIRONMENTAL SCIENCES↗

Soil microbial ecology and microbiome-metabolite linkages improve understanding of ecosystem states along terrestrial-aquatic interfaces

These data are from Bandopadhyay et al., "Soil microbial ecology and microbiome-metabolite linkages improve understanding of ecosystem states along terrestrial-aquatic interfaces". This study aims to understand the soil microbial ecology along terrestrial-aquatic interfaces of a freshwater and estuarine region and how it relates to organic matter. We analyzed soil microbial (16S rRNA gene) and organic matter (Fourier-transform ion cyclotron resonance mass spectrometry, FTICR-MS) composition from upland (forested), transition (stressed forest), and wetland positions at three sites in each of the Lake Erie (freshwater) and Chesapeake Bay (estuarine) regions. This dataset includes 16S rRNA gene amplicon data (only processed file types included here) and organic matter composition from FTICR-MS data (raw and processed files included here) from upland (forested), transition (stressed forest), and wetland positions at three sites in each of the Lake Erie and Chesapeake Bay regions. These sites are part of the COMPASS-FME project (https://compass.pnnl.gov/FME/COMPASSFME). File formats and software needed to access files: 16S rRNA gene amplicon data: These files follow the format reported here https://ess-dive.gitbook.io/amplicon-sequencing-reporting-format#updates-in-v1.0.1. As per this format, there are four file types reported: 1. Taxon tables (also called sequence-by-sample or OTU (operational taxonomic unit)/ESV (exact sequence variant) tables) : available in a .txt file format and accessible using TextEdit or MS Excel. 2. Representative sequences (also called consensus sequences) : available in a .fasta format and accessible using TextEdit. 3. Sequencing metadata : available in a MS Excel workbook file format and CSV file format 4. Bioinformatic metadata : available in a MS Excel workbook file format and CSV file format FTICR-MS data: 1. Raw data converted to a processed file with intensities of the peaks in the given samples : available in a MS Excel CSV file format 2. Processed file used in analyses and visualizations (appended as icr_long_) : available in a MS Excel CSV file format 3. Metadata file for ICR features (appended as icr_meta) : available in a MS Excel CSV file format

54 ENVIRONMENTAL SCIENCES↗

Manuscript model outputs, model source code, and figure scripts: The role of geomorphology in mediating biomass allocation impacts on salt-marsh resilience and carbon accumulation

This data package provides model source code (C++), model outputs, and figure generation scripts (R) needed to reproduce the following manuscript: Bruns, Nicholas E., Genevieve L. Noyce, and Matthew L. Kirwan. "The Role of Geomorphology in Mediating Biomass Allocation Impacts on Salt-Marsh Resilience and Carbon Accumulation." Estuarine, Coastal and Shelf Science 327 (December 2025): 109549. https://doi.org/10.1016/j.ecss.2025.109549. This manuscript investigates how geomorphology mediates the impact of biomass allocation shifts on salt marsh persistence and carbon (C) sequestration under sea level rise. We use a 1-D soil-column model (Kirwan and Mudd 2012) to perform experiments across a range of static root:shoot ratios (RSR = 1-4) spanning observed values. The model explicitly simulates interactions between tidal inundation, productivity, inorganic sediment deposition, and organic matter accumulation. Experiments determine whether geomorphic feedbacks amplify, dampen, or leave unchanged the ecosystem response to biomass allocation shifts. A first experiment uses constant sea level rise (2.5 mm/yr) to examine equilibrium dynamics and their influence on carbon accumulation rates across different suspended sediment concentrations (SSC). A second set of experiments calculates threshold sea level rise rates for marsh drowning across RSR and SSC combinations. Final experiments apply accelerating sea level rise scenarios (2000-2200) derived from NOAA projections (Sweet et al. 2022) to generate an envelope of expected responses, quantifying the importance of biomass allocation shifts on marsh carbon accumulation and persistence. All experiments are repeated across SSC ranging from 5-50 mg/L to investigate how these interactions vary in micro-tidal marshes with different sediment supplies. Package contents: * README.txt with detailed description of package contents and instructions for reproducing manuscript figures and model outputs * R scripts for generating all manuscript figures * C++ baseline model code from Kirwan and Mudd (2012) * C++ source code for 4 experimental model variants used in the manuscript, extending above baseline code * Model input files (.csv, .txt) including sea level rise scenarios * Model output files (.txt) used in manuscript analyses Temporal coverage: Model simulations span years 1900-2200, with accelerating sea level rise scenarios for 2000-2200. Key variables: Root:shoot ratio, suspended sediment concentration, carbon accumulation rate, vertical accretion rate, marsh elevation, inundation depth, threshold sea level rise rate.

54 ENVIRONMENTAL SCIENCES↗

ELM model simulations of Plum Island Ecosystems LTER low marsh site 2018-2020

Model simulations using the Department of Energy's Energy Exascale Earth System Model (E3SM) land model (ELM) with improved capabilities to represent vegetation response to salinity and inundation. The simulations were conducted for a tidal salt marsh at Plum Island Ecosystems Long Term Ecological Research (LTER) site near Rowley, Massachusetts, USA; the site is a low marsh dominated by Spartina alterniflora. The model was forced with site-specific meteorology, salinity and tidal cycles from 2018-2020. Four sets of model simulations are included and described below:1. Parameterization of the salinity response function. These simulations tested different combinations of values for optimal salinity and salinity tolerance.2. Model evaluation. This comparison conducted simulations using the default model, the salinity function only, the submergence function only, and both the salinity and submergence functions. 3. Salinity scenarios. These simulations used the 2018 salinity input data varied by -5 to +10 ppt salinity.4. Water level scenarios. These simulations used the tide height varied by -10 to +50 cm. These simulations were used to demonstrate how the salinity and submergence functions better represent carbon uptake by tidal salt marshes.The data package includes netCDF files used as forcing files for tide height and salinity, one for each year 2018-2020 at observed salinity concentrations, and an additional three forcing files in which salinity concentrations were varied 5 ppt lower, 5 ppt higher, and 10 ppt higher than the measured 2018 time series. Also included are python scripts for creating forcing files, plain text parameter and command files for running simulations, model outputs in netCDF format, and python scripts for visualizing outputs. Code for the modified E3SM model is archived in Sulman et al 2023 at doi:10.15485/1991625. More detail about files is provided in the README.md file.

54 ENVIRONMENTAL SCIENCES↗