USGS Science⌕ Search

SEARCH · USGS Science

Results for “Carbon”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

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

At least 757 records · Page 42Linked to original sources

Linking soil organic matter dynamics and erosion-induced terrestrial carbon sequestration at different landform positions

Recently, the potential for terrestrial carbon (C) sequestration by soil erosion and deposition has received increased interest. Erosion and deposition constitute a sink for atmospheric carbon dioxide relative to a preerosional state or a noneroding scenario, if the posterosion watershed C balance is increased due to (1) partial replacement of eroded C by new photosynthate in the eroded site; and (2) preservation from decomposition of at least some eroded soil organic carbon (SOC) arriving in depositional settings. Little is known, however, about differences in C dynamics at different erosional and depositional landform positions within the same eroding system. We determined the contribution of different landform positions to erosion-induced terrestrial C sequestration by measuring rates of net primary productivity (NPP), replacement of eroded C, and decomposition of organic matter (OM) at four categorically different landform positions within a naturally eroding toposequence in northern California. We found that eroded C is replaced by NPP 15 times over in the summit of the site studied and 5 times over in the slope. Profile-averaged, long-term rate constant for SOM decomposition was 2 to 14 times slower in the depositional settings compared with that in eroding slopes. As a result, the inventory of C in the depositional settings was 2 to 3 times larger than that of the eroding positions. Owing to both C replacement at eroding sites and reduced rates of OM decomposition in depositional sites, soil erosion constitutes a C sink from the atmosphere at our study site.

California↗

Controls on mangrove forest‐atmosphere carbon dioxide exchanges in western Everglades National Park

We report on net ecosystem production (NEP) and key environmental controls on net ecosystem exchange (NEE) of carbon dioxide (CO 2 ) between a mangrove forest and the atmosphere in the coastal Florida Everglades. An eddy covariance system deployed above the canopy was used to determine NEE during January 2004 through August 2005. Maximum daytime NEE ranged from −20 to −25 μ mol (CO 2 ) m −2 s −1 between March and May. Respiration (R d ) was highly variable (2.81 ± 2.41 μ mol (CO 2 ) m −2 s −1 ), reaching peak values during the summer wet season. During the winter dry season, forest CO 2 assimilation increased with the proportion of diffuse solar irradiance in response to greater radiative transfer in the forest canopy. Surface water salinity and tidal activity were also important controls on NEE. Daily light use efficiency was reduced at high (>34 parts per thousand (ppt)) compared to low (<17 ppt) salinity by 46%. Tidal inundation lowered daytime R d by ∼0.9 μ mol (CO 2 ) m −2 s −1 and nighttime R d by ∼0.5 μ mol (CO 2 ) m −2 s −1 . The forest was a sink for atmospheric CO 2 , with an annual NEP of 1170 ± 127 g C m −2 during 2004. This unusually high NEP was attributed to year‐round productivity and low ecosystem respiration which reached a maximum of only 3 g C m −2 d −1 . Tidal export of dissolved inorganic carbon derived from belowground respiration likely lowered the estimates of mangrove forest respiration. These results suggest that carbon balance in mangrove coastal systems will change in response to variable salinity and inundation patterns, possibly resulting from secular sea level rise and climate change.

Florida↗

A model‐data intercomparison of CO2 exchange across North America: Results from the North American Carbon Program site synthesis

Our current understanding of terrestrial carbon processes is represented in various models used to integrate and scale measurements of CO 2 exchange from remote sensing and other spatiotemporal data. Yet assessments are rarely conducted to determine how well models simulate carbon processes across vegetation types and environmental conditions. Using standardized data from the North American Carbon Program we compare observed and simulated monthly CO 2 exchange from 44 eddy covariance flux towers in North America and 22 terrestrial biosphere models. The analysis period spans ∼220 site‐years, 10 biomes, and includes two large‐scale drought events, providing a natural experiment to evaluate model skill as a function of drought and seasonality. We evaluate models' ability to simulate the seasonal cycle of CO 2 exchange using multiple model skill metrics and analyze links between model characteristics, site history, and model skill. Overall model performance was poor; the difference between observations and simulations was ∼10 times observational uncertainty, with forested ecosystems better predicted than nonforested. Model‐data agreement was highest in summer and in temperate evergreen forests. In contrast, model performance declined in spring and fall, especially in ecosystems with large deciduous components, and in dry periods during the growing season. Models used across multiple biomes and sites, the mean model ensemble, and a model using assimilated parameter values showed high consistency with observations. Models with the highest skill across all biomes all used prescribed canopy phenology, calculated NEE as the difference between GPP and ecosystem respiration, and did not use a daily time step.

