Wetland carbon storage controlled by millennial-scale variation in relative sea-level rise.

Coastal wetlands (mangrove, tidal marsh and seagrass) sustain the highest rates of carbon sequestration per unit area of all natural systems1,2, primarily because of their comparatively high productivity and preservation of organic carbon within sedimentary substrates3. Climate change and associated relative sea-level rise (RSLR) have been proposed to increase the rate of organic-carbon burial in coastal wetlands in the first half of the twenty-first century4, but these carbon-climate feedback effects have been modelled to diminish over time as wetlands are increasingly submerged and carbon stores become compromised by erosion4,5. Here we show that tidal marshes on coastlines that experienced rapid RSLR over the past few millennia (in the late Holocene, from about 4,200 years ago to the present) have on average 1.7 to 3.7 times higher soil carbon concentrations within 20 centimetres of the surface than those subject to a long period of sea-level stability. This disparity increases with depth, with soil carbon concentrations reduced by a factor of 4.9 to 9.1 at depths of 50 to 100 centimetres. We analyse the response of a wetland exposed to recent rapid RSLR following subsidence associated with pillar collapse in an underlying mine and demonstrate that the gain in carbon accumulation and elevation is proportional to the accommodation space (that is, the space available for mineral and organic material accumulation) created by RSLR. Our results suggest that coastal wetlands characteristic of tectonically stable coastlines have lower carbon storage owing to a lack of accommodation space and that carbon sequestration increases according to the vertical and lateral accommodation space6 created by RSLR. Such wetlands will provide long-term mitigating feedback effects that are relevant to global climate-carbon modelling.

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

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.

Methane fluxes in tidal marshes of the conterminous United States.

