Cenozoic oxygen isotopic values for benthic foraminifera from DSDP and ODP low and high latitude marine sediment cores

The climate during the Cenozoic era changed in several steps from ice-free poles and warm conditions to ice-covered poles and cold conditions. Since the 1950s, a body of information on ice volume and temperature changes has been built up predominantly on the basis of measurements of the oxygen isotopic composition of shells of benthic foraminifera collected from marine sediment cores. The statistical methodology of time series analysis has also evolved, allowing more information to be extracted from these records. Here we provide a comprehensive view of Cenozoic climate evolution by means of a coherent and systematic application of time series analytical tools to each record from a compilation spanning the interval from 4 to 61 Myr ago. We quantitatively describe several prominent features of the oxygen isotope record, taking into account the various sources of uncertainty (including measurement, proxy noise, and dating errors). The estimated transition times and amplitudes allow us to assess causal climatological-tectonic influences on the following known features of the Cenozoic oxygen isotopic record: Paleocene-Eocene Thermal Maximum, Eocene-Oligocene Transition, Oligocene-Miocene Boundary, and the Middle Miocene Climate Optimum. We further describe and causally interpret the following features: Paleocene-Eocene warming trend, the two-step, long-term Eocene cooling, and the changes within the most recent interval (Miocene-Pliocene). We review the scope and methods of constructing Cenozoic stacks of benthic oxygen isotope records and present two new latitudinal stacks, which capture besides global ice volume also bottom water temperatures at low (less than 30°) and high latitudes. This review concludes with an identification of future directions for data collection, statistical method development, and climate modeling.

Geochemistry of ODP Site 145-882

Measurements of benthic foraminiferal cadmium:calcium (Cd/Ca) have indicated that the glacial-interglacial change in deep North Pacific phosphate (PO4) concentration was minimal, which has been taken by some workers as a sign that the biological pump did not store more carbon in the deep glacial ocean. Here we present sedimentary redox-sensitive trace metal records from Ocean Drilling Program (ODP) Site 882 (NW subarctic Pacific, water depth 3244 m) to make inferences about changes in deep North Pacific oxygenation - and thus respired carbon storage - over the past 150,000 yr. These observations are complemented with biogenic barium and opal measurements as indicators for past organic carbon export to separate the influences of deep-water oxygen concentration and sedimentary organic carbon respiration on the redox state of the sediment. Our results suggest that the deep subarctic Pacific water mass was depleted in oxygen during glacial maxima, though it was not anoxic. We reconcile our results with the existing benthic foraminiferal Cd/Ca by invoking a decrease in the fraction of the deep ocean nutrient inventory that was preformed, rather than remineralized. This change would have corresponded to an increase in the deep Pacific storage of respired carbon, which would have lowered atmospheric carbon dioxide (CO2) by sequestering CO2 away from the atmosphere and by increasing ocean alkalinity through a transient dissolution event in the deep sea. The magnitude of change in preformed nutrients suggested by the North Pacific data would have accounted for a majority of the observed decrease in glacial atmospheric pCO2.

Grain-size distributions and bulk chemistry of sediment core MD00-2361

The terrigenous fraction of sediments from a deep-sea sediment core recovered from the northwestern Western Australian continental slope offshore North West Cape, SE Indian Ocean, reveals a history of Western Australian climate throughout the last 550 ka. End-member modelling of a data set of grain-size distributions (n = 438) of the terrigenous sediment fraction allows to interpret the record in terms of aeolian and fluvial sediment deposition, both related to palaeo-environmental conditions in the North West Cape area. The data set can be best described by two aeolian end members and one fluvial one. Changes in the ratio of the two aeolian end members over the fluvial one are interpreted as aridity variations in northwestern Western Australia. These grain-size data are compared with bulk geochemical data obtained by XRF scans of the core. Log-ratios of the elements Zr/Fe and Ti/Ca, which indicate a terrigenous origin, corroborate the grain-size data. We postulate that the mid- to late Quaternary near North West Cape climate was relatively arid during the glacial and relatively humid during the interglacial stages, owing to meridional shifts in the atmospheric circulation system. Opposite to published palaeo-environmental records from the same latitude (20°S) offshore Chile and offshore Namibia, the Australian aridity record does not show the typical southern hemisphere climate variability of humid glacials and dry interglacials, which we interpret to be the result of the relatively large land mass of the Australian continent, which emphasises a strong monsoonal climatic system.

