Nd isotope data for ODP Leg 208 holes

The flow of deep-water masses is a key component of heat transport in the modern climate system, yet the role of deep-ocean heat transport during periods of extreme warmth is poorly understood. The present mode of meridional overturning circulation is characterized by deep-water formation in both the North Atlantic and the Southern Ocean. However, a different mode of meridional overturning circulation operated during the extreme greenhouse warmth of the early Cenozoic, during which time the Southern Ocean was the dominant region of deep-water formation. The combination of general global cooling and tectonic evolution of the Atlantic basins over the past ~55 m.y. ultimately led to the development of a mode of overturning circulation characterized by both Southern Ocean and North Atlantic deep-water sources. The change in deep-water circulation mode may, in turn, have affected global climate; however, unraveling the causes and consequences of this transition requires a better understanding of the timing of the transition. New Nd isotope data from the southeastern Atlantic Ocean indicate that the initial transition to a bipolar mode of deep-water circulation occurred in the early Oligocene, ca. 33 Ma. The likely cause of significant deep-water production in the North Atlantic was tectonic deepening of the sill separating the Greenland-Norwegian Sea from the North Atlantic.

Stable oxygen isotope data of Globorotalia inflata and radiocarbon dates for sediment cores GeoB6211-2, GeoB13862-1, and GeoB6308-3

The Southern Westerly Winds (SWW) exert a crucial influence over the world ocean and climate. Nevertheless, a comprehensive understanding of the Holocene temporal and spatial evolution of the SWW remains a significant challenge due to the sparsity of high-resolution marine archives and appropriate SWW proxies. Here, we present a north-south transect of high-resolution planktonic foraminiferal oxygen isotope records from the western South Atlantic. Our proxy records reveal Holocene migrations of the Brazil- Malvinas Confluence (BMC), a highly sensitive feature for changes in the position and strength of the northern portion of the SWW. Through the tight coupling of the BMC position to the large-scale wind field, the records allow a quantitative reconstruction of Holocene latitudinal displacements of the SWW across the South Atlantic. Our data reveal a gradual poleward movement of the SWW by about 1-1.5° from the early to the mid-Holocene. Afterwards variability in the SWW is dominated by millennial-scale displacements in the order of 1° in latitude with no recognizable longer-term trend. These findings are confronted with results from a state-of-the-art transient Holocene climate simulation using a comprehensive coupled atmosphere-ocean general circulation model. Proxy-inferred and modeled SWW shifts compare qualitatively, but the model underestimates both orbitally forced multi-millennial and internal millennial SWW variability by almost an order of magnitude. The underestimated natural variability implies a substantial uncertainty in model projections of future SWW shifts.

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.

Nitrogen isotope variations along the northeast Pacific margin

Glacial-interglacial changes in sedimentary d15N over the last 120 kyr display a remarkably similar pattern in timing and amplitude in core records extending from the denitrification zone in the eastern tropical North Pacific (ETNP), where subsurface denitrification is active, to the Oregon margin, where no denitrification occurs today. Low d15N values (4-6 per mil) generally characterize glacial stages 2 and 4, and higher d15N values (7-10 per mil) are representative of the Holocene, millennial-scale periods within stage 3, and stage 5. The inferred synchroneity of d15N variations along the entire margin implies that the nitrate isotopic signal produced in the oxygen-poor subsurface waters in the ETNP is rapidly advected northward and recorded at sites far beyond the boundaries of the modern denitrification zone. Similar to d15N, primary production indicators (percent Corg, Ba/Al, and percent opal) show glacial-interglacial as well as millennial-scale variations along the NE Pacific margin, with higher primary production during warm periods. However, the relative phasing between d15N and paleoproduction tracers within individual records changes latitudinally. Whereas d15N and primary production vary approximately synchronously in the midlatitudes, production lags d15N in the ETNP by several kiloyears. This lag calls for a new understanding of the processes driving denitrification in the ETNP. We suggest that oxygen input by the Equatorial Undercurrent as well as local organic matter flux controls denitrification rates in the ETNP. Moreover, the differences in relative timing point to a time-transgressive development of upwelling-favorable winds along the NE Pacific margin after the last glaciation, with those in the north developing several kiloyears earlier.

Paleogene and early Neogene stable isotope record of foraminifera from ODP Leg 119 samples

