Age model and Plio-Pleistocene benthic stable oxygen isotope ratios of ODP Site 184-1143 in the Sourh China Sea

Based on benthic foraminiferal delta18O from ODP Site 1143, a 5-Myr astronomical timescale for the West Pacific Plio-Pleistocene was established using an automatic orbital tuning method. The tuned Brunhes/Matuyama paleomagnetic polarity reversal age agrees well with the previously published age of 0.78 Ma. The tuned ages for several planktonic foraminifer bio-events also agree well with published dates, and new ages for some other bio-events in the South China Sea were also estimated. The benthic delta18O from Site 1143 is highly coherent with the Earth's orbit (ETP) both at the obliquity and precession bands for the last 5 Myr, and at the eccentricity band for the last 2 Myr. In general, the 41-kyr cycle was dominant through the Plio-Pleistocene although the 23-kyr cycle was also very strong. The 100-kyr cycle became dominant only during the last 1 Myr. A comparison of the benthic delta18O between the Atlantic (ODP 659) and the East and West Pacific (846 and 1143) reveals that the Atlantic-Pacific benthic oxygen isotope difference ratio (Delta delta18OAtl-Pac) displays an increasing trend in three time intervals: 3.6-2.7 Ma, 2.7-2.1 Ma and 1.5-0.25 Ma. Each of the intervals begins with a rapid negative shift in Delta delta18OAtl-Pac, followed by a long period with an increasing trend, corresponding to the growth of the Northern Hemisphere ice sheet. This means that all three intervals of ice sheet growth in the Northern Hemisphere were accompanied at the beginning by a rapid relative warming of deep water in the Atlantic as compared to that of the Pacific, followed by its gradual relative cooling. This general trend, superimposed on the frequent fluctuations with glacial cycles, should yield insights into the processes leading to the boreal glaciation. Cross-spectral analyses of the Delta delta18OAtl-Pac with the Earth's orbit suggests that after the initiation of Northern Hemisphere glaciation at about 2.5 Ma, obliquity rather than precession had become the dominant force controlling the vertical structure or thermohaline circulation in the paleo-ocean.

Stable oxygen isotope record and lithogenic composition of ODP Hole 117-722B

