Pore water and sediment geochemistry of the Bothian Sea

Methane is a powerful greenhouse gas and its biological conversion in marine sediments, largely controlled by anaerobic oxidation of methane (AOM), is a crucial part of the global carbon cycle. However, little is known about the role of iron oxides as an oxidant for AOM. Here we provide the first field evidence for iron-dependent AOM in brackish coastal surface sediments and show that methane produced in Bothnian Sea sediments is oxidized in distinct zones of iron- and sulfate-dependent AOM. At our study site, anthropogenic eutrophication over recent decades has led to an upward migration of the sulfate/methane transition zone in the sediment. Abundant iron oxides and high dissolved ferrous iron indicate iron reduction in the methanogenic sediments below the newly established sulfate/methane transition. Laboratory incubation studies of these sediments strongly suggest that the in situ microbial community is capable of linking methane oxidation to iron oxide reduction. Eutrophication of coastal environments may therefore create geochemical conditions favorable for iron-mediated AOM and thus increase the relevance of iron-dependent methane oxidation in the future. Besides its role in mitigating methane emissions, iron-dependent AOM strongly impacts sedimentary iron cycling and related biogeochemical processes through the reduction of large quantities of iron oxides.

Gas composition, concentration and stable isotope ratios (Table 2, 3, 4)

The deployment of CCS (carbon capture and storage) at industrial scale implies the development of effective monitoring tools. Noble gases are tracers usually proposed to track CO2. This methodology, combined with the geochemistry of carbon isotopes, has been tested on available analogues. At first, gases from natural analogues were sampled in the Colorado Plateau and in the French carbogaseous provinces, in both well-confined and leaking-sites. Second, we performed a 2-years tracing experience on an underground natural gas storage, sampling gas each month during injection and withdrawal periods. In natural analogues, the geochemical fingerprints are dependent on the containment criterion and on the geological context, giving tools to detect a leakage of deep-CO2 toward surface. This study also provides information on the origin of CO2, as well as residence time of fluids within the crust and clues on the physico-chemical processes occurring during the geological story. The study on the industrial analogue demonstrates the feasibility of using noble gases as tracers of CO2. Withdrawn gases follow geochemical trends coherent with mixing processes between injected gas end-members. Physico-chemical processes revealed by the tracing occur at transient state. These two complementary studies proved the interest of geochemical monitoring to survey the CO2 behaviour, and gave information on its use.

Fig. S1. Loci of possible event locations based on differential ts - tp arrival times at the single functioning geophone

Geological storage of CO2 that has been captured at large, point source emitters represents a key potential method for reduction of anthropogenic greenhouse gas emissions. However, this technology will only be viable if it can be guaranteed that injected CO2 will remain trapped in the subsurface for thousands of years or more. A signi?cant issue for storage security is the geomechanical response of the reservoir. Concerns have been raised that geomechanical deformation induced by CO2 injection will create or reactivate fracture networks in the sealing caprocks, providing a pathway for CO2 leakage. In this paper, we examine three large-scale sites where CO2 is injected at rates of ab. 1 megatonne/y or more: Sleipner, Weyburn, and In Salah. We compare and contrast the observed geomechanical behavior of each site, with particular focus on the risks to storage security posed by geomechanical deformation. At Sleipner, the large, high-permeability storage aquifer has experienced little pore pressure increase over 15 y of injection, implying little possibility of geomechanical deformation. At Weyburn, 45 y of oil production has depleted pore pressures before increases associated with CO2 injection. The long history of the ?eld has led to complicated, sometimes nonintuitive geomechanical deformation. At In Salah, injection into the water leg of a gas reservoir has increased pore pressures, leading to uplift and substantial microseismic activity. The differences in the geomechanical responses of these sites emphasize the need for systematic geomechanical appraisal before injection in any potential storage site.

Geochemistry of sediment cores GeoB13820-1 and GeoB13863-1 from the western South Atlantic

