(Table S1) Age determination of sediment core RAPiD-21-03K

A decadal resolution time series of sea surface temperature (SST) spanning the last two millennia is reconstructed by combining a proxy record from a new sediment sequence with previously published data from core MD99-2275, north of Iceland. The alkenone based SST reconstruction is validated with historic observational data and compared to a new similar temporal resolution reconstruction obtained from sediment core RAPiD21-3K, in the subpolar North Atlantic. The two SST paleorecords show consistent multidecadal scale coolings throughout the interval and similar expressions during the contrasted climatic periods 'colloquially known' as the Medieval Climatic Anomaly (MCA) and Little Ice Age (LIA). In order to further understand the temporal and spatial SST variations and investigate the influence of natural forcings on the observed SST changes during the last millennium, we compare our time series to simulations using the Institut Pierre-Simon Laplace IPSLCM4-v2 climate model. This comparison highlights the potential importance of volcanism as a natural forcing driving coherent abrupt cooling events captured in the subpolar North Atlantic records.

(Table S1) Age determination of ODP Site 165-1002

A series of 14C measurements in Ocean Drilling Program cores from the tropical Cariaco Basin, which have been correlated to the annual-layer counted chronology for the Greenland Ice Sheet Project 2 (GISP2) ice core, provides a high-resolution calibration of the radiocarbon time scale back to 50,000 years before the present. Independent radiometric dating of events correlated to GISP2 suggests that the calibration is accurate. Reconstructed 14C activities varied substantially during the last glacial period, including sharp peaks synchronous with the Laschamp and Mono Lake geomagnetic field intensity minimal and cosmogenic nuclide peaks in ice cores and marine sediments. Simulations with a geochemical box model suggest that much of the variability can be explained by geomagnetically modulated changes in 14C production rate together with plausible changes in deep-ocean ventilation and the global carbon cycle during glaciation.

(Table S1) Relative abundances of coccolithophore taxa in surface sediment samples from the Equatorial and Southeastern Pacific

This study aims to contribute to a more detailed knowledge of the biogeography of coccolithophores in the Equatorial and Southeastern Pacific Ocean. Census data of fossil coccoliths are presented in a suite of core-top sediment samples from 15°N to 50.6°S and from 71°W to 93°W. Following standard preparation of smear slides, a total of 19 taxa are recognized in light microscopy and their relative abundances are determined for 134 surface sediment samples. Considering the multivariate character of oceanic conditions and their effects on phytoplankton, a Factor Analysis was performed and three factors were retained. Factor 1, dominated by Florisphaera profunda and Gephyrocapsa oceanica, includes samples located under warm water masses and indicates the occurrence of calcite dissolution in the water column in the area offshore Chile. Factor 2 is related to cold, low-salinity surface-water masses from the Chilean upwelling, and is dominated by Emiliania huxleyi, Gephyrocapsa sp. < 3 µm, Coccolithus pelagicus and Gephyrocapsa muellerae. Factor 3 is linked to more saline, coastal upwelling areas where Calcidiscus leptoporus and Helicosphaera carteri are the dominant species.