SEARCH · Engineering Papers

Results for “Process sensitivity”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

DOE OSTI2021

NRAP-Open-IAM Analytical Reservoir Model: Development and Testing

Geological carbon sequestration (GCS) is a key technology for reducing global carbon dioxide (CO 2 ) emissions. Over the last decade, the U.S. Department of Energy has invested in understanding the science base, developing practical implementation methods, and demonstrating secure GCS technologies to mitigate the environmental impacts associated with the atmospheric release of CO 2 . As part of the National Risk Assessment Partnership, a systems-level risk assessment tool, called the NRAP-Open-IAM, has been developed to conduct risk assessment and enable safe operations at a GCS site. The current NRAP-Open-IAM contains a simple reservoir model component that calculates the evolution of CO 2 saturation and fluid pressure in a storage reservoir during CO 2 injection operations. This report presents the development and testing of a new analytical reservoir reduced-order model (ROM), which is extended from an existing semi-analytical model for estimation of CO 2 and brine leakage along legacy wells, and enhances the capability of the NRAP-Open-IAM to simulate more types of reservoir conditions. The developed model is validated against three reference studies, and the results indicate that the new ROM predicts the behavior of the two-phase fluids (brine and injected CO 2 ) well and is applicable to different reservoir simulation boundary conditions (i.e., constant pressure boundary and infinite-acting boundary) without a priori user specification of the boundary type. Sensitivity analysis for a set of model parameters is performed using 4,000 synthetic cases prepared via a fully automated process and using machine-learning-based feature selection. The stochastic analysis identifies gravitational number (i.e., ratio of gravitational forces to viscous force) and distance between the injection well and observation location as the most impactful parameters for matching the pressure and CO 2 saturation, respectively, between the numerical simulations and the ROM. This report details the possible ROM uncertainties and serves as a guide for users to understand the use and limitations of this ROM. The code implementation of the model will be released as a module within the NRAP-Open-IAM.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Comparative Study for Pressure Gain Combustion-Gas Turbine System Performance in NGCC Configurations

The Department of Energy (DOE) has a goal to increase natural gas combined cycle (NGCC) plant efficiency to greater than 70 percent on a lower heating value (LHV) basis. Improvement in gas turbine (GT) efficiency is key achieving this goal. While many of the traditional improvement methods (higher firing temperature, advanced materials, improved cooling methods, etc.) may increase current GT and NGCC efficiencies incrementally, combining these traditional advancements with breakthrough technologies such as pressure gain combustion (PGC) may provide significant system efficiency gains. PGC is a combustion process that results in a net pressure increase across the combustor (up to pressure ratio of 2.2) rather than the traditional (<5%) pressure drop. There are multiple PGC technologies being developed that create the pressure rise phenomena, including resonant pulse combustors (RPCs), pulse detonation engines (PDEs), wave rotor combustor (WRC) engines and rotating detonation engines (RDEs). In this study, scoping analyses of PGC concepts applied to NGCC power plants were performed by combining simplified analytical models for the RPC, PDE, WRC, and RDE types of pressure gain processes with a GT model to evaluate impacts of using each of the PGC technologies on key GT performance parameters. Sensitivity analyses were also performed. Ultimately, the PGC technologies were found to have the potential to increase simple cycle efficiency by 4–6 percent (absolute) and NGCC efficiency by 2–3 percent (absolute) based on the results of the current study.

03 NATURAL GAS↗

DOE OSTI2024

Measuring Muon Antineutrino Charged-Current Interactions Without Mesons in the Final State, in the NOvA Near Detector

NOvA is a long-baseline neutrino experiment based at Fermilab in the US, with the primary aim of measuring neutrino and antineutrino oscillations. This will enhance our understanding of electroweak interactions by measuring the neutrino mixing angles, CP-violating phase and neutrino mass ordering. To measure these oscillations, we first need to have a deep understanding of how neutrinos and antineutrinos interact with matter. Antineutrino interaction cross sections are, at present, particularly poorly constrained, and processes such as meson exchange currents are not well understood in the antineutrino sector. This analysis will develop a cross-section measurement of muon antineutrino interactions without mesons (e.g. pions or kaons) in the final state, in the NOvA near detector. A high-statistics, high-purity sample is obtained through a cut-based selection process implementing machine learning techniques. The sample is dominated by quasi-elastic and meson exchange current interactions which are sensitive to nuclear effects such as Final-State Interactions. The cross section will be extracted as a function of the incoming neutrino energy and the kinematics of the outgoing particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Leveraging the signature of heterotrophic respiration on atmospheric CO 2 for model benchmarking