Trioctahedral micas in xenolithic ejecta from recent volcanism of the Somma-Vesuvius

This study reports the first crystal chemical database resulting from a detailed structural investigation of trioctahedral micas found in xenolithic ejecta produced during the AD 1631, 1872 and 1944 eruptions, three explosive episodes of recent volcanic period of Vesuvius volcano (Southern Italy). Three xenolith types were selected: metamorphic/metasomatic skarns, pyrometamorphic/hydrothermally altered nodules and mafic cumulates. They are related to different magma chemistry and effusive styles: from sub-plinian and most evolved (AD 1631 eruption) to violent strombolian with medium evolution degree (AD 1872 eruption) to vulcanian-effusive, least evolved (AD 1944 eruption) event, respectively. Both xenoliths and micas were investigated employing multiple techniques: the xenoliths were characterized by X-ray fluorescence, inductively-coupled plasma-mass spectrometry, optical microscopy, X-ray powder diffraction, and quantitative energy-dispersive microanalysis; the micas were studied by electron probe microanalysis and single crystal X-ray diffraction. The mica-bearing xenoliths show variable texture and mineralogical assemblage, clearly related to their different origin. Based on the major oxide chemistry, only one xenolithic sample falls in the skarn compositional field from the Somma-Vesuvius literature, some fall close to the skarns and cumulate fields, others plot close to the syenite/foidolite/essexite field. A subgroup of the selected ejecta does not fall or approach any of the compositional fields. Trace and rare earth element patterns show some petrological affinity between studied xenoliths and erupted magmas with typical Eu, Ta and Nb negative anomalies. Strongly depleted patterns were detected for the 1631 metamorphic/metasomatic skarns xenoliths. Three distinct mica groups were distinguished: 1) Mg-, Al-rich, low Ti-bearing, low to moderate F-bearing varieties (1631 xenolith), 2) Al-moderate, F- and Mg-rich, Ti-, Fe-poor varieties (1872 xenolith), and 3) Al-, Ti- and Fe-rich, F-poor phases (1944 xenolith). All the analysed mica crystals are 1M polytypes with the expected space group C2/m. Micas from xenoliths of the 1631 Vesuvius eruption are phlogopites characterized by a combination of low extent of oxy-type and variable extent OH-F-substitutions, as testified by the range of F concentration (from ~ 0.20 to 0.80 apfu). Micas from xenoliths of the 1872 Vesuvius eruption exhibit structural peculiarities typical of fluorophlogopites, i.e. OH-F-substitution is predominant. Micas from the xenolith of the 1944 Vesuvius eruption display features typical of oxy-substituted micas. The variability of the crystal chemical features of the studied micas are consistent with the remarkable variation of their host rocks. Micas from 1631 nodules are related to metasomatic, skarn-type environment, deriving from the metamorphosed wall-rocks hosting the magma reservoir. The fluorophlogopites from the 1872 xenoliths testify for strongly dehydrated environmental conditions compared to those of the 1631 and 1944 hosts. Finally, magma storage condition at depth, associated to a decreasing aH2O may have promoted major oxy-type substitutions in 1944 biotites.

Temperature observation of CORK installation during ODP Leg 168

Four boreholes, drilled a few tens of meters into igneous basement on the eastern flank of the Juan de Fuca Ridge during ODP Leg 168, were sealed and instrumented for long-term monitoring to observe the hydrologic state of young sediment-sealed oceanic crust. The thermal regime is dominated by the effects of rapid fluid circulation in uppermost igneous basement driven by very small non-hydrostatic pressure gradients. Upper basement temperatures are uniform laterally between pairs of holes over distances of hundreds of meters to kilometers. In the case of two holes drilled into a sediment-buried basement ridge and adjacent valley, basement temperatures differ by less than 2 K despite the 2.2 km lateral separation of the sites and the 2.5:1 contrast in sediment cover thickness. Under conductive conditions, upper basement temperatures would differ by roughly 50 K. By comparison with modeling results, the observed degree of isothermality suggests a fluid flux of at least 10**-6 m/s (30 m/yr), and an effective permeability in the range of 10**-10 to 10**-9 m**2 in the uppermost igneous crust. The pressure difference available to drive this rapid flux between the ridge and valley, estimated by comparing the observed pressures via the isothermal upper basement hydrostat that is inferred to connect the two sites, is small (~2 kPa) and also suggests high permeability. Relative to the hydrostats defined by the local conductive sediment geotherms, substantial super-hydrostatic pressure (+18 kPa) is present within the buried basement ridge, and sub-hydrostatic pressure is present in the adjacent valley (-26 kPa). Such pressure differentials are the direct consequence of the advection-dominated thermal regime and small pressure losses in high-permeability basement, and are available to drive fluid seepage through sediment sections vertically up above and horizontally away from buried ridges, and down above valleys. No constraints are provided by any of the observations on the depth in the crust to which thermally or chemically significant flow might extend, although just as in the overlying sediments, the pattern of deep flow may be affected by the near-isothermal and near-hydrostatic conditions present in the permeable uppermost crustal section.

