Constraints on the magnitude and patterns of ocean cooling at the Last Glacial Maximum - Supplementary material

Observation-based reconstructions of sea surface temperature from relatively stable periods in the past, such as the Last Glacial Maximum, represent an important means of constraining climate sensitivity and evaluating model simulations. The first quantitative global reconstruction of sea surface temperatures during the Last Glacial Maximum was developed by the Climate Long-Range Investigation, Mapping and Prediction (CLIMAP) project in the 1970s and 1980s. Since that time, several shortcomings of that earlier effort have become apparent. Here we present an updated synthesis of sea surface temperatures during the Last Glacial Maximum, rigorously defined as the period between 23 and 19 thousand years before present, from the Multiproxy Approach for the Reconstruction of the Glacial Ocean Surface (MARGO) project. We integrate microfossil and geochemical reconstructions of surface temperatures and include assessments of the reliability of individual records. Our reconstruction reveals the presence of large longitudinal gradients in sea surface temperature in all of the ocean basins, in contrast to the simulations of the Last Glacial Maximum climate available at present.

Stable isotope ratios and foraminiferal abundance of sediment cores from the Western Iberian Margin

Rapid climate changes at the onset of the last deglaciation and during Heinrich Event H4 were studied in detail at IMAGES cores MD95-2039 and MD95-2040 from the Western Iberian margin. A major reorganisation of surface water hydrography, benthic foraminiferal community structure, and deepwater isotopic composition commenced already 540 years before the Last Isotopic Maximum (LIM) at 17.43 cal. ka and within 670 years affected all environments. Changes were initiated by meltwater spill in the Nordic Seas and northern North Atlantic that commenced 100 years before concomitant changes were felt off western Iberia. Benthic foraminiferal associations record the drawdown of deepwater oxygenation during meltwater and subsequent Heinrich Events H1 and H4 with a bloom of dysoxic species. At a water depth of 3380 m, benthic oxygen isotopes depict the influence of brines from sea ice formation during ice-rafting pulses and meltwater spill. The brines conceivably were a source of ventilation and provided oxygen to the deeper water masses. Some if not most of the lower deep water came from the South Atlantic. Benthic foraminiferal assemblages display a multi-centennial, approximately 300-year periodicity of oxygen supply at 2470-m water depth. This pattern suggests a probable influence of atmospheric oscillations on the thermohaline convection with frequencies similar to Holocene climate variations. For Heinrich Events H1 and H4, response times of surface water properties off western Iberia to meltwater injection to the Nordic Seas were extremely short, in the range of a few decades only. The ensuing reduction of deepwater ventilation commenced within 500-600 years after the first onset of meltwater spill. These fast temporal responses lend credence to numerical simulations that indicate ocean-climate responses on similar and even faster time scales.

Isotopic composition of terrestrial plant waxes (dD and d13C of n-alkanes) of sediment core SO188_342 from the Northern Bay of Bengal for the last 18 ka

The Indian Summer Monsoon (ISM) is a major global climatic phenomenon. Long-term precipitation proxy records of the ISM, however, are often fragmented and discontinuous, impeding an estimation of the magnitude of precipitation variability from the Last Glacial to the present. To improve our understanding of past ISM variability, we provide a continuous reconstructed record of precipitation and continental vegetation changes from the lower Ganges-Brahmaputra-Meghna catchment and the Indo-Burman ranges over the last 18,000 years (18 ka). The records derive from a marine sediment core from the northern Bay of Bengal (NBoB), and are complemented by numerical model results of spatial moisture transport and precipitation distribution over the Bengal region. The isotopic composition of terrestrial plant waxes (dD and d13C of n-alkanes) are compared to results from an isotope-enabled general atmospheric circulation model (IsoCAM) for selected time slices (pre-industrial, mid-Holocene and Heinrich Stadial 1). Comparison of proxy and model results indicate that past changes in the dD of precipitation and plant waxes were mainly driven by the amount effect, and strongly influenced by ISM rainfall. Maximum precipitation is detected for the Early Holocene Climatic Optimum (EHCO; 10.5-6 ka BP), whereas minimum precipitation occurred during the Heinrich Stadial 1 (HS1; 16.9-15.4 ka BP). The IsoCAM model results support the hypothesis of a constant moisture source (i.e. the NBoB) throughout the study period. Relative to the pre-industrial period the model reconstructions show 20% more rain during the mid-Holocene (6 ka BP) and 20% less rain during the Heinrich Stadial 1 (HS1), respectively. A shift from C4-plant dominated ecosystems during the glacial to subsequent C3/C4-mixed ones during the interglacial took place. Vegetation changes were predominantly driven by precipitation variability, as evidenced by the significant correlation between the dD and d13C alkane records. When compared to other records across the ISM domain, precipitation and vegetation changes inferred from our records and the numerical model results provide evidence for a coherent regional variability of the ISM from the Last Glacial to the present.