Spatial and temporal variations in atmospheric carbon dioxide (CO 2 ) reflect large-scale net carbon exchange between the atmosphere and terrestrial ecosystems. Soil heterotrophic respiration (HR) is one of the component fluxes that drive this net exchange, but, given observational limitations, it is difficult to quantify this flux or to evaluate global-scale model simulations thereof. Here, we show that atmospheric CO 2 can provide a useful constraint on large-scale patterns of soil heterotrophic respiration. We analyze three soil model configurations (CASA-CNP, MIMICS, and CORPSE) that simulate HR fluxes within a biogeochemical test bed that provides each model with identical net primary productivity (NPP) and climate forcings. We subsequently quantify the effects of variation in simulated terrestrial carbon fluxes (NPP and HR from the three soil test-bed models) on atmospheric CO 2 distributions using a three-dimensional atmospheric tracer transport model. Our results show that atmospheric CO 2 observations can be used to identify deficiencies in model simulations of the seasonal cycle and interannual variability in HR relative to NPP. In particular, the two models that explicitly simulated microbial processes (MIMICS and CORPSE) were more variable than observations at interannual timescales and showed a stronger-than-observed temperature sensitivity. Our results prompt future research directions to use atmospheric CO 2 , in combination with additional constraints on terrestrial productivity or soil carbon stocks, for evaluating HR fluxes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Measuring Muon Antineutrino Charged-Current Interactions without Mesons in the Final State, in the NOvA Near Detector

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

A sensitive and specific genetically-encoded potassium ion biosensor for in vivo applications across the tree of life

Potassium ion (K + ) plays a critical role as an essential electrolyte in all biological systems. Genetically-encoded fluorescent K + biosensors are promising tools to further improve our understanding of K + -dependent processes under normal and pathological conditions. Here, we report the crystal structure of a previously reported genetically-encoded fluorescent K + biosensor, GINKO1, in the K + -bound state. Using structure-guided optimization and directed evolution, we have engineered an improved K + biosensor, designated GINKO2, with higher sensitivity and specificity. We have demonstrated the utility of GINKO2 for in vivo detection and imaging of K + dynamics in multiple model organisms, including bacteria, plants, and mice.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Sequential Learning of Active Subspaces

This study shows that in recent years, active subspace methods (ASMs) have become a popular means of performing subspace sensitivity analysis on black-box functions. Naively applied, however, ASMs require gradient evaluations of the target function. In the event of noisy, expensive, or stochastic simulators, evaluating gradients via finite differencing may be infeasible. In such cases, often a surrogate model is employed, on which finite differencing is performed. When the surrogate model is a Gaussian process (GP), we show that the ASM estimator is available in closed form, rendering the finite-difference approximation unnecessary. We use our closed-form solution to develop acquisition functions focused on sequential learning tailored to sensitivity analysis on top of ASMs. We also show that the traditional ASM estimator may be viewed as a method of moments estimator for a certain class of GPs. We demonstrate how uncertainty on GP hyperparameters may be propagated to uncertainty on the sensitivity analysis, allowing model-based confidence intervals on the active subspace. Our methodological developments are illustrated on several examples.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Jet charge with global event shapes: Probing quark flavor dynamics

