Foraminiferal isotope ratios from IODP Site U1313

Large-amplitude millennial-scale climate oscillations have been identified in late Pleistocene climate archives from around the world. These oscillations appear to be of larger amplitude during times of enlarged ice sheets. This observation suggests the existence of a relationship between large-amplitude millennial variations in climate and extreme glacial conditions and therefore that the emergence of millennial-scale climate variability may be linked to the Pliocene intensification of northern hemisphere glaciation (iNHG). Here we test this hypothesis using new late Pliocene high-resolution (ab. 400 year) records of ice-rafted debris deposition and stable isotopes in planktic foraminiferal calcite (Globigerinoides ruber) generated from Integrated Ocean Drilling Program Site U1313 in the subpolar North Atlantic (a reoccupation of the classic Deep Sea Drilling Project Site 607). Our records span marine oxygen isotope stages (MIS) 103-95 (ab. 2600 to 2400 ka), the first interval during iNHG (ab. 3.5 to 2.5 Ma) in which large-amplitude glacial-interglacial cycles and inferred sea level changes occur. Our records reveal small-amplitude variability at periodicities of ab. 1.8 to 6.2 kyr that prevails regardless of (inter)glacial state with no significant amplification during the glacials MIS 100, 98, and 96. These findings imply that the threshold for the amplification of such variability to the proportions seen in the marine archive of the last glacial was not crossed during the late Pliocene and, in view of all available data, likely not until the Mid-Pleistocene Transition.

Geochemistry and physics on sediment cores from the Peru Basin

We studied preservation/dissolution cycles and paleoproductivity in eight sediment cores from the Peru Basin south of the highly productive surface waters of the eastern equatorial Pacific. Stratigraphy is based on stable oxygen isotopes and on combined magnetostratigraphy and biostratigraphy. Sediment cores which span the last 8 m.y., were retrieved during cruise 79 with RV SONNE close to the carbonate compensation depth (CCD). In general, sediments show Pacific-type carbonate cycles. We interpret a pronounced carbonate peak between 6 and 7 Ma as the result of a western and northern extension of the highly productive Peru Current. Decreased carbonate contents from the late Miocene to the late Pliocene might be associated with a slow contraction of the latitudinal extent of the high-productivity belt north of the study areas. During the Pliocene, carbonate variations showed 400 kyr cycles indicating the growth and decay of ice sheets, which should have been associated with pulsations of the Antarctic ice cap. An abrupt collapse of the carbonate system occurred at 2.4 Ma. Higher frequency variations of the carbonate record indicate the major increase of the northern hemisphere glaciation. During the Quaternary, carbonate fluxes are high during glacials and low during interglacials. Large amplitude variations with long broad minima and maxima, associated with small migrations of the lysocline and the CCD (< 200 m), are indicative of the preservation/dissolution history in the Peru Basin. During the early Pleistocene, climatic forcing by the 41 kyr obliquity cycle is not observed in the carbonate record. During the last 800 kyr, variability in the carbonate record was dominated by the 100 kyr eccentricity cycle. Fluxes of biogenic material (calcium carbonate, organic carbon, opal, and barium) were greatest during glacials, which imply higher productivity and export production of the Peru Current during cold climatic periods. Dissolution was greatest during interglacials as inferred from the relatively poor preservation of planktonic foraminifera and from the low accumulation rate of carbonate. After the Mid-Brunhes Event (400 ka), we observe a plateaulike shift to enhanced dissolution and to intensified productivity.

(Table S3) Benthic foraminiferal Mg/Ca, Li/Ca and stable isotope ratios of ODP Site 122-761

