skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Benthic macrofaunal carbon fluxes and environmental drivers of spatial variability in a large coastal-plain estuary
Abstract. While the importance of carbon cycling in estuaries is increasingly recognized, the role of benthic macrofauna remains poorly quantified due to limited spatial and temporal resolution in biomass measurements. Here, we ask: (1) To what extent do benthic macrofauna contribute to estuarine carbon cycling via respiration and calcification? and (2) How well can routinely collected environmental variables predict their biomass? We analyzed data from 8128 benthic samples collected from the Chesapeake Bay between 1995 and 2022 and estimated associated carbon fluxes using empirical relationships. We then used generalized additive models to relate observed and modeled environmental variables to the biomass. Biomass was highest in the upper mainstem of the Bay (Upper Bay) and upper Potomac River Estuary, the largest tidal tributary of the Bay. In the Upper Bay, benthic macrofauna respired 18 %–45 % of the estimated organic carbon supply. Calcification-driven alkalinity reduction reached 6.31 ± 2.84 mol m−2 yr−1 in the Potomac River Estuary, aligning with prior estimates of alkalinity sinks in the tributary and highlighting the potential importance of calcifying fauna in alkalinity dynamics. Estimated CO2 production in the Upper Bay from benthic respiration and calcification (151 g C m−2 yr−1) also exceeded observed air–sea CO2 fluxes (74.5 g C m−2 yr−1). Generalized additive models revealed that low salinity, moderate dissolved oxygen, and elevated nitrate best predicted high-biomass zones, with the three predictors explaining 52 % of biomass deviance. These predictive relationships offer a pathway to estimate macrofaunal biomass and associated carbon fluxes in systems where direct biomass measurements are sparse. Our findings demonstrate that benthic macrofauna play a substantial and spatially structured role in estuarine carbon cycling. Incorporating their contributions into estuarine biogeochemical models will improve predictions of ecosystem responses to environmental and anthropogenic change.  more » « less
Award ID(s):
1852428
PAR ID:
10655276
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
European Geophysical Union
Date Published:
Journal Name:
Biogeosciences
Volume:
22
Issue:
23
ISSN:
1726-4189
Page Range / eLocation ID:
7769 to 7795
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Variations in estuarine carbonate chemistry can have critical impacts on marine calcifying organisms, yet the drivers of this variability are difficult to quantify from observations alone, due to the strong spatiotemporal variability of these systems. Terrestrial runoff and wetland processes vary year to year based on local precipitation, and estuarine processes are often strongly modulated by tides. In this study, a 3D-coupled hydrodynamic-biogeochemical model is used to quantify the controls on the carbonate system of a coastal plain estuary, specifically the York River estuary. Experiments were conducted both with and without tidal wetlands. Results show that on average, wetlands account for 20–30% of total alkalinity (TA) and dissolved inorganic carbon (DIC) fluxes into the estuary, and double-estuarine CO2outgassing. Strong quasi-monthly variability is driven by the tides and causes fluctuations between net heterotrophy and net autotrophy. On longer time scales, model results show that in wetter years, lower light availability decreases primary production relative to biological respiration (i.e., greater net heterotrophy) resulting in substantial increases in CO2outgassing. Additionally, in wetter years, advective exports of DIC and TA to the Chesapeake Bay increase by a factor of three to four, resulting in lower concentrations of DIC and TA within the estuary. Quantifying the impacts of these complex drivers is not only essential for a better understanding of coastal carbon and alkalinity cycling, but also leads to an improved assessment of the health and functioning of coastal ecosystems both in the present day and under future climate change. 
    more » « less
  2. ABSTRACT Meiofauna, organisms smaller than 1 mm, are the most abundant and diverse invertebrates inhabiting the world's ocean floor but their contribution to benthic oxygen demand is still poorly constrained. This knowledge is crucial for understanding seabed respiration, global marine carbon, and oxygen cycles, which are relevant to all nutrient cycling and energy flows in the ecosystem. It is common to predict meiofauna respiration based on their biomass or volume, which are difficult to quantify, and thus meiofauna are rarely included in biogeochemistry studies. In addition, it is still unknown how well the generalized allometric relations describe all meiofauna respiration. Therefore, we used a novel approach specially developed for single meiofauna respiration measurements to derive the respiration rates of 10 meiofauna groups in two marine and one brackish coastal muddy environments under oxic and hypoxic conditions, representing natural sediment conditions. Our estimates suggest that large ostracods and juveniles of macrofauna (e.g., bivalves, trumpet worms, and priapulids) had the highest individual respiration rates. Meiofauna community as a whole contributed 3–33% to sediment oxygen uptake. However, the most important contributors to the overall sediment oxygen uptake were nematodes and foraminifera which had lower respiration rates but were highly abundant. Therefore, out of more than 22 meiofauna phyla, we recommend that nematode and foraminifera respiration, which contributes 3–30% (total 3–33%) to sediment oxygen uptake, should be taken into consideration in any estimations of benthic oxygen and carbon cycles. 
    more » « less
  3. Understanding the relative contributions of aboveground and belowground processes to soil accretion and carbon density may explain carbon sequestration rates in mangroves across different coastal environmental settings. We reformulated the nutrient mangrove model (NUMAN) by improving algorithms and uncertainty analysis using literature values and site-specific observations to evaluate the relative contributions of organic and inorganic sedimentation for three mangrove sites with marked soil fertility gradients reflected by nitrogen-to-phosphorus (N:P) ratios including Shark River (N:P = 28), Rookery Bay (N:P = 54–78), and Taylor Slough (N:P = 102) in south Florida. NUMAN 2.0 considers cellulose as a refractory organic-matter sub-pool and simultaneously incorporates coarse-root inputs to soil formation. The model simulation also captures root necromass accumulation. Monte Carlo (MC) simulations (N = 1000 per site) were conducted to capture uncertainty by treating five key parameters as random variables: lignin content in fine, coarse, and large roots; inorganic sediment loading; and root biomass at the surface. With robust mass balancing of organic matter, NUMAN 2.0 generates precise predictions of surface accretion and carbon density. NUMAN 2.0 simulations estimated mean (standard deviation) soil carbon sequestration rates at 130.1 (55.4) for Shark River, 72.5 (3.7) for Rookery Bay, and 130.0 (83.9) g m-2 yr-1 for Taylor Slough, compared to field values of 123.0, 86.0, and 108.8 (8.7) g m-2 yr-1 , respectively. Simulation experiments with NUMAN 2.0 suggest that belowground organic matter dominates soil formation and carbon sequestration generally in coastal environmental settings with little allochthonous input such as carbonate settings, while wood litterfall should dominate soil organic matter in top 10 cm in estuaries, and bays. 
    more » « less
  4. Different sources and processes contribute to pCO2 and CO2 exchange with the atmosphere in the rivers and floodplains of the Amazon basin. We measured or estimated pCO2, CO2 fluxes with the atmosphere, planktonic community respiration (PCR), and environmental and landscape variables along the Negro and Amazon-Solimões rivers during different periods of the fluvial hydrological cycle. Values of pCO2 ranged from 307 to 7,527 μatm, while CO2 fluxes ranged from -9.3 to 1,128 mmol m-2 d-1 in the Amazon-Solimões basin. In the Negro basin, pCO2 values ranged from 648 to 6,526 μatm, and CO2 fluxes from 35 to 1,025 mmol m-2 d-1. In a general linear model including data from Negro and Amazon-Solimões basins, seasonal and spatial variation in flooded vegetated habitat area, dissolved oxygen, depth and water temperature explained 85% of surface pCO2 variation. Levels of pCO2 varied with inundation extent, with higher pCO2 values occurring in periods with greater water depth and inundation area, and lower dissolved oxygen concentrations and water temperatures. In a separate analysis for the Amazon-Solimões river and floodplains, ecosystem type (lotic or lentic), hydrological period, water temperature, dissolved oxygen, depth and dissolved phosphorus explained 83% of pCO2 variation. Our results demonstrate the influence of alluvial floodplains and seasonal variations in their limnological characteristics on the pCO2 levels in river channels of the lowland Amazon. 
    more » « less
  5. Rates of ecosystem metabolic properties, such as plankton community respiration, can be used as an assessment of the eutrophication state of a waterbody and are the primary biogeochemical rates causing oxygen depletion in coastal waters. However, given the additional labor involved in measuring biogeochemical rate processes, few monitoring programs regularly measure these properties and thus few long-term monitoring records of plankton respiration exist. An eight-year, biweekly plankton community respiration rate time series was analyzed as part of a monitoring program situated in the lower Patuxent River estuary, a tributary of Chesapeake Bay. We found that particulate nutrients (nitrogen and phosphorus) were the most highly correlated co-variates with respiration rate. Additionally, statistical and kinetic models including variables both water temperature and particulate nitrogen were able to explain 74% of the variability in respiration. Over the long-term record, both particulate nutrients and respiration rate were elevated when measured at higher tides. Separate measurements of respiration rate during ten consecutive days and during high and low tide on three separate days also support the enhancement of respiration with high tide. The enhancement was likely due to the import of particulate nutrients from the highly productive mid-bay region. This analysis of the longest consistently measured community respiration rate dataset in Chesapeake Bay has implications for how to interpret long-term records of measurements made at fixed locations in estuaries. 
    more » « less