Benthic foraminiferal biofacies and stable isotopic record of St. Stephens Quarry, Alabama

We integrate upper Eocene-lower Oligocene lithostratigraphic, magnetostratigraphic, biostratigraphic, stable isotopic, benthic foraminiferal faunal, downhole log, and sequence stratigraphic studies from the Alabama St. Stephens Quarry (SSQ) core hole, linking global ice volume, sea level, and temperature changes through the greenhouse to icehouse transition of the Cenozoic. We show that the SSQ succession is dissected by hiatuses associated with sequence boundaries. Three previously reported sequence boundaries are well dated here: North Twistwood Creek-Cocoa (35.4-35.9 Ma), Mint Spring-Red Bluff (33.0 Ma), and Bucatunna-Chickasawhay (the mid-Oligocene fall, ca. 30.2 Ma). In addition, we document three previously undetected or controversial sequences: mid-Pachuta (33.9-35.0 Ma), Shubuta-Bumpnose (lowermost Oligocene, ca. 33.6 Ma), and Byram-Glendon (30.5-31.7 Ma). An ~0.9 per mil d18O increase in the SSQ core hole is correlated to the global earliest Oligocene (Oi1) event using magnetobiostratigraphy; this increase is associated with the Shubuta-Bumpnose contact, an erosional surface, and a biofacies shift in the core hole, providing a first-order correlation between ice growth and a sequence boundary that indicates a sea-level fall. The d18O increase is associated with a eustatic fall of ~55 m, indicating that ~0.4 per mil of the increase at Oi1 time was due to temperature. Maximum d18O values of Oi1 occur above the sequence boundary, requiring that deposition resumed during the lowest eustatic lowstand. A precursor d18O increase of 0.5 per mil (33.8 Ma, midchron C13r) at SSQ correlates with a 0.5 per mil increase in the deep Pacific Ocean; the lack of evidence for a sea-level change with the precursor suggests that this was primarily a cooling event, not an ice-volume event. Eocene-Oligocene shelf water temperatures of ~17-19 °C at SSQ are similar to modern values for 100 m water depth in this region. Our study establishes the relationships among ice volume, d18O, and sequences: a latest Eocene cooling event was followed by an earliest Oligocene ice volume and cooling event that lowered sea level and formed a sequence boundary during the early stages of eustatic fall.

Age models for sediment cores in the eastern tropical Pacific

We present new high-resolution N isotope records from the Gulf of Tehuantepec and the Nicaragua Basin spanning the last 50-70 ka. The Tehuantepec site is situated within the core of the north subtropical denitrification zone while the Nicaragua site is at the southern boundary. The d15N record from Nicaragua shows an 'Antarctic' timing similar to denitrification changes observed off Peru-Chile but is radically different from the northern records. We attribute this to the leakage of isotopically heavy nitrate from the South Pacific oxygen minimum zone (OMZ) into the Nicaragua Basin. The Nicaragua record leads the other eastern tropical North Pacific (ETNP) records by about 1000 years because denitrification peaks in the eastern tropical South Pacific (ETSP) before denitrification starts to increase in the Northern Hemisphere OMZ, i.e., during warming episodes in Antarctica. We find that the influence of the heavy nitrate leakage from the ETSP is still noticeable, although attenuated, in the Gulf of Tehuantepec record, particularly at the end of the Heinrich events, and tends to alter the recording of millennial timescale denitrification changes in the ETNP. This implies (1) that sedimentary d15N records from the southern parts of the ETNP cannot be used straightforwardly as a proxy for local denitrification and (2) that denitrification history in the ETNP, like in the Arabian Sea, is synchronous with Greenland temperature changes. These observations reinforce the conclusion that on millennial timescales during the last ice age, denitrification in the ETNP is strongly influenced by climatic variations that originated in the high-latitude North Atlantic region, while commensurate changes in Southern Ocean hydrography more directly, and slightly earlier, affected oxygen concentrations in the ETSP. Furthermore, the d15N records imply ongoing physical communication across the equator in the shallow subsurface continuously over the last 50-70 ka.

