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.

Stable isotope composition of early Danian sediments from the Atlantic Ocean

A long-standing question in Paleogene climate concerns the frequency and mechanism of transient greenhouse gas-driven climate shifts (hyperthermals). The discovery of the greenhouse gas-driven Paleocene-Eocene Thermal Maximum (PETM; ~55 Ma) has spawned a search for analogous events in other parts of the Paleogene record. On the basis of high-resolution bulk sediment and foraminiferal stable isotope analyses performed on three lower Danian sections of the Atlantic Ocean, we report the discovery of a possible greenhouse gas-driven climatic event in the earliest Paleogene. This event - that we term the Dan-C2 event - is characterized by a conspicuous double negative excursion in delta13C and delta18O, associated with a double spike in increased clay content and decreased carbonate content. This suggests a double period of transient greenhouse gas-driven warming and dissolution of carbonates on the seafloor analogous to the PETMin the early Paleocene at ~65.2 Ma. However, the shape of the two negative carbon isotope excursions that make up the Dan-C2 event is different from the PETM carbon isotope profile. In the Dan-C2 event, these excursions are fairly symmetrical and each persisted for about ~40 ky and are separated by a short plateau that brings the combined duration to ~100 ky, suggesting a possible orbital control on the event. Because of the absence of a long recovery phase, we interpret the Dan-C2 event to have been associated with a redistribution of carbon that was already in the biosphere. The Dan-C2 event and other early Paleogene hyperthermals such as the short-lived early Eocene ELMO eventmay reflect amplification of a regular cycle in the size and productivity of the marine biosphere and the balance between burial of organic and carbonate carbon.

Sea surface temperture reconstruction for ODP Site 189-1172

Relative to the present day, meridional temperature gradients in the Early Eocene age (~56-53 Myr ago) were unusually low, with slightly warmer equatorial regions (Pearson et al., 2007, doi:10.1130/G23175A.1 ) but with much warmer subtropical Arctic (Sluijs et al., 2008, doi:10.1029/2007PA001495) and mid-latitude (Sluijs et al., 2007, doi:10.1038/nature06400) climates. By the end of the Eocene epoch (~34 Myr ago), the first major Antarctic ice sheets had appeared (Zachos et al., 1992, doi:10.1130/0091-7613(1992)020<0569:EOISEO>2.3.CO;2; Barker et al., 2007, doi:10.1016/j.dsr2.2007.07.027), suggesting that major cooling had taken place. Yet the global transition into this icehouse climate remains poorly constrained, as only a few temperature records are available portraying the Cenozoic climatic evolution of the high southern latitudes. Here we present a uniquely continuous and chronostratigraphically well-calibrated TEX86 record of sea surface temperature (SST) from an ocean sediment core in the East Tasman Plateau (palaeolatitude ~65° S). We show that southwest Pacific SSTs rose above present-day tropical values (to ~34° C) during the Early Eocene age (~53 Myr ago) and had gradually decreased to about 21° C by the early Late Eocene age (~36 Myr ago). Our results imply that there was almost no latitudinal SST gradient between subequatorial and subpolar regions during the Early Eocene age (55-50 Myr ago). Thereafter, the latitudinal gradient markedly increased. In theory, if Eocene cooling was largely driven by a decrease in atmospheric greenhouse gas concentration Zachos et al. (2008, doi:10.1038/nature06588), additional processes are required to explain the relative stability of tropical SSTs given that there was more significant cooling at higher latitudes.

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.

Hydrate abundances at the high-flux Batumi seep area in the southeastern Black Sea

