The early Eocene global warming events: Implications from ODP Leg 208 Walvis Ridge depth transect

Five sections drilled in multiple holes over a depth transect of more than 2200 m at the Walvis Ridge (SE Atlantic) during Ocean Drilling Program (ODP) Leg 208 resulted in the first complete early Paleogene deep-sea record. Here we present high-resolution stratigraphic records spanning a ~4.3 million yearlong interval of the late Paleocene to early Eocene. This interval includes the Paleocene-Eocene thermal maximum (PETM) as well as the Eocene thermal maximum (ETM) 2 event. A detailed chronology was developed with nondestructive X-ray fluorescence (XRF) core scanning records and shipboard color data. These records were used to refine the shipboard-derived spliced composite depth for each site and with a record from ODP Site 1051 were then used to establish a continuous time series over this interval. Extensive spectral analysis reveals that the early Paleogene sedimentary cyclicity is dominated by precession modulated by the short (100 kyr) and long (405 kyr) eccentricity cycles. Counting of precession-related cycles at multiple sites results in revised estimates for the duration of magnetochrons C24r and C25n. Direct comparison between the amplitude modulation of the precession component derived from XRF data and recent models of Earth's orbital eccentricity suggests that the onset of the PETM and ETM2 are related to a 100-kyr eccentricity maximum. Both events are approximately a quarter of a period offset from a maximum in the 405-kyr eccentricity cycle, with the major difference that the PETM is lagging and ETM2 is leading a 405-kyr eccentricity maximum. Absolute age estimates for the PETM, ETM2, and the magnetochron boundaries that are consistent with recalibrated radiometric ages and recent models of Earth's orbital eccentricity cannot be precisely determined at present because of too large uncertainties in these methods. Nevertheless, we provide two possible tuning options, which demonstrate the potential for the development of a cyclostratigraphic framework based on the stable 405-kyr eccentricity cycle for the entire Paleogene.

Sedimentation rate, bulk density, accumulation rate of terrigenous material and grain sizes of ODP Leg 167 sites

The terrigenous mineral fraction of sediments recovered by drilling during Ocean Drilling Program Leg 167 at Sites 1018 and 1020 is used to evaluate changes in the source and transport of fine-grained terrigenous sediment and its relation to regional climates and the paleoceanographic evolution of the California Current system during the late Pleistocene. Preliminary time scales developed by correlation of oxygen isotope stratigraphies with the global SPECMAP record show average linear sedimentation rates in excess of 100 m/m.y., which provide an opportunity for high-resolution studies of terrigenous flux, grain size, and mineralogy. The mass flux of terrigenous minerals at Site 1018 varies from 5 to 30 g/(cm**2 x k.y.) and displays a general trend toward increased flux during glacials. The terrigenous record at Site 1020 shows a similar pattern of increased glacial input, but overall accumulation rates are significantly lower. Spectral analysis demonstrates that most of this variability is concentrated in frequency bands related to orbital cycles of eccentricity, tilt, and precession. Detailed grain-size analysis performed on the isolated terrigenous mineral fraction shows that sediments from Site 1018 are associated with higher energy transport and depositional regimes than those found at Site 1020. Grain-size data are remarkably uniform throughout the last 500 k.y., with no discernible difference observed between glacial and interglacial size distributions within each site. X-ray diffraction analysis of the <2-µm clay component suggests that the deposition of minerals found at Site 1020 is consistent with transport from a southern source during intervals of increased terrigenous input.

(Table 1) Sapropel analysis from ODP Sites 160-964, 160-966, 160-967 and 160-969

Seventeen eastern Mediterranean Pliocene sapropels from ODP Sites 964, 966, 967 and 969, some of which are coeval, have been analysed for their geochemistry. The sapropels are characterized by very high organic carbon contents (up to 30%) which are reported to be the result of both increased productivity and improved preservation. Although the organic matter in the sapropels is mainly of marine origin, the d13Corg values and C/N ratios appear "terrestrial". This is the result of anaerobic organic matter degradation which preferentially removed nitrogen- and 13C-rich organic components. A comparison with Ti/Al profiles, which mimic the precession index, and a calculation of organic carbon accumulation rates indicate that sedimentation rates were at most 30% lower or at most 50% higher during sapropel formation. Thus, sapropel formation lasted from between 2000 and 10,000 years at Site 964 to between 4500 and 12,000 years at Site 967. A synthesis of new data and a comparison with existing models indicates that productivity, which increased due to extra nutrients supplied as a result of winter mixing and as a result of enhanced input by the Nile, was the driving mechanism behind sapropel formation. The resulting sapropel formation was simultaneous at different depths, but lasted longer in the part of the basin closest to the Nile.