Journal of Geophysical Research: Biogeosciences↗

Characterizing the performance of ecosystem models across time scales: A spectral analysis of the North American Carbon Program site‐level synthesis

[1] Ecosystem models are important tools for diagnosing the carbon cycle and projecting its behavior across space and time. Despite the fact that ecosystems respond to drivers at multiple time scales, most assessments of model performance do not discriminate different time scales. Spectral methods, such as wavelet analyses, present an alternative approach that enables the identification of the dominant time scales contributing to model performance in the frequency domain. In this study we used wavelet analyses to synthesize the performance of 21 ecosystem models at 9 eddy covariance towers as part of the North American Carbon Program's site‐level intercomparison. This study expands upon previous single‐site and single‐model analyses to determine what patterns of model error are consistent across a diverse range of models and sites. To assess the significance of model error at different time scales, a novel Monte Carlo approach was developed to incorporate flux observation error. Failing to account for observation error leads to a misidentification of the time scales that dominate model error. These analyses show that model error (1) is largest at the annual and 20–120 day scales, (2) has a clear peak at the diurnal scale, and (3) shows large variability among models in the 2–20 day scales. Errors at the annual scale were consistent across time, diurnal errors were predominantly during the growing season, and intermediate‐scale errors were largely event driven. Breaking spectra into discrete temporal bands revealed a significant model‐by‐band effect but also a nonsignificant model‐by‐site effect, which together suggest that individual models show consistency in their error patterns. Differences among models were related to model time step, soil hydrology, and the representation of photosynthesis and phenology but not the soil carbon or nitrogen cycles. These factors had the greatest impact on diurnal errors, were less important at annual scales, and had the least impact at intermediate time scales.

Journal of Geophysical Research: Biogeosciences↗

Carbon dioxide and methane emissions from the Yukon River system

Carbon dioxide (CO 2 ) and methane (CH 4 ) emissions are important, but poorly quantified, components of riverine carbon (C) budgets. This is largely because the data needed for gas flux calculations are sparse and are spatially and temporally variable. Additionally, the importance of C gas emissions relative to lateral C exports is not well known because gaseous and aqueous fluxes are not commonly measured on the same rivers. We couple measurements of aqueous CO 2 and CH 4 partial pressures ( p CO 2 , p CH 4 ) and flux across the water-air interface with gas transfer models to calculate subbasin distributions of gas flux density. We then combine those flux densities with remote and direct observations of stream and river water surface area and ice duration, to calculate C gas emissions from flowing waters throughout the Yukon River basin. CO 2 emissions were 7.68 Tg C yr −1 (95% CI: 5.84 −10.46), averaging 750 g C m −2 yr −1 normalized to water surface area, and 9.0 g C m −2 yr −1 normalized to river basin area. River CH 4 emissions totaled 55 Gg C yr −1 or 0.7% of the total mass of C emitted as CO 2 plus CH 4 and ∼6.4% of their combined radiative forcing. When combined with lateral inorganic plus organic C exports to below head of tide, C gas emissions comprised 50% of total C exported by the Yukon River and its tributaries. River CO 2 and CH 4 derive from multiple sources, including groundwater, surface water runoff, carbonate equilibrium reactions, and benthic and water column microbial processing of organic C. The exact role of each of these processes is not yet quantified in the overall river C budget.

Yukon River Basin↗

Enhanced transfer of terrestrially derived carbon to the atmosphere in a flooding event

Rising CO 2 concentration in the atmosphere, global climate change, and the sustainability of the Earth's biosphere are great societal concerns for the 21st century. Global climate change has, in part, resulted in a higher frequency of flooding events, which allow for greater exchange between soil/plant litter and aquatic carbon pools. Here we demonstrate that the summer 2011 flood in the Mississippi River basin, caused by extreme precipitation events, resulted in a “flushing” of terrestrially derived dissolved organic carbon (TDOC) to the northern Gulf of Mexico. Data from the lower Atchafalaya and Mississippi rivers showed that the DOC flux to the northern Gulf of Mexico during this flood was significantly higher than in previous years. We also show that consumption of radiocarbon-modern TDOC by bacteria in floodwaters in the lower Atchafalaya River and along the adjacent shelf contributed to northern Gulf shelf waters changing from a net sink to a net source of CO 2 to the atmosphere in June and August 2011. This work shows that enhanced flooding, which may or may not be caused by climate change, can result in rapid losses of stored carbon in soils to the atmosphere via processes in aquatic ecosystems.

