Data S1. Results of the analyses considering developmental time, supplementary figures and tables & Video S1

Our ability to project the impact of global change on marine ecosystem is limited by our poor understanding on how to predict species sensitivity. For example, the impact of ocean acidification is highly species-specific, even in closely related taxa. The aim of this study was to test the hypothesis that the tolerance range of a given species to decreased pH corresponds to their natural range of exposure. Larvae of the green sea urchin Strongylocentrotus droebachiensis were cultured from fertilization to metamorphic competence (29 days) under a wide range of pH (from pHT = 8.0/pCO2 ~ 480 ?atm to pHT = 6.5/pCO2 ~ 20 000 ?atm) covering present (from pHT 8.7 to 7.6), projected near-future variability (from pHT 8.3 to 7.2) and beyond. Decreasing pH impacted all tested parameters (mortality, symmetry, growth, morphometry and respiration). Development of normal, although showing morphological plasticity, swimming larvae was possible as low as pHT >= 7.0. Within that range, decreasing pH increased mortality and asymmetry and decreased body length (BL) growth rate. Larvae raised at lowered pH and with similar BL had shorter arms and a wider body. Relative to a given BL, respiration rates and stomach volume both increased with decreasing pH suggesting changes in energy budget. At the lowest pHs (pHT <= 6.5), all the tested parameters were strongly negatively affected and no larva survived past 13 days post fertilization. In conclusion, sea urchin larvae appeared to be highly plastic when exposed to decreased pH until a physiological tipping point at pHT = 7.0. However, this plasticity was associated with direct (increased mortality) and indirect (decreased growth) consequences for fitness.

(Table S1) Globigerinoides sacculifer Boron isotope ratios of ODP Hole 165-999A

everal hypotheses have been put forward to explain the onset of intensive glaciations on Greenland, Scandinavia, and North America during the Pliocene epoch between 3.6 and 2.7 million years ago (Ma). A decrease in atmospheric CO2 may have played a role during the onset of glaciations, but other tectonic and oceanic events occurring at the same time may have played a part as well. Here we present detailed atmospheric CO2 estimates from boron isotopes in planktic foraminifer shells spanning 4.6-2.0 Ma. Maximal Pliocene atmospheric CO2 estimates gradually declined from values around 410 µatm to early Pleistocene values of 300 ?atm at 2.0 Ma. After the onset of large-scale ice sheets in the Northern Hemisphere, maximal pCO2 estimates were still at 2.5 Ma +90 µatm higher than values characteristic of the early Pleistocene interglacials. By contrast, Pliocene minimal atmospheric CO2 gradually decreased from 310 to 245 µatm at 3.2 Ma, coinciding with the start of transient glaciations on Greenland. Values characteristic of early Pleistocene glacial atmospheric CO2 of 200 ?atm were abruptly reached after 2.7 Ma during the late Pliocene transition. This trend is consistent with the suggestion that ocean stratification and iron fertilization increased after 2.7 Ma in the North Pacific and Southern Ocean and may have led to increased glacial CO2 storage in the oceanic abyss after 2.7 Ma onward.

(Table S1) Trace element and isotopic compositions of anhydrites sampled in the Manus Basin, Papua New Guinea