Here, we present results from sediments collected in the Argentine Basin, a non-steady state depositional marine system characterized by abundant oxidized iron within methane-rich layers due to sediment reworking followed by rapid deposition. Our comprehensive inorganic data set shows that iron reduction in these sulfate and sulfide-depleted sediments is best explained by a microbially mediated process-implicating anaerobic oxidation of methane coupled to iron reduction (Fe-AOM) as the most likely major mechanism. Although important in many modern marine environments, iron-driven AOM may not consume similar amounts of methane compared with sulfate-dependent AOM. Nevertheless, it may have broad impact on the deep biosphere and dominate both iron and methane cycling in sulfate-lean marine settings. Fe-AOM might have been particularly relevant in the Archean ocean, >2.5 billion years ago, known for its production and accumulation of iron oxides (in iron formations) in a biosphere likely replete with methane but low in sulfate. Methane at that time was a critical greenhouse gas capable of sustaining a habitable climate under relatively low solar luminosity, and relationships to iron cycling may have impacted if not dominated methane loss from the biosphere.

Development of bacterial and archaeal communities in erupted subsurface muds at the Håkon Mosby mud volcano

The microbial oxidation of methane controls the emission of the greenhouse gas methane from the ocean floor. However, some seabed structures such as mud volcanoes have leaky microbial methane filters and can be important sources of methane. We investigated the disturbance and recovery of a methanotrophic mud volcano microbiome (Håkon Mosby mud volcano, 1250 m water depth), to assess time scales of community succession and function in the natural deep-sea environment. We analyzed 10 surface and 5 subsurface sediment samples across HMMV mud flows from most recently discharged subsurface muds towards old consolidated muds as well as one reference site (REF) located approximately 0.5 km outside of the HMMV. Surface samples were obtained in 2003, 2009 and 2010. The surface of the new mud flows at the geographical center was sampled in 2009 and 2010. Around 100 m south of the center, we sampled more consolidated aged muds in 2003 and 2010. Old mud flows were sampled around 300 m southeast and 100 m north of the geographical center in 2003, 2009 and 2010. Surface sediment samples (0-20 cm) were recovered either by TV-guided Multicorer or by push cores using the remotely operated vehicle Quest (Marum, University Bremen). Subsurface sediments of all zones (>2 m below sea floor) were obtained in 2003 by gravity corer. After recovery, sediments were immediately subsampled in a refrigerated container (0°C) and further processed for biogeochemical analyses or preserved at -20°C for later DNA analyses. Our study show that freshly erupted muds hosted heterotrophic deep subsurface communities, which were replaced by surface communities within a few years of exposure. Aerobic methanotrophy was established at the top surface layer within less than a year, followed by anaerobic methanotrophy, sulfate reduction and finally thiotrophy. Our data indicate that it takes decades in cold environments before efficient methanotrophic communities establish to control methane emission. The observed succession provides insights to the response time of complex deep-sea communities to seafloor disturbances.

Relative abundance of anaerobic methanotrophic archaea in bacterial mat Cascadia margin off Oregon