Sea surface temperature reconstruction of ODP Site 175-1087

Successful application of the alkenone palaeothermometer, the UK'37 index, relies upon the assumption that fossil alkenone synthesisers responded to growth-temperature changes in a similar manner to the modern producers, chiefly the coccolithophores Emiliania huxleyi and Gephyrocapsa oceanica. We compare coccolith and UK'37 data from ODP Site 1087 in the south-east Atlantic between 1500 and 500 ka, and show that evolutionary events and changes in species dominance within the coccolithophore populations had little impact on the UK'37 record. The relative abundances of the C37 and C38 alkenones also closely resembled those found in modern populations, and suggest a similar temperature sensitivity of UK'37 during the early and mid-Pleistocene to that found at present. These results support the application of the UK'37 index to reconstruct sea-surface temperatures (SSTs) throughout the Quaternary. The UK'37 record at ODP Site 1087 contains an SST signal that documents the emergence of the 100-kyr cycles that characterise the late Quaternary ice volume records. This is preceded by significant cooling at ODP Site 1087, marked by a negative shift in SSTs and a positive shift in the planktonic delta18O some 250-kyr earlier, at ca 1150-1000 ka. This results in a permanent fall in average SSTs of around 1.5 °C. The predicted increase in aridity onshore as a result of this cooling can be identified in a number of published records from southern Africa, and may have played a role in some important evolutionary events of the mid-Pleistocene.

Age models and sea-surface paleotemperature reconstruction of sediment cores from the eastern equatorial Atlantic

Based on the faunal record of planktonic foraminifers in three long gravity sediment cores from the eastern equatorial Atlantic, the sea-surface temperature history ove the last 750,000 years was studied at a resolution of 3,000 to 10,000 years. Detailed oxygen-isotope and paleomagnetic stratigraphy helped to identify the following major faunal events: Globorotaloides hexagonus and Globorotalia tumida flexuosa became extinct in the eastern tropical Atlantic at the isotope stage 4/5 boundary, now dated at 68,000 years B.P. The persistent occurrence of the pink variety of Globigerinoides ruber started during the late stage 12 at 410,000 years B.P. CARTUNE-age. This datum may provide an easily detectible faunal stratigraphic marker for the mid-Brunhes Chron. The updated scheme of the Ericson zones helped the recognition of a hiatus at the northwestern slope of the Sierra Leone Basin covering oxygen-isotope stages 10 to 12. Classifying the planktonic foraminifer counts into six faunal assemblages, according to the factor analysis derived model of Pflaumann (1985), the tropical and the tropical-upwelling communities account for 57 % at Site 16415, and 86 % at Site 13519, respectively of the variance of the faunal record. A largely continuous paleotemperature record for both winter and summer seasons was obtained from the top of the Sierra Leone Rise with the winter temperatures ranging between 20 and 25 °C, and the summer ones between 24 and 30 °C. The record of cores from greater water depths is frequently interrupted by samples with no-analogue faunal communities and/or poor preservation. Based on the seasonality signal, during cold periods the termal equator shifted to a geographically mnore asymmetrical northern position. Dissolution altering the faunal communities becomes stronger with greater water depth, the estimated mean minimum loss of specimens increases from 70 % to 80 % between 2,860 and 3,850 water depth although some species will be more susceptible than others. Enhanced dissolution occured during stage 4 but also during cold phases in the warm stage 7 and 9. Correlations between the Foraminiferal Dissolution Index and the estimated sea-surface temperatures are significant. Foraminiferal flux rates, negatively correlated to the flux rates of organic carbon and of diatoms, may be a result of enhanced dissolution during cold stages, destroying still more of the faunal signal than indicated by the calculated minimum loss. The fluctuations of the oxygen-isotope curves and the hibernal sea-surfave temperatures are fairly coherent. During warm oxygen-isotope stages the temperature maxima lag often by 5 to 15 ka behind the respective sotope minima. During cold stages, sea-surface temperature changes are partly out of phase and contain additional fluctuations.