Methane (CH4) is a potent greenhouse gas (GHG) with atmospheric concentrations that have nearly tripled since pre-industrial times. Wetlands account for a large share of global CH4 emissions, yet the magnitude and factors controlling CH4 fluxes in tidal wetlands remain uncertain. We synthesized CH4 flux data from 100 chamber and 9 eddy covariance (EC) sites across tidal marshes in the conterminous United States to assess controlling factors and improve predictions of CH4 emissions. This effort included creating an open-source database of chamber-based GHG fluxes (https://doi.org/10.25573/serc.14227085). Annual fluxes across chamber and EC sites averaged 26 ± 53 g CH4 m-2 year-1, with a median of 3.9 g CH4 m-2 year-1, and only 25% of sites exceeding 18 g CH4 m-2 year-1. The highest fluxes were observed at fresh-oligohaline sites with daily maximum temperature normals (MATmax) above 25.6°C. These were followed by frequently inundated low and mid-fresh-oligohaline marshes with MATmax ≤25.6°C, and mesohaline sites with MATmax >19°C. Quantile regressions of paired chamber CH4 flux and porewater biogeochemistry revealed that the 90th percentile of fluxes fell below 5 ± 3 nmol m-2 s-1 at sulfate concentrations >4.7 ± 0.6 mM, porewater salinity >21 ± 2 psu, or surface water salinity >15 ± 3 psu. Across sites, salinity was the dominant predictor of annual CH4 fluxes, while within sites, temperature, gross primary productivity (GPP), and tidal height controlled variability at diel and seasonal scales. At the diel scale, GPP preceded temperature in importance for predicting CH4 flux changes, while the opposite was observed at the seasonal scale. Water levels influenced the timing and pathway of diel CH4 fluxes, with pulsed releases of stored CH4 at low to rising tide. This study provides data and methods to improve tidal marsh CH4 emission estimates, support blue carbon assessments, and refine national and global GHG inventories.

Effects of spatial variability in vegetation phenology, climate, landcover, biodiversity, topography, and soil property on soil respiration across a coastal ecosystem.

Coastal terrestrial-aquatic interfaces (TAIs) are crucial contributors to global biogeochemical cycles and carbon exchange. The soil carbon dioxide (CO2) efflux in these transition zones is however poorly understood due to the high spatiotemporal dynamics of TAIs, as various sub-ecosystems in this region are compressed and expanded by complex influences of tides, changes in river levels, climate, and land use. We focus on the Chesapeake Bay region to (i) investigate the spatial heterogeneity of the coastal ecosystem and identify spatial zones with similar environmental characteristics based on the spatial data layers, including vegetation phenology, climate, landcover, diversity, topography, soil property, and relative tidal elevation; (ii) understand the primary driving factors affecting soil respiration within sub-ecosystems of the coastal ecosystem. Specifically, we employed hierarchical clustering analysis to identify spatial regions with distinct environmental characteristics, followed by the determination of main driving factors using Random Forest regression and SHapley Additive exPlanations. Maximum and minimum temperature are the main drivers common to all sub-ecosystems, while each region also has additional unique major drivers that differentiate them from one another. Precipitation exerts an influence on vegetated lands, while soil pH value holds importance specifically in forested lands. In croplands characterized by high clay content and low sand content, the significant role is attributed to bulk density. Wetlands demonstrate the importance of both elevation and sand content, with clay content being more relevant in non-inundated wetlands than in inundated wetlands. The topographic wetness index significantly contributes to the mixed vegetation areas, including shrub, grass, pasture, and forest. Additionally, our research reveals that dense vegetation land covers and urban/developed areas exhibit distinct soil property drivers. Overall, our research demonstrates an efficient method of employing various open-source remote sensing and GIS datasets to comprehend the spatial variability and soil respiration mechanisms in coastal TAI. There is no one-size-fits-all approach to modeling carbon fluxes released by soil respiration in coastal TAIs, and our study highlights the importance of further research and monitoring practices to improve our understanding of carbon dynamics and promote the sustainable management of coastal TAIs.

Controls on spatial variation in porewater methane concentrations across United States tidal wetlands.

Tidal wetlands can be a substantial sink of greenhouse gases, which can be offset by variable methane (CH4) emissions under certain environmental conditions and anthropogenic interventions. Land managers and policymakers need maps of tidal wetland CH4 properties to make restoration decisions and inventory greenhouse gases. However, there is a mismatch in spatial scale between point-based sampling of porewater CH4 concentration and its predictors, and the coarser resolution mapping products used to upscale these data. We sampled porewater CH4 concentrations, salinity, sulfate (SO42-), ammonium (NH4+), and total Fe using a spatially stratified sampling at 27 tidal wetlands in the United States. We measured porewater CH4 concentrations across four orders of magnitude (0.05 to 852.9 µM). The relative contribution of spatial scale to variance in CH4 was highest between- and within-sites. Porewater CH4 concentration was best explained by SO42- concentration with segmented linear regression (p < 0.01, R2 = 0.54) indicating lesser sensitivity of CH4 to SO42- below 0.62 mM SO42-. Salinity was a significant proxy for CH4 concentration, because it was highly correlated with SO42- (p < 0.01, R2 = 0.909). However, salinity was less predictive of CH4 with segmented linear regression (p < 0.01, R2 = 0.319) relative to SO42-. Neither NH4+, total Fe, nor relative tidal elevation correlated significantly with porewater CH4; however, NH4+ was positively and significantly correlated with SO42- after detrending CH4 for its relationship with SO42- (p < 0.01, R2 = 0.194). Future sampling should focus on within- and between-site environmental gradients to accurately map CH4 variation. Mapping salinity at sub-watershed scales has some potential for mapping SO42-, and by proxy, constraining spatial variation in porewater CH4 concentrations. Additional work is needed to explain site-level deviations from the salinity-sulfate relationship and elucidate other predictors of methanogenesis. This work demonstrates a unique approach to remote team science and the potential to strengthen collaborative research networks.

Global dataset of soil organic carbon in tidal marshes.

Tidal marshes store large amounts of organic carbon in their soils. Field data quantifying soil organic carbon (SOC) stocks provide an important resource for researchers, natural resource managers, and policy-makers working towards the protection, restoration, and valuation of these ecosystems. We collated a global dataset of tidal marsh soil organic carbon (MarSOC) from 99 studies that includes location, soil depth, site name, dry bulk density, SOC, and/or soil organic matter (SOM). The MarSOC dataset includes 17,454 data points from 2,329 unique locations, and 29 countries. We generated a general transfer function for the conversion of SOM to SOC. Using this data we estimated a median (± median absolute deviation) value of 79.2 ± 38.1 Mg SOC ha-1 in the top 30 cm and 231 ± 134 Mg SOC ha-1 in the top 1 m of tidal marsh soils globally. This data can serve as a basis for future work, and may contribute to incorporation of tidal marsh ecosystems into climate change mitigation and adaptation strategies and policies.

A Conterminous USA-Scale Map of Relative Tidal Marsh Elevation.

Tidal wetlands provide myriad ecosystem services across local to global scales. With their uncertain vulnerability or resilience to rising sea levels, there is a need for mapping flooding drivers and vulnerability proxies for these ecosystems at a national scale. However, tidal wetlands in the conterminous USA are diverse with differing elevation gradients, and tidal amplitudes, making broad geographic comparisons difficult. To address this, a national-scale map of relative tidal elevation (Z*MHW), a physical metric that normalizes elevation to tidal amplitude at mean high water (MHW), was constructed for the first time at 30 × 30-m resolution spanning the conterminous USA. Contrary to two study hypotheses, watershed-level median Z*MHW and its variability generally increased from north to south as a function of tidal amplitude and relative sea-level rise. These trends were also observed in a reanalysis of ground elevation data from the Pacific Coast by Janousek et al. (Estuaries and Coasts 42 (1): 85-98, 2019). Supporting a third hypothesis, propagated uncertainty in Z*MHW increased from north to south as light detection and ranging (LiDAR) errors had an outsized effect under narrowing tidal amplitudes. The drivers of Z*MHW and its variability are difficult to determine because several potential causal variables are correlated with latitude, but future studies could investigate highest astronomical tide and diurnal high tide inequality as drivers of median Z*MHW and Z*MHW variability, respectively. Watersheds of the Gulf Coast often had propagated Z*MHW uncertainty greater than the tidal amplitude itself emphasizing the diminished practicality of applying Z*MHW as a flooding proxy to microtidal wetlands. Future studies could focus on validating and improving these physical map products and using them for synoptic modeling of tidal wetland carbon dynamics and sea-level rise vulnerability analyses. Supplementary Information: The online version contains supplementary material available at 10.1007/s12237-021-01027-9.

Author Correction: Accuracy and Precision of Tidal Wetland Soil Carbon Mapping in the Conterminous United States.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has not been fixed in the paper.

Accuracy and Precision of Tidal Wetland Soil Carbon Mapping in the Conterminous United States.

Tidal wetlands produce long-term soil organic carbon (C) stocks. Thus for carbon accounting purposes, we need accurate and precise information on the magnitude and spatial distribution of those stocks. We assembled and analyzed an unprecedented soil core dataset, and tested three strategies for mapping carbon stocks: applying the average value from the synthesis to mapped tidal wetlands, applying models fit using empirical data and applied using soil, vegetation and salinity maps, and relying on independently generated soil carbon maps. Soil carbon stocks were far lower on average and varied less spatially and with depth than stocks calculated from available soils maps. Further, variation in carbon density was not well-predicted based on climate, salinity, vegetation, or soil classes. Instead, the assembled dataset showed that carbon density across the conterminous united states (CONUS) was normally distributed, with a predictable range of observations. We identified the simplest strategy, applying mean carbon density (27.0 kg C m-3), as the best performing strategy, and conservatively estimated that the top meter of CONUS tidal wetland soil contains 0.72 petagrams C. This strategy could provide standardization in CONUS tidal carbon accounting until such a time as modeling and mapping advancements can quantitatively improve accuracy and precision.

Prolonged California aridity linked to climate warming and Pacific sea surface temperature.

California has experienced a dry 21(st) century capped by severe drought from 2012 through 2015 prompting questions about hydroclimatic sensitivity to anthropogenic climate change and implications for the future. We address these questions using a Holocene lake sediment record of hydrologic change from the Sierra Nevada Mountains coupled with marine sediment records from the Pacific. These data provide evidence of a persistent relationship between past climate warming, Pacific sea surface temperature (SST) shifts and centennial to millennial episodes of California aridity. The link is most evident during the thermal-maximum of the mid-Holocene (~8 to 3 ka; ka = 1,000 calendar years before present) and during the Medieval Climate Anomaly (MCA) (~1 ka to 0.7 ka). In both cases, climate warming corresponded with cooling of the eastern tropical Pacific despite differences in the factors producing increased radiative forcing. The magnitude of prolonged eastern Pacific cooling was modest, similar to observed La Niña excursions of 1(o) to 2 °C. Given differences with current radiative forcing it remains uncertain if the Pacific will react in a similar manner in the 21st century, but should it follow apparent past behavior more intense and prolonged aridity in California would result.