Geochemistry of Hawaiian volcanics

The "Ko'olau" component of the Hawaiian mantle plume represents an extreme (EM1-type) end member of Hawaiian shield lavas in radiogenic isotope space, and was defined on the basis of the composition of subaerial lavas exposed in the Makapu'u section of Ko'olau Volcano. The 679 m-deep Ko'olau Scientific Drilling Project (KSDP) allows the long-term evolution of Ko'olau Volcano to be reconstructed and the longevity of the "Ko'olau" component in the Hawaiian plume to be tested. Here, we report triple spike Pb isotope and Sr and Nd isotope data on KSDP core samples, and rejuvenation stage Honolulu Volcanics (HV) (together spanning ~2.8 m.y.), and from ~110 Ma basalts from ODP Site 843, thought to be representative of the Pacific lithosphere under Hawai'i. Despite overlapping ranges in Pb isotope ratios, KSDP and HV lavas form two distinct linear arrays in 208Pb/204Pb-206Pb/204Pb isotope space. These arrays intersect at the radiogenic end indicating they share a common component. This "Kalihi" component has more radiogenic Pb, Nd, Hf, but less radiogenic Sr isotope ratios than the "Makapu'u" component. The mixing proportions of these two components in the lavas oscillated through time with a net increase in the "Makapu'u" component upsection. Thus, the "Makapu'u" enriched component is a long-lived feature of the Hawaiian plume, since it is present in the main shield-building stage KSDP lavas. We interpret the changes in mixing proportions of the Makapu'u and Kalihi components as related to changes in both the extent of melting as well as the lithology (eclogite vs. peridotite) of the material melting as the volcano moves away from the plume center. The long-term Nd isotope trend and short-term Pb isotope fluctuations seen in the KSDP record cannot be ascribed to a radial zonation of the Hawaiian plume: rather, they reflect the short length-scale heterogeneities in the Hawaiian mantle plume. Linear Pb isotope regressions through the HV, recent East Pacific Rise MORB and ODP Site 843 datasets are clearly distinct, implying that no simple genetic relationship exists between the HV and the Pacific lithosphere. This observation provides strong evidence against generation of HV as melts derived from the Pacific lithosphere, whether this be recent or old (100 Ma). The depleted component present in the HV is unlike any MORB-type mantle and most likely represents material thermally entrained by the upwelling Hawaiian plume and sampled only during the rejuvenated stage. The "Kalihi" component is predominant in the main shield building stage lavas but is also present in the rejuvenated HV. Thus this material is sampled throughout the evolution of the volcano as it moves from the center (main shield-building stage) to the periphery (rejuvenated stage) of the plume. The presence of a plume-derived material in the rejuvenated stage has significant implications for Hawaiian mantle plume melting models.

Analysis of Mn-nodules from Breitgrund, Baltic Sea

The usually high concentrations of Zn, Pb, Cd, and Cu in the most recently accreted portions of ferromanganese nodules from the western Baltic Sea are thought to reflect increased metal input due to anthropogenic mobilization. If so, the point of increase represents a time horizon within the structure of the nodule. Similar trace metal distributions of radiometrically dated sediments from the same area suggest that the ferromanganese nodules have grown in thickness between 0.02 and 0.16 mm yr-1. From this growth rate anthropogenic Zn flux to the nodule surface was calculated to be 80 mg m-2 yr-1.

Fossil chironomids of Lake Temje