Oxygen and carbon isotopic records of monogeneric and monospecific benthic and planktonic foraminifer samples from Sites 744 and 738 drilled on the southern end of the Kerguelen Plateau during ODP Leg 119 reveal the evolution of polar Indian Ocean water masses from the early Paleocene to the middle Miocene. Results from Site 738 are from sediments of early Paleocene to late Eocene age and those from Site 744 are late Eocene to middle Miocene. They suggest that intermediate waters at this location did not originate in the high latitudes during the early Eocene. Surface and near-surface waters cooled gradually after the maximum warming at 56 Ma, when surface waters were about 18°C. Intermediate waters cooled after 52 Ma. The highest temperatures (lowest d18O values) of the Cenozoic occurred from 56 to 52 Ma. The records of equatorial Pacific Site 577 and Weddell Sea Site 690 resemble that of the polar Indian Ocean in this interval. The well-documented d13C excursions toward positive values in the late Paleocene and negative values in the early Eocene are represented by foraminifers increases of 1.5 per mil and following decreases of about 3 per mil. Most of the cooling in the Paleogene occurred in the middle and late Eocene. A 2°C decrease of surface water at about 38.4 Ma heralded the beginning of extensive glacial conditions in Antarctica in the early Oligocene. At Site 744, the global d18O shift just above the Eocene/Oligocene boundary is 1.15 per mil, and occurred gradually in sediments dated at 36.5-35.9 Ma. Ice-rafted debris was deposited beginning at 36.1 Ma for about the next 2 m.y. This simultaneous occurrence of the global d18O shift with ice-rafted debris is evidence for early Oligocene glaciation in East Antarctica. Moreover, early and late Oligocene Cibicidoides d18O values between 2 and 2.2 per mil indicate intermediate water cooling and a small ice-volume effect. Production of cold dense bottom water in Antarctica was intensified with continental cooling and glaciation in the early Oligocene. Comparison of Oligocene and early Miocene isotopic data from high-latitude and low-latitude deepsea sites indicates that there were probably at least two sources of bottom waters at this time.

Age determination of sediment cores from the Skagerrak

Paleoceanographic and paleoenvironmental interpretations based on foraminifera, sedimentary data, radiocarbon dates, and stable isotope measurements were derived from two sections in the Skagerrak: a 115-m-thick Holocene marine section drilled onshore at Skagen near the northernmost tip of Jutland, Denmark, and a 9-m piston core from the Skagerrak, north of Skagen. The foraminiferal data show that arctic-subarctic environments in the deep Skagerrak-Kattegat area were succeeded by boreal conditions at 9.6 ka. This was a result of northward migration of the Atlantic polar front and inflow of warm Atlantic water into the area through the Norwegian Channel. A gradual warming of the water masses after 9.6 ka is indicated by the data. Rare foraminifera and high sedimentation rates are found between approximately 8.6 ka and 7.6 ka at both core locations. The modern foraminiferal assemblages of the area were fully established at 7.6 ka indicating that the modern circulation pattern in the Skagerrak-Kattegat after the opening of the English Channel and the Danish Straits was not established before this date. At 5.5 ka a sudden change to coarser sediments (higher-energy environments) and the appearance of the foraminifer Eoeponidella laesoeensis is recorded in the Skagen core. This indicates a rapid change in the hydrography reflecting altered meteorological and hydrographic conditions in the Skagerrak-Kattegat, including a strengthening of the Jutland Current and increased inflow of North Sea water into the Kattegat. The event is interpreted as a response to cooling at the end of the Holocene climatic optimum in late Atlantic time and possibly reflects a rapid cooling event of North Atlantic surface water masses.

Miocene foraminiferal, bolboformid, dinoflagellate and diatom age calibrations for the northeastern North Atlantic

