Dust fluxes and iron fertilization in Holocene and Last Glacial Maximum climates

Mineral dust aerosols play a major role in present and past climates. To date, we rely on climate models for estimates of dust fluxes to calculate the impact of airborne micronutrients on biogeochemical cycles. Here we provide a new global dust flux data set for Holocene and Last Glacial Maximum (LGM) conditions based on observational data. A comparison with dust flux simulations highlights regional differences between observations and models. By forcing a biogeochemical model with our new data set and using this model's results to guide a millennial-scale Earth System Model simulation, we calculate the impact of enhanced glacial oceanic iron deposition on the LGM-Holocene carbon cycle. On centennial timescales, the higher LGM dust deposition results in a weak reduction of <10?ppm in atmospheric CO2 due to enhanced efficiency of the biological pump. This is followed by a further ~10?ppm reduction over millennial timescales due to greater carbon burial and carbonate compensation.

Physical properties, grain size distribution and clay mineralogy of ODP Leg 117 sites

The fabric of sediments recovered at sites drilled on the Indus Fan, Owen Ridge, and Oman margin during Ocean Drilling Program Leg 117 was examined by scanning electron microscopy to document changes that accompany sediment burial. Two sediment types were studied: (1) biogenic sediments consisting of a variety of marly nannofossil and nannofossil oozes and chalks and (2) terrigenous sediments consisting of fine-grained turbidites deposited in association with the Indus Fan. Biogenic sediments were examined with samples from the seafloor to depths of 306 m below seafloor (mbsf) on the Owen Ridge (Site 722) and 368 mbsf on the Oman margin (Sites 723 and 728). Over these depth ranges the biogenic sediments are characterized by a random arrangement of microfossils and display little chemical diagenetic alteration. The microfossils are dispersed within a fine-grained matrix that is predominantly microcrystalline carbonate particles on the Owen Ridge and clay and organic matter on the Oman margin. Sediments with abundant siliceous microfossils display distinct, open fabrics with high porosity. Porosity reduction resulting from gravitational compaction appears to be the primary process affecting fabric change in the biogenic sediment sections. Fabric of illite-rich clayey silts and silty claystones from the Indus Fan (Site 720) and Owen Ridge (Sites 722 and 731) was examined for a composite section extending from 45 to 985 mbsf. In this section fabric of the fine-grained turbidites changes from one with small flocculated clay domains, random particle arrangement, and high porosity to a fabric with larger domains, strong preferred particle orientation roughly parallel to bedding, and lower porosity. These changes are accomplished by a growth in domain size, primarily through increasing face-to-face contacts, and by particle reorientation which is characterized by a sharp increase in alignment with bedding between 200 and 400 mbsf. Despite extensive particle reorientation, flocculated clay fabric persists in the deepest samples examined, particularly adjacent to silt grains, and the sediments lack fissility. Fabric changes over the 45-985 mbsf interval occur in response to gravitational compaction. Porosity reduction and development of preferred particle orientation in the Indus Fan and Owen Ridge sections occur at greater depths than outlined in previous fabric models for terrigenous sediments as a consequence of a greater abundance of silt and a greater abundance of illite and chlorite clays.

Age determination of sediment cores from the Panama Basin

High-resolution percent Corg and delta18Oforam records obtained from Panama Basin core Atlantis II 54-25PC and additional data from nearby core P7 show that enhanced burial of organic carbon has characterized every major glacial period for the last 500 kyr in that area. Both Corg concentration and mass accumulation rate profiles exhibit a sawtooth pattern with maxima occurring typically in the later stages of glacial periods. Comparison with dust records suggests that the carbon accumulation rate profile reflects both the upwelling history and a variable rate of iron input during the late Quaternary. The sawtooth character may derive from increased wind velocities and rates of upwelling during glacials which are indirectly related to ice volume (Sarnthein et al., 1988). The rapid decline in export production at the end of glacials in the equatorial Pacific may be attributed to the retreat of ice sheets (thus reduced wind velocities and upwelling) coupled with a coincident decline in atmospheric dust load and/or delivery rate. The Corg accumulation rate profiles do not correlate well with atmospheric CO2 records. For example, atmospheric CO2 was already at a minimum 40 kyr ago when production in the Panama Basin began increasing dramatically, commensurate with an increase in global dust levels. Using the relationship between the degree of photosynthetic fractionation and the concentration of free CO2 in the surface ocean postulated by Popp et al. (1989), delta13Corg measurements made on core P7 show that Panama Basin surface waters have been supplying CO2 to the atmosphere continually for at least the last 50 kyr. There is no evidence for a flux of CO2 into the surface ocean in this area at any time during this period despite the higher production. If the Panama Basin cores are representative of the eastern and central equatorial Pacific, then these observations weaken the influence on CO2 drawdown postulated for increased glacial productivity at low latitudes.

