Sea surface temperature reconstructions of sediment cores from the Atlantic sector of the Southern Ocean

Late Quaternary summer sea surface temperatures (SSTs) have been derived from radiolarian assemblages in the East Atlantic sector of the Southern Ocean. In the subantarctic and the polar frontal zone, glacial SSTs (oxygen isotope stages 2, 4, 6, and 8) were 3°-5°C cooler than today, indicating northward displacements of the isotherms about 2°-4° of latitudes. During interglacials, SSTs almost reached modern levels (oxygen isotope stages 7 and 9) or exceeded them by 2°-3°C (oxygen isotope stages 1 and 5.5). In the subantarctic Atlantic Ocean, changes in SST and calcium carbonate content of the sediment precede variations in global ice volume in the range of the main Milankovitch frequencies. Comparisons with the timing of North Atlantic Deep Water (NADW) proxy records suggests that this early response in the subantarctic Atlantic Ocean is not triggered by the flux of NADW to the Southern Ocean.

Sea-surface temperature reconstructions for the Quaternary western Atlantic

Paleotemperature estimates calculated by the SIMMAX Modern Analog Technique are presented for two gravity cores from the Rio Grande Rise, one from the Brazil Slope, and one from the Ceara Rise. The estimates are based on comparisons between modern and fossil planktonic foraminiferal assemblages and were carried out on samples from Quaternary sediments. Estimated warm-season temperatures from the Rio Grande Rise (at approx. 30° S) range from around 19°C to 24°C, with some coincidence of warm peaks with interglacial stages. The temperature estimates (also warm-season) from the more tropical Brazil Slope (at approx. 8° S) and Ceara Rise (at approx. 4° N) cores are more stable, remaining between 26°C and 28°C throughout most of their lengths. This fairly stable situation in the tropical western Atlantic is interrupted in oxygen isotope stage 6 by a significant drop of 2-3°C in both of these cores. Temperature estimates from the uppermost samples in all cores compare very well to the modern-day measured values. Affinities of some foraminiferal species for warmer or cooler surface temperatures are identified within the temperature range of the examined samples based on their abundance values. Especially notable among the warmer species are, Globorotalia menardii, Globigerinita glutinata, Globigerinoides ruber, and Globigerinoides sacculifer. Species indicative of cooler surface temperatures include Globorotalia inflata, Globigerina bulloides, Neogloboquadrina pachyderma, and Globigerina falconensis. A cluster analysis was carried out to assist in understanding the degree of variation which occurs in the foraminiferal assemblages, and how temperature differences influence the faunal compositions of the samples. It is demonstrated that fairly similar samples may have unexpectedly different estimated temperatures due to small differences in key species and, conversely, quite different assemblages can result in similar or identical temperature estimates which confirms that other parameters than just temperature affect faunal content.

Sea-surface paleotemperature reconstructions for the middle Miocene-Holocene of ODP Site 133-811 in the Coral Sea off northeast Australia

The derivation of a detailed sea-surface paleotemperature curve for the middle Miocene-Holocene (10-0 Ma) from ODP Site 811 on the Queensland Plateau, northeast Australia, has clarified the role of sea-surface temperature fluctuations as a control on the initiation and development of the extensive carbonate platforms of this region. This curve was derived from isotopic analyses of the planktonic foraminifer Globigerinoides ruber, and converted to temperature using the surface-water paleotemperature equation accounting for variations in global ice volume. The accuracy of these data were confirmed by derivation of paleotemperatures using the water column isotopic gradient (Delta delta18O), corrected for salinity and variations in seafloor water mass temperature. Results indicate that during this period surface-water temperatures were, on average, greater than the minimum required for tropical reef growth (20°C; Veron, 1986), with the exception of the late Miocene and earliest early Pliocene (10-4.9 Ma), when there were repeated intervals of temperatures between 18-20°C. Tropical reef growth on the Queensland Plateau was extensive from the early to early middle Miocene (~21-13 Ma), after which reef development began to decline. A lowstand near 11 Ma probably exposed shallower portions of the plateau; after re-immersion near 7 Ma, the areal extent of reef development was greatly reduced (~ 50%). Paleotemperature data from Site 811 indicate that decreased sea-surface temperatures were likely to have been instrumental in reducing the area of active reef growth on the Queensland Plateau. Reduced reefal growth rates continued until the late Pliocene or Quaternary, despite the increase of average sea-surface paleotemperatures to 22-23°C. Studies on modern corals show that when sea-surface temperatures are below ~24°C, as they were from the late Miocene to the Pleistocene off northeast Australia, corals are stressed and growth rates are greatly reduced. Consequently, when temperatures are in this range, corals have difficulty keeping pace with subsidence and changing environmental factors. In the late Pliocene, sedimentation rates increased due to increases in non-reefal carbonate production and falling sea levels. It was not until the mid-Quaternary (0.6-0.7 Ma) that sea-surface paleotemperatures increased above 24°C as a result of the formation of a western Coral Sea warm water pool. Because of age discrepancies, it is unclear exactly when an effective barrier developed on the central Great Barrier Reef; the formation of the warm water pool was likely to have either assisted the formation of this barrier and/or permitted increased coral growth rates. Fluctuations in sea-surface temperature can account for much of the observed spatial and temporal variations of reef growth and carbonate platform distribution off northeast Australia, and therefore we conclude that paleotemperature variations are a critical control on the development of carbonate platforms, and must be considered an important cause of ancient platform "drowning".

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

