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 record and planktonic foraminiferal fragmentation of the late Pleistocene section of ODP Hole 134-828A

A high-resolution record of foraminiferal fragmentation (a dissolution indicator) for the last 250 k.y. (isotopic Stages 1 to 7) is identified in the upper 61.9 m of Ocean Drilling Program (ODP) Hole 828A, west Vanuatu. This record is comparable in detail to the atmospheric CO2 record and the d18O stack. Phase shifts between preservation spikes and maximum ice volumes (d18O of Globigerinoides sacculifer) are analogous to those on Ontong Java Plateau. Mass spectrometer (AMS14C) dating of a sample taken at the base of dissolution cycle B1 and the position of the last glacial maximum indicates a lag in time of ~8 k.y. in the Vanuatu region for the last glacial termination. When dissolution spikes are compared with minimum ice volumes there is no phase shift for the last two glacial terminations. The difference between Vanuatu and Ontong Java Plateau may be explained by local CO2 sinks and the interplay between intermediate and deep water masses. Terrigenous input increasingly affected sediment of Hole 828A on the North d'Entrecasteaux Ridge (NDR) as it approached Espiritu Santo Island. Mud and silt suspended in mid-water flows become important after 125 ka, while turbidites bypass the New Hebrides Trench only towards the last glacial maximum (LGM). Terrigenous supply seems to affect the lysocline profile that changed from an "open ocean" to a "near continent" type, thus favoring dissolution. Fragmentation of planktonic foraminifers is a more sensitive indicator of lysocline variations than is foraminiferal susceptibility to dissolution, the foraminiferal dissolution index, the abundance of benthic foraminifers, or CaCO3 content. A modern foraminiferal lysocline for the neighboring area (between 10°S and 30°S, and 160°E and 180°E) is found at 3.1 km below sea level, compared to west Vanuatu where it is shallower. The past lysocline level was deeper than 3086 m during intervals of dissolution minima, and ranged from ~2550 to 3000 m during intervals of dissolution maxima. The high sedimentation rates (in the order of 10 to 50 cm/k.y.) found in Hole 828A offer a great potential for future high-resolution studies either in this hole or other western localities along the NDR. Areas of high sedimentation near continental regions have been discarded for paleoceanographic and/or paleoclimatic studies. Nonetheless, conditions analogous to those found in Hole 828A are expected to occur in many trench areas around the world where mid-water flows have preserved as yet undiscovered fine high-resolution sedimentary records.

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.

Quaternary carbonate and stable isotope record of ODP Holes 133-817A and 133-818B

