Comparison of TEX86 and UK'37 temperature proxies in sinking particles in the Cariaco Basin

The Cariaco Basin, a silled, permanently anoxic basin on the continental shelf of Venezuela with a dynamic chemocline (-240-350 m), has been subject of > 20 years of oceanographic observation and sediment trap studies. We evaluated UK'37 and the TEX86 temperature proxies using sinking particles collected in shallow sediment trap samples at 275 m (Trap A) and 455 m (Trap B) (within and below the chemocline). The organic geochemical temperature proxies, UK'37. (based on coccolithophorid alkenone lipids) and TEX86 (based on archaeal glycerol dialkyl glycerol tetraether (GDGT) lipids), use observed relationships between the ratio of specific lipids and measured sea surface temperature to hindcast past sea surface temperatures. In this study, both UK'37 and TEX86 temperature proxies record seasonal temperature variations, including the cooling associated with upwelling events. UK'37-based temperatures are colder than measured sea surface temperatures, and better correlated temperature at the chlorophyll maximum. In sediment trap material collected below the chemocline (Trap B), UK'37 values are higher than those in Trap A. Warmer subchemocline UK'37 based temperatures may be related to autooxidation of sinking particles, either by small amounts of available oxygen or by alternate electron acceptors concentrated in the biologically dynamic chemocline (e.g. intermediate sulfur compounds). The absolute flux weighted TEX86 temperature values measured in sinking particles from Trap A match the measured SST well. The differences in the TEX86 values between Traps A and B are small and reflect less impact of degradation. Overall, the TEX86 temperatures in sinking particles in the Cariaco Basin reflect annual SST.

Temperature proxies of sediment cores from the Sea of Marmara

Reconstructing ocean temperature values is a major target in paleoceanography and climate research. However, most temperature proxies are organism-based and thus suffer from an "ecological bias". Multiproxy approaches can potentially overcome this bias but typically require more investment in time and resources, while being susceptible to errors induced by sample preparation steps necessary before analysis. Three lipid-based temperature proxies are widely used: UK'37 (based on long chain alkenones from phytoplanktonic haptophytes), TEX86 [based on glycerol dialkyl glycerol tetraethers (GDGTs) from pelagic archaea] and LDI (based on long chain diols from phytoplanktonic eustigmatophytes). So far, separate analytical methods, including gas chromatography (GC) and liquid chromatography (LC), have been used to determine these proxies. Here we present a sensitive method for determining all three in a single normal phase high performance LC-atmospheric pressure chemical ionization mass spectrometry (NP-HPLC-APCI-MS) analysis. Each of the long chain alkenones and long chain diols was separated and unambiguously identified from the accurate masses and characteristic fragmentation during multiple stage MS analysis (MS2). Comparison of conventional GC and HPLC-MS methods showed similar results for UK'37 and LDI, respectively, using diverse environmental samples and an Emiliania huxleyi culture. Including the three sea surface temperature (SST) proxies; the NP-HPLC-APCI-MS method in fact allows simultaneous determination of nine paleoenvironmental proxies. The extent to which the ecology of the source organisms (ecological bias) influences lipid composition and thereby the reconstructed temperature values was demonstrated by applying the new method to a sediment core from the Sea of Marmara, covering the last 21 kyr BP. Reconstructed SST values differed considerably between the proxies for the Last Glacial Maximum (LGM) and the period of Sapropel S1 formation at ca. 10 kyr BP, whereas the trends during the late Holocene were similar. Changes in the composition of alkenone-producing species at the transition from the LGM to the Bølling/Allerød (B/A) were inferred from unreasonably high UK'37-derived SST values (ca. 20 °C) during the LGM. We ascribe discrepancies between the reconstructed temperature records during S1 deposition to habitat change, e.g. a different depth due to changes in nutrient availability.

Sea surface temperature reconstruction for ODP Hole 178-1098B

