Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Tidal”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Exploring the origin of the distance bimodality of stars in the periphery of the Small Magellanic Cloud with APOGEE and Gaia

ABSTRACT The Magellanic Cloud system represents a unique laboratory for study of both interacting dwarf galaxies and the ongoing process of the formation of the Milky Way and its halo. We focus on one aspect of this complex, three-body interaction – the dynamical perturbation of the Small Magellanic Cloud (SMC) by the Large Magellanic Cloud (LMC), and specifically potential tidal effects on the SMC’s eastern side. Using Gaia astrometry and the precise radial velocities (RVs) and multielement chemical abundances from Apache Point Observatory Galactic Evolution Experiment (APOGEE-2) Data Release 17, we explore the well-known distance bimodality on the eastern side of the SMC. Through estimated stellar distances, proper motions, and RVs, we characterize the kinematics of the two populations in the bimodality and compare their properties with those of SMC populations elsewhere. Moreover, while all regions explored by APOGEE seem to show a single chemical enrichment history, the metallicity distribution function (MDF), of the ‘far’ stars on the eastern periphery of the SMC is found to resemble that for the more metal-poor fields of the western periphery, whereas the MDF for the ‘near’ stars on the eastern periphery resembles that for stars in the SMC Centre. The closer eastern periphery stars also show RVs (corrected for SMC rotation and bulk motion) that are, on average, approaching us relative to all other SMC populations sampled. We interpret these trends as evidence that the near stars on the eastern side of the SMC represent material pulled out of the central SMC as part of its tidal interaction with the LMC.

Almeida, Andres (ORCID:0009000007332479)↗
DOE OSTI2024

Properties of a fading AGN from SDSS-IV MaNGA

ABSTRACT We identify a fading AGN SDSS J220141.64+115124.3 from the internal Product Launch-11 (MPL-11) in Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. The central region with a projected radius of $\sim$2.4 kpc is characterized as LINER-like line ratios while the outskirts extended to $\sim$15 kpc show Seyfert-like line ratios. The ${\rm{[O {\small III}]}}$$\lambda$5007 luminosity of the Seyfert regions is a factor of 37 (2) higher than the LINER regions without (with) dust attenuation correction, suggesting that the AGN activity decreases at least $\sim$8 $\times$ 10$^3$ yr ($\sim$2.4 kpc/light-speed) ago. We model the emission line spectra in the central region with double Gaussian components (a narrow core and a broad wing) and analyse the properties of each component. The narrow core component mostly co-rotates with the stellar disc, whereas the broad wing component with a median of the velocity dispersion $\sim$300 km s$^{-1}$ is related to a wind outflow. The kinematic position angle (PA) of the ionized gas shows an $\sim 20^{\circ }$ twist from the galaxy centre to 1.5 effective radius. The median of the PA difference between the gas and stellar components is as large as $\sim 50^{\circ }$ within 0.4 effective radius. The tidal feature in DESI image and star–gas misalignment suggest this galaxy is a merger remnant. Combining all these observational results as well as public available X-ray and MIR luminosities, we confirm this is a fading AGN, the merger process kick-started the central engine to quasar phase which ionized gas composed of tidal debris, and now the activity of the central black hole decreases. The discontinuity in ${\rm{[O {\small III}]}}$$\lambda$5007 flux and EQW maps is due to multiple AGN outbursts triggered by merger remnant gas inflows.

Mo, Hao (ORCID:0000000234430768)↗
DOE OSTI2025

Auriga Streams – I: disrupting satellites surrounding Milky Way-mass haloes at multiple resolutions