Stable carbon and oxygen isotope ratios of carbonates of ODP Site 103-639

A Tithonian sequence of shallow-water limestones, intercalated with siliciclastics and overlain by dolomite, was recovered during drilling at ODP Site 639 on the edge of a tilted fault block. The carbonates were strongly affected by fracturing, dolomitization, dedolomitization, and compaction. The chronology and nature of the fractures, fracture infilling, and diagenesis of the host rock are established and correlated for both the limestone and the dolomite. A first phase of dolomitization affected limestone that was already, at least partially, indurated. In the limestone unit, fractures were filled by calcite and dolomite; most of the dolomite was recrystallized into calcite, except for the upper part. In the dolomitic unit, the first-formed dolomite was progressively recrystallized into saddle dolomite, as fractures were simultaneously activated. The dolomitic textures become less magnesian (the molar ratio mMg/mCa goes from 1.04-0.98 to 0.80), and the d18O (PDB) ranges from -10 per mil to -8 per mil. The varying pores and fissures are either cemented by a calcic saddle dolomite (mMg/mCa ranging from 0.95 to 0.80) or filled with diverse internal sediments of detrital calcic dolomite, consisting of detrital dolomite silt (d18O from -9 per mil to -7 per mil) and laminated yellow filling (with different d18O values that range from -4 per mil to +3 per mil). These internal sediments clearly contain elements of the host rock and fragments of saddle crystals. They are covered by marls with calpionellids of early Valanginian age, which permits dating of most of the diagenetic phases as pre-Valanginian. The dolomitization appears to be related to fracturing resulting from extensional tectonics; it is also partially related to an erosional episode. Two models of dolomitization can be proposed from the petrographic characteristics and isotopic data. Early replacement of aragonite bioclasts by sparite, dissolution linked to dolomitization, and negative d18O values of dolomite suggest a freshwater influence and 'mixing zone' model. On the other hand, the significant presence of saddle dolomite and repeated negative d18O values suggest a temperature effect; because we can dismiss deep burial, hydrothermal formation of dolomite would be the most probable model. For both of these hypotheses, the vadose filling of cavities and fractures by silt suggests emersion, and the different, and even positive, d18O values of the last-formed yellow internal sediment could suggest dolomitization of the top of the sequence under saline to hypersaline conditions. Fracturing resulting in the reopening of porosity and the draining of dolomitizing fluids was linked to extensional tectonics prior to the tilting of the block. These features indicate an earlier beginning to the rifting of the Iberian margin than previously known. Dolomitization, emersion, and erosion correspond to eustatic sea-level lowering at the Berriasian/Valanginian boundary. Diagenesis, rather than sedimentation, seems to mark this global event and to provide a record of the regional tectonic history.

Organic carbon analyses of ODP Site 160-971 in the Mediterranean Sea

Mud volcanism on the Mediterranean Ridge is caused by extrusion of overpressured sediments, with consequent formation of spectacular dome-shaped features composed of mud breccias at the seafloor. The organic material in the mud breccia of the Napoli mud volcano is a mixture of different facies, stratigraphic origin and thermal maturities. One portion is synsedimentary organic material with only minor diagenetic alterations and represents sedimenting material that was embedded into the mud volcano during its extrusion. The mud breccia also contains thermally mature organic material of mainly terrestrial provenance with algae of fresh- and brackish-water origin. Vitrinite reflectance data of this maturity generation range from 0.65 to 0.90% R(oil) and thus characterize thermally mature source rocks, a rank which is corroborated by fluorescence and molecular characteristics. The predominance of vitrinite in the maceral assemblages and the occurrence of biomarkers of terrigenous origin suggest that the major part of the mud matrix derives from a lacustrine or riverine sedimentary unit in the subsurface, possibly from the Messinian stage. A third generation of organic material includes inertinites and vitrinites of high reflectance, which represent recycled organic matter present in any marine sediment. By use of the Lopatin method for modelling the thermal maturation of hydrocarbon source rocks from the vitrinite reflectance data, we calculated that the depth of mobilization ranges from 4900 m to 7500 m, depending upon the temperature gradient used.