The sensitivities of benthic foraminiferal Mg/Ca and Li/Ca to bottom water temperature and carbonate saturation state have recently been assessed. Here we present a new approach that uses paired Mg/Ca and Li/Ca records to calculate simultaneous changes in temperature and saturation state. Using previously published records, we first use this approach to document a cooling of deep ocean waters associated with the establishment of the Antarctic ice sheet at the Eocene-Oligocene climate transition. We then apply this approach to new records of the Middle Miocene Climate Transition from ODP Site 761 to estimate variations in bottom water temperature and the oxygen isotopic composition of seawater. We estimate that the oxygen isotopic composition of seawater varied by ~1 per mil between the deglacial extreme of the Miocene Climatic Optimum and the glacial maximum following the Middle Miocene Climate Transition, indicating large amplitude variations in ice volume. However, the longer-term change between 15.3 and 12.5 Ma is marked by a ~1°C cooling of deep waters, and an increase in the oxygen isotopic composition of seawater of ~0.6 per mil. We find that bottom water saturation state increased in the lead up to the Middle Miocene Climate Transition and decreased shortly after. This supports decreasing pCO2 as a driver for global cooling and ice sheet expansion, in agreement with existing boron isotope and leaf stomatal index CO2 records but in contrast to the published alkenone CO2 records.

Siliceous components of Pliocene-Holocene sediments of the South Atlantic

The Pliocene-Holocene sediments recovered on ODP Leg 114 from Holes 699A, 701C, and 704B are the subject of a detailed investigation to interpret changes in the Oceanographic environment of the South Atlantic in the vicinity of the Polar Front Zone (PFZ). The cores sample sediments at shallow (Hole 704B, 2532 m), intermediate (Hole 699A, 3716 m), and basinal (Hole 701C, 4647 m) depths. Sites 699 and 704 come under the influence of the Antarctic Circumpolar Current (ACC) and Circumpolar Deep Water. It is possible that the upper reaches of Antarctic Bottom Water (AABW) may also affect Hole 699A. Site 701 is influenced by AABW. Closely spaced samples were analyzed for grain-size distribution, sand fraction components, biosiliceous microfossils, organic carbon, and water content. PFZ migrations are traced using changes in bulk sedimentaccumulation rates and the abundance of the diatoms Actiniscus ssp. and Genus et species indet. 1 Fenner (1991), as well as changes in sediment grain size and composition. Diatomaceous sediments of Gilbert age in Hole 699A indicate that the PFZ was positioned over this site, but during the Gauss it migrated north, bringing in less productive Antarctic Surface Water. All cores document a very gradual southerly movement of the PFZ throughout the Matuyama (with some sharp fluctuations of the northen PFZ border over Site 704 between 1.45 and 1.83 m.y.). This regressive shift culminated in the late Matuyama. The latest Matuyama to earliest Brunhes record in Hole 699A has been removed by a hiatus lasting from 1.0 to 0.6 m.y., which was probably caused by intensification of the deep-reaching ACC. The corresponding interval in Hole 704B, the shallowest core, contains evidence of winnowing. Sharp fluctuations of large amplitude and high frequency in the lithology of the sediments from Hole 704B in the eastern South Atlantic, starting at about 0.75 m.y. and characterizing the whole Brunhes Epoch, record the rapid movement of the northern border of the PFZ over the site. These reflect strong glacial/interglacial alternations in climate. To a lesser extent, lithologic fluctuations in Hole 701C reflect the same phenomenon, whereas in Hole 699A the lithology does not vary as dramatically.

Age model and grain size distribution of the fine fraction of ODP Site 114-704

The physical properties of sediments beneath an upwelling area in the southern part of the Atlantic Ocean (ODP Hole 704A) were investigated. Highly significant correlations characterize the relationship of carbonate content to bulk density (R = 0.85), carbonate content to porosity (R = 0.84), and carbonate content to impedance (R = 0.84). No relationship exists between carbonate content and compressional-wave velocity (R = 0.24), indicating that amplitude variations in impedance are primarily controlled by variations in bulk density, which, in turn, are controlled by climatically driven biogenic opal and carbonate deposition. In general, maxima in impedance correspond to maxima in carbonate content (minima in opal content). The impedance record exhibits its most drastic change at about 2.4 Ma, marking dramatic increases in the average content of biogenic opal and the beginning of large-amplitude fluctuations. Between 0.7 and 0.4 Ma carbonate content, bulk density, and grain density decrease while opal content drastically increases. Similar changes have been observed in sediments beneath an upwelling cell off northwest Africa, indicating an oceanwide enhancement in upwelling or in the calcite corrosiveness of bottom water that appears to be isochronous.

