Compilation of 220 sediment core d13C data from the glacial Atlantic Ocean

We compare a compilation of 220 sediment core d13C data from the glacial Atlantic Ocean with three-dimensional ocean circulation simulations including a marine carbon cycle model. The carbon cycle model employs circulation fields which were derived from previous climate simulations. All sediment data have been thoroughly quality controlled, focusing on epibenthic foraminiferal species (such as Cibicidoides wuellerstorfi or Planulina ariminensis) to improve the comparability of model and sediment core carbon isotopes. The model captures the general d13C pattern indicated by present-day water column data and Late Holocene sediment cores but underestimates intermediate and deep water values in the South Atlantic. The best agreement with glacial reconstructions is obtained for a model scenario with an altered freshwater balance in the Southern Ocean that mimics enhanced northward sea ice export and melting away from the zone of sea ice production. This results in a shoaled and weakened North Atlantic Deep Water flow and intensified Antarctic Bottom Water export, hence confirming previous reconstructions from paleoproxy records. Moreover, the modeled abyssal ocean is very cold and very saline, which is in line with other proxy data evidence.

Distribution of planktic foraminifera in surface sediments of the Atlantic Ocean

We present a data set of 738 planktonic foraminiferal species counts from sediment surface samples of the eastern North Atlantic and the South Atlantic between 87°N and 40°S, 35°E and 60°W including published Climate: Long-Range Investigation, Mapping, and Prediction (CLIMAP) data. These species counts are linked to Levitus's [1982] modern water temperature data for the four caloric seasons, four depth ranges (0, 30, 50, and 75 m), and the combined means of those depth ranges. The relation between planktonic foraminiferal assemblages and sea surface temperature (SST) data is estimated using the newly developed SIMMAX technique, which is an acronym for a modern analog technique (MAT) with a similarity index, based on (1) the scalar product of the normalized faunal percentages and (2) a weighting procedure of the modern analog's SSTs according to the inverse geographical distances of the most similar samples. Compared to the classical CLIMAP transfer technique and conventional MAT techniques, SIMMAX provides a more confident reconstruction of paleo-SSTs (correlation coefficient is 0.994 for the caloric winter and 0.993 for caloric summer). The standard deviation of the residuals is 0.90°C for caloric winter and 0.96°C for caloric summer at 0-m water depth. The SST estimates reach optimum stability (standard deviation of the residuals is 0.88°C) at the average 0- to 75-m water depth. Our extensive database provides SST estimates over a range of -1.4 to 27.2°C for caloric winter and 0.4 to 28.6°C for caloric summer, allowing SST estimates which are especially valuable for the high-latitude Atlantic during glacial times.

Glacial North Atlantic: Sea surface conditions reconstructed by GLAMAP 2000

The response of the tropical ocean to global climate change and the extent of sea ice in the glacial nordic seas belong to the great controversies in paleoclimatology. Our new reconstruction of peak glacial sea surface temperatures (SSTs) in the Atlantic is based on census counts of planktic foraminifera, using the Maximum Similarity Technique Version 28 (SIMMAX-28) modern analog technique with 947 modern analog samples and 119 well-dated sediment cores. Our study compares two slightly different scenarios of the Last Glacial Maximum (LGM), the Environmental Processes of the Ice Age: Land, Oceans, Glaciers (EPILOG), and Glacial Atlantic Ocean Mapping (GLAMAP 2000) time slices. The comparison shows that the maximum LGM cooling in the Southern Hemisphere slightly preceeded that in the north. In both time slices sea ice was restricted to the north western margin of the nordic seas during glacial northern summer, while the central and eastern parts were ice-free. During northern glacial winter, sea ice advanced to the south of Iceland and Faeroe. In the central northern North Atlantic an anticyclonic gyre formed between 45° and 60°N, with a cool water mass centered west of Ireland, where glacial cooling reached a maximum of >12°C. In the subtropical ocean gyres the new reconstruction supports the glacial-to-interglacial stability of SST as shown by CLIMAP Project Members (CLIMAP) [1981]. The zonal belt of minimum SST seasonality between 2° and 6°N suggests that the LGM caloric equator occupied the same latitude as today. In contrast to the CLIMAP reconstruction, the glacial cooling of the tropical east Atlantic upwelling belt reached up to 6°-8°C during Northern Hemisphere summer. Differences between these SIMMAX-based and published U37[k]- and Mg/Ca-based equatorial SST records are ascribed to strong SST seasonalities and SST signals that were produced by different planktic species groups during different seasons.

