Estimation of sea surface temperature from sediment core PS1768-8

The quantitative diatom analysis of 218 surface sediment samples recovered in the Atlantic and western Indian sector of the Southern Ocean is used to define a base of reference data for paleotemperature estimations from diatom assemblages using the Imbrie and Kipp transfer function method. The criteria which justify the exclusion of samples and species out of the raw data set in order to define a reference database are outlined and discussed. Sensitivity tests with eight data sets were achieved evaluating the effects of overall dominance of single species, different methods of species abundance ranking, and no-analog conditions (e.g., Eucampia Antarctica) on the estimated paleotemperatures. The defined transfer functions were applied on a sediment core from the northern Antarctic zone. Overall dominance of Fragilariopsis kerguelensis in the diatom assemblages resulted in a close affinity between paleotemperature curve and relative abundance pattern of this species downcore. Logarithmic conversion of counting data applied with other ranking methods in order to compensate the dominance of F. kerguelensis revealed the best statistical results. A reliable diatom transfer function for future paleotemperature estimations is presented.

Diatom abundance in surface sediments of the Southern Ocean

Based on the quantitative analysis of diatom assemblages preserved in 274 surface sediment samples recovered in the Pacific, Atlantic and western Indian sectors of the Southern Ocean we have defined a new reference database for quantitative estimation of late-middle Pleistocene Antarctic sea ice fields using the transfer function technique. The Detrended Canonical Analysis (DCA) of the diatom data set points to a unimodal distribution of the diatom assemblages. Canonical Correspondence Analysis (CCA) indicates that winter sea ice (WSI) but also summer sea surface temperature (SSST) represent the most prominent environmental variables that control the spatial species distribution. To test the applicability of transfer functions for sea ice reconstruction in terms of concentration and occurrence probability we applied four different methods, the Imbrie and Kipp Method (IKM), the Modern Analog Technique (MAT), Weighted Averaging (WA), and Weighted Averaging Partial Least Squares (WAPLS), using logarithm-transformed diatom data and satellite-derived (1981-2010) sea ice data as a reference. The best performance for IKM results was obtained using a subset of 172 samples with 28 diatom taxa/taxa groups, quadratic regression and a three-factor model (IKM-D172/28/3q) resulting in root mean square errors of prediction (RMSEP) of 7.27% and 11.4% for WSI and summer sea ice (SSI) concentration, respectively. MAT estimates were calculated with different numbers of analogs (4, 6) using a 274-sample/28-taxa reference data set (MAT-D274/28/4an, -6an) resulting in RMSEP's ranging from 5.52% (4an) to 5.91% (6an) for WSI as well as 8.93% (4an) to 9.05% (6an) for SSI. WA and WAPLS performed less well with the D274 data set, compared to MAT, achieving WSI concentration RMSEP's of 9.91% with WA and 11.29% with WAPLS, recommending the use of IKM and MAT. The application of IKM and MAT to surface sediment data revealed strong relations to the satellite-derived winter and summer sea ice field. Sea ice reconstructions performed on an Atlantic- and a Pacific Southern Ocean sediment core, both documenting sea ice variability over the past 150,000 years (MIS 1 - MIS 6), resulted in similar glacial/interglacial trends of IKM and MAT-based sea-ice estimates. On the average, however, IKM estimates display smaller WSI and slightly higher SSI concentration and probability at lower variability in comparison with MAT. This pattern is a result of different estimation techniques with integration of WSI and SSI signals in one single factor assemblage by applying IKM and selecting specific single samples, thus keeping close to the original diatom database and included variability, by MAT. In contrast to the estimation of WSI, reconstructions of past SSI variability remains weaker. Combined with diatom-based estimates, the abundance and flux pattern of biogenic opal represents an additional indication for the WSI and SSI extent.

(Table S1) Diatom abundance according to ice cycles in ANDRILL AND1-1B drill core

The influence of Antarctica and the Southern Ocean on Late Pliocene global climate reconstructions has remained ambiguous due to a lack of well-dated Antarctic-proximal, paleoenvironmental records. Here we present ice sheet, sea-surface temperature, and sea ice reconstructions from the ANDRILL AND-1B sediment core recovered from beneath the Ross Ice Shelf. We provide evidence for a major expansion of an ice sheet in the Ross Sea that began at ~3.3 Ma, followed by a coastal sea surface temperature cooling of ~2.5°C, a stepwise expansion of sea ice, and polynya-style deep mixing in the Ross Sea between 3.3 and 2.5 Ma. The intensification of Antarctic cooling resulted in strengthened westerly winds and invigorated ocean circulation. The associated northward migration of Southern Ocean fronts has been linked with reduced Atlantic Meridional Overturning Circulation by restricting surface water connectivity between the ocean basins, with implications for heat transport to the high latitudes of the North Atlantic. While our results do not exclude low-latitude mechanisms as drivers for Pliocene cooling, they indicate an additional role played by southern high-latitude cooling during development of the bipolar world.