Geochemistry of Neogene-Quaternary sediments in ODP Site 105-645

Subcontinuously cored early(?)-middle Miocene to recently deposited sediments from ODP Site 645 were studied texturally, mineralogically, and geochemically. The entire sequence contains minerals and associated chemical elements that are chiefly of detrital origin. In particular, the clay minerals, which include smectite, kaolinite, chlorite, and illite, are detrital. No obvious evidence of diagenesis with depth, of burial, of volcanism, or of hydrothermal alteration was observed. The sedimentary textures, clay mineralogy, and <2-µm fraction geochemistry of the early middle Miocene sediments (630 to 1147 mbsf) suggest the pronounced but variable influence of a southward bottom current. Two clay facies are defined. The lower one, Cj (780 to 1147 mbsf), is characterized by the great abundance of discrete smectite (with less than 15% illite interlayers), probably detrital in origin, and reworked older, discrete, smectite-rich sediments. The upper clay facies, C2 (630 to 780 mbsf), shows a net decrease of the fully expandable clay abundances, with a great abundance of mixed-layer, illite-smectite clays (60 to 80% of illite interlayers). Such clay assemblages can be inherited from paleosoils or older sedimentary rocks. An important change occurs at 630 mbsf (clay fraction) or 600 mbsf (sedimentary texture), which may be explained by the beginning of continental glaciation (630 mbsf, ~9 Ma) and the onset of ice rafting in Baffin Bay (600 mbsf, ~8 Ma). Above this level, the characteristics and modifications of the clay assemblages are controlled climatically and can be explained by the fluctuations of (1) ice-rafting, (2) speed of weak bottom currents, and (3) some supply by mud turbiditic currents. Three clay facies (C3, C4, and C5) can be defined by the abrupt increases of the inherited chlorite and illite clays.

Helium isotope concentrations and ratios in sediments from ODP Hole 130-806B

We have determined the helium abundance and isotopic composition of seafloor carbonate sediments from the flanks of the Ontong Java Plateau, western equatorial Pacific Ocean (ODP Site 806). These results provide a two million year record of the burial flux of extraterrestrial 3He, which we believe is a proxy for the terrestrial accretion rate of interplanetary dust particles. The 3He burial flux prior to ~700 ka was relatively low, ~0.5 pcc/cm**2/kyr, but from 700 ka to the present, the burial flux gradually increased to a value of ~1.0 pcc/cm**2/kyr. 100 kyr periodicity in the 3He burial flux is apparent over the last 700 kyr and correlates with the oxygen isotope record of global climate, with high 3He burial fluxes associated with interglacial periods. This periodicity and phase are consistent with previous 3He measurements in North Atlantic sediments. Although 100 kyr periodicity in 3He burial flux is in agreement with recent predictions of the accretion rate of interplanetary dust based on a model of the orbital evolution of asteroidal debris, the measurements and predictions differ by one half cycle in phase. Nevertheless, our observations suggest the terrestrial accretion rate of interplanetary dust is controlled by orbital eccentricity and/or inclination relative to the solar-system invariable plane. Such control is a necessary but not sufficient condition for the hypothesis of that variations in extraterrestrial dust accretion modulates terrestrial climate with a 100 kyr period. We also identify several brief (<25 kyr) intervals of strongly enhanced 3He burial, possibly related to random and transient fluctuations in the accretion rate of asteroidal or cometary dust particles.

Oxygen-isotope composition, temperature and depth of formation of chert at DSDP Leg 62 Holes