This study presents a new Miocene biostratigraphic synthesis for the high-latitude northeastern North Atlantic region. Via correlations to the bio-magnetostratigraphy and oxygen isotope records of Ocean Drilling Program and Deep Sea Drilling Project Sites, the ages of shallower North Sea deposits have been better constrained. The result has been an improved precision and documentation of the age designations of the existing North Sea foraminiferal zonal boundaries of King (1989) and Gradstein and Bäckström (1996). All calibrations have been updated to the Astronomically Tuned Neogene Time Scale (ATNTS) of Lourens et al. (2004). This improved Miocene biozonation has been achieved through: the updating of age calibrations for key microfossil bioevents, identification of new events, and integration of new biostratigraphic data from a foraminiferal analysis of commercial wells in the North Sea and Norwegian Sea. The new zonation has been successfully applied to two commercial wells and an onshore research borehole. At these high latitudes, where standard zonal markers are often absent, integration of microfossil groups significantly improves temporal resolution. The new zonation comprises 11 Nordic Miocene (NM) Zones with an average duration of 1 to 2 million years. This multi-group combination of a total of 92 bioevents (70 foraminifers and bolboformids; 16 dinoflagellate cysts and acritarchs; 6 marine diatoms) facilitates zonal identification throughout the Nordic Atlantic region. With the highest proportion of events being of calcareous walled microfossils, this zonation is primarily suited to micropaleontologists. A correlation of this Miocene biostratigraphy with a re-calibrated oxygen isotope record for DSDP Site 608 suggests a strong correlation between Miocene planktonic microfossil turnover rates and the inferred paleoclimatic trends. Benthic foraminifera zonal boundaries appear to often coincide with Miocene global sequence boundaries. The biostratigraphic record is punctuated by four main stratigraphic hiati which show variation in their geographic and temporal extent. These are related to the following regional unconformities: basal Neogene, Lower/Middle Miocene ("mid-Miocene unconformity"), basal Upper Miocene and basal Messinian unconformities. Further coring of Neogene sections in the North Sea and Norwegian Sea may better constrain their extent and their effect on the biostratigraphic record.

Uranium contents and isotope ratios of uranium and strontium in sediments, leachates and pore fluids from ODP Hole 162-984A

Bulk dissolution rates for sediment from ODP Site 984A in the North Atlantic are determined using the 234U/238U activity ratios of pore water, bulk sediment, and leachates. Site 984A is one of only several sites where closely spaced pore water samples were obtained from the upper 60 meters of the core; the sedimentation rate is high (11-15 cm/ka), hence the sediments in the upper 60 meters are less than 500 ka old. The sediment is clayey silt and composed mostly of detritus derived from Iceland with a significant component of biogenic carbonate (up to 30%). The pore water 234U/238U activity ratios are higher than seawater values, in the range of 1.2 to 1.6, while the bulk sediment 234U/238U activity ratios are close to 1.0. The 234U/238U of the pore water reflects a balance between the mineral dissolution rate and the supply rate of excess 234U to the pore fluid by a-recoil injection of 234Th. The fraction of 238U decays that result in a-recoil injection of 234U to pore fluid is estimated to be 0.10 to 0.20 based on the 234U/238U of insoluble residue fractions. The calculated bulk dissolution rates, in units of g/g/yr are in the range of 0.0000004 to 0.000002 1/yr. There is significant down-hole variability in pore water 234U/238U activity ratios (and hence dissolution rates) on a scale of ca. 10 m. The inferred bulk dissolution rate constants are 100 to 1000 times slower than laboratory-determined rates, 100 times faster than rates inferred for older sediments based on Sr isotopes, and similar to weathering rates determined for terrestrial soils of similar age. The results of this study suggest that U isotopes can be used to measure in situ dissolution rates in fine-grained clastic materials. The rate estimates for sediments from ODP Site 984 confirm the strong dependence of reactivity on the age of the solid material: the bulk dissolution rate (R_d) of soils and deep-sea sediments can be approximately described by the expression R_d ~ 0.1 1/age for ages spanning 1000 to 500,000,000 yr. The age of the material, which encompasses the grain size, surface area, and other chemical factors that contribute to the rate of dissolution, appears to be a much stronger determinant of dissolution rate than any single physical or chemical property of the system.

Volcanic rocks from the submarine Vityaz Ridge: their age, chemical and isotopic compositions

