Magnesium/Calcium ratios and sea surface temperature estimation for sediments of the Timor Sea

We use a multiproxy approach to monitor changes in the vertical profile of the Indonesian Throughflow as well as monsoonal wind and precipitation patterns in the Timor Sea on glacial-interglacial, precessional, and suborbital timescales. We focus on an interval of extreme climate change and sea level variation: marine isotope (MIS) 6 to MIS 5e. Paleoproductivity fluctuations in the Timor Sea follow a precessional beat related to the intensity of the Australian (NW) monsoon. Paired Mg/Ca and d18O measurements of surface- and thermocline-dwelling planktonic foraminifers (G. ruber and P. obliquiloculata) indicate an increase of >4°C in both surface and thermocline water temperatures during Termination II. Tropical sea surface temperature changed synchronously with ice volume (benthic d18O) during deglaciation, implying a direct coupling of high- and low-latitude climate via atmospheric and/or upper ocean circulation. Substantial cooling and freshening of thermocline waters occurred toward the end of Termination II and during MIS 5e, indicating a change in the vertical profile of the Indonesian Throughflow from surface- to thermocline-dominated flow.

Stable isotope ratios of foraminifera and sea surface temperature estimation of sediment core W2804A-14

Carbon isotopically based estimates of CO2 levels have been generated from a record of the photosynthetic fractionation of 13C (epsilon p) in a central equatorial Pacific sediment core that spans the last ~255 ka. Contents of 13C in phytoplanktonic biomass were determined by analysis of C37 alkadienones. These compounds are exclusive products of Prymnesiophyte algae which at present grow most abundantly at depths of 70-90 m in the central equatorial Pacific. A record of the isotopic compostion of dissolved CO2 was constructed from isotopic analyses of the planktonic foraminifera Neogloboquadrina dutertrei, which calcifies at 70-90 m in the same region. Values of epsilon p, derived by comparison of the organic and inorganic delta values, were transformed to yield concentrations of dissolved CO2 (c e) based on a new, site-specific calibration of the relationship between epsilon p and c e. The calibration was based on reassessment of existing epsilon p versus c e data, which support a physiologically based model in which epsilon p is inversely related to c e. Values of PCO2, the partial pressure of CO2 that would be in equilibrium with the estimated concentrations of dissolved CO2, were calculated using Henry's law and the temperature determined from the alkenone-unsaturation index UK 37. Uncertainties in these values arise mainly from uncertainties about the appropriateness (particularly over time) of the site-specific relationship between epsilon p and 1/c e. These are discussed in detail and it is concluded that the observed record of epsilon p most probably reflects significant variations in Delta pCO2, the ocean-atmosphere disequilibrium, which appears to have ranged from ~110 µatm during glacial intervals (ocean > atmosphere) to ~60 µatm during interglacials. Fluxes of CO2 to the atmosphere would thus have been significantly larger during glacial intervals. If this were characteristic of large areas of the equatorial Pacific, then greater glacial sinks for the equatorially evaded CO2 must have existed elsewhere. Statistical analysis of air-sea pCO2 differences and other parameters revealed significant (p < 0.01) inverse correlations of Delta pCO2 with sea surface temperature and with the mass accumulation rate of opal. The former suggests response to the strength of upwelling, the latter may indicate either drawdown of CO2 by siliceous phytoplankton or variation of [CO2]/[Si(OH)4] ratios in upwelling waters.

Sea surface temperature estimation for the last 70 ka from IODP Site 306-U1313

During the six Heinrich Events of the last 70 ka episodic calving from the circum-Atlantic ice sheets released large numbers of icebergs into the North Atlantic. These icebergs and associated melt-water flux are hypothesized to have led to a shutdown of Atlantic Meridional Overturning Circulation (AMOC) and severe cooling in large parts of the Northern Hemisphere. However, due to the limited availability of high-resolution records the magnitude sea surface temperature (SST) changes related to the impact of Heinrich Events on the mid-latitude North Atlantic is poorly constrained. Here we present a record of UK37'-based SSTs derived from sediments of Integrated Ocean Drilling Project (IODP) Site U1313, located at the southern end of the ice-rafted debris (IRD)-belt in the mid-latitude North Atlantic (41°N). We demonstrate that all six Heinrich Events are associated with a rapid warming of surface waters by 2 to 4°C in a few thousand years. The presence of IRD leaves no doubt about the simultaneous timing and correlation between rapid surface water warming and Heinrich Events. We argue that this warming in the mid-latitude North Atlantic is related to a northward expansion of the subtropical gyre during Heinrich Events. As a wide-range of studies demonstrated that in the central IRD-belt Heinrich Events are associated with low SSTs, these results thus identify an anti-phased (seesaw) pattern in SSTs during Heinrich Events between the mid-latitude (warm) and northern North Atlantic (cold). This highlights the complex response of surface water characteristics in the North Atlantic to Heinrich Events that is poorly reproduced by fresh water hosing experiments and challenges the widely accepted view that within the IRD-belt of the North Atlantic Heinrich Events coincide with periods of low SSTs.