Relative abundance of diatom species in surface samples from the Southern Ocean, Sea ice species

Diatoms from 228 Southern Ocean core-top sediment samples were examined to determine the geographic distributions of 32 major diatom species/taxa preserved in the sediments of three zonally-distinct regions; Sea Ice, Open Ocean and the Tropical/Subtropical. In the first of three papers, 14 species/taxa occurring in the region where sea ice covers the ocean surface on an annual basis are geographically documented. Comparisons are drawn between the diatom abundances on the sea floor, sea ice parameters (annual duration and concentration in February and September) and February sea-surface temperature. Such parameters are commonly used in reconstructing past oceanographic conditions in the Sea Ice and Open Ocean zones. Analysis of the geographic patterns and sea-surface parameter correlations reveals species-specific distributions regulated primarily by sea ice coverage and sea-surface temperature, which support the use of diatom remains for the estimation of these past sea-surface environmental parameters. Comparison with reliable accounts of the 14 species from the sediments or plankton also provides the first glimpses into species-specific ecology and habitat linkages.

(Table 1) Diatom mats from ODP Leg 177 Sites in the Southern Ocean, Atlantic Sector

The interaction between biogenic silica export and burial, paleoceanography, diatom species succession and mats formation was examined based on relative abundances data of Plio/Pleistocene diatoms from six cores recovered during ODP Leg 177 on a transect across the Antarctic Circumpolar Current (ACC) in the Atlantic sector of the Southern Ocean. Fragilariopsis kerguelensis, Actinocyclus ingens and species of the genus Thalassiothrix were the main contributors to the diatom assemblages. Three main steps marked the development of the silica system in the Southern Ocean: Step 1 (at ca. 2.77 Ma), establishment of increased biogenic silica burial in the Antarctic Circumpolar Current area, following the large-scale oceanic reorganization connected to the increased northern hemisphere glaciation; Step 2 (at ca. 1.93 Ma), the Antarctic Polar Front becomes the main biogenic silica sink, diatom mats are widespread, and are also found slightly to the north and south of the APF; Step 3 (at ca. 0.63 Ma), with the strong drop in abundance (and later extinction at 0.38 Ma) of A. ingens and the rise to dominance of F. kerguelensis, the system enters a glacial-interglacial mode, with diatom mats occurring during interglacials at the APF and in the Antarctic Zone, but disappearing north of it.

Nd-Sr isotopes, diatom and biogenic opal content of Pliocene sediments from IODP Site 318-U1361

Warm intervals within the Pliocene epoch (5.33-2.58 million years ago) were characterized by global temperatures comparable to those predicted for the end of this century (Haywood and Valdes, doi:10.1016/S0012-821X(03)00685-X) and atmospheric CO2 concentrations similar to today (Seki et al., 2010, doi:10.1016/j.epsl.2010.01.037; Bartoli et al., 2011, doi:10.1029/2010PA002055; Pagani et al., 2010, doi:10.1038/ngeo724). Estimates for global sea level highstands during these times (Miller et al., 2012, doi:10.1130/G32869.1) imply possible retreat of the East Antarctic ice sheet, but ice-proximal evidence from the Antarctic margin is scarce. Here we present new data from Pliocene marine sediments recovered offshore of Adélie Land, East Antarctica, that reveal dynamic behaviour of the East Antarctic ice sheet in the vicinity of the low-lying Wilkes Subglacial Basin during times of past climatic warmth. Sedimentary sequences deposited between 5.3 and 3.3 million years ago indicate increases in Southern Ocean surface water productivity, associated with elevated circum-Antarctic temperatures. The geochemical provenance of detrital material deposited during these warm intervals suggests active erosion of continental bedrock from within the Wilkes Subglacial Basin, an area today buried beneath the East Antarctic ice sheet. We interpret this erosion to be associated with retreat of the ice sheet margin several hundreds of kilometres inland and conclude that the East Antarctic ice sheet was sensitive to climatic warmth during the Pliocene.