Sea surface temperature and export productivity reconstructions of sediments off the western Iberian margin

Present day hydrographic conditions along the western Iberian margin are characterized by seasonal upwelling with filaments that can penetrate more than 200 km into the open ocean and constitute areas of cold and highly productive waters. In order to investigate spatial and temporal gradients in temperature and productivity conditions during the last 150 ky, high-resolution proxy records were generated in 3 cores (SU92-03, MD95-2040, MD95-2042), located along the Iberian coast between 43°12'N and 37°48'N and forming a N-S profile. In all cores, planktonic foraminifera census counts are used to reconstruct summer sea surface temperature (SSTsu) and export productivity (Pexpsu) using the modern analog technique SIMMAX 28. SSTsu and Pexpsu values similar to the present are observed throughout the Holocene and MIS 5e periods for each site, respectively, indicating fairly stable conditions equivalent to the modern ones. On glacial/interglacial timescales, SSTsu increases by 2-3 °C from the northern to southernmost site. Pexpsu, on the other hand, shows a decrease of 30-40 gC/m**2/yr from North to South at present time and during interglacial periods, and no significant variation (90-100 gC/m**2/yr) during glacial periods. The northernmost core SU92-03 reveals the coldest conditions with records more similar to MD95-2040 than to MD95-2042, the later of which is, as at present, more affected by subtropical waters. Core SU92-03 shows higher interglacial productivity similar to open ocean mid- to high latitude sites, while the other two cores monitor higher glacial productivity conform with other upwelling sites off NW Africa. A boundary between differences in glacial/interglacial productivity appears to be present in our study between 43°12'N and 40°35'N. Especially north of 40°N, coldest SSTsu and lowest Pexpsu are found during Heinrich events (H)1-H8 and H10-H11. In contrast, lowest Pexpsu do not coincide with these events at site MD95-2042, but appear to be related to the presence of relatively warm and nutrient-poor subtropical Eastern North Atlantic Central Water advected with the Azores Current.

Stable isotope record, calcium carbonate conentrations, and sea surface temperture reconstructions of sediment cores from the North Atlantic

We analyze five high-resolution time series spanning the last 1.65 m.y.: benthic foraminiferal delta18O and delta13O, percent CaCO3, and estimated sea surface temperature (SST) at North Atlantic Deep Sea Drilling Project site 607 and percent CaCO3 at site 609. Each record is a multicore composite verified for continuity by splicing among multiple holes. These climatic indices portray changes in northern hemisphere ice sheet size and in North Atlantic surface and deep circulation. By tuning obliquity and precession components in the delta18O record to orbital variations, we have devised a time scale (TP607) for the entire Pleistocene that agrees in age with all K/Ar-dated magnetic reversals to within 1.5%. The Brunhes time scale is taken from Imbrie et al. [1984], except for differences near the stage 17/16 transition (0.70 to 0.64 Ma). All indicators show a similar evolution from the Matuyama to the Brunhes chrons: orbital eccentricity and precession responses increased in amplitude; those at orbital obliquity decreased. The change in dominance from obliquity to eccentricity occurred over several hundred thousand years, with fastest changes around 0.7 to 0.6 Ma. The coherent, in-phase responses of delta18O, delta13O, CaCO3 and SST at these rhythms indicate that northern hemisphere ice volume changes have controlled most of the North Atlantic surface-ocean and deep-ocean responses for the last 1.6 m.y. The delta13O, percent CaCO3, and SST records at site 607 also show prominent changes at low frequencies, including a prominent long-wavelength oscillation toward glacial conditions that is centered between 0.9 and 0.6 Ma. These changes appear to be associated neither with orbital forcing nor with changes in ice volume.

