Stable isotope analysis and temperature reconstruction data from DSDP Hole 94-609B and ODP Hole 162-984B

The Greenland ice sheet is accepted as a key factor controlling the Quaternary glacial scenario. However, the origin and mechanisms of major Arctic glaciation starting at 3.15 Ma and culminating at 2.74 Ma are still controversial. For this phase of intense cooling Ravelo et al. proposed a complex gradual forcing mechanism. In contrast, our new submillennial-scale paleoceanographic records from the Pliocene North Atlantic suggest a far more precise timing and forcing for the initiation of northern hemisphere glaciation (NHG), since it was linked to a 2-3 °C surface water warming during warm stages from 2.95 to 2.82 Ma. These records support previous models, claiming that the final closure of the Panama Isthmus (3.0- ~2.5 Ma induced an increased poleward salt and heat transport. Associated strengthening of North Atlantic Thermohaline Circulation and in turn, an intensified moisture supply to northern high latitudes resulted in the build-up of NHG, finally culminating in the great, irreversible climate crash at marine isotope stage G6 (2.74 Ma). In summary, there was a two-step threshold mechanism that marked the onset of NHG with glacial-to-interglacial cycles quasi-persistent until today.

Neogene calcareous nannofossil abundance and datums of ODP Site 186-1150 and 186-1151 sediments

Calcareous nannofossil assemblages were studied from Ocean Drilling Program Holes 1150A, 1150B, 1151A, 1151C, and 1151D in order to estimate the age of sediments drilled in the Japan Trench of the western Pacific Ocean. The abundance and species diversity of nannofossil flora are generally low but are sufficient to show that the sedimentary sequences range from Quaternary to Miocene in age (nannofossil Zones CN15-CN3). The abundance of Coccolithus pelagicus, a cold-water indicator, was studied from sediments younger than 3.83 Ma from both Holes 1150A and 1151A in order to elucidate past climate conditions. Between 3.83 and 2.82 Ma, the abundance of C. pelagicus was generally low, but abundance increased significantly after 2.82 Ma. In agreement with previous studies, this increase appears to be related to a change in the current system around the western Pacific Ocean and eastern Atlantic Ocean that occurred in response to the final elevation of the Isthmus of Panama.

Model-based changes in global annual mean surface temperature change (Delta T_g) and radiative forcing due to land ice albedo changes (Delta R_[LI]) over the last 5 Myr, supplementary material

It is still an open question how equilibrium warming in response to increasing radiative forcing - the specific equilibrium climate sensitivity S - depends on background climate. We here present palaeodata-based evidence on the state dependency of S, by using CO2 proxy data together with a 3-D ice-sheet-model-based reconstruction of land ice albedo over the last 5 million years (Myr). We find that the land ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcing depends on the CO2 data set used. This non-linearity has not, so far, been accounted for in similar approaches due to previously more simplistic approximations, in which land ice albedo radiative forcing was a linear function of sea level change. The latitudinal dependency of ice-sheet area changes is important for the non-linearity between land ice albedo and sea level. In our set-up, in which the radiative forcing of CO2 and of the land ice albedo (LI) is combined, we find a state dependence in the calculated specific equilibrium climate sensitivity, S[CO2,LI], for most of the Pleistocene (last 2.1 Myr). During Pleistocene intermediate glaciated climates and interglacial periods, S[CO2,LI] is on average ~ 45 % larger than during Pleistocene full glacial conditions. In the Pliocene part of our analysis (2.6-5 Myr BP) the CO2 data uncertainties prevent a well-supported calculation for S[CO2,LI], but our analysis suggests that during times without a large land ice area in the Northern Hemisphere (e.g. before 2.82 Myr BP), the specific equilibrium climate sensitivity, S[CO2,LI], was smaller than during interglacials of the Pleistocene. We thus find support for a previously proposed state change in the climate system with the widespread appearance of northern hemispheric ice sheets. This study points for the first time to a so far overlooked non-linearity in the land ice albedo radiative forcing, which is important for similar palaeodata-based approaches to calculate climate sensitivity. However, the implications of this study for a suggested warming under CO2 doubling are not yet entirely clear since the details of necessary corrections for other slow feedbacks are not fully known and the uncertainties that exist in the ice-sheet simulations and global temperature reconstructions are large.

Stable isotope analysis of benthic foraminifera and sedimentology on ODP Site 162-984

Sediments recovered at ODP Site 984 on the Reykjanes Ridge provided multicentennial-scale records (SST, planktic and benthic delta18O, IRD and magnetic susceptibility) of Late Pliocene climate change over the onset of Northern Hemisphere glaciation (NHG), 2.95-2.82 Ma. Short-term climate variations prior and after the onset of continent-wide glaciation were compared to test the hypothesis, whether Dansgaard-Oeschger (DO) cycles may have been triggered by continental ice breakouts. By means of spectral analyses for two selected interglacial stages prior to and after NHG (G15 and G1), we found that climate variability resembled that of the Holocene and the mid-Pliocene warm period. In contrast, DO-like periodicities near 1470, 2900, and 4400 yr indeed only occurred in glacial stages after the onset of NHG (G14, G6, and 104), but hardly in stage G20 prior to the onset. These results suggest a causal link between DO cycles and the Late Pliocene onset of major NHG and ice breakouts in the North Atlantic. This data set provides all primary data and spectral analysis related to this scientific work.

Isotopic ratios, concentration of noble gases and mineral and bulk composition of extraterrestrial dust particles in Dome Fuji ice core, East Antarctica

Two silicate-rich dust layers were found in the Dome Fuji ice core in East Antarctica, at Marine Isotope Stages 12 and 13. Morphologies, textures, and chemical compositions of constituent particles reveal that they are high-temperature melting products and are of extraterrestrial origin. Because similar layers were found ~2000 km east of Dome Fuji, at EPICA (European Project for Ice Coring in Antarctica)-Dome C, particles must have rained down over a wide area 434 and 481 ka. The strewn fields occurred over an area of at least 3 × 10**6 km**2. Chemical compositions of constituent phases and oxygen isotopic composition of olivines suggest that the upper dust layer was produced by a high-temperature interaction between silicate-rich melt and water vapor due to an impact explosion or an aerial burst of a chondritic meteoroid on the inland East Antarctic ice sheet. An estimated total mass of the impactor, on the basis of particle flux and distribution area, is at least 3 × 10**9 kg. A possible parent material of the lower dust layer is a fragment of friable primitive asteroid or comet. A hypervelocity impact of asteroidal/cometary material on the upper atmosphere and an explosion might have produced aggregates of sub-µm to µm-sized spherules. Total mass of the parent material of the lower layer must exceed 1 × 10**9 kg. The two extraterrestrial horizons, each a few millimeters in thickness, represent regional or global meteoritic events not identified previously in the Southern Hemisphere.

Composition of Eocene claystones of ODP Hole 105-647A

Authigenic carbonates were recovered in lower to middle Eocene claystones at Ocean Drilling Program Site 647 in the Labrador Sea. Detailed chemical, petrographic, and X-ray investigations reveal that these diagenetic carbonates have a complex mineralogical composition. At least five different carbonate phases are identified: calcium-rich rhodochrosite, rhodochrosite, manganosiderite, siderite, and calcite. Manganese carbonates are the dominant carbonate phases formed throughout the section. Textural analyses show two major generations of carbonate formation. Early cementation of micritic carbonate in burrow structures was followed by carbonate cementation forming microsparry to sparry crystals. At approximately 620 meters below seafloor (mbsf), three concretions of iron carbonates occur, which indicates a special pore-water chemistry. Thin section analyses from this level show (1) several generations of diagenetic carbonates, (2) widespread secondary cavity formation in burrow structures, and (3) various cement precipitations in voids. We suggest that this level represents a hiatus or highly condensed sequence, as indicated by (1) the low carbonate content in host sediments, (2) carbonate dissolution reflected by the high ratio of benthic to planktonic foraminifers, and (3) complex diagenetic alteration in the carbonate concretions. Iron and manganese enrichments observed in lithologic Unit IV may have been derived from a hydrothermal source at the adjacent, then active, Labrador Sea mid-ocean ridge. Authigenic smectites forming numerous pseudomorphs of siliceous microfossils are precipitated in burrow structures. We propose that diagenetic smectite formation from biogenic opal and iron oxyhydroxide (analogous to smectite formation in surface sediments of the East Pacific area) occurred in the Labrador Sea during the early and middle Eocene.

Effective stress, porosity, p-wave velocity and mineral composition of ODP Hole 174A-1073A sediments

Porosity, permeability, and compressional (P-wave) velocity were measured as a function of stress on sediments from Ocean Drilling Program Site 1073, U.S. Mid-Atlantic continental slope. Thin sections, scanning electron microscopy, and X-ray diffraction analyses provided mineralogical characteristics of the samples. Uniaxial strain boundary conditions were imposed on the samples during consolidation tests with the maximum effective axial stress reaching 13 MPa. The maximum effective radial stress necessary to maintain uniaxial strain was 7.6 MPa. Over an effective axial stress interval of 0 to 5.2 MPa, Sample 174A-1073A-26X-2, 82-89 cm (226.65 meters below seafloor [mbsf]), exhibited the largest decrease in porosity (51% to 41%), whereas Sample 71X-1, 2-8 cm (644.70 mbsf), exhibited the smallest decrease in porosity (48% to 45%). All samples showed negligible porosity increases during unloading. The permeability (on the order of 1 x 10-17 m**2) of Sample 174A-1073A-71X-1, 2-8 cm, was twice that measured on Sample 8H-1, 23-26 cm (63.75 mbsf), even though the former was considerably deeper and older. The differences in porosity-stress behavior and permeability between shallow and deep samples is related to lithologic, mineralogic, and diagenetic differences between the sediments above and below the Pliocene-Pleistocene to Miocene unconformity. P-wave velocity for Samples 174A-1073A-41X-5, 97-103 cm (372.35 mbsf), and 71X-1, 2-8 cm, increased with decreasing porosity, but did not change significantly during unloading.