Marine Isotope Stage (MIS) 11 results of model data with Community Climate System Model version 3 including the Dynamic Global Vegetation Model (CCSM3-DGVM) in NetCDF format at global and regional scale

The climate of Marine Isotope Stage (MIS) 11, the interglacial roughly 400,000 years ago, is investigated for four time slices, 416, 410, 400, and 394 ka. The overall picture is that MIS 11 was a relatively warm interglacial in comparison to preindustrial, with Northern Hemisphere (NH) summer temperatures early in MIS 11 (416-410 ka) warmer than preindustrial, though winters were cooler. Later in MIS 11, especially around 400 ka, conditions were cooler in the NH summer, mainly in the high latitudes. Climate changes simulated by the models were mainly driven by insolation changes, with the exception of two local feedbacks that amplify climate changes. Here, the NH high latitudes, where reductions in sea ice cover lead to a winter warming early in MIS 11, as well as the tropics, where monsoon changes lead to stronger climate variations than one would expect on the basis of latitudinal mean insolation change alone, are especially prominent. The results support a northward expansion of trees at the expense of grasses in the high northern latitudes early during MIS 11, especially in northern Asia and North America.

Sea surface temperatures during Marine Isotope Stage (MIS) 11

The Marine Isotope Stage (MIS) 11 (424-374 ka) was characterized by a protracted deglaciation and an unusually long climatic optimum. It remains unclear to what degree the climate development during this interglacial reflects the unusually weak orbital forcing or greenhouse gas trends. Previously, arguments about the duration and timing of the MIS11 climatic optimum and about the pace of the deglacial warming were based on a small number of key records, which appear to show regional differences. In order to obtain a global signal of climate evolution during MIS11, we compiled a database of 78 sea surface temperature (SST) records from 57 sites spanning MIS11, aligned these individually on the basis of benthic (N = 28) or planktonic (N = 31) stable oxygen isotope curves to a common time frame and subjected 48 of them to an empirical orthogonal function (EOF) analysis. The analysis revealed a high commonality among all records, with the principal SST trend explaining almost 49% of the variability. This trend indicates that on the global scale, the surface ocean underwent rapid deglacial warming during Termination V, in pace with carbon dioxide rise, followed by a broad SST optimum centered at ~410 kyr. The second EOF, which explained ~18% of the variability, revealed the existence of a different SST trend, characterized by a delayed onset of the temperature optimum during MIS11 at ~398 kyr, followed by a prolonged warm period lasting beyond 380 kyr. This trend is most consistently manifested in the mid-latitude North Atlantic and Mediterranean Sea and is here attributed to the strength of the Atlantic meridional overturning circulation. A sensitivity analysis indicates that these results are robust to record selection and to age-model uncertainties of up to 3-6 kyr, but more sensitive to SST seasonal attribution and SST uncertainties >1 °C. In order to validate the CCSM3 (Community Climate System Model, version 3) predictive potential, the annual and seasonal SST anomalies recorded in a total of 74 proxy records were compared with runs for three time slices representing orbital configuration extremes during the peak interglacial of MIS11. The modeled SST anomalies are characterized by a significantly lower variance compared to the reconstructions. Nevertheless, significant correlations between proxy and model data are found in comparisons on the seasonal basis, indicating that the model captures part of the long-term variability induced by astronomical forcing, which appears to have left a detectable signature in SST trends.

Paleoceanography of sediment cores SO202-27-6 during MIS 3

The modern subarctic Pacific is characterized by a steep salinity-driven surface water stratification, which hampers the supply of saline and nutrient-rich deeper waters into the euphotic zone, limiting productivity. However, the strength of the halocline might have varied in the past. Here, we present diatom oxygen (d18Odiat) and silicon (d30Sidiat) stable isotope data from the open subarctic North-East (NE) Pacific (SO202-27-6; Gulf of Alaska), in combination with other proxy data (Neogloboquadrina pachydermasin d18O, biogenic opal, Ca and Fe intensities, IRD), to evaluate changes in surface water hydrography and productivity during Marine Isotope Stage (MIS) 3, characterized by millennial-scale temperature changes (Dansgaard-Oeschger (D-O) cycles) documented in Greenland ice cores.