Plant-wax, alkenone and major element analyses of ODP Site 659 for the Last Glacial cycle and associated coretop data

We present a hydrologic reconstruction of the Sahara-Sahel transition, covering the complete last glacial cycle (130 ka), based on a combination of plant-wax-specific hydrogen (dD) and carbon isotopes (d13C). The dD and d13C signatures of long-chain n-alkanes from ODP Site 659 off NW Africa reveal a significant anti-correlation. Complementary to published pollen data, we infer that this plant-wax signal reflects sensitive responses of the vegetation cover to precipitation changes in the Sahel region, as well as varying contributions from biomes north of the Sahara (C3 domain) by North-East Trade Winds (NETW). During arid phases, especially the northern parts of the Sahel likely experienced crucial water stress, which resulted in a pronounced contraction of the vegetation cover, thus reducing the amount of C4 plant waxes from the region. The increase in NETW strength during dry periods further promoted a more pronounced C3-plant-wax signal derived from the North African C3 plant domain. During humid periods, the C4-dominated Sahelian environments spread northward into the Saharan realm, in association with lower NETW inputs of C3 plant waxes. Arid-humid cycles deduced from plant-wax dD are in accordance with concomitant changes in weathering intensity reflected in varying major element distributions. Environmental shifts are generally linked to periods with large fluctuations in Northern Hemisphere summer insolation. During Marine Isotope Stages 2 and 3, when insolation variability was low, coupling of the hydrologic regime to alkenone-based estimates of NE Atlantic sea-surface temperatures becomes apparent.

Planktic foraminiferal Mg/Ca and oxygen isotopes, plus estimated d18Oseawater and ice volume corrected d18Oseawater for the past 35000 years of sediment core GeoB10069-3

Evidence from geologic archives suggests that there were large changes in the tropical hydrologic cycle associated with the two prominent northern hemisphere deglacial cooling events, Heinrich Stadial 1 (HS1; ~19 to 15 kyr BP; kyr BP = 1000 yr before present) and the Younger Dryas (~12.9 to 11.7 kyr BP). These hydrologic shifts have been alternatively attributed to high and low latitude origin. Here, we present a new record of hydrologic variability based on planktic foraminifera-derived d18O of seawater (d18Osw) estimates from a sediment core from the tropical Eastern Indian Ocean, and using 12 additional d18Osw records, construct a single record of the dominant mode of tropical Eastern Equatorial Pacific and Indo-Pacific Warm Pool (IPWP) hydrologic variability. We show that deglacial hydrologic shifts parallel variations in the reconstructed interhemispheric temperature gradient, suggesting a strong response to variations in the Atlantic Meridional Overturning Circulation and the attendant heat redistribution. A transient model simulation of the last deglaciation suggests that hydrologic changes, including a southward shift in the Intertropical Convergence Zone (ITCZ) which likely occurred during these northern hemisphere cold events, coupled with oceanic advection and mixing, resulted in increased salinity in the Indonesian region of the IPWP and the eastern tropical Pacific, which is recorded by the d18Osw proxy. Based on our observations and modeling results we suggest the interhemispheric temperature gradient directly controls the tropical hydrologic cycle on these time scales, which in turn mediates poleward atmospheric heat transport.

Stable isotope ratios of foraminifera from the continental shelf off Troms, Northern Norway

A stable oxygen and carbon isotope stratigraphy is established for a Late Weichselian/Holocene glaciomarine/marine seguence in Andfjorden and Malangsdjupet on the continental shelf off Troms, Northern Norway. The stratigraphy demonstrates that the global signals, Termination I B and possibly also I A (upper parts), are present and radiocarbon date to 10.3-9.7 kyr B.P. and >14-13.5 kyr B.P., respectively. A temperature increase of 5°-6°C and possibly a small salinity increase occurred during Term. I. A near-glacial environment between 13 and 14 kyr B.P. was characterized by poorly ventilated bottom waters followed by a meltwater pulse at circa 13 kyr B.P. During the beginning intrusion of Atlantic Water between 13 and 10 kyr B.P., the bottom water was characterized by somewhat fluctuating temperatures and salinities. Temperatures close to those of the present were established around 9.7 kyr B.P. and seem to have been rather stable since.

Composition of bottom sediments from the Isfjord, Western Svalbard

Lithology, heavy mineral associations, and chemical composition of bottom sediments studied in two gravity cores from Isfjord, Western Spitsbergen (Svalbard) accompanied by high-frequency seismic records, provide a new insight on provenance and glaciomarine sedimentation in the fjord from the last deglaciation through Holocene.