Stable isotope record of foraminifera from sediment core EW9209-1JPC

Stable isotopic measurements of G. sacculifer and C. wuellerstorfi in a core from the western equatorial Atlantic imply that there are parallel, suborbital oscillations in surface water hydrography and deep water circulation occurring during oxygen isotope stages 2 and 3. Low values of G. sacculifer delta18O accompany high values of C. wuellerstorfi delta13C, linking warmer sea surface temperatures (SSTs) in the tropics with increased production of lower North Atlantic Deep Water (NADW). The amplitude of the delta18O oscillations is 0.6 per mil (or 2°-3°C), which is superimposed on a glacial/interglacial amplitude of about 2.1per mil. Using the G. sacculifer delta18O data, we calculate that surface waters were colder during stage 2 than calculated by CLIMAP [1976, 1981]. The longer-period (>2 kyr) oscillations in air temperature recorded in the Greenland and Antarctic ice cores appear to correlate with oscillations in sea surface temperature in the equatorial Atlantic. The magnitude of these oscillations in tropical SST is too large to have resulted from changes in meridional heat transport caused by the global conveyor alone. The apparent synchroneity of equatorial SST and polar air temperature changes, as well as the amplitude of the SST changes at the equator, are consistent with the climate effects expected from changes in the atmosphere's greenhouse gas content (H2Ovapor, CO2, and CH4).

Sea surface temperature estimation of the South Tasman Rise

The Subtropical Front (STF) marking the northern boundary of the Southern Ocean has a steep gradient in sea surface temperature (SST) of approximately 4°C over 0.5° of latitude. Presently, in the region south of Tasmania, the STF lies nominally at 47°S in the summer and 45°S in the winter. We present here SST reconstructions in a latitudinal transect of cores across the South Tasman Rise, southeast of Australia, during the late Quaternary. SST reconstructions are based on two paleotemperature proxies, alkenones and faunal assemblages, which are used to assess past changes in SST in spring and summer. The north-south alignment in core locations allows reconstruction of movement of the STF over the last 100 ka. Surface water temperatures during the last glaciation in this region were ~4°C colder than today. Additional temperature changes greater in magnitude than 4°C seen in individual cores can be attributed to changes in the water mass overlying the core site caused by the movement of the front across that location. During the penultimate interglacial, SST was ~2°C warmer and the STF was largely positioned south of 47°S. Movement of the STF to the north occurred during cool climate periods such as the last marine isotope stages 3 and 4. In the last glaciation, the front was at its farthest north position, becoming pinned against the Tasmanian landmass. It moved south by 4° latitude to 47°S in summer during the deglaciation but remained north of 45°S in spring throughout the early deglaciation. After 11 ka B.P. inferred invigoration of the East Australia Current appears to have pushed the STF seasonally south of the East Tasman Plateau, until after 6 ka B.P. when it achieved its present configuration.

Radiolarian assemblages of the late Neogene North Pacific

Radiolarian census and abundance data were collected from three deep-sea cores drilled by the Ocean Drilling Program Sites 884, 887 and 1151 to investigate patterns of ecologic changes in space and time during the last 16 million years for the mid-latitude to subarctic North Pacific. High concentrations of radiolarians occurred between 9.0 and 2.7 Ma. Radiolarian species richness was highest in the early middle Miocene at each site and gradually decreased up to about 7 Ma, coinciding with a well-established global cooling trend. A degree of overlap index calculated for radiolarian assemblages revealed 11 faunal change events, of which 8 corresponded to global cooling events and expansions of polar ice sheets. Three of the faunal change events were observed within the peak of radiolarian accumulation rate and were ascribed to changes in primary productivity in the North Pacific rather than global climatic changes. Our assemblage analyses revealed that north-south differentiation in radiolarian assemblages in the northwestern Pacific has existed since 16 Ma and became more distinct via major steps at 6.8 Ma and 2.7 Ma, coinciding with major glaciation events, and that east-west faunal contrasts in the subarctic region became obvious beginning at 11.7 Ma and changed to a different mode around 6.8 Ma. The observed east-west faunal differences possibly reflect east to west climate differences that were characterized by cooler temperatures in the east than the west during the late Miocene (11.7-6.8 Ma) and then by the opposite temperature trend (6.8 Ma-Recent). A severe glaciation at 2.7 Ma played a large role, particularly in temporal changes in radiolarian accumulation rate and assemblage composition.

