Relative abundances of twenty-eight dominant benthic foraminifera species of DSDP Hole 23-219 (Table 1s)

Tropical climate is variable on astronomical time scale, driving changes in surface and deep-sea fauna during the Pliocene-Pleistocene. To understand these changes in the tropical Indian Ocean over the past 2.36 Myr, we quantitatively analyzed deep-sea benthic foraminifera and selected planktic foraminifera from >125 µm size fraction from Deep Sea Drilling Project Site 219. The data from Site 219 was combined with published foraminiferal and isotope data from Site 214, eastern Indian Ocean to determine the nature of changes. Factor and cluster analyses of the 28 highest-ranked species distinguished four biofacies, characterizing distinct deep-sea environmental settings. These biofacies have been named after their most dominant species such as Stilostomella lepidula-Pleurostomella alternans (Sl-Pa), Nuttallides umbonifer-Globocassidulina subglobosa (Nu-Gs), Oridorsalis umbonatus-Gavelinopsis lobatulus (Ou-Gl) and Epistominella exigua-Uvigerina hispido-costata (Ee-Uh) biofacies. Biofacies Sl-Pa ranges from ~2.36 to 0.55 Myr, biofacies Nu-Gs ranges from ~1.9 to 0.65 Myr, biofacies Ou-Gl ranges from ~1 to 0.35 Myr and biofacies Ee-Uh ranges from 1.1 to 0.25 Myr. The proxy record indicates fluctuating tropical environmental conditions such as oxygenation, surface productivity and organic food supply. These changes appear to have been driven by changes in monsoonal wind intensity related to glacial-interglacial cycles. A shift at ~1.2-0.9 Myr is observed in both the faunal and isotope records at Site 219, indicating a major increase in monsoon-induced productivity. This coincides with increased amplitude of glacial cycles, which appear to have influenced low latitude monsoonal climate as well as deep-sea conditions in the tropical Indian Ocean.

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.

Dating analysis from DSDP Hole 90-594

Date-32 is a fast and easily used computer program developed to date Quaternary deep-sea cores by associating variations in the earth's orbit with recurring oscillations in core properties, such as carbonate content or isotope composition. Starting with known top and bottom dates, distortions in the periodicities of the core properties due to varying sedimentation rates are realigned by fast Fourier analysis so as to maximise the spectral energy density at the orbital frequencies. This allows age interpolation to all parts of the core to an accuracy of 10 kyrs, or about 1.5% of the record duration for a typical Brunhes sequence. The influence of astronomical forcing is examined and the method is applied to provide preliminary dates in a high-resolution Brunhes record from DSDP Site 594 off southeastern New Zealand.

Paleomagnetic characterization and age of several ODP Leg 198 sites, Shatsky Rise

Shatsky Rise, a medium-sized large igneous province in the west Central Pacific Ocean, has three main topographic highs that preserve a thick sedimentary record from Cretaceous through Cenozoic. During Ocean Drilling Program (ODP) Leg 198 to Shatsky Rise, a total of ~768 m of late Miocene-Holocene sediments was recovered from six sites. Sites 1207 and 1208 were drilled on the Northern and Central Highs, respectively, and yielded expanded late Miocene-Holocene sequences. Sites 1209, 1210, 1211, and 1212 were drilled on the Southern High and yielded shorter sequences of similar age. Clearly interpretable magnetic stratigraphies were obtained from all sites using the shipboard pass-through magnetometer. These results were augmented using discrete sample cubes (7 cm**3) collected shipboard and measured postcruise. Miocene age sediments are separated by a hiatus from Oligocene, Eocene, and Cretaceous age sediments beneath. An astrochronological age model was developed for the six sites based on cycles observed in reflectance data, measured shipboard. This age model is in good agreement with published astrochronological polarity chron ages in the 1 to 6 Ma interval.

Natural remanent magnetization (NRM) from ODP Site 1263 covering magnetochron C20r and C21n and high-resolution bulk carbon isotope (d13C) records from Sites 702 and 1263

To explore cause and consequences of past climate change, very accurate age models such as those provided by the astronomical timescale (ATS) are needed. Beyond 40 million years the accuracy of the ATS critically depends on the correctness of orbital models and radioisotopic dating techniques. Discrepancies in the age dating of sedimentary successions and the lack of suitable records spanning the middle Eocene have prevented development of a continuous astronomically calibrated geological timescale for the entire Cenozoic Era. We now solve this problem by constructing an independent astrochronological stratigraphy based on Earth's stable 405 kyr eccentricity cycle between 41 and 48 million years ago (Ma) with new data from deep-sea sedimentary sequences in the South Atlantic Ocean. This new link completes the Paleogene astronomical timescale and confirms the intercalibration of radioisotopic and astronomical dating methods back through the Paleocene-Eocene Thermal Maximum (PETM, 55.930 Ma) and the Cretaceous-Paleogene boundary (66.022 Ma). Coupling of the Paleogene 405 kyr cyclostratigraphic frameworks across the middle Eocene further paves the way for extending the ATS into the Mesozoic.