Revision of the early-middle Pleistocene calcareous nannofossil biochronology

The extant nannofossil biostratigraphic and biochronologic framework for the early-middle Pleistocene time interval has been tested through the micropaleontological analysis of globally distributed high-quality low- to mid-latitude deep-sea successions. The quantitative temporal distribution patterns of relative abundances of selected taxa were reconstructed in critical intervals, and the following biohorizons were defined: first occurrence of medium-sized Gephyrocapsa spp. (bmG); last occurrence of Calcidiscus macintyrei (tCm); first occurrence of large Gephyrocapsa spp. (blG); last occurrence of large Gephyrocapsa spp. (tlG); first occurrence of Reticulofenestra asanoi (bRa); re-entrance of medium-sized Gephyrocapsa spp. (reemG) and last occurrence of Reticulofenestra asanoi (tRa). The detailed patterns of abundance change at these biohorizons were used to generate a detailed biostratigraphy, and the biostratigraphic data were transformed into a precise biochronology by means of correlation to isotope stratigraphies and astronomical timescales. The degree of isochrony or diachrony of the biohorizons was evaluated. Biohorizons tlG and tRa are isochronous occurring close to marine isotope stages (MIS)55 and MIS 22, respectively, and bmG and blG are slightly diachronous on the order of 30-40 kyr, whereas biohorizons tCm, reemG and bRa are confirmed as diachronous on the order of 100, 80 and 60 kyr, respectively. Some of the events are clearly controlled by environmental conditions, e.g. the last occurrence of R. asanoi, related to significant environmental changes associated with the first large-amplitude glaciation of the late Quaternary, MIS 22.

An integrated biostratigraphy (conodonts and foraminifers) and chronostratigraphy (paleomagnetic reversals, magnetic susceptibility, elemental chemistry, carbon isotopes and geochronology) for the Permian-Triassic strata of Guandao section

The chronostratigraphy of Guandao section has served as the foundation for numerous studies of the end-Permian extinction and biotic recovery in south China. Guandao section is continuous from the Permian-Triassic boundary to the Upper Triassic.Conodonts enable broad delineation of stage and substage boundaries and calibration of foraminifer biostratigraphy as follows. Changhsingian- Griesbachian: first Hindeodus parvus, and first appearance of foraminifers Postcladella kalhori and Earlandia sp. Griesbachian-Dienerian: first Neospathodus dieneri, and last appearance of foraminifer P. grandis. Dienerian-Smithian: first Novispathodus waageni and late Dienerian first appearance of foraminifer Hoyenella ex gr. sinensis. Smithian-Spathian: first Nv? crassatus and last appearance of foraminifers Arenovidalina n. sp. and Glomospirella cf. vulgaris. Spathian-Aegean: first Chiosella timorensis and first appearance of foraminifer Meandrospira dinarica. Aegean-Bithynian: first Nicoraella germanica and first appearance of foraminifer Pilammina densa. Bithynian-Pelsonian: after last Neogondolella regalis, prior to first Paragondolella bulgarica and first appearance of foraminifer Aulotortus eotriasicus. Pelsonian-Illyrian: first Pg. excelsa and last appearance of foraminifers Meandrospira ? deformata and Pilamminella grandis. Illyrian-Fassanian: first Budurovignathus truempyi, and first appearance of foraminifers Abriolina mediterranea and Paleolituonella meridionalis. Fassanian-Longobardian: first Bv. mungoensis and last appearance of foraminifer A. mediterranea. Longobardian-Cordevolian: first Quadralella polygnathiformis and last appearance of foraminifers Turriglomina mesotriasica and Endotriadella wirzi. The section contains primary magnetic signature with frequent reversals occurring around the Permian-Triassic, Olenekian-Anisian, and Anisian-Ladinian boundaries. Predominantly normal polarity occurs in the lower Smithian, Bithynian, and Longobardian-Cordevolian. Predominantly reversed polarity occurs in the upper Griesbachian, Induan-Olenekian, Pelsonian and lower Illyrian. Reversals match well with the GPTS. Large amplitude carbon isotope excursions, attaining values as low as -2.9 per mil d13C and high as +5.7 per mil d13C, characterize the Lower Triassic and basal Anisian. Values stabilize around +2 per mil d13C through the Anisian to Carnian. Similar signatures have been reported globally. Magnetic susceptibility and synthetic gamma ray logs show large fluctuations in the Lower Triassic and an overall decline in magnitude of fluctuation through the Middle and Upper Triassic. The largest spikes in magnetic susceptibility and gamma ray, indicating greater terrestrial lithogenic flux, correspond to positive d13C excursions. Several volcanic ash horizons occur in the Lower Triassic and Olenekian-Anisian boundary. High resolution U-Pb analysis of zircons provide a robust age of 247.2 Ma for the Olenekian-Anisian boundary.

