Water temperature reconstructions based on paired Uk'37 and Tex86 records, and paired surface and subsurface Mg/Ca records

This dataset contains the collection of available published paired Uk'37 and Tex86 records spanning multi-millennial to multi-million year time scales, as well as a collection of Mg/Ca-derived temperatures measured in parallel on surface and subsurface dwelling foraminifera, both used in the analyses of Ho and Laepple, Nature Geoscience 2016. As the signal-to-noise ratios of proxy-derived Holocene temperatures are relatively low, we selected records that contain at least the last deglaciation (oldest sample >18kyr BP).

Age models, accumulation rates, and grain size distributions of sediment cores from the northern subtropical Atlantic

The terrigenous sediment proportion of the deep sea sediments from off Northwest Africa has been studied in order to distinguish between the aeolian and the fluvial sediment supply. The present and fossil Saharan dust trajectories were recognized from the distribution patterns of the aeolian sediment. The following timeslices have been investigated: Present, 6,000, 12,000 and 18,000 y. B. P. Furthermore, the quantity of dust deposited off the Saharan coast has been estimated. For this purpose, 80 surface sediment samples and 34 sediment cores have been analysed. The stratigraphy of the cores has been achieved from oxygen isotopic curves, 14C-dating, foraminiferal transfer temperatures, and carbonate contents. Silt sized biogenic opal generally accounts for less than 2 % of the total insoluble sediment proportion. Only under productive upwelling waters and off river mouths, the opal proportion exceeds 2 % significantly. The modern terrigenous sediment from off the Saharan coast is generally characterized by intensely stained quartz grains. They indicate an origin from southern Saharan and Sahelian laterites, and a zonal aeolian transport in midtropospheric levels, between 1.5 an 5.5 km, by 'Harmattan' Winds. The dust particles follow large outbreaks of Saharan air across the African coast between 15° and 21° N. Their trajectories are centered at about 18° N and continue further into a clockwise gyre situated south of the Canary Islands. This course is indicated by a sickle-shaped tongue of coarser grain sizes in the deep-sea sediment. Such loess-sized terrigenous particles only settle within a zone extending to 700 km offshore. Fine silt and clay sized particles, with grain sizes smaller than 10- 15 µm, drift still further west and can be traced up to more than 4,000 km distance from their source areas. Additional terrigenous silt which is poor in stained quartz occurs within a narrow zone off the western Sahara between 20° and 27° N only. It depicts the present dust supply by the trade winds close to the surface. The dust load originates from the northwestern Sahara, the Atlas Mountains and coastal areas, which contain a particularly low amount of stained quartz. The distribution pattern of these pale quartz sediments reveals a SSW-dispersal of dust being consistent with the present trade wind direction from the NNE. In comparison to the sediments from off the Sahara and the deeper subtropical Atlantic, the sediments off river mouths, in particular off the Senegal river, are characterized by an additional input of fine grained terrigenous particles (< 6 µm). This is due to fluvial suspension load. The fluvial discharge leads to a relative excess of fine grained particles and is observed in a correlation diagram of the modal grain sizes of terrigenous silt with the proportion of fine fraction (< 6 µm). The aeolian sediment contribution by the Harmattan Winds strongly decreased during the Climatic Optimum at 6,000 y. B. P. The dust discharge of the trade winds is hardly detectable in the deep-sea sediments. This probably indicates a weakened atmospheric circulation. In contrast, the fluvial sediment supply reached a maximum, and can be traced to beyond Cape Blanc. Thus, the Saharan climate was more humid at 6,000 y B. P. A latitudinal shift of the Harmattan driven dust outbreaks cannot be observed. Also during the Glacial, 18,000 y. B. P., Harmattan dust transport crossed the African coast at latitudes of 15°-20° N. Its sediment load increased intensively, and markedly coarser grains spread further into the Atlantic Ocean. An expanded zone of pale-quart sediments indicates an enhanced dust supply by the trade winds blowing from the NE. No synglacial fluvial sediment contribution can be recognized between 12° and 30° N. This indicates a dry glacial climate and a strengthened stmospheric circulation over the Sahelian and Saharan region. The climatic transition pahes, at 12, 000 y. B. P., between the last Glacial and the Intergalcial, which is compareable to the Alerod in Europe, is characterized by an intermediate supply of terrigenous particles. The Harmattan dust transport wa weaker than during the Glacial. The northeasterly trade winds were still intensive. River supply reached a first postglacial maximum seaward of the Senegal river mouth. This indicates increasing humidity over the southern Sahara and a weaker atmospheric circulation as compared to the glacial. The accumulation rates of the terrigenous silt proportion (> 6 µm) decrcase exponentially with increasing distance from the Saharan coast. Those of the terrigenous fine fraction (< 6 µm) follow the same trend and show almost similar gradients. Accordingly, also the terrigenous fine fraction is believed to result predominantly from aeolian transport. In the Atlantic deep-sea sediments, the annual terrigenous sediment accumulation has fluctuated, from about 60 million tons p. a. during the Late Glacial (13,500-18,000 y. B. P, aeolian supply only) to about 33 million tons p. a. during the Holocene Climatic Optimum (6,000-9,000 y. B. P, mainly fluvial supply), when the river supply has reached a maximum, and to about 45 million tons p. a. during the last 4,000 years B. P. (fluvial supply only south of 18° N).