Changes in bottom water corrosiveness and carbonate ion concentrations in the sub-Antarctic Atlantic during the last glacial period

The glacial climate system transitioned rapidly between cold (stadial) and warm (interstadial) conditions in the Northern Hemisphere. This variability, referred to as Dansgaard-Oeschger variability, is widely believed to arise from perturbations of the Atlantic Meridional Overturning Circulation. Evidence for such changes during the longer Heinrich stadials has been identified, but direct evidence for overturning circulation changes during Dansgaard-Oeschger events has proven elusive. Here we reconstruct bottom water [CO3]2- variability from B/Ca ratios of benthic foraminifera and indicators of sedimentary dissolution, and use these reconstructions to infer the flow of northern-sourced deep water to the deep central sub-Antarctic Atlantic Ocean. We find that nearly every Dansgaard-Oeschger interstadial is accompanied by a rapid incursion of North Atlantic Deep Water into the deep South Atlantic. Based on these results and transient climate model simulations, we conclude that North Atlantic stadial-interstadial climate variability was associated with significant Atlantic overturning circulation changes that were rapidly transmitted across the Atlantic. However, by demonstrating the persistent role of Atlantic overturning circulation changes in past abrupt climate variability, our reconstructions of carbonate chemistry further indicate that the carbon cycle response to abrupt climate change was not a simple function of North Atlantic overturning.

Modeling the sensitivity of precipitation oxygen isotopes to the Last Glacial Maximum boundary conditions: A study using Community Atmosphere Model (CAM3.0)

To understand the validity of d18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (d18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and d18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and d18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in d18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the d18Oprecip distribution among the simulations.

Geochemistry of foraminifera from sediment cores of the central California margin

Under present climate conditions, convection at high latitudes of the North Pacific is restricted to shallower depths than in the North Atlantic. To what extent this asymmetry between the two ocean basins was maintained over the past 20 kyr is poorly known because there are few unambiguous proxy records of ventilation from the North Pacific. We present new data for two sediment cores from the California margin at 800 and 1600 m depth to argue that the depth of ventilation shifted repeatedly in the northeast Pacific over the course of deglaciation. The evidence includes benthic foraminiferal Cd/Ca, 18O/16O, and 13C/12C data as well as radiocarbon age differences between benthic and planktonic foraminifera. A number of features in the shallower of the two cores, including an interval of laminated sediments, are consistent with changes in ventilation over the past 20 kyr suggested by alternations between laminated and bioturbated sediments in the Santa Barbara Basin and the Gulf of California [Keigwin and Jones, 1990 doi:10.1029/PA005i006p01009; Kennett and Ingram, 1995 doi:10.1038/377510a0; Behl and Kennett, 1996 doi:10.1038/379243a0]. Data from the deeper of the two California margin cores suggest that during times of reduced ventilation at 800 m, ventilation was enhanced at 1600 m depth, and vice versa. This pronounced depth dependence of ventilation needs to be taken into account when exploring potential teleconnections between the North Pacific and the North Atlantic.

Separating ITCZ- and ENSO-related rainfall changes in the Galápagos over the last 3 kyr using D/H ratios of multiple lipid biomarkers

We present a 3000-yr rainfall reconstruction from the Galápagos Islands that is based on paired biomarker records from the sediment of El Junco Lake. Located in the eastern equatorial Pacific, the climate of the Galápagos Islands is governed by movements of the Intertropical Convergence Zone (ITCZ) and the El Niño-Southern Oscillation (ENSO). We use a novel method for reconstructing past ENSO- and ITCZ-related rainfall changes through analysis of molecular and isotopic biomarker records representing several types of plants and algae that grow under differing climatic conditions. We propose that ?D values of dinosterol, a sterol produced by dinoflagellates, record changes in mean rainfall in El Junco Lake, while dD values of C34 botryococcene, a hydrocarbon unique to the green alga Botryococcus braunii, record changes in rainfall associated with moderate-to-strong El Niño events. We use these proxies to infer changes in mean rainfall and El Niño-related rainfall over the past 3000 yr. During periods in which the inferred change in El Niño-related rainfall opposed the change in mean rainfall, we infer changes in the amount of ITCZ-related rainfall. Simulations with an idealized isotope hydrology model of El Junco Lake help illustrate the interpretation of these proxy reconstructions. Opposing changes in El Niño- and ITCZ-related rainfall appear to account for several of the largest inferred hydrologic changes in El Junco Lake. We propose that these reconstructions can be used to infer changes in frequency and/or intensity of El Niño events and changes in the position of the ITCZ in the eastern equatorial Pacific over the past 3000 yr. Comparison with El Junco Lake sediment grain size records indicates general agreement of inferred rainfall changes over the late Holocene.