The disintegration of ice shelves, reduced sea-ice and glacier extent, and shifting ecological zones observed around Antarctica (Cook et al., 2005, doi:10.1126/science.1104235; Stammerjohn et al., 2008, doi:10.1016/j.dsr2.2008.04.026) highlight the impact of recent atmospheric (Steig et al., 2009, doi:10.1038/nature07669) and oceanic warming (Gille, 2002, doi:10.1126/science.1065863) on the cryosphere. Observations (Cook et al., 2005, doi:10.1126/science.1104235; Stammerjohn et al., 2008, doi:10.1016/j.dsr2.2008.04.026) and models (Pollard and DeConto, 2009, doi:10.1038/nature07809) suggest that oceanic and atmospheric temperature variations at Antarctica's margins affect global cryosphere stability, ocean circulation, sea levels and carbon cycling. In particular, recent climate changes on the Antarctic Peninsula have been dramatic, yet the Holocene climate variability of this region is largely unknown, limiting our ability to evaluate ongoing changes within the context of historical variability and underlying forcing mechanisms. Here we show that surface ocean temperatures at the continental margin of the western Antarctic Peninsula cooled by 3-4 °C over the past 12,000?years, tracking the Holocene decline of local (65° S) spring insolation. Our results, based on TEX86 sea surface temperature (SST) proxy evidence from a marine sediment core, indicate the importance of regional summer duration as a driver of Antarctic seasonal sea-ice fluctuations (Huybers and Denton, 2008, doi:10.1038/ngeo311). On millennial timescales, abrupt SST fluctuations of 2-4 °C coincide with globally recognized climate variability (Mayewski et al., 2004, doi:10.1016/j.yqres.2004.07.001). Similarities between our SSTs, Southern Hemisphere westerly wind reconstructions (Moreno et al., 2010, doi:10.1130/G30962.1) and El Niño/Southern Oscillation variability (Conroy et al., 2008, doi:10.1016/j.quascirev.2008.02.015) indicate that present climate teleconnections between the tropical Pacific Ocean and the western Antarctic Peninsula (Yuan et al., 2004, doi:10.1017/S0954102004002238) strengthened late in the Holocene epoch. We conclude that during the Holocene, Southern Ocean temperatures at the western Antarctic Peninsula margin were tied to changes in the position of the westerlies, which have a critical role in global carbon cycling (Moreno et al., 2010, doi:10.1130/G30962.1; Anderson et al., 2009, doi:10.1126/science.1167441).

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.

Sea surface temperture reconstruction for ODP Site 189-1172

Relative to the present day, meridional temperature gradients in the Early Eocene age (~56-53 Myr ago) were unusually low, with slightly warmer equatorial regions (Pearson et al., 2007, doi:10.1130/G23175A.1 ) but with much warmer subtropical Arctic (Sluijs et al., 2008, doi:10.1029/2007PA001495) and mid-latitude (Sluijs et al., 2007, doi:10.1038/nature06400) climates. By the end of the Eocene epoch (~34 Myr ago), the first major Antarctic ice sheets had appeared (Zachos et al., 1992, doi:10.1130/0091-7613(1992)020<0569:EOISEO>2.3.CO;2; Barker et al., 2007, doi:10.1016/j.dsr2.2007.07.027), suggesting that major cooling had taken place. Yet the global transition into this icehouse climate remains poorly constrained, as only a few temperature records are available portraying the Cenozoic climatic evolution of the high southern latitudes. Here we present a uniquely continuous and chronostratigraphically well-calibrated TEX86 record of sea surface temperature (SST) from an ocean sediment core in the East Tasman Plateau (palaeolatitude ~65° S). We show that southwest Pacific SSTs rose above present-day tropical values (to ~34° C) during the Early Eocene age (~53 Myr ago) and had gradually decreased to about 21° C by the early Late Eocene age (~36 Myr ago). Our results imply that there was almost no latitudinal SST gradient between subequatorial and subpolar regions during the Early Eocene age (55-50 Myr ago). Thereafter, the latitudinal gradient markedly increased. In theory, if Eocene cooling was largely driven by a decrease in atmospheric greenhouse gas concentration Zachos et al. (2008, doi:10.1038/nature06588), additional processes are required to explain the relative stability of tropical SSTs given that there was more significant cooling at higher latitudes.

Geochemical and palynological measurements from the ODP Leg 189 at Site 1172 on the East Tasman Plateau