Stable isotope record of sediment cores from the North Atlantic

Abrupt and short climate changes, such as the Younger Dryas, punctuated the last glacial-to-interglacial transition (Ruddiman and McIntyre, 1981 doi:10.1016/0031-0182(81)90097-3; Duplessy et al., 1981 doi:10.1016/0031-0182(81)90096-1; Oeschger et al. 1984; Broecker et al., 1985 doi:10.1038/315021a0). Broecker et al. (1988 doi:10.1029/PA003i001p00001) proposed that these may have been caused by an interruption of thermohaline circulation as inputs of glacial meltwater freshened the surface waters of the North Atlantic. The finding (Fairbanks, 1989 doi:10.1038/342637a0) that meltwater discharge was minimal during the Younger Dryas, however, led to the suggestion that the surface-water salinity drop might have been caused instead by changes in the freshwater budget (the difference between precipitation and evaporation), accompanied by a reduction in poleward advection of saline subtropical water. Here we use micropalaeontological and stable-isotope records from foraminifera in two cores from the North Atlantic to generate two continuous, high-resolution records of sea surface temperature and salinity changes over the past 18,000 years. Despite the injection of glacial meltwater during warm episodes, we find that sea surface salinity and temperature remain positively correlated during deglaciation. Cold, low-salinity events occurred during the early stages of deglaciation (14,500-13,000 years ago) and the Younger Dryas, but the minor injections of meltwater at high latitudes during these events are insufficient to account for the observed salinity changes. We conclude that an additional feedback from changes in the hydrological cycle and in advection was necessary to trigger changes in thermohaline circulation and thus in climate. This feedback did not act when the meltwater injection occurred at low latitude.

Age determination and stable isotope composition of sediment cores from the Arctic Ocean

On the basis of 52 sediment cores, analyzed and dated at high resolution, the paleoceanography and climate of the Last Glacial Maximum (LGM) were reconstructed in detail for the Fram Strait and the eastern and central Arctic Ocean. Sediment composition and stable isotope data suggest three distinct paleoenvironments: (1) a productive region in the eastern to central Fram Strait and along the northern Barents Sea continental margin characterized by Atlantic Water advection, frequent open water conditions, and occasional local meltwater supply and iceberg calving from the Barents Sea Ice Sheet; (2) an intermediate region in the southwestern Eurasian Basin (up to 84-85°N) and the western Fram Strait characterized by subsurface Atlantic Water advection and recirculation, a moderately high planktic productivity, and a perennial ice cover that breaks up only occasionally; and (3) a central Arctic region (north of 85°N in the Eurasian Basin) characterized by a low-salinity surface water layer and a thick ice cover that strongly reduces bioproduction and bulk sedimentation rates. Although the total inflow of Atlantic Water into the Arctic Ocean may have been reduced during the LGM, its impact on ice coverage and halocline structure in the Fram Strait and southwestern Eurasian Basin was strong.

(Table S1) Stable carbon and oxygen isotope ratios of benthic foraminifera of the western Atlantic Ocean

Oxygen and carbon isotopic data were produced on the benthic foraminiferal taxa Cibicidoides and Planulina from 25 new piston cores, gravity cores, and multicores from the Brazil margin. The cores span water depths from about 400 to 3000 m and intersect the major water masses in this region. These new data fill a critical gap in the South Atlantic Ocean and provide the motivation for updating the classic glacial western Atlantic d13C transect of Duplessy et al. (1988). The distribution of 13C of SumCO2 requires the presence of three distinct water masses in the glacial Atlantic Ocean: a shallow (~1000 m), southern source water mass with an end-member d13C value of about 0.3-0.5 per mil VPDB, a middepth (~1500 m), northern source water mass with an end-member value of about 1.5 per mil, and a deep (>2000 m), southern source water with an end-member value of less than -0.2 per mil, and perhaps as low as the -0.9 per mil values observed in the South Atlantic sector of the Southern Ocean (Ninnemann and Charles, 2002, doi:10.1016/S0012-821X(02)00708-2). The origins of the water masses are supported by the meridional gradients in benthic foraminiferal d18O. A revised glacial section of deep water d13C documents the positions and gradients among these end-member intermediate and deep water masses. The large property gradients in the presence of strong vertical mixing can only be maintained by a vigorous overturning circulation.