In a hierarchically formed Universe, galaxies accrete smaller systems that tidally disrupt as they evolve in the host’s potential. We present a complete catalogue of disrupting galaxies accreted onto Milky Way-mass haloes from the Auriga suite of cosmological magnetohydrodynamic zoom-in simulations. We classify accretion events as intact satellites, stellar streams, or phase-mixed systems based on automated criteria calibrated to a visually classified sample, and match accretions to their counterparts in haloes re-simulated at higher resolution. Most satellites at the present day have lost substantial amounts of stellar mass – 67 per cent have $f_\text{bound} < 0.97$ (our threshold of lost stellar mass to no longer be considered intact), while 53 per cent satisfy a more stringent $f_\text{bound} < 0.8$. Streams typically outnumber intact systems, contribute a smaller fraction of overall accreted stars, and are substantial contributors at intermediate distances from the host centre ($\sim$0.1 to $\sim 0.7R_\text{200m}$, or $\sim$35 to $\sim$250 kpc for the Milky Way). We also identify accretion events that disrupt to form streams around massive intact satellites instead of the main host. Streams are more likely than intact or phase-mixed systems to have experienced pre-processing, suggesting this mechanism is important for setting disruption rates around Milky Way-mass haloes. All of these results are preserved across different simulation resolutions, though we do find some hints that satellites disrupt more readily at lower resolution. The Auriga haloes suggest that disrupting satellites surrounding Milky Way-mass galaxies are the norm and that a wealth of tidal features waits to be uncovered in upcoming surveys.

galaxies: haloes↗
DOE OSTI2025

Direct nonparametric multimessenger constraints on the equation of state of cold dense nuclear matter

We utilize the now substantial amount of astrophysical observations of neutron stars (NSs), along with perturbative quantum chromodynamics (pQCD) calculations at high density, to directly constrain the NS equation of state (EOS). To this end, we construct nonparametric EOS priors by using Gaussian processes trained on 75 EOSs, which include models with either hadrons, hyperons, or quarks at high densities. We create a prior using the full EOS sample (model agnostic), and one prior for each EOS family to test model discrimination. We introduce a novel inference approach, which allows the simultaneous sampling of intrinsic and extrinsic parameters of binary NS mergers, as well as a nonparametric equation of state. We showcase this method in a Bayesian updating scheme by first performing a complete analysis of the binary NS merger event GW170817 with minimal assumptions, and sequentially adding information from x-ray and radio NS observations, along with pQCD calculations. Besides providing standard constraints, such as the pressure at twice nuclear saturation density 𝑝⁡(2⁢𝜌 sat ) = 4.3$^{+0.6}_{−0.6}$ × 10 34 dyne/cm 2 , at 95% confidence level, for the model agnostic prior, our methodology shows how the choice of EOS families used in conditioning changes the inferred astrophysical properties of the EOS, namely tidal deformability and maximum supported NS mass. We find hyperonic priors predicting higher tidal deformabilities for a 1.4⁢𝑀 ⊙ NS, and hadronic priors being preferred by the considered astrophysical data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Interface Modes in Inspiralling Neutron Stars: A Gravitational-Wave Probe of First-Order Phase Transitions

At the extreme densities in neutron stars, a phase transition to deconfined quark matter is anticipated. Yet masses, radii, and tidal deformabilities offer only indirect measures of a first-order phase transition, requiring many detections to resolve or being ineffective observables if the discontinuity exists at lower densities. Here, we report on a smoking-gun gravitational-wave signature of a first-order transition: the resonant tidal excitation of an interface mode. Using relativistic perturbation theory with an equation-of-state family informed by chiral effective field theory, we show that such a resonance may be detectable with next-generation interferometers and possibly already with LIGO A+ for sufficiently loud events.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Modifying PyUltraLight to model scalar dark matter with self-interactions

