(Fig. 4) Multiproxy data set of box core MSM5/5-712-1

A multiproxy data set of an AMS radiocarbon dated 46 cm long sediment core from the continental margin off western Svalbard reveals multidecadal climatic variability during the past two millennia. Investigation of planktic and benthic stable isotopes, planktic foraminiferal fauna, and lithogenic parameters aims to unveil the Atlantic Water advection to the eastern Fram Strait by intensity, temperatures, and salinities. Atlantic Water has been continuously present at the site over the last 2,000 years. Superimposed on the increase in sea ice/icebergs, a strengthened intensity of Atlantic Water inflow and seasonal ice-free conditions were detected at ~ 1000 to 1200 AD, during the well-known Medieval Climate Anomaly (MCA). However, temperatures of the MCA never exceeded those of the 20th century. Since ~ 1400 AD significantly higher portions of ice rafted debris and high planktic foraminifer fluxes suggest that the site was located in the region of a seasonal highly fluctuating sea ice margin. A sharp reduction in planktic foraminifer fluxes around 800 AD and after 1730 AD indicates cool summer conditions with major influence of sea ice/icebergs. High amounts of the subpolar planktic foraminifer species Turborotalia quinqueloba in size fraction 150-250 µm indicate strengthened Atlantic Water inflow to the eastern Fram Strait already after ~ 1860 AD. Nevertheless surface conditions stayed cold well into the 20th century indicated by low planktic foraminiferal fluxes. Most likely at the beginning of the 20th century, cold conditions of the terminating Little Ice Age period persisted at the surface whereas warm and saline Atlantic Water already strengthened, hereby subsiding below the cold upper mixed layer. Surface sediments with high abundances of subpolar planktic foraminifers indicate a strong inflow of Atlantic Water providing seasonal ice-free conditions in the eastern Fram Strait during the last few decades.