The Quaternary history of metastable CaCO3 input and preservation within Antarctic Intermediate Water (AAIW) was examined by studying sediments from ODP Holes 818B (745 mbsl) and 817A (1015 mbsl) drilled in the Townsville Trough on the southern slope of the Queensland Plateau. These sites lie within the core of modern AAIW, and near the aragonite saturation depth (~1000 m). Thus, they are well positioned to monitor chemical changes that may have occurred within this watermass during the past 1.6 m.y. The percent of fine aragonite content, percent of fine magnesian calcite content, and percent of whole pteropods (>355 µm) were used to separate the fine aragonite input signal from the CaCO3 preservation signal. Stable d18O and d13C isotopic ratios were determined for the planktonic foraminifer Globigerinoides sacculifer and, in Hole 818B, for the benthic foraminifer Cibicidoides spp. to establish the oxygen isotope stratigraphy and to study the relationship between intermediate and shallow water d13C of Sum CO2 and the relationship between benthic foraminiferal d13C and CaCO3 preservation within intermediate waters of the Townsville Trough. Data were converted from depth to age using oxygen isotope stratigraphy, nannostratigraphy, and foraminiferal biostratigraphy. Several long hiatuses and the absence of magnetostratigraphy did not permit time series analysis. The principal results of the CaCO3 preservation study include the following (1) a general increase in CaCO3 preservation between 0.9 and 1.6 Ma; (2) a CaCO3 dissolution maximum near 0.9 Ma, primarily expressed in the Hole 818B fine aragonite record; (3) an abrupt and permanent increase of fine aragonite content between 0.86 and 0.875 Ma in both Holes 818B and 817A probably reflecting a dramatic increase of fine carbonate sediment production on the Queensland Plateau; (4) an improvement in CaCO3 preservation near 0.87 Ma, which accompanied the increase of sediment input, indicated by the first appearance of whole pteropods in the deeper Hole 817A and a "spike" in the percent whole pteropods in Hole 818B; (5) a period of strong CaCO3 dissolution during the mid-Brunhes Chron from 0.36 to 0.41 Ma; and (6) a complex CaCO3 preservation pattern between 0.36 Ma and the present characterized by a general increase in CaCO3 preservation through time with good preservation during interglacial stages and poor preservation during glacial stages. The long-term aragonite preservation histories for Holes 818B and 817A appear to be similar in general shape, although different in detail, to CaCO3 preservation records from the deep Indian and central equatorial Pacific oceans as well as from intermediate water sites in the Bahamas and the Maldives. All of these areas have experienced CaCO3 dissolution at about 0.9 Ma and during the mid-Brunhes Chron. However, the late Quaternary (0 to 0.36 Ma) glacial to interglacial preservation pattern in Holes 818B and 817A is out of phase with CaCO3 preservation records for sediments deposited in Pacific deep and bottom waters. The sharp increase in bank production and export from the Queensland Plateau and the coincident improvement of CaCO3 preservation between 0.86 and 0.875 Ma may have been synchronous with the initiation of the Great Barrier Reef and roughly coincides with an increase in carbonate accumulation on the Bahama banks, in the western North Atlantic Ocean, and on Mururoa atoll, in the central South Pacific Ocean. The development of these reef systems during the middle Quaternary may be related to the transition in the frequency and amplitude of global sea level change from 41 k.y. low amplitude cycles prior to 0.9 Ma to 100 k.y. high amplitude cycles after 0.73 Ma. Carbon isotopic analyses show that benthic foraminiferal d13C values (Cibicidoides spp.) have been heavier than planktonic foraminiferal d13C values (G. sacculifer) throughout most of the last 0.54 m.y., which may indicate that 13C-enriched intermediate water (AAIW) occupied the Townsville Trough during much of the late Quaternary. Furthermore, both planktonic and benthic foraminiferal d13C values are often observed to be heaviest during interglacial to glacial transitions, and lightest during glacial to interglacial transitions. We suggest that this pattern is the result of changes in the preformed d13C of Sum CO2 of AAIW and may reflect changes in nutrient utilization by primary producers in Antarctic surface waters, changes in the d13C of upwelled Circumpolar Deep Water, or changes in the extent and/or temperature of equilibration between surface water and atmospheric CO2 within the Antarctic Polar Frontal Zone (the source area for AAIW). Finally, the poor correlation between percent of whole pteropods (aragonite preservation) and d13C of Cibicidoides spp. may be the result of a decoupling of d13C from CO2 due to the numerous and complex variables that combine to produce the preformed d13C of AAIW.

(Table 3) Stable oxagen isotope record of planktonic foraminifera from Central Walvish Ridge

Above the Walvis Ridge, in the SE Atlantic Ocean, we collected living plantkic foraminifera from the upper water column using depth stratified plantkon tows. The oxygen isotope composition (d18Oc) in shells of foraminifera and shell concentration profiles show seasonal and depth habitats of individual species. The tow results are compared with the average annual deposition d18Oc from sediment traps and the interannual average d18Oc of fossil specimens in top sediments at the same site. The species Globigerinita glutinata best reflects the austral winter/spring sea surface temperature (SST). Its d18Oc signal in top sediments remains pristine. In contrast, tow results also show that Globigerinoides ruber continues to calcify below the surface mixed layer (SML), i.e., down to the deep chlorophyll maximum (DCM); hence its d18Oc signature of exported specimens reflects the SST only when SML incorporates the DCM. Deep tow and sediment trap results show that both Globorotalia truncatulinoides and Globorotalia inflata record the temperature between 150 and 350 m, depending on the season and the shell size. However, for all fossil taxa in sediments apart from Globigerinita glutinata, we observe a positive d18Oc shift with respect to the sediment trap and plankton tow values, likely related to the interannual flux changes and deep encrustation.