Carbon and oxygen isotopic composition of benthic foraminifers from ODP Site 167-1014 and Hole 167-1020C

Benthic foraminiferal stable isotope data are presented for Sites 1014 (Tanner Basin, 1176 m) and 1020 (Gorda Ridge, 3040 m) to constrain past changes in Pacific deep- and intermediate-water nutrient chemistry associated with the onset of large-amplitude 100-k.y. climate cycles after ~900 ka. The Site 1014 data were based on analyses of separate species of Cibicidoides, whereas only Cibicidoides wuellerstorfi was used to generate the Site 1020 record. The present data span 380-920 and 620-950 ka at Sites 1014 and 1020, respectively.

Geochemistry of tephra from the Taupo Volcanic Zone

The Taupo Volcanic Zone (TVZ), central North Island, New Zealand, is the most frequently active Quaternary rhyolitic system in the world. Silicic tephras recovered from Ocean Drilling Programme Site 1123 (41°47.16'S, 171°29.94'W; 3290 m water depth) in the southwest Pacific Ocean provide a well-dated record of explosive TVZ volcanism since ~1.65 Ma. We present major, minor and trace element data for 70 Quaternary tephra layers from Site 1123 determined by electron probe microanalysis (1314 analyses) and laser ablation inductively coupled plasma mass spectrometry (654 analyses). Trace element data allow for the discrimination of different tephras with similar major element chemistries and the establishment of isochronous tie-lines between three sediment cores (1123A, 1123B and 1123C) recovered from Site 1123. These tephra tie-lines are used to evaluate the stratigraphy and orbitally tuned stable isotope age model of the Site 1123 composite record. Trace element fingerprinting of tephras identifies ~4.5 m and ~7.9 m thick sections of repeated sediments in 1123A (49.0-53.5 mbsf [metres below seafloor]) and 1123C (48.1-56.0 mbsf), respectively. These previously unrecognised repeated sections have resulted in significant errors in the Site 1123 composite stratigraphy and age model for the interval 1.15-1.38 Ma and can explain the poor correspondence between d18O profiles for Site 1123 and Site 849 (equatorial Pacific) during this interval. The revised composite stratigraphy for Site 1123 shows that the 70 tephra layers, when correlated between cores, correspond to ~37-38 individual eruptive events (tephras), 7 of which can be correlated to onshore TVZ deposits. The frequency of large-volume TVZ-derived silicic eruptions, as recorded by the deposition of tephras at Site 1123, has not been uniform through time. Rather it has been typified by short periods (25-50 ka) of intense activity bracketed by longer periods (100-130 ka) of quiescence. The most active period (at least 1 event per 7 ka) occurred between ~1.53 and 1.66 Ma, corresponding to the first ~130 ka of TVZ rhyolitic magmatism. Since 1.2 Ma, ~80% of tephras preserved at Site 1123 and the more proximal Site 1124 were erupted and deposited during glacial periods. This feature may reflect either enhanced atmospheric transport of volcanic ash to these sites (up to 1000 km from source) during glacial conditions or, more speculatively, that these events are triggered by changes in crustal stress accumulation associated with large amplitude sea-level changes. Only 8 of the ~37-38 Site 1123 tephra units (~20%) can be found in all three cores, and 22 tephra units (~60%) are only present in one of the three cores. Whether a tephra is preserved in all three cores does not have any direct relationship to eruptive volume. Instead it is postulated that tephra preservation at Site 1123 is 'patchy' and influenced by the vigorous nature of their deposition to the deep ocean floor as vertical density currents. At this site, at least 5 cores would need to have been drilled within a proximity of 10's to 100's of metres of each other to yield a >99% chance of recovering all the silicic tephras deposited on the ocean surface above it in the past 1.65 Ma.

