(Supplemental Table 2) Plant characteristics, GHG fluxes, and soil chemical and microbial properties in experimental treatments in Alexandra Fiord

We evaluated above- and belowground ecosystem changes in a 16 year, combined fertilization and warming experiment in a High Arctic tundra deciduous shrub heath (Alexandra Fiord, Ellesmere Island, NU, Canada). Soil emissions of the three key greenhouse gases (GHGs) (carbon dioxide, methane, and nitrous oxide) were measured in mid-July 2009 using soil respiration chambers attached to a FTIR system. Soil chemical and biochemical properties including Q10 values for CO2, CH4, and N2O, Bacteria and Archaea assemblage composition, and the diversity and prevalence of key nitrogen cycling genes including bacterial amoA, crenarchaeal amoA, and nosZ were measured. Warming and fertilization caused strong increases in plant community cover and height but had limited effects on GHG fluxes and no substantial effect on soil chemistry or biochemistry. Similarly, there was a surprising lack of directional shifts in the soil microbial community as a whole or any change at all in microbial functional groups associated with CH4 consumption or N2O cycling in any treatment. Thus, it appears that while warming and increased nutrient availability have strongly affected the plant community over the last 16 years, the belowground ecosystem has not yet responded. This resistance of the soil ecosystem has resulted in limited changes in GHG fluxes in response to the experimental treatments.

Mesocosm experiment Cape Verde 2012: Data

Two 7-day mesocosm experiments were conducted in October 2012 at the Instituto Nacional de Desenvolvimento das Pescas (INDP), Mindelo, Cape Verde. Surface water was collected at night before the start of the respective experiment with RV Islândia south of São Vicente (16°44.4'N, 25°09.4'W) and transported to shore using four 600L food safe intermediate bulk containers. Sixteen mesocosm bags were distributed in four flow-through water baths and shaded with blue, transparent lids to approximately 20% of surface irradiation. Mesocosm bags were filled from the containers by gravity, using a submerged hose to minimize bubbles. The accurate volume inside the individual bags was calculated after addition of 1.5 mmol silicate and measuring the resulting silicate concentration. The volume ranged from 105.5 to 145 L. The experimental manipulation comprised addition of different amounts of inorganic N and P. In the first experiment, the P supply was changed at constant N supply in thirteen of the sixteen units, while in the second experiment the N supply was changed at constant P supply in twelve of the sixteen units. In addition to this, "cornerpoints" were chosen that were repeated during both experiments. Four cornerpoints should have been repeated, but setting the nutrient levels in one mesocosm was not succesfull and therefore this mesocosm also was set at the center point conditions. Experimental treatments were evenly distributed between the four water baths. Initial sampling of the mesocosms on day 1 of each run was conducted between 9:45 and 11:30. After nutrient manipulation, sampling was conducted on a daily basis between 09:00 and 10:30 for days 2 to 8.

EPOCA Svalbard 2009 benthic experiment: Serripes study, 2009

Ocean acidification influences sediment/water nitrogen fluxes, possibly by impacting on the microbial process of ammonia oxidation. To investigate this further, undisturbed sediment cores collected from Ny Alesund harbour (Svalbard) were incubated with seawater adjusted to CO2 concentrations of 380, 540, 760, 1,120 and 3,000 µatm. DNA and RNA were extracted from the sediment surface after 14 days' exposure and the abundance of bacterial and archaeal ammonia oxidising (amoA) genes and transcripts quantified using quantitative polymerase chain reaction. While there was no change to the abundance of bacterial amoA genes, an increase to 760 µatm pCO2 reduced the abundance of bacterial amoA transcripts by 65 %, and this was accompanied by a shift in the composition of the active community. In contrast, archaeal amoA gene and transcript abundance both doubled at 3,000 µatm, with an increase in species richness also apparent. This suggests that ammonia oxidising bacteria and archaea in marine sediments have different pH optima, and the impact of elevated CO2 on N cycling may be dependent on the relative abundances of these two major microbial groups. Further evidence of a shift in the balance of key N cycling groups was also evident: the abundance of nirS-type denitrifier transcripts decreased alongside bacterial amoA transcripts, indicating that NO3 ? produced by bacterial nitrification fuelled denitrification. An increase in the abundance of Planctomycete-specific 16S rRNA, the vast majority of which grouped with known anammox bacteria, was also apparent at 3,000 µatm pCO2. This could indicate a possible shift from coupled nitrification-denitrification to anammox activity at elevated CO2.