Submarine gas hydrates are a major global reservoir of the potent greenhouse gas methane. Since current assessments of worldwide hydrate-bound carbon vary by one order of magnitude, new technical efforts are required for improved and accurate hydrate quantifications. Here we present hydrate abundances determined for surface sediments at the high-flux Batumi seep area in the southeastern Black Sea at 840 m water depth using state-of-the art autoclave technology. Pressure sediment cores of up to 2.65 m in length were recovered with an autoclave piston corer backed by conventional gravity cores. Quantitative core degassing yielded volumetric gas/bulk sediment ratios of up to 20.3 proving hydrate presence. The cores represented late glacial to Holocene hemipelagic sediments with the shallowest hydrates found at 90 cmbsf. Calculated methane concentrations in the different cores surpassed methane equilibrium concentrations in the two lowermost lithological Black Sea units sampled. The results indicated hydrate fractions of 5.2% of pore volume in the sapropelic Unit 2 and mean values of 21% pore volume in the lacustrine Unit 3. We calculate that the studied area of ~ 0.5 km**2 currently contains about 11.3 kt of methane bound in shallow hydrates. Episodic detachment and rafting of such hydrates is suggested by a rugged seafloor topography along with variable thicknesses in lithologies. We propose that sealing by hydrate precipitation in coarse-grained deposits and gas accumulation beneath induces detachment of hydrate/sediment chunks. Floating hydrates will rapidly transport methane into shallower waters and potentially to the sea-atmosphere boundary. In contrast, persistent in situ dissociation of shallow hydrates appears unlikely in the near future as deep water warming by about 1.6 °C and/or decrease in hydrostatic pressure corresponding to a sea level drop of about 130 m would be required. Because hydrate detachment should be primarily controlled by internal factors in this area and in similar hydrated settings, it serves as source of methane in shallow waters and the atmosphere which is mainly decoupled from external forcing.

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.

Table 6 Generalised properties of Utsira caprocks, based on analysis of cuttings (page 79)

Executive Summary: Carbon dioxide capture and storage (CCS) is one option for mitigating atmospheric emissions of carbon dioxide and thereby contributes in actions for stabilization of atmospheric greenhouse gas concentrations. The Bellona Foundation is striving to achieve wide implementation of carbon dioxide (CO2) capture and storage both in Norway and internationally. Bellona considers CCS as the only viable large scale option to close the gap between energy production and demand in an environmentally sound way, thereby ensuring that climate changes and acidification of the oceans due to increased CO2 concentrations in the atmosphere will be stabilised. ff

Database of Ice-Rich Yedoma Permafrost (IRYP)