Surface wind stress modeled for the Panama Seaway with the Community Climate System Model 3

The early Pliocene warm phase was characterized by high sea surface temperatures and a deep thermocline in the eastern equatorial Pacific. A new hypothesis suggests that the progressive closure of the Panamanian seaway contributed substantially to the termination of this zonally symmetric state in the equatorial Pacific. According to this hypothesis, intensification of the Atlantic meridional overturning circulation (AMOC) - induced by the closure of the gateway - was the principal cause of equatorial Pacific thermocline shoaling during the Pliocene. In this study, twelve Panama seaway sensitivity experiments from eight ocean/climate models of different complexity are analyzed to examine the effect of an open gateway on AMOC strength and thermocline depth. All models show an eastward Panamanian net throughflow, leading to a reduction in AMOC strength compared to the corresponding closed-Panama case. In those models that do not include a dynamic atmosphere, deepening of the equatorial Pacific thermocline appears to scale almost linearly with the throughflow-induced reduction in AMOC strength. Models with dynamic atmosphere do not follow this simple relation. There are indications that in four out of five models equatorial wind-stress anomalies amplify the tropical Pacific thermocline deepening. In summary, the models provide strong support for the hypothesized relationship between Panama closure and equatorial Pacific thermocline shoaling.

Analyses of benthic and planktonic foraminifera during Marine Isotope Stage 5 of ODP Hole 172-1058C

Subtropical Gyres are an important constituent of the ocean-atmosphere system due to their capacity to store vast amounts of warm and saline waters. Here we decipher the sensitivity of the (sub)surface North Atlantic Subtropical Gyre with respect to orbital and millennial scale climate variability between ~140 and 70 ka, Marine Isotope Stage (MIS) 5. Using (isotope)geochemical proxy data from surface and thermocline dwelling foraminifers from Blake Ridge off the west coast of North America (ODP Site 1058) we show that the oceanographic development at subsurface (thermocline) level is substantially different from the surface ocean. Most notably, surface temperatures and salinities peak during the penultimate deglaciation (Termination II) and early MIS 5e, implying that subtropical surface ocean heat and salt accumulation might have resulted from a sluggish northward heat transport. In contrast, maximum thermocline temperatures are reached during late MIS 5e when surface temperatures are already declining. We argue that the subsurface warming originated from intensified Ekman downwelling in the Subtropical Gyre due to enhanced wind stress. During MIS 5a-d a tight interplay of the subtropical upper ocean hydrography to high latitude millennial-scale cold events can be observed. At Blake Ridge, the most pronounced of these high latitude cold events are related to surface warming and salt accumulation in the (sub)surface. Similar to Termination II, heat accumulated in the Subtropical Gyre probably due to a reduced Atlantic Meridional Overturning Circulation. Additionally, a southward shift and intensification of the subtropical wind belts lead to a decrease of on-site precipitation and enhanced evaporation, coupled to intensified gyre circulation. Subsequently, the northward advection of these warm and saline water likely contributed to the fast resumption of the overturning circulation at the end of these high latitude cold events.