A brief (~150 kyr) period of widespread global average surface warming marks the transition between the Paleocene and Eocene epochs, ~56 million years ago. This so-called "Paleocene-Eocene thermal maximum" (PETM) is associated with the massive injection of 13C-depleted carbon, reflected in a negative carbon isotope excursion (CIE). Biotic responses include a global abundance peak (acme) of the subtropical dinoflagellate Apectodinium. Here we identify the PETM in a marine sedimentary sequence deposited on the East Tasman Plateau at Ocean Drilling Program (ODP) Site 1172 and show, based on the organic paleothermometer TEX86, that southwest Pacific sea surface temperatures increased from ~26 °C to ~33°C during the PETM. Such temperatures before, during and after the PETM are >10 °C warmer than predicted by paleoclimate model simulations for this latitude. In part, this discrepancy may be explained by potential seasonal biases in the TEX86 proxy in polar oceans. Additionally, the data suggest that not only Arctic, but also Antarctic temperatures may be underestimated in simulations of ancient greenhouse climates by current generation fully coupled climate models. An early influx of abundant Apectodinium confirms that environmental change preceded the CIE on a global scale. Organic dinoflagellate cyst assemblages suggest a local decrease in the amount of river run off reaching the core site during the PETM, possibly in concert with eustatic rise. Moreover, the assemblages suggest changes in seasonality of the regional hydrological system and storm activity. Finally, significant variation in dinoflagellate cyst assemblages during the PETM indicates that southwest Pacific climates varied significantly over time scales of 103 - 104 years during this event, a finding comparable to similar studies of PETM successions from the New Jersey Shelf.

Age model, lipids, and iron/calcium-ratios of sediment core GeoB9307-3

Hide Intense debate persists about the climatic mechanisms governing hydrologic changes in tropical and subtropical southeast Africa since the Last Glacial Maximum, about 20,000 years ago. In particular, the relative importance of atmospheric and oceanic processes is not firmly established. Southward shifts of the intertropical convergence zone (ITCZ) driven by high-latitude climate changes have been suggested as a primary forcing, whereas other studies infer a predominant influence of Indian Ocean sea surface temperatures on regional rainfall changes. To address this question, a continuous record representing an integrated signal of regional climate variability is required, but has until now been missing. Here we show that remote atmospheric forcing by cold events in the northern high latitudes appears to have been the main driver of hydro-climatology in southeast Africa during rapid climate changes over the past 17,000 years. Our results are based on a reconstruction of precipitation and river discharge changes, as recorded in a marine sediment core off the mouth of the Zambezi River, near the southern boundary of the modern seasonal ITCZ migration. Indian Ocean sea surface temperatures did not exert a primary control over southeast African hydrologic variability. Instead, phases of high precipitation and terrestrial discharge occurred when the ITCZ was forced southwards during Northern Hemisphere cold events, such as Heinrich stadial 1 (around 16,000 years ago) and the Younger Dryas (around 12,000 years ago), or when local summer insolation was high in the late Holocene, i.e., during the last 4,000 years.

Palynology and sea surface temperature reconstruction for the mid-Cretaceous Oceanic Anoxic Event 2 at ODP Hole 210-1276A

The mid-Cretaceous is thought to be a greenhouse world with significantly higher atmospheric pCO2 and sea-surface temperatures as well as a much flatter latitudinal thermal gradient compared to the present. This time interval was punctuated by the Cenomanian/Turonian Oceanic Anoxic Event (OAE-2, ~ 93.5 Myr ago), an episode of global, massive organic carbon burial that likely resulted in a large and abrupt pCO2 decline. However, the climatic consequences of this pCO2 drop are yet poorly constrained. We determined the first, high-resolution sea-surface temperature (SST) record across OAE-2 from a deep-marine sedimentary sequence at Ocean Drilling Program (ODP) Site 1276 in the mid-latitudinal Newfoundland Basin, NW Atlantic. By employing the organic palaeothermometer TEX86, we found that SSTs across the OAE-2 interval were extremely high, but were punctuated by a remarkably large cooling (5-11 °C), which is synchronous with the 2.5-5.5 °C cooling in SST records from equatorial Atlantic sites, and the "Plenus Cold Event". Because this global cooling event is concurrent with increased organic carbon burial, it likely acted in response to the associated pCO2 drop. Our findings imply a substantial increase in the latitudinal SST gradient in the proto-North Atlantic during this period of global cooling and reduced atmospheric pCO2, suggesting a strong coupling between pCO2 and latitudinal thermal gradients under greenhouse climate conditions.

Revised calibration of the MBT-CBT paleotemperature proxy

