Response of planktic foraminiferal size to late Quaternary climate ...

17 downloads 0 Views 749KB Size Report
and interglacial (three out of stage 3, 5.1, 5.5, 7.3) samples ..... their relationship to orbital forcing recorded in ... response to orbital forcing in the late Quatern-.
PALEOCEANOGRAPHY, VOL. 18, NO. 2, 1039, doi:10.1029/2002PA000831, 2003

Response of planktic foraminiferal size to late Quaternary climate change Daniela N. Schmidt,1 Sabrina Renaud,2 and Jo¨rg Bollmann Geological Institute, ETH Zu¨rich, Zu¨rich, Switzerland Received 17 July 2002; revised 15 January 2003; accepted 11 February 2003; published 22 May 2003.

[1] Geographical size distribution within entire Holocene foraminiferal assemblages is related to global environmental gradients such as temperature, primary productivity, and environmental variability. This study demonstrates that these correlations are also recognizable in late Quaternary assemblages from three locations in the South Atlantic on temporal and latitudinal scales. The size response to temporal paleoenvironmental changes during glacial-interglacial cycles mimics the geographic Holocene size variability. The amplitude of size variability is directly related to the amplitude of the climatic fluctuations as shown by the stable size-temperature relationship over time. The documented changes in the assemblage size are caused by species replacement and intraspecific size variability. The relative importance of these processes depends on the environmental setting. Species have been shown to reach their maximum size and abundance under certain optimum conditions and decrease in size if environmental conditions differ from these optima. We confirm that late Quaternary species INDEX TERMS: 3030 Marine Geology and sizes were largest at paleotemperatures identical to Holocene ones. Geophysics: Micropaleontology; 4267 Oceanography: General: Paleoceanography; 4815 Oceanography: Biological and Chemical: Ecosystems, structure and dynamics; 4855 Oceanography: Biological and Chemical: Plankton; 9325 Information Related to Geographic Region: Atlantic Ocean; KEYWORDS: planktic foraminifers, morphometry, paleoecology, climate change, Quaternary Citation: Schmidt, D. N., S. Renaud, and J. Bollmann, Response of planktic foraminiferal size to late Quaternary climate change, Paleoceanography, 18(2), 1039, doi:10.1029/2002PA000831, 2003.