Geophysical Research Letters↗

Mechanisms of earthquake‐induced chemical and fluid transport to carbonate groundwater springs after earthquakes

Mechanisms by which hydrochemical changes occur after earthquakes are not well documented. We use the 2016–2017 central Italy seismic sequence, which caused notable hydrochemical transient variations in groundwater springs to address this topic, with special reference to effects on fractured carbonate aquifers. Hydrochemistry measured before and after the earthquakes at four springs at varying distances from the epicenters all showed immediate postmainshock peaks in trace element concentrations but little change in major elements. Most parameters returned to preearthquake values before the last events of the seismic sequence. The source of solutes, particularly trace elements, is longer residence time pore water stored in slow‐moving fractures or abandoned karstic flow paths. These fluids were expelled into the main flow paths after an increase in pore pressure, hydraulic conductivity, and shaking from coseismic aquifer stress. The weak response to the later earthquakes is explained by progressive depletion of high solute fluids as earlier shocks flushed out the stored fluids in the fractures. Spring δ 13 C DIC values closest to a deep magma source to the west became enriched relative to preearthquake values following the 24 August event. This enrichment indicates input from deeply sourced dissolved CO 2 gas after dilation of specific fault conduits. Differences in carbon isotopic responses between springs are attributed to proximity to the deep gaseous CO 2 source. Most of the transient chemical changes seen in the three fractured carbonate aquifers are attributed to local shaking and emptying of isolated pores and fractures and are not from rapid upward movement of deep fluids.

Water Resources Research↗

The role of the upper tidal estuary in wetland blue carbon storage and flux

Carbon (C) standing stocks, C mass balance, and soil C burial in tidal freshwater forested wetlands (TFFW) and TFFW transitioning to low‐salinity marshes along the upper estuary are not typically included in “blue carbon” accounting, but may represent a significant C sink. Results from two salinity transects along the tidal Waccamaw and Savannah rivers of the US Atlantic Coast show total C standing stocks were 321‐1264 Mg C ha ‐1 among all sites, generally shifting to greater soil storage as salinity increased. Carbon mass balance inputs (litterfall, woody growth, herbaceous growth, root growth, surface accumulation) minus C outputs (surface litter and root decomposition, gaseous C) over a period of up to 11 years were 340‐900 g C m ‐2 yr ‐1 . Soil C burial was variable (7‐337 g C m ‐2 yr ‐1 ), and lateral C export was estimated as C mass balance minus soil C burial as 267‐849 g C m ‐2 yr ‐1 . This represents a large amount of C export to support aquatic biogeochemical transformations. Despite reduced C persistence within emergent vegetation, decomposition of organic matter, and higher lateral C export, total C storage increased as forests converted to marsh with salinization. These tidal river wetlands exhibited high N mineralization in salinity‐stressed forested sites and considerable P mineralization in low salinity marshes. Large C standing stocks and rates of C sequestration suggest that TFFW and oligohaline marshes are considerably important globally to coastal C dynamics and in facilitating energy transformations in areas of the world in which they occur.

Georgia, South Carolina↗

Carbon chemistry of intact versus chronically drained peatlands in the southeastern USA

The Great Dismal Swamp (GDS) is a large temperate swamp in Virginia/North Carolina with peat soils historically resistant to microbial decomposition. However, this peatland has been subject to ~200 years of disturbance during which extensive drainage, fire suppression, and wide-spread logging have increased decomposition and dramatically decreased the distribution of Atlantic white cedar (AWC). The purpose of this study was to determine the impact of long-term drainage and AWC loss on the carbon chemistry of GDS peats. Peat cores were collected from three drained GDS vegetation communities (pocosin, AWC, and red maple-black gum) and compared to cores collected from an intact, undrained AWC peatland at the Alligator River National Wildlife Refuge (AR) in North Carolina, USA. The AR peats had higher lignin content in the deeper peat intervals, and lignin content and % organic carbon were largely invariant with depth compared to the GDS peats. The concentrations of syringyl group phenols were greater in the surface layers of GDS peats, likely reflecting the selective removal of AWC and transition from gymnosperms to angiosperms. Acid to aldehyde ratios for vanillyl and syringyl group phenols indicated that the GDS peats were more decomposed, particularly at depth, and that this occurred under aerobic conditions. Moreover, solid state 13C NMR confirmed a coincident loss of carbohydrates and increase in recalcitrant byproducts of carbohydrate degradation with depth. These data indicate that long-term drainage has accelerated the decomposition of peat at the GDS, reducing the capacity and stability of the carbon sink.

