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.

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.

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.

(Table S1) Age determination of Atlantic Ocean sediment cores

Well-dated benthic foraminifer oxygen isotopic records (d18O) from different water depths and locations within the Atlantic Ocean exhibit distinct patterns and significant differences in timing over the last deglaciation. This has two implications: on the one hand, it confirms that benthic d18O cannot be used as a global correlation tool with millennial-scale precision, but on the other hand, the combination of benthic isotopic records with independent dating provides a wealth of information on past circulation changes. Comparing new South Atlantic benthic isotopic data with published benthic isotopic records, we show that (1) circulation changes first affected benthic d18O in the 1000-2200 m range, with marked decreases in benthic d18O taking place at ~17.5 cal. kyr B.P. (ka) due to the southward propagation of brine waters generated in the Nordic Seas during Heinrich Stadial 1 (HS1) cold period; (2) the arrival of d18O-depleted deglacial meltwater took place later at deeper North Atlantic sites; (3) hydrographic changes recorded in North Atlantic cores below 3000 m during HS1 do not correspond to simple alternations between northern- and southern-sourced water but likely reflect instead the incursion of brine-generated deep water of northern as well as southern origin; and (4) South Atlantic waters at ~44°S and ~3800 m depth remained isolated from better-ventilated northern-sourced water masses until after the resumption of North Atlantic Deep Water (NADW) formation at the onset of the Bølling-Allerod, which led to the propagation of NADW into the South Atlantic.

Relative abundance of dinoflagellate cysts in surface sediments and dinoflagellate counts from three sediment cores from the Atlantic sector of the Southern Ocean

In this study we investigate the potential of organic-walled dinoflagellate cysts (dinocysts) as tools for quantifying past sea-surface temperatures (SST) in the Southern Ocean. For this purpose, a dinocyst reference dataset has been formed, based on 138 surface sediment samples from different circum-Antarctic environments. The dinocyst assemblages of these samples are composed of phototrophic (gonyaulacoid) and heterotrophic (protoperidinioid) species that provide a broad spectrum of palaeoenvironmental information. The relationship between the environmental parameters in the upper water column and the dinocyst distribution patterns of individual species has been established using the statistical method of Canonical Correspondence Analysis (CCA). Among the variables tested, summer SST appeared to correspond to the maximum variance represented in the dataset. To establish quantitative summer SST reconstructions, a Modern Analogue Technique (MAT) has been performed on data from three Late Quaternary dinocyst records recovered from locations adjacent to prominent oceanic fronts in the Atlantic sector of the Southern Ocean. These dinocyst time series exhibit periodic changes in the dinocyst assemblage during the last two glacial/interglacial-cycles. During glacial conditions the relative abundance of protoperidinioid cysts was highest, whereas interglacial conditions are characterised by generally lower cyst concentrations and increased relative abundance of gonyaulacoid cysts. The MAT palaeotemperature estimates show trends in summer SST changes following the global oxygen isotope signal and a strong correlation with past temperatures of the last 140,000 years based on other proxies. However, by comparing the dinocyst results to quantitative estimates of summer SSTs based on diatoms, radiolarians and foraminifer-derived stable isotope records it can be shown that in several core intervals the dinocyst-based summer SSTs appeared to be extremely high. In these intervals the dinocyst record seems to be highly influenced by selective degradation, leading to unusual temperature ranges and to unrealistic palaeotemperatures. We used the selective degradation index (kt-index) to determine those intervals that have been biased by selective degradation in order to correct the palaeotemperature estimates. We show that after correction the dinocyst based SSTs correspond reasonably well with other palaeotemperature estimates for this region, supporting the great potential of dinoflagellate cysts as a basis for quantitative palaeoenvironmental studies.

Morphology of Cycladophora davisiana davisiana from the Atlantic Ocean