(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.

X-ray fluorescent (XRF) ratios and diatom groups of sediment core GeoB13813-4

The continental shelf adjacent to the Río de la Plata (RdlP) exhibits extremely complex hydrographic and ecological characteristics which are of great socioeconomic importance. Since the long-term environmental variations related to the atmospheric (wind fields), hydrologic (freshwater plume), and oceanographic (currents and fronts) regimes are little known, the aim of this study is to reconstruct the changes in the terrigenous input into the inner continental shelf during the late Holocene period (associated with the RdlP sediment discharge) and to unravel the climatic forcing mechanisms behind them. To achieve this, we retrieved a 10 m long sediment core from the RdlP mud depocenter at 57 m water depth (GeoB 13813-4). The radiocarbon age control indicated an extremely high sedimentation rate of 0.8 cm per year, encompassing the past 1200 years (AD 750-2000). We used element ratios (Ti / Ca, Fe / Ca, Ti / Al, Fe / K) as regional proxies for the fluvial input signal and the variations in relative abundance of salinity-indicative diatom groups (freshwater versus marine-brackish) to assess the variability in terrigenous freshwater and sediment discharges. Ti / Ca, Fe / Ca, Ti / Al, Fe / K and the freshwater diatom group showed the lowest values between AD 850 and 1300, while the highest values occurred between AD 1300 and 1850. The variations in the sedimentary record can be attributed to the Medieval Climatic Anomaly (MCA) and the Little Ice Age (LIA), both of which had a significant impact on rainfall and wind patterns over the region. During the MCA, a weakening of the South American summer monsoon system (SAMS) and the South Atlantic Convergence Zone (SACZ), could explain the lowest element ratios (indicative of a lower terrigenous input) and a marine-dominated diatom record, both indicative of a reduced RdlP freshwater plume. In contrast, during the LIA, a strengthening of SAMS and SACZ may have led to an expansion of the RdlP river plume to the far north, as indicated by higher element ratios and a marked freshwater diatom signal. Furthermore, a possible multidecadal oscillation probably associated with Atlantic Multidecadal Oscillation (AMO) since AD 1300 reflects the variability in both the SAMS and SACZ systems.

High resolution stable isotope and Mg/Ca record of sediment core MD02-2550C2

An early Holocene record from the Gulf of Mexico (GOM) reveals climatic and hydrologic changes during the interval from 10.5 to 7 thousand calendar years before present from paired analyses of Mg/Ca and d18O on foraminiferal calcite. The sea surface temperature record based on foraminiferal Mg/Ca contains six oscillations and an overall ~1.5°C warming that appears to be similar to the September-March insolation difference. The d18O of seawater in the GOM (d18OGOM) record contains six oscillations, including a -0.8 per mil excursion that may be associated with the "8.2 ka climate event" or a broader climate anomaly. Faunal census records from three GOM cores exhibit similar changes, suggesting subcentennial-scale variability in the incursions of Caribbean waters into the GOM. Overall, our results provide evidence that the subtropics were characterized by decadal- to centennial-scale climatic and hydrologic variability during the early Holocene.

Radiocarbon age, Mg/Ca and d18O measurements on planktonic foraminifera of sediment core GeoB12605-3

The sea surface temperature (SST) of the tropical Indian Ocean is a major component of global climate teleconnections. While the Holocene SST history is documented for regions affected by the Indian and Arabian monsoons, data from the near-equatorial western Indian Ocean are sparse. Reconstructing past zonal and meridional SST gradients requires additional information on past temperatures from the western boundary current region. We present a unique record of Holocene SST and thermocline depth variations in the tropical western Indian Ocean as documented in foraminiferal Mg/Ca ratios and d18O from a sediment core off northern Tanzania. For Mg/Ca and thermocline d18O, most variance is concentrated in the centennial to bicentennial periodicity band. On the millennial time scale, an early to mid-Holocene (~7.8-5.6 ka) warm phase is followed by a temperature drop by up to 2°C, leading to a mid-Holocene cool interval (5.6-4.2 ka). The shift is accompanied by an initial reduction in the difference between surface and thermocline foraminiferal d18O, consistent with the thickening of the mixed layer and suggestions of a strengthened Walker circulation. However, we cannot confirm the expected enhanced zonal SST gradient, as the cooling of similar magnitude had previously been found in SSTs from the upwelling region off Sumatra and in Flores air temperatures. The SST pattern probably reflects the tropical Indian Ocean expression of a large-scale climate anomaly rather than a positive Indian Ocean Dipole-like mean state.

Model results in NetCDF format of North African vegetation-precipitation in early and mid-Holocene climate simulated by CCSM3-DGVM//

The present study analyses the sign, strength, and working mechanism of the vegetation-precipitation feedback over North Africa in middle (6 ka BP) and early Holocene (9 ka BP) simulations using the comprehensive coupled climate-vegetation model CCSM3-DGVM (Community Climate System Model version 3 and a dynamic global vegetation model). The coupled model simulates enhanced summer rainfall and a northward migration of the West African monsoon trough along with an expansion of the vegetation cover for the early and middle Holocene compared to the pre-industrial period. It is shown that dynamic vegetation enhances the orbitally triggered summer precipitation anomaly by approximately 20% in the Sahara-Sahel region (10-25° N, 20° W-30° E) in both the early and mid-Holocene experiments compared to their fixed-vegetation counterparts. The primary vegetation-rainfall feedback identified here operates through surface latent heat flux anomalies by canopy evaporation and transpiration and their effect on the mid-tropospheric African easterly jet, whereas the effects of vegetation changes on surface albedo and local water recycling play a negligible role. Even though CCSM3-DGVM simulates a positive vegetation-precipitation feedback in the North African region, this feedback is not strong enough to produce multiple equilibrium climate-ecosystem states on a regional scale.

Results of calculated surface air temperature anomaly of COSMOS Pliocene experiment 2 in NetCDF format

In this manuscript we describe the experimental procedure employed at the Alfred Wegener Institute in Germany in the preparation of the simulations for the Pliocene Model Intercomparison Project (PlioMIP). We present a description of the utilized Community Earth System Models (COSMOS, version: COSMOS-landveg r2413, 2009) and document the procedures that we applied to transfer the Pliocene Research, Interpretation and Synoptic Mapping (PRISM) Project mid-Pliocene reconstruction into model forcing fields. The model setup and spin-up procedure are described for both the paleo- and preindustrial (PI) time slices of PlioMIP experiments 1 and 2, and general results that depict the performance of our model setup for mid-Pliocene conditions are presented. The mid-Pliocene, as simulated with our COSMOS setup and PRISM boundary conditions, is both warmer and wetter in the global mean than the PI. The globally averaged annual mean surface air temperature in the mid-Pliocene standalone atmosphere (fully coupled atmosphere-ocean) simulation is 17.35 °C (17.82 °C), which implies a warming of 2.23 °C (3.40 °C) relative to the respective PI control simulation.

Seafloor spreading magnetic anomaly isochron map compilation for the Weddell Sea and Scotia Sea

The deep sea sedimentary record is an archive of the pre-glacial to glacial development of Antarctica and changes in climate, tectonics and ocean circulation. Identification of the pre-glacial, transitional and full glacial components in the sedimentary record is necessary for ice sheet reconstruction and to build circum-Antarctic sediment thickness grids for past topography and bathymetry reconstructions, which constrain paleoclimate models. A ~3300 km long Weddell Sea to Scotia Sea transect consisting of multichannel seismic reflection data from various organisations, were used to interpret new horizons to define the initial basin-wide seismostratigraphy and to identify the pre-glacial to glacial components. We mapped seven main units of which three are in the inferred Cretaceous-Paleocene pre-glacial regime, one in the Eocene-Oligocene transitional regime and three units in the Miocene-Pleistocene full glacial climate regime. Sparse borehole data from ODP leg 113 and SHALDRIL constrain the ages of the upper three units. Compiled seafloor spreading magnetic anomalies constrain the basement ages and the hypothetical age model. In many cases, the new horizons and stratigraphy contradict the interpretations in local studies. Each seismic sedimentary unit and its associated base horizon are continuous and traceable for the entire transect length, but reflect a lateral change in age whilst representing the same deposition process. The up to 1240 m thick pre-glacial seismic units form a mound in the central Weddell Sea basin and, in conjunction with the eroded flank geometry, support the interpretation of a Cretaceous proto-Weddell Gyre. The base reflector of the transitional seismic unit, which marks the initial ice sheet advances to the outer shelf, has a lateral model age of 26.6-15.5 Ma from southeast to northwest. The Pliocene-Pleistocene glacial deposits reveals lower sedimentations rates, indicating a reduced sediment supply. Sedimentation rates for the pre-glacial, transitional and full glacial components are highest around the Antarctic Peninsula, indicating higher erosion and sediment supply of a younger basement. We interpret an Eocene East Antarctic Ice Sheet expansion, Oligocene grounding of the West Antarctic Ice Sheet and Early Miocene grounding of the Antarctic Peninsula Ice Sheet.