Effects of elevated CO2 and temperature on an intertidal meiobenthic community

In the near future, the marine environment is likely to be subjected to simultaneous increases in temperature and decreased pH. The potential effects of these changes on intertidal, meiofaunal assemblages were investigated using a mesocosm experiment. Artificial Substrate Units containing meiofauna from the extreme low intertidal zone were exposed for 60 days to eight experimental treatments (four replicates for each treatment) comprising four pH levels: 8.0 (ambient control), 7.7 & 7.3 (predicted changes associated with ocean acidification), and 6.7 (CO2 point-source leakage from geological storage), crossed with two temperatures: 12 °C (ambient control) and 16 °C (predicted). Community structure, measured using major meiofauna taxa was significantly affected by pH and temperature. Copepods and copepodites showed the greatest decline in abundance in response to low pH and elevated temperature. Nematodes increased in abundance in response to low pH and temperature rise, possibly caused by decreased predation and competition for food owing to the declining macrofauna density. Nematode species composition changed significantly between the different treatments, and was affected by both seawater acidification and warming. Estimated nematode species diversity, species evenness, and the maturity index, were substantially lower at 16 °C, whereas trophic diversity was slightly higher at 16 °C except at pH 6.7. This study has demonstrated that the combination of elevated levels of CO2 and ocean warming may have substantial effects on structural and functional characteristics of meiofaunal and nematode communities, and that single stressor experiments are unlikely to encompass the complexity of abiotic and biotic interactions. At the same time, ecological interactions may lead to complex community responses to pH and temperature changes in the interstitial environment.

Sea surface temperature reconstruction of sediment core GeoB6518-1

Past hydrological changes in Africa have been linked to various climatic processes, depending on region and timescale. Long-term precipitation changes in the regions of northern and southern Africa influenced by the monsoons are thought to have been governed by precessional variations in summer insolation (Kutzbach and Liu, 1997, doi:10.1126/science.278.5337.440; Partridge et al., 1997, doi:10.1016/S0277-3791(97)00005-X). Conversely, short-term precipitation changes in the northern African tropics have been linked to North Atlantic sea surface temperature anomalies, affecting the northward extension of the Intertropical Convergence Zone and its associated rainbelt (Hastenrath, 1990, doi:10.1002/joc.3370100504, Street-Perrott and Perrott, 1990, doi:10.1038/343607a0). Our knowledge of large-scale hydrological changes in equatorial Africa and their forcing factors is, however, limited (Gasse, 2000, doi:10.1016/S0277-3791(99)00061-X). Here we analyse the isotopic composition of terrigenous plant lipids, extracted from a marine sediment core close to the Congo River mouth, in order to reconstruct past central African rainfall variations and compare this record to sea surface temperature changes in the South Atlantic Ocean. We find that central African precipitation during the past 20,000 years was mainly controlled by the difference in sea surface temperatures between the tropics and subtropics of the South Atlantic Ocean, whereas we find no evidence that changes in the position of the Intertropical Convergence Zone had a significant influence on the overall moisture availability in central Africa. We conclude that changes in ocean circulation, and hence sea surface temperature patterns, were important in modulating atmospheric moisture transport onto the central African continent.

High resolution benthic stable-isotope records from ODP Site 1263 and 1262