Microchemical analyses of rare earth element (REE) concentrations and Sr and S isotope ratios of anhydrite are used to identify sub-seafloor processes governing the formation of hydrothermal fluids in the convergent margin Manus Basin, Papua New Guinea. Samples comprise drill-core vein anhydrite and seafloor massive anhydrite from the PACMANUS (Roman Ruins, Snowcap and Fenway) and SuSu Knolls (North Su) active hydrothermal fields. Chondrite-normalized REE patterns in anhydrite show remarkable heterogeneity on the scale of individual grains, different from the near uniform REEN patterns measured in anhydrite from mid-ocean ridge deposits. The REEN patterns in anhydrite are correlated with REE distributions measured in hydrothermal fluids venting at the seafloor at these vent fields and are interpreted to record episodes of hydrothermal fluid formation affected by magmatic volatile degassing. 87Sr/86Sr ratios vary dramatically within individual grains between that of contemporary seawater and that of endmember hydrothermal fluid. Anhydrite was precipitated from a highly variable mixture of the two. The intra-grain heterogeneity implies that anhydrite preserves periods of contrasting hydrothermal versus seawater dominant near-seafloor fluid circulation. Most sulfate d34S values of anhydrite cluster around that of contemporary seawater, consistent with anhydrite precipitating from hydrothermal fluid mixed with locally entrained seawater. Sulfate d34S isotope ratios in some anhydrites are, however, lighter than that of seawater, which are interpreted as recording a source of sulfate derived from magmatic SO2 degassed from underlying felsic magmas in the Manus Basin. The range of elemental and isotopic signatures observed in anhydrite records a range of sub-seafloor processes including high-temperature hydrothermal fluid circulation, varying extents of magmatic volatile degassing, seawater entrainment and fluid mixing. The chemical and isotopic heterogeneity recorded in anhydrite at the inter- and intra-grain scale captures the dynamics of hydrothermal fluid formation and sub-seafloor circulation that is highly variable both spatially and temporally on timescales over which hydrothermal deposits are formed. Microchemical analysis of hydrothermal minerals can provide information about the temporal history of submarine hydrothermal systems that are variable over time and cannot necessarily be inferred only from the study of vent fluids.

(Table S1) Tie-points used in age model construction for ODP Leg 172 sites

ix Ocean Drilling Program (ODP) sites, in the Northwest Atlantic have been used to investigate kinematic and chemical changes in the "Western Boundary Undercurrent" (WBUC) during the development of full glacial conditions across the Marine Isotope Stage 5a/4 boundary (~70,000 years ago). Sortable silt mean grain size(sort s) measurements are employed to examine changes in near bottom flow speeds, together with carbon isotopes measured in benthic foraminifera and % planktic foraminiferal fragmentation as proxies for changes in water-mass chemistry. A depth transect of cores, spanning 1.8-4.6 km depth, allows changes in both the strength and depth of the WBUC to be constrained across millennial scale events. Sort s measurements reveal that the flow speed structure of the WBUC during warm intervals ("interstadials") was comparable to modern (Holocene) conditions. However, significant differences are observed during cold intervals, with higher relative flow speeds inferred for the shallow component of the WBUC (~2 km depth) during all cold "stadial" intervals (including Heinrich Stadial 6), and a substantial weakening of the deep component (~3-4 km) during full glacial conditions. Our results therefore reveal that the onset of full glacial conditions was associated with a regime shift to a shallower mode of circulation (involving Glacial North Atlantic Intermediate Water) that was quantitatively distinct from preceding cold stadial events. Furthermore, our chemical proxy data show that the physical response of the WBUC during the last glacial inception was probably coupled to basin-wide changes in the water-mass composition of the deep Northwest Atlantic.

(Table S1) Geochemistry for sediment core MD02-2515

A high-resolution, 55-kyr long record of chalcophile and redox-sensitive trace element accumulation (Ag, Cd, Re, Mo) from MD02-2515, western Guaymas Basin, is investigated in conjunction with patterns in stratigraphy and productivity. High opal concentrations (~58 wt. %), representing increased diatom production, coincide with laminated sediments, and dilute the concentrations of organic carbon (Corg) and metals. A similarity between opal and normalized Corg, Ag and Cd concentrations suggests delivery to the sediments by diatom export production, while patterns in normalized Re and Mo accumulation suggest a different emplacement mechanism. Although Mo enrichment in organic-rich, laminated sediments typically represents anoxic conditions at other locations, Mo (and Re) in Guaymas Basin is enriched in nonlaminated and bioturbated sediments that are representative of oxygenated conditions. Adsorption onto Fe- and/or Mn-oxyhydroxide surfaces during oxygenation inadequately explains both the Re and Mo enrichments. Thus, recently published mechanisms invoking direct Re and Mo removal from the water column and bioturbation-assisted irrigation of Re into the sediments are used to explain the counterintuitive observations in Guaymas Basin. The MD02-2515 stratigraphic and proxy records are also different from other records in the northeast Pacific in that there is little correspondence with Greenland Dansgaard-Oeschger interstadials. There is some correlation with Heinrich events, suggesting that ventilation of intermediate waters and/or reduced productivity may be important in controlling stratigraphy and trace element accumulation. The results question whether MD02-2515 records can be compared to northeast Pacific open-margin records, especially before 17 kyr BP.

