Late Pliocene sea-surface temperature record of DSDP Hole 80-548 from the Northeast Atlantic Ocean (Table 1)

Application of quantitative micropaleontologic methods of sea-surface temperature (SST) estimation to the planktic foraminiferal faunas at Hole 548 in the northeast Atlantic between 3.5 and 2.0 Ma provide a detailed record of Pliocene surface temperature changes during the onset of major Northern Hemisphere glaciation. SST estimates indicate a relatively warm interval of "enhanced interglacials (warm peaks)" between 3.2 and 2.8 Ma followed by steady cooling with a significant step-like change toward colder temperatures occurring about 2.4 Ma. This record compares well with faunal and SST records from other North Atlantic sites.

(Table 1) Calculated sea surface temperature of sediment core RC11-210

Proxy indicators of sea surface temperature and equatorial divergence based on radiolarian assemblage data, and of trade wind intensity based on eolian grain size data show similar aspects of variability during the late Pleistocene: All indicators fluctuate at higher frequencies than the 100,000-year glacial-interglacial cycle, display reduced amplitude variations since 300,000 years ago, exhibit a change in the record character at about 300,000 years ago (the mid-Brunhes climatic event), and have higher amplitude variations in sediments 300,000-850,000 years old. Time series analyses were conducted to determine the spectral character of each record (delta18O of planktonic foraminifer, sea surface temperature values, equatorial divergence indicators, and wind intensity indicators) and to quantify interrecord coherence and phase relationships. The record was divided at the 300,000-year clear change in climatic variability (nonstationarity). The delta18O-based time scale is better lower in the core so our spectral analyses concentrated on the interval from 402,000-774,000 years. The delta18O spectra show 100,000- and 41,000-year power in the younger portion, 0-300,000 years, and 100,000-, 41,000- and 23,000-year power in the older interval, all highly coherent and in phase with the SPECMAP average stacked isotope record. Unlike the isotope record the dominant period in both the eolian grain size and equatorial divergence indicators is 31,000 years. This period is also important in the sea surface temperature signal where the dominant spectral peak is 100,000 years. The 31,000-year spectral component is coherent and in phase between the eolian and divergence records, confirming the link between atmospheric and ocean surface circulation for the first time in the paleoclimate record. Since the 31,000-year power appears in independent data sets within this core and also appears in other equatorial records [J. Imbrie personal communication, 1987], we assume it to be real and representative of both a nonlinear response to orbital forcing, possibly a combination of orbital tilt and eccentricity, and some resonance phenomenon required to amplify the response at this period so that it appears as a dominant frequency component. The mid-Brunhes climatic event is an important aspect of these records, but its cause remains unknown.

Occurrence of planktic foraminifera in Pliocene to Holocene sediments of DSDP Hole 81-552A in the North Atlantic (Appendix 2)

A new planktic foraminifer transfer function (GSF18) related 5 North Atlantic assemblages to winter and summer sea surface temperature. GSF18, based on recombined and simplified core top census data, preserves most environmental information and reproduces modern North Atlantic conditions with approximately the same accuracy as previous transfer functions, but can be more readily applied to faunal samples ranging in age from Pliocene to Holocene. Transfer function GSF18 has been applied to faunal data from Deep Sea Drilling Project Hole 552A to produce a 2.5 m.y. sea-surface temperature (SST) time series. Estimates show several periods between 2.3 and 4.6 Ma during which mean SST's were both several degrees warmer and several degrees cooler than modern conditions. Between 2.9 and 4.0 Ma SST was generally warmer than modern except for a 250 k.y. interval centered at 3.3 Ma. Maximum SST, with respect to modern conditions, occurred after the cool interval near 3.1 Ma when SST was approximately 3.6° C warmer than present conditions. Comparison of SST estimates with stable isotope data suggest that after peak warming at 3.1 Ma, there was an overall surface water cooling with concomitant build up of global ice volume, culminating in Northern Hemisphere glaciation. This event is also indicated by the presence of ice rafted detritus in 552A sediments at about 2.45 Ma.

Paleotemperature reconstructions of the Atlantic Sector of the Southern Ocean