Sea surface temperature estimation for MIS 16 of IODP Site 306-U1313

Hudson Strait (HS) Heinrich Events, ice-rafting events in the North Atlantic originating from the Laurentide ice sheet (LIS), are among the most dramatic examples of millennial-scale climate variability and have a large influence on global climate. However, it is debated as to whether the occurrence of HS Heinrich Events in the (eastern) North Atlantic in the geological record depends on greater ice discharge, or simply from the longer survival of icebergs in cold waters. Using sediments from Integrated Ocean Drilling Program (IODP) Site U1313 in the North Atlantic spanning the period between 960 and 320 ka, we show that sea surface temperatures (SSTs) did not control the first occurrence of HS Heinrich(-like) Events in the sedimentary record. Using mineralogy and organic geochemistry to determine the characteristics of ice-rafting debris (IRD), we detect the first HS Heinrich(-like) Event in our record around 643 ka (Marine Isotope Stage (MIS) 16), which is similar as previously reported for Site U1308. However, the accompanying high-resolution alkenone-based SST record demonstrates that the first HS Heinrich(-like) Event did not coincide with low SSTs. Thus, the HS Heinrich(-like) Events do indicate enhanced ice discharge from the LIS at the end of the Mid-Pleistocene Transition, not simply the survivability of icebergs due to cold conditions in the North Atlantic.

Sea surface temperature estimation of the South Tasman Rise

The Subtropical Front (STF) marking the northern boundary of the Southern Ocean has a steep gradient in sea surface temperature (SST) of approximately 4°C over 0.5° of latitude. Presently, in the region south of Tasmania, the STF lies nominally at 47°S in the summer and 45°S in the winter. We present here SST reconstructions in a latitudinal transect of cores across the South Tasman Rise, southeast of Australia, during the late Quaternary. SST reconstructions are based on two paleotemperature proxies, alkenones and faunal assemblages, which are used to assess past changes in SST in spring and summer. The north-south alignment in core locations allows reconstruction of movement of the STF over the last 100 ka. Surface water temperatures during the last glaciation in this region were ~4°C colder than today. Additional temperature changes greater in magnitude than 4°C seen in individual cores can be attributed to changes in the water mass overlying the core site caused by the movement of the front across that location. During the penultimate interglacial, SST was ~2°C warmer and the STF was largely positioned south of 47°S. Movement of the STF to the north occurred during cool climate periods such as the last marine isotope stages 3 and 4. In the last glaciation, the front was at its farthest north position, becoming pinned against the Tasmanian landmass. It moved south by 4° latitude to 47°S in summer during the deglaciation but remained north of 45°S in spring throughout the early deglaciation. After 11 ka B.P. inferred invigoration of the East Australia Current appears to have pushed the STF seasonally south of the East Tasman Plateau, until after 6 ka B.P. when it achieved its present configuration.

TEX86 values and temperature estimation for sediment trap samples

The newly introduced temperature proxy, the tetraether index of archaeal lipids with 86 carbon atoms (TEX86), is based on the number of cyclopentane moieties in the glycerol dialkyl glycerol tetraether (GDGT) lipids of marine Crenarchaeota. The composition of sedimentary GDGTs used for TEX86 paleothermometry is thought to reflect sea surface temperature (SST). However, marine Crenarchaeota occur ubiquitously in the world oceans over the entire depth range and not just in surface waters. We analyzed the GDGT distribution in settling particulate organic matter collected in sediment traps from the northeastern Pacific Ocean and the Arabian Sea to investigate the seasonal and spatial distribution of the fluxes of crenarchaeotal GDGTs and the origin of the TEX86 signal transported to the sediment. In both settings the TEX86 measured at all trap deployment depths reflects SST. In the Arabian Sea, analysis of an annual time series showed that the SST estimate based on TEX86 in the shallowest trap at 500 m followed the in situ SST with a 1 to 3 week time delay, likely caused by the relatively low settling speed of sinking particles. This revealed that the GDGT signal that reaches deeper water is derived from the upper water column rather than in situ production of GDGTs. The GDGT temperature signal in deeper traps at 1500 m and 3000 m did not show a seasonal cyclicity observed in the 500 m trap but rather reflected the annual mean SST. This is probably due to a homogenization of the TEX86 SST signal carried by particles as they ultimately reach the interior of the ocean. Our data confirm the use of TEX86 as a temperature proxy of surface ocean waters.