The upper 38 m of Hole 722B sediments (Owen Ridge, northwest Arabian Sea) was sampled at 20 cm intervals and used to develop records of lithogenic percent, mass accumulation rate, and grain size spanning the past 1 m.y. Over this interval, the lithogenic component of Owen Ridge sediments can be used to infer variability in the strength of Arabian Sea summer monsoon winds (median grain size) and the aridity of surrounding dust source-areas (mass accumulation rate; MAR in g/cm**2/k.y). The lithogenic MAR has strong 100, 41, and 23 k.y. cyclicities and is forced primarily by changes in source-area aridity associated with glacial-interglacial cycles. The lithogenic grain size, on the other hand, exhibits higher frequency variability (23 k.y.) and is forced by the strength of summer monsoon winds which, in turn, are forced by the effective sensible heating of the Indian-Asian landmass and by the availability of latent heat from the Southern Hemisphere Indian Ocean. These forcing mechanisms combine to produce a wind-strength record which has no strong relationship to glacial-interglacial cycles. Discussion of the mechanisms responsible for production of primary Milankovitch cyclicities in lithogenic records from the Owen Ridge is presented elsewhere (Clemens and Prell, 1990, doi:10.1029/PA005i002p00109). Here we examine the 1 m.y. record from Hole 722B focusing on different aspects of the lithogenic components including an abrupt change in the monsoon wind-strength record at 500 k.y., core-to-core reproducibility, comparison with magnetic susceptibility, coherency with a wind-strength record from the Pacific Ocean, and combination frequencies in the wind-strength record. The Hole 722B lithogenic grain-size record shows an abrupt change at 500 k.y. possibly indicating decreased monsoon wind-strength over the interval from 500 k.y. to present. The grain-size decrease appears to be coincident with a loss of spectral power near the 41 k.y. periodicity. However, the grain-size decrease is not paralleled in the Globigerina bulloides upwelling record, an independent record of summer monsoon wind-strength (Prell, this volume). These observations leave us with competing hypotheses possibly involving: (1) a decrease in the sensitivity of monsoon windstrength to obliquity forcing, (2) decoupling of the grain size and G. bulloides records via a decoupling of the nutrient supply from wind-driven upwelling, and/or (3) a change in dust source-area or the patterns of dust transporting winds. Comparison of the lithogenic grain size and weight percent records from Hole 722B with those from a nearby core shows that the major and most minor events are well replicated. These close matches establish our confidence in the lithogenic extraction techniques and measurements. Further, reproducibility on a core-to-core scale indicates that the eolian depositional signal is regionally strong, coherent, and well preserved. The lithogenic weight percent and magnetic susceptibility are extremely well correlated in both the time and frequency domains. From this we infer that the magnetically susceptible component of Owen Ridge sediments is of terrestrial origin and transported to the Owen Ridge via summer monsoon winds. Because of the high correlation with the lithogenic percent record, the magnetic susceptibility record can be cast in terms of lithogenic MAR and used as a high resolution proxy for continental aridity. In addition to primary Milankovitch periodicities, the Hole 722B grain-size record exhibits periodicity at 52 k.y. and at 29 k.y. Both periodicities are also found in the grain-size record from piston core RC11-210 in the equatorial Pacific Ocean. Comparison of the two grain-size records shows significant coherence and zero phase relationships over both the 52 and 29 k.y. periodicities suggesting that the strengths of the Indian Ocean monsoon and the Pacific southeasterly trade winds share common forcing mechanisms. Two possible origins for the 52 and 29 k.y. periodicities in the Hole 722B wind-strength record are (1) direct Milankovitch forcing (54 and 29 k.y. components of obliquity) and (2) combination periodicities resulting from nonlinear interactions within the climate system. We find that the 52 and 29 k.y. periodicities show stronger coherency with crossproducts of eccentricity and obliquity (29 k.y.) and precession and obliquity (52 k.y.) than with direct obliquity forcing. Our working hypothesis attributes these periodicities to nonlinear interaction between external insolation forcing and internal climatic feedback mechanisms involving an interdependence of continental snow/ice-mass (albedo) and the hydrological cycle (latent heat availability).

Stable oxygen isotope record and relative abundances of planktonic foraminifera of ODP Hole 117-728A