The microbially mediated anaerobic oxidation of methane (AOM) is the major biological sink of the greenhouse gas methane in marine sediments (doi:10.1007/978-94-009-0213-8_44) and serves as an important control for emission of methane into the hydrosphere. The AOM metabolic process is assumed to be a reversal of methanogenesis coupled to the reduction of sulfate to sulfide involving methanotrophic archaea (ANME) and sulfate-reducing bacteria (SRB) as syntrophic partners which were describes amongst others in Boetius et al. (2000; doi:10.1038/35036572). In this study, 16S rRNA-based methods were used to investigate the distribution and biomass of archaea in samples from sediments above outcropping methane hydrate at Hydrate Ridge (Cascadia margin off Oregon) and (ii) massive microbial mats enclosing carbonate reefs (Crimea area, Black Sea). Sediment samples from Hydrate Ridge were obtained during R/V SONNE cruises SO143-2 in August 1999 and SO148-1 in August 2000 at the crest of southern Hydrate Ridge at the Cascadia convergent margin off the coast of Oregon. The second study area is located in the Black Sea and represents a field in which there is active seepage of free gas on the slope of the northwestern Crimea area. Here, a field of conspicuous microbial reefs forming chimney-like structures was discovered at a water depth of 230 m in anoxic waters. The microbial mats were sampled by using the manned submersible JAGO during the R/V Prof. LOGACHEV cruise in July 2001. At Hydrate Ridge the surface sediments were dominated by aggregates consisting of ANME-2 and members of the Desulfosarcina-Desulfococcus branch (DSS) (ANME-2/DSS aggregates), which accounted for >90% of the total cell biomass. The numbers of ANME-1 cells increased strongly with depth; these cells accounted 1% of all single cells at the surface and more than 30% of all single cells (5% of the total cells) in 7- to 10-cm sediment horizons that were directly above layers of gas hydrate. In the Black Sea microbial mats ANME-1 accounted for about 50% of all cells. ANME-2/DSS aggregates occurred in microenvironments within the mat but accounted for only 1% of the total cells. FISH probes for the ANME-2a and ANME-2c subclusters were designed based on a comparative 16S rRNA analysis. In Hydrate Ridge sediments ANME-2a/DSS and ANME-2c/DSS aggregates differed significantly in morphology and abundance. The relative abundance values for these subgroups were remarkably different at Beggiatoa sites (80% ANME-2a, 20% ANME-2c) and Calyptogena sites (20% ANME-2a, 80% ANME-2c), indicating that there was preferential selection of the groups in the two habitats.

Tab. 1 Mole % CO2 and 13C values of CO2 at baseline and at the final monitoring event for the eight observation wells

During CO2 storage operations in mature oilfields or saline aquifers it is desirable to trace the movement of injected CO2 for verification and safety purposes. We demonstrate the successful use of carbon isotope abundance ratios for tracing the movement of CO2 injected at the Cardium CO2 Storage Monitoring project in Alberta between 2005 and 2007. Injected CO2 had a d13C value of -4.6±1.1 per mil that was more than 10 per mil higher than the carbon isotope ratios of casing gas CO2 prior to CO2 injection with average d13C values ranging from -15.9 to -23.5 per mil. After commencement of CO2 injection, d13C values of casing gas CO2 increased in all observation wells towards those of the injected CO2 consistent with a two-source end-member mixing model. At four wells located in a NE-SW trend with respect to the injection wells, breakthrough of injected CO2 was registered chemically (>50 mol % CO2) and isotopically 1-6 months after commencement of CO2 injection resulting in cumulative CO2 fluxes exceeding 100000 m**3 during the observation period. At four other wells, casing gas CO2 contents remained below 5 mol % resulting in low cumulative CO2 fluxes (<2000 m**3) throughout the entire observation period, but carbon isotope ratios indicated contributions between <30 and 80% of injected CO2. Therefore, we conclude that monitoring the movement of CO2 in the injection reservoir with geochemical and isotopic techniques is an effective approach to determine plume expansion and to identify potential preferential flow paths provided that the isotopic composition of injected CO2 is constant and distinct from that of baseline CO2.

Physical properties of 39 sediment cores from the Bengal Fan and northeast Indian Ocean

