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.

Sedimentation rates and geochemical analysis of surface sediments from the South Atlantic

Extensive investigations of sedimentary barium were performed in the southern South Atlantic in order to assess the reliability of the barium signal in Antarctic sediments as a proxy for paleoproductivity. Maximum accumulation rates of excess barium were calculated for the Antarctic zone south of the polar front where silica accumulates at high rates. The correspondence between barium and opal supports the applicability of barium as a proxy for productivity. Within the Antarctic zone north of today's average sea ice maximum, interglacial vertical rain rates of excess barium are high, with a maximum occurring during the last deglaciation and early Holocene and during oxygen isotope chronozone 5.5. During these periods, the maximum silica accumulation was supposedly located south of the polar front. Glacial paleoproductivity, instead, was low within the Antarctic zone. North of the polar front, significantly higher barium accumulation occurs during glacial times. The vertical rain rates, however, are as high as in the glacial Antarctic zone. Therefore there was no evidence for an increased productivity in the glacial Southern Ocean.

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.

Estimation of sea surface temperatures and salinity of ODP Site 177-1089 (Appendix A)

ODP Site 1089 is optimally located in order to monitor the occurrence of maxima in Agulhas heat and salt spillage from the Indian to the Atlantic Ocean. Radiolarian-based paleotemperature transfer functions allowed to reconstruct the climatic history for the last 450 kyr at this location. A warm sea surface temperature anomaly during Marine Isotope Stage (MIS) 10 was recognized and traced to other oceanic records along the surface branch of the global thermohaline (THC) circulation system, and is particularly marked at locations where a strong interaction between oceanic and atmospheric overturning cells and fronts occurs. This anomaly is absent in the Vostok ice core deuterium, and in oceanic records from the Antarctic Zone. However, it is present in the deuterium excess record from the Vostok ice core, interpreted as reflecting the temperature at the moisture source site for the snow precipitated at Vostok Station. As atmospheric models predict a subtropical Indian source for such moisture, this provides the necessary teleconnection between East Antarctica and ODP Site 1089, as the subtropical Indian is also the source area of the Agulhas Current, the main climate agent at our study location. The presence of the MIS 10 anomaly in the delta13C foraminiferal records from the same core supports its connection to oceanic mechanisms, linking stronger Agulhas spillover intensity to increased productivity in the study area. We suggest, in analogy to modern oceanographic observations, this to be a consequence of a shallow nutricline, induced by eddy mixing and baroclinic tide generation, which are in turn connected to the flow geometry, and intensity, of the Agulhas Current as it flows past the Agulhas Bank. We interpret the intensified inflow of Agulhas Current to the South Atlantic as responding to the switch between lower and higher amplitude in the insolation forcing in the Agulhas Current source area. This would result in higher SSTs in the Cape Basin during the glacial MIS 10, due to the release into the South Atlantic of the heat previously accumulating in the subtropical and equatorial Indian and Pacific Ocean. If our explanation for the MIS 10 anomaly in terms of an insolation variability switch is correct, we might expect that a future Agulhas SSST anomaly event will further delay the onset of next glacial age. In fact, the insolation forcing conditions for the Holocene (the current interglacial) are very similar to those present during MIS 11 (the interglacial preceding MIS 10), as both periods are characterized by a low insolation variability for the Agulhas Current source area. Natural climatic variability will force the Earth system in the same direction as the anthropogenic global warming trend, and will thus lead to even warmer than expected global temperatures in the near future.

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.

Global data compilation of benthic data sets II

In this study we present a global distribution pattern and budget of the minimum flux of particulate organic carbon to the sea floor (J POC alpha). The estimations are based on regionally specific correlations between the diffusive oxygen flux across the sediment-water interface, the total organic carbon content in surface sediments, and the oxygen concentration in bottom waters. For this, we modified the principal equation of Cai and Reimers [1995] as a basic monod reaction rate, applied within 11 regions where in situ measurements of diffusive oxygen uptake exist. By application of the resulting transfer functions to other regions with similar sedimentary conditions and areal interpolation, we calculated a minimum global budget of particulate organic carbon that actually reaches the sea floor of ~0.5 GtC yr**-1 (>1000 m water depth (wd)), whereas approximately 0.002-0.12 GtC yr**-1 is buried in the sediments (0.01-0.4% of surface primary production). Despite the fact that our global budget is in good agreement with previous studies, we found conspicuous differences among the distribution patterns of primary production, calculations based on particle trap collections of the POC flux, and J POC alpha of this study. These deviations, especially located at the southeastern and southwestern Atlantic Ocean, the Greenland and Norwegian Sea and the entire equatorial Pacific Ocean, strongly indicate a considerable influence of lateral particle transport on the vertical link between surface waters and underlying sediments. This observation is supported by sediment trap data. Furthermore, local differences in the availability and quality of the organic matter as well as different transport mechanisms through the water column are discussed.

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.

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).

Organic carbon accumulation in the South Atlantic Ocean

