Supplementary material: XRF core scanning data, CaCO3 values, foraminiferal carbon and oxygen isotope data, Mg/Ca data, clay mineralogy and helium isotope results

The Paleocene-Eocene Thermal Maximum (PETM; ~56 Ma) is associated with abrupt climate change, carbon cycle perturbation, ocean acidification, as well as biogeographic shifts in marine and terrestrial biota that were largely reversed as the climatic transient waned. We report a clear exception to the behavior of the PETM as a reversing climatic transient in the eastern North Atlantic (Deep-Sea Drilling Project Site 401, Bay of Biscay) where the PETM initiates a greatly prolonged environmental change compared to other places on Earth where records exist. The observed environmental perturbation extended well past the d13C recovery phase and up to 650 kyr after the PETM onset according to our extraterrestrial 3He-based age-model. We observe a strong decoupling of planktic foraminiferal d18O and Mg/Ca values during the PETM d13C recovery phase, which in combination with results from helium isotopes and clay mineralogy, suggests that the PETM triggered a hydrologic change in western Europe that increased freshwater flux and the delivery of weathering products to the eastern North Atlantic. This state change persisted long after the carbon-cycle perturbation had stopped. We hypothesize that either long-lived continental drainage patterns were altered by enhanced hydrological cycling induced by the PETM, or alternatively that the climate system in the hinterland area of Site 401 was forced into a new climate state that was not easily reversed in the aftermath of the PETM.

Measurements of formation, degradation, and sinking velocity of diatom aggregates from laboratory studies performed at 15 °C and 4 °C

Most deep ocean carbon flux profiles show low and almost constant fluxes of particulate organic carbon (POC) in the deep ocean. However, the reason for the non-changing POC fluxes at depths is unknown. This study presents direct measurements of formation, degradation, and sinking velocity of diatom aggregates from laboratory studies performed at 15 °C and 4 °C during a three-week experiment. The average carbon-specific respiration rate during the experiment was 0.12 ± 0.03 at 15 °C, and decreased 3.5-fold when the temperature was lowered to 4 °C. No direct influence of temperature on aggregate sinking speed was observed. Using the remineralisation rate measured at 4 °C and an average particle sinking speed of 150 m d**-1, calculated carbon fluxes were similar to those collected in deep ocean sediment traps from a global data set, indicating that temperature plays a major role for deep ocean fluxes of POC.

Data base on heterotrophic dinoflagellates grazing and growth rate as a function of size, prey size and type compiled from the literature

Microzooplankton (the 20 to 200 µm size class of zooplankton) is recognised as an important part of marine pelagic ecosystems. In terms of biomass and abundance heterotrophic dinoflagellates are one of the important groups of organism in microzooplankton. However, their rates - grazing and growth - , feeding behaviour and prey preferences are poorly known and understood. A set of data was assembled in order to derive a better understanding of heterotrophic dinoflagellates rates, in response to parameters such as prey concentration, prey type (size and species), temperature and their own size. With these objectives, literature was searched for laboratory experiments with information on one or more of these parameters effect studied. The criteria for selection and inclusion in the database included: (i) controlled laboratory experiment with a known dinoflagellate feeding on a known prey; (ii) presence of ancillary information about experimental conditions, used organisms - cell volume, cell dimensions, and carbon content. Rates and ancillary information were measured in units that meet the experimenter need, creating a need to harmonize the data units after collection. In addition different units can link to different mechanisms (carbon to nutritive quality of the prey, volume to size limits). As a result, grazing rates are thus available as pg C dinoflagellate-1 h-1, µm3 dinoflagellate-1 h-1 and prey cell dinoflagellate-1 h-1; clearance rate was calculated if not given and growth rate is expressed as the growth rate per day.

Processed data from CTD casts and fluorospectrometer measurements from the SONNE SO226/3 cruise track

The Sonne transit cruise SO226-3 DipFIP took place from March 4th (Wellington, New Zealand) to March 28th (Kaohsiung, Taiwan) in 2013. CTD data for 16 stations along the cruise track were recorded using the onboard a SEABIRD SBE 9 plus CTD down to depth of 800m. Obtained hydrographic data were binned to 1 m intervals with the available SBE software. Obvious outliers in the readings of the oxygen sensor close to the sea surface have been manually removed. Fluorospectrometer (bbe Moldaenke) pigment data measured for 24 depth intervals are available for 10 stations. Measurements were conducted in the shipboard laboratory on water samples from the CTD rosette. Data are averages of a least 30 readings per sample.

Radiocarbon ages and diatom isotope data for Lake El'gygytgyn sediment core Lz1029

Determining the response of sites within the Arctic Circle to long-term climatic change remains an essential pre-requisite for assessing the susceptibility of these regions to future global warming and Arctic amplification. To date, existing records from North East Russia have demonstrated significant spatial variability across the region during the late Quaternary. Here we present diatom d18O and d30Si data from Lake El'gygytgyn, Russia, and suggest environmental changes that would have impacted across West Beringia from the Last Glacial Maximum to the modern day. In combination with other records, the results raise the potential for climatic teleconnections to exist between the region and sites in the North Atlantic. The presence of a series of 2-3 per mil decreases in d18Odiatom during both the Last Glacial and the Holocene indicates the sensitivity of the region to perturbations in the global climate system. Evidence of an unusually long Holocene thermal maximum from 11.4 ka BP to 7.6 ka BP is followed by a cooling trend through the remainder of the Holocene in response to changes in solar insolation. This is culminated over the last 900 years by a significant decrease in d18Odiatom of 2.3 per mil, which may be related to a strengthening and easterly shift of the Aleutian Low in addition to possible changes in precipitation seasonality.