Vast portions of Arctic and sub-Arctic Siberia, Alaska and the Yukon Territory are covered by ice-rich silty to sandy deposits that are containing large ice wedges, resulting from syngenetic sedimentation and freezing. Accompanied by wedge-ice growth in polygonal landscapes, the sedimentation process was driven by cold continental climatic and environmental conditions in unglaciated regions during the late Pleistocene, inducing the accumulation of the unique Yedoma deposits up to >50 meters thick. Because of fast incorporation of organic material into syngenetic permafrost during its formation, Yedoma deposits include well-preserved organic matter. Ice-rich deposits like Yedoma are especially prone to degradation triggered by climate changes or human activity. When Yedoma deposits degrade, large amounts of sequestered organic carbon as well as other nutrients are released and become part of active biogeochemical cycling. This could be of global significance for future climate warming as increased permafrost thaw is likely to lead to a positive feedback through enhanced greenhouse gas fluxes. Therefore, a detailed assessment of the current Yedoma deposit coverage and its volume is of importance to estimate its potential response to future climate changes. We synthesized the map of the coverage and thickness estimation, which will provide critical data needed for further research. In particular, this preliminary Yedoma map is a great step forward to understand the spatial heterogeneity of Yedoma deposits and its regional coverage. There will be further applications in the context of reconstructing paleo-environmental dynamics and past ecosystems like the mammoth-steppe-tundra, or ground ice distribution including future thermokarst vulnerability. Moreover, the map will be a crucial improvement of the data basis needed to refine the present-day Yedoma permafrost organic carbon inventory, which is assumed to be between 83±12 (Strauss et al., 2013, doi:10.1002/2013GL058088) and 129±30 (Walter Anthony et al., 2014, doi:10.1038/nature13560) gigatonnes (Gt) of organic carbon in perennially-frozen archives. Hence, here we synthesize data on the circum-Arctic and sub-Arctic distribution and thickness of Yedoma for compiling a preliminary circum-polar Yedoma map. For compiling this map, we used (1) maps of the previous Yedoma coverage estimates, (2) included the digitized areas from Grosse et al. (2013) as well as extracted areas of potential Yedoma distribution from additional surface geological and Quaternary geological maps (1.: 1:500,000: Q-51-V,G; P-51-A,B; P-52-A,B; Q-52-V,G; P-52-V,G; Q-51-A,B; R-51-V,G; R-52-V,G; R-52-A,B; 2.: 1:1,000,000: P-50-51; P-52-53; P-58-59; Q-42-43; Q-44-45; Q-50-51; Q-52-53; Q-54-55; Q-56-57; Q-58-59; Q-60-1; R-(40)-42; R-43-(45); R-(45)-47; R-48-(50); R-51; R-53-(55); R-(55)-57; R-58-(60); S-44-46; S-47-49; S-50-52; S-53-55; 3.: 1:2,500,000: Quaternary map of the territory of Russian Federation, 4.: Alaska Permafrost Map). The digitalization was done using GIS techniques (ArcGIS) and vectorization of raster Images (Adobe Photoshop and Illustrator). Data on Yedoma thickness are obtained from boreholes and exposures reported in the scientific literature. The map and database are still preliminary and will have to undergo a technical and scientific vetting and review process. In their current form, we included a range of attributes for Yedoma area polygons based on lithological and stratigraphical information from the original source maps as well as a confidence level for our classification of an area as Yedoma (3 stages: confirmed, likely, or uncertain). In its current version, our database includes more than 365 boreholes and exposures and more than 2000 digitized Yedoma areas. We expect that the database will continue to grow. In this preliminary stage, we estimate the Northern Hemisphere Yedoma deposit area to cover approximately 625,000 km². We estimate that 53% of the total Yedoma area today is located in the tundra zone, 47% in the taiga zone. Separated from west to east, 29% of the Yedoma area is found in North America and 71 % in North Asia. The latter include 9% in West Siberia, 11% in Central Siberia, 44% in East Siberia and 7% in Far East Russia. Adding the recent maximum Yedoma region (including all Yedoma uplands, thermokarst lakes and basins, and river valleys) of 1.4 million km² (Strauss et al., 2013, doi:10.1002/2013GL058088) and postulating that Yedoma occupied up to 80% of the adjacent formerly exposed and now flooded Beringia shelves (1.9 million km², down to 125 m below modern sea level, between 105°E - 128°W and >68°N), we assume that the Last Glacial Maximum Yedoma region likely covered more than 3 million km² of Beringia. Acknowledgements: This project is part of the Action Group "The Yedoma Region: A Synthesis of Circum-Arctic Distribution and Thickness" (funded by the International Permafrost Association (IPA) to J. Strauss) and is embedded into the Permafrost Carbon Network (working group Yedoma Carbon Stocks). We acknowledge the support by the European Research Council (Starting Grant #338335), the German Federal Ministry of Education and Research (Grant 01DM12011 and "CarboPerm" (03G0836A)), the Initiative and Networking Fund of the Helmholtz Association (#ERC-0013) and the German Federal Environment Agency (UBA, project UFOPLAN FKZ 3712 41 106).

Stable carbon isotope ratios of alkane of IODP Hole 302-M0004A

The Palaeocene/Eocene thermal maximum represents a period of rapid, extreme global warming approx ~55 million years ago, superimposed on an already warm world (Zachos et al., 2003, doi:10.1126/science.1090110; Bowen et al., 2004, doi:10.1038/nature03115; Thomas et al., 2002, doi:10.1130/0091-7613(2002)030<1067:WTFFTF>2.0.CO;2). This warming is associated with a severe shoaling of the ocean calcite compensation depth **4 and a >2.5 per mil negative carbon isotope excursion in marine and soil carbonates (Zachos et al., 2003, doi:10.1126/science.1090110; Bowen et al., 2004, doi:10.1038/nature03115; Thomas et al., 2002, doi:10.1130/0091-7613(2002)030<1067:WTFFTF>2.0.CO;2; Zachos et al., doi:10.1126/science.1109004). Together these observations indicate a massive release of 13C-depleted carbon (Zachos et al., doi:10.1126/science.1109004) and greenhouse-gas-induced warming. Recently, sediments were recovered from the central Arctic Ocean (Backman et al., 2006, doi:10.2204/iodp.proc.302.2006), providing the first opportunity to evaluate the environmental response at the North Pole at this time. Here we present stable hydrogen and carbon isotope measurements of terrestrial-plant- and aquatic-derived n-alkanes that record changes in hydrology, including surface water salinity and precipitation, and the global carbon cycle. Hydrogen isotope records are interpreted as documenting decreased rainout during moisture transport from lower latitudes and increased moisture delivery to the Arctic at the onset of the Palaeocene/Eocene thermal maximum, consistent with predictions of poleward storm track migrations during global warming (Backman et al., 2006, doi:10.2204/iodp.proc.302.2006). The terrestrial-plant carbon isotope excursion (about ~4.5 to ~6 per mil) is substantially larger than those of marine carbonates. Previously, this offset was explained by the physiological response of plants to increases in surface humidity (Bowen et al., 2004, doi:10.1038/nature03115). But this mechanism is not an effective explanation in this wet Arctic setting, leading us to hypothesize that the true magnitude of the excursion - and associated carbon input - was greater than originally surmised. Greater carbon release and strong hydrological cycle feedbacks may help explain the maintenance of this unprecedented warmth.of this unprecedented warmth.

(Table 1) Gravity and pressure within the Sleipner area relative to ones at the reference station in 2002 and 2005

At Sleipner, CO2 is being separated from natural gas and injected into an underground saline aquifer for environmental purposes. Uncertainty in the aquifer temperature leads to uncertainty in the in situ density of CO2. In this study, gravity measurements were made over the injection site in 2002 and 2005 on top of 30 concrete benchmarks on the seafloor in order to constrain the in situ CO2 density. The gravity measurements have a repeatability of 4.3 µGal for 2003 and 3.5 µGal for 2005. The resulting time-lapse uncertainty is 5.3 µGal. Unexpected benchmark motions due to local sediment scouring contribute to the uncertainty. Forward gravity models are calculated based on both 3D seismic data and reservoir simulation models. The time-lapse gravity observations best fit a high temperature forward model based on the time-lapse 3D seismics, suggesting that the average in situ CO2 density is about to 530kg/m**3. Uncertainty in determining the average density is estimated to be ±65 kg/m**3 (95% confidence), however, this does not include uncertainties in the modeling. Additional seismic surveys and future gravity measurements will put better constraints on the CO2 density and continue to map out the CO2 flow.

(Table S1) Stable carbon and oxygen isotopes and Mg/Ca ratios of planktonic foraminifera from ODP Hole 198-1209B

The Paleocene-Eocene Thermal Maximum (PETM) has been attributed to a rapid rise in greenhouse gas levels. If so, warming should have occurred at all latitudes, although amplified toward the poles. Existing records reveal an increase in high-latitude sea surface temperatures (SSTs) (8° to 10°C) and in bottom water temperatures (4° to 5°C). To date, however, the character of the tropical SST response during this event remains unconstrained. Here we address this deficiency by using paired oxygen isotope and minor element (magnesium/calcium) ratios of planktonic foraminifera from a tropical Pacific core to estimate changes in SST. Using mixed-layer foraminifera, we found that the combined proxies imply a 4° to 5°C rise in Pacific SST during the PETM. These results would necessitate a rise in atmospheric pCO2 to levels three to four times as high as those estimated for the late Paleocene.

Temperature, salinity, methane concentration, oxidation rates and methanotrophic cell number at Svalbard seeps in summer 2011 and 2012 (cruises POS419 and MSM21/4)

Large amounts of the greenhouse gas methane are released from the seabed to the water column where it may be consumed by aerobic methanotrophic bacteria. This microbial filter is consequently the last marine sink for methane before its liberation to the atmosphere. The size and activity of methanotrophic communities, which determine the capacity of the water column methane filter, are thought to be mainly controlled by nutrient and redox dynamics, but little is known about the effects of ocean currents. Here, we report measurements of methanotrophic activity and biomass (CARD-FISH) at methane seeps west of Svalbard, and related them to physical water mass properties (CTD) and modelled current dynamics. We show that cold bottom water containing a large number of aerobic methanotrophs was rapidly displaced by warmer water with a considerably smaller methanotrophic community. This water mass exchange, caused by short-term variations of the West Spitsbergen Current, constitutes a rapid oceanographic switch severely reducing methanotrophic activity in the water column. Strong and fluctuating currents are widespread oceanographic features common at many methane seep systems and are thus likely to globally affect methane oxidation in the ocean water column.

Half-degree gridded nitrogen and phosphorus fertilizer use for global agriculture production during 1900-2013

In addition to enhance agricultural productivity, synthetic nitrogen (N) and phosphorous (P) fertilizer application in croplands dramatically altered global nutrient budget, water quality, greenhouse gas balance, and their feedbacks to the climate system. However, due to the lack of geospatial fertilizer input data, current Earth system/land surface modeling studies have to ignore or use over-simplified data (e.g., static, spatially uniform fertilizer use) to characterize agricultural N and P input over decadal or century-long period. We therefore develop a global time-series gridded data of annual synthetic N and P fertilizer use rate in croplands, matched with HYDE 3,2 historical land use maps, at a resolution of 0.5º latitude by longitude during 1900-2013. Our data indicate N and P fertilizer use rates increased by approximately 8 times and 3 times, respectively, since the year 1961, when IFA (International Fertilizer Industry Association) and FAO (Food and Agricultural Organization) survey of country-level fertilizer input were available. Considering cropland expansion, increase of total fertilizer consumption amount is even larger. Hotspots of agricultural N fertilizer use shifted from the U.S. and Western Europe in the 1960s to East Asia in the early 21st century. P fertilizer input show the similar pattern with additional hotspot in Brazil. We find a global increase of fertilizer N/P ratio by 0.8 g N/g P per decade (p< 0.05) during 1961-2013, which may have important global implication of human impacts on agroecosystem functions in the long run. Our data can serve as one of critical input drivers for regional and global assessment on agricultural productivity, crop yield, agriculture-derived greenhouse gas balance, global nutrient budget, land-to-aquatic nutrient loss, and ecosystem feedback to the climate system.

Isotopic data of atmospheric methane from EDML and Vostok ice cores over a full glacial cycle

The response of natural CH4 sources to climate changes will be an important factor to consider as concentrations of this potent greenhouse gas continue to increase. Polar ice cores provide the means to assess this sensitivity in the past and have shown a close connection between CH4 levels and northern hemisphere temperature variability over the last glacial cycle. However, the contribution of the various CH4 sources and sinks to these changes is still a matter of debate. Contemporaneous stable CH4 isotope records in ice cores provide additional boundary conditions for assessing changes in the CH4 sources and sinks. Here we present new ice core CH4 isotope data covering the last 160,000 years, showing a clear decoupling between CH4 loading and carbon isotopic variations over most of the record. We suggest that d13CH4 variations were not dominated by a change in the source mix but rather by climate- and CO2-related ecosystem control on the isotopic composition of the methane precursor material, especially in seasonally inundated wetlands in the tropics. In contrast, relatively stable d13CH4 intervals occurred during large CH4 loading changes concurrently with past climate changes implying that most CH4 sources (most notably tropical wetlands) responded simultaneously.

(Table S1) Radiocarbon dates of organic matter from Lena Delta ice wedges

Relative to the past 2,000 years, the Arctic region has warmed significantly over the past few decades. However, the evolution of Arctic temperatures during the rest of the Holocene is less clear. Proxy reconstructions, suggest a long-term cooling trend throughout the mid- to late Holocene, whereas climate model simulations show only minor changes or even warming. Here we present a record of the oxygen isotope composition of permafrost ice wedges from the Lena River Delta in the Siberian Arctic. The isotope values, which reflect winter season temperatures, became progressively more enriched over the past 7,000 years, reaching unprecedented levels in the past five decades. This warming trend during the mid- to late Holocene is in opposition to the cooling seen in other proxy records. However, most of these existing proxy records are biased towards summer temperatures. We argue that the opposing trends are related to the seasonally different orbital forcing over this interval. Furthermore, our reconstructed trend as well as the recent maximum are consistent with the greenhouse gas forcing and climate model simulations, thus reconciling differing estimates of Arctic and northern high-latitude temperature evolution during the Holocene.