Here we introduce a modification of the pysiultralight code that models the dynamical evolution of ultralight axionlike scalar dark matter fields. Our modified code, pysiultralight, adds a quartic, self-interaction term to reflect the one which arises naturally in axionlike particle models. Using a particle mass of 10−22 eV/c 2 , we show that pysiultralight produces spatially oscillating solitons, exploding solitons, and collapsing solitons which prior analytic work shows will occur with attractive self-interactions. Using our code we calculate the oscillation frequency as a function of soliton mass and equilibrium radius in the presence of attractive self-interactions. We show that when the soliton mass is below the critical mass ($M_c$ = $\frac{\sqrt{3}}{2} M_{max}$) described by Chavanis and the initial radius is within a specific range, solitons are unstable and explode. We test the maximum mass criteria described by Chavanis and Chavanis and Delfini for a soliton to collapse when attractive self-interactions are included. We also analyze both binary soliton collisions and a soliton rotating around a central mass with attractive and repulsive self-interactions. We find that when attractive self-interactions are included, the density profiles get distorted after a binary collision. We also find that a soliton is less susceptible to tidal stripping when attractive self-interactions are included. We find that the opposite is true for repulsive self-interactions in that solitons would be more easily tidally stripped. Including self-interactions might therefore influence the survival timescales of infalling solitons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Numerical relativity simulations of prompt collapse mergers: Threshold mass and phenomenological constraints on neutron star properties after GW170817

We determine the threshold mass for prompt (no bounce) black hole formation in equal-mass neutron star (NS) mergers using a new set of 227 numerical relativity simulations. We consider 23 phenomenological and microphysical finite-temperature equations of state (EOS), including models with hyperons and first-order phase transitions to deconfined quarks. We confirm the existence of EOS-insensitive relations between the threshold mass, binary tidal parameter at the threshold (Λth), maximum mass of nonrotating NSs, and radii of reference mass NSs. We combine the EOS-insensitive relations, phenomenological constraints on NS properties, and observational data from GW170817 to derive an improved lower limit on radii of maximum mass and a 1.6 M ⊙ NS of 9.81 and 10.90 km, respectively. We also constrain the radius and quadrupolar tidal deformability (Λ) of a 1.4 M ⊙ NS to be larger than 10.74 km and 172, respectively. We consider uncertainties in all independent parameters—fitting coefficients as well as GW170817 masses while reporting the range of radii constraints. We discuss an approach to constrain the upper as well as lower limit of NS maximum mass using future binary NS detections and their identification as prompt or delayed collapse. With future observations, it will be possible to derive even tighter constraints on the properties of matter at and above nuclear density using the method proposed in this work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Reconstructing the neutron star equation of state from observational data via automatic differentiation

Neutron star observables like masses, radii, and tidal deformability are direct probes to the dense matter equation of state (EoS). Here, a novel deep learning method that optimizes an EoS in the automatic differentiation framework of solving inverse problems is presented. The trained neural network EoS yields narrow bands for the relationship between the pressure and speed of sound as a function of the mass density. The results are consistent with those obtained from conventional approaches and the observational bound on the tidal deformability inferred from the gravitational wave event, GW170817.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Detectability of subsolar mass neutron stars through a template bank search

Here, we study the detectability of gravitational-wave signals from subsolar-mass binary neutron star systems by the current generation of ground-based gravitational-wave detectors. We find that finite size effects from large tidal deformabilities of the neutron stars and lower merger frequencies can significantly impact the sensitivity of the detectors to these sources. By simulating a matched-filter based search using injected binary neutron star signals with tidal deformabilities derived from physically motivated equations of state, we calculate the reduction in sensitivity of the detectors. We conclude that the loss in sensitive volume can be as high as 78.4% for an equal mass binary system of chirp mass 0.17M ⊙ , in a search conducted using binary black hole template banks. We use this loss in sensitive volume, in combination with the results from the search for subsolar-mass binaries conducted on data collected by the LIGO-Virgo observatories during their first three observing runs, to obtain a conservative upper limit on the merger rate of subsolar-mass binary neutron stars. Since the discovery of a low-mass neutron star would provide new insight into formation mechanisms of neutron stars and further constrain the equation of state of dense nuclear matter, our result merits a dedicated search for subsolar-mass binary neutron star signals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Impact of O4 future detections on the determination of the dense matter equations of state