Organic carbon content, stable isotopes and planktonic foraminifera in sediments of tghe norther Arabian Sea oxygen minimum zone

The northern Arabian Sea is one of the few regions in the open ocean where thermocline water is severely depleted in oxygen. The intensity of this oxygen minimum zone (OMZ) has been reconstructed over the past 225,000 years using proxies for surface water productivity, water column denitrification, winter mixing, and the aragonite compensation depth (ACD). Changes in OMZ intensity occurred on orbital and suborbital timescales. Lowest O2 levels correlate with productivity maxima and shallow winter mixing. Precession-related productivity maxima lag early summer insolation maxima by ~6 kyr, which we attribute to a prolonged summer monsoon season related to higher insolation at the end of the summer. Periods with a weakened or even non-existent OMZ are characterized by low productivity conditions and deep winter mixing attributed to strong and cold winter monsoonal winds. The timing of deep winter mixing events corresponds with that of periods of climatic cooling in the North Atlantic region.

Petrologic data from Site U1308 over the last three glaciations

A composite North Atlantic record from DSDP Site 609 and IODP Site U1308 spans the past 300,000 years and shows that variability within the penultimate glaciation differed substantially from that of the surrounding two glaciations. Hematite stained grains exhibit similar repetitive down-core variations within the Marine Isotope Stage (MIS) 8 and 4-2 intervals, but little cyclic variability within the MIS 6 section. There is also no petrologic evidence, in terms of detrital carbonate-rich (Heinrich) layers, for surging of the Laurentide Ice Sheet through the Hudson Strait during MIS 6. Rather, very high background concentration of ice-rafted debris (IRD) indicates near continuous glacial meltwater input that likely increased thermohaline disruption sensitivity to relatively weak forcing events, such as expanded sea ice over deepwater formation sites. Altered (sub)tropical precipitation patterns and Antarctic warming during high orbital precession and low 65° N summer insolation appears related to high abundance of Icelandic glass shards and southward sea ice expansion. Differing European and North American ice sheet configurations, perhaps aided by larger variations in eccentricity leading to cooler summers, may have contributed to the relative stability of the Laurentide Ice Sheet in the Hudson Strait region during MIS 6.

Stable isotope record and age determination of ODP Site 121-758 in the northeastern Indian Ocean

Numerous studies have shown that delta18O records from benthic and planktonic foraminifera, primarily a proxy of global ice volume variations, reflect Milankovitch periodicities. To study climatic response to orbital forcing at Ocean Drilling Program site 758, we have generated continuous delta18O and delta13C records from a single benthic foraminiferal species Cibicides wuellerstorfi for the last 3.6 m.y. and extended the planktonic foraminiferal isotope records of Farrell and Janecek (1991, doi:10.2973/odp.proc.sr.121.124.1991) (0-2.5 Ma, based on Globigerinoides sacculifer) to 3.6 Ma (Chen, 1994). We then constructed an age model by matching, correlating and tuning the benthic delta18O record to a model simulation of ice volume (Imbrie and Imbrie, 1980, doi:10.1126/science.207.4434.943). The filtered 41- and 23-kyr signals based on the resultant astronomically tuned age model are highly correlated to obliquity (r=0.83) and precession (r=0.75), respectively. Although derived with methodology different from Shackleton et al. (1990) and Hilgen (1991, doi:10.1016/0012-821X(91)90206-W, 1991, doi:10.1016/0012-821X(91)90082-S), our results generally agree with their published astronomical timescales for the time interval from 0 to 3.0 Ma, providing additional support for the newly emerging chronology based on orbital tuning. Slight discrepancies exist in the time interval from 3.0 to 3.6 Ma, suggesting several possibilities, including differences in the approaches of orbital tuning and the relatively low amplitude of delta18O variations in our record. However, even if the discrepancies are due to the relatively low amplitude of the isotope signals in our record at 3.0-3.6 Ma, our resultant timescale as a whole does not adversely affect our evaluation of the paleoclimatology and paleoceanography of the Indian Ocean, such as the evolution of the 100-, 41- and 23-kyr cycles, and variation of global ice volume and deepwater temperature during the past 3.6 m.y.

Barium intensities and tie points between holes of Leg 208