Diatom assemblages from ODP Leg 177 sites 1093, 1094 and core PS2089-2, from the present Antarctic sea ice free zone and close to the Polar Front, were analyzed in order to reconstruct the climate development around the Mid-Brunhes Event 400 000 yr ago, as reflected by summer sea surface temperature (SSST) and sea ice distribution. Dense sample spacing allows a mean temporal resolution during Marine Isotope Stage (MIS) 11 (423-362 ka) of 300-400 yr. SSST values were estimated from diatom assemblages using a transfer function technique. The distribution pattern of sea ice diatoms indicates that the present-day ice free Antarctic Zone was seasonally covered by sea ice during the cold MIS 12 and MIS 10. These glacial periods are characterized by sea ice fluctuations with a periodicity of 3 and 1.85 kyr, suggesting the occurrence of Dansgaard-Oeschger-style millennial-scale oscillations in the Atlantic sector of the Southern Ocean during the glacial stages MIS 12 and MIS 10. Termination V (MIS 12/11) is characterized by a distinct temperature increase of 4-6°C, intersected especially at the southern site 1094 and core PS2089-2 by two distinct cooling events reminiscent of the Younger Dryas, which are associated with a northward shift of the winter sea ice edge in the Antarctic Zone. The SSST record is characterized by distinct temperature intervals bounded by stepwise, rapid changes. Maximum temperatures were reached during Termination V and the early MIS 11, exceeding modern values by 2°C over a period of 8 kyr. This pattern indicates a very early response of the Southern Ocean to global climate on Milankovitch-driven climate variability. The SSST optimum is marked by millennial-scale temperature oscillations with an amplitude of ca. 1°C and periodicities of ca. 1.85 and 1.47 kyr, probably reflecting changes in the ocean circulation system. The SSSTs during the MIS 11 temperature optimum do not exceed values obtained from other interglacial optima such as the early periods of MIS 5 or MIS 1 from the Antarctic Zone. However, the total duration of the warmest period was distinctly longer than observed from other interglacials. The comparison of the South Atlantic climate record with a high-resolution record from ODP Leg 162, site 980from the North Atlantic shows a strong conformity in the climate development during the studied time interval.

(Table 1) Radiocarbon ages of ODP Hole 146-893A

High-resolution planktonic foraminiferal census data from Santa Barbara Basin (Ocean Drilling Program hole 893A) demonstrate major assemblage switches between 25 and 60 ka that were associated with Dansgaard-Oeschger cycles. Stadials dominated by Neogloboquadrina pachyderma (sinistral), and Globigerinoides glutinata suggest a strong subpolar California Current influence, while interstadials marked by abundant N. pachyderma (dextral) and G. bulloides indicate a relative increase in subtropical countercurrent influence. Modern analog technique and transfer function (F-20RSC) temperature reconstructions support d18O evidence of large rapid (70 years or less) sea surface temperature shifts (3° to 5°C) between stadials and interstadials. Changes in the vertical temperature gradient and water column structure (thermocline depth) are recorded by planktonic faunal oscillations suggest bimodal stability in the organization of North Pacific surface ocean circulation. Santa Barbara Basin surface water demonstrates the rapid response of the California Current System to reorganization of North Pacific atmospheric circulation during rapid climate change.

Foraminiferal faunal estimates of paleotemperature: Circumventing the no-analog problem yields cool ice age tropics

The sensitivity of the tropics to climate change, particularly the amplitude of glacial-to-interglacial changes in sea surface temperature (SST), is one of the great controversies in paleoclimatology. Here we reassess faunal estimates of ice age SSTs, focusing on the problem of no-analog planktonic foraminiferal assemblages in the equatorial oceans that confounds both classical transfer function and modern analog methods. A new calibration strategy developed here, which uses past variability of species to define robust faunal assemblages, solves the no-analog problem and reveals ice age cooling of 5° to 6°C in the equatorial current systems of the Atlantic and eastern Pacific Oceans. Classical transfer functions underestimated temperature changes in some areas of the tropical oceans because core-top assemblages misrepresented the ice age faunal assemblages. Our finding is consistent with some geochemical estimates and model predictions of greater ice age cooling in the tropics than was inferred by Climate: Long-Range Investigation, Mapping, and Prediction (CLIMAP) [1981] and thus may help to resolve a long-standing controversy. Our new foraminiferal transfer function suggests that such cooling was limited to the equatorial current systems, however, and supports CLIMAP's inference of stability of the subtropical gyre centers.

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.

Planktonic foraminiferal d18O and d13C values, derived sea-surface temperatures from d18O and selected planktonic foraminiferal species for core P69 of Hawke Bay off eastern North Island