Here, in view of the next LIGO-Virgo-KAGRA Observing period O4 (to start in Spring 2023), we address the question of the ability of the interferometers network to discriminate among different neutron stars equation of states better than what was possible with the observation of the binary neutron stars merger GW170817. We show that the observation of an event similar to GW170817 during O4 would allow us to resolve the dimensionless effective tidal deformability $\tilde{Λ}$ within an uncertainty 7 times better than the one obtained in O2. Thanks to the expected increase in sensitivities, we show that any GW170817-like single-event within a distance of 100 Mpc would imply significantly improved constraints of the neutron stars equations of state. We also illustrate the important impact of the noise in the analysis of the signal, showing how it can impact the effective tidal deformability probability density function for large signal-to-noise ratio.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Coastal Carbon Library and Atlas: Open source soil data and tools supporting blue carbon research and policy

Abstract Quantifying carbon fluxes into and out of coastal soils is critical to meeting greenhouse gas reduction and coastal resiliency goals. Numerous ‘blue carbon’ studies have generated, or benefitted from, synthetic datasets. However, the community those efforts inspired does not have a centralized, standardized database of disaggregated data used to estimate carbon stocks and fluxes. In this paper, we describe a data structure designed to standardize data reporting, maximize reuse, and maintain a chain of credit from synthesis to original source. We introduce version 1.0.0. of the Coastal Carbon Library, a global database of 6723 soil profiles representing blue carbon‐storing systems including marshes, mangroves, tidal freshwater forests, and seagrasses. We also present the Coastal Carbon Atlas, an R‐shiny application that can be used to visualize, query, and download portions of the Coastal Carbon Library. The majority (4815) of entries in the database can be used for carbon stock assessments without the need for interpolating missing soil variables, 533 are available for estimating carbon burial rate, and 326 are useful for fitting dynamic soil formation models. Organic matter density significantly varied by habitat with tidal freshwater forests having the highest density, and seagrasses having the lowest. Future work could involve expansion of the synthesis to include more deep stock assessments, increasing the representation of data outside of the U.S., and increasing the amount of data available for mangroves and seagrasses, especially carbon burial rate data. We present proposed best practices for blue carbon data including an emphasis on disaggregation, data publication, dataset documentation, and use of standardized vocabulary and templates whenever appropriate. To conclude, the Coastal Carbon Library and Atlas serve as a general example of a grassroots F.A.I.R. (Findable, Accessible, Interoperable, and Reusable) data effort demonstrating how data producers can coordinate to develop tools relevant to policy and decision‐making.

Holmquist, James R.↗
DOE OSTI2023

Salinity–induced limits to mangrove canopy height

Aim: Mangrove canopy height is a key metric to assess tidal forests' resilience in the face of climate change. In terrestrial forests, tree height is primarily determined by water availability, plant hydraulic design, and disturbance regime. However, the role of water stress remains elusive in tidal environments, where saturated soils are prevalent, and salinity can substantially affect the soil water potential. Location: Global. Time Period: The canopy height dataset provides a global snapshot of the maximum mangrove height geographical distribution for the year 2000. Climate and environmental variables extend over the period 1970–2018. Major Taxa Studied: Mangroves. Methods: We use global observations of maximum canopy height, species richness, air temperature, and seawater salinity—a proxy of soil water salt concentration—to explore the causal link between salinity and mangrove stature. Results: Our findings suggest that salt stress limits mangrove height. High salinity favours more salt-tolerant species, narrowing the spectrum of viable traits. Highly salt-tolerant mangroves have evolved to cope with high salt concentrations in the soil, but this adaptation comes at a cost. They typically have lower rates of photosynthesis and growth, resulting in reduced productivity and smaller stature compared to more salt-sensitive mangrove species. This suggests a causal link between salinity, biodiversity, and tree height, where high salinity selects for more salt-tolerant species that tend to be less productive and shorter. Conclusions: We hypothesize that the salinity-induced limit to mangrove canopy height is the direct result of a reduction of primary productivity, an increment in the risk of xylem cavitation, and an indirect consequence of the decrease in biodiversity. As sea-level rise enhances coastal salinisation, failure to account for these effects can lead to incorrect estimates of future carbon stocks in Tropical coastal ecosystems and endanger preservation efforts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Lab-Scale, Closed-Loop Experimental Characterization, Model Refinement, and Validation of a Hydrokinetic Energy-Harvesting Ocean Kite