Diatom oxygen and silicon isotopes data from IODP Hole 323-U1341B in the Bering Sea through the late Pliocene and over the onset of major Northern Hemisphere Glaciation

The Bering Sea is one of the most biologically productive regions in the marine system and plays a key role in regulating the flow of waters to the Arctic Ocean and into the subarctic North Pacific Ocean. Cores from Integrated Ocean Drilling Program (IODP) Expedition 323 to the Bering Sea provide the first opportunity to obtain reconstructions from the region that extend back to the Pliocene. Previous research at Bowers Ridge, south Bering Sea, has revealed stable levels of siliceous productivity over the onset of major Northern Hemisphere Glaciation (NHG) (circa 2.85-2.73 Ma). However, diatom silica isotope records of oxygen (d18Odiatom) and silicon (d30Sidiatom) presented here demonstrate that this interval was associated with a progressive increase in the supply of silicic acid to the region, superimposed on shift to a more dynamic environment characterized by colder temperatures and increased sea ice. This concluded at 2.58 Ma with a sharp increase in diatom productivity, further increases in photic zone nutrient availability and a permanent shift to colder sea surface conditions. These transitions are suggested to reflect a gradually more intense nutrient leakage from the subarctic northwest Pacific Ocean, with increases in productivity further aided by increased sea ice- and wind-driven mixing in the Bering Sea. In suggesting a linkage in biogeochemical cycling between the south Bering Sea and subarctic Northwest Pacific Ocean, mainly via the Kamchatka Strait, this work highlights the need to consider the interconnectivity of these two systems when future reconstructions are carried out in the region.

Planktonic foraminiferal biostratigraphy and stable carbon isotope ratios of late middle Eocene sediments from Blake Plateau, West Atlantic Ocean

The muricate planktonic foraminiferal genera Morozovella and Acarinina were abundant and diverse during the upper Palaeocene to middle Eocene and dominated the tropical and subtropical assemblages. A significant biotic turnover in planktonic foraminifera occurred in the latest middle Eocene with a notable reduction in the acarininid lineage and the extinction of the morozovellids. These genera are extensively employed as palaeoclimatic and biostratigraphic markers and, therefore, this turnover episode is an important event in the record of the Cenozoic planktonic foraminifera. Sediments from the western North Atlantic (Ocean Drilling Program Site 1052) were examined in order to investigate these extinction events, in terms of both timing and mechanisms. Biostratigraphic events of the middle and late Eocene have been examined with a sampling resoluti on of approximately 3 kyr. These have been calibrated to the magneto- and astrochronology to accurately define the timing of key biostratigraphic events, particularly the extinction of Morozovella spinulosa which is a distinct biomarker for late middle Eocene sediments. High-resolution biostratigraphy reveals that the extinctions in the muricate group occurred in a stepwise form. The large acarininids (Acarinina praetopilensis) terminate 10 kyr prior to the extinction of M. spinulosa and small acarininids (Acarinina medizzai and Acarinina echinata) continue into the upper Eocene. High-resolution stable isotope analyses have been conducted on planktonic and benthic foraminifera from the western North Atlantic to reconstruct sea surface temperatures (SSTs) and deep water temperatures and the structure of the water column around this major biotic turnover. Whilst the extinctions of M. spinulosa and A. praetopilensis occur during a long-term cooling trend, the biotic turnover in the muricate group does not appear to be related to significant climatic change. Sea surface temperatures decrease slowly prior to the extinction events, and there is no evidence for a large-temperature shift associated with the faunal changes. The turnover event was therefore probably related to the increased surface water productivity and the deterioration of photosymbiotic partnerships with algae.

