Stable carbon and oxygen isotope compositions of Eocene benthic foraminifera from ODP Site 207-1258

'Hyperthermals' are intervals of rapid, pronounced global warming known from six episodes within the Palaeocene and Eocene epochs (~65-34 million years (Myr) ago) (Zachos et al., 2005, doi:10.1126/science.1109004; 2008, doi:10.1038/nature06588; Roehl et al., 2007, doi:10.1029/2007GC001784; Thomas et al., 2000; Cramer et al., 2003, doi:10.1029/2003PA000909; Lourens et al., 2005, doi:10.1038/nature03814; Petrizzo, 2005, doi:10.2973/odp.proc.sr.198.102.2005; Sexton et al., 2006, doi:10.1029/2005PA001253; Westerhold et al., 2007, doi:10.1029/2006PA001322; Edgar et al., 2007, doi:10.1038/nature06053; Nicolo et al., 2007, doi:10.1130/G23648A.1; Quillévéré et al., 2008, doi:10.1016/j.epsl.2007.10.040; Stap et al., 2010, doi:10.1130/G30777.1). The most extreme hyperthermal was the 170 thousand year (kyr) interval (Roehl et al., 2007) of 5-7 °C global warming (Zachos et al., 2008) during the Palaeocene-Eocene Thermal Maximum (PETM, 56 Myr ago). The PETM is widely attributed to massive release of greenhouse gases from buried sedimentary carbon reservoirs (Zachos et al., 2005; 2008; Lourenbs et al., 2005; Nicolo et al., 2007; Dickens et al., 1995, doi:10.1029/95PA02087; Dickens, 2000; 2003, doi:10.1016/S0012-821X(03)00325-X; Panchuk et al., 2008, doi:10.1130/G24474A.1) and other, comparatively modest, hyperthermals have also been linked to the release of sedimentary carbon (Zachos et al., 2008, Lourens et al., 2005; Nicolo et al., 2007; Dickens, 2003; Panchuk et al., 2003). Here we show, using new 2.4-Myr-long Eocene deep ocean records, that the comparatively modest hyperthermals are much more numerous than previously documented, paced by the eccentricity of Earth's orbit and have shorter durations (~40 kyr) and more rapid recovery phases than the PETM. These findings point to the operation of fundamentally different forcing and feedback mechanisms than for the PETM, involving redistribution of carbon among Earth's readily exchangeable surface reservoirs rather than carbon exhumation from, and subsequent burial back into, the sedimentary reservoir. Specifically, we interpret our records to indicate repeated, large-scale releases of dissolved organic carbon (at least 1,600 gigatonnes) from the ocean by ventilation (strengthened oxidation) of the ocean interior. The rapid recovery of the carbon cycle following each Eocene hyperthermal strongly suggests that carbon was resequestered by the ocean, rather than the much slower process of silicate rock weathering proposed for the PETM (Zachos et al., 2005; 2003). Our findings suggest that these pronounced climate warming events were driven not by repeated releases of carbon from buried sedimentary sources (Zachos et al., 2008, Lourens et al., 2005; Nicolo et al., 2007; Dickens, 2003; Panchuk et al., 2003) but, rather, by patterns of surficial carbon redistribution familiar from younger intervals of Earth history.

Alkenone analyses of sediment cores from the western Arabian Sea