We propose measuring the jet electric charge of jet regions, defined within the framework of global event shapes, as a probe of quark flavor dynamics within the nucleon and the hadronization process. In particular, we consider a measurement of the jet region charge while simultaneously keeping track of the energy flow throughout the event, as characterized by the global event shape. As a concrete example, we focus on the measurement of the (Q), the jet charge of the jet region (𝐽) defined within the framework of the 1-Jettiness global event shape (𝜏 1 ) for the Deep Inelastic Scattering (DIS) process, 𝑒 − +𝑝→𝑒 − +𝐽+𝑋, with unpolarized or longitudinally polarized protons. The 1-Jettiness distribution, binned according to jet charge, allows for enhanced quark flavor separation of the initial state unpolarized or polarized PDFs. On the other hand, the jet charge distribution binned by 1-Jettiness can serve as a probe of quark flavor dynamics in the final state hadronization process. We derive a factorization theorem for simultaneous measurements of 𝜏 1 and Q in the resummation region, 𝜏 1 ≪𝑃 𝐽 𝑇 , where 𝑃 𝐽 𝑇 denotes the transverse momentum of the jet region. The correlation between the jet region charge and the charge of the struck quark is expected to be the strongest in the resummation region. The factorization theorem contains a new universal , generalizing the standard jet function to include a jet charge measurement. Therefore these universal functions can be extracted from a global analysis of N-jettiness and thrust at 𝑒 + ⁢𝑒 − colliders. We provide simulation studies to demonstrate the sensitivity of the 1-Jettiness jet charge observable to quark flavor dynamics in nucleon structure and explore the possibility of probing the final state hadronization process. This observable is well-suited for applications with existing HERA data and the future Electron-Ion Collider (EIC).