This paper presents a study wherein we experimentally characterize the dynamics and control system of a lab-scale ocean kite, and then refine, validate, and extrapolate this model for use in a full-scale system. Ocean kite systems, which harvest tidal and ocean current resources through high-efficiency cross-current motion, enable energy extraction with an order of magnitude less material (and cost) than stationary systems with the same rated power output. However, an ocean kite represents a nascent technology that is characterized by relatively complex dynamics and requires sophisticated control algorithms. In order to characterize the dynamics and control of ocean kite systems rapidly, at a relatively low cost, the authors have developed a lab-scale, closed-loop prototyping environment for characterizing tethered systems, whereby 3D printed systems are tethered and flown in a water channel environment. While this system has been shown to be capable of yielding similar dynamic characteristics to some full-scale systems, there are also fundamental limitations to the geometric scales and flow speeds within the water channel environment, making many other real-world scenarios impossible to replicate from the standpoint of dynamic similarity. Here, to address these scenarios, we show how the lab-scale framework is used to refine and validate a scalable dynamic model of a tethered system, which can then be extrapolated to full-scale operation. In this work, we present an extensive case study of this model refinement, validation, and extrapolation on an ocean kite system intended for operation in the Gulf Stream or similar current environments.

42 ENGINEERING↗
DOE OSTI2024

Multiple microbial guilds mediate soil methane cycling along a wetland salinity gradient

ABSTRACT Estuarine wetlands harbor considerable carbon stocks, but rising sea levels could affect their ability to sequester soil carbon as well as their potential to emit methane (CH 4 ). While sulfate loading from seawater intrusion may reduce CH 4 production due to the higher energy yield of microbial sulfate reduction, existing studies suggest other factors are likely at play. Our study of 11 wetland complexes spanning a natural salinity and productivity gradient across the San Francisco Bay and Delta found that while CH 4 fluxes generally declined with salinity, they were highest in oligohaline wetlands (ca. 3-ppt salinity). Methanogens and methanogenesis genes were weakly correlated with CH 4 fluxes but alone did not explain the highest rates observed. Taxonomic and functional gene data suggested that other microbial guilds that influence carbon and nitrogen cycling need to be accounted for to better predict CH 4 fluxes at landscape scales. Higher methane production occurring near the freshwater boundary with slight salinization (and sulfate incursion) might result from increased sulfate-reducing fermenter and syntrophic populations, which can produce substrates used by methanogens. Moreover, higher salinities can solubilize ionically bound ammonium abundant in the lower salinity wetland soils examined here, which could inhibit methanotrophs and potentially contribute to greater CH 4 fluxes observed in oligohaline sediments. IMPORTANCE Low-level salinity intrusion could increase CH4 flux in tidal freshwater wetlands, while higher levels of salinization might instead decrease CH4 fluxes. High CH4 emissions in oligohaline sites are concerning because seawater intrusion will cause tidal freshwater wetlands to become oligohaline. Methanogenesis genes alone did not account for landscape patterns of CH4 fluxes, suggesting mechanisms altering methanogenesis, methanotrophy, nitrogen cycling, and ammonium release, and increasing decomposition and syntrophic bacterial populations could contribute to increases in net CH4 flux at oligohaline salinities. Improved understanding of these influences on net CH4 emissions could improve restoration efforts and accounting of carbon sequestration in estuarine wetlands. More pristine reference sites may have older and more abundant organic matter with higher carbon:nitrogen compared to wetlands impacted by agricultural activity and may present different interactions between salinity and CH4. This distinction might be critical for modeling efforts to scale up biogeochemical process interactions in estuarine wetlands.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Neutrino emission from binary neutron star mergers: characterising light curves and mean energies