Records of total organic carbon (TOC) and C37 alkenones were used as indicators for past primary productivity in the western and eastern Arabian Sea. Data from GeoB 3005, an open ocean site in the western Arabian Sea upwelling area, are compared with similar records of GeoB 3007 from the Owen Ridge, Ocean Drilling Program (ODP) Site 723 from the continental margin off Oman and MD 900963 from the eastern Arabian Sea. TOC/C37 alkenone records together with other proxies used to reconstruct upwelling intensity, indicate periods of high productivity in tune with precessional forcing all over the Arabian Sea. Based on their phase-relationship to variations in boreal summer insolation they can be divided into three groups. In the western Arabian Sea the precession-related phasing is different between productivity proxies and those for summer monsoon wind strength and upwelling intensity. TOC and C37 alkenone records from the western Arabian Sea lag the other monsoonal indicators by about 5 kyr, but lead productivity indicators from the eastern Arabian Sea by 3 kyr. Based on the differences in phase relationships associated with the precessional cycling between productivity and monsoonal proxies in the western Arabian Sea it is proposed that the TOC/C37 alkenone signal in the western Arabian Sea document a combined signal of moderate SW monsoon winds and of strengthened and prolonged NE monsoon winds. In the eastern Arabian Sea the phasing hints to coincidence between maximum productivity and stronger NE monsoon winds associated with precession-related maxima in ice volume. In contrast, variations in paleoproductivity at site GeoB 3007 from the Owen Ridge indicate productivity maxima during glacial substages 8.2, 6.2 and 2.2, whereas precessionrelated changes are of only minor importance at this location. The results of frequency analyses confirm that productivity at site GeoB 3007 responds predominantly to glacialinterglacial climate changes, while site GeoB 3005 from the open ocean upwelling region near the Gulf of Aden is dominated by precessional insolation. A possible explanation for the pattern revealed at the Owen Ridge is the periodic NW-SE displacement of the Findlater Jet axis, which separates the region of open ocean upwelling to the northwest from downwelling to the southeast ofthe jet. The carbon isotopes of planktic foraminifera reflect nutrient related d13C variations of dissolved inorganic carbon. The difference between the planktic foraminifera Globigerinoides ruber (w), living in the upper 50 m of the water column, and the deeper Iiving Neogloboquadrina dutertrei (Delta d13Cr-d) of core GeoB 3005 displays nutrient variations in the upwelling area near the Gulf of Aden. The results of cross-spectral analyses between Deltad13Cr-d of GeoB 3005 and proxies for SW monsoon intensity indicate, too, a dissociation of productivity from monsoonal upwelling intensity. Instead, productivity depends mainly on the availability of nutrients, while upwelling intensity of sub-surface water masses seems to be of only secondary importance. Additionally, sea surface temperatures (SSTs) were reconstructed using the unsaturation ratio of C37 alkenones. Alkenone SSTs reflect annual mean temperatures rather than explicitly the season of upwelling. This is evident from alkenone SSTs in a transect of surface sediments extending from the inner Gulf of Aden into the western Arabian Sea. The alkenone-derived SST records of GeoB 3005 from the open ocean upwelling region near the Gulf of Aden and GeoB 3007 from the Owen Ridge reveal similar variations with high SSTs during interglacial and low SSTs during glacial periods. The glacial marine oxygen isotope stage (MIS) 6 remains relatively warm and was not as cold as MIS 3 to 4 and 8 according to the alkenone SST. Similar variation-patterns were reconstructed in the coastal upwelling area off Oman for ODP Site 723 as weIl as in the eastern Arabian Sea for MD 900963, where upwelling is not as pronounced as in the western Arabian Sea. Spectral-analyses indicate that SST changes are in good agreement with the modulation of low-latitude precessional insolation changes by eccentricity. However, they do not show the pronounced cydicity in the precessional frequency band, which is characteristic for variations in paleoproductivity. Although the overall variation pattern is very similar, a dose comparison between the western (GeoB 3005) and the eastern Arabian Sea (MD 900963) shows larger differences between both sites during cold intervals than during periods of warm SSTs. This is attributed to a more effective cooling of surface waters in the western Arabian Sea by prolonged NE monsoon winds during times of expanded Northern Hemisphere ice-sheets, thereby lowering the annual mean SSTs stronger than in the eastern Arabian Sea.

Geochemistry of Eocene/Oligocene sediments of ODP Hole 154-925A