Stable carbon and oxygen isotope compositions of bulk rock samples analyzed from the Zumaia section

Astronomical tuning of sedimentary records to precise orbital solutions has led to unprecedented resolution in the geological time scale. However, the construction of a consistent astronomical time scale for the Paleocene is controversial due to uncertainties in the recognition of the exact number of 405-kyr eccentricity cycles and accurate correlation between key records. Here, we present a new Danian integrated stratigraphic framework using the land-based Zumaia and Sopelana hemipelagic sections from the Basque Basin and deep-sea records drilled during Ocean Drilling Program (ODP) Legs 198 (Shatsky Rise, North Pacific) and 208 (Walvis Ridge, South Atlantic) that solves previous discrepancies. The new coherent stratigraphy utilises composite images from ODP cores, a new whole-rock d13C isotope record at Zumaia and new magnetostratigraphic data from Sopelana. We consistently observe 11 405-kyr eccentricity cycles in all studied Danian successions. We achieve a robust correlation of bioevents and stable isotope events between all studied sections at the ~100-kyr short-eccentricity level, a prerequisite for paleoclimatic interpretations. Comparison with and subsequent tuning of the records to the latest orbital solution La2011 provides astronomically calibrated ages of 66.022 ± 0.040 Ma and 61.607 ± 0.040 Ma for the Cretaceous-Paleogene (K-Pg) and Danian-Selandian 105 (D-S) boundaries respectively. Low sedimentation rates appear common in all records in the mid-Danian interval, including conspicuous condensed intervals in the oceanic records that in the past have hampered the proper identification of cycles. The comprehensive interbasinal approach applied here reveals pitfalls in time scale construction, filtering techniques in particular, and indicates that some caution and scrutiny has to be applied when building orbital chronologies. Finally, the Zumaia section, already hosting the Selandian Global Boundary Stratotype Section and Point (GSSP), could serve as the global Danian unit stratotype in the future.

Calculated sea level records and deep ocean temperatures from DSDP and ODP sites

During the mid-Pleistocene transition the dominant 41 ka periodicity of glacial cycles transitioned to a quasi-100 ka periodicity for reasons not yet known. This study investigates the potential role of deep ocean hydrography by examining oxygen isotope ratios in benthic foraminifera. Oxygen isotope records from the Atlantic, Pacific and Indian Ocean basins are separated into their ice volume and local temperature/hydrography components using a piece-wise linear transfer function and a temperature calibration. Although our method has certain limitations, the deep ocean hydrography reconstructions show that glacial deep ocean temperatures approached freezing point as the mid-Pleistocene transition progressed. Further analysis suggests that water mass reorganisation could have been responsible for these temperature changes, leading to such stable conditions in the deep ocean that some obliquity cycles were skipped until precessional forcing triggered deglaciation, creating the apparent quasi-100 ka pattern. This study supports previous work that suggests multiples of obliquity cycles dominate the quasi-100 ka glacial cycles with precession components driving deglaciations.

Astronomically modulated Neogene sediment records of ODP Site 184-1148

General global cooling over the Neogene has been modulated by changes in Earth's orbital parameters. Investigations of deep-sea sediment sequences show that various orbital cycles can dominate climate records for different latitudes or for different time intervals. However, a comprehensive understanding of astronomical imprints over the entire Neogene has been elusive because of the general absence of long, continuous records extending beyond the Pliocene. We present benthic foraminiferal d18O and d13C records over the past 23 Ma at Ocean Drilling Program Site 1148 in the northern South China Sea and construct an astronomically tuned timescale (TJ08) for these records based on natural gamma radiation and color reflectance data at this site. Our results show that a 41 ka cycle has dominated sediment records at this location over the Neogene, displaying a linear response to orbital forcing. A 100 ka cycle has also been significant. However, it is correlated nonlinearly with Earth's orbital variations at the 100 ka band. The sediment records also display a prominent 405 ka cycle. Although this cycle was coherent with orbital forcing during the Oligocene and the early Miocene, it was not coherent with Earth's orbital variations at the 405 ka band over the whole Neogene. Amplification of Northern Hemisphere and Southern Hemisphere glaciation since the middle Miocene may be responsible for this change in sedimentary response. Our benthic foraminifera d18O and d13C records further exhibit amplitude variations with longer periods of 600, 1000, 1200, and 2400 ka. Apparently, these cycles are nonlinear responses to insolation forcing.