Age models, color measurements, and calculated sea surface temperatures of sediment cores from the Bering Sea and the Sea of Okhotsk

Past changes in North Pacific sea surface temperatures and sea-ice conditions are proposed to play a crucial role in deglacial climate development and ocean circulation but are less well known than from the North Atlantic. Here, we present new alkenone-based sea surface temperature records from the subarctic northwest Pacific and its marginal seas (Bering Sea and Sea of Okhotsk) for the time interval of the last 15 kyr, indicating millennial-scale sea surface temperature fluctuations similar to short-term deglacial climate oscillations known from Greenland ice-core records. Past changes in sea-ice distribution are derived from relative percentage of specific diatom groups and qualitative assessment of the IP25 biomarker related to sea-ice diatoms. The deglacial variability in sea-ice extent matches the sea surface temperature fluctuations. These fluctuations suggest a linkage to deglacial variations in Atlantic meridional overturning circulation and a close atmospheric coupling between the North Pacific and North Atlantic. During the Holocene the subarctic North Pacific is marked by complex sea surface temperature trends, which do not support the hypothesis of a Holocene seesaw in temperature development between the North Atlantic and the North Pacific.

Age determination and sea-surface reconstruction of sediment cores from the Norwegian Sea

Polycystine radiolarians are used to reconstruct summer sea surface temperatures (SSSTs) for the Late Pleistocene-Holocene (600-13,400 14C years BP) in the Norwegian Sea. At 13,200 14C years BP, the SSST was close to the average Holocene SSST (~12°C). It then gradually dropped to 7.1°C in the Younger Dryas. Near the Younger Dryas-Holocene transition (~10,000 14C years BP), the SSST increased 5°C in about 530 years. Four abrupt cooling events, with temperature drops of up to 2.1°C, are recognized during the Holocene: at 9340, 7100 ("8200 calendar years event"), 6400 and 1650 14C years BP. Radiolarian SSSTs and the isotopic signal from the GISP2 ice core are strongly coupled, stressing the importance of the Norwegian Sea as a mediator of heat/precipitation exchange between the North Atlantic, the atmosphere, and the Greenland ice sheet. Radiolarian and diatom-derived SSSTs display similarities, with the former not showing the recently reported Holocene cooling trend.

Model-based changes in global annual mean surface temperature change (Delta T_g) and radiative forcing due to land ice albedo changes (Delta R_[LI]) over the last 5 Myr, supplementary material

It is still an open question how equilibrium warming in response to increasing radiative forcing - the specific equilibrium climate sensitivity S - depends on background climate. We here present palaeodata-based evidence on the state dependency of S, by using CO2 proxy data together with a 3-D ice-sheet-model-based reconstruction of land ice albedo over the last 5 million years (Myr). We find that the land ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcing depends on the CO2 data set used. This non-linearity has not, so far, been accounted for in similar approaches due to previously more simplistic approximations, in which land ice albedo radiative forcing was a linear function of sea level change. The latitudinal dependency of ice-sheet area changes is important for the non-linearity between land ice albedo and sea level. In our set-up, in which the radiative forcing of CO2 and of the land ice albedo (LI) is combined, we find a state dependence in the calculated specific equilibrium climate sensitivity, S[CO2,LI], for most of the Pleistocene (last 2.1 Myr). During Pleistocene intermediate glaciated climates and interglacial periods, S[CO2,LI] is on average ~ 45 % larger than during Pleistocene full glacial conditions. In the Pliocene part of our analysis (2.6-5 Myr BP) the CO2 data uncertainties prevent a well-supported calculation for S[CO2,LI], but our analysis suggests that during times without a large land ice area in the Northern Hemisphere (e.g. before 2.82 Myr BP), the specific equilibrium climate sensitivity, S[CO2,LI], was smaller than during interglacials of the Pleistocene. We thus find support for a previously proposed state change in the climate system with the widespread appearance of northern hemispheric ice sheets. This study points for the first time to a so far overlooked non-linearity in the land ice albedo radiative forcing, which is important for similar palaeodata-based approaches to calculate climate sensitivity. However, the implications of this study for a suggested warming under CO2 doubling are not yet entirely clear since the details of necessary corrections for other slow feedbacks are not fully known and the uncertainties that exist in the ice-sheet simulations and global temperature reconstructions are large.

Plio-Pleistocene Mg/Ca and TEX86 sea surface temperature records of ODP Sites 184-1143 and 165-999,