Virginia, North Carolina↗

Hydrothermal manganese-oxide mineralization of a carbonate ooze, Samoan hotspot region, South Pacific Ocean

Low-temperature hydrothermal manganese oxides occur throughout the global oceans. However, the hydrothermal replacement of a carbonate ooze by manganese oxides is described here for the first time. The 24 samples dredged from three locations around the Territory of American Samoa in the South Pacific Ocean include six samples with remnant carbonate and volcaniclastic sediments, which we refer to as “low Mn,” and 18 samples in which the mineralization is pervasive and has replaced most original sediment, termed “high Mn.” The 18 high-Mn samples exhibit a mean Mn content of 51 wt.%. Higher Li and Mn contents and lower Fe contents in the high-Mn samples indicate a hydrothermal origin and distinguish these samples from hydrogenetic ferromanganese crusts. Mn-oxide layers are up to 90 mm thick, some with columns to 44 mm long and 10 mm wide, the magnitude of which has not been described previously. Samples are composed of birnessite and 10 Å phyllomanganate minerals. Textures of the thickest Mn layers indicate mineralization below the seabed from ascending fluids during multiple hydrothermal pulses. Mineralization took place by complete to partial replacement and cementation of foraminiferal sediments intermixed with volcaniclastic sediments in varying amounts. Our results highlight the production of carbonate sediment-hosted hydrothermal Mn oxides from multiple hydrothermal sources within the Samoan volcanic chain. The potential extensive distribution on the regional scale of this newly described mineralization process and unique element enrichments raise questions about its broader distribution globally and potential importance to hydrothermal processes, element mass balance, and seabed mineral resources.

Geochemistry, Geophysics, Geosystems↗

On the influence of biomass burning on the seasonal CO2 signal as observed at monitoring stations

We investigated the role of biomass burning in simulating the seasonal signal in both prognostic and diagnostic analyses. The prognostic analysis involved the High-Resolution Biosphere Model, a prognostic terrestrial biosphere model, and the coupled vegetation fire module, which together produce a prognostic data set of biomass burning. The diagnostic analysis involved the Simple Diagnostic Biosphere Model (SDBM) and the Hao and Liu [1994] diagnostic data set of biomass burning, which have been scaled to global 2 and 4 Pg C yr −1 , respectively. The monthly carbon exchange fields between the atmosphere and the biosphere with a spatial resolution of 0.5° × 0.5°, the seasonal atmosphere-ocean exchange fields, and the emissions from fossil fuels have been coupled to the three-dimensional atmospheric transport model TM2. We have chosen eight monitoring stations of the National Oceanic and Atmospheric Administration network to compare the predicted seasonal atmospheric CO 2 signals with those deduced from atmosphere-biosphere carbon exchange fluxes without any contribution from biomass burning. The prognostic analysis and the diagnostic analysis with global burning emissions of 4 Pg C yr −1 agree with respect to the change in the amplitude of the seasonal CO 2 concentration introduced through biomass burning. We find that the seasonal CO 2 signal at stations in higher northern latitudes (north of 30°N) is marginally influenced by biomass burning. For stations in tropical regions an increase in the CO 2 amplitude of more than 1 ppmv (up to 50% with respect to the observed trough to peak amplitude) has been calculated. Biomass burning at stations farther south accounts for an increase in the CO 2 amplitude of up to 59% (0.6 ppmv). A change in the phase of the seasonal CO 2 signal at tropical and southern stations has been shown to be strongly influenced by the onset of biomass burning in southern tropical Africa and America. Comparing simulated and observed seasonal CO 2 signals, we find higher discrepancies at southern tropical stations if biomass burning emissions are included. This is caused by the additional increase in the amplitude in the prognostic analysis and a phase shift in a diagnostic analysis. In contrast, at the northern tropical stations biomass burning tends to improve the estimates of the seasonal CO 2 signal in the prognostic analysis because of strengthening of the amplitude. Since the SDBM predicts the seasonal CO 2 signal reasonably well for the northern hemisphere tropical stations, no general improvement of the fit occurs if biomass burning emissions are considered.