The Paleocene - Eocene thermal maximum (PETM) is one of the best known examples of a transient climate perturbation, associated with a brief, but intense, interval of global warming and a massive perturbation of the global carbon cycle from injection of isotopically light carbon into the ocean-atmosphere system. One key to quantifying the mass of carbon released, identifying the source(s), and understanding the ultimate fate of this carbon is to develop high-resolution age models. Two independent strategies have been employed, cycle stratigraphy and analysis of extraterrestrial Helium (HeET), both of which were first tested on Ocean Drilling Program (ODP) Site 690. Both methods are in agreement for the onset of the PETM and initial recovery, or the clay layer ("main body"), but seem to differ in the final recovery phase of the event above the clay layer, where the carbonate contents rise and carbon isotope values return toward background values. Here we present a state-of-the-art age model for the PETM derived from a new orbital chronology developed with cycle stratigraphic records from sites drilled during ODP Leg 208 (Walvis Ridge, Southeastern Atlantic) integrated with published records from Site 690 (Weddell Sea, Southern Ocean, ODP Leg 113). During Leg 208, five Paleocene - Eocene (P-E) boundary sections (Sites 1262 to 1267) were recovered in multiple holes over a depth transect of more than 2200 m at the Walvis Ridge yielding the first stratigraphically complete P-E deep-sea sequence with moderate to relatively high sedimentation rates (1 to 3 cm/kyr). A detailed chronology was developed with non-destructive X-ray fluorescence (XRF) core scanning records on the scale of precession cycles, with a total duration of the PETM now estimated to be ~ 170 kyr. The revised cycle stratigraphic record confirms original estimates for the duration of the onset and initial recovery, but suggests a new duration for the final recovery that is intermediate to the previous estimates by cycle stratigraphy and HeET.

Stable isotopes, CaCO3, opal and terrigeneous of ODP Site 108-663 (Table A1)

High- and low-latitude forcing of terrestrial African paleoclimate variability is demonstrated using 900 ka eolian and biogenic component records from Ocean Drilling Program site 663 in the eastern equatorial Atlantic. Terrigenous (eolian dust) and phytolith (savannah grass cuticle) accumulation rate records vary predominantly at 100 and 41 kyr periodicities and spectral phase estimates implicate high-latitude forcing. The abundance of freshwater diatoms (Melosira) transported from dry African lake beds varies coherently at 23-19 kyr orbital periodicities and at a phasing which implicates low-latitude precessional monsoon forcing. Modeling studies demonstrate that African climate is sensitive to both high- and low-latitude boundary conditions. African monsoon intensity is modulated by direct insolation variations due to orbital precession, whereas remote high-latitude forcing can be related to cool North Atlantic sea surface temperatures (SSTs) which promote African aridity and enhance dust-transporting wind speeds. The site 663 terrigenous and phytolith records covary with North Atlantic SST variability at 41 °N (site 607). We suggest that Pleistocene African climate has responded to both high-latitude North Atlantic SST variability as well as low-latitude precessional monsoon forcing; the high-latitude influence dominates the sedimentary record. Prior to circa 2.4 Ma, terrigenous variations occurred primarily at precessional periodicities (23-19 kyr), indicating that African climate was largely controlled by low-latitude insolation variations prior to the onset of high-amplitude glacial-interglacial climate change.

Radiocarbon datings and temperature reconstructions of sediment core GeoB10038-4 off SW Sumatra

Surface and thermocline conditions of the eastern tropical Indian Ocean were reconstructed through the past glacial-interglacial cycle by using Mg/Ca and alkenone-paleothermometry, stable oxygen isotopes of calcite and seawater, and terrigenous fraction performed on sediment core GeoB 10038-4 off SW Sumatra (~6°S, 103°E, 1819 m water depth). Results show that annual mean surface and thermocline temperatures varied differently and independently, and suggest that surface temperatures have been responding to southern high-latitude climate, whereas the more variable thermocline temperatures were remotely controlled by changes in the thermocline temperatures of the North Indian Ocean. Except for glacial terminations, salinity proxies indicate that changing intensities of the boreal summer monsoon did not considerably affect annual mean conditions off Sumatra during the past 133,000 years. Our results do not show a glacial-interglacial pattern in the thermocline conditions and reject a linear response of the tropical Indian Ocean thermocline to mid- and high-latitude climate change. Alkenone-based surface temperature estimates varied in line with the terrigenous fraction of the sediment and the East Asian winter monsoon proxy records at the precession band suggestive of monsoon (sea level) to be the dominant control on alkenone temperatures in the eastern tropical Indian Ocean on sub-orbital (glacial-interglacial) timescales.

(Table 1) Stable carbon isotope record and total nitrogen concentration in ODP Hole 175-1085B sediments