A 380 cm long sediment core from Lake Temje (central Yakutia, Eastern Siberia) was studied to infer Holocene palaeoenvironmental change in the extreme periglacial setting of eastern Siberia during the last 10,000 years. Data on sediment composition were used to characterize changes in the depositional environment during the ontogenetic development of the Lake Temje. The analysis of fossil chironomid remains and statistical treatment of chironomid data by the application of a newly developed regional Russian transfer functions provided inferences of mean July air temperatures (T_July) and water depths (WD). Reconstructed WDs show minor changes throughout the core and range between 80 and 120 cm. All the fluctuations in reconstructed water depth lie within the mean error of prediction of the inference model (RMSEP = 0.35) so it is not possible to draw conclusions from the reconstructions. A qualitative and quantitative reconstruction of Holocene climate in central Yakutia recognized three stages of palaeoenvironmental changes. The early Holocene between 10 and 8 ka BP was characterized by colder-than-today and moist summer conditions. Cryotextures in the lake sediments document full freezing of the lake water during the winter time. A general warming trend started around 8.0 ka BP in concert with enhanced biological productivity. Reconstructed mean T_July were equal or up to 1.5 °C higher than today between 6.0 ka and 5.0 ka BP. During the entire late Holocene after 4.8 ka BP, reconstructed mean T_July remained below modern value. Limnological conditions did not change significantly. The inference of a mid-Holocene climate optimum supports scenarios of Holocene climatic changes in the subpolar part of eastern Siberia and indicates climate teleconnections to the North Atlantic realm.

Late Cretaceous planktonic foraminiferal assemblages from the Exmouth Plateau, ODP Sites 122-762 and 122-763, eastern Indian Ocean

The evolution of planktonic foraminifera during the Late Cretaceous is marked in the Santonian by the disappearance of complex morphotypes (the marginotruncanids), and the contemporary increasing importance and diversification of another group of complex taxa, the globotruncanids. Upper Turonian to lower Campanian planktonic foraminiferal assemblages from Holes 762C and 763B (Ocean Drilling Program, Leg 122, Exmouth Plateau, 47°S palaeolatitude) were studied in detail to evaluate the compositional variations at the genus and species level based on the assumption that, in the Cretaceous oceans as in the modern, any faunal change was associated with changes in the characteristics and the degree of stability of the oceanic surface waters. Three major groups were recognised based on gross morphology, and following the assumption that Cretaceous planktonic foraminifera, although extinct, had life-history strategies comparable to those of modern planktonics: 1 - r-selected opportunists; 2 - k-selected specialists; 3 - r/k intermediate morphotypes which include all genera that display a range of trophic strategies in-between opportunist and specialist taxa. Although planktonic foraminiferal assemblages are characterised by a progressive appearance of complex taxa, this trend is discontinuous. Variation in number of species and specimens within genera has allowed recognition of five discrete intervals each of them reflecting different oceanic conditions based on fluctuations in diversity and abundance of the major morphotypes. Planktonic forms show cyclical fluctuations in diversity and abundance of cold (r-strategists) and warm taxa (k-strategists), perhaps representing alternating phases of unstable conditions (suggesting a weakly stratified upper water column in a mesotrophic environment), and well-stratified surface and near-surface waters (indicating a more oligotrophic environment). Interval 1, middle Turonian to early Coniacian in age, is dominated by the r/k intermediate morphotypes which alternate with r-strategists. These cyclical alternations are used to identify three additional subintervals. Interval 2, aged middle to late Coniacian, is characterised by the increasing number of species and relative abundance of k-strategists. After this maximum diversification the k-strategists show a progressive decrease reaching a minimum value in Interval 3 (early to late Santonian), which corresponds to the extinction of the genus Marginotruncana. In the Interval 4, latest Santonian in age, the k-strategists, represented mainly by the genera Globotruncana, increase again in diversity and abundance. The last Interval 5 (early Campanian) is dominated by juvenile globotruncanids and r-strategists which fluctuate in opposite phase. The positive peak (Interval 2) related to the maximum diversification of warm taxa (k-strategists) in the Coniacian seems to correspond to a warmer episode. It is followed by a marked decrease in the relative abundance of warm taxa (k-strategists crisis) with a minimum in the late Santonian (Interval 3), reflecting a decrease in temperature. Detailed analysis of faunal variations allows the Santonian faunal turnover to be ascribed to a cooling event strong enough to cause the extinction of the marginotruncanids.

Silt grain-size analysis of ODP holes from the western equatorial Atlantic