We obtained sediment physical properties and geochemical data from 47 piston and gravity cores located in the Bay of Bengal, to study the complex history of the Late Pleistocene run-off from the Ganges and Brahmaputra rivers and its imprint on the Bengal Fan. Grain-size parameters were predicted from core logs of density and velocity to infer sediment transport energy and to distinguish different environments along the 3000-km-long transport path from the delta platform to the lower fan. On the shelf, 27 cores indicate rapidly prograding delta foresets today that contain primarily mud, whereas outer shelf sediment has 25% higher silt contents, indicative of stronger and more stable transport regime, which prevent deposition and expose a Late Pleistocene relic surface. Deposition is currently directed towards the shelf canyon 'Swatch of No Ground', where turbidites are released to the only channel-levee system that is active on the fan during the Holocene. Active growth of the channel-levee system occurred throughout sea-level rise and highstand with a distinct growth phase at the end of the Younger Dryas. Coarse-grained material bypasses the upper fan and upper parts of the middle fan, where particle flow is enhanced as a result of flow-restriction in well-defined channels. Sandier material is deposited mainly as sheet-flow deposits on turbidite-dominated plains at the lower fan. The currently most active part of the fan with 10-40 cm thick turbidites is documented for the central channel including inner levees (e.g., site 40). Site 47 from the lower fan far to the east of the active channel-levee system indicates the end of turbidite sedimentation at 300 ka for that location. That time corresponds to the sea-level lowering during late isotopic stage 9 when sediment supply to the fan increased and led to channel avulsion farther upstream, probably indicating a close relation of climate variability and fan activity. Pelagic deep-sea sites 22 and 28 contain a 630-kyear record of climate response to orbital forcing with dominant 21- and 41-kyear cycles for carbonate and magnetic susceptibility, respectively, pointing to teleconnections of low-latitude monsoonal forcing on the precession band to high-latitude obliquity forcing. Upper slope sites 115, 124, and 126 contain a record of the response to high-frequency climate change in the Dansgaard-Oeschger bands during the last glacial cycle with shared frequencies between 0.75 and 2.5 kyear. Correlation of highs in Bengal Fan physical properties to lows in the d18O record of the GISP2 ice-core suggests that times of greater sediment transport energy in the Bay of Bengal are associated with cooler air temperatures over Greenland. Teleconnections were probably established through moisture and other greenhouse-gas forcing that could have been initiated by instabilities in the methane hydrate reservoir in the oceans.

Planktonic foraminiferal oxygen isotopes of ODP Site 143-865 samples

Cool tropical sea surface temperatures (SSTs) are reported for warm Paleogene greenhouse climates based on the d18O of planktonic foraminiferal tests. These results are difficult to reconcile with models of greenhouse gas-forced climate. It has been suggested that this "cool tropics paradox" arises from postdepositional alteration of foraminiferal calcite, yielding erroneously high d18O values. Recrystallization of foraminiferal tests is cryptic and difficult to quantify, and the compilation of robust d18O records from moderately altered material remains challenging. Scanning electron microscopy of planktonic foraminiferal chamber-wall cross sections reveals that the basal area of muricae, pustular outgrowths on the chamber walls of species belonging to the genus Morozovella, contain no mural pores and may be less susceptible to postdepositional alteration. We analyzed the d18O in muricae bases of morozovellids from the central Pacific (Ocean Drilling Program Site 865) by ion microprobe using 10 ?m pits with an analytical reproducibility of ±0.34 per mil (2 standard deviations). In situ measurements of d18O in these domains yield consistently lower values than those published for conventional multispecimen analyses. Assuming that the original d18O is largely preserved in the basal areas of muricae, this new d18O record indicates Early Paleogene (~49-56 Ma) tropical SSTs in the central Pacific were 4°-8°C higher than inferred from the previously published d18O record and that SSTs reached at least ~33°C during the Paleocene-Eocene thermal maximum. This study demonstrates the utility of ion microprobe analysis for generating more reliable paleoclimate records from moderately altered foraminiferal tests preserved in deep-sea sediments.

Planktonic foraminifera, nannofossils and stable isotope record across tghe Paleocene-Eocene thermal maximum of ODP Hole 113-690B