In recent years, temporal fluctuations in the abundance of C. d. davisiana have been used frequently as a highresolution stratigraphic and paleoenvironmental tool. The modern ecology and morphologic variation (temporal and geographic) of this radiolarian species is evaluated to ascertain its potential stratigraphic and paleoenvironmental significance. Statistics were obtained on the width and height of all C. d. davisiana segments from Pleistocene populations of differing ages from the Northern Hemisphere (Labrador Sea and Iceland-Faeroe Ridge) and Southern Hemisphere (Namibian shelf and Meteor Rise). Results reveal that segment height variations between and within populations are more conservative than segment width. The mean sizes of the thorax and first abdominal segment have distinguishable differences between C. d. davisiana found in the North and South Atlantic. All populations have no significant difference in first abdominal segment width, however, mean heights of this segment differ greatly between populations of the North and South Atlantic. Second abdominal segment sizes show no clear population grouping. Size differences in post-cephalic segment size of these populations would appear to be related to some isolation of gene pools and possibly unknown paleoenvironmental factors. Temporal changes in the postcephalic size of C. d. davisiana may be used to: (1) identify temporally equivalent peaks in abundance of the species in a given region, (2) possibly evaluate the degree of mixing of water'masses between regions, and (3) trace the initial spread of the species from its area of origin. Cleve's 1887 plankton samples, between Greenland and Spitzsbergen, were studied and used in conjunction with other data to make the following conclusions on the modern ecology of C. d. davisiana in the Arctic and Greenland-Norwegian Seas. (1) It is presently absent in surface water plankton samples, (2) it currently lives at depths below 500 m, where it is rare, (3) it does not live in the upper 200 m under Arctic ice but is rare at greater depths, (4) it is absent in the upper 200 m near permanent Greenland Sea ice where normal oceanic salinity prevails, and (5) it is most common in deep marginal fjord environments which may serve as a refuge for the species during interglacial periods. In the Atlantic Ocean, the abundance of C. d. davisiana does not exceed 1% of the assemblage between the Subtropical Convergence of each hemisphere. In the Norwegian and Labrador Seas the species may occasionally be in the range of 1-5% of the modern radiolarian assemblage and never more than 5% in the southern high latitudes. Apparently only in the modern Sea of Okhotsk, does the species presently occur in high abundance. We concur with Morley and Hays (1983) that increased abundances are likely caused by the development of a strong low-salinity surface layer associated with seasonal sea ice melting and a strong temperature minimum above warmer and higher salinity intermediate waters. Similar conditions were frequent during the Pleistocene in the high latitudes and its modern scarcity outside the Sea of Okhotsk must be related to the absence of the presently unique conditions in the latter region.

Assessing consistency of fish survey data : uncertainties in the estimation of mackerel icefish (Champsocephalus gunnari) abundance at South Georgia

Acknowledgments The authors wish to thank the crews, fishermen and scientists who conducted the various surveys from which data were obtained, and Mark Belchier and Simeon Hill for their contributions. This work was supported by the Government of South Georgia and South Sandwich Islands. Additional logistical support provided by The South Atlantic Environmental Research Institute with thanks to Paul Brickle. Thanks to Stephen Smith of Fisheries and Oceans Canada (DFO) for help in constructing bootstrap confidence limits. Paul Fernandes receives funding from the MASTS pooling initiative (The Marine Alliance for Science and Technology for Scotland), and their support is gratefully acknowledged. MASTS is funded by the Scottish Funding Council (grant reference HR09011) and contributing institutions. We also wish to thank two anonymous referees for their helpful suggestions on earlier versions of this manuscript.

Nd isotopes an geochemistry of bulk deep sea sediments of the Atlantic Ocean