Recent studies have shown that the Early Eocene Climatic Optimum (EECO) was preceded by a series of short-lived global warming events, known as hyperthermals. Here we present high-resolution benthic stable carbon and oxygen isotope records from ODP Sites 1262 and 1263 (Walvis Ridge, SE Atlantic) between ~54 and ~52 million years ago, tightly constraining the character, timing, and magnitude of six prominent hyperthermal events. These events, which include Eocene Thermal Maximum (ETM) 2 and 3, are studied in relation to orbital forcing and long-term trends. Our findings reveal an almost linear relationship between d13C and d18O for all these hyperthermals, indicating that the eccentricity-paced co-variance between deep-sea temperature changes and extreme perturbations in the exogenic carbon pool persisted during these events towards the onset of the EECO, in accord with previous observations for the Paleocene Eocene Thermal Maximum (PETM) and ETM2. The covariance of d13C and d18O during H2 and I2, which are the second pulses of the "paired" hyperthermal events ETM2-H2 and I1-I2, deviates with respect to the other events. We hypothesize that this could relate to a relatively higher contribution of an isotopically heavier source of carbon, such as peat or permafrost, and/or to climate feedbacks/local changes in circulation. Finally, the d18O records of the two sites show a systematic offset with on average 0.2 per mil heavier values for the shallower Site 1263, which we link to a slightly heavier isotopic composition of the intermediate water mass reaching the northeastern flank of the Walvis Ridge compared to that of the deeper northwestern water mass at Site 1262.

Oliogocene to Pleistocene benthic foraminifera in ODP Leg 121 sites

Oligocene to Pleistocene bathyal benthic foraminifers at Broken Ridge (Site 754) and Ninetyeast Ridge (Site 756), eastern Indian Ocean, were investigated for then- stratigraphic distribution and their response to paleoceanographic changes. Q-mode factor analysis was applied to relative abundance data of the most abundant benthic foraminifers. At Site 754, seven varimax assemblages were recognized from the upper Oligocene to the Pleistocene: the Gyroidina orbicularis-Rectuvigerina striata Assemblage in the uppermost Oligocene; the Lenticulina spp. Assemblage in the upper Oligocene to lower Miocene, and in lower Miocene to lowermost middle Miocene; the Burseolina cf. pacifica-Cibicidoides mundulus Assemblage in the lower Miocene; the Planulina wuellerstorfi Assemblage in the upper middle Miocene; the Globocassidulina spp. Assemblage in the upper Miocene; the Gavelinopsis lobatulus-Uvigerina proboscidea Assemblage in the Pliocene; and the Ehrenbergina spp. Assemblage in the Pleistocene. The major faunal changes are complex, but exist between the Lenticulina spp. Assemblage and the P. wuellerstorfi Assemblage at ~13.8 Ma, and between the Ehrenbergina spp. Assemblage and the G. lobatulus Assemblage at ~5 Ma. The development of the P. wuellerstorfi and Globocassidulina spp. Assemblages after 13.8 Ma is correlated with the decrease in temperature of the intermediate waters of the ocean, in turn related to Antarctic glacial expansion. The faunal changes at ~5 Ma are related to the development of low oxygen intermediate water, formed in the presence of a strong thermocline. At Site 756, six varimax assemblages are distributed as follows: the Cibicidoides cf. mundulus-Oridorsalis umbonatus Assemblage in the lower Oligocene; the Epistominella umbonifera-Cibicidoides mundulus Assemblage from the upper Oligocene to the lower Miocene; the Cibicidoides mundulus-Burseolinapacifica Assemblage from lower Miocene to the lower middle Miocene; the Globocassidulina spp. Assemblage from the upper lower Miocene to the Pliocene; the Uvigerina proboscidea Assemblage in the upper Miocene and the Pliocene; and the Globocassidulina sp. D Assemblage in the Pliocene. The main faunal change at this site is between the E. umbonifera Assemblage and the Globocassidulina spp. Assemblage, at ~17.1 Ma. The timing of this faunal change is coeval with faunal changes in the North Atlantic and the Pacific. The change is related to a change in bottom water characteristics caused by an increased influence of carbonate corrosive water from the Antarctic source region, and a change in surface productivity. A low oxygen event at Site 756, which started at about 7.3 Ma, occurred about 2.3 m.y. before that at Site 754. The different response to global paleoceanographic changes is not yet explained, but may be due to the difference of marine topography and the degree of upwelling