The carbonate saturation profile of the oceans shoaled markedly during a transient global warming event known as the Paleocene-Eocene thermal maximum (PETM) (circa 55 Ma). The rapid release of large quantities of carbon into the ocean-atmosphere system is believed to have triggered this intense episode of dissolution along with a negative carbon isotope excursion (CIE). The brevity (120-220 kyr) of the PETM reflects the rapid enhancement of negative feedback mechanisms within Earth's exogenic carbon cycle that served the dual function of buffering ocean pH and reducing atmospheric greenhouse gas levels. Detailed study of the PETM stratigraphy from Ocean Drilling Program Site 690 (Weddell Sea) reveals that the CIE recovery period, which postdates the CIE onset by ~80 kyr, is represented by an expanded (~2.5 m thick) interval containing a unique planktic foraminiferal assemblage strongly diluted by coccolithophore carbonate. Collectively, the micropaleontological and sedimentological changes preserved within the CIE recovery interval reflect a transient state when ocean-atmosphere chemistry fostered prolific coccolithophore blooms that suppressed the local lysocline to relatively deeper depths. A prominent peak in the abundance of the clay mineral kaolinite is associated with the CIE recovery interval, indicating that continental weathering/runoff intensified at this time as well (Robert and Kennett, 1994). Such parallel stratigraphic changes are generally consonant with the hypothesis that enhanced continental weathering/runoff and carbonate precipitation helped sequester carbon during the PETM recovery period (e.g., Dickens et al., 1997, doi:10.1130/0091-7613(1997)025<0259:ABOGIT>2.3.CO;2 ; Zachos et al., 2005, doi:10.1126/science.1109004).

Sea surface temperatures during Marine Isotope Stage (MIS) 11

The Marine Isotope Stage (MIS) 11 (424-374 ka) was characterized by a protracted deglaciation and an unusually long climatic optimum. It remains unclear to what degree the climate development during this interglacial reflects the unusually weak orbital forcing or greenhouse gas trends. Previously, arguments about the duration and timing of the MIS11 climatic optimum and about the pace of the deglacial warming were based on a small number of key records, which appear to show regional differences. In order to obtain a global signal of climate evolution during MIS11, we compiled a database of 78 sea surface temperature (SST) records from 57 sites spanning MIS11, aligned these individually on the basis of benthic (N = 28) or planktonic (N = 31) stable oxygen isotope curves to a common time frame and subjected 48 of them to an empirical orthogonal function (EOF) analysis. The analysis revealed a high commonality among all records, with the principal SST trend explaining almost 49% of the variability. This trend indicates that on the global scale, the surface ocean underwent rapid deglacial warming during Termination V, in pace with carbon dioxide rise, followed by a broad SST optimum centered at ~410 kyr. The second EOF, which explained ~18% of the variability, revealed the existence of a different SST trend, characterized by a delayed onset of the temperature optimum during MIS11 at ~398 kyr, followed by a prolonged warm period lasting beyond 380 kyr. This trend is most consistently manifested in the mid-latitude North Atlantic and Mediterranean Sea and is here attributed to the strength of the Atlantic meridional overturning circulation. A sensitivity analysis indicates that these results are robust to record selection and to age-model uncertainties of up to 3-6 kyr, but more sensitive to SST seasonal attribution and SST uncertainties >1 °C. In order to validate the CCSM3 (Community Climate System Model, version 3) predictive potential, the annual and seasonal SST anomalies recorded in a total of 74 proxy records were compared with runs for three time slices representing orbital configuration extremes during the peak interglacial of MIS11. The modeled SST anomalies are characterized by a significantly lower variance compared to the reconstructions. Nevertheless, significant correlations between proxy and model data are found in comparisons on the seasonal basis, indicating that the model captures part of the long-term variability induced by astronomical forcing, which appears to have left a detectable signature in SST trends.

Simulated change in atmospheric CO2 and D14C around 14.6 kyr BP, at the onset of the northern hemispheric warming into the Bølling/Allerød