Stable isotope record and age model for sapropel stratigraphy of ODP Site 160-967

Stable oxygen isotope data from four holes drilled at the Ocean Drilling Program Site 967, which is located on the lower northern slope of the Eratosthenes Seamount, provide a continuous record of Eastern Mediterranean surface-water conditions during the last 3.2 Ma. A high-resolution stratigraphy for the Pliocene-Pleistocene sequence was established by using a combination of astronomical calibration of sedimentary cycles, nannofossil stratigraphy, and stable oxygen isotope fluctuations. Sapropels and color cycles are present throughout the last 3.2 Ma at Site 967, and their ages, as determined by calibration against the precessional component of the astronomical record, are consistent with those estimated for the sapropels of the classical land-based marine sequences of the Punta Piccola, San Nicola, Singa, and Vrica sections (southern Italy). The Site 967 oxygen isotope record shows large amplitude fluctuations mainly caused by variations in surface water salinity throughout the entire period. Spectral analysis shows that fluctuations in the d18O record were predominantly influenced by orbital obliquity and precessional forcing from 3.2 to 1 Ma, and all main orbital frequencies characterize the d18O record for the last million years. The start of sapropel formation at 3.2 Ma indicates a possible link between sapropel formation and the build up of northern hemisphere ice sheets. The dominance of the obliquity cycle in the interval from 3.2-1 Ma further points to the sensitivity of Eastern Mediterranean climate to the fluctuations in the volume of Arctic ice sheets. An intensification of negative isotope anomalies at Site 967, relative to the open ocean, supports a link between high run-off (during warm periods) and sapropel formation. freshwater input would have inhibited deep-water formation, which led to stagnation of deeper waters. Comparison with the land sections also confirms that differential preservation and diagenesis play a key role in sapropel occurrence.

Pliocene to Pleistocene calcium carbonate and organic carbon record of ODP Hole 117-722B

Site 722 provides high resolution records of percent CaCO3, magnetic susceptibility, d18O, organic carbon, and coarse fraction for the past 3.4 m.y. from the crest of the Owen Ridge, northwestern Arabian Sea. Within this time interval, most of the carbonate percent variations can be attributed to terrigenous dilution and do not reflect changes in the carbonate system. From the late Pliocene to Present, the average rate of calcium carbonate accumulation increases from 1 to 3 g/cm**2/k.y. and the average accumulation of organic carbon decreases from 75 to 30 mg/cm**2/k.y. The carbonate component is more dissolved in the older interval. The long-term variations in carbonate accumulation may reflect a greater input of organic matter in the late Pliocene, which decomposes to produce CO2 and dissolve carbonate. Magnetic susceptibility and % noncarbonate (100 - CaCO3%) reflect changes in the amount of the lithogenic component in the sediments. The period of variation of lithogenic material is the same period as the original forcing of the regional summer monsoon, however, the timing matches global aridity patterns and global ice volume (sea level) changes. This preliminary analysis suggests that the high frequency variation of lithogenic material persists for at least the last 3.4 m.y. Within the last million years, calcium carbonate accumulation has a large amplitude signal that covaries with major changes in ice volume. Both calcium carbonate and noncarbonate (mostly terrigenous) accumulation are greatest during glacial stages. Interglacial intervals are characterized by low mass accumulation rates, increased foraminifer fragmentation, and increased opal concentration. The accumulation of organic carbon matches the high frequency changes in sedimentation rates. We attribute this high correlation to enhanced preservation of organic carbon by increased sedimentation rate. Of the three major biological components studied, only opal exhibits the variations expected for a biological productivity system forced by monsoonal upwelling driven by changes in northern hemisphere summer radiation.