The MBT-CBT proxy for the reconstruction of paleotemperatures and past soil pH is based on the distribution of branched glycerol dialkyl glycerol tetraether (brGDGT) membrane lipids. The Methylation of Branched Tetraether (MBT) and the Cyclisation of Branched Tetraether (CBT) indices were developed to quantify these distributions, and significant empirical relations between these indices and annual mean air temperature (MAT) and/or soil pH were found in a large data set of soils. In this study, we extended this soil dataset to 278 globally distributed surface soils. Of these soils, 26% contains all nine brGDGTs, while in 63% of the soils the seven most common brGDGTs were detected, and the latter were selected for calibration purposes. This resulted in new transfer functions for the reconstruction of pH based on the CBT index: pH = 7.90-1.97 × CBT (r**2 = 0.70; RMSE = 0.8; n = 176), as well as for MAT based on the CBT index and methylation index based on the seven most abundant GDGTs (defined as MBT'): MAT = 0.81-5.67 × CBT + 31.0 × MBT' (r**2 = 0.59; RMSE = 5.0 °C; n = 176). The new transfer function for MAT has a substantially lower correlation coefficient than the original equation (r**2 = 0.77). To investigate possible improvement of the correlation, we used our extended global surface soil dataset to statistically derive the indices that best describe the relations of brGDGT composition with MAT and soil pH. These new indices, however, resulted in only a relatively minor increase in correlation coefficients, while they cannot be explained straightforwardly by physiological mechanisms. The large scatter in the calibration cannot be fully explained by local factors or by seasonality, but MAT for soils from arid regions are generally substantially (up to 20 °C) underestimated, suggesting that absolute brGDGT-based temperature records for these areas should be interpreted with caution. The applicability of the new MBT'-CBT calibration function was tested using previously published MBT-CBT-derived paleotemperature records covering the last deglaciation in Central Africa and East Asia, the Eocene-Oligocene boundary and the Paleocene-Eocene thermal maximum. The results show that trends remain similar in all records, but that absolute temperature estimates and the amplitude of temperature changes are lower for most records, and generally in better agreement with independent proxy data.

Abrupt vegetation change after the Late Quaternary megafaunal extinction in southeastern Australia

A substantial extinction of megafauna occurred in Australia between 50 and 45 kyr ago, a period that coincides with human colonization of Australia. Large shifts in vegetation also occurred around this time, but it is unclear whether the vegetation changes were driven by the human use of fire-and thus contributed to the extinction event-or were a consequence of the loss of megafaunal grazers. Here we reconstruct past vegetation changes in southeastern Australia using the stable carbon isotopic composition of higher plant wax n-alkanes and levels of biomass burning from the accumulation rates of the biomarker levoglucosan from a well-dated sediment core offshore from the Murray-Darling Basin. We find that from 58 to 44 kyr ago, the abundance of plants with the C-4 carbon fixation pathway was generally high-between 60 and 70%. By 43 kyr ago, the abundance of C-4 plants dropped to 30% and biomass burning increased. This transient shift lasted for about 3,000 years and came after the period of human arrival and directly followed megafauna extinction at 48.9-43.6 kyr ago. We conclude that the vegetation shift was not the cause of the megafaunal extinction in this region. Instead, our data are consistent with the hypothesis that vegetation change was the consequence of the extinction of large browsers and led to the build-up of fire-prone vegetation in the Australian landscape.

Paleoclimate record combinatio of speleothems and of sediment core GeoB3910-2

A substantial strengthening of the South American monsoon system (SAMS) during Heinrich Stadials (HS) points toward decreased cross-equatorial heat transport as the main driver of monsoonal hydroclimate variability at millennial time-scales. In order to better constrain the exact timing and internal structure of HS1 over tropical South America we assessed two precisely dated speleothem records from central-eastern and northeastern Brazil in combination with two marine records of terrestrial organic and inorganic matter input into the western equatorial Atlantic. During HS1 we recognize at least two events of widespread intensification of the SAMS across the entire region influenced by the South Atlantic Convergence Zone (SACZ) at 16.11-14.69 kyr BP and 18.1-16.66 kyr BP (labeled as HS1a and HS1c, respectively), separated by a dry excursion from 16.66-16.11 kyr BP (HS1b). In view of the spatial structure of precipitation anomalies, the widespread increase of monsoon precipitation over the SACZ domain was termed 'Mega-SACZ'.