Response of spring diatoms to CO2 availability in the Western North Pacific as determined by next-generation sequencing

Next-generation sequencing (NGS) technologies have enabled us to determine phytoplankton community compositions at high resolution. However, few studies have adopted this approach to assess the responses of natural phytoplankton communities to environmental change. Here, we report the impact of different CO2 levels on spring diatoms in the Oyashio region of the western North Pacific as estimated by NGS of the diatom-specific rbcL gene (DNA), which encodes the large subunit of RubisCO. We also examined the abundance and composition of rbcL transcripts (cDNA) in diatoms to assess their physiological responses to changing CO2 levels. A short-term (3-day) incubation experiment was carried out on-deck using surface Oyashio waters under different pCO2 levels (180, 350, 750, and 1000 µatm) in May 2011. During the incubation, the transcript abundance of the diatom-specific rbcL gene decreased with an increase in seawater pCO2 levels. These results suggest that CO2 fixation capacity of diatoms decreased rapidly under elevated CO2 levels. In the high CO2 treatments (750 and 1000 µatm), diversity of diatom-specific rbcL gene and its transcripts decreased relative to the control treatment (350µatm), as well as contributions of Chaetocerataceae, Thalassiosiraceae, and Fragilariaceae to the total population, but the contributions of Bacillariaceae increased. In the low CO2 treatment, contributions of Bacillariaceae also increased together with other eukaryotes. These suggest that changes in CO2 levels can alter the community composition of spring diatoms in the Oyashio region. Overall, the NGS technology provided us a deeper understanding of the response of diatoms to changes in CO2 levels in terms of their community composition, diversity, and photosynthetic physiology.

Phospholipid fatty acid and quinone composition of sediments from ODP Leg 205 sites

Assessing the habitability of deep-sea sediments undergoing compaction, compression, and subduction at convergent margins adds to our understanding of the limits of the terrestrial biosphere. In this work, we report exploratory biomarker data on sediments obtained at Ocean Drilling Program (ODP) Sites 1253, 1254, and 1255 during drilling at the Costa Rica subduction trench and forearc sedimentary wedge. The samples selected for postcruise biomarker analyses were located within intervals of potentially enhanced fluid flow within the décollement and sedimentary wedge fault zones (Sites 1254 and 1255) and within basal carbonates at the reference site (Site 1253). The passage of fluids that are geochemically distinct from ambient interstitial water provides a disequilibrium setting that may enhance habitability. Biomarker data show low levels of microbial biomass in subseafloor sediments sampled at the Costa Rica convergent margin as deep as ~370 meters below seafloor.

(Table S2) Relative abundance of archaeal 16S rRNA genes in sediment cores