Paleoproductivity, nutrient burial, and carbon cycling were investigated across the Eocene/Oligocene (E/O) boundary (begin to end; 36.9-32.7 Ma at ~40 kyr resolution, timescale of Shackleton et al. (1999, doi:10.1098/rsta.1999.0407) at Ocean Drilling Program Site 925 on the Ceara Rise in the western equatorial Atlantic (3040 m present water depth; 748.26-850.70 mbsf). Downcore bulk sediment records of biogenic barium, total reactive phosphorus, biogenic silica, and calcium carbonate are interpreted to represent export production, net nutrient burial, biogenic opal production, and inorganic carbon burial, respectively. The global positive excursion in d13C subsequent to the E/O boundary is recorded at Site 925. Export production appears to have been externally forced by orbital parameters at eccentricity frequencies during the study interval, based on spectral analysis of the biogenic barium and reactive phosphorus records. Biogenic silica production or preservation increased after the Eocene/Oligocene boundary to a higher baseline, although overall productivity and nutrient burial did not increase, based on barium and reactive phosphorus records. Thus, although absolute production did not increase at this site, a shift in relative abundance of siliceous versus carbonate productivity may have resulted in a change in relative organic carbon burial. This may have contributed to the positive excursion in global oceanic d13C subsequent to the Eocene/Oligocene boundary, although the silica maximum persists after the carbon isotope excursion ends.

Planktonic foraminiferal preservation in Plio/Pleistocene sediments from the western equatorial Atlantic and Caribbean

Records of mean sortable silt and planktonic foraminiferal preservation from the Ceará Rise (western equatorial Atlantic) and from the Caribbean are presented to analyze the Pliocene (3.5-2.2 Ma) to Pleistocene (1.6-0.3 Ma) evolution of near-bottom current strength and the carbonate corrosiveness of deep water. During the mid-Pleistocene climate transition (~1 Ma) a drastic decrease in glacial bottom current strength and an increase in carbonate corrosiveness is registered, demonstrating a substantial decrease in the glacial contribution of the Lower North Atlantic Deep Water (LNADW) to the Atlantic Ocean. Also, an increased sensitivity to eccentricity orbital forcing is registered after the MPT. By contrast, carbonate preservation increases considerably in the deep Caribbean in response to a strong and persistent stable contribution of Upper North Atlantic Deep Water (UNADW). We found evidence for the strongest and most stable circulation within the LNADW cell during the Northern Hemisphere cooling period between ~3.2 and 2.75 Ma. This is in agreement with the 'superconveyor model' which postulates that the highest NADW production took place prior to ~2.7 Ma. A considerable decrease in bottom current strength and planktonic foraminiferal preservation is observed synchronous with the first occurrence of large-scale continental ice sheets in the Northern Hemisphere. This documents the final termination of the 'superconveyor' at ca. 2.75 Ma. However, our data do not support a 'superconveyor' in the interval between 3.5 and 3.2 Ma when high-amplitude fluctuations in bottom current flow and preservation in planktonic foraminifera are observed. Because of the great sensitivity of NADW production to changes in surface water salinity, we assume that the high-amplitude fluctuations of LNADW circulation prior to ~3.2 Ma are linked to changes in the Atlantic salinity budget. After 2.75 Ma they are primarily controlled by ice-sheet forcing. In contrast to the stepwise deterioration of planktonic foraminiferal preservation in the western deep Atlantic, a trend toward better preservation from the Pliocene to Pleistocene is observed in the deep Caribbean. This indicates a long-term increase in the contribution of UNADW to the Atlantic Ocean.

Tropical Atlantic SST during the mid-Pleistocene transition (MPT) of ODP Hole 175-1077B

We compare a new mid-Pleistocene sea surface temperature (SST) record from the eastern tropical Atlantic to changes in continental ice volume, orbital insolation, Atlantic deepwater ventilation, and Southern Ocean front positions to resolve forcing mechanisms of tropical Atlantic SST during the mid-Pleistocene transition (MPT). At the onset of the MPT, a strong tropical cooling occurred. The change from a obliquity- to a eccentricity-dominated cyclicity in the tropical SST took place at about 650 kyr BP. In orbital cycles, tropical SST changes significantly preceded continental ice-volume changes but were in phase with movements of Southern Ocean fronts. After the onset of large-amplitude 100-kyr variations, additional late glacial warming in the eastern tropical Atlantic was caused by enhanced return flow of warm waters from the western Atlantic driven by strong trade winds. Pronounced 80-kyr variations in tropical SST occurred during the MPT, in phase with and likely directly forced by transitional continental ice-volume variations. During the MPT, a prominent anomalous long-term tropical warming occurred, likely generated by extremely northward displaced Southern Ocean fronts. While the overall pattern of global climate variability during the MPT was determined by changes in mean state and frequency of continental ice volume variations, tropical Atlantic SST variations were primarily driven by early changes in Subantarctic sea-ice extent and coupled Southern Ocean frontal positions.

Parameters of the Earths orbit for the last 5 Million years in 1 kyr resolution

New values for the astronomical parameters of the Earth's orbit and rotation (eccentricity, obliquity and precession) are proposed for paleoclimatic research related to the Late Miocene, the Pliocene and the Quaternary. They have been obtained from a numerical solution of the Lagrangian system of the planetary point masses and from an analytical solution of the Poisson equations of the Earth-Moon system. The analytical expansion developed in this paper allows the direct determination of the main frequencies with their phase and amplitude. Numerical and analytical comparisons with the former astronomical solution BER78 are performed so that the accuracy and the interval of time over which the new solution is valid can be estimated. The corresponding insolation values have also been computed and compared to the former ones. This analysis leads to the conclusion that the new values are expected to be reliable over the last 5 Ma in the time domain and at least over the last 10 Ma in the frequency domain.

Astronomical time scale for Paleocene and early Eocene sediments

The astronomical-tuned time scale is rapidly extended into the Paleogene but, due to the existence of an Eocene gap, different tuning options had to be presented for the Paleocene. These options differ both in number and tuning of ~405-kyr eccentricity related cycles and are only partially consistent with recalculated 40Ar/39Ar constraints for the Cretaceous/Paleogene (K/Pg) and Paleocene/Eocene (P/E) boundaries. In this paper, we evaluate the cyclostratigraphic interpretation of records from ODP Leg 198 and 208 sites, and the Zumaia section to solve the problem of the different tuning options. We found that the interval between the K/Pg boundary and the early Late Paleocene biotic event (ELPE) comprises 17 instead of 16 * ~405-kyr eccentricity related cycles as previously proposed, while the entire Paleocene contains 25 * ~405-kyr cycles. Starting from 40Ar/39Ar age constraints for the K/Pg boundary, a new tuning to 405-kyr eccentricity is presented for the Paleocene and earliest Eocene, which results in ages of ~66.0 and ~ 56.0 Ma for the K/Pg and P/E boundaries, respectively. This tuning introduces considerable differences in age for a number of nannofossil events at ODP Sites 1209 and 1262 in the interval between 61 and 63 Ma, but eliminates large and abrupt changes in the seafloor spreading rate. The tuning seems further consistent with recalculated 40Ar/39Ar ages for ash layer -17 of early Eocene age. However, despite this apparent consistency with existing radio-isotopic constraints, an alternative 405-kyr younger or, less likely, older tuning cannot be excluded at this stage.

Short-term spectral analysis of ODP Hole 114-704B

Short-term spectral analysis was carried out on geochemical logging data from ODP Site 704. The FFT was used to compute the amplitude spectra of short-term overlapping segments to produce depth-period-amplitude spectrograms of the logging data. The spectrograms provided a means of evaluating the significance of the observed periodic components. The periodic components that were consistently present and prominent across a given record interval were considered to be significant. Changes in the spectrogram characteristics seem to reflect changes in either lithology, sedimentation rates, or hiatuses and may therefore provide useful information to aid in stratigraphic and paleoenvironmental studies. The dominant periodicity during the late Pleistocene and Brunhes Chron (0.97 to 0.47 Ma) was determined to be > 100,000 yr whereas the upper Matuyama Chron was dominated by the 41,000-yr periodicity. These periodicities suggest that the sedimentation patterns within the upper Matuyama Chron (0.98-1.78 Ma) were influenced by the Milankovitch obliquity cycle and those within the latest Matuyama-Brunhes Chron (<0.98 Ma) by the eccentricity cycle. The Brunhes/Matuyama boundary therefore represents a major discontinuity. Periodicities observed within the lower Matuyama and the upper Gauss Chron did not correlate with any of the periodicities within the Milankovitch frequency bands.

Age model of sediment core V20-120 (Table 2)

Time series analyses of atmospheric and oceanic variables in a late Pleistocene record from the northwest Pacific show the complex relationship of the response of various segments of the climate system to changes in the earth's orbit. Most variance spectra of time series from this subarctic record contain frequency peaks with periods corresponding to at least one of the major orbital components of eccentricity, obliquity, or precession. Although the radiolarian faunal (water mass) assemblages have prominent spectral peaks with 41,000-year periods which are coherent with obliquity at this frequency, only the Transitional faunal assemblage contains variance focused at a frequency corresponding to the 100,000-year period of eccentricity. Three of these faunal time series also show variance concentrated at a frequency with a 20,000-year period. These three time series are not coherent at a 20,000-year frequency with either of the dominant spectral peaks of precession. They are coherent, however, with variations in the second harmonic of the obliquity cycle. Changes in obliquity apparently affect siliceous faunal abundances in the northwest Pacific region of this high-latitude site more than variations in eccentricity or precession. Maxima in the time series of quartz abundance occur during low values of eccentricity and high glacial ice volume. Because atmospheric winds serve as the major source of supply of quartz to the sediments at this site, these high quartz values reflect increased aridity at the source region. Except for short periods during interglacials, the characteristics of the surface waters in this region of the subarctic Pacific during much of the last 460,000 years were similar to those which exist today in the Sea of Okhotsk. The spectrum of winter sea surface temperature estimates, derived from siliceous microfaunal abundances, contains dominant frequency peaks at periods of 100,000, 41,000, and 23,000 years which are coherent with eccentricity, obliquity, and precession, respectively. Based upon the relationship of the Subarctic Front with the dominance of specific faunal asemblages, the front was positioned south of its present-day location throughout much of the late Pleistocene.

Milankovitch frequency peaks of ODP Site 113-693 (Table 3)

The first well logs collected below the Antarctic circle were obtained during Leg 113 at Site 693 on the Dronning Maud Land Margin (Antarctica) in the Weddell Sea. Gamma-ray, resistivity, and sonic logs were collected between 108.0 and 439.0 mbsf. The downhole logs show good agreement with the data collected from cores and provide a continuous measurement of the sedimentary record. These continuous log records show that the rather uniform Tertiary lithology seen in cores is characterized by high-frequency variability in the log data. Several thin hard streaks are identified, the largest of which coincides with a major Miocene hiatus. Associated with this hiatus is a change to lower illite content (and correspondingly lower gamma-ray counts) and to a significant increase in diatom content. Spectral analysis of the logs was performed on the lower Pliocene through upper Oligocene interval (108.0-343.0 mbsf). Between 108.0 and 245.0 mbsf, average sedimentation rates (50 and 26 m/m.y.) are high enough to show that variance is present in the orbital eccentricity (~95 k.y.) and obliquity (~41 k.y.) bands. Between 253.0 and 343.0 mbsf, the sedimentation rate (8 m/m.y.) is too low to resolve high frequency variations. The Milankovitch frequencies are best developed in the resistivity logs. Resistivity is responding to changes in porosity, which in these sediments is controlled by the abundance of biosiliceous sediments, particularly diatoms. The orbital forcing suggested by the Milankovitch frequencies may be influencing diatom productivity by inducing oscillations in upwelling, ice coverage, pack ice, and/or polynya. Although variations in diatom abundance were observed in the cores, they were not attributed to a Milankovitch signal, and therefore in this environment, downhole logs are an important contribution to the detection and understanding of orbitally influenced changes in sedimentation.

Massic magnetic susceptibility and CaCO3 content at the P. ohmi and T. hugii zones from the Rio Argos X.Ag-1 section, southeastern Spain

An orbital floating time scale of the Hauterivian-Barremian transition (Early Cretaceous) is proposed using high-resolution magnetic susceptibility measurements. Orbital tuning was performed on the Río Argos section (southeast Spain), the candidate for a Global boundary Stratotype Section and Point (GSSP) for the Hauterivian-Barremian transition. Spectral analyses of MS variations, coupled with the frequency ratio method, allow the recognition of precession, obliquity and eccentricity frequency bands. Orbitally-tuned magnetic susceptibility provides minimum durations for ammonite biozones. The durations of well-constrained ammonite zones are assessed at 0.78 myr for Pseudothurmannia ohmi (Late Hauterivian) and 0.57 myr for Taveraidiscus hugii (Early Barremian). These results are consistent with previous estimates from the other reference section (Angles, southeast France) and tend to show that the Río Argos section displays a complete succession for this time interval. They differ significantly from those proposed in the Geologic Time Scale 2008 and may help to improve the next compilation. The Faraoni Oceanic Anoxic Event, a key Early Cretaceous oceanographic perturbation occurring at the P. ohmi/P. catulloi subzone boundary has a duration estimated at 0.10-0.15 myr, which is similar to previous assessments.

Age models and sedimentation rates of sites in the tropical Atlantic

Late Pleistocene signals of calcium carbonate, organic carbon, and opaline silica concentration and accumulation are documented in a series of cores from a zonal/meridional/depth transect in the equatorial Atlantic Ocean to reconstruct the regional sedimentary history. Spectral analysis reveals that maxima and minima in biogenous sedimentation occur with glacial-interglacial cyclicity as a function of both (1) primary production at the sea surface modulated by orbitally forced variation in trade wind zonality and (2) destruction at the seafloor by variation in the chemical character of advected intermediate and deep water from high latitudes modulated by high-latitude ice volume. From these results a pattern emerges in which the relative proportion of signal variance from the productivity signal centered on the precessional (23 kyr) band decreases while that of the destruction signal centered on the obliquity (41 kyr) and eccentricity (100 kyr) periods increases below ~3600-m ocean depth.

Stable isotopes measured on Miocene benthic foraminifera from ODP sites 202-1237 and 184-1146

One of the most enigmatic features of Cenozoic long-term climate evolution is the long-lasting positive carbon-isotope excursion or "Monterey Excursion", which started during a period of global warmth after 16.9 Ma and ended at not, vert, similar 13.5 Ma, approximately 400 kyr after major expansion of the Antarctic ice-sheet. We present high-resolution (1-9 kyr) astronomically-tuned climate proxy records in two complete sedimentary successions from the northwestern and southeastern Pacific (ODP Sites 1146 and 1237), which shed new light on the middle Miocene carbon-isotope excursion and associated climatic transition over the interval 17.1-12.7 Ma. We recognize three distinct climate phases with different imprints of orbital variations into the climatic signals (1146 and 1237 d18O, d13C; 1237 XRF Fe, fraction > 63 µm): (1) climate optimum prior to 14.7 Ma characterized by minimum ice volume and prominent 100 and 400 kyr variability, (2) long-term cooling from 14.7 to 13.9 Ma, principally driven by obliquity and culminating with rapid cryosphere expansion and global cooling at the onset of the last and most pronounced d13C increase, (3) "Icehouse" mode after 13.9 Ma with distinct 100 kyr variability and improved ventilation of the deep Pacific. The "Monterey" carbon-isotope excursion (16.9-13.5 Ma) consists overall of nine 400 kyr cycles, which show high coherence with the long eccentricity period. Superposed on these low-frequency oscillations are high-frequency variations (100 kyr), which closely track the amplitude modulation of the short eccentricity period. In contrast to d13C, the d18O signal additionally shows significant power in the 41 kyr band, and the 1.2 Myr amplitude modulation of the obliquity cycle is clearly imprinted in the 1146 d18O signal. Our results suggest that eccentricity was a prime pacemaker of middle Miocene climate evolution through the modulation of long-term carbon budgets and that obliquity-paced changes in high-latitude seasonality favored the transition into the "Icehouse" climate.