Results of geological studies at the submarine Vityaz Ridge carried out during cruises 37 and 41 of R/V Akademik Lavrent'ev in 2005 and 2006 are reported. The studied area is located at an near-island trench of the slope in the central part of the Kuril Island arc. Morphologically it consists of two parts: an inner volcanic arc represented by the Great Kuril Range and an outer arc corresponding to the submarine Vityaz Ridge. Diverse rocks composing the basement and the sedimentary cover of the ridge were recovered by dredging. Based on K-Ar dating and geochemistry, volcanics were divided into Paleocene, Eocene, late Oligocene, and Pliocene-Pleistocene complexes. Each of the complexes reflects a tectonomagmatic stage in the ridge evolution. Geochemical and isotope data on the volcanics indicate contribution of ancient crustal material in the magma source and, correspondingly, formation of this structure on the continental basement. Two-stage model ages (TDM2) vary in a wide range from zero values in mafic rocks to 0.77 Ga in felsic varieties, pointing to presence of Precambrian protolith in the source of the felsic rocks of the Vityaz Ridge. The Pliocene-Pleistocene volcanics are classed with tholeiitic, calc-alkaline, and subalkaline series, which differ in alkali contents and REE fractionation. Values of (La/Sm)_n and (La/Yb)_n ratios vary from 0.74 and 0.84 in the tholeiitic varieties to 1.19 and 1.44 in the calc-alkaline and 2.32 and 3.73 in the subalkaline rocks. All three varieties occur within the same volcanic edifices and formed during differentiation of magmatic melts that were channeled along fault zones from the mantle source slightly enriched in crustal component.

Stable isotope composition and ice rafted debris in early Late Pliocene sediments from Maud Rise, South Atlantic

We have conducted an integrated study of ice-rafted debris (IRD) and oxygen isotopes (measured on Cibicides, Globigerina bulloides, and Neogloboquadrina pachyderma, using identical samples). We used samples from the early Late Pliocene Gauss Chron from ODP Site 114-704 on the Meteor Rise in the subantarctic South Atlantic. During the early Gauss Chron, the oxygen isotopic ratios are generally up to 0.5?-0.6? less than their respective Holocene values. The lowest values in this record can accommodate a warming of about 2.5°C or a sea-level rise of about 50 m, but not both, and probably result from some warming and a small reduction in global ice volume. Starting with isotope stage MG2 [ 3.23 Ma on the Berggren et al. ( 1985) time scale; 3.38 on the Shackleton et al. ( 1995b) time scale] oxygen-isotopic values generally increase (and oscillate about a Holocene mean). The first significant IRD appears at the same time. There is a subsequent increase in IRD amounts upsection. In order to reach the site, this material must have been transported by large, tabular icebergs derived from Antarctic ice shelves or ice tongues, similar to occasional, large modern icebergs. This combined record suggests strongly that the Antarctic ice sheet was essentially intact; some warming at the drill site is indicated, but not a major reduction in ice-volume on Antarctica.

Age determination and faunal composition of sediments off northern Greenland

We present a record encompassing marine isotope stages 7-1 from a hitherto unexplored and heavily ice-covered area of the Arctic Ocean, the Lomonosov Ridge off the northern Greenland-Canada continental margin, using nannofossil and benthic foraminifera stratigraphy. Planktic foraminifera assemblages are used as a key paleoceanographic proxy, and a surprisingly large variability is found for an interior Arctic Ocean site. Abundant small (63-125 µm) subpolar Turborotalita quinqueloba occur in two sections, possibly representing substages 5e (last interglacial) and 5a (warm interstadial). However, the present-day circulation pattern and the very distant location of high productive regions cannot explain such high abundances of subpolar specimens in the interior, perennially sea ice-covered Arctic Ocean. Hence our proxy record indicates that last interglacial sea ice concentrations were reduced off some areas of northern Greenland-Canada. Whether this was part of a larger regional pattern or it represents the influence of polynya areas with locally increased productivity remains to be solved. With respect to glacial conditions, increased ice-rafted debris (IRD) deposition in the area appears to be associated with glacial stages 6, 4, and late 3. Stage 2 sediments (including the Last Glacial Maximum) are condensed with a sparse IRD content only.

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.

Triennial skeletal stable oxygen isotope ratios of Diploria labyrinthiformis

The reconstruction of the climatic history during the past several hundred years requires a sufficient geographical coverage of combined climate proxy series. Especially in order to identify causal connections between the atmosphere and the ocean, inclusion of marine records into composite climate time series is of fundamental importance. We present two skeletal delta18O chronologies of coral skeletons of Diploria labyrinthiformis from Bermuda fore-reef sites covering periods in the nineteenth and twentieth centuries and compare them with instrumental temperature data. Both time series are demonstrated to display sea-surface temperature (SST) variability on inter-annual to decadal time scales. On the basis of a specific modern delta18O vs instrumental SST calibration we reconstruct a time series of SST anomalies between AD 1350 and 1630 covering periods during the Little Ice Age. The application of the coral delta18O vs temperature relationship leads to estimates of past SST variability which are comparable to the magnitude of modern variations. Parallel to delta18O chronologies we present time series of skeletal bulk density. Coral delta18O and skeletal density reveal a strong similarity during Little Ice Age, confirming the reliability of both proxy climate indicators. The past coral records, presented in this study, share features with a previously published climate proxy record from Bermuda and a composite time series of reconstructed Northern Hemisphere summer temperatures. The coral proxy data presented here represent a valuable contribution to elucidate northern Atlantic subtropical climate variation during the past several centuries.