(Table S1) Carbon and oxygen isotope values, and Mn and Sr concentrations of Early Jurassic sediments

For much of the Mesozoic record there has been an inconclusive debate on the possible global significance of isotopic proxies for environmental change and of sequence stratigraphic depositional sequences. We present a carbon and oxygen isotope and elemental record for part of the Early Jurassic based on marine benthic and nektobenthic molluscs and brachiopods from the shallow marine succession of the Cleveland Basin, UK. The invertebrate isotope record is supplemented with carbon isotope data from fossil wood, which samples atmospheric carbon. New data elucidate two major global carbon isotope events, a negative excursion of ~2 per mil at the Sinemurian-Pliensbachian boundary, and a positive excursion of ~2 per mil in the Late Pliensbachian. The Sinemurian-Pliensbachian boundary event is similar to the slightly younger Toarcian Oceanic Anoxic Event and is characterized by deposition of relatively deepwater organic-rich shale. The Late Pliensbachian strata by contrast are characterized by shallow marine deposition. Oxygen isotope data imply cooling locally for both events. However, because deeper water conditions characterize the Sinemurian-Pliensbachian boundary in the Cleveland Basin the temperature drop is likely of local significance; in contrast a cool Late Pliensbachian shallow seafloor agrees with previous inference of partial icehouse conditions. Both the large-scale, long-term and small-scale, short-duration isotopic cycles occurred in concert with relative sea level changes documented previously from sequence stratigraphy. Isotope events and the sea level cycles are concluded to reflect processes of global significance, supporting the idea of an Early Jurassic in which cyclic swings from icehouse to greenhouse and super greenhouse conditions occurred at timescales from 1 to 10 Ma.

(Tables S1,2) Stable isotopes in muscle tissue of polar bears (Ursus maritimus) from various regions of the Arctic

The relative contribution of regional contamination versus dietary differences to geographic variation in polar bear (Ursus maritimus) contaminant levels is unknown. Dietary variation between Alaska, Canada, East Greenland, and Svalbard subpopulations was assessed by muscle nitrogen and carbon stable isotope (d15N, d13C) and adipose fatty acid (FA) signatures relative to their main prey (ringed seals). Western and southern Hudson Bay signatures were characterized by depleted d15N and d13C, lower proportions of C20 and C22 monounsaturated FAs and higher proportions of C18 and longer chain polyunsaturated FAs. East Greenland and Svalbard signatures were reversed relative to Hudson Bay. Alaskan and Canadian Arctic signatures were intermediate. Between-subpopulation dietary differences predominated over interannual, seasonal, sex, or age variation. Among various brominated and chlorinated contaminants, diet signatures significantly explained variation in adipose levels of polybrominated diphenyl ether (PBDE) flame retardants (14-15%) and legacy PCBs (18-21%). However, dietary influence was contaminant class-specific, since only low or nonsignificant proportions of variation in organochlorine pesticide (e.g., chlordane) levels were explained by diet. Hudson Bay diet signatures were associated with lower PCB and PBDE levels, whereas East Greenland and Svalbard signatures were associated with higher levels. Understanding diet/food web factors is important to accurately interpret contaminant trends, particularly in a changing Arctic.

(Table S1) Stable isotope record and growth rates of cultured Mytilus edulis