Global Biogeochemical Cycles↗

Rhythmic bedding produced in Cretaceous pelagic carbonate environments: Sensitive recorders of climatic cycles

Various types of rhythmic bedding are prominent features of Cretaceous pelagic carbonates. These bedding rhythms are the primary depositional result of variations in carbonate productivity, terrigenous dilution, redox conditions and/or the energy of bottom currents. Each bedding cycle type is different in its expression in the stratigraphic record but ultimately was caused by rhythmic climatic variations related to variations in receipt of solar insolation as the result of the earth's orbital characteristics, the so‐called Milankovitch cycles. Stratigraphic sequences may contain a single bedding cycle type or a composite of several types. Although sediment type and sedimentary structures were influenced by variations in local climatic, oceanographic, tectonic, and diagenetic parameters, pelagic carbonate environments during the Cretaceous were particularly sensitive recorders of orbitally induced changes in worldwide climate.

Paleoceanography↗

Origin and distribution of carbon dioxide in the unsaturated zone of the southern High Plains of Texas

Partial pressures of CO 2 , O 2 , N 2 , and Ar were monitored at two locations in the Ogallala aquifer system on the Southern High Plains of Texas. Samples were collected monthly during parts of 1980–1981 from nine depths ranging from 0.6 to 36 meters below land surface. P CO 2 was observed to be greater at depth than in the active soil zone and thus appears to contradict the normal process in which CO 2 is generated in the soil zone and diffuses upward to the atmosphere and downward to the water table. The δ 13 C of the CO 2 gas was quite uniform and averaged −17.9 per mil. P O 2 declined with depth, suggesting in situ generation of CO 2 by the oxidation of carbon. Several hypotheses were considered to explain the origin of the CO 2 at depth. It was concluded that the most probable hypothesis was that dissolved and particulate organic carbon introduced by recharging water was oxidized to CO 2 by the aerobic microbial community that utilized oxygen diffusing in from the atmosphere. This hypothesis is consistent with the CO 2 concentration profile, calculated production profile of CO 2 , δ 13 C values of CO 2 gas, caliche, soil humic acid fraction, and dissolved carbonate in groundwater. The abundance of CO 2 , its concentration profile, and its probable origin provide information for evaluating the observed complex sequence of caliche dissolution and precipitation known to occur in the aquifer.

Texas↗

Liquid carbon dioxide of magmatic origin and its role in volcanic eruptions

Natural liquid carbon dioxide is produced commercially from a 2.5-km-deep well near the 4,500-yr-old maar volcano, Mount Gambier, South Australia. The carbon dioxide has accumulated in a dome that is located on the extension of a linear chain of volcanic activity. A magmatic origin for the fluid is suggested by the geological setting, δ 13 C PDB of –4.0‰, for the CO 2 (where PDB represents the carbon-isotope standard), and a relatively high 3 He component of the contained helium and high 3 He/C ratio (6.4 x 10 −10 ). The 3 He/ 4 He and He/Ne ratios are 3.0 and > 1,370 times those of air, respectively. The CO 2 , as collected at the Earth's surface at 29.5 °C and 75 bar, expands more than 300-fold to form a gas at 1 atm and 22 °C. We suggest that liquid CO 2 or high-density CO 2 fluid (the critical point is 31.1 °C, 73.9 bar) of volcanic origin that expands explosively from shallow levels in the Earth's crust may be a major contributor to 'phreatic' volcanic eruptions and maar formation. Less violent release of magmatic CO 2 into crater lakes may cause gas bursts with equally disastrous consequences such as occurred at Lake Nyos, Cameroon, in August 1986.

Nature↗

Permafrost collapse is accelerating carbon release