Nd isotopes preserved in fossil fish teeth and ferromanganese crusts have become a common tool for tracking variations in water mass composition and circulation through time. Studies of Nd isotopes extracted from Pleistocene to Holocene bulk sediments using hydroxylamine hydrochloride (HH) solution yield high resolution records of Nd isotopes that can be interpreted in terms of deep water circulation, but concerns about diagenesis and potential contamination of the seawater signal limit application of this technique to geologically young samples. In this study we demonstrate that Nd extracted from the > 63 µm, decarbonated fraction of older Ocean Drilling Program (ODP) sediments using a 0.02 M HH solution produces Nd isotopic ratios that are within error of values from cleaned fossil fish teeth collected from the same samples, indicating that the HH-extractions are robust recorders of deep sea Nd isotopes. This excellent correlation was achieved for 94 paired fish teeth and HH-extraction samples ranging in age from the Miocene to Cretaceous, distributed throughout the north, tropical and south Atlantic, and composed of a range of lithologies including carbonate-rich oozes/chalks and black shales. The strong Nd signal recovered from Cretaceous anoxic black shale sequences is unlikely to be associated with ferromanganese oxide coatings, but may be derived from abundant phosphatic fish teeth and debris or organic matter in these samples. In contrast to the deep water Nd isotopic signal, Sr isotopes from HH-extractions are often offset from seawater values, suggesting that evaluation of Sr isotopes is a conservative test for the integrity of Nd isotopes in the HH fraction. However, rare earth elements (REE) from the HH-extractions and fish teeth produce distinctive middle REE bulge patterns that may prove useful for evaluating whether the Nd isotopic signal represents uncontaminated seawater. Alternatively, a few paired HH-extraction and cleaned fish teeth samples from each site of interest can be used to verify the seawater composition of the HH-extractions. The similarity between isotopic values for the HH-extraction and fish teeth illustrates that the extensive cleaning protocol applied to fish teeth samples is not necessary in typical, carbonate-rich, deep sea sediments.

Stable oxygen isotope record of benthic foraminifera from the tropical Atlantic

Benthic foraminiferal oxygen isotope ratios from two sediment cores recovered at 426 and 1299 m water depth in the eastern and western tropical Atlantic show that a slowdown of the thermohaline circulation (THC) during Heinrich event H1 and the Younger Dryas was accompanied by rapid and intense warming of intermediate depth waters. Millennial-scale covariations of low paleosalinities in the subpolar North Atlantic with decreased benthic oxygen isotope ratios in the eastern tropical Atlantic throughout the past 10,000 years suggest that THC weakening might be related to middepth warming during the Holocene period as well. Climate model experiments simulating a strong reduction of the THC in the Atlantic Ocean under present-day and glacial conditions reveal that the increase of temperature in the middepth tropical and South Atlantic is a common feature for both climatic states, caused by a reduced ventilation of cold intermediate and deep waters in conjunction with downward mixing of heat from the thermocline. From the similarity of the paleoclimatic records with the model simulations, we infer that the characteristic pattern of temperature change in the Atlantic Ocean related to weakened thermohaline circulation can serve as an indicator of present-day and future THC slowdown.

Abundance of organic-walled dinoflagellate cysts in surface sediments from the Atlantic sector of the Southern Ocean (Appendix 2)

Thirty-two surface sediment samples from the Southern Ocean (eastern Atlantic sector), between the Subtropical Front and the Weddell Gyre, were investigated to provide information on the distribution of modern organic-walled dinoflagellate cysts in relation to the oceanic fronts of the Antarctic Circumpolar Current (ACC). A clearly distinguishable distribution pattern was observed in relation to the water masses and fronts of the ACC. The dinoflagellate cysts of species characteristic of open oceanic environments, such as Impagidinium species, are highly abundant around the Subtropical Front, whereas south of this front, cosmopolitan species such as Nematosphaeropsis labyrinthus and the cysts of Protoceratium reticulatum characterise the transition from subtropical to subantarctic surface waters. The subantarctic surface waters are dominated by the cysts of heterotrophic dinoflagellates, such as Protoperidinium spp. and Selenopemphix antarctica. The cysts of Protoperidinium spp. form the dominant part of the assemblages around the Antarctic Polar Front, whereas S. antarctica concentrations increase further to the south. The presence of S. antarctica in sediments of the Maud Rise, a region of seasonal sea-ice cover, reflects its tolerance for low temperatures and sea-ice cover. A previously undescribed species, Cryodinium meridianum gen. nov. sp. nov., has a restricted distribution pattern between the Antarctic Polar Front and the ACC-Weddell Gyre Boundary.

Carbon isotope composition of lipids from bones and host sediments of the Atlantic and Pacific Oceans

Carbon in lipids separated from organic matter of fish and marine mammal bones from bottom of the Pacific and Atlantic oceans has d13C values ranging from -21.6 to -25.8 per mil and is isotopically lighter than that in lipids and total organic matter of host sediments. During fossilization of organic phosphate carbon isotope composition of bound lipids of fish bone becomes lighter and that of bones of mammals becomes heavier, possibly as a result of metabolisms of these organisms and composition of phospholipids in them.