Variations in the strength of coastal upwelling in the South East Atlantic Ocean and summer monsoonal rains over South Africa are controlled by the regional atmospheric circulation regime. Although information about these parameters exists for the last glacial period, little detailed information exists for older time periods. New information from ODP Site 1085 for Marine Isotope Stages (MIS) 12-10 shows that glacial-interglacial productivity trends linked to upwelling variability followed a pattern similar to the last glacial cycle, with maximums shortly before glacial maxima, and minimums shortly before glacial terminations. During the MIS-11/10 transition, several periodic oscillations in productivity and monsoonal proxies are best explained by southwards shifts in the southern sub-tropical high-pressure cells followed by abrupt northwards shifts. Comparison to coeval sea-surface temperature measurements suggests that these monsoonal cycles were tightly coupled to anti-phased hemispheric climate change, with an intensified summer monsoon during periods of Northern (Southern) Hemisphere cooling (warming). The timing of these events suggests a pacing by insolation over precession periods. A lack of similar regional circulation shifts during the MIS-13/12 transition is likely due to the large equatorwards shift in the tropical convection zone that occurred during this extreme glaciation.

Benthic foraminiferal stable isotope stratigraphy of ODP Site 138-846 in the tropical Pacific Ocean

A stable-isotope stratigraphy at Site 846 (tropical Pacific, 3°06'S, 90°49'W, 3307 m water depth), based on the benthic foraminifers Cibicides wuellerstorfi and Uvigerina peregrina, yields a high-resolution record of deep-sea delta18O and delta13C over the past 1.8 Ma, with an average sampling interval of 3 k.y. Variance in the delta18O and delta13C records is concentrated in the well-known orbital periods of 100, 41, and 23 k.y. In the 100-k.y. band, both isotopic signals grow from relatively low amplitudes prior to 1.2 Ma, to high amplitudes in the late Quaternary since 0.7 Ma. The amplitude of delta18O and especially of delta13C decreases in the 41-k.y. band as it grows in the 100-k.y. band, consistent with a transfer of energy into an orbitally-paced internal oscillation. A weak 30-k.y. rhythm, present in both delta18O and delta13C, may reflect nonlinear interaction between the 41-k.y. and 100-k.y. bands in the evolving climate system. In the 23-k.y. and 19-k.y. bands associated with orbital precession, delta18O and delta13C are not coherent with each other on long time scales, and do not evolve like the 100-k.y. and 41-k.y. bands. This suggests that the source of the growing 100-k.y. oscillation is not a nonlinear response to precession, in contrast to predictions of some climate models. Sedimentation rates at this site also vary with a strong 100-k.y. cycle. Unlike the isotope records, the amplitude of 100-k.y. variations in sedimentation rate is relatively constant over the past 1.8 Ma, ranging from about 15 to 70 m/m.y. Prior to 0.9 Ma, sedimentation rates co-vary with orbital eccentricity, rather than with global climate as reflected by delta18O or delta13C. A source of this 100-k.y. cycle of sedimentation rate in the absence of similar ice volume fluctuations may be precessional heating of equatorial land masses, which in an energy balance climate model drives variations of monsoonal climates with a 100-k.y. rhythm. For the interval younger than 0.9 Ma, high sedimentation rates in the 100-k.y. band are consistently associated with glacial stages. This change of pattern suggests that when the amplitude of glacial cycles become large enough, their global effects overpower a local monsoon-driven variation in sedimentation rate at Site 846.

Physical properties of four profiles from the Triassic in Germany

The Norian Steinmergel-Keuper (SMK) represents a low-latitude cyclically-bedded playa system of the Mid-German Basin. We investigated a drilling site (core Morsleben) and sections from marginal positions. Dolomite/red mudstone beds form rhythmic alternations that were associated with varying monsoon activity. Hence, low K/Al ratios of dolomite beds suggest increased chemical weathering of the crystalline hinterland and therefore increased monsoonal rainfall. High K/Al ratios in red mudstone beds reflect increased physical weathering of the hinterlands during dryer periods. Dolomite layers reflect the lake stage (maximum monsoon) while red mudstones indicate the dry phase (minimum monsoon) of the playa cycle. We distinguished five major types of cyclic facies alternations, representing specific facies zones in the playa system. We have implemented spectrophotometry as a tool for high-resolution cyclostratigraphy. The dense sampling increment (up to 1 cm) allows for the recognition of all orbital frequencies. Sediment colour profiles reveal striking hierarchical cycles from semi-precession (SP, 99 kyr) over precession (P, 19.8 kyr) and obliquity (O, 36 kyr) to eccentricity (E1-2 109 kyr; E3, 413 kyr). A significant about 2 Myr-signal is attributed to the longer-term eccentricity E4. One monsoonal (precession) cycle includes two carbonate precipitation events. We propose that stratified mudstone and red mudstone are associated with maximum and minimum monsoon during the transition of the solstices in perihelion and aphelion, respectively. The two carbonate precipitation events were most likely created when equinoxes were in perihelion and aphelion, respectively. A sedimentary semi-precession response cycle is a novel finding for the Norian strata. The obliquity signal is attributed to incoming atmospheric moisture from the northeast of the SMK basin. The E4 cycle controls lake-level changes over long times. Apparently, E4 is responsible whether or not a threshold value is crossed. Bundles of 109 kyr and 413 kyr in red mudstones suggest a dry system with reduced monsoonal activity. In contrast, humid periods reveal thick layers of dolomite beds, indicating that during those intervals the monsoonal activity was strong enough to prevent the playa system from drying out completely.