High resolution stratigraphy based on oxygen isotope ratios of the planktonic foraminifers Neogloboquadrina dutertrei (d'Orbigny), Globigeriniodes ruber (d'Orbigny), and Globigerina bulloides (d'Orbigny), magnetic susceptibility, and calcium carbonate content covers the sedimentary record of ODP Hole 728A drilled on the Oman Margin from approximately 10 k.y. to 525 k.y., comprising isotopic stages 1-13. Below stage 13 isotopic stage boundaries cannot be defined with certainty in our data. Sediment accumulation rates were calculated from the isotopic record of N. dutertrei by matching it with the age model SPECMAP curve. During the glacial periods sediment accumulation rates were higher than during the interglacial periods, reflecting increased input from the shelf during low-stands of sea level and increased eolian input. Periodograms for the past 524 k.y. on oxygen isotope records of N. dutertrei, G. ruber, and G. bulloides, on calcium carbonate content, magnetic susceptibility, and on a foraminiferal fragmentation record show powers matching the Milankovitch periodicities. High powers are concentrated around 103 k.y. In the spectra of oxygen isotope ratios of N. dutertrei, magnetic susceptibility, and foraminiferal fragmentation these are significant at the 80% confidence level with respect to a first order autoregressive model. Power concentrations near 43 k.y., matching obliquity, are present but subdued in all spectra. Power concentrations near 23 k.y., matching precession, are significant in the spectra of the oxygen isotope record of N. dutertrei, magnetic susceptibility, and calcium carbonate content record. Fragmentation of planktonic foraminifers increased during the interglacial periods. This is attributed to dissolution of the tests in an expanded oxygen minimum zone (OMZ), where undersaturation of calcium carbonate is caused by enhanced production in the euphotic zone, which would suggest stronger monsoonal induced upwelling during interglacial periods. Extension of the OMZ could also be increased by outflow of low oxygen marginal basin bottom water.

Pleistocene sedimentation in the Benguela Current, southeastern Atlantic

We infer variations in paleoproductivity and eolian input at ODP Site 1082 in the Walvis Basin from stable oxygen isotope compositions of the planktonic foraminifera Globorotalia inflata, total organic carbon mass accumulation rates (TOC MAR), and X-ray fluorescence analyses of Fe content. The most pronounced paleoclimatic changes correspond to the time at about 0.9 Ma, when glacial conditions in the northern hemisphere (NH) led to the onset of pronounced 100-kyr glacial-interglacial cycles. We used Fe intensity as a proxy for eolian terrigenous input, and TOC MAR as a paleoproductivity indicator. Paleoproductivity and eolian input show generally higher-amplitude variations of glacial-interglacial cyclicity from 1.5 to 0.58 Ma, indicating pronounced variations in upwellingfavorable winds in this area. At 0.58 Ma, paleoproductivity and eolian input shifted abruptly to lower-amplitude variations with a periodicity of 100 kyr while delta18O values show a trend toward more negative isotope values for the past 0.65 Myr. Especially during glacial periods, oxygen isotope values indicate increasingly warmer sea-surface temperatures toward the end of the Pleistocene. To evaluate the relative influences of NH glaciation and southern hemisphere (SH) insolation as potential forcing mechanisms for variations of eolian input and productivity in the northern Benguela system, we filtered our proxy records at orbital frequencies. The filtered records of Fe intensity and TOC MAR indicate a strong influence of the 100-kyr and 41-kyr frequency bands, supporting our assumption that strong ice buildup in the NH is the dominant trigger for climate changes on the continent and probably in trade-wind intensity. SH insolation and low-latitude precession-related insolation changes were important for paleoproductivity variations in the northern Benguela system, modifying the nutrient supply by southern ocean intermediate waters and the zonal direction of upwelling-inducing trades by the African monsoon system, respectively.

Constructing a stacked benthic d18O record

Composite stacks were constructed by superimposing 6 to 13 benthic foraminiferal d18O records covering the period 0-850 ka. An initial timescale for each core was established using radioisotopic age control points and assuming constant sedimentation rates between these points. The average of these records is our 13-core "untuned" stack. Next, we matched the 41 kyr component of each record individually to variations in Earth's obliquity. Four of the 13 records produced timescales that were inconsistent with one or more of the known radioisotopic ages. The nine remaining cores were averaged to create a "minimally tuned" stack. Six of the minimally tuned cores were assembled into a "tropical" stack. For each stack we estimated the uncertainty envelope from the standard deviation of the constituents. Spectral analysis of the three stacks indicates that benthic d18O is dominated by a 100 kyr oscillation that has a narrow spectral peak. The contribution of precession to the total variance is small when compared to prior results from planktic stacks.

Proxy data and model simulations of sediment core GeoB7920-2

The evolution of the northwest African hydrological balance throughout the Pleistocene epoch influenced the migration of prehistoric humans**1. The hydrological balance is also thought to be important to global teleconnection mechanisms during Dansgaard-Oeschger and Heinrich events**2. However, most high-resolution African climate records do not span the millennial-scale climate changes of the last glacial-interglacial cycle**1, 3, 4, 5, or lack an accurate chronology**6. Here, we use grain-size analyses of siliciclastic marine sediments from off the coast of Mauritania to reconstruct changes in northwest African humidity over the past 120,000 years. We compare this reconstruction to simulations of palaeo-humidity from a coupled atmosphere-ocean-vegetation model. These records are in good agreement, and indicate the reoccurrence of precession-forced humid periods during the last interglacial period similar to the Holocene African Humid Period. We suggest that millennial-scale arid events are associated with a reduction of the North Atlantic meridional overturning circulation and that millennial-scale humid events are linked to a regional increase of winter rainfall over the coastal regions of northwest Africa.

XRF iron (Fe) counts of ODP Hole 171-1050A and 171-1051A

Upper Paleocene to lower Eocene sediments drilled at Ocean Drilling Program (ODP) Site 1051 (Blake Nose, off Florida) display well-defined orbital cycles, a detailed magnetic stratigraphy, and a suite of planktonic foraminiferal datums. We derived a cyclostratigraphy by using spectral analysis of high-resolution records of elemental concentrations obtained by an X-ray fluorescence (XRF) Core Scanner. XRF counts of iron serve as a proxy for the relative amount of terrestrial material. Sliding-window spectral analysis, bandpass filtering, and direct counting of precession and obliquity cycles yield minimum durations for magnetic polarity chrons C22 to C26 (~49 to ~61 Ma), calculations of sediment accumulation rates, as well as constraints on the timing of biostratigraphic and climatological events in the vicinity of the Initial Eocene Thermal Maximum (IETM). Durations of polarity chrons as represented in sediments drilled at Site 1051 were estimated using a conservative assignment of 41 k.y. for obliquity cycles and 21 k.y. for precession cycles. Combined polarity chrons C26r and C26n span 3.61 m.y., and chron C25r spans 1.07 m.y. Polarity chron C24r is estimated as 2.877 m.y. The interpretation of polarity chron C24n is ambiguous, but its duration is probably <1.23 m.y. Polarity chron C23r spans 0.53 m.y., chron C23n is 0.74 m.y., and chron C22r is 0.9 m.y. Spectral analysis through this interval indicates that spectral peaks shift through time and are related to changes in sedimentation rate in Site 1051. The sedimentation rates dramatically increased ~200 k.y. after the IETM and remained high for most of chron C24r.

Age determination of Southeast Atlantic sediment cores

Angola Basin and Cape Basin (southeast Atlantic) surface sediments and sediment cores show that maxima in the abundance of taraxerol (relative to other land-derived lipids) covary with maxima in the relative abundance of pollen from the mangrove tree genus Rhizophora and that in the surface sediments offshore maxima in the relative abundance of taraxerol occur at latitudes with abundant coastal mangrove forests. Together with the observation that Rhizophora mangle and Rhizophora racemosa leaves are extraordinarily rich in taraxerol, this strongly indicates that taraxerol can be used as a lipid biomarker for mangrove input to the SE Atlantic. The proxy-environment relations for taraxerol and Rhizophora pollen down-core show that increased taraxerol and Rhizophora pollen abundances occur during transgressions and periods with a humid climate. These environmental changes modify the coastal erosion and sedimentation patterns, enhancing the extent of the mangrove ecosystem and/or the transport of mangrove organic matter offshore. Analyses of mid-Pleistocene sediments show that interruption of the pattern of taraxerol maxima during precession minima occurs almost only during periods of low obliquity. This demonstrates the complex environmental response of the interaction between precession-related humidity cycles and obliquity-related sea-level changes on mangrove input.

Age models for ODP Leg 172 sites

Ten ODP sites drilled in a depth transect (2164-4775 m water depth) during Leg 172 recovered high-deposition rate (>20 cm/kyr) sedimentary sections from sediment drifts in the western North Atlantic. For each site an age model covering the past 0.8-0.9 Ma has been developed. The time scales have a resolution of 10-20 kyr and are derived by tuning variations of estimated carbonate content to the orbital parameters precession and obliquity. Based on the similarity in the signature of proxy records and the spectral character of the time series, the sites are divided into two groups: precession cycles are better developed in carbonate records from a group of shallow sites (2164-2975 m water depth, Sites 1055-1058) while the deeper sites (2995-4775 m water depth, Sites 1060-1063) are characterized by higher spectral density in the obliquity band. The resulting time scales show excellent coherence with other dated carbonate and isotope records from low latitudes. Besides the typical Milankovitch cyclicity significant variance of the resulting carbonate time series is concentrated at millennial-scale changes with periods of about 12, 6, 4, 2.5, and 1.5 kyr. Comparisons of carbonate records from the Blake Bahama Outer Ridge and the Bermuda Rise reveal a remarkable similarity in the time and frequency domain indicating a basin-wide uniform sedimentation pattern during the last 0.9 Ma.

Delta 66Zn ratios and age of carbonates from ODP Site 138-849 (Table 1)

The carbonate fraction of sediment core ODP 849, leg 138, located in the eastern equatorial Pacific, mostly consisting of coccoliths, was separated and analyzed for its Zn isotopic composition. The overall variation in Zn isotopic composition, as determined by multiple-collector, magnetic-sector, inductively coupled plasma mass spectrometry, was found to be on the order of 1? (expressed in delta66Zn, where deltaxZn=[(xZn/64Zn)sample/(xZn/64Zn)standard -1]*10**3 and x=66, 67 or 68) over the last 175 ka. The analytical precision was 0.04 per mil and the overall reproducibility was usually better than 0.07 per mil. The Zn isotopic composition signal exhibits several marked peaks and a high-frequency variability. A periodogram of the delta66Zn signal showed two periodicities of 35.2 and 21.2 ka. We suggest that the latter is caused by the precession of the Earth's axis of rotation. The periodogram exhibits a minimum at 41.1 ka, thus showing that the Zn isotopic composition is independent of the obliquity in the eastern equatorial Pacific. The range of delta66Zn values observed for the carbonate fraction of ODP 849 overlaps with the range observed for Fe-Mn nodules in the world's oceans, which suggests that seawater/carbonate Zn isotope fractionation is weak. We therefore assume that most of the Zn isotope variability is a result of the selective entrainment of the light isotopes by organic matter in the surface ocean. The ODP 849 delta66Zn record seems to follow the changes in the insolation cycles. Changes in the late summer/fall equatorial insolation modulate the intensity of the equatorial upwelling, hence the mixing between deep and surface waters. We propose that during decreased summer/fall equatorial insolation, when a steep thermocline can develop (El Niño-like conditions), the surface waters cannot be replenished by deep waters and become depleted in the lighter Zn isotopes by biological activity, thus resulting in the progressive increase of the delta66Zn values of the carbonate shells presumably in equilibrium with surface seawater.

Stable oxygen isotope. Mg/Ca and Temperature record of Caribbean and Pacific sediments

The final phase of the closure of the Panamanian Gateway and the intensification of Northern Hemisphere Glaciation (NHG) both occurred during the Late Pliocene. Glacial-interglacial (G-IG) variations in sea level might, therefore, have had a significant impact on the remaining connections between the East Pacific and the Caribbean. Here, we present combined foraminiferal Mg/Ca and d18O measurements from Ocean Drilling Program (ODP) Site 1241 from the East Pacific and ODP Site 999 from the Caribbean. The studied time interval covers the first three major G-IG Marine Isotope Stages (MIS 95-100, ~2.5 Ma) after the intensification of NHG. Analyses were performed on the planktonic foraminifera Neogloboquadrina dutertrei and Globigerinoides sacculifer, representing water mass properties in the thermocline and the mixed-layer, respectively. Changes in sea water temperature, relative salinity, and water column stratification strongly suggest that the Panamanian Gateway temporarily closed during glacial MIS 98 and 100, as a result of changes in ice volume equivalent to a drop in sea level of 60-90 m. Reconstructed sea surface temperatures (SST) from G. sacculifer show a glacial decrease of 2.5°C at Site 1241, but increases of up to 3°C at Site 999 during glacial MIS 98 and 100 suggesting that the Panamanian Gateway closed during these glacial periods. The Mg/Ca-temperatures of N. dutertrei remain relatively stable in the East Pacific, but do show a 3°C warming in the Caribbean at the onset of these glacial periods suggesting that the closing of the gateway also changed the water column stratification. We infer that the glacial closure of the gateway allowed the Western Atlantic Warm Pool to extend into the southern Caribbean, increasing SST (G. sacculifer) and deepening the thermocline (N. dutertrei). Additionally, ice volume appears to have become large enough during MIS 100 to survive the relatively short lasting interglacial MIS 99 so that the gateway remained closed. Towards the end of MIS 98, during MIS 97 and into MIS 96 temperatures on both sides are mostly similar suggesting water masses exchanged again. Additionally, Caribbean variations in SST and d18Owater follow a precession-like cyclicity rather than the obliquity-controlled variations characteristic of the East-Pacific and many other tropical areas, suggesting that regional atmospheric processes related to the trade winds and the Intertropical Convergence Zone (ITCZ) had a dominant impact in the Caribbean.

Sediment, isotopic and geochemical record during Marine Isotope Stages 29 to 34 of IODP Site 339-U1387

Centennial-to-millennial scale records from IODP Site U1387, drilled during IODP Expedition 339 into the Faro Drift at 558 m water depth, now allow evaluating the climatic history of the upper core of the Mediterranean Outflow (MOW) and of the surface waters in the northern Gulf of Cadiz during the early Pleistocene. This study focuses on the period from Marine Isotope Stage (MIS) 29 to 34, i.e. the interval surrounding extreme interglacial MIS 31. Conditions in the upper MOW reflect obliquity, precession and millennial-scale variations. The benthic d18O signal follows obliquity with the exception of an additional, smaller d18O peak that marks the MIS 32/31 transition. Insolation maxima (precession minima) led to poor ventilation and a sluggish upper MOW core, whereas insolation minima were associated with enhanced ventilation and often also increased bottom current velocity. Millennial-scale periods of colder sea-surface temperatures (SST) were associated with short-term maxima in flow velocity and better ventilation, reminiscent of conditions known from MIS 3. A prominent contourite layer, coinciding with insolation cycle 100, was formed during MIS 31 and represents one of the few contourites developing within an interglacial period. MIS 31 surface water conditions were characterized by an extended period (1065-1091 ka) of warm SST, but SST were not much warmer than during MIS 33. Interglacial to glacial transitions experienced 2 to 3 stadial/interstadial cycles, just like their mid-to-late Pleistocene counterparts. Glacial MIS 30 and 32 recorded periods of extremely cold (< 12°C) SST that in their climatic impact were comparable to the Heinrich events of the mid and late Pleistocene. Glacial MIS 34, on the other hand, was a relative warm glacial period off southern Portugal. Overall, surface water and MOW conditions at Site U1387 show strong congruence with Mediterranean climate, whereas millennial-scale variations are closely linked to North Atlantic circulation changes.

The early Eocene global warming events: Implications from ODP Leg 208 Walvis Ridge depth transect

Five sections drilled in multiple holes over a depth transect of more than 2200 m at the Walvis Ridge (SE Atlantic) during Ocean Drilling Program (ODP) Leg 208 resulted in the first complete early Paleogene deep-sea record. Here we present high-resolution stratigraphic records spanning a ~4.3 million yearlong interval of the late Paleocene to early Eocene. This interval includes the Paleocene-Eocene thermal maximum (PETM) as well as the Eocene thermal maximum (ETM) 2 event. A detailed chronology was developed with nondestructive X-ray fluorescence (XRF) core scanning records and shipboard color data. These records were used to refine the shipboard-derived spliced composite depth for each site and with a record from ODP Site 1051 were then used to establish a continuous time series over this interval. Extensive spectral analysis reveals that the early Paleogene sedimentary cyclicity is dominated by precession modulated by the short (100 kyr) and long (405 kyr) eccentricity cycles. Counting of precession-related cycles at multiple sites results in revised estimates for the duration of magnetochrons C24r and C25n. Direct comparison between the amplitude modulation of the precession component derived from XRF data and recent models of Earth's orbital eccentricity suggests that the onset of the PETM and ETM2 are related to a 100-kyr eccentricity maximum. Both events are approximately a quarter of a period offset from a maximum in the 405-kyr eccentricity cycle, with the major difference that the PETM is lagging and ETM2 is leading a 405-kyr eccentricity maximum. Absolute age estimates for the PETM, ETM2, and the magnetochron boundaries that are consistent with recalibrated radiometric ages and recent models of Earth's orbital eccentricity cannot be precisely determined at present because of too large uncertainties in these methods. Nevertheless, we provide two possible tuning options, which demonstrate the potential for the development of a cyclostratigraphic framework based on the stable 405-kyr eccentricity cycle for the entire Paleogene.

Sedimentation rate, bulk density, accumulation rate of terrigenous material and grain sizes of ODP Leg 167 sites

The terrigenous mineral fraction of sediments recovered by drilling during Ocean Drilling Program Leg 167 at Sites 1018 and 1020 is used to evaluate changes in the source and transport of fine-grained terrigenous sediment and its relation to regional climates and the paleoceanographic evolution of the California Current system during the late Pleistocene. Preliminary time scales developed by correlation of oxygen isotope stratigraphies with the global SPECMAP record show average linear sedimentation rates in excess of 100 m/m.y., which provide an opportunity for high-resolution studies of terrigenous flux, grain size, and mineralogy. The mass flux of terrigenous minerals at Site 1018 varies from 5 to 30 g/(cm**2 x k.y.) and displays a general trend toward increased flux during glacials. The terrigenous record at Site 1020 shows a similar pattern of increased glacial input, but overall accumulation rates are significantly lower. Spectral analysis demonstrates that most of this variability is concentrated in frequency bands related to orbital cycles of eccentricity, tilt, and precession. Detailed grain-size analysis performed on the isolated terrigenous mineral fraction shows that sediments from Site 1018 are associated with higher energy transport and depositional regimes than those found at Site 1020. Grain-size data are remarkably uniform throughout the last 500 k.y., with no discernible difference observed between glacial and interglacial size distributions within each site. X-ray diffraction analysis of the <2-µm clay component suggests that the deposition of minerals found at Site 1020 is consistent with transport from a southern source during intervals of increased terrigenous input.

(Table 1) Sapropel analysis from ODP Sites 160-964, 160-966, 160-967 and 160-969

Seventeen eastern Mediterranean Pliocene sapropels from ODP Sites 964, 966, 967 and 969, some of which are coeval, have been analysed for their geochemistry. The sapropels are characterized by very high organic carbon contents (up to 30%) which are reported to be the result of both increased productivity and improved preservation. Although the organic matter in the sapropels is mainly of marine origin, the d13Corg values and C/N ratios appear "terrestrial". This is the result of anaerobic organic matter degradation which preferentially removed nitrogen- and 13C-rich organic components. A comparison with Ti/Al profiles, which mimic the precession index, and a calculation of organic carbon accumulation rates indicate that sedimentation rates were at most 30% lower or at most 50% higher during sapropel formation. Thus, sapropel formation lasted from between 2000 and 10,000 years at Site 964 to between 4500 and 12,000 years at Site 967. A synthesis of new data and a comparison with existing models indicates that productivity, which increased due to extra nutrients supplied as a result of winter mixing and as a result of enhanced input by the Nile, was the driving mechanism behind sapropel formation. The resulting sapropel formation was simultaneous at different depths, but lasted longer in the part of the basin closest to the Nile.