1. Introduction [2] Planktic foraminifers have, since their origin 140 million years ago, undergone three evolutionary radiations [Bolli et al., 1985] in both species richness and morphology [Cifelli, 1969; Norris, 1991]. Over this period, planktic foraminifers have grown from a negligible contribution to a major sink for pelagic carbonate. The excellent preservation, global occurrence and high abundance of planktic foraminifers within Cenozoic marine deep sea sediments are prime reasons for their extensive application in paleoceanographic and paleoclimatic studies. The analysis of species compositions [Imbrie and Kipp, 1971; Imbrie et al., 1973], stable isotopes [Shackleton and Opdyke, 1973; Hays et al., 1976], and other chemical tracers in the foraminiferal tests, e.g., Mg/ Ca [Nu¨rnberg et al., 1996; Elderfield and Ganssen, 2000], Cd/Ca ratios [Delaney and Boyle, 1987], and Boron isotopes [Sanyal et al., 1995] have proven to be valuable tools for paleoclimatic and paleoceanographic reconstructions. [3] In addition, morphological variability within individual taxa, such as size, shape and coiling direction of the test were shown to be correlated to fluctuations in temperature [Ericson, 1959; Be´ et al., 1973; Hecht et al., 1976; Malmg1 Now at Department of Geology, Royal Holloway, University of London, London, UK. 2 Now at Pale´oEnvironnements et Pale´oBiosphe`re, UMR S125 CNRS UCB Lyon 1, Villeurbanne, France.

Copyright 2003 by the American Geophysical Union. 0883-8305/03/2002PA000831$12.00

ren and Kennett, 1978a, 1978b; Naidu and Malmgren, 1996] or fertility [Naidu and Malmgren, 1995]. Different adult sizes are mainly influenced by environmentally controlled growth rates [Caron et al., 1987a, 1987b; Bijma et al., 1990a] since the lifespan is determined and reproduction triggered by the synodic lunar cycles [Spindler et al., 1979; Bijma et al., 1990a; Schiebel et al., 1997]. [4] The response of size and morphology to paleoceanographic changes has been variable, depending on species and location. Distinct size fluctuations related to late Quaternary climatic cycles seem to have occurred in the vicinity of frontal systems or in upwelling-areas [Be´ and Duplessy, 1976; Malmgren and Healy-Williams, 1978; Malmgren and Kennett, 1978b; Naidu and Malmgren, 1995]. In more stable environments observed stratigraphic size variability is seemingly uncorrelated to climatic cycles [Malmgren and Healy-Williams, 1978; Huber et al., 2000]. Despite successful application of size change for paleoceanographic reconstructions, these ambiguous correlations with environmental parameters have restricted the use of foraminiferal size variability as a global paleoclimatic index. [5] Most paleoceanographic and paleoclimatic reconstructions are based on the assumption that evolutionary adaptations of individual taxa do not change over the time periods considered, a hypothesis that has yet to be tested. Our size analyses of planktic foraminifera in globally distributed Holocene assemblages have allowed us to identify distinct size and abundance optima related to specific environmental conditions (D. N. Schmidt et al., Size distribution of Holocene planktic foraminifer assemblages:

17 - 1

17 - 2

SCHMIDT ET AL.: FORAMINIFERAL SIZE AND CLIMATE CHANGE

Figure 1. Core locations, main current systems, and fronts in the South Atlantic. Currents after Peterson and Stramma [1991]. Shaded area indicates depth less than 4000 m. Cores investigated for foraminiferal size data are indicated in bold; cores with reference data in italics.

biogeography, ecology and adaptation, submitted to Paleobiology, 2003; hereinafter referred to as Schmidt et al., submitted manuscript, 2003). [6] We therefore can test the stability of these relationships through the late Quaternary climatic cycles in the selected cores. The maximum size of entire planktic foraminifera assemblages were found to generally increase from the poles to the tropics parallel to surface temperatures. Decreased sizes relative to this general trend were observed in assemblages from upwelling areas and frontal systems. [7] On the basis of these new insights, we attempt here to test the temporal stability of the observed correlation of size and environmental gradients. Specifically, the present study of assemblage size variability in the Quaternary will allow us to test the following hypotheses: [8] 1. Foraminiferal size adaptation to specific environmental optima has been identified in the Holocene. Globally, Holocene planktic foraminiferal assemblage sizes depend on temperature, primary productivity, and environmental variability. By comparison with paleoproxy estimates in Pleistocene cores, the temporal stability of these species’ optima can be investigated. Quaternary climatic changes, e.g., glacial-interglacial cycles, should have led to temporal size fluctuations. [9] 2. On the basis of the linear correlations of size and temperature in the Holocene, we expect that the size response was proportional to the extent of physical environmental change in any particular area. [10] 3. Size changes of foraminiferal assemblages may be caused by changes in species composition and/or by intraspecific size variability. These two potential causes should be distinct and identifiable in specific biogeographic areas.

[11] To test these hypotheses, we focus our study on the last 300 kyr using cores from three different environmental settings in the South Atlantic (Figure 1): (1) the equatorial upwelling (core GeoB1105) with high-amplitude paleoproxy variability of temperature and productivity influencing a diverse assemblage, (2) the subtropical gyre (core GeoB1413) with minimal past variability, and (3) the subpolar frontal system (core PS2498) of intermediate variability and low species richness.

2. Environmental Settings and Faunal Composition 2.1. Modern Environmental Settings [12] The modern South Atlantic is dominated by a large subtropical gyre, which is bordered by the equatorial current system in the north, the subtropical front in the south, the Benguela Current in the east and the Brazil current in the west (Figure 1). The three cores studied represent a transect along the southern Mid-Atlantic Ridge well above the lysocline (Table 1), with core GeoB1105 in the equatorial upwelling zone, core GeoB1413 at the edge of the subtropical gyre and core PS2498 near the subantarctic front (Figure 1). Age models, sedimentology and paleoceanography of these cores were described by Wefer et al. [1996] and Mackensen et al. [2001]. [13] The equatorial upwelling core (Table 2), GeoB1105, is from a location which today is characterized by strong seasonal upwelling leading to a seasonal sea-surface temperature variation from 23.1C to 28.3C [Levitus et al., 1994]. It is characterized by strong seasonal changes in the depth of thermocline [Wolff et al., 1999] and a high primary

17 - 3

SCHMIDT ET AL.: FORAMINIFERAL SIZE AND CLIMATE CHANGE Table 1. Core Location, Water Depth, Time Interval Investigated, Average Sedimentation Rate and Reference for the Age Modelsa

Core

Latitude, 

Longitude, 

GeoB1105-4 GeoB1413-4 V22-174 PS2498-1 PS2082

139.9 S 1540.8 S 101S 4409.2 S 4315 S

1225.7 W 927.3 W 1252 W 1409.2 W 1145 E

Depth, m

Age, ka

3225 3789 2630 3783 4610

300 400

Average Sedimentation Rate, cm/ka

130

4.80b 1.07b 1.97c 8.94d 3.55d

a Cores investigated for foraminiferal size are in bold, and nearby cores with additional paleoproxy data in italic. b From Wefer et al. [1996]. c From Imbrie et al. [1984]. d From Mackensen et al. [2001].

productivity of 100 gCm 2a 1 [Antoine et al., 1996]. The subtropical gyre core, GeoB1413, is located in an area of more stable environmental conditions, with seasonal temperature changes between 21.2C to 24.8C, a deep thermocline, and low primary productivity of 45 –60 gCm 2a 1. Core PS2498 is from the subpolar frontal system, characterized by temperatures ranging from 12.3C to 7.2C, a shallow thermocline and intermediate primary productivity of 60 gCm 2a 1 (Table 2). 2.2. Late Quaternary Environmental Variability [14] Glacial paleotemperatures were colder at all three core locations. In the equatorial upwelling region (GeoB1105), glacial paleotemperatures using various proxies were estimated to have been 4 – 5C colder during summer and 7 – 8C colder during winter [Wefer et al., 1996], suggesting increased seasonality. In the subtropical gyre (GeoB1413), a temperature-decrease of 2.5C at Termination I was calculated (S. Mulitza, University of Bremen, personal communication, 2001). In the subpolar environment, the temperature difference between the Holocene and the last glacial maximum is estimated to have been 7C [Brathauer and Abelmann, 1999]. [15] At the equator (GeoB1105) upwelling was intensified during glacial periods, evidenced by increased seasonal temperature contrasts [Wefer et al., 1996], and higher total organic carbon contents [Schneider et al., 1996]. To infer paleoproductivity we used a transfer function based on the accumulation rate of organic matter (Corg) according to the method of Stein [1991]. On the basis of unpublished Corg, porosity and wet bulk density data from Peter Mu¨ller for

GeoB1105 http://www.pangaea.de), we calculate a glacial paleoproductivity of 160 gCm 2a 1, and 40 gCm 2a 1 for the interglacials 5 and 7. For subtropical core GeoB1413, no paleoproductivity estimates are available. However, the carbon isotopic composition of planktic foraminifers was interpreted as being indicative of nutrient depletion of thermocline waters during glacial times [Wefer et al., 1996]. Paleoproductivity at PS2498 was estimated as 61 gCm 2a 1 during interglacials and 110 gCm 2a 1 during the last glacial maximum by Mackensen et al. [2001], indicating a northward shift of the polar front during glacials. [16] In the cores, the different environmental conditions are reflected in distinct faunal compositions throughout the late Quaternary. The equatorial upwelling faunas of GeoB1105 are mixtures of tropical species, such as Globigerinoides ruber (white and pink), Globigerinoides sacculifer, Globorotalia menardii, and of species associated with upwelling, such as Globorotalia inflata, Globigerina bulloides, and Globigerinella calida. The subtropical gyre assemblages of core GeoB1413 are mainly composed of G. ruber with minor contributions of G. sacculifer, whereas the subpolar assemblages of core PS2498 are dominated by G. bulloides with minor proportions of Neogloboquadrina pachyderma (dextral and sinistral) and G. inflata. In all three cores the species contributing to the Holocene assemblages are also present during glacials, except G. menardii, although relative abundances vary.

3. Material and Methods 3.1. Material [17] Three deep sea cores from the South Atlantic were used in this study. All three cores are from water depths above the present lysocline [Broecker and Takahashi, 1978]. They contain undisturbed sediments and have well-calibrated age models based on oxygen isotope measurements [Wefer et al., 1996; Mackensen et al., 2001]. Very different sedimentation rates result in variable age coverage in the three cores (Table 1): i.e., 300 kyr for GeoB1105, 400 kyr for GeoB1413, and 130 kyr for PS2498. Samples for the analysis of test size of all specimens in an assemblage have been taken at average intervals of about 3.2 kyr in core GeoB1105 (98 samples), 4.7 kyr in core GeoB1413 (81 samples), and 2.3 kyr in core PS2498 (56 samples). Intraspecific size variability was investigated in two different data sets: First in three glacial and interglacial samples for the largest 40

Table 2. Environmental Characteristics for the Site Locations of Studied Coresa Core

SST, C

DTs, C

DTgl, C

Dof Th, m

GeoB1105

25.9

4.5

4 – 5 summer 7 – 8 winter

20 – 30 summer 75 winter

GeoB1413 PS2498

23.2 9.7

2.5 3.0

7–8 7

150 mm to allow comparison with the faunal composition data from Wefer et al. [1996] and Mix et al. [1999]. Additionally, this size fraction allows one to exclude juvenile forms and to focus on the adult stage within the ontogeny of the foraminifers. 3.2. Size Analysis [18] In the present study, three different estimators for size changes are considered: (1) a variable of the size distribution within the entire assemblage, (2) the maximum test size attained for single species and (3) the average size of the species Globorotalia truncatulinoides. 3.2.1. Assemblage Size [19] To characterize the size distribution within an entire assemblage, a complete split of each sample was analyzed, containing on average 2012 individuals (minimum 779, maximum 4414). The entire data set consists of 462961 measurements. Using an automated image analysis system [Bollmann et al., 2003], the maximum diameter of each randomly oriented foraminifer was measured since this size parameter is least affected by orientation. We have calculated the accuracy of diameter determination at the applied magnification of 160 to be 3.31 mm. In order to test the accuracy and the reproducibility of measurements, standard glass microbeads have been analyzed (mean diameter with standard deviation of 168 mm ± 7.4 mm, 331 mm ± 14.1 mm, and 655 mm ± 29.0 mm, Duke Scientific cooperation). Repeated preparations and measurements of these microbeads showed a good reproducibility and an accuracy (172 mm ± 5.1 mm; 330 mm ± 3.51 mm; 636 mm ± 24.2 mm) of mean diameter for all microbeads. The reproducibility of 10-times repeated measurements of a natural foraminiferal sample of the mean size of the distribution is ±4.1mm, that of the median is ±3.9 mm and that of the 95-percentile is ±9.71 mm, which represents an error of 1.5 to 2.2%. Confidence intervals for the mean range from 1.5 mm (minimum value in PS2498) to 6.8 mm (maximum value in GeoB1413). The use of the 95th percentile prevents the use of the confidence interval for statistical precision. [20] Because the minimum diameter is given by the sieve size (>150 mm), it represents an artificial cut-off of the natural size distribution and therefore has no biological significance. Since the distributions are highly skewed toward larger sizes, the most suitable descriptor is the measured size value separating the largest five percent of the assemblage from the smaller 95% (sizeassemblage5). Details of the measuring technique and the choice of descriptors to characterize the assemblages are by Schmidt et al. (submitted manuscript, 2003). 3.2.2. Maximum Size of Single Species [21] The results of the Holocene survey indicate that the maximum test size measured for individual foraminiferal

Table 3. Overview of Paleoenvironmental Proxies Considered in the Present Studya Label GeoB1105

GeoB1413 V22-174 PS2498 PS2082

Proxy carbon/oxygen isotopes of G. ruber, G. inflata, G. crassaformis, carbon/oxygen isotopes of F. wuellerstorfi species richness diversity temperature reconstruction (forams) alkenone temperature Corg, density, porosity carbon/oxygen isotopes for G. ruber, G. crassaformis SST reconstruction species richness carbon/oxygen isotopes of F. wuellerstorfi primary productivity temperature (radiolaria)

Reference Wefer et al. [1996] Bickert and Wefer [1996] Wefer et al. [1996] Wefer et al. [1996] Wefer et al. [1996] Schneider et al. [1996] Wefer et al. [1996]; P. Mu¨ller (unpublished data, 1999) Wefer et al. [1996] Imbrie et al. [1984] Mix et al. [1999] Mackensen et al. [2001] Mackensen et al. [2001] Brathauer and Abelmann [1999]

a Paleotemperature estimates are based on changes in the composition of radiolarian or foraminiferal faunas [Wefer et al., 1996; Brathauer and Abelmann, 1999; Mix et al., 1999] or on alkenone measurements [Schneider et al., 1996]. The mass accumulation of total organic carbon in the sediments (Corg), density, and porosity have been used to calculate paleoproductivity.

species is characteristic of its ecological optimum. To trace these optima back in time, we identified the largest size for all species represented among the 40 largest individuals encountered. This was done in three representative glacial and interglacial samples from each core. The environmental conditions leading to maximum size were determined by comparison of measured size and environmental proxies. 3.2.3. Average Size of Globorotalia Truncatulinoides [22] To test the impact of environmental change at the scale of a species, the size variability of G. truncatulinoides was analyzed. This taxon was selected because it is present in all cores throughout the investigated time span. Its abundance in the total foraminiferal assemblages ranged from 0.1 to 6%. Because of these low abundances, G. truncatulinoides measurements may provide independent evidence from that of the entire assemblage. A representative split of the fraction >150 mm was picked and oriented, including on average 40 specimens (minimum 8, maximum 76). The maximum diameter measured for the spiral side was chosen as the size estimator (sizetrunc), because it is most comparable to the diameters determined in unoriented tests of the other size measurements in this study (assemblage size and maximum size per species). 3.3. Paleoenvironmental Proxies [23] Quantification of paleoenvironmental change is largely based on the use of ‘‘proxy’’ evidence. Proxies are measurable characteristics of sediments and fossil biota, which, via transfer functions, can be used to reconstruct the paleoenvironment. To calibrate our observed foraminiferal assemblage size changes in the past, we compiled data from several published paleo-proxy analyses (Table 3). These include oxygen and carbon isotopes studies of various planktic and benthic foraminifers, temperature

SCHMIDT ET AL.: FORAMINIFERAL SIZE AND CLIMATE CHANGE

17 - 5

Figure 2. Plot of sizeassemblage5 changes compared to paleotemperature [Wefer et al., 1996; Brathauer and Abelmann, 1999; Mix et al., 1999] and paleoproductivity [Wefer et al., 1996; Mackensen et al., 2001; P. Mu¨ller, unpublished data, 1999] for the last 300 kyr at the three locations. Numbers represent the isotope stages, and gray bars indicate glacial stages.

and primary productivity reconstructions, and estimators of species richness. [24] In equatorial core GeoB1105, paleoenvironmental reconstructions include isotopic analysis of benthic and planktic foraminifera, paleotemperatures based on foraminiferal transfer functions, modern analogue technique, and alkenones. Paleoproductivity estimates are based on organic carbon contents. Additionally, changes in functional groups of foraminifers, such as upwelling indicators, tropical species, symbiont-bearing and nonsymbiont bearing ones (Table 3), have been calculated on the basis of faunal counts of Wefer et al. [1996]. [25] In subtropical gyre core GeoB1413 only the isotopic composition of planktic foraminiferal shells are available. This is used for correlation with nearby core V22-174 (Figure 1), from which paleotemperature estimates are available on the basis of faunal transfer functions [Imbrie et al., 1984, 1989; Mix et al., 1999]. [26] In subpolar core PS2498, paleoproductivity estimates based on organic carbon accumulation rate are taken from Mackensen et al. [2001]. Paleotemperature estimates were taken from core PS2082 (Figure 1) and are based on transfer functions compiled from radiolarian assemblages [Brathauer and Abelmann, 1999]. Since both cores are at similar latitudes within the circumpolar current and today’s environmental settings are very similar, we assume that this relationship prevailed back in time. 3.4. Statistical Analysis [27] The relationship of sizeassemblage5 and sizetrunc with environmental proxies was tested using a simple linear

regression. The size difference between glacial and interglacial measurements was examined with a Student’s t-test with the null hypothesis that mean glacial sizes are different from interglacial ones. The hypothesis was rejected based on the t-value and the probability (p).

4. Results and Interpretation 4.1. Assemblage Size Changes and Paleotemperature Variability [28] The smallest sizeassemblage5 values (300 –492 mm) are found in the subpolar core PS2498, the largest sizes (472– 568 mm) in the subtropical gyre core GeoB1413, and intermediate sizes (378 – 527 mm) in the equatorial core GeoB1105 (Figure 2). Linear correlation (Table 4) indicates a positive correlation of size and temperature on a global scale (Figure 3), when all three cores are considered together (r = 0.909, p < 0.001). A similar relationship (Figure 3) prevails on a local scale in the subpolar (r = 0.798, p < 0.001) and the equatorial upwelling samples (r = 0.629, p < 0.001). The largest range in size is observed in subpolar core PS2498, associated with a corresponding paleotemperature range from 14.1C to 5.8C. Equatorial core GeoB1105 is also characterized by large size fluctuations in accord with a paleotemperature variability from 24.6C to 13C. Little size change is found in subtropical gyre core GeoB1413 with associated paleotemperature changes of 5C derived from the nearby core V22174. Since size and temperature changes are minimal, no significant correlation between size and paleotemperature is found in core GeoB1413 in the subtropics. Despite this lack

17 - 6

SCHMIDT ET AL.: FORAMINIFERAL SIZE AND CLIMATE CHANGE

Table 4. Linear Correlation of Environmental Factors With Sizeassemblage5 for All Coresa Proxy

All Cores GeoB1105 GeoB1413 PS2498

Oceanography d18O F.wuellerstorfi d13C F.wuellerstorfi d18O G.ruber d13C G. ruber 18 d O G.inflata d13C G.inflata d18O G.crassaformis d13C G.crassaformis

0.577

SST warm MAT SST cold MAT SST warm TF SST cold TF SST warm SST cold SST radiolaria SST alkenone

0.741

0.55 0.552

0.205

0.359 0.186 0.142 0.134

0.114 0.101

0.399

0.18

0.716 0.798 0.446 0.655

0.163 0.017

Temperature

0.909 0.802 0.629 0.622

Productivity TOC accumulation

0.049

0.641

0.582

are observed in subpolar core PS2498 and in the equatorial upwelling core GeoB1105 (Figure 2), although no universal size-paleoproductivity correlation could be found within our data set (Table 4). This is consistent with the results of the Holocene study (Schmidt et al., submitted manuscript, 2003), where productivity leads to only minor departures from the global size-temperature trend. However, the effect of productivity on size seems to differ according to the local setting, as represented by the discrepancy in the slopes of the size-productivity regressions at the sites PS2498 and GeoB1105 (Figure 4). 4.3. Assemblage Size Changes and Species Composition [31] Changes in temperature, productivity and frontal dynamics, characterized by high turbulence, frequently appearing eddies [Beckmann et al., 1987], and storm events [Schiebel et al., 1995], can affect assemblage size by altering the faunal composition. We tested the effect of species replacement in the equatorial upwelling core GeoB1105, where changes of faunal composition related to a strengthening of the upwelling during glacials have previously been documented [Wefer et al., 1996]. Species replacement occurs within this core according to the intensity of the upwelling situation. An upwelling assemblage is characterized by increased abundances of Neogloboquadrina dutertrei, G. bulloides and N. pachyderma

Diversity S

0.608

0.092

a

Values indicated are correlation coefficients (r). Bold values indicate significant correlation (p < 0.05); values in italic in ‘‘all’’ are based on only two cores (SST cold for GeoB1105 and GeoB1413, primary productivity for GeoB1105 and PS2498). For references, see Table 3. Abbreviations are as follows: MAT, Modern analogue technique; TF, transfer functions refer to paleotemperature estimates based on faunal composition (see Wefer et al. [1996] for explanation); TOC, total organic carbon; S, species richness.

of correlation, the results of GeoB1413 clearly fit into the regional pattern (Figure 3). Therefore both the geographical and temporal size variability in the Quaternary are consistent with today’s biogeographic size changes (Schmidt et al., submitted manuscript, 2003). [29] Throughout the late Quaternary, the sizeassemblage5 variability in all three cores parallels fluctuations in paleotemperature estimates (Figure 2). The linear relationship between size and temperature throughout the last 300 kyr shows that the amplitude of size variability is directly related to the amplitude of climatic fluctuations. Our results are therefore in good agreement with the general sizetemperature relationship evident from the Holocene study (Schmidt et al., submitted manuscript, 2003). 4.2. Assemblage Size Changes and Paleoproductivity Variability [30] In Holocene assemblages, negative deviations from the general size-temperature relationship could be related to high primary productivity (Schmidt et al., submitted manuscript, 2003). In the late Quaternary, significant negative correlations between paleoproductivity and sizeassemblage5

Figure 3. Sizeassemblage5 (mm) compared to paleotemperature (C) of all samples with paleotemperature estimates (for references, see Figure 2). Lines represent regression lines. For correlation coefficients, see Table 4. Additionally, the mean sizeassemblage5 per biogeographic zone of the Holocene are indicated (D. N. Schmidt et al., submitted manuscript, 2003). Error bars indicate the standard deviations for temperature, respectively sizeassemblage5.

SCHMIDT ET AL.: FORAMINIFERAL SIZE AND CLIMATE CHANGE

Figure 4. Plot of paleoproductivity (gCm-2a-1) versus sizeassemblage5 (mm) for core GeoB1105 (gray squares) and core PS2498 (black circles). Primary productivity estimates for PS2498 from Mackensen et al. [2001] and for GeoB1105 based on data from P. Mu¨ller (unpublished data, 1999) and Wefer et al. [1996]. For correlation coefficients, see Table 4.

17 - 7

ronmental conditions [Hecht, 1976; Schmidt et al., submitted manuscript, 2003]. Maximum test size measured for individual foraminiferal species is characteristic for its ecological optimum. To trace these optima back in time, we identified the largest size for all species represented among the 40 largest individual encountered in representative glacial (three out of LGM, stages 4, 6 and 8, interstadial 5.2) and interglacial (three out of stage 3, 5.1, 5.5, 7.3) samples from each core. [33] The maximum size of the cosmopolitan species Globorotalia truncatulinoides, G. inflata, and G. crassaformis increased with increasing paleotemperature from 7C to 26C paleotemperature. This trend is parallel to the size increase in Holocene sediments (Figure 6a) and confirms the general temperature/size relationship found in Holocene sediments. In contrast, maximum sizes of Globigerinoides sacculifer, Globorotalia menardii and G. tumida, were found to be consistently smaller than in the Holocene (Figure 6b). This reduced maximum size might be caused by consistently lower Pleistocene than Holocene temperature estimates for the analyzed cores. Consequently largest Pleistocene test sizes are smaller than their optimum size (Figure 6b). Nevertheless the general species-specific

compared to a tropical assemblage, which is dominated by G. menardii, G. ruber, G. sacculifer, Globigerinella siphonifera, G. calida, Globigerina digitata, and Pulleniatina obliquiloculata [Hemleben et al., 1989]. On the basis of these differences, we used the ratio of tropical versus upwelling species as an indicator for species replacement. Accordingly, we also calculated the ratio of symbiont bearing (G. menardii, G. ruber, G. sacculifer, G. siphonifera, Orbulina universa, G. calida, Globigerina falconensis) versus the nonsymbiont bearing species (G. truncatulinoides, N. dutertrei, Globorotalia scitula, Globigerinoidestenellus, P. obliquiloculata, G. inflata, Globorotaloideshexagonus). We have not grouped any of the species with an unclear assignment to one of the groups [Hemleben et al., 1989; R. Schiebel, ETH, Zurich, personal communication, 2002]. The ratio of tropical versus upwelling species covaries (r = 0.761, p < 0.001) with assemblage sizes in core GeoB1105 (Figure 5), as does the ratio of symbiont bearing to nonsymbiont bearing ones (r = 0.815, p < 0.001). During times of intense upwelling, the relative abundance of tropical species or symbiont bearing ones is reduced by increasing numbers of upwelling indicators and nonsymbiont bearing species. The sizeassemblage5 is smaller during times with high abundances of upwelling species and nonsymbiont bearing ones. The observed size ranges in GeoB1105 (378– 527 mm) fluctuate between values characteristic for Holocene upwelling assemblages (451 mm) and those characteristic for subtropical/tropical assemblages (509/625 mm) (Schmidt et al., submitted manuscript, 2003). 4.4. Size Increases Within Species: Indicator of Optimum Ecological Conditions [32] Previous studies demonstrated that test sizes within individual species vary systematically with changing envi-

Figure 5. Plot of sizeassemblage5 (mm) (black circles) versus ratio of tropical to upwelling species (dark gray triangles) and symbiont bearing to nonsymbiont bearing species (gray rectangles) in core GeoB1105. Ratios based on data from Wefer et al. [1996].

17 - 8

SCHMIDT ET AL.: FORAMINIFERAL SIZE AND CLIMATE CHANGE

Figure 6. (a) Maximum size (mm) for G. truncatulinoides, G. inflata, and G. crassaformis. Open symbols, Holocene geographic data (Schmidt et al., submitted manuscript, 2003); filled symbols (late Quaternary based on isotope stages 2 to 8) versus paleotemperature. (b) Maximum sizes for investigated species in the Holocene compared to maximum sizes in the late Quaternary. Line indicates 1:1 relationship. (c) Optimum sea surface temperature in the Holocene versus optimum (paleo) temperatures in late Quaternary (isotope stages 2 to 8) for the investigated species. Ecological optimum is estimated as the temperature where each species reaches its largest size and maximum abundance (for references of temperature estimates, see Figure 2). To obtain consistency of temperature estimates from all three sites, optimum paleotemperatures are summer estimates, in contrast to winter estimates in Figure 2. Line indicates 1:1 relationship.

relationship of maximum test size and water temperature is constant over the analyzed time period (stages 2 to 8) and the data show a good general agreement (Figure 6c), indicating a stable relationship. [34] To increase our understanding of the influence of the intraspecific variability on sizeassemblage5, we studied in detail the reaction of G. truncatulinoides to the environmental change in the last 130 kyr. Since determinations of maximum size are affected by a random bias, we have picked on average 40 specimens and calculated the average size of the species per sample. Additionally, we increased the resolution in comparison to the maximum size data set of three glacial interglacial samples up to 32 samples a core. [35] However, the intraspecific size variability of G. truncatulinoides does not show clear glacial-interglacial size variability over the past 130 kyr in cores PS2498 and GeoB1413 (Student’s t-test: null hypothesis sizes are different for glacials and interglacials: p = 0.576, p = 0.422 respectively). None of the cores (Figure 7) record a significant correlation between sizetrunc and sizeassemblage5 (linear correlation GeoB1105: r = 0.610, p = 0.061; GeoB1413: r = 0.184, p = 0.359; PS2498: r = 0.062, p =

0.807). Either the correlations of size and paleotemperature may only prevail at large scales, or alternatively the biology and ecology of G. truncatulinoides are more complicated than those of other taxa.

5. Discussion 5.1. Stability of Abiotic Size Controls During the Late Quaternary [36] Our data strongly suggest that the Holocene relationship between planktic foraminiferal assemblage size and temperature, productivity and frontal dynamics (Schmidt et al., submitted manuscript, 2003) persisted over the past 300 kyr. As in the Holocene, the main trend emerging at both regional and local temporal scales is the positive correlation of size and temperature. Size-temperature relationships have been documented for various planktic foraminiferal species in the Holocene in laboratory cultures [Caron et al., 1987a, 1987b; Bijma et al., 1990b; Spero et al., 1991], plankton tows [Ortiz et al., 1995] and sediments [Kennett, 1968; Be´ et al., 1973; Hecht, 1976; Malmgren and Kennett, 1976]. This makes sense from a physiological standpoint since metabolic rates are

SCHMIDT ET AL.: FORAMINIFERAL SIZE AND CLIMATE CHANGE

17 - 9

Figure 7. Sizeassemblage5 (mm) and sizetrunc (mm) versus time for the last 130 kyr in PS2498, GeoB1413, and GeoB1105. Shaded area represents glacial stages.

known to increase with temperature [Caron et al., 1987a, 1987b; Bijma et al., 1990b; Spero et al., 1991]. Additionally, higher species richness at warm temperatures [Rutherford et al., 1999] would increase the likelihood of the presence of larger taxa in the assemblages (see Schmidt et al., submitted manuscript, 2003). [37] In the Holocene and the late Quaternary, secondary effects influencing size are productivity and variability of the environment, especially in upwelling areas and frontal systems. An excellent example of the local dominance of primary productivity changes on foraminiferal species’ sizes was presented by Naidu and Malmgren [1995] in the Arabian Sea. They recorded positive correlation between upwelling strength, indicated by the accumulation of G. bulloides, and size for several distinct species. However, within the paleoproductivity range of this study (30 – 160 gCm 2a 1), size and productivity display a negative correlation. Our own Holocene analyses have demonstrated that upwelling areas and frontal systems, apparently because of their high seasonal variability, prevent single species from attaining their size and abundance optima (Schmidt et

al., submitted manuscript, 2003), that they are able to reach in more stable areas. The dynamics of the frontal systems lead to expatriation [Berger, 1970; Weyl, 1978] and vertical displacement of biota. Fronts appear to function as environmental barriers [Schiebel et al., 2002] between stable ecosystems, each harboring groups of well-adapted species able to grow to large sizes. Based on this, the negative size-productivity relationships in equatorial upwelling core GeoB1105 and the subpolar front data set (PS2498) are most likely the result of increased environmental variability. On the basis of our Holocene data, we infer that increasing seasonality of temperatures during glacials lead to a size reduction despite higher paleoproductivities and consequently resulting in the observed negative correlation between paleoproductivity and size in two or our cores (Figure 4). 5.2. Biotic Processes: Species Replacement and Intraspecific Size Variability [38] In general, assemblage size differences are a combined result of modifications in faunal composition coupled

17 - 10

SCHMIDT ET AL.: FORAMINIFERAL SIZE AND CLIMATE CHANGE

with size changes within single species as a response to environmental variability. In addition, different species have specific size spectra [Hecht, 1976] and specific environmental optima. The size response to a change in a particular environmental parameter will show a positive correlation below and a negative correlation above the optimum conditions. For instance, Globigerina bulloides is known to have replaced subtropical and tropical species in the equatorial upwelling areas during glacial times. Since G. bulloides is generally smaller then the species it replaces, e.g., G. menardii, G. sacculifer, G. ruber, the total assemblage size decreases. [39] In addition to species replacement, there is evidence for size reduction of single species in upwelling areas [Ortiz et al., 1995; Schmidt et al., submitted manuscript, 2003]. Growth of foraminifera under nonoptimum conditions may be limited by spatial heterogeneity of, e.g., temperature, nutrients, or turbidity. Increasing turbidity was found to favor asymbiotic species [Ortiz et al., 1995] which are smaller than symbiont-bearing ones [Hecht, 1976]. In addition, symbiotic species grow to smaller sizes under light limited conditions [Be´ et al., 1982]. Both processes may be influential in our size record of equatorial core GeoB1105 (Figure 5). During intensive glacial upwelling the proportion of symbiont-bearing species decreases to 20% relative to 60% in interglacial periods, causing a decrease in sizeassemblage5. In contrast, in the subtropical gyre, the same species dominate cold and warm intervals in the record. Assemblage size changes in these cores are therefore not the result of species replacement and may be representative of intraspecific size variability. [40] We have attempted to characterize intraspecific size changes in G. truncatulinoides. The size of G. truncatulinoides increases among the cores in parallel with their general paleotemperature estimates (Figure 7). Despite these biogeographic patterns, the fine-scale stratigraphic size variability is not very closely correlated to paleotemperature estimates [Renaud and Schmidt, 2003]. The absence of the expected correlation underlines the ecological and taxonomic complexity of this species. First, G. truncatulinoides lives in deep waters down to 2000 m [Hemleben et al., 1985, 1989]. Paleotemperature estimates of surface waters may not be relevant to the past habitat changes experienced by this species [Lohmann and Schweitzer, 1990]. Second, genetic analyses by de Vargas et al. [2001] indicate that G. truncatulinoides consists of four different genetic species, three of which have been found segregated along a temperature gradient. Species replacement is therefore likely to occur. Hence the down core sizetrunc measurements may be a combined result of intraspecific size variability and replacement by species of differing sizes [Renaud and Schmidt, 2003]. 5.3. Differential Response to Environmental Change of Various Amplitudes [41] Our results show that the relative importance of species replacement and intraspecific size changes depends on the environmental setting. Warm water species have narrow ecological tolerances, resulting in temperature optima which are very close to each other (Figure 6c).

Relatively minor temperature changes will therefore lead to species replacement. In contrast, few species are adapted to cold waters, limiting the role of species replacement in polar and subpolar areas. Observed size changes are thus much more likely to be due to intraspecific size variability. Our cores are suitable for testing these hypotheses. 5.3.1. GeoB1105: Species Replacement in Equatorial Upwelling Systems [42] Temperatures were lower and primary productivity was higher during the glacial and stadial phases [Wefer et al., 1996] than today. The glacial/interglacial summertemperature variability ranged from 18 to 27C. Several tropical and subtropical species have their optimum in this temperature range. Most of these species show a strong decrease in size outside their temperature optima, which eliminates their influence on the size of the largest 5% of the assemblage. Possible intraspecific size variability is attenuated by the dramatic turnover in species composition as expressed by order of magnitude changes in the ratio of tropical versus upwelling species (Figure 5). Upwellingspecies (N. dutertrei, G. bulloides, and N. pachyderma) dominate the glacial assemblages while tropical species (e.g., G. sacculifer, G. ruber, G. menardii) prevail during interglacial periods (Figure 5). In cold intervals, with vigorous upwelling and decreased temperatures, large tropical symbiont bearing species disappear. A conflation of these three processes (decrease in the number of (large) symbiont-bearing species and of (large) tropical species, coupled with observed size decrease of the persisting symbiont-bearing species) contributes to smaller size during cold intervals. 5.3.2. GeoB1413: Environmental and Size Stability in the Central Gyre [43] The position of this core in the northern part of the subtropical gyre suggests relatively stable Quaternary environmental conditions as indicated by the limited range in sea surface temperatures [Imbrie et al., 1989] and relatively constant assemblage characteristics [Mix et al., 1999]. This signal may have been additionally dampened by the low sedimentation rates of GeoB1413 (1 cm/kyr) since bioturbation reduces the glacial-interglacial amplitudes. The restricted size variability cannot solely be attributed to the low sedimentations rates because even in long ranging glacials such as stage 6 no size reduction in comparison to the interglacials can be measured. Accordingly, the small variability of sizeassemblage5 through time is related to more stable environmental conditions in the subtropical gyre than in the two other environmental settings. 5.3.3. PS2498: Subpolar Frontal Dynamics and Intraspecific Size Variability [44] The polar front today is about 6 to the south of core PS2498 [Peterson and Whitworth, 1989] and shifted relative to the core location during glacials [Mackensen et al., 1994; Brathauer and Abelmann, 1999], when temperatures were lower, primary productivity was higher and frontal dynamics increased. Sizeassemblage5 closely matched these changes. The largest species encountered are G. bulloides, G. truncatulinoides and N. pachyderma, which all grow to similar sizes under low-temperature conditions.

SCHMIDT ET AL.: FORAMINIFERAL SIZE AND CLIMATE CHANGE

Thus intraspecific size variability seems to be more important than the size effects of species replacement.

6. Conclusions [45] We tested hypotheses about the processes controlling foraminiferal assemblage size variability in the late Quaternary with the following results: [46] 1. The size response to temporal environmental changes during glacial-interglacial cycles mimics the spatial Holocene size variability. In the late Quaternary, as in the Holocene, planktic foraminiferal assemblage size variability is primarily affected by paleotemperature. Paleoproductivity and frontal dynamics also seem to have affected assemblage sizes. [47] 2. The amplitude of size change matches the extent of the inferred paleoclimatic fluctuations. [48] 3. Changes in abiotic parameters also affect biotic responses such as species replacement and intraspecific size variability. The relative importance of these depends on the environmental setting of the core. In the subtropics and tropics, represented by GeoB1413 and GeoB1105 respectively, the species have narrow ecological niches and species replacement is most likely to occur. In the subpolar environment, species have wider environmental

17 - 11

preferences and grow to similar sizes, but they still display intraspecific size decrease if moved away from their ecological optimum. Therefore intraspecific size change without species replacement is more likely to occur in this environment. [49] The similarity of the results observed in both the Holocene and late Quaternary records provide evidence that the ecological preferences of the investigated species did not change much during the last 300 kyr. Therefore the assumption of an analogous situation for paleoceanographic reconstructions based on planktic foraminiferal temperature optima is supported by our results. Expanding our approach to a longer time record might permit the recognition of incipient nonanalogue situations, where the influence of evolutionary change would conceivably become the dominant factor.

[50] Acknowledgments. Barbara Donner, Andreas Mackensen, and Gerold Wefer generously provided material for this study. We would like to thank Ralf Schiebel for valuable taxonomic discussions as well as Gerrit Meinecke, Stefan Mulitza, and Hans Thierstein. Peter Mu¨ller generously provided unpublished data. This paper has benefited from the comments of Judy McKenzie. The research was funded by the Swiss National Fund Projects: 2053-053676 Plankton Ecology and 2100-063461 Taxonomy and Oceanic nutrient cycling, ecology and climate: from oceanic to marginal environments.

References Antoine, D., J.-M. Andre´, and A. Morel, Oceanic primary production, 2: Estimation at global scale from satellite (coastal zone color scanner) chlorophyll, Global Biogeochem. Cycles, 10, 57 – 69, 1996. Be´, A. W. H., and J. C. Duplessy, Subtropical convergence fluctuations and Quaternary climates in the middle latitudes of the Indian Ocean, Science, 194, 419 – 422, 1976. Be´, A. W. H., S. M. Harrison, and L. Lott, Orbulina universa d’Orbigny in the Indian Ocean, Micropaleontology, 19, 150 – 192, 1973. Be´, A. W. H., H. J. Spero, and O. R. Anderson, Effects of symbiont elimination and reinfection on the life processes of the planktonic foraminifer Globigerinoides sacculifer, Mar. Biol., 70, 73 – 86, 1982. Beckmann, W., A. Auras, and C. Hemleben, Cyclonic coldcore eddy in the eastern North Atlantic, III: Zooplankton, Mar. Ecol. Prog. Ser., 39, 165 – 173, 1987. Berger, W. H., Planktonic foraminifera selective solution and paleoclimatic interpretation, Deep Sea Res., 15(1), 31 – 43, 1968. Berger, W. H., Planktonic foraminifera: Differential production and expatriation off Baja California, Limnol. Oceanogr., 15, 183 – 204, 1970. Bickert, T., and G. Wefer, Late Quaternary deep water circulation in the South Atlantic: Reconstruction from carbonate dissolution and benthic stable isotopes, in The South Atlantic: Present and Past Circulation, edited by G. Wefer et al., pp. 599 – 620, SpringerVerlag, New York, 1996. Bijma, J., J. Erez, and C. Hemleben, Lunar and semi-lunar reproductive cycles in some spinose planktonic foraminifers, J. Foraminiferal Res., 20, 117 – 127, 1990a. Bijma, J., W. W. Faber, and C. Hemleben, Temperature and salinity limits for growth and sur-

vival of some planktonic foraminifers in laboratory cultures, J. Foraminiferal Res., 20, 95 – 116, 1990b. Bolli, H. M., J. B. Saunders, and K. Perch-Nielsen, Plankton Stratigraphy, Cambridge Univ. Press, New York, 1985. Bollmann, J., P. Quinn, M. Vela, B. Brabec, S. Brechner, D. N. Schmidt, R. Schiebel, and H. R. Thierstein, Automated analysis of microfossils, in Image Analysis, Sediments and Paleoenvironments, edited by P. Francus, Kluwer Acad., Norwell, Mass., in press, 2003. Brathauer, U., and A. Abelmann, Late Quaternary variations in sea surface temperatures and their relationship to orbital forcing recorded in the Southern Ocean (Atlantic sector), Paleoceanography, 14, 135 – 148, 1999. Broecker, W. S., and T. Takahashi, The relationship between lysocline depth and in situ carbonate ion concentration, Deep Sea Res. Part I, 25, 65 – 95, 1978. Caron, D. A., W. W. J. Faber, and A. W. H. Be´, Effects of temperature and salinity on the growth and survival of the planktonic foraminifer Globigerinoides sacculifer, J. Mar. Biol. Assoc. UK, 67, 323 – 342, 1987a. Caron, D. A., W. W. J. Faber, and A. W. H. Be´, Growth of the spinose planktonic foraminifer Orbulina universa in laboratory culture and the effect of temperature on life processes, J. Mar. Biol. Assoc. UK, 67, 343 – 358, 1987b. Cifelli, R., Radiation of Cenozoic planktonic foraminifera, Syst. Zool., 18, 154 – 168, 1969. de Vargas, C., S. Renaud, H. Hilbrecht, and J. Pawlowski, Pleistocene adaptive radiation in Globorotalia truncatulinoides: Genetic, morphologic and environmental evidence, Paleobiology, 27, 104 – 125, 2001. Delaney, M. L., and E. A. Boyle, Cd/Ca in late Miocene benthic foraminifera and changes in the global organic carbon budget, Nature, 330, 156 – 159, 1987.

Elderfield, H., and G. Ganssen, Past temperature and d18O of surface ocean waters inferred from foraminiferal Mg/Ca ratios, Nature, 405, 442 – 445, 2000. Ericson, D. B., Coiling direction of Globigerina pachyderma as a climatic index, Science, 130, 219 – 230, 1959. Hays, J. D., J. Imbrie, and N. J. Shackleton, Variations in the Earth’s orbit: Pacemaker of the ice ages, Science, 194, 1121 – 1132, 1976. Hecht, A. D., An ecologic model for test size variation recent planktonic foraminifera: Application to the fossil record, J. Foraminiferal Res., 6, 295 – 311, 1976. Hecht, A. D., A. W. Be´, and L. Lott, Ecologic and paleoclimatic implications of morphologic variation of Orbulina universa in the Indian Ocean, Science, 194, 422 – 424, 1976. Hemleben, C., M. Spindler, I. Breitinger, and W. G. Deuser, Field and laboratory studies on the ontogeny and ecology of some Globorotaliid species from the Sargasso Sea off Bermuda, J. Foraminiferal Res., 15, 254 – 272, 1985. Hemleben, C., M. Spindler, and O. R. Anderson, Modern Planktonic Foraminifera, 363 pp., Springer-Verlag, New York, 1989. Huber, R., H. Meggers, K.-H. Baumann, M. E. Raymo, and R. Henrich, Shell size variation of the planktonic foraminifer Neogloboquadrina pachyderma sin. in the Norwegian-Greenland Sea during the last 1.3 Myr: Implications for paleoceanographic reconstructions, Palaeogeogr. Palaeoclimatol. Palaeoecol., 160, 193 – 212, 2000. Imbrie, J., and N. Kipp, A new micropaleontological method for quantitative paleoclimatology: Application to a late Pleistocene Caribbean core, in The Late Cenozoic Glacial Ages, edited by K. K. Turekian, pp. 71 – 181, Yale Univ. Press, New York, 1971. Imbrie, J., J.-V. Donk, and N.-G. Kipp, Paleoclimatic investigation of a Late Pleistocene Car-

17 - 12

SCHMIDT ET AL.: FORAMINIFERAL SIZE AND CLIMATE CHANGE

ibbean deep-sea core; Comparison of isotopic and faunal methods, Quat. Res., 3, 10 – 38, 1973. Imbrie, J., J. D. Hays, D. G. Martinson, A. McIntyre, A. C. Mix, J. J. Morley, N. G. Pisias, W. L. Prell, and N. J. Shackleton, The orbital theory of Pleistocene climate: Support from a revised chronology of the marine d18O record, in Milankovitch and Climate, edited by A. L. Berger et al., pp. 269 – 305, D. Reidel, Norwell, Mass., 1984. Imbrie, J., A. McIntyre, and A. C. Mix, Oceanic response to orbital forcing in the late Quaternary: Observational and experimental strategy, in Climate and Geosciences, a Challenge for Science and Society in the 21st Century, edited by A. Berger, S. H. Schneider, and J. C. Duplessy, pp. 121 – 164, Kluwer Acad., Norwell, Mass., 1989. Kennett, J. P., Latitudinal variation in Globigerina pachyderma (Ehrenberg) in surface sediments of the southwest Pacific Ocean, Micropaleontology, 14, 305 – 318, 1968. Levitus, S., T. Boyer, R. Burgett, and M. Conkright, World Ocean Atlas 1994 and CD-ROM Data Sets, 604 pp., Natl. Oceanogr. Data Cent., Washington, D. C., 1994. Lohmann, G. P., and P. N. Schweitzer, Globorotalia truncatulinoides’ growth and chemistry as probes of the past thermocline, 1: Shell size, Paleoceanography, 5, 55 – 75, 1990. Mackensen, A., H. Grobe, H.-W. Hubberten, and G. Kuhn, Benthic foraminiferal assemblages and the d13C-signal in the Atlantic sector of the southern ocean: Glacial-to-interglacial contrasts, in Carbon Cycling in the Glacial Ocean: Constraints on the Ocean’s Role in Global Change, edited by R. Zahn et al., pp. 105 – 144, Springer-Verlag, New York, 1994. Mackensen, A., M. Rudolph, and G. Kuhn, Late Pleistocene deep-water circulation in the subantarctic eastern Atlantic, Global Planet. Change, 30, 197 – 229, 2001. Malmgren, B. A., and N. Healy-Williams, Variation in test diameter of Orbulina universa in the paleoclimatology of the late Quaternary of the Gulf of Mexico, Palaeogeogr. Palaeoclimatol. Palaeoecol, 25, 235 – 240, 1978. Malmgren, B. A., and J. P. Kennett, Biometric analysis of phenotypic variation in recent Globigerina bulloides d’Orbigny in the Southern Indian Ocean, Mar. Micropaleontol., 1, 3 – 25, 1976. Malmgren, B. A., and J. P. Kennett, Test size variation in Globigerina bulloides in response to Quaternary palaeoceanographic changes, Nature, 275, 123 – 124, 1978a. Malmgren, B. A., and J. P. Kennett, Late Quaternary paleoclimatic applications of mean size variations in Globigerina bulloides d’Orbigny

in the Southern Indian Ocean, J. Paleontol., 52, 1195 – 1207, 1978b. Mix, A. C., A. E. Morey, N. G. Pisias, and S. W. Hostetler, Foraminiferal faunal estimates of paleotemperature: Circumventing the no-analogue problem yields cool ice age tropics, Paleoceanography, 14, 350 – 359, 1999. Naidu, P. D., and B. A. Malmgren, Monsoon upwelling effects on test size of some planktonic foraminifera species from the Oman Margin, Arabian Sea, Paleoceanography, 10, 117 – 122, 1995. Naidu, P. D., and B. A. Malmgren, Relationship between late Quaternary upwelling history and coiling properties of Neogloboquadrina pachyderma and Globigerina bulloides in the Arabian Sea, J. Foraminiferal Res., 26, 64 – 70, 1996. Norris, R. D., Biased extinction and evolutionary trends, Paleobiology, 17, 388 – 399, 1991. Nu¨rnberg, D., J. Bijma, and C. Hemleben, Assessing the reliability of magnesium in foraminiferal calcite as a proxy for water mass temperatures, Geochim. Cosmochim. Acta, 60, 803 – 814, 1996. Ortiz, J. D., A. C. Mix, and R. W. Collier, Environmental control of living symbiotic and asymbiotic foraminifera of the California Current, Paleoceanography, 10, 987 – 1009, 1995. Peterson, R. G., and L. Stramma, Upper-level circulation in the South Atlantic Ocean, Prog. Oceanogr., 26, 1 – 73, 1991. Peterson, R. G., and T. I. Whitworth, The subantarctic and polar fronts in relation to deep water masses through the southwestern Atlantic, J. Geophys. Res., 94, 10,817 – 10,838, 1989. Renaud, S., and D. N. Schmidt, Habitat tracking as a response of the planktic foraminifer Globorotalia truncatulinoides to environmental fluctuations during the last 140 kyr, Mar. Micropaleontol, in press, 2003. Rutherford, S., S. d’ Hondt, and W. Prell, Environmental controls on the geographic distribution of zooplankton diversity, Nature, 400, 749 – 753, 1999. Sanyal, A., N. G. Hemming, G. N. Hanson, and W. Broecker, Evidence for a higher pH in the glacial ocean from boron isotopes in foraminifera, Nature, 373, 234 – 236, 1995. Schiebel, R., B. Hiller, and C. Hemleben, Impacts of storms on Recent planktic foraminiferal test production and CaCO3 flux in the North Atlantic at 47 degrees N, 20 degrees W (JGOFS), Mar. Micropaleontol., 26, 115 – 129, 1995. Schiebel, R., J. Bijma, and C. Hemleben, Population dynamics of the planktic foraminifer Globigerina bulloides from the eastern North Atlantic, Deep Sea Res. Part II, 44, 1701 – 1713, 1997.

Schiebel, R., B. Schmuker, C. Hemleben, and M. Alves, Tracking the Recent and Late Pleistocene Azores Front by the distribution of Planktic Foraminifers, J. Mar. Syst., 37, 213 – 227, 2002. Schneider, R. R., P. J. Mu¨ller, G. Ruhland, G. Meineke, H. Schmidt, and G. Wefer, Late Quaternary surface temperatures and productivity in the east-equatorial South Atlantic: Response to changes in trade/monsoon wind forcing and surface water advection, in The South Atlantic: Present and Past Circulation, edited by G. Wefer et al., pp. 527 – 551, Springer-Verlag, New York, 1996. Shackleton, N. J., and N. D. Opdyke, Oxygen isotope and palaeomagnetic stratigraphy of equatorial Pacific core V28-238: Oxygen isotope temperatures and ice volumes on a 105 and 106 year scale, Quat. Res., 3, 39 – 55, 1973. Spero, H. J., I. Lerche, and D. F. Williams, Opening the carbon isotope ‘‘vital effect’’ box, 2: Quantitative model for interpreting foraminiferal carbon isotope data, Paleoceanography, 6, 639 – 655, 1991. Spindler, M., C. Hemleben, U. Bayer, A. W. H. Be´, and O. R. Anderson, Lunar periodicity of reproduction in the planktonic foraminifer Hastigerina pelagica, Mar. Ecol. Prog. Ser., 1, 61 – 64, 1979. Stein, R., Accumulation of organic carbon in marine sediments, in Lecture Notes in Earth Sciences, edited by S. Bhattacharji et al., 217 pp., Springer-Verlag, New York, 1991. Wefer, G., et al., Late Quaternary surface circulation of the South Atlantic: The stable isotope record and implications for heat transport and productivity, in The South Atlantic: Present and Past Circulation, edited by G. Wefer et al., pp. 461 – 502, Springer-Verlag, New York, 1996. Weyl, P. K., Micropaleontology and ocean surface climate, Science, 202, 475 – 481, 1978. Wolff, T., S. Mulitza, C. Ru¨ hlemann, and G. Wefer, Response of the tropical Atlantic thermocline to late Quaternary trade wind changes, Paleoceanography, 14, 374 – 383, 1999.

J. Bollmann, Geological Institute, ETH Zurich, Sonnegstr. 5, CH-8092 Zurich, Switzerland. ([email protected]) S. Renaud, Pale´oEnvironnements et Pale´oBiosphere, UMR S125 CNRS UCB Lyon 1, Baˆt Ge´ode, 2 rue, Dubois, Campus de la Donna, F69622 Villeurbanne, France. (renaud@cismsun. univ-lyon1.fr) D. N. Schmidt, Department of Geology, Royal Holloway, University of London, Egham, Surrey TW20 OEX, UK. ([email protected])