To study the effects of temperature, salinity, and life processes (growth rates, size, metabolic effects, and physiological/genetic effects) on newly precipitated bivalve carbonate, we quantified shell isotopic chemistry of adult and juvenile animals of the intertidal bivalve Mytilus edulis (Blue mussel) collected alive from western Greenland and the central Gulf of Maine and cultured them under controlled conditions. Data for juvenile and adult M. edulis bivalves cultured in this study, and previously by Wanamaker et al. (2006, doi:10.1029/2005GC001189), yielded statistically identical paleotemperature relationships. On the basis of these experiments we have developed a species-specific paleotemperature equation for the bivalve M. edulis [T °C = 16.28 (±0.10) - 4.57 (±0.15) {d18Oc VPBD - d18Ow VSMOW} + 0.06 (±0.06) {d18Oc VPBD - d18Ow VSMOW}**2; r**2 = 0.99; N = 323; p < 0.0001]. Compared to the Kim and O'Neil (1997) inorganic calcite equation, M. edulis deposits its shell in isotope equilibrium (d18Ocalcite) with ambient water. Carbon isotopes (d13Ccalcite) from sampled shells were substantially more negative than predicted values, indicating an uptake of metabolic carbon into shell carbonate, and d13Ccalcite disequilibrium increased with increasing salinity. Sampled shells of M. edulis showed no significant trends in d18Ocalcite based on size, cultured growth rates, or geographic collection location, suggesting that vital effects do not affect d18Ocalcite in M. edulis. The broad modern and paleogeographic distribution of this bivalve, its abundance during the Holocene, and the lack of an intraspecies physiologic isotope effect demonstrated here make it an ideal nearshore paleoceanographic proxy throughout much of the North Atlantic Ocean.

(Table S1) Stable oxygen isotopic values and Mg/Ca ratios of Globigerinoides ruber w from sediment core AAS9_21

Variations in sea surface temperature (SST), d18O of sea water (?18Ow), and salinity were reconstructed for the past 68 ka using a sediment core (AAS9/21) from the eastern Arabian Sea (EAS) in order to understand the changes in evaporation and precipitation associated with the monsoon system. The Mg/Ca-derived SST record varies by ~4°C; it shows that marine isotope stage (MIS) 4 was warmer than MIS 3, that the Last Glacial Maximum was 4°C cooler than the present, and that there was a 2°C increase within the Holocene. MIS 4 records higher d18Ow and salinity values than MIS 2, suggesting variable flow of low-salinity Bay of Bengal flow into the EAS during glacial periods. The transition from MIS 4 to MIS 3 was marked with a conspicuous shift from higher to lower d18Ow values, which reflects a decrease in the evaporation-precipitation budget in the EAS, perhaps due to the strengthening of southwest monsoon. Monsoon reconstructions based on d18Ow reveal that monsoon-driven precipitation was higher during MIS 3 and MIS 1 and was lower during MIS 2 and MIS 4. This is consistent with earlier monsoon reconstructions based on upwelling indices from the western Arabian Sea. However, the amplitude of monsoon fluctuations derived through upwelling indices and d18Ow varies significantly, which may indicate spatial variability of monsoon rainfall.

(Table S1) Age determination ofs ediment cores of the Azores Current System

In order to test the sensitivity of marine primary productivity in the midlatitude open ocean North Atlantic to changes in the Atlantic Meridional Overturning Circulation (AMOC), we investigated two spliced sediment cores from a site south of the Azores Islands at the northern rim of the North Atlantic subtropical gyre. For this purpose we analyzed coccolithophore assemblages, diatom abundances, alkenones and conducted X-ray fluorescence (XRF) core scanning. During times of reduced AMOC, especially during Heinrich event 1 (H1) and the Younger Dryas, we observe a strong increase in productivity as evidenced by high coccolith accumulation rates, high alkenone concentrations/accumulation rates, high Ba/Ti-ratios, high abundances of diatoms and low abundances ofF. profunda. The increased productivity is partly caused by a more southern position of the Azores Front (AzF), and hence by a less northward extension of the subtropical gyre, as deduced from high abundances of the temperate coccolithophore species G. muellerae and low abundances of subtropical species (Oolithotus spp., Umbellosphaera spp., Umbilicosphaeraspp.). However, to explain the full range of the observed productivity increase, other factors like increased westerly winds and advection of nutrient-rich surface waters have also to be considered. Because this pattern can also be observed in other sediment cores from the midlatitude North Atlantic, we propose that during times of reduced AMOC there has been a band of strongly increased productivity across the North Atlantic at the northern rim of the contracted subtropical gyre, which partly counteracts the decreased organic carbon pump in the high northern latitudes.