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.

Magnesium/Calcium ratios measured on planktonic foraminifera of the LGM

As part of the Multi-proxy Approach for the Reconstruction of the Glacial Ocean (MARGO) incentive, published and unpublished temperature reconstructions for the Last Glacial Maximum (LGM) based on planktonic foraminiferal Mg/Ca ratios have been synthesised and made available in an online database. Development and applications of Mg/Ca thermometry are described in order to illustrate the current state of the method. Various attempts to calibrate foraminiferal Mg/Ca ratios with temperature, including culture, trap and core-top approaches have given very consistent results although differences in methodological techniques can produce offsets between laboratories which need to be assessed and accounted for where possible. Dissolution of foraminiferal calcite at the sea-floor generally causes a lowering of Mg/Ca ratios. This effect requires further study in order to account and potentially correct for it if dissolution has occurred. Mg/Ca thermometry has advantages over other paleotemperature proxies including its use to investigate changes in the oxygen isotopic composition of seawater and the ability to reconstruct changes in the thermal structure of the water column by use of multiple species from different depth and or seasonal habitats. Presently available data are somewhat limited to low latitudes where they give fairly consistent values for the temperature difference between Late Holocene and the LGM (2-3.5 °C). Data from higher latitudes are more sparse, and suggest there may be complicating factors when comparing between multi-proxy reconstructions.

Mg/Ca ratios in planktonic foraminifera

The primary Mg/Ca ratio of foraminiferal shells is a potentially valuable paleoproxy for sea surface temperature (SST) reconstructions. However, the reliable extraction of this ratio from sedimentary calcite assumes that we can overcome artifacts related to foraminiferal ecology and partial dissolution, as well as contamination by secondary calcite and clay. The standard batch method for Mg/Ca analysis involves cracking, sonicating, and rinsing the tests to remove clay, followed by chemical cleaning, and finally acid-digestion and single-point measurement. This laborious procedure often results in substantial loss of sample (typically 30-60%). We find that even the earliest steps of this procedure can fractionate Mg from Ca, thus biasing the result toward a more variable and often anomalously low Mg/Ca ratio. Moreover, the more rigorous the cleaning, the more calcite is lost, and the more likely it becomes that any residual clay that has not been removed by physical cleaning will increase the ratio. These potentially significant sources of error can be overcome with a flow-through (FT) sequential leaching method that makes time- and labor-intensive pretreatments unnecessary. When combined with time-resolved analysis (FT-TRA) flow-through, performed with a gradually increasing and highly regulated acid strength, produces continuous records of Mg, Sr, Al, and Ca concentrations in the leachate sorted by dissolution susceptibility of the reacting material. Flow-through separates secondary calcite from less susceptible biogenic calcite and clay, and further resolves the biogenic component into primary and more resistant fractions. FT-TRA reliably separates secondary calcite (which is not representative of original life habitats) from the more resistant biogenic calcite (the desired signal) and clay (a contaminant of high Mg/Ca, which also contains Al), and further resolves the biogenic component into primary and more resistant fractions that may reflect habitat or other changes during ontogeny. We find that the most susceptible fraction of biogenic calcite in surface dwelling foraminifera gives the most accurate value for SST and therefore best represents primary calcite. Sequential dissolution curves can be used to correct the primary Mg/Ca ratio for clay, if necessary. However, the temporal separation of calcite from clay in FT-TRA is so complete that this correction is typically <=2%, even in clay-rich sediments. Unlike hands-on batch methods, that are difficult to reproduce exactly, flow-through lends itself to automation, providing precise replication of treatment for every sample. Our automated flow-through system can process 22 samples, two system blanks, and 48 mixed standards in <12 hours of unattended operation. FT-TRA thus represents a faster, cheaper, and better way to determine Mg/Ca ratios in foraminiferal calcite.

Stable organic carbon isotope record and total organic carbon content from the Cenomanian-Turonian boundary

Ocean Drilling Program (ODP) Leg 207 recovered expanded sections of organic-carbon-rich laminated shales on Demerara Rise (western tropical Atlantic). High-resolution organic carbon isotope and total organic carbon (TOC) records are presented, which span the Cenomanian-Turonian boundary interval (CTBI), including the Oceanic Anoxic Event (OAE) 2, from four sites oriented along a NW striking depth transect. These records represent the first high-resolution carbon isotope records across OAE 2 from the South American margin of the tropical Atlantic. Due to the scarcity of age significant fossils, the main purpose of this study was to develop a detailed carbon isotope stratigraphy in order to correlate the CTBI across the depth transect and to tie this to biostratigraphically well-defined sections in the Western Interior Basin (Pueblo, USA), boreal shelf seas (Eastbourne, England), and western Tethys (Oued Mellegue, Tunisia). All four sections studied document a 6 per mil increase of d13Corg values at the base of the CTBI, which is followed by an interval of elevated d13Corg values and a subsequent decrease. Our results supply an important stratigraphic base for subsequent paleoceanographic studies on Late Cenomanian to Early Turonian sediments from Demerara Rise and elsewhere.

Calcite saturation and yield and Mg/Ca, Mn/Ca, Al/Ca, Fe/Ca for four species of planktic foraminifera cleaned using Mg-cleaning and Cd-cleaning methods

Four species of planktic foraminifera from core-tops spanning a depth transect on the Ontong Java Plateau were prepared for Mg/Ca analysis both with (Cd-cleaning) and without (Mg-cleaning) a reductive cleaning step. Reductive cleaning caused etching of foraminiferal calcite, focused on Mg-rich inner calcite, even on tests which had already been partially dissolved at the seafloor. Despite corrosion, there was no difference in Mg/Ca of Pulleniatina obliquiloculata between cleaning methods. Reductive cleaning decreased Mg/Ca by an average (all depths) of ~ 4% for Globigerinoides ruber white and ~ 10% for Neogloboquadrina dutertrei. Mg/Ca of Globigerinoides sacculifer (above the calcite saturation horizon only) was 5% lower after reductive cleaning. The decrease in Mg/Ca due to reductive cleaning appeared insensitive to preservation state for G. ruber, N. dutertrei and P. obliquiloculata. Mg/Ca of Cd-cleaned G. sacculifer appeared less sensitive to dissolution than that of Mg-cleaned. Mg-cleaning is adequate, but SEM and contaminants (Al/Ca, Fe/Ca and Mn/Ca) show that Cd-cleaning is more effective for porous species. A second aspect of the study addressed sample loss during cleaning. Lower yield after Cd-cleaning for G. ruber, G. sacculifer and N. dutertrei confirmed this to be the more aggressive method. Strongest correlations between yield and Delta[CO3^2-] in core-top samples were for Cd-cleaned G. ruber (r = 0.88, p = 0.020) and Cd-cleaned P. obliquiloculata (r = 0.68, p = 0.030). In a down-core record (WIND28K) correlation, r, between yield values > 30% and dissolution index, XDX, was -0.61 (p = 0.002). Where cleaning yield < 30% most Mg-cleaned Mg/Ca values were biased by dissolution.