Recent evidence suggests that the Subtropical Convergence (STC) zone east of New Zealand shifted little from its modern position along Chatham Rise during the last glaciation, and that offshore surface waters north of the STC zone cooled only slightly. However, at nearshore core site P69 (2195 m depth), 115 km off the east coast of North Island and ca 300 km north of the modern STC zone, planktonic foraminiferal species, transfer function data and stable oxygen and carbon isotope records suggest that surface waters were colder by up to 6°C during the late last glacial period compared to the Holocene, and included a strong upwelling signature. Presently site P69 is bathed by south-flowing subtropical waters in the East Cape Current. The nearshore western end of Chatham Rise supports a major bathymetric depression, the Mernoo Saddle, through which some exchange between northern subtropical and southern subantarctic water presently occurs. It is proposed that as a result of much intensified current flows south of the Rise during the last glaciation, a consequence of more compressed subantarctic water masses, lowered sea level, and an expanded and stronger Westerly Wind system, there was accelerated leakage northwards of both Australasian Subantarctic Water and upwelled Antarctic Intermediate Water over Mernoo Saddle in a modified and intensified Southland Current. The expanded cold water masses displaced the south-flowing warm East Cape Current off southeastern North Island, and offshore divergence was accompanied by wind-assisted upwelling of nutrient-rich waters in the vicinity of P69. A comparable kind of inshore cold water jetting possibly characterised most glacial periods since the latest Miocene, and may account for the occasional occurrence of subantarctic marine fossils in onland late Cenozoic deposits north of the STC zone, rather than invoking wholesale major oscillations of the oceanic STC itself.

Diatom abundance in surface sediments of the northern North Pacific

In order to map the modern distribution of diatoms and to establish a reliable reference data set for paleoenvironmental reconstruction in the northern North Pacific, a new data set including the relative abundance of diatom species preserved in a total of 422 surface sediments was generated, which covers a broad range of environmental variables characteristic of the subarctic North Pacific, the Sea of Okhotsk and the Bering Sea between 30° and 70°N. The biogeographic distribution patterns as well as the preferences in sea surface temperature of 38 diatom species and species groups are documented. A Q-mode factor analysis yields a three-factor model representing assemblages associated with the Arctic, Subarctic and Subtropical water mass, indicating a close relationship between the diatom composition and the sea surface temperatures. The relative abundance pattern of 38 diatom species and species groups was statistically compared with nine environmental variables, i.e. the summer sea surface temperature and salinity, annual surface nutrient concentration (nitrate, phosphate, silicate), summer and winter mixed layer depth and summer and winter sea ice concentrations. Canonical Correspondence Analysis (CCA) indicates 32 species and species groups have strong correspondence with the pattern of summer sea surface temperature. In addition, the total diatom flux data compiled from ten sediment traps reveal that the seasonal signals preserved in the surface sediments are mostly from spring through autumn. This close relationship between diatom composition and the summer sea surface temperature will be useful in deriving a transfer function in the subarctic North Pacific for the quantitative paleoceanographic and paleoenvironmental studies. The relative abundance of the sea-ice indicator diatoms Fragilariopsis cylindrus and F. oceanica of >20% in the diatom composition is used to represent the winter sea ice edge in the Bering Sea. The northern boundary of the distribution of F. doliolus in the open ocean is suggested to be an indicator of the Subarctic Front, while the abundance of Chaetoceros resting spores may indicate iron input from nearby continents and shelves and induced productivity events in the study area.

Planktonic foraminiferal faunas from surface-sediment samples from the South China Sea and the western Pacific

We present 30 new planktonic foraminiferal census data of surface sediment samples from the South China Sea, recovered between 630 and 2883 m water depth. These new data, together with the 131 earlier published data sets from the western Pacific, are used for calibrating the SIMMAX-28 transfer function to estimate past sea-surface temperatures. This regional SIMMAX method offers a slightly better understanding of the marginal sea conditions of the South China Sea than the linear transfer function FP-12E, which is based only on open-ocean data. However, both methods are biased toward the tropical temperature regime because of the very limited data from temperate to subpolar regions. The SIMMAX formula was applied to sediment core 17940 from the northeastern South China Sea, with sedimentation rates of 20-80 cm/ka. Results revealed nearly unchanged summer temperatures around 28°C for the last 30 ky, while winter temperatures varied between 19.5°C in the last glacial maximum and 26°C during the Holocene. During Termination 1A, the winter estimates show a Younger Dryas cooling by 3°C subsequent to a temperature optimum of 24°C during the Bölling=Alleröd. Estimates of winter temperature differences between 0 and 100 m water depth document the seasonal variations in the thickness of the mixed layer and provide a new proxy for estimating past changes in the strength of the winter monsoon.

LGM, planktic foraminifera (EPILOG version)