Age determination of sediment cores of the California Current System

Alkenone sea surface temperature (SST) records were generated from the Ocean Drilling Program's (ODP) Sites 1014 and 1016 to examine the response of the California Current System to global climate change during the last 136 ka. The temperature differences between these sites (Delta SST(NEP)=SST(ODP1014)-SST(ODP1016)) reflected the intensity of the California Current and varied between 0.4 and 6.1 °C. A high Delta SST(NEP) (weaker California Current) was found for late marine isotope stage (MIS) 2 and early MIS 5e, while a low Delta SST(NEP) (stronger California Current) was detected for mid-MIS 5e and MIS 1. Spectral analysis indicated that this variation pattern dominated 23- (precession) and 30-ka periods. Comparison of the Delta SST(NEP) and SST based on data from core MD01-2421 at the Japan margin revealed anti-phase variation; the high Delta SST(NEP) (weakening of the California Current) corresponded to the low SST at the Japan margin (the southward displacement of the NW Pacific subarctic boundary), and vice versa. This variation was synchronous with a model prediction of the tropical El Niño-Southern Oscillation behavior. These findings suggest that the intensity of the North Pacific High varied in response to precessional forcing, and also that the response has been linked with the changes of tropical ocean-atmosphere interactions.

Iberian Margin 420 kyr geochemical and color variations record

Here we report 420 kyr long records of sediment geochemical and color variations from the southwestern Iberian Margin. We synchronized the Iberian Margin sediment record to Antarctic ice cores and speleothem records on millennial time scales and investigated the phase responses relative to orbital forcing of multiple proxy records available from these cores. Iberian Margin sediments contain strong precession power. Sediment "redness" (a* and 570-560 nm) and the ratio of long-chain alcohols to n-alkanes (C26OH/(C26OH + C29)) are highly coherent and in-phase with precession. Redder layers and more oxidizing conditions (low alcohol ratio) occur near precession minima (summer insolation maxima). We suggest these proxies respond rapidly to low-latitude insolation forcing by wind-driven processes (e.g., dust transport, upwelling, precipitation). Most Iberian Margin sediment parameters lag obliquity maxima by 7-8 ka, indicating a consistent linear response to insolation forcing at obliquity frequencies driven mainly by high-latitude processes. Although the lengths of the time series are short (420 ka) for detecting 100 kyr eccentricity cycles, the phase relationships support those obtained by Shackleton []. Antarctic temperature and the Iberian Margin alcohol ratios (C26OH/(C26OH + C29)) lead eccentricity maxima by 6 kyr, with lower ratios (increased oxygenation) occurring at eccentricity maxima. CO2, CH4, and Iberian SST are nearly in phase with eccentricity, and minimum ice volume (as inferred from Pacific d18Oseawater) lags eccentricity maxima by 10 kyr. The phase relationships derived in this study continue to support a potential role of the Earth's carbon cycle in contributing to the 100 kyr cycle.

Alkenone variations, total organic carbon and stable oxygen isotope ratios of the western Arabian Sea

We here present records of total organic carbon (TOC) and C37 alkenones, used as indicators for past primary productivity, from the western (WAS) and eastern Arabian Sea (EAS). New data from an open ocean site of the WAS upwelling area are compared with similar records from Ocean Drilling Program (ODP) Site 723 from the continental margin off Oman and MD 900963 from the EAS. These records together with other proxies used to reconstruct upwelling intensity, indicate periods of high productivity in tune with precessional forcing. On the basis of their phase relationship to boreal summer insolation they can be divided into three groups: in the WAS differences between monsoonal proxies (1) and productivity (2) document a combined signal of moderate SW monsoon winds and of strengthened and prolonged NE monsoon winds, whereas in the EAS phasing indicates maximum productivity (3) at times of stronger NE monsoon winds associated with precession-related maxima in ice volume.