Deep drilling into the marine sea floor has uncovered a vast sedimentary ecosystem of microbial cells (Parkes et al., 1994, doi:10.1038/371410a0; D'Hondt et al., 2004, doi:10.1126/science.1101155). Extrapolation of direct counts of stained microbial cells to the total volume of habitable marine subsurface sediments suggests that between 56 Pg (Parkes et al., 1994, doi:10.1038/371410a0) and 303 Pg (Whitman et al., 1998) of cellular carbon could be stored in this largely unexplored habitat. From recent studies using various culture-independent techniques, no clear picture has yet emerged as to whether Archaea or Bacteria are more abundant in this extensive ecosystem (Schippers et al., doi:10.1038/nature03302; Inagaki et al., doi:10.1073/pnas.0511033103 ; Mauclaire et al., doi:10.1111/j.1472-4677.2004.00035.x; Biddle et al., doi:10.1073/pnas.0600035103). Here we show that in subsurface sediments buried deeper than 1 m in a wide range of oceanographic settings at least 87% of intact polar membrane lipids, biomarkers for the presence of live cells (Biddle et al., doi:10.1073/pnas.0600035103; Sturt et al., 2004, doi:10.1002/rcm.1378), are attributable to archaeal membranes, suggesting that Archaea constitute a major fraction of the biomass. Results obtained from modified quantitative polymerase chain reaction and slot-blot hybridization protocols support the lipid-based evidence and indicate that these techniques have previously underestimated archaeal biomass. The lipid concentrations are proportional to those of total organic carbon. On the basis of this relationship, we derived an independent estimate of amounts of cellular carbon in the global marine subsurface biosphere. Our estimate of 90 Pg of cellular carbon is consistent, within an order of magnitude, with previous estimates, and underscores the importance of marine subsurface habitats for global biomass budgets.

(Table 1) Basin area, average late-summer speedup, number of sinks and runoff of West Greenland Ice Sheet catchments

We use interferometric synthetic aperture radar observations recorded in a land-terminating sector of western Greenland to characterise the ice sheet surface hydrology and to quantify spatial variations in the seasonality of ice sheet flow. Our data reveal a non-uniform pattern of late-summer ice speedup that, in places, extends over 100 km inland. We show that the degree of late-summer speedup is positively correlated with modelled runoff within the 10 glacier catchments of our survey, and that the pattern of late-summer speedup follows that of water routed at the ice sheet surface. In late-summer, ice within the largest catchment flows on average 48% faster than during winter, whereas changes in smaller catchments are less pronounced. Our observations show that the routing of seasonal runoff at the ice sheet surface plays an important role in shaping the magnitude and extent of seasonal ice sheet speedup.

(Appendix) Vertebral number of topsmelt, Atherinops affinis, adults and offspring from different latitudes along the North-American Pacific coast

Organisms that are distributed across spatial climate gradients often exhibit adaptive local variations in morphological and physiological traits, but to what extent such gradients shape evolutionary responses is still unclear. Given the strong natural contrast in latitudinal temperature gradients between the North-American Pacific and Atlantic coast, we asked how increases in vertebral number (VN, known as Jordan's Rule) with latitude would differ between Pacific (Atherinops affinis) and Atlantic Silversides (Menidia menidia), two ecologically equivalent and taxonomically similar fishes with similar latitudinal distributions. VN was determined from radiographs of wild-caught adults (genetic + environmental differences) and its genetic basis confirmed by rearing offspring in common garden experiments. Compared to published data on VN variation in M. menidia (a mean increase of 7.0 vertebrae from 32 to 46°N, VN slope = 0.42/lat), the latitudinal VN increase in Pacific Silversides was approximately half as strong (a mean increase of 3.3 vertebrae from 28 to 43°N, VN slope = 0.23/lat). This mimicked the strong Atlantic (1.11°C/lat) versus weak Pacific latitudinal gradient (0.40°C/lat) in median annual sea surface temperature (SST). Importantly, the relationship of VN to SST was not significantly different between the two species (average slope = -0.39 vertebrae/°C), thus suggesting a common thermal dependency of VN in silverside fishes. Our findings provide novel support for the hypothesis that temperature gradients are the ultimate cause of Jordan's Rule, even though its exact adaptive significance remains speculative. A second investigated trait, the mode of sex determination in Atlantic versus Pacific Silversides, revealed patterns that were inconsistent with our expectation: M. menidia displays temperature-dependent sex determination (TSD) at low latitudes, where growing seasons are long or unconstrained, but also a gradual shift to genetic sex determination (GSD) with increasing latitude due to more and more curtailed growing seasons. Sex ratios in A. affinis, on the other hand, were independent of latitude and rearing temperature (indicating GSD), even though growing seasons are thermally unconstrained across most of the geographical distribution of A. affinis. This suggests that additional factors (e.g., longevity) play an important role in shaping the mode of sex determination in silverside fishes.

Siliceous microfossil assemblage and accumulation rate variations across the Middle Eocene Climatic Optimum at ODP Site 171-1051

The Middle Eocene Climatic Optimum (MECO; ~ 40 million years ago [Ma]) is one of the most prominent transient global warming events in the Paleogene. Although the event is well documented in geochemical and isotopic proxy records at many locations, the marine biotic response to the MECO remains poorly constrained. We present new high-resolution, quantitative records of siliceous microplankton assemblages from the MECO interval of Ocean Drilling Program (ODP) Site 1051 in the subtropical western North Atlantic Ocean, which are interpreted in the context of published foraminiferal and bulk carbonate stable isotope (d18O and d13C) records. High diatom, radiolarian and silicoflagellate accumulation rates between 40.5 and 40.0 Ma are interpreted to reflect an ~ 500 thousand year (kyr) interval of increased nutrient supply and resultant surface-water eutrophication that was associated with elevated sea-surface temperatures during the prolonged onset of the MECO. Relatively low pelagic siliceous phytoplankton sedimentation accompanied the peak MECO warming interval and the termination of the MECO during an ~ 70 kyr interval centered at ~ 40.0 Ma. Following the termination of the MECO, an ~ 200-kyr episode of increased siliceous plankton abundance indicates enhanced nutrient levels between ~ 39.9 and 39.7 Ma. Throughout the Site 1051 record, abundance and accumulation rate fluctuations in neritic diatom taxa are similar to the trends observed in pelagic taxa, implying either similar controls on diatom production in the neritic and pelagic zones of the western North Atlantic or fluctuations in sea level and/or shelf accommodation on the North American continental margin to the west of Site 1051. These results, combined with published records based on multiple proxies, indicate a geographically diverse pattern of surface ocean primary production changes across the MECO. Notably, however, increased biosiliceous accumulation is recorded at both ODP Sites 1051 and 748 (Southern Ocean) in response to MECO warming. This may suggest that increased biosiliceous sediment accumulation, if indeed a widespread phenomenon, resulted from higher continental silicate weathering rates and an increase in silicic acid supply to the oceans over several 100 kyr during the MECO.

Porewater concentrations, pH, single cell number, anaerobic oxidation of methane (AOM) rate, sulfate reduction (SR) rate of samples during METEOR cruise M76/3b at the REGAB cold seep site from the West African margin in 2008

The giant pockmark REGAB (West African margin, 3160 m water depth) is an active methane-emitting cold seep ecosystem, where the energy derived from microbially mediated oxidation of methane supports high biomass and diversity of chemosynthetic communities. Bare sediments interspersed with heterogeneous chemosynthetic assemblages of mytilid mussels, vesicomyid clams and siboglinid tubeworms form a complex seep ecosystem. To better understand if benthic bacterial communities reflect the patchy distribution of chemosynthetic fauna, all major chemosynthetic habitats at REGAB were investigated using an interdisciplinary approach combining porewater geochemistry, in situ quantification of fluxes and consumption of methane, as well bacterial community fingerprinting. This study revealed that sediments populated by different fauna assemblages show distinct biogeochemical activities and are associated with distinct sediment bacterial communities. The methane consumption and methane effluxes ranged over one to two orders of magnitude across habitats, and reached highest values at the mussel habitat, which hosted a different bacterial community compared to the other habitats. Clam assemblages had a profound impact on the sediment geochemistry, but less so on the bacterial community structure. Moreover, all clam assemblages at REGAB were restricted to sediments characterized by complete methane consumption in the seafloor, and intermediate biogeochemical activity. Overall, variations in the sediment geochemistry were reflected in the distribution of both fauna and microbial communities; and were mostly determined by methane flux.