The western warm pools of the Atlantic and Pacific oceans are a critical source of heat and moisture for the tropical climate system. Over the past five million years, global mean temperatures have cooled by 3-4 °C. Yet, current reconstructions of sea surface temperatures indicate that temperature in the warm pools has remained stable during this time. This stability has been used to suggest that tropical sea-surface temperatures are controlled by some sort of thermostat-like regulation. Here we reconstruct sea surface temperatures in the South China Sea, Caribbean Sea and western equatorial Pacific Ocean for the past five million years, using a combination of the Mg/Ca, TEXH86-and Uk'37 surface temperature proxies. Our data indicate that during the period of Pliocene warmth from about 5 to 2.6 million years ago, the western Pacific and western Atlantic warm pools were about 2 °C warmer than today. We suggest that the apparent lack of warming seen in the previous reconstructions was an artefact of low seawater Mg/Ca ratios in the Pliocene oceans. Taking this bias into account, our data indicate that tropical sea surface temperatures did change in conjunction with global mean temperatures. We therefore conclude that the temperature of the warm pools of the equatorial oceans during the Pliocene was not limited by a thermostat-like mechanism.

Uk'37, derived sea surface temperatures and coccolith abundances in North Atlantic sediment cores

The UK37' index has proven to be a robust proxy to estimate past sea surface temperatures (SSTs) over a range of time scales, but like any other proxy, it has uncertainties. For instance, in reconstructions of the Last Glacial Maximum (LGM) in the northern North Atlantic, UK37' indicates higher temperatures than those derived from foraminiferal proxies. Here we evaluate whether such warm glacial estimates are caused by the advection of reworked alkenones in ice-rafted debris (IRD) to deep-sea sediments. We have quantified both coccolith assemblages and alkenones in sediments from glaciogenic debris flows in the continental margins of the northern North Atlantic, and from a deep-sea core from the Reykjanes Ridge. Certain debris flow deposits in the North Atlantic were generated by the presence of massive ice-sheets in the past, and their associated ice streams. Such deposits are composed of the same materials that were present in the IRD at the time they were generated. We conclude that ice rafting from some locations was a transport pathway to the deep sea floor of reworked alkenones and pre-Quaternary coccolith species during glacial stages, but that not all of the IRD contained alkenones, even when reworked coccoliths were present. We speculate that the ratio of reworked coccoliths to alkenone concentration might be useful to infer whether significant reworked alkenone inputs from IRD did occur at a particular site in the glacial North Atlantic. We also observe that alkenones in some of the debris flows contain a colder signal than estimated for LGM sediments in the northern North Atlantic. This is also clear in the deep-sea core studied where the warmest intervals do not correspond to the intervals with large inputs of reworked coccoliths or IRD. We conclude that any possible bias to UK37' estimates associated with reworked alkenones is not necessarily towards higher values, and that the high SST anomalies for the LGM are unlikely to be the result of a bias caused by IRD inputs.

Paleoclimate reconstructions from sediment core MD01-2378

On the basis of the radiocarbon (14C) plateau-tuning method a new age model for Timor Sea Core MD01-2378 was established. It revealed a precise centennial-scale phasing of climate events in the ocean, cryo-, and atmosphere during the last deglacial and provides important new insights into causal linkages controlling events of global climate change. At Site MD01-2378 reservoir ages of surface waters dropped from 1600 yr prior to 20 cal ka to 250-500 yr after 18.8 cal ka. This evidence was crucial for generating a high-resolution age model for deglacial events in the Indo-Pacific Warm Pool. Sea-surface temperatures (SST) started to change near 18.8 cal ka, that is ~500 yr after the start of, presumably northern hemispheric, deglacial melt and sea level rise as shown by the benthic foraminiferal oxygen isotope ratio (d18O). However, the SST rise occurred 500-1000 yr prior to the onset of deglacial Antarctic warming and the first major rise in atmospheric carbon dioxide at about 18 ka. The increase in SST may partly reflect reduced seasonal upwelling of cold subsurface waters along the eastern margin of the Indian Ocean, which is reflected by a doubling of the thermal gradient between the sea surface and the thermocline, a halving of chlorin productivity from 19 to 18.5 cal ka, and in particular, by the strong decrease in surface water reservoir ages. Two significant increases in deglacial Timor Sea surface salinities from 19-18.5 and 15.5-14.5 cal ka, may partly reflect the deglacial increase in the distance of local river mouths, partly an inter-hemispheric millennial-scale see-saw in tropical monsoon intensity, possibly linked to a deglacial increase in the dominance of Pacific El Niño regimes over Heinrich stadial 1.