We measured oxygen-isotope compositions of 16 siliceous rocks from Deep Sea Drilling Project Sites 463, 464, 465, and 466 (Leg 62). Samples are from deposits that range in age from about 40 to 103 m.y. and that occur at sub-bottom depths of 9 to 461 meters. Mean d18O values range from 28.4 to 36.8 per mil and 36.0 ± 0.3 per mil for quartz-rich and opal-CTrich rocks, respectively. d18O values in chert decrease with increasing sub-bottom depth; the slope of the d18O/depth curve is less steep for Site 464 than for the other sites which indicates that chert at Site 464 formed at higher temperatures than chert at Sites 463, 465, and 466. Temperatures of formation of cherts were 7 to 42°C, using the silica-water fractionation factor of Knauth and Epstein (1976), or 19 to 56°C, using the equation of Clayton et al. (1972). Temperatures in the sediment where the cherts now occur are lower than their isotopically determined temperatures of formation, which means that the cherts record an earlier history when temperatures in the sediment section were greater. Estimated sediment temperatures when the cherts formed are comparable to, but generally slightly lower than, those calculated from Knauth and Epstein's equation. The isotopic composition of cherts is more closely related to environment of formation (diagenetic environment) or paleogeothermal gradients, than to paleoclimates (bottom-water temperatures). Opal-CT-rich rocks may better record paleo-bottom-water temperature. In Leg 62 cherts, better crystallinity of quartz corresponds to lower d18O values; this implies progressively higher temperatures of equilibration between quartz and water during maturation of quartz. The interrelationship of d18O and crystallinity is noted also in continental-margin deposits such as the Monterey Formation - but for higher temperatures. The apparent temperature difference between open-ocean and continental-margin deposits can be explained by the dominant control of temperature on silica transformation in the rapidly deposited continental-margin deposits, whereas time, as well as temperature, has a strong influence on the transformations in open-ocean deposits. Comparisons between the chemistry and d18O values of cherts reveal two apparent trends: both boron and SiO2 increase as d18O increases. However, the correspondence between SiO2 and d18O is only apparent, because the two cherts lowest in SiO2 are also the most deeply buried, so the trend actually reflects depth of burial. The correspondence between boron and d18O supports the conclusion that boron is incorporated in the quartz crystal structure during precipitation

Concentrations and accumulation rates of platinum-group elements and rhenium in Pacific and Atlantic sediments, DSDP and ODP data

The nature of Re-platinum-group element (PGE; Pt, Pd, Ir, Os, Ru) transport in the marine environment was investigated by means of marine sediments at and across the Cretaceous-Tertiary boundary (KTB) at two hemipelagic sites in Europe and two pelagic sites in the North and South Pacific. A traverse across the KTB in the South Pacific pelagic clay core found elevated levels of Re, Pt, Ir, Os, and Ru, each of which is approximately symmetrically distributed over a distance of ~1.8 m across the KTB. The Re-PGE abundance patterns are fractionated from chondritic relative abundances: Ru, Pt, Pd, and Re contents are slightly subchondritic relative to Ir, and Os is depleted by ~95% relative to chondritic Ir proportions. A similar depletion in Os (~90%) was found in a sample of the pelagic KTB in the North Pacific, but it is enriched in Ru, Pt, Pd, and Re relative to Ir. The two hemipelagic KTB clays have near-chondritic abundance patterns. The ~1.8-m-wide Re-PGE peak in the pelagic South Pacific section cannot be reconciled with the fallout of a single impactor, indicating that postdepositional redistribution has occurred. The elemental profiles appear to fit diffusion profiles, although bioturbation could have also played a role. If diffusion had occurred over ~65 Ma, the effective diffusivities are ~10**?13 cm**2/s, much smaller than that of soluble cations in pore waters (~10**?6 cm**2/s). The coupling of Re and the PGEs during redistribution indicates that postdepositional processes did not significantly fractionate their relative abundances. If redistribution was caused by diffusion, then the effective diffusivities are the same. Fractionation of Os from Ir during the KTB interval must therefore have occurred during aqueous transport in the marine environment. Distinctly subchondritic Os/Ir ratios throughout the Cenozoic in the South Pacific core further suggest that fractionation of Os from Ir in the marine environment is a general process throughout geologic time because most of the inputs of Os and Ir into the ocean have Os/Ir ratios >/=1. Mass balance calculations show that Os and Re burial fluxes in pelagic sediments account for only a small fraction of the riverine Os (<10%) and Re (<0.1%) inputs into the oceans. In contrast, burial of Ir in pelagic sediments is similar to the riverine Ir input, indicating that pelagic sediments are a much larger repository for Ir than for Os and Re. If all of the missing Os and Re is assumed to reside in anoxic sediments in oceanic margins, the calculated burial fluxes in anoxic sediments are similar to observed burial fluxes. However, putting all of the missing Os and Re into estuarine sediments would require high concentrations to balance the riverine input and would also fail to explain the depletion of Os at pelagic KTB sites, where at most ~25% of the K-T impactor's Os could have passed through estuaries. If Os is preferentially sequestered in anoxic marine environments, it follows that the Os/Ir ratio of pelagic sediments should be sensitive to changes in the rates of anoxic sediment deposition. There is thus a clear fractionation of Os and Re from Ir in precipitation out of sea water in pelagic sections. Accordingly, it is inferred here that Re and Os are removed from sea water in anoxic marine depositional regimes.