Diatom assemblages, HBIs and TEXL86 in sediment core NBP99-03_10

The West Antarctic ice sheet is particularly sensitive to global warming and its evolution and impact on global climate over the next few decades remains difficult to predict. In this context, investigating past sea ice conditions around Antarctica is of primary importance. Here, we document changes in sea ice presence, upper water column temperatures (0-200 m) and primary productivity over the last 9000 yr BP (before present) in the western Antarctic Peninsula (WAP) margin from a sedimentary core collected in the Palmer Deep Basin. Employing a multi-proxy approach, based on the combination of two biomarkers proxies (highly branched isoprenoid (HBI) alkenes for sea ice and TEXL86 for temperature) and micropaleontological data (diatom assemblages), we derived new Holocene records of sea ice conditions and upper water column temperatures. The early Holocene (9000-7000 yr BP) was characterized by a cooling phase with a short sea ice season. During the mid-Holocene (~7000-3800 yr BP), local climate evolved towards slightly colder conditions and a prominent extension of the sea ice season occurred, promoting a favorable environment for intensive diatom growth. The late Holocene (the last ~2100 yr) was characterized by warmer temperatures and increased sea ice presence, accompanied by reduced local primary productivity, likely in response to a shorter growing season compared to the early or mid-Holocene. The gradual increase in annual sea ice duration over the last 7000 yr might have been influenced by decreasing mean annual and spring insolation, despite increasing summer insolation. We postulate that, in addition to precessional changes in insolation, seasonal variability, via changes in the strength of the circumpolar Westerlies and upwelling activity, was further amplified by the increasing frequency/amplitude of the El Nino-Southern Oscillation (ENSO). However, between 3800 and 2100 yr BP, the lack of correlation between ENSO and climate variability in the WAP suggests that other climatic factors might have been more important in controlling WAP climate at this time.

Daily polynya area and ice production in the North Water polynya region

The North Water (NOW) Polynya is a regularly-forming area of open-water and thin-ice, located between northwestern Greenland and Ellesmere Island (Canada) at the northern tip of Baffin Bay. Due to its large spatial extent, it is of high importance for a variety of physical and biological processes, especially in wintertime. Here, we present a long-term remote sensing study for the winter seasons 1978/1979 to 2014/2015. Polynya characteristics are inferred from (1) sea ice concentrations and brightness temperatures from passive microwave satellite sensors (Advanced Microwave Scanning Radiometer (AMSR-E and AMSR2), Scanning Multichannel Microwave Radiometer (SMMR), Special Sensor Microwave Imager/Sounder (SSM/I-SSMIS)) and (2) thin-ice thickness distributions, which are calculated using MODIS ice-surface temperatures and European Center for Medium-Range Weather Forecasts (ECMWF) atmospheric reanalysis data in a 1D thermodynamic energy-balance model. Daily ice production rates are retrieved for each winter season from 2002/2003 to 2014/2015, assuming that all heat loss at the ice surface is balanced by ice growth. Two different cloud-cover correction schemes are applied on daily polynya area and ice production values to account for cloud gaps in the MODIS composites. Our results indicate that the NOW polynya experienced significant seasonal changes over the last three decades considering the overall frequency of polynya occurrences, as well as their spatial extent. In the 1980s, there were prolonged periods of a more or less closed ice cover in northern Baffin Bay in winter. This changed towards an average opening on more than 85% of the days between November and March during the last decade. Noticeably, the sea ice cover in the NOW polynya region shows signs of a later-appearing fall freeze-up, starting in the late 1990s. Different methods to obtain daily polynya area using passive microwave AMSR-E/AMSR2 data and SSM/I-SSMIS data were applied. A comparison with MODIS data (thin-ice thickness < 20 cm) shows that the wintertime polynya area estimates derived by MODIS are about 30 to 40% higher than those derived using the polynya signature simulation method (PSSM) with AMSR-E data. In turn, the difference in polynya area between PSSM and a sea ice concentration (SIC) threshold of 70% is fairly low (approximately 10%) when applied to AMSR-E data. For the coarse-resolution SSM/I-SSMIS data, this difference is much larger, particularly in November and December. Instead of a sea ice concentration threshold, the PSSM method should be used for SSM/I-SSMIS data. Depending on the type of cloud-cover correction, the calculated ice production based on MODIS data reaches an average value of 264.4 ± 65.1 km**3 to 275.7 ± 67.4 km**3 (2002/2003 to 2014/2015) and shows a high interannual variability. Our achieved long-term results underline the major importance of the NOW polynya considering its influence on Arctic ice production and associated atmosphere/ocean processes.