This much is clear: the Arctic is warming fast, and frozen soils are starting to thaw, often for the first time in thousands of years. But how this happens is as murky as the mud that oozes from permafrost when ice melts. As the temperature of the ground rises above freezing, microorganisms break down organic matter in the soil. Greenhouse gases — including carbon dioxide, methane and nitrous oxide — are released into the atmosphere, accelerating global warming. Soils in the permafrost region hold twice as much carbon as the atmosphere does — almost 1,600 billion tonnes 1 . What fraction of that will decompose? Will it be released suddenly, or seep out slowly? We need to find out. Current models of greenhouse-gas release and climate assume that permafrost thaws gradually from the surface downwards. Deeper layers of organic matter are exposed over decades or even centuries, and some models are beginning to track these slow changes. But models are ignoring an even more troubling problem. Frozen soil doesn’t just lock up carbon — it physically holds the landscape together. Across the Arctic and Boreal regions, permafrost is collapsing suddenly as pockets of ice within it melt. Instead of a few centimetres of soil thawing each year, several metres of soil can become destabilized within days or weeks. The land can sink and be inundated by swelling lakes and wetlands. Abrupt thawing of permafrost is dramatic to watch. Returning to field sites in Alaska, for example, we often find that lands that were forested a year ago are now covered with lakes 2 . Rivers that once ran clear are thick with sediment. Hillsides can liquefy, sometimes taking sensitive scientific equipment with them. This type of thawing is a serious problem for communities living around the Arctic (see ‘Arctic permafrost’). Roads buckle, houses become unstable. Access to traditional foods is changing, because it is becoming dangerous to travel across the land to hunt. Families cannot reach lines of game traps that have supported them for generations.

Nature↗

Global carbon dioxide emissions from inland waters

Carbon dioxide (CO 2 ) transfer from inland waters to the atmosphere, known as CO 2 evasion, is a component of the global carbon cycle. Global estimates of CO 2 evasion have been hampered, however, by the lack of a framework for estimating the inland water surface area and gas transfer velocity and by the absence of a global CO 2 database. Here we report regional variations in global inland water surface area, dissolved CO 2 and gas transfer velocity. We obtain global CO 2 evasion rates of 1.8 &thinsp; petagrams of carbon (Pg &thinsp; C) per year from streams and rivers and 0.32 &thinsp; Pg &thinsp; C &thinsp; yr &minus;1 from lakes and reservoirs, where the upper and lower limits are respectively the 5th and 95th confidence interval percentiles. The resulting global evasion rate of 2.1 &thinsp; Pg &thinsp; C &thinsp; yr &minus;1 is higher than previous estimates owing to a larger stream and river evasion rate. Our analysis predicts global hotspots in stream and river evasion, with about 70 per cent of the flux occurring over just 20 per cent of the land surface. The source of inland water CO 2 is still not known with certainty and new studies are needed to research the mechanisms controlling CO 2 evasion globally.

Nature↗

Differentiating moss from higher plants is critical in studying the carbon cycle of the boreal biome

The satellite-derived normalized difference vegetation index (NDVI), which is used for estimating gross primary production (GPP), often includes contributions from both mosses and vascular plants in boreal ecosystems. For the same NDVI, moss can generate only about one-third of the GPP that vascular plants can because of its much lower photosynthetic capacity. Here, based on eddy covariance measurements, we show that the difference in photosynthetic capacity between these two plant functional types has never been explicitly included when estimating regional GPP in the boreal region, resulting in a substantial overestimation. The magnitude of this overestimation could have important implications regarding a change from a current carbon sink to a carbon source in the boreal region. Moss abundance, associated with ecosystem disturbances, needs to be mapped and incorporated into GPP estimates in order to adequately assess the role of the boreal region in the global carbon cycle.

Nature Communications↗

Carbon isotope equilibration during sulphate-limited anaerobic oxidation of methane

Collectively, marine sediments comprise the largest reservoir of methane on Earth. The flux of methane from the sea bed to the overlying water column is mitigated by the sulphate-dependent anaerobic oxidation of methane by marine microbes within a discrete sedimentary horizon termed the sulphate–methane transition zone. According to conventional isotope systematics, the biological consumption of methane leaves a residue of methane enriched in 13 C (refs 1–3). However, in many instances the methane within sulphate–methane transition zones is depleted in 13 C, consistent with the production of methane, and interpreted as evidence for the intertwined anaerobic oxidation and production of methane 4–6 . Here, we report results from experiments in which we incubated cultures of microbial methane consumers with methane and low levels of sulphate, and monitored the stable isotope composition of the methane and dissolved inorganic carbon pools over time. Residual methane became progressively enriched in 13 C at sulphate concentrations above 0.5 mM, and progressively depleted in 13 C below this threshold. We attribute the shift to 13 C depletion during the anaerobic oxidation of methane at low sulphate concentrations to the microbially mediated carbon isotope equilibration between methane and carbon dioxide. We suggest that this isotopic e ect could help to explain the 13 C-depletion of methane in subseafloor sulphate–methane transition zones.

Nature Geoscience↗