A compilation of 1118 surface sediment samples from the South Atlantic was used to map modern seafloor distribution of organic carbon content in this ocean basin. Using new data on Holocene sedimentation rates, we estimated the annual organic carbon accumulation in the pelagic realm (>3000 m water depth) to be approximately 1.8*10**12 g C/year. In the sediments underlying the divergence zone in the Eastern Equatorial Atlantic (EEA), only small amounts of organic carbon accumulate in spite of the high surface water productivity observed in that area. This implies that in the Eastern Equatorial Atlantic, organic carbon accumulation is strongly reduced by efficient degradation of organic matter prior to its burial. During the Last Glacial Maximum (LGM), accumulation of organic carbon was higher than during the mid-Holocene along the continental margins of Africa and South America (Brazil) as well as in the equatorial region. In the Eastern Equatorial Atlantic in particular, large relative differences between LGM and mid-Holocene accumulation rates are found. This is probably to a great extent due to better preservation of organic matter related to changes in bottom water circulation and not just a result of strongly enhanced export productivity during the glacial period. On average, a two- to three-fold increase in organic carbon accumulation during the LGM compared to mid-Holocene conditions can be deduced from our cores. However, for the deep-sea sediments this cannot be solely attributed to a glacial productivity increase, as changes in South Atlantic deep-water circulation seem to result in better organic carbon preservation during the LGM.

Multivariate statistical analysis on the kaolinite/chlorite ratios from 20 South Atlantic sediment cores

Multivariate statistical analysis on the kaolinite/chlorite ratios from 20 South Atlantic sediment cores allowed for the extraction of two processes controlling the fluctuations of the kaolinite/chlorite ratio during the last 130,000 yrs, (1) the relative strength of North Atlantic Deep Water (NADW) inflow into the South Atlantic Ocean and (2) the influx of aeolian sediments from the south African continent. The NADW fluctuation can be traced in the entire deep South Atlantic while the dust signal is restricted to the vicinity of South Africa. Our data indicate that NADW formation underwent significant changes in response to glacial/interglacial climate changes with enhanced export to the Southern Hemisphere during interglacials. The most pronounced phases with Enhanced South African Dust Export (ESADE) occurred during cold Marine Isotope Stage (MIS) 5d and across the Late Glacial/Holocene transition from 16 ka to 4 ka (MIS 2 to 1). This particular pattern is attributed to the interaction of Antarctic Sea Ice extent, the position of the westerlies and the South African monsoon system.

Stable oxygen isotope and Mg/Ca ratios on Globorotalia inflata

We present a species-specific Mg/Ca-calcification temperature calibration for Globorotalia inflata from a suite of 38 core top samples from the South Atlantic (from 8° to 49°S). G. inflata is a deep-dwelling planktonic foraminifer commonly occurring in subtropical to subpolar conditions, which qualifies it for reconstructions of the permanent thermocline. Apparent calcification depths and calcification temperatures were determined by comparing measured d18O with equilibrium d18O of calcite based on water column properties. Based on our core top samples, G. inflata apparent calcification depth is constant throughout the South Atlantic mid-latitudes with a depth of 350-400 m within the permanent thermocline. The resulting Mg/Ca-calcification temperature calibration is Mg/Ca = 0.72 +/-0.045/0.042 exp (0.076 +0.006 calcification 2 temperature) (r2 = 0.81) and covers the temperature range 3.1-16.5°C. We applied our Mg/Ca calibration to gravity core PS2495-3 from the Mid-Atlantic Ridge at ca. 41°S to test its validity by reconstructing a low-resolution record covering the last two glacial-interglacial cycles. Our paleotemperature record reveals large changes in temperature for Terminations I and II, when permanent thermocline temperature increased by as much as 8°C. The G. inflata paleotemperature record suggests that oceanic fronts repeatedly migrated over the location of site PS2495-3 during the last 160 kyr. This study shows the potential of G. inflata Mg/Ca to reconstruct paleotemperatures in the permanent thermocline.

Degradation rates and productivity signals of organic-walled dinoflagellate cyst species in surface sediment samples

To understand the role of the ocean within the global carbon cycle, detailed information is required on key-processes within the marine carbon cycle; bio-production in the upper ocean, export of the produced material to the deep ocean and the storage of carbon in oceanic sediments. Quantification of these processes requires the separation of signals of net primary production and the rate of organic matter decay as reflected in fossil sediments. This study examines the large differences in degradation rates of organic-walled dinoflagellate cyst species to separate these degradation and productivity signals. For this, accumulation rates of cyst species known to be resistant (R-cysts) or sensitive (S-cysts) to aerobic degradation of 62 sites are compared to mean annual chlorophyll-a, sea-surface temperature, sea-surface salinity, nitrate and phosphate concentrations of the upper waters and deep-water oxygen concentrations. Furthermore, the degradation of sensitive cysts, as expressed by the degradation constant k and reaction time t, has been related to bottom water [O2]. The studied sediments were taken from the Arabian Sea, north-western African Margin (North Atlantic), western-equatorial Atlantic Ocean/Caraibic, south-western African margin (South Atlantic) and Southern Ocean (Atlantic sector). Significant relationships are observed between (a) accumulation rates of R-cysts and upper water chlorophyll-a concentrations, (b) accumulation rates of S-cysts and bottom water [O2] and (c) degradation rates of S-cysts (kt) and bottom water [O2]. Relationships that are extremely weak or are clearly insignificant on all confidence intervals are between (1) S-cyst accumulation rates and chlorophyll-a concentrations, sea-surface temperature (SST), sea-surface salinity (SSS), phosphate concentrations (P) and nitrate concentrations (N), (2) between R-cyst accumulation rates and bottom water [O2], SST, SSS, P and N, and between (3) kt and water depth. Co-variance is present between the parameters N and P, N, P and chlorophyll-a, oxygen and water depth. Correcting for this co-variance does not influence the significance of the relationship given above. The possible applicability of dinoflagellate cyst degradation to estimate past net primary production and deep ocean ventilation is discussed.