Stable carbon isotope ratios of Cibicidoides wuellerstorfi of MIS 2 time slices

The interval of time represented by marine isotope stages 11 and 12 (~360-470 ka) contains what may be the most extreme glacial and interglacial climate conditions of the Late Pleistocene. It has been suggested that sea level rose by ~160 m at the termination of glacial stage 12. This is 30% greater than the sea level rise that followed the most recent glacial maximum. There have been few detailed studies of the unique conditions that existed during the stage 11-12 time period because of the lack of high-quality core material. This problem has been addressed by the collection of high deposition rate cores from sediment drifts in the western North Atlantic during Ocean Drilling Project Leg 172. Benthic foraminiferal d13C data from cores collected between ~4600 and 1800 m were used to reconstruct bathymetric gradients in deep and intermediate water properties for selected time slices during this glacial-interglacial cycle. During glacial stage 12, the deep western North Atlantic was filled by a water mass that was more nutrient-enriched than modern Antarctic Bottom Water. Above 2000 m, a more nutrient-depleted water mass existed during this glacial stage. Such an intermediate water mass has been described for more recent glacial periods and presumably forms in a more proximate region of the North Atlantic. Interglacial stage 11 water mass properties closely resemble those of the present-day western North Atlantic. A nutrient-depleted water mass (d13C of 0.75-1.0 per mil), similar to modern North Atlantic Deep Water existed between 3500 and 2000 m. This was underlain by a water mass with lower d13C values (<0.75 per mil) that probably was derived from a southern source. Using Leg 172 data, along with previously published results from the Atlantic and Pacific oceans, we estimate a mean global d13C change of 0.95 per mil from stage 12 to stage 11. This is twice the whole ocean ?13C change reported for the transition from the last glacial maximum to the Holocene.

Sea surface temperature estimation for MIS 16 of IODP Site 306-U1313

Hudson Strait (HS) Heinrich Events, ice-rafting events in the North Atlantic originating from the Laurentide ice sheet (LIS), are among the most dramatic examples of millennial-scale climate variability and have a large influence on global climate. However, it is debated as to whether the occurrence of HS Heinrich Events in the (eastern) North Atlantic in the geological record depends on greater ice discharge, or simply from the longer survival of icebergs in cold waters. Using sediments from Integrated Ocean Drilling Program (IODP) Site U1313 in the North Atlantic spanning the period between 960 and 320 ka, we show that sea surface temperatures (SSTs) did not control the first occurrence of HS Heinrich(-like) Events in the sedimentary record. Using mineralogy and organic geochemistry to determine the characteristics of ice-rafting debris (IRD), we detect the first HS Heinrich(-like) Event in our record around 643 ka (Marine Isotope Stage (MIS) 16), which is similar as previously reported for Site U1308. However, the accompanying high-resolution alkenone-based SST record demonstrates that the first HS Heinrich(-like) Event did not coincide with low SSTs. Thus, the HS Heinrich(-like) Events do indicate enhanced ice discharge from the LIS at the end of the Mid-Pleistocene Transition, not simply the survivability of icebergs due to cold conditions in the North Atlantic.

A multiproxy analysis of late MIS 6 to early MIS 5d OSP Site 161-975 sediments