Mantry, Sonny [University of North Georgia, Dahlon↗

DOE OSTI2026

Real-Time GPU-Accelerated OFDR With an Integrated Auxiliary Interferometer

A GPU-accelerated optical frequency domain reflectometry (OFDR) system with an improved integrated auxiliary interferometer is proposed. Unlike conventional approaches that require separate auxiliary interferometers and multiple detection channels, the proposed OFDR system embeds this functionality directly into the signal via an intentional beat component. This enables self-calibration of laser nonlinearity while maintaining a cost-effective hardware configuration. Building on this simplified configuration, the system leverages GPU acceleration with an NVIDIA RTX 4070 Ti to achieve real-time performance, delivering high-throughput signal processing for continuous OFDR interrogation. The signal processing pipeline comprises signal capture, resampling for nonlinearity compensation, and frequency shift computation, all optimized for parallel execution. Hardware benchmarking demonstrates substantial acceleration over CPU implementations, achieving up to a 45× speedup for resampling and frequency shift computations and enabling processing latencies below 30 ms. Thermal response validation is conducted under two complementary scenarios: localized heating using a water bath and cryogenic-temperature conditions using liquid nitrogen. Under localized heating, the system achieves an accuracy of 0.249 °C with a thermal sensitivity of 5.971 GHz/°C, while cryogenic-temperature validation demonstrates a frequency shift response with a sensitivity of 2.383 GHz/°C and an accuracy of 2.04 °C. The high acceleration of the proposed GPU-accelerated OFDR system and its accuracy are achieved by exploiting CUDA-based stride indexing, enabling efficient parallel segmentation and processing of large datasets without additional memory copies. The benchmarking results confirm the robustness, accuracy, and deployability of the proposed OFDR system across a wide temperature range, establishing it as a practical platform for real-time distributed fiber sensing in structurally dynamic environments.

Harb, Salah [Lawrence Berkeley National Laboratory↗

DOE OSTI2025

Resolving Mesoscale Convective Systems: Grid Spacing Sensitivity in the Tropics and Midlatitudes

Abstract Mesoscale convective systems (MCSs) are a critical global water cycle component and drive extreme precipitation events in tropical and midlatitude regions. However, simulating deep convection remains challenging for modern numerical weather and climate models due to the complex interactions of processes from microscales to synoptic scales. Recent models with kilometer‐scale horizontal grid spacings offer notable improvements in simulating deep convection compared to coarser‐resolution models. Still, deficiencies in representing key physical processes, such as entrainment, lead to systematic biases. Additionally, evaluating model outputs using process‐oriented observational data remain difficult. This study presents an ensemble of MCS simulations with spanning the deep convective gray zone ( from 12 km to 125 m) in the Southern Great Plains of the U.S. and the Amazon Basin. Comparing these simulations with Atmospheric Radiation Measurement (ARM) wind profiler observations, we find greater sensitivity in the Amazon Basin compared to the Great Plains. Convective drafts converge structurally at sub‐kilometer scales, but some deficiencies remain. In both regions, simulated up and downdrafts are too deep and extreme downdrafts are not strong enough. Furthermore, Amazonian updrafts are too strong. Overall, we observe higher sensitivity in the tropics, including an artificial buildup in vertical kinetic energy at scales of , suggesting a need for 250 m in this region. Nevertheless, bulk convergence—agreement of storm‐average statistics—is achievable with kilometer‐scale simulations within a 10% error margin with 1 km providing a good balance between accuracy and computational cost.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Probing the Structure and Dynamics of Fluid Mixtures in Porous Materials Through Ultrafast Vibrational Spectro-Microscopy and Many-Body Molecular Dynamics

This research program focused on the characterization of both molecular structure and dynamics of aqueous solutions in various metal-organic frameworks (MOFs) and organic polymers that have recently been proposed for applications in water treatment technologies. This was accomplished by integrating ultrafast vibrational spectroscopy and microscopy, which are sensitive to local environments, molecular orientation, and dynamical couplings, with many-body molecular dynamics (MB-MD) simulations, which enable realistic modeling of aqueous systems and porous materials under different thermodynamic conditions. Our studies allowed for identifying the physical mechanisms and characterizing the underlying molecular interactions between water and host materials that determine the adsorption and transport processes of water in various MOFs and β-cyclodextrin polymers – prototypical examples of porous materials and polymers, that have recently been proposed for technological applications in water purification. By combining bulk- and surface-sensitive, spatially resolved ultrafast vibrational spectroscopy with MB-MD simulations, it was possible to gain broad insight into the structure and mobility of water under heterogeneous confinement at various length scales. In particular, measurements of surface- and bulk-sensitive vibrational spectra at single crystal level, integrated with MB-MD simulations were used to characterize domain-specific, structure–binding affinity relationships. Besides providing specific information about the behavior of aqueous solutions in prototypical MOFs and organic polymers for potential applications in water treatment technologies, our studies also provided fundamental insights into the properties of aqueous solutions in confined heterogeneous environments, which have broad implications in many areas, ranging from heterogeneous catalysis to oil recovery, ion transport processes, and reaction-diffusion processes in crowded systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Sensitivity evaluation of the Kudryavtsev permafrost model

Modeling is an important way to assess current and future permafrost spatial distribution and dynamics, especially in data poor areas like the Arctic region. Here, we evaluate a physics-based analytical model, Kudryavtsev's active layer model, which is widely used because it has relatively few data requirements. This model was recently incorporated into a component modeling toolbox, allowing for coupled modeling of permafrost and geomorphic processes over geological timescales. However, systematic quantitative assessment of the influence of its controlling parameters on permafrost temperature and active layer thickness predictions has not been undertaken before. In this work, we investigate the sensitivity of the Kudryavtsev's active layer model by Monte Carlo simulations to generate probability distributions for input parameters and compare predictions with a comprehensive benchmark dataset of in-situ permafrost observations over entire Alaska. Predicted permafrost surface temperature is highly dependent on mean annual air temperature (r = 0.78 on average), annual temperature amplitude (–0.41), and winter-averaged snow thickness (0.30). Uncertainty of predicted permafrost temperature is relatively small (RMSE = 1 °C), when air temperature and snow depth are well constrained. Similarly, RMSE between simulated and observed ALT at stations is ~0.08 m. However, under given air temperature and snow conditions, soil water content bias can significantly affect modeled active layer thickness (RMSE = 0.1 m or 40% of the observed active layer thickness). If soil water content has a large bias, improvements in other parameters may not significantly improve the active layer predictions of the Kudryavtsev's model.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

The NtrYX Two-Component System Regulates the Bacterial Cell Envelope

ABSTRACT Activity of the NtrYX two-component system has been associated with important processes in diverse bacteria, ranging from symbiosis to nitrogen and energy metabolism. In the facultative alphaproteobacterium Rhodobacter sphaeroides , loss of the two-component system NtrYX results in increased lipid production and sensitivity to some known cell envelope-active compounds. In this study, we show that NtrYX directly controls multiple properties of the cell envelope. We find that the response regulator NtrX binds upstream of cell envelope genes, including those involved in peptidoglycan biosynthesis and modification and in cell division. We show that loss of NtrYX impacts the cellular levels of peptidoglycan precursors and lipopolysaccharide and alters cell envelope structure, increasing cell length and the thickness of the periplasm. Cell envelope function is also disrupted in the absence of NtrYX, resulting in increased outer membrane permeability. Based on the properties of R. sphaeroides cells lacking NtrYX and the target genes under direct control of this two-component system, we propose that NtrYX plays a previously undescribed, and potentially conserved, role in the assembly, structure, and function of the cell envelope in a variety of bacteria. IMPORTANCE The bacterial cell envelope provides many important functions. It protects cells from harsh environments, serves as a selective permeability barrier, houses bioenergetic functions, defines sensitivity to antibacterial agents, and plays a crucial role in biofilm formation, symbiosis, and virulence. Despite the important roles of this cellular compartment, we lack a detailed understanding of the biosynthesis and remodeling of the cell envelope. Here, we report that the R. sphaeroides two-component signaling system NtrYX is a previously undescribed regulator of cell envelope processes, providing evidence that it is directly involved in controlling transcription of genes involved in cell envelope assembly, structure, and function in this and possibly other bacteria. Thus, our data report on a newly discovered process used by bacteria to assemble and remodel the cell envelope.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

PETSc TSAdjoint: A Discrete Adjoint ODE Solver for First-Order and Second-Order Sensitivity Analysis

Here, we present a new software system PETSc TSAdjoint for first-order and second order adjoint sensitivity analysis of time-dependent nonlinear differential equations. The derivative calculation in PETSc TSAdjoint is essentially a high-level algorithmic differentiation process. The adjoint models are derived by differentiating the timestepping algorithms and implementing them based on the parallel infrastructure in PETSc. Full differentiation of the library code, including MPI routines, is avoided, and users do not need to derive their own adjoint models for their specific applications. PETSc TSAdjoint can compute the first-order derivative, that is, the gradient of a scalar functional, and the Hessian-vector product, which carries second-order derivative information, while requiring minimal input (a few callbacks) from the users. The adjoint model employs optimal checkpointing schemes in a manner that is transparent to users. Finally, usability, efficiency, and scalability are demonstrated through examples from a variety of applications.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Synchrotron-based techniques for characterizing STCH water-splitting materials

Understanding the role of oxygen vacancy–induced atomic and electronic structural changes to complex metal oxides during water-splitting processes is paramount to advancing the field of solar thermochemical hydrogen production (STCH). The formulation and confirmation of a mechanism for these types of chemical reactions necessitate a multifaceted experimental approach, featuring advanced structural characterization methods. Synchrotron X-ray techniques are essential to the rapidly advancing field of STCH in part due to properties such as high brilliance, high coherence, and variable energy that provide sensitivity, resolution, and rapid data acquisition times required for the characterization of complex metal oxides during water-splitting cycles. X-ray diffraction (XRD) is commonly used for determining the structures and phase purity of new materials synthesized by solid-state techniques and monitoring the structural integrity of oxides during water-splitting processes (e.g., oxygen vacancy–induced lattice expansion). X-ray absorption spectroscopy (XAS) is an element-specific technique and is sensitive to local atomic and electronic changes encountered around metal coordination centers during redox. While in operando measurements are desirable, the experimental conditions required for such measurements (high temperatures, controlled oxygen partial pressures, and H 2 O) practically necessitate in situ measurements that do not meet all operating conditions or ex situ measurements. Here, we highlight the application of synchrotron X-ray scattering and spectroscopic techniques using both in situ and ex situ measurements, emphasizing the advantages and limitations of each method as they relate to water-splitting processes. The best practices are discussed for preparing quenched states of reduction and performing synchrotron measurements, which focus on XRD and XAS at soft (e.g., oxygen K-edge, transition metal L-edges, and lanthanide M-edges) and hard (e.g., transition metal K-edges and lanthanide L-edges) X-ray energies. The X-ray absorption spectra of these complex oxides are a convolution of multiple contributions with accurate interpretation being contingent on computational methods. The state-of-the-art methods are discussed that enable peak positions and intensities to be related to material electronic and structural properties. Through careful experimental design, these studies can elucidate complex structure–property relationships as they pertain to nonstoichiometric water splitting. A survey of modern approaches for the evaluation of water-splitting materials at synchrotron sources under various experimental conditions is provided, and available software for data analysis is discussed.

08 HYDROGEN↗

DOE OSTI2023

Data-Driven Affinely Adjustable Robust Volt/VAr Control

Recent years have seen the increasing proliferation of distributed energy resources with intermittent power outputs, posing new challenges to the voltage management in distribution networks. To this end, this paper proposes a data-driven affinely adjustable robust Volt/VAr control (AARVVC) scheme, which modulates the smart inverter’s reactive power in an affine function of its active power, based on the voltage sensitivities with respect to real/reactive power injections. To achieve a fast and accurate estimation of voltage sensitivities, we propose a data-driven method based on deep neural network (DNN), together with a rule-based bus-selection process using the bidirectional search method. Our method only uses the operating statuses of selected buses as inputs to DNN, thus significantly improving the training efficiency and reducing information redundancy. Finally, a distributed consensus-based solution, based on the alternating direction method of multipliers (ADMM), for the AARVVC is applied to decide the inverter’s reactive power adjustment rule with respect to its active power. Only limited information exchange is required between each local agent and the central agent to obtain the slope of the reactive power adjustment rule, and there is no need for the central agent to solve any (sub)optimization problems. Finally, numerical results on the modified IEEE-123 bus system validate the effectiveness and superiority of the proposed data-driven AARVVC method.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2025

Mono‐Materials Created by Engineering a Continuum of P3HB Stereomicrostructures in a One‐Step Catalytic Process

Multi-material products that combine multiple complementary polymers can create products with desired performance but present challenges to end-of-life (EoL) management. The emerging mono-material product design based on a single polymer type addresses the fundamental EoL issue, but challenges of delivering vastly tunable material properties by the single polymer still remain. Here, we introduce a simple strategy to produce biodegradable poly(3-hydroxybutyrate) (P3HB) materials with a wide range of material properties by engineering a stereomicrostructure continuum, achieved through polymerizing diastereomeric mixtures of racemic and meso-dimethyl diolides at various feed ratios with a single catalyst. This one-step, one-pot process produces biodegradable P3HB mono-materials ranging from rigid to flexible thermoplastics, to tough thermoplastic elastomers, to a pressure-sensitive adhesive (PSA), which have been then combined to fabricate prototype all-P3HB PSA tapes, demonstrating the feasibility of designing mono-material products through engineering polymer stereomicrostructures.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Evaluating sulfurization as a blue carbon sink in a southern California salt marsh

Blue carbon ecosystems such as seagrass meadows, mangrove forests, and salt marshes are important carbon sinks that can store carbon for millennia. Recently, organic matter sulfurization and pyritization have been proposed as mechanisms of net carbon storage in blue carbon ecosystems. At our study site, organic sulfur that is resistant to acid hydrolysis (protokerogen) is an order of magnitude less abundant than pyrite sulfur, suggesting a dominance of pyritization over sulfurization. The C/N ratios and carbon isotope compositions suggest that nearly half of total organic carbon and ≥ 80% of protokerogen is composed of marsh plant material. Sediment protokerogen appears to be sulfurized based on its low δ 34 S values (− 10‰), abundance of disulfides, and higher S/C ratio (~ 1.0%) relative to potential biogenic sulfur sources. However, the interpretation of protokerogen δ 34 S values is complicated by the wide range in sulfur isotope compositions of marsh plants. Evidence for sulfurization occurs within the shallowest sediments across different vegetation zones, yielding consistent products, while pyritization appears to be more sensitive to alterations in sediment redox conditions. Based on organic sulfur and pyrite content, sulfurization may be a more spatially consistent process than pyritization, with implications for carbon storage. The relative abundance of pyrite and protokerogen organic sulfur indicates that pyritization is favored at our study site, but this is likely to vary across the spectrum of blue carbon ecosystems.

54 ENVIRONMENTAL SCIENCES↗