(Table S1) Age determination of sediment cores from the Labrador Sea

High resolution flow speed reconstructions of two core sites located on Gardar Drift in the northeast Atlantic Basin and Orphan Knoll in the northwest Atlantic Basin reveal a long-term decrease in flow speed of Northeast Atlantic Deep Water (NEADW) after 6,500 years. Benthic foraminiferal oxygen isotopes of sites currently bathed in NEADW show a 0.2per mil depletion after 6,500 years, shortly after the start of the development of a carbon isotope gradient between NEADW and Norwegian Sea Deep Water. We consider these changes in near-bottom flow vigor and benthic foraminiferal isotope records to mark a significant reorganization of the Holocene deep ocean circulation, and attribute the changes to a weakening of NEADW flow during the mid to late Holocene that allowed the shoaling of Lower Deep Water and deeper eastward advection of Labrador Sea Water into the northeast Atlantic Basin.

(Table S1) TEX86 and BIT record from sediment core PBBC-1

A molecular organic geochemical proxy (TEX86) for sea surface temperature (SST) is compared with a foraminifera-based SST proxy (Mg/Ca) in a decadal-resolution marine sedimentary record spanning the last 1000 years from the Gulf of Mexico. We assess the relative strengths of the organic and inorganic paleoceanographic techniques for reconstructing high-resolution SST variability during recent climate events, including the Little Ice Age (LIA) and the Medieval Warm Period (MWP). SST estimates based on the molecular organic proxy TEX86 show a similar magnitude and pattern of SST variability to foraminiferal Mg/Ca-SST estimates but with some important differences. For instance, both proxies show a cooling (1°C-2°C) of Gulf of Mexico SSTs during the LIA. During the MWP, however, Mg/Ca-SSTs are similar to near-modern SSTs, while TEX86 indicates SSTs that were cooler than modern. Using the respective SST calibrations for each proxy results in TEX86-SST estimates that are 2°C-4°C warmer than Mg/Ca-SST throughout the 1000 year record. We interpret the TEX86-SST as a summer-weighted SST signal from the upper mixed layer, whereas the Mg/Ca-SST better reflects the mean annual SST. Downcore differences in the SST estimates between the two proxies (DeltaT = TEX86 - Mg/Ca) are interpreted in the context of varying seasonality and/or changing water column temperature gradients.

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.