[1] We used planktic foraminiferal assemblages in 70 sediment cores from the tropical and subtropical South Atlantic Ocean (10°N-37°S) to estimate annual mean sea surface temperatures (SSTs) and seasonality for the Last Glacial Maximum with a modified version of the Imbrie-Kipp transfer function method (IKTF) that takes into account the abundance of rare but temperature sensitive species. In contrast to CLIMAP Project Members [1981], the reconstructed SSTs indicate cooler glacial SSTs in the entire tropical/subtropical South Atlantic with strongest cooling in the upwelling region off Namibia (7-10°C) and smallest cooling (1-2°C) in the western subtropical gyre. In the western Atlantic, our data support recent temperature estimates from other proxies. In the upwelling regions in the eastern Atlantic, our data conflict with SST reconstructions from alkenones, which may be due to an environmental preference of the alkenone-producing algae or to an underestimation of foraminiferal SSTs due to anomalous high abundances of N. pachyderma (sinistral).

Estimates of annual mean surface temperatures in the late Quaternary of the tropical Atlantic

At least two modes of glacial-interglacial climate change have existed within the tropical Atlantic Ocean during the last 20,000 years. The first mode (defined by cold glacial and warm interglacial conditions) occurred symmetrically north and south of the equator and dominated the eastern boundary currents and tropical upwelling areas. This pattern suggests that mode 1 is driven by a glacial modification of surface winds in both hemispheres. The second mode of oceanic climate change, defined by temperature extremes centered on the deglaciation, was hemispherically asymmetrical, with the northern tropical Atlantic relatively cold and the southern tropical Atlantic relatively warm during deglaciation. A likely cause for this pattern of variation is a reduction of the presently northward cross-equatorial heat flux during deglaciation. No single mechanism accounts for all the data. Potential contributors to oceanic climate changes are linkage to high-latitude climates, modification of monsoonal winds by ice sheet and/or insolation changes, atmospheric CO2 and greenhouse effects, indirect effects of glacial meltwater, and variations in thermohaline overturn of the oceans.

Sea-surface temperature and sea ice distribution of the Southern Ocean at the EPILOG Last Glacial Maximum

Based on the quantitative study of diatoms and radiolarians, summer sea-surface temperature (SSST) and sea ice distribution were estimated from 122 sediment core localities in the Atlantic, Indian and Pacific sectors of the Southern Ocean to reconstruct the last glacial environment at the EPILOG (19.5-16.0 ka or 23 000-19 000 cal yr. B.P.) time-slice. The statistical methods applied include the Imbrie and Kipp Method, the Modern Analog Technique and the General Additive Model. Summer SSTs reveal greater surface-water cooling than reconstructed by CLIMAP (Geol. Soc. Am. Map Chart. Ser. MC-36 (1981) 1), reaching a maximum (4-5 °C) in the present Subantarctic Zone of the Atlantic and Indian sector. The reconstruction of maximum winter sea ice (WSI) extent is in accordance with CLIMAP, showing an expansion of the WSI field by around 100% compared to the present. Although only limited information is available, the data clearly show that CLIMAP strongly overestimated the glacial summer sea ice extent. As a result of the northward expansion of Antarctic cold waters by 5-10° in latitude and a relatively small displacement of the Subtropical Front, thermal gradients were steepened during the last glacial in the northern zone of the Southern Ocean. Such reconstruction may, however, be inapposite for the Pacific sector. The few data available indicate reduced cooling in the southern Pacific and give suggestion for a non-uniform cooling of the glacial Southern Ocean.

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.

Paleotemperature and productivity reconstruction of sediment core D11957P on the Iberian Margin

Foraminifera counts and climatic assemblages from the Tore Seamount are used to approach the glacial and interglacial changes in temperature and productivity on the Iberian Margin over the last 225 kyr. Chronostratigraphy is based on Globigerinoides ruber and Globigerina bulloides oxygen isotopes and supported by foraminifera and carbonate stadial fluctuations. Foraminifera indicate cooling from late interglacial stage 5 to the beginning of Termination I (TI). Neogloboquadnna pachyderma-s reflects cold conditions during glacial stages 4-2. In contrast, glacial stage 6 is dominated by warmer N. pachyderma-d and dutertrei and a restricted arctic assemblage. Past sea surface temperatures confirm the general cooling, reaching 4.3°C (SIMMAX.28) during stage 2. Multiple productivity proxies such as organic carbon, productivity-related foraminifera, and delta13C constrain the changes observed. A productivity increase occurs after interglacial stage 5, enhanced from late glacial stage 3 to TI Present-day satellite-detected phytoplankton plumes off Portugal would have accounted in the past glacial stages for the general productivity increase over the Tore. On top of this, welldefined peaks of organic carbon and productivity-related foraminifera correspond with Heinrich events 1-4.