Abstract Neutrinos are copiously emitted by neutron star mergers, due to the high temperatures reached by dense matter during the merger and its aftermath. Neutrinos influence the merger dynamics and shape the properties of the ejecta, including the resulting r -process nucleosynthesis and kilonova emission. In this work, we analyse neutrino emission from a large sample of binary neutron star merger simulations in Numerical Relativity, covering a broad range of initial masses, nuclear equation of state and viscosity treatments. We extract neutrino luminosities and mean energies, and compute quantities of interest such as the peak values, peak broadnesses, time averages and decrease time scales. We provide a systematic description of such quantities, including their dependence on the initial parameters of the system. We find that for equal-mass systems the total neutrino luminosity (several $$10^{53}{\hbox {erg}~{\hbox {s}}^{-1}}$$ 10 53 erg s - 1 ) decreases as the reduced tidal deformability increases, as a consequence of the less violent merger dynamics. Similarly, tidal disruption in asymmetric mergers leads to systematically smaller luminosities. Peak luminosities can be twice as large as the average ones. Electron antineutrino luminosities dominate (initially by a factor of 2-3) over electron neutrino ones, while electron neutrinos and heavy flavour neutrinos have similar luminosities. Mean energies are nearly constant in time and independent on the binary parameters. Their values reflect the different decoupling temperature inside the merger remnant. Despite present uncertainties in neutrino modelling, our results provide a broad and physically grounded characterisation of neutrino emission, and they can serve as a reference point to develop more sophisticated neutrino transport schemes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Global Technical, Economic, and Feasible Potential of Renewable Electricity

Renewable electricity generation will need to be rapidly scaled to address climate change and other environmental challenges. Doing so effectively will require an understanding of resource availability. We review estimates for renewable electricity of the global technical potential, defined as the amount of electricity that could be produced with current technologies when accounting for geographical and technical limitations as well as conversion efficiencies; economic potential, which also includes cost; and feasible potential, which accounts for societal and environmental constraints. We consider utility-scale and rooftop solar photovoltaics, concentrated solar power, onshore and offshore wind, hydropower, geothermal electricity, and ocean (wave, tidal, ocean thermal energy conversion, and salinity gradient energy) technologies. We find that the reported technical potential for each energy resource ranges over several orders of magnitude across and often within technologies. Therefore, we also discuss the main factors explaining why authors find such different results. According to this review and on the basis of the most robust studies, we find that technical potentials for utility-scale solar photovoltaic, concentrated solar power, onshore wind, and offshore wind are above 100 PWh/year. Hydropower, geothermal electricity, and ocean thermal energy conversion have technical potentials above 10 PWh/year. Rooftop solar photovoltaic, wave, and tidal have technical potentials above 1 PWh/year. Salinity gradient has a technical potential above 0.1 PWh/year. The literature assessing the global economic potential of renewables, which considers the cost of each renewable resource, shows that the economic potential is higher than current and near-future electricity demand. Fewer studies have calculated the global feasible potential, which considers societal and environmental constraints. While these ranges are useful for assessing the magnitude of available energy sources, they may omit challenges for large-scale renewable portfolios.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Modeled sub-hourly nitrate concentrations in subsurface water across a salt marsh system in Elkhorn Slough, California

Salt marshes are dynamic hydrologic systems where terrestrial groundwater, terrestrial surface water, and seawater mix due to bi-directional flows and pressure gradients. Salt marshes are known to remove excess nitrate from water thereby improving water quality in coastal systems. However, our knowledge of the short-term processes (e.g., intra-tidal) in salt marshes is limited. Further, we lack knowledge about how these short-term dynamics may change as a function of seasonality. To address this knowledge gap, we develope a sampling methodology that allowed us to continuously measure nitrate in subsurface water at three different depths across three marsh elevations at sub-hourly time steps. Here we present a data table with the results of nitrate in porewater across the three salt marsh elevations (information in the data set). We analyzed nitrate with a s::can spectrophotometer. We used the UV-Vis absorbance of the samples to calculate nitrate concentrations using a local calibration library. More information about the methodology can be found in Grande et al (2023). Seasonal and Tidal Variations in Hydrologic Inputs Drive Salt Marsh Porewater Nitrate Dynamics, Hydrological Processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

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↗