A multiproxy analysis based on planktic foraminiferal abundances, derived SSTs, and stable planktic isotopes measurements together with alkenone abundances and Uk'37 SSTs was performed on late MIS 6 to early MIS 5d sediment recovered from Site 975 (ODP Leg 161) in the South Balearic Islands basin (Western Mediterranean) with emphasis on reconstructing the climate progression of the last interglacial period. A number of abrupt climate changes related to alternative influence of nutrient rich northern and oligotrophic southern water masses were revealed. Heinrich event 11 and cooling events C27, C26, C25, C24, C23, which have been previously described in the North Atlantic, were recognized. However, in comparison to the eastern North Atlantic mid-latitude region, events C27 and C26 at Site 975 seem to be significantly more pronounced. Together with evidence of a two-phase climate optimum with maximum SSTs reached during its later phase, this implies a close similarity in climate dynamics between the Western Mediterranean and the Nordic seas. We propose that postglacial effects in the Nordic seas had an influence on the western Mediterranean climate via atmospheric circulation and that these effects competed with the insolation force.

Pb isotopes in individual sand-sized (>150 µm) ice-rafted feldspars deposited at subpolar North Atlantic Ocean Deep Sea Drilling Project Site 611 during the late Pliocene (Marine Isotope Stages, MIS, G6 through 100)

The onset of abundant ice-rafted debris (IRD) deposition in the Nordic Seas and subpolar North Atlantic Ocean 2.72 millions of years ago (Ma) is thought to record the Pliocene onset of major northern hemisphere glaciation (NHG) due to a synchronous advance of North American Laurentide, Scandinavian and Greenland ice-sheets to their marine calving margins during marine isotope stage (MIS) G6. Numerous marine and terrestrial records from the Nordic Seas region indicate that extensive ice sheets on Greenland and Scandinavia increased IRD inputs to these seas from 2.72 Ma. The timing of ice-sheet expansion on North America as tracked by IRD deposition in the subpolar North Atlantic Ocean, however, is less clear because both Europe and North America are potential sources for icebergs in this region. Moreover, cosmogenic-dating of terrestrial tills on North America indicate that the Laurentide Ice Sheet did not extend to ~39°N until 2.4 ±0.14 Ma, at least 180 ka after the onset of major IRD deposition at 2.72 Ma. To address this problem,we present the first detailed analysis of the geochemical provenance of individual sand-sized IRD deposited in the subpolar North Atlantic Ocean between MIS G6 and 100 (~2.72-2.52 Ma). IRD provenance is assessed using laser ablation lead (Pb) isotope analyses of single ice-rafted (>150 mm) feldspar grains. To track when an ice-rafting setting consistent with major NHG first occurred in the North Atlantic Ocean during the Pliocene intensification of NHG (iNHG), we investigate when the Pb-isotope composition (206Pb/204Pb, 207Pb/204Pb, 208Pb/204Pb) of feldspars deposited at DSDP Site 611 first resembles that determined for IRD deposited at this site during MIS 100, the oldest glacial for which there exists convincing evidence for widespread glaciation of North America. Whilst Quaternary-magnitude IRD fluxes exist at Site 611 during glacials from 2.72 Ma, we find that the provenance of this IRD is not constant. Instead, we find that the Pb isotope composition of IRD at our study site is not consistent with major NHG until MIS G2 (2.64 Ma). We hypothesise that IRD deposition in the North Atlantic Ocean prior to MIS G2 was dominated by iceberg calving from Greenland and Scandinavia. We further suggest that the grounding line of continental ice on Northeast America may not have extended onto the continental shelf and calved significant numbers of icebergs to the North Atlantic Ocean during glacials until 2.64 Ma.

Biomarker and x-ray diffraction data from Site U1313 (MIS 16 - 9)