Stable oxygen isotope, UK'37, productivity, CaCO3 and XRF records from sediment core MD06-3075

Proxy records of hydrologic variability in the West Pacific Warm Pool (WPWP) have revealed wide-scale changes in past convective activity in response to orbital and sub-orbital climate forcings. However, attributing proxy responses to regional changes in WPWP hydrology versus local variations in precipitation requires independent records linking the terrestrial and marine realms. We present high-resolution stable isotope, UK'37 sea-surface temperature, X-ray fluorescence (XRF) core scanning and coccolithophore-derived paleoproductivity records covering the past 120 ka from International Marine Global Change (IMAGES) Program Core MD06-3075 (6°29' N, 125°50' E, water depth 1878 m), situated in the Davao Gulf on the southern side of Mindanao. XRF-derived log(Fe/Ca) records provide a robust proxy for runoff-driven sedimentary discharge from Mindanao, whilst past changes in local productivity are associated with variable freshwater runoff and stratification of the surface layer. Significant precessional-scale variability in sedimentary discharge occurred during Marine Isotope Stage (MIS) 5, with peaks in discharge contemporaneous with Northern Hemisphere summer insolation minima. We attribute these changes to the latitudinal migration of the Intertropical Convergence Zone (ITCZ) over the WPWP together with variability in the strength of the Walker circulation acting on precessional timescales. Between 60 and 15 ka sedimentary discharge at Mindanao was muted, displaying little orbital- or millennial-scale variability, likely in response to weakened precessional insolation forcing and lower sea level driving increased subsidence of air masses over the exposed Sunda Shelf. These results highlight the high degree of local variability in the precipitation response to past climate changes in the WPWP.

Middle to late Miocene oxygen isotope stratigraphy of ODP site 1085 (SE Atlantic): new constrains on Miocene climate variability and sea-level fluctuations

The middle Miocene delta18O increase represents a fundamental change in earth's climate system due to a major expansion and permanent establishment of the East Antarctic Ice Sheet accompanied by some effect of deepwater cooling. The long-term cooling trend in the middle to late Miocene was superimposed by several punctuated periods of glaciations (Mi-Events) characterized by oxygen isotopic shifts that have been related to the waxing and waning of the Antarctic ice-sheet and bottom water cooling. Here, we present a high-resolution benthic stable oxygen isotope record from ODP Site 1085 located at the southwestern African continental margin that provides a detailed chronology for the middle to late Miocene (13.9-7.3 Ma) climate transition in the eastern South Atlantic. A composite Fe intensity record obtained by XRF core scanning ODP Sites 1085 and 1087 was used to construct an astronomically calibrated chronology based on orbital tuning. The oxygen isotope data exhibit four distinct delta18O excursions, which have astronomical ages of 13.8, 13.2, 11.7, and 10.4 Ma and correspond to the Mi3, Mi4, Mi5, and Mi6 events. A global climate record was extracted from the oxygen isotopic composition. Both long- and short-term variabilities in the climate record are discussed in terms of sea-level and deep-water temperature changes. The oxygen isotope data support a causal link between sequence boundaries traced from the shelf and glacioeustatic changes due to ice-sheet growth. Spectral analysis of the benthic delta18O record shows strong power in the 400-kyr and 100-kyr bands documenting a paleoceanographic response to eccentricity-modulated variations in precession. A spectral peak around 180-kyr might be related to the asymmetry of the obliquity cycle indicating that the response of the dominantly unipolar Antarctic ice-sheet to obliquityinduced variations probably controlled the middle to late Miocene climate system. Maxima in the delta18O record, interpreted as glacial periods, correspond to minima in 100-kyr eccentricity cycle and minima in the 174-kyr obliquity modulation. Strong middle to late Miocene glacial events are associated with 400-kyr eccentricity minima and obliquity modulation minima. Thus, fluctuations in the amplitude of obliquity and eccentricity seem to be the driving force for the middle to late Miocene climate variability.