One of the most abrupt and yet unexplained past rises in atmospheric CO2 (10 p.p.m.v. in two centuries) occurred in quasi-synchrony with abrupt northern hemispheric warming into the Bølling/Allerød, 14,600 years ago. Here we use a U/Th-dated record of atmospheric D14C from Tahiti corals to provide an independent and precise age control for this CO2 rise. We also use model simulations to show that the release of old (nearly 14C-free) carbon can explain these changes in CO2 and D14C. The D14C record provides an independent constraint on the amount of carbon released (125 Pg C). We suggest, in line with observations of atmospheric CH4 and terrigenous biomarkers, that thawing permafrost in high northern latitudes could have been the source of carbon, possibly with contribution from flooding of the Siberian continental shelf during meltwater pulse 1A. Our findings highlight the potential of the permafrost carbon reservoir to modulate abrupt climate changes via greenhouse-gas feedbacks.

Stable carbon isotopes, Methane Index, and distributions of different GDGTs

Gas hydrates represent one of the largest pools of readily exchangeable carbon on Earth's surface. Releases of the greenhouse gas methane from hydrates are proposed to be responsible for climate change at numerous events in geological history. Many of these inferred events, however, were based on carbonate carbon isotopes which are susceptible to diagenetic alterations. Here we propose a molecular fossil proxy, i.e., the "Methane Index (MI)", to detect and document the destabilization and dissociation of marine gas hydrates. MI consists of the relative distribution of glycerol dibiphytanyl glycerol tetraethers (GDGTs), the core membrane lipids of archaea. The rational behind MI is that in hydrate-impacted environments, the pool of archaeal tetraether lipids is dominated by GDGT-1, -2 and -3 due to the large contribution of signals from the methanotrophic archaeal community. Our study in the Gulf of Mexico cold-seep sediments demonstrates a correlation between MI and the compound-specific carbon isotope of GDGTs, which is strong evidence supporting the MI-methane consumption relationship. Preliminary applications of MI in a number of hydrate-impacted and/or methane-rich environments show diagnostic MI values, corroborating the idea that MI may serve as a robust indicator for hydrate dissociation that is useful for studies of global carbon cycling and paleoclimate change.

Mg/Ca and stable isotope record of benthic foraminifera of the Paleocene-Eocene Thermal Maximum

A rapid increase in greenhouse gas levels is thought to have fueled global warming at the Paleocene-Eocene Thermal Maximum (PETM). Foraminiferal magnesium/calcium ratios indicate that bottom waters warmed by 4° to 5°C, similar to tropical and subtropical surface ocean waters, implying no amplification of warming in high-latitude regions of deep-water formation under ice-free conditions. Intermediate waters warmed before the carbon isotope excursion, in association with downwelling in the North Pacific and reduced Southern Ocean convection, supporting changing circulation as the trigger for methane hydrate release. A switch to deep convection in the North Pacific at the PETM onset could have amplified and sustained warming.

Stable isotope record of foraminifera from sediment core EW9209-1JPC

Stable isotopic measurements of G. sacculifer and C. wuellerstorfi in a core from the western equatorial Atlantic imply that there are parallel, suborbital oscillations in surface water hydrography and deep water circulation occurring during oxygen isotope stages 2 and 3. Low values of G. sacculifer delta18O accompany high values of C. wuellerstorfi delta13C, linking warmer sea surface temperatures (SSTs) in the tropics with increased production of lower North Atlantic Deep Water (NADW). The amplitude of the delta18O oscillations is 0.6 per mil (or 2°-3°C), which is superimposed on a glacial/interglacial amplitude of about 2.1per mil. Using the G. sacculifer delta18O data, we calculate that surface waters were colder during stage 2 than calculated by CLIMAP [1976, 1981]. The longer-period (>2 kyr) oscillations in air temperature recorded in the Greenland and Antarctic ice cores appear to correlate with oscillations in sea surface temperature in the equatorial Atlantic. The magnitude of these oscillations in tropical SST is too large to have resulted from changes in meridional heat transport caused by the global conveyor alone. The apparent synchroneity of equatorial SST and polar air temperature changes, as well as the amplitude of the SST changes at the equator, are consistent with the climate effects expected from changes in the atmosphere's greenhouse gas content (H2Ovapor, CO2, and CH4).