Geochemistry and accumulation rates of the early Paleogene in the Weddel Sea, Maud Rise and Kerguelen Plateau

Marine biological productivity has been invoked as a possible climate driver during the early Paleogene through its potential influence on atmospheric carbon dioxide concentrations. However, the relationship of export productivity (the flux of organic carbon (C) from the surface ocean to the deep ocean) to organic C burial flux (the flux of organic C from the deep ocean that is buried in marine sediments) is not well understood. We examine the various components involved with atmosphere-to-ocean C transfer by reconstructing early Paleogene carbonate and silica production (using carbonate and silica mass accumulation rates (MARs)); export productivity (using biogenic barium (bio-Ba) MARs); organic C burial flux (using reactive phosphorus (P) MARs); redox conditions (using uranium and manganese contents); and the fraction of organic C buried relative to export productivity (using reactive P to bio-Ba ratios). Our investigations concentrate on Paleocene/Eocene sections of Sites 689/690 from Maud Rise and Site 738 from Kerguelen Plateau. In both regions, export productivity, organic C burial flux, and the fraction of organic C buried relative to export productivity decreased from the Paleocene/early Eocene to the middle Eocene. A shift is indicated from an early Paleogene two-gyre circulation in which nutrients were not efficiently recycled to the surface via upwelling in these regions, to a circulation more like the present day with efficient recycling of nutrients to the surface ocean. Export productivity was enhanced for Kerguelen Plateau relative to Maud Rise throughout the early Paleogene, possibly due to internal waves generated by the plateau regardless of gyre circulation.

Ba distribution in surface sediments of the Southern Ocean

We present excess Ba (Baxs) data (i.e., total Ba corrected for lithogenic Ba) for surface sediments from a north-south transect between the Polar Front Zone and the northern Weddell Gyre in the Atlantic sector and between the Polar Front Zone and the Antarctic continent in the Indian sector. Focus is on two different processes that affect excess Ba accumulation in the sediments: sediment redistribution and excess Ba dissolution. The effect of these processes needs to be corrected for in order to convert accumulation rate into vertical rain rate, the flux component that can be linked to export production. In the Southern Ocean a major process affecting Ba accumulation rate is sediment focusing, which is corrected for using excess 230Th. This correction, however, may not always be straightforward because of boundary scavenging effects. A further major process affecting excess Ba accumulation is barite dissolution during exposure at the sediment-water column interface. Export production estimates derived from excess 230Th and barite dissolution corrected Baxs accumulation rates (i.e., excess Ba vertical rain rates) are of the same magnitude but generally larger than export production estimates based on water column proxies (234Th-deficit in the upper water column; particulate excess Ba enrichment in the mesopelagic water column). We believe export production values based on excess Ba vertical rain rate might be overestimated due to inaccurate assessment of the Baxs preservation rate. Barite dissolution has, in general, been taken into account by relating it to exposure time before burial depending on the rate of sediment accumulation. However, the observed decrease of excess Ba content with increasing water column depth (or increasing hydrostatic pressure) illustrates the dependence of barite preservation on degree of saturation in the deep water column in accordance with available thermodynamic data. Therefore correction for barite dissolution would not be appropriate by considering only exposure time of the barite to some uniformly undersaturated deep water but requires also that regional differences in degree of undersatuation be taken into account.

Planktonic foraminifers from Oligocene to Pleistocene sediments of DSDP Leg 92 sites

Leg 92 of the Deep Sea Drilling Project cored sediments containing calcareous microfossils at six sites along 19°S latitude in the South Pacific Ocean. Shipboard examination of these sediments revealed planktonic foraminifers of uppermost Oligocene through Pleistocene age that were identified and assigned to biostratigraphic zones according to the tropical zonation scheme of Blow (1969). Preservation of planktonic foraminifers in the sites from Leg 92 has been affected by the position of each site with respect to the lysocline and calcium carbonate compensation depth (CCD) at the time of deposition, depth of burial, and sediment accumulation rate (rate of burial). An additional factor may also be important, especially in the sediments deposited immediately above basement. Evidence of poor preservation in basal sediments of Holes 600C and 601, which have always been shallower than both the lysocline and the CCD, suggests that hydrothermal solutions circulating within young oceanic crust may penetrate the overlying sediments and affect the preservation of calcareous microfossils deposited there.

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.