A reconstruction of Milankovitch to millennial-scale variability of sea-surface temperature (SST) and sea-surface productivity in the Pleistocene mid-latitude North Atlantic Ocean (MIS 16-9) and its relationship to ice sheet instability was carried out on sediments from IODP Site U1313. This reconstruction is based on alkenone and n-alkane concentrations, Uk37' index, total organic carbon (TOC) and carbonate contents, X-Ray diffraction (XRD) data, magnetic susceptibility, and accumulation rates. Increased input of ice-rafted debris (IRD) occurred during MIS 16, 12, and 10, characterized by high concentrations of dolomite, quartz, and feldspars and elevated accumulation rates of terrigenous matter. Minimum input values of terrigenous matter, on the other hand, were determined for MIS 13 and 11. Peak values of dolomite, coinciding with quartz, plagioclase, and kalifeldspar peaks and maxima in long-chain n-alkanes indicative for land plants, are interpreted as Heinrich-like Events related to sudden instability of the Laurentide Ice Sheet during early and late (deglacial) phases of the glacials. The coincidence of increased TOC values with elevated absolute concentrations of alkenones suggest increased glacial productivity, probably due to a more southern position of the Polar Front. Alkenone-based SST reached absolute maxima of about 19°C during MIS 11.3 and absolute minima of <10°C during MIS 12 and 10. Within MIS 11, prominent cooling events (MIS 11.22 and 11.24) occurred. The absolute SST minima recorded directly before and after the glacial maxima MIS 10.2 and 12.2, are related to Heinrich-like Event meltwater pulses, as supported by the coincidence of SST minima and maxima in C37:4 alkenones and dolomite. These sudden meltwater pulses - especially during Terminations IV and V - probably caused a collapse of phytoplankton productivity as indicated by the distinct drop in alkenone concentrations. Ice-sheet disintegration and subsequent surges and outbursts of icebergs and meltwater discharge may have been triggered by increased insolation in the Northern High Latitudes.

Stable carbon and oxygen isotope record of IODP Site 306-U1313, MIS 23-19

Since the seminal work by Hays et al. (1976), a plethora of studies has demonstrated a correlation between orbital variations and climatic change. However, information on how changes in orbital boundary conditions affected the frequency and amplitude of millennial-scale climate variability is still fragmentary. The Marine Isotope Stage (MIS) 19, an interglacial centred at around 785 ka, provides an opportunity to pursue this question and test the hypothesis that the long-term processes set up the boundary conditions within which the short-term processes operate. Similarly to the current interglacial, MIS 19 is characterised by a minimum of the 400-kyr eccentricity cycle, subdued amplitude of precessional changes, and small amplitude variations in insolation. Here we examine the record of climatic conditions during MIS 19 using high-resolution stable isotope records from benthic and planktonic foraminifera from a sedimentary sequence in the North Atlantic (Integrated Ocean Drilling Program Expedition 306, Site U1313) in order to assess the stability and duration of this interglacial, and evaluate the climate system's response in the millennial band to known orbitally induced insolation changes. Benthic and planktonic foraminiferal d18O values indicate relatively stable conditions during the peak warmth of MIS 19, but sea-surface and deep-water reconstructions start diverging during the transition towards the glacial MIS 18, when large, cold excursions disrupt the surface waters whereas low amplitude millennial scale fluctuations persist in the deep waters as recorded by the oxygen isotope signal. The glacial inception occurred at ~779 ka, in agreement with an increased abundance of tetra-unsaturated alkenones, reflecting the influence of icebergs and associated meltwater pulses and high-latitude waters at the study site. After having combined the new results with previous data from the same site, and using a variety of time series analysis techniques, we evaluate the evolution of millennial climate variability in response to changing orbital boundary conditions during the Early-Middle Pleistocene. Suborbital variability in both surface- and deep-water records is mainly concentrated at a period of ~11 kyr and, additionally, at ~5.8 and ~3.9 kyr in the deep ocean; these periods are equal to harmonics of precession band oscillations. The fact that the response at the 11 kyr period increased over the same interval during which the amplitude of the response to the precessional cycle increased supports the notion that most of the variance in the 11 kyr band in the sedimentary record is nonlinearly transferred from precession band oscillations. Considering that these periodicities are important features in the equatorial and intertropical insolation, these observations are in line with the view that the low-latitude regions play an important role in the response of the climate system to the astronomical forcing. We conclude that the effect of the orbitally induced insolation is of fundamental importance in regulating the timing and amplitude of millennial scale climate variability.