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.

Salinity estimation from Gulf of Guinea, sediment core MD03-2707

The isotopic composition of surface seawater is widely used to infer past changes in sea surface salinity using paired foraminiferal Mg/Ca and d18O from marine sediments. At low latitudes, paleosalinity reconstructions using this method have largely been used to document changes in the hydrological cycle. This method usually assumes that the modern seawater d18O (d18Osw)/salinity relationship remained constant through time. Modelling studies have shown that such assumptions may not be valid because large-scale atmospheric circulation patterns linked to global climate changes can alter the seawater d18Osw/salinity relationship locally. Such processes have not been evidenced by paleo-data so far because there is presently no way to reconstruct past changes in the seawater d18Osw/salinity relationship. We have addressed this issue by applying a multi-proxy salinity reconstruction from a marine sediment core collected in the Gulf of Guinea. We measured hydrogen isotopes in C37:2 alkenones (dDa) to estimate changes in seawater dD. We find a smooth, long-term increase of ~10 per mil in dDa between 10 and 3 kyr BP, followed by a rapid decrease of ~10 per mil in dDa between 3 kyr BP and core top to values slightly lighter than during the early Holocene. Those features are inconsistent with published salinity estimations based on d18Osw and foraminiferal Ba/Ca, as well as nearby continental rainfall history derived from pollen analysis. We combined dDa and d18Osw values to reconstruct a Holocene record of salinity and compared it to a Ba/Ca-derived salinity record from the same sedimentary sequence. This combined method provides salinity trends that are in better agreement with both the Ba/Ca-derived salinity and the regional precipitation changes as inferred from pollen records. Our results illustrate that changes in atmospheric circulation can trigger changes in precipitation isotopes in a counter-intuitive manner that ultimately impacts surface salinity estimates based on seawater isotopic values. Our data suggest that the trends in Holocene rainfall isotopic values at low latitudes may not uniquely result from changes in local precipitation associated with the amount effect.