Bulk geochemical analysis and planktonic foraminifera of sediment core GeoB7112-5

A multivariable approach utilising bulk sediment, planktonic Foraminifera and siliceous phytoplankton has been used to reconstruct rapid variations in palaeoproductivity in the Peru-Chile Current System off northern Chile for the past 19000 cal. yr. During the early deglaciation (19000-16000 cal. yr BP), our data point to strongest upwelling intensity and highest productivity of the past 19 000 cal. yr. The late deglaciation (16000-13000 cal. yr BP) is characterised by a major change in the oceanographic setting, warmer water masses and weaker upwelling at the study site. Lowest productivity and weakest upwelling intensity are observed from the early to the middle Holocene (13000-4000 cal. yr BP), and the beginning of the late Holocene (<4000 cal. yr BP) is marked by increasing productivity, mainly driven by silicate-producing organisms. Changes in the productivity and upwelling intensity in our record may have resulted from a large-scale compression and/or displacement of the South Pacific subtropical gyre during more productive periods, in line with a northward extension of the Antarctic Circumpolar Current and increased advection of Antarctic water masses with the Peru-Chile Current. The corresponding increase in hemispheric thermal gradient and wind stress induced stronger upwelling. During the periods of lower productivity, this scenario probably reversed.

Benthic and planktic foraminiferal stable isotope records combined with Nd isotope ratios of ODP Hole 198-1210B

The Pacific Ocean is the largest water body on Earth, and circulation in the Pacific contributed significantly to climate evolution in the latest Cretaceous, the culmination of a period of long-term cooling. Here, we present new high-resolution late Campanian to Maastrichtian benthic and planktic foraminiferal stable isotope data and a neodymium (Nd) isotope record obtained from sedimentary ferromanganese oxide coatings of Ocean Drilling Program Hole 1210B from the tropical Pacific Ocean (Shatsky Rise). These new records resolve 13 million years in the latest Cretaceous, providing insights into changes in surface and bottom water temperatures and source regions of deep to intermediate waters covering the carbon isotope excursions of the Campanian-Maastrichtian Boundary Event (CMBE) and the Mid-Maastrichtian event (MME). Our new benthic foraminiferal d18O and Nd isotope records together with published Nd isotope data show markedly parallel trends across the studied interval over a broad range of bathyal to abyssal water depths interpreted to reflect changes in the intensity of deep-ocean circulation in the tropical Pacific. In particular, we observe a three-million-year-long period of cooler conditions in the early Maastrichtian (72.5 to 69.5 Ma) when a concomitant change toward less radiogenic seawater Nd isotope signatures probably marks a period of enhanced admixture and northward flow of deep waters with Southern Ocean provenance. We suggest this change to have been triggered by intensified formation and convection of deep waters in the high southern latitudes, a process that weakened during the MME (69.5 to 68.5 Ma). The early Maastrichtian cold interval is closely related to the negative and positive carbon isotope trends of the CMBE and MME. The millions-of-years long duration of these carbon cycle perturbations suggests a tectonic forcing of climatic cooling, possibly related to changes in ocean basin geometry and bathymetry.