Orbital-resolution Pliocene-Pleistocene simulations of CO2, sea level, global temperature and benthic d18O, and d11B-based proxy-CO2 data

During the past five million yrs, benthic d18O records indicate a large range of climates, from warmer than today during the Pliocene Warm Period to considerably colder during glacials. Antarctic ice cores have revealed Pleistocene glacial-interglacial CO2 variability of 60-100 ppm, while sea level fluctuations of typically 125 m are documented by proxy data. However, in the pre-ice core period, CO2 and sea level proxy data are scarce and there is disagreement between different proxies and different records of the same proxy. This hampers comprehensive understanding of the long-term relations between CO2, sea level and climate. Here, we drive a coupled climate-ice sheet model over the past five million years, inversely forced by a stacked benthic d18O record. We obtain continuous simulations of benthic d18O, sea level and CO2 that are mutually consistent. Our model shows CO2 concentrations of 300 to 470 ppm during the Early Pliocene. Furthermore, we simulate strong CO2 variability during the Pliocene and Early Pleistocene. These features are broadly supported by existing and new d11B-based proxy CO2 data, but less by alkenone-based records. The simulated concentrations and variations therein are larger than expected from global mean temperature changes. Our findings thus suggest a smaller Earth System Sensitivity than previously thought. This is explained by a more restricted role of land ice variability in the Pliocene. The largest uncertainty in our simulation arises from the mass balance formulation of East Antarctica, which governs the variability in sea level, but only modestly affects the modeled CO2 concentrations.

Multi-proxy analyses of sediment core GeoB13801-2

Over the Uruguayan shelf and uppermost slope the coalescence of northward flowing Subantarctic Shelf Water and southward flowing Subtropical Shelf Water forms a distinct thermohaline front termed the Subtropical Shelf Front (STSF). Running in a SW direction diagonally across the shelf from the coastal waters at 32°S towards the shelf break at ca. 36°S, the STSF represents the shelf-ward extension of the Brazil-Malvinas Confluence zone. This study reconstructs latitudinal STSF shifts during the Holocene based on benthic foraminifera d18O and d13C, total organic carbon, carbonate contents, Ti/Ca, and grain-size distribution from a high-accumulation sedimentary record located at an uppermost continental-slope terrace. Our data provide direct evidence for: (1) a southern STSF position (to the South of the core site) at the beginning of the early Holocene (>9.4 cal ka BP) linked to a more southerly position of the Southern Westerly Winds in combination with restricted shelf circulation intensity due to lower sea level; (2) a gradual STSF northward migration (bypassing the core site towards the North) primarily forced by the northward migration of the Southern Westerly Winds from 9.4 cal ka BP onwards; (3) a relatively stable position of the front in the interval between 7.2 and 4.0 cal ka BP; (4) millennial-scale latitudinal oscillations close to 36°S of the STSF after 4.0 cal ka BP probably linked to the intensification in El Niño Southern Oscillation; and (5) a southward migration of the STSF during the last 200 years possibly linked to anthropogenic influences on the atmosphere.

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.

A Late Pleistocene sea level stack

Late Pleistocene sea level has been reconstructed from ocean sediment core data using a wide variety of proxies and models. However, the accuracy of individual reconstructions is limited by measurement error, local variations in salinity and temperature, and assumptions particular to each technique. Here we present a sea level stack (average) which increases the signal-to-noise ratio of individual reconstructions. Specifically, we perform principal component analysis (PCA) on seven records from 0-430 ka and five records from 0-798 ka. The first principal component, which we use as the stack, describes ~80 % of the variance in the data and is similar using either five or seven records. After scaling the stack based on Holocene and Last Glacial Maximum (LGM) sea level estimates, the stack agrees to within 5 m with isostatically adjusted coral sea level estimates for Marine Isotope Stages 5e and 11 (125 and 400 ka, respectively). When we compare the sea level stack with the d18O of benthic foraminifera, we find that sea level change accounts for about ~40 % of the total orbital-band variance in benthic d18O, compared to a 65 % contribution during the LGM-to-Holocene transition. Additionally, the second and third principal components of our analyses reflect differences between proxy records associated with spatial variations in the d18O of seawater.