Grain-size records of the terrigenous and calcareous silt fraction, preservation of planktic foraminifera, and benthic foraminiferal stable-isotope data (delta13C, delta18O values of C. wuellerstorfi) at ODP Site 927 on the Ceará Rise (5°27.7'N, 44°28.8'W), are used to reconstruct variations in the history of bottom current strength, ventilation, and carbonate corrosiveness of deep waters during the time interval from 0.8 to 0.3 Ma. Glacial periods are characterized by generally smaller mean sizes of the terrigenous sortable silt fraction (mean(SS)), lower delta13C values, and poorer preservation of planktic foraminifera compared to interglacials. This indicates lower bottom current speeds, larger nutrient contents and more corrosive deep water. By contrast, larger mean(SS) sizes, higher delta13C values, and well preserved planktic foraminifera indicate strong circulation and a well ventilated deep-water mass during interglacials. The observed changes are most likely related to the weakening and strengthening of circulation of Lower North Atlantic Deep Water (LNADW). Cross-spectral analysis between the mean(SS) and benthic delta18O records reveals that minima in mean(SS) occur about 7.6 k.y. after the maximum in ice volume. This indicates a considerable lag time between ice-shield induced changes in LNADW production and subsequent changes in the velocity of LNADW flow in the western equatorial Atlantic. Striking changes in bottom current speed occur regularly during glacial to interglacial transitions. Extremely fine mean(SS) minima point to an almost complete shutdown of bottom current vigor in response to a cessation of LNADW production caused by an enhanced melt water release during the initial phases of deglaciation. However, each of the fine minima extremes is followed by a rapid shift to very high mean(SS) values that indicate strong bottom currents, and hence, vigorous LNADW flow during the early interglacials. After the onset of glacial Stage 12, generally poorer carbonate preservation and higher variability is registered. This coincides with a global decrease in carbonate preservation during the mid-Brunhes (mid-Brunhes dissolution event). Detailed grain-size analysis of the calcareous fine fraction (<63 µm) revealed a considerable reduction of particles in the fraction from 7 to 63 µm during periods of enhanced dissolution. This indicates a preferential dissolution of larger planktic foraminiferal fragments which leads to an enrichment of coccoliths in the calcareous fine fraction.

Calcareous nannofossil stratigraphy, mass accumulation rates and bulk density of sediments from ODP Leg 173 sites

Drilling on the Iberia Abyssal Plain during Ocean Drilling Program Leg 173 allowed us to recover Upper Cretaceous through Paleocene sediments at Sites 1068 and 1069 and only upper Paleocene sediments at Site 1067, which expands considerably the Upper Cretaceous to Paleocene record for this region. Of these three sites, Site 1068 recovered uppermost Cretaceous sediments as well as the most complete Paleocene record, whereas Site 1067 yielded only uppermost Paleocene sediments (Zone CP8). Site 1069 provided a rather complete upper Campanian through Maastrichtian section but a discontinuous Paleocene record. After a detailed calcareous nannofossil biostratigraphy was documented in distribution charts, we calculated mass accumulation rates for Holes 1068A and 1069A. Sediments in Hole 1068A apparently record the final stages of burial of a high basement block by turbidity flows. Accumulation rates through the Upper Cretaceous indicate relatively high rates, 0.95 g/cm**2/k.y., but may be unreliable because of the lack of datum points and/or possible hiatuses. Accumulation rates in the Paleocene section of Hole 1068A fluctuated every few million years from lower (~0.35 g/cm**2/k.y.) to higher rates (~0.85 g/cm**2/k.y.) until the latest Paleocene, when rates increased to an average of ~2.0 g/cm**2/k.y. Mass accumulation rates for the Upper Cretaceous in Hole 1069A indicate a steady rate of ~0.60 g/cm**2/k.y. from 75 to 72 Ma. There may have been one or more hiatuses between 72 and 68 Ma (combined Zone CC24 through Subzone CC25b), as indicated by the very low accumulation rate of 0.15 g/cm**2/k.y. The Paleocene section of Hole 1069A does not show the same continuous record, which may result from fluctuations in the carbonate compensation depth and poor recovery (average = 40%). Zones CP4 and CP5 are missing within a barren interval; this and numerous other barren intervals affect the precision of the nannofossil zonation and calculation of mass accumulation rates. However, in spite of these missing zones, mass accumulation rates do not seem to indicate the presence of hiatuses as the rates for this barren interval average ~1.0 g/cm**2/k.y. This study set out to test the hypothesis that a reliable biostratigraphic record could be constructed from sediments derived from turbidity flows deposited below the carbonate compensation depth. As illustrated here, not only could a reliable biostratigraphic record be determined from these sediments, but sedimentation and mass accumulation rates could also be determined, allowing inferences to be drawn concerning the sedimentary history of this passive margin. The reliability of this record is confirmed by independent verification by the establishment of a magnetostratigraphy for the same cores.