Stable isotope analysis of Porites coral core from the Andaman Islands

North Atlantic climate variations are reflected in sedimentary records from the northern Indian Ocean in which two basins, the Arabian Sea and the Bay of Bengal, are strongly affected by the monsoon. Contrary to the Bay of Bengal the Arabian Sea plays an important role in the global marine nitrogen cycle. In its mid-water oxygen minimum zone (OMZ) bioavailable fixed nitrogen is reduced to nitrogen gas (NO3- - > N2), whereas oxygen concentrations are slightly above the threshold of nitrate reduction in the OMZ of the Bay of Bengal. A coral colony (Porites lutea) growing south of Port Blair on the Andaman Islands in the Bay of Bengal was studied for its response to changes in the monsoon system and its link to temperature changes in the North Atlantic Ocean, between 1975 and 2006. Its linear extension rates, d13C and d18O values measured within the coral skeleton reveal a strong seasonality, which seems to be caused by the monsoon-driven reversal of the surface ocean circulation. The sampling site appears to be influenced by low salinity Bay of Bengal Water during the NE monsoon (boreal winter) and by the high salinity Arabian Sea Water during the SW monsoon in summer. The high salinity Arabian Sea Water circulates along with the Summer Monsoon Current (S-MC) from the Arabia Sea into the Bay of Bengal. Decreasing d18O and reconstructed salinity values correlate to the increasing SSTs in the North Atlantic Ocean indicating a reduced influence of the S-MC at the sampling site in the course of northern hemispheric warming. During such periods oxygen-depletion became stronger in the OMZ of the Arabian Sea as indicated by the sedimentary records. A reduced propagation of oxygen-depleted high salinity Arabian Sea Water into the Bay of Bengal could be a mechanism maintaining oxygen concentration above the threshold of nitrate reduction in the OMZ of the Bay of Bengal in times of global warming.

Revised calibration of the MBT-CBT paleotemperature proxy

The MBT-CBT proxy for the reconstruction of paleotemperatures and past soil pH is based on the distribution of branched glycerol dialkyl glycerol tetraether (brGDGT) membrane lipids. The Methylation of Branched Tetraether (MBT) and the Cyclisation of Branched Tetraether (CBT) indices were developed to quantify these distributions, and significant empirical relations between these indices and annual mean air temperature (MAT) and/or soil pH were found in a large data set of soils. In this study, we extended this soil dataset to 278 globally distributed surface soils. Of these soils, 26% contains all nine brGDGTs, while in 63% of the soils the seven most common brGDGTs were detected, and the latter were selected for calibration purposes. This resulted in new transfer functions for the reconstruction of pH based on the CBT index: pH = 7.90-1.97 × CBT (r**2 = 0.70; RMSE = 0.8; n = 176), as well as for MAT based on the CBT index and methylation index based on the seven most abundant GDGTs (defined as MBT'): MAT = 0.81-5.67 × CBT + 31.0 × MBT' (r**2 = 0.59; RMSE = 5.0 °C; n = 176). The new transfer function for MAT has a substantially lower correlation coefficient than the original equation (r**2 = 0.77). To investigate possible improvement of the correlation, we used our extended global surface soil dataset to statistically derive the indices that best describe the relations of brGDGT composition with MAT and soil pH. These new indices, however, resulted in only a relatively minor increase in correlation coefficients, while they cannot be explained straightforwardly by physiological mechanisms. The large scatter in the calibration cannot be fully explained by local factors or by seasonality, but MAT for soils from arid regions are generally substantially (up to 20 °C) underestimated, suggesting that absolute brGDGT-based temperature records for these areas should be interpreted with caution. The applicability of the new MBT'-CBT calibration function was tested using previously published MBT-CBT-derived paleotemperature records covering the last deglaciation in Central Africa and East Asia, the Eocene-Oligocene boundary and the Paleocene-Eocene thermal maximum. The results show that trends remain similar in all records, but that absolute temperature estimates and the amplitude of temperature changes are lower for most records, and generally in better agreement with independent proxy data.

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.