A multiproxy analysis of late MIS 6 to early MIS 5d OSP Site 161-975 sediments

A multiproxy analysis based on planktic foraminiferal abundances, derived SSTs, and stable planktic isotopes measurements together with alkenone abundances and Uk'37 SSTs was performed on late MIS 6 to early MIS 5d sediment recovered from Site 975 (ODP Leg 161) in the South Balearic Islands basin (Western Mediterranean) with emphasis on reconstructing the climate progression of the last interglacial period. A number of abrupt climate changes related to alternative influence of nutrient rich northern and oligotrophic southern water masses were revealed. Heinrich event 11 and cooling events C27, C26, C25, C24, C23, which have been previously described in the North Atlantic, were recognized. However, in comparison to the eastern North Atlantic mid-latitude region, events C27 and C26 at Site 975 seem to be significantly more pronounced. Together with evidence of a two-phase climate optimum with maximum SSTs reached during its later phase, this implies a close similarity in climate dynamics between the Western Mediterranean and the Nordic seas. We propose that postglacial effects in the Nordic seas had an influence on the western Mediterranean climate via atmospheric circulation and that these effects competed with the insolation force.

Stable isotope ratios of foraminifera from last glacial-interglacial sediments of the southern high latitudes

A set of numerical equations is developed to estimate past sea surface temperatures (SST) from fossil Antarctic diatoms. These equations take into account both the biogeographic distribution and experimentally derived silica dissolution. The data represent a revision and expansion of a floral data base used previously and includes samples resulting from progressive opal dissolution experiments. Factor analysis of 166 samples (124 Holocene core top and 42 artificial samples) resolved four factors. Three of these factors depend on the water mass distribution (one Subantarctic and two Antarctic assemblages); factor 4 corresponds to a 'dissolution assemblage'. Inclusion of this factor in the data analysis minimizes the effect of opal dissolution on the assemblages and gives accurate estimates of SST over a wide range of biosiliceous dissolution. A transfer function (DTF 166/34/4) is derived from the distribution of these factors versus summer SST. Its standard error is +/- 1°C in the -1 to +10 °C summer temperature range. This transfer function is used to estimate SST changes in two southern ocean cores (43°S and 55°S) which cover the last climatic cycle. The time scale is derived from the changes in foraminiferal oxygen and carbon isotopic ratios. The reconstructed SST records present strong analogies with the air temperature record over Antarctica at the Vostok site, derived from changes in the isotopic ratio of the ice. This similarity may be used to compare the oceanic isotope stratigraphy and the Vostok time scale derived from ice flow model. The oceanic time scale, if taken at face value, would indicate that large changes in ice accumulation rates occurred between warm and cold periods.

Calculated sea surface temperatures of two sediment cores

The reconstruction of ocean history employs a large variety of methods with origins in the biological, chemical, and physical sciences, and uses modern statistical techniques for the interpretation of extensive and complex data sets. Various sediment properties deliver useful information for reconstructing environmental parameters. Those properties that have a close relationship to environmental parameters are called ''proxy variables'' (''proxies'' for short). Proxies are measurable descriptors for desired (but unobservable) variables. Surface water temperature is probably the most important parameter for describing the conditions of past oceans and is crucial for climate modelling. Proxies for temperature are: abundance of microfossils dwelling in surface waters, oxygen isotope composition of planktic foraminifers, the ratio of magnesium or strontium to calcium in calcareous shells or the ratio of certain organic molecules (e.g. alkenones produced by coccolithophorids). Surface water salinity, which is important in modelling of ocean circulation, is much more difficult to reconstruct. At present there is no established method for a direct determination of this parameter. Measurements associated with the paleochemistry of bottom waters to reconstruct bottom water age and flow are made on benthic foraminifers, ostracodes, and deep-sea corals. Important geochemical tracers are d13C and Cd/Ca ratios. When using benthic foraminifers, knowledge of the sediment depth habitat of species is crucial. Reconstructions of productivity patterns are of great interest because of important links to current patterns, mixing of water masses, wind, the global carbon cycle, and biogeography. Productivity is reflected in the flux of carbon into the sediment. There are a number of fluxes other than those of organic carbon that can be useful in assessing productivity fluctuations. Among others, carbonate and opal flux have been used, as well as particulate barite. Furthermore, microfossil assemblages contain clues to the intensity of production as some species occur preferentially in high-productivity regions while others avoid these. One marker for the fertility of sub-surface waters (that is, nutrient availability) is the carbon isotope ratio within that water (13C/12C, expressed as d13C). Carbon isotope ratios in today's ocean are negatively correlated with nitrate and phosphate contents. Another tracer of phosphate content in ocean waters is the Cd/Ca ratio. The correlation between this ratio and phosphate concentrations is quite well documented. A rather new development to obtain clues on ocean fertility (nitrate utilization) is the analysis of the 15N/14N ratio in organic matter. The fractionation dynamics are analogous to those of carbon isotopes. These various ratios are captured within the organisms growing within the tagged water. A number of reconstructions of the partial pressure of CO2 have been attempted using d13C differences between planktic and benthic foraminifers and d13C values of bulk organic material or individual organic components. To define the carbon system in sea water, two elements of the system have to be known in addition to temperature. These can be any combination of total CO2 , alkalinity, or pH. To reconstruct pH, the boron isotope composition of carbonates has been used. Ba patterns have been used to infer the distribution of alkalinity in past oceans. Information relating to atmospheric circulationand climate is transported to the ocean by wind or rivers, in the form of minerals or as plant andanimal remains. The most useful tracers in this respect are silt-sized particles and pollen.

Shell size variation of Neogloboquadrina pachyderma sin. in the Norwegian-Greenland Sea during the last 1.3 Myr

We present measurements of the maximum diameter of the planktonic foraminifer Neogloboquadrina pachyderma sin. from six sediment cores (Ocean Drilling Program sites 643, 644, 907, 909, 985 and 987) from the Norwegian-Greenland Sea. Our data show a distinct net increase in mean shell size of N. pachyderma sin. at all sites during the last 1.3 Ma, with largest shell sizes reached after 0.4 Ma. External factors such as glacial-interglacial variability and carbonate dissolution alone cannot account for the observed variation in mean shell size of N. pachyderma sin. We consider the observed shell size increase to mirror an evolutionary trend towards better adaptation of N. pachyderma sin. to the cold water environment after 1.1-1.0 Ma. Probably, the Mid Pleistocene climate shift and the associated change of amplitude and frequency of glacial-interglacial fluctuations have triggered the evolution of this planktonic foraminifer. Oxygen and carbon stable isotope analyses of different shell size classes indicate that the observed shell size increase could not be explained by the functional concept that larger shells promote increasing sinking velocities during gametogenesis. For paleoceanographic reconstructions, the evolutionary adaptation of Neogloboquadrina pachyderma sin. to the cold water habitat has significant implications. Carbonate sedimentation in highest latitudes is highly dependent on the presence of this species. In the Norwegian-Greenland Sea, carbonate-poor intervals before 1.1 Ma are, therefore, not necessarily related to severe glacial conditions. They are probably attributed to the absence of this not yet polar-adapted species. Further, transfer function and modern analog techniques used for the reconstruction of surface water conditions in high latitudes could, therefore, contain a large range of errors if they were applied to samples older than 1.1-1.0 Myrs.

Planktonic foraminifera in the Arabian Sea

Variations in primary productivity (PP) have been reconstructed in eutrophic, mesotrophic and oligotrophic parts of the Arabian Sea over the past 135 000 years applying principal component analysis and transfer function to planktic foraminiferal assemblages. Temporal variation in paleoproductivity is most pronounced in the mesotrophic northern (NAST site) and oligotrophic eastern (EAST site) Arabian Sea, and comparatively weak in the western eutrophic GeoB 3011-1 site in the upwelling area off Oman. Higher PP during interglacials (250-320 g C/m**2 year) than during cold stages (210-270 g C/m**2 year) at GeoB 3011-1 could have been caused by a strengthened upwelling during intensified summer monsoons and increased wind velocities. At NAST, during interglacials, PP is estimated to exceed g C/m**2 year 1, and during glacials to be as low as 140-180 g C/m**2 year. These fluctuations may result from a (1) varying impact of filaments that are associated to the Oman coastal upwelling, and (2) from open-ocean upwelling associated to the Findlater Jet. At EAST, highest productivity of about 380 g C/m**2 year is documented for the transition from isotope stage 5 to 4. We suggest that during isotope stages 2, 4, 5.2, the transition 5/4, and the end of stage 6, deep mixing of surface waters was caused by moderate to strong winter monsoons, and induced an injection of nutrients into the euphotic layer leading to enhanced primary production. The deepening of the mixed layer during these intervals is confirmed by an increased concentration of deep-dwelling planktic foraminiferal species. A high-productivity event in stage 3, displayed by estimated PP values, and by planktic foraminifera and radiolaria flux and accumulation rate, likely resulted from a combination of intensified SW monsoons with moderate to strong NE monsoons. Differential response of Globigerina bulloides, Globigerinita glutinata and mixed layer species to the availability of food is suited to subdivide productivity regimes on a temporal and spatial scale.

Calculated sea level records and deep ocean temperatures from DSDP and ODP sites

During the mid-Pleistocene transition the dominant 41 ka periodicity of glacial cycles transitioned to a quasi-100 ka periodicity for reasons not yet known. This study investigates the potential role of deep ocean hydrography by examining oxygen isotope ratios in benthic foraminifera. Oxygen isotope records from the Atlantic, Pacific and Indian Ocean basins are separated into their ice volume and local temperature/hydrography components using a piece-wise linear transfer function and a temperature calibration. Although our method has certain limitations, the deep ocean hydrography reconstructions show that glacial deep ocean temperatures approached freezing point as the mid-Pleistocene transition progressed. Further analysis suggests that water mass reorganisation could have been responsible for these temperature changes, leading to such stable conditions in the deep ocean that some obliquity cycles were skipped until precessional forcing triggered deglaciation, creating the apparent quasi-100 ka pattern. This study supports previous work that suggests multiples of obliquity cycles dominate the quasi-100 ka glacial cycles with precession components driving deglaciations.

Planktic foraminifer-based sea-surface temperature reconstruction for sediment cores JM97-948/2A and MD95-2011 from the Vøring Plateau, eastern Norwegian Sea

The early to mid-Holocene thermal optimum is a well-known feature in a wide variety of paleoclimate archives from the Northern Hemisphere. Reconstructed summer temperature anomalies from across northern Europe show a clear maximum around 6000 years before present (6 ka). For the marine realm, Holocene trends in sea-surface temperature reconstructions for the North Atlantic and Norwegian Sea do not exhibit a consistent pattern of early to mid- Holocene warmth. Sea-surface temperature records based on alkenones and diatoms generally show the existence of a warm early to mid-Holocene optimum. In contrast, several foraminifer and radiolarian based temperature records from the North Atlantic and Norwegian Sea show a cool mid- Holocene anomaly and a trend towards warmer temperatures in the late Holocene. In this paper, we revisit the foraminifer record from the Vøring Plateau in the Norwegian Sea. We also compare this record with published foraminifer based temperature reconstructions from the North Atlantic and with modelled (CCSM3) upper ocean temperatures. Model results indicate that while the seasonal summer warming of the seasurface was stronger during the mid-Holocene, sub-surface depths experienced a cooling. This hydrographic setting can explain the discrepancies between the Holocene trends exhibited by phytoplankton and zooplankton based temperature proxy records.

Global compilation of last glacial maximum oxygen isotopic ratios for planktonic and benthic foraminifera

We use the fully coupled atmosphere-ocean three-dimensional model of intermediate complexity iLOVECLIM to simulate the climate and oxygen stable isotopic signal during the Last Glacial Maximum (LGM, 21 000 yr). By using a model that is able to explicitly simulate the sensor (d18O), results can be directly compared with data from climatic archives in the different realms. Our results indicate that iLOVECLIM reproduces well the main feature of the LGM climate in the atmospheric and oceanic components. The annual mean d18O in precipitation shows more depleted values in the northern and southern high latitudes during the LGM. The model reproduces very well the spatial gradient observed in ice core records over the Greenland ice-sheet. We observe a general pattern toward more enriched values for continental calcite d18O in the model at the LGM, in agreement with speleothem data. This can be explained by both a general atmospheric cooling in the tropical and subtropical regions and a reduction in precipitation as confirmed by reconstruction derived from pollens and plant macrofossils. Data-model comparison for sea surface temperature indicates that iLOVECLIM is capable to satisfyingly simulate the change in oceanic surface conditions between the LGM and present. Our data-model comparison for calcite d18O allows investigating the large discrepancies with respect to glacial temperatures recorded by different microfossil proxies in the North Atlantic region. The results argue for a trong mean annual cooling between the LGM and present (>6°C), supporting the foraminifera transfer function reconstruction but in disagreement with alkenones and dinocyst reconstructions. The data-model comparison also reveals that large positive calcite d18O anomaly in the Southern Ocean may be explained by an important cooling, although the driver of this pattern is unclear. We deduce a large positive d18Osw anomaly for the north Indian Ocean that contrasts with a large negative d18Osw anomaly in the China Sea between the LGM and present. This pattern may be linked to changes in the hydrological cycle over these regions. Our simulation of the deep ocean suggests that changes in d18Osw between the LGM and present are not spatially homogenous. This is supported by reconstructions derived from pore fluids in deep-sea sediments. The model underestimates the deep ocean cooling thus biasing the comparison with benthic calcite d18O data. Nonetheless, our data-model comparison support a heterogeneous cooling of few degrees (2-4°C) in the LGM Ocean.

Planktic foraminiferal abundance and faunal-derived sea surface temperature of eastern equatorial Pacific ODP Site 202-1240

Glacial cooling (~1-5°C) in the eastern equatorial Pacific (EEP) cold tongue is often attributed to increased equatorial upwelling, stronger advection from the Peru-Chile Current (PCC), and to the more remote subpolar southeastern Pacific water mass. However, evidence is scarce for identifying unambiguously which process plays a more important role in driving the large glacial cooling in the EEP. To address this question, here we adopt a faunal calibration approach using planktic foraminifers with a new compilation of coretop data from the eastern Pacific, and present new downcore variation data of fauna assemblage and estimated sea surface temperatures (SSTs) for the past 160 ka (Marine Isotope Stage (MIS) 6) from ODP Site 1240 in the EEP. With significant improvement achieved by adding more coretop data from the eastern boundary current, our downcore calibration results indicate that most of the glacial cooling episodes over the past 160 ka in the EEP are attributable to increased influence from the subpolar water mass from high latitudes of the southern Pacific. By applying this new calibration of the fauna SST transfer function to a latitudinal transect of eastern Pacific (EP) cores, we find that the subpolar water mass has been a major dynamic contributor to EEP cold tongue cooling since MIS 6.

Planktonic foraminifera of sediment core MD01-2390

Sea-surface temperature (SST) estimates in the sediment core MD01-2390 based on planktonic foraminiferal species abundances using five different transfer function techniques suggest nearly unchanged or unusually higher temperatures in the tropical southern South China Sea (SCS) during the Last Glacial Maximum (LGM) relative to modern temperatures. These results are in contrast to substantial cooling of 2-5 °C inferred by geochemical (Uk'37, Mg/Ca ratios) and terrestrial proxies from the western tropical Pacific region. Using multivariate statistics we show that the glacial southern SCS harboured unique planktonic foraminiferal assemblages that have no modern analogs. Analyses of faunal variation through the core reveal that planktonic foraminiferal assemblages responded to temperature changes inferred from Mg/Ca data but that this signal is subdued by superimposed variations in the relative abundance of Pulleniatina obliquiloculata and Neogloboquadrina pachyderma (dextral). These species occur in glacial samples at proportions that are not observed in the calibration data set. The glacial high abundance of N. pachyderma (dextral) are interpreted to reflect a seasonal (winter) inflow of cold surface water from the northeast via the Bashi Strait due to the combined effects of an intensified winter monsoon, a southward shift of the polar front and the eastward migration of the Kuroshio Current. In contrast, processes controlling the high relative abundances of P. obliquiloculata during the LGM may be unique to the southern SCS. We propose a scenario involving a stronger (winter) mixing or enhanced upwelling due to an intensified winter monsoon that prevented shallow-dwelling, warm indicators to establish larger populations during the LGM. Our results indicate that a no-analog behaviour of planktonic foraminifera faunas is responsible for the warm glacial conditions in this part of the western Pacific warm pool as implied by foraminiferal transfer functions and that a more significant surface cooling in the region as implied by terrestrial and geochemical (Mg/Ca ratios; alkenone unsaturation index) marine proxies is a more likely scenario.

PaleoTools - extended software package for quantitative paleoenvironmental reconstructions

The Imbrie and Kipp transfer function method (IKM) and the modern analog technique (MAT) are accepted tools for quantitative paleoenvironmental reconstructions. However, no uncomplicated, flexible software has been available to apply these methods on modern computer devices. For this reason the software packages PaleoToolBox, MacTransfer, WinTransfer, MacMAT, and PanPlot have been developed. The PaleoToolBox package provides a flexible tool for the preprocessing of microfossil reference and downcore data as well as hydrographic reference parameters. It includes procedures to randomize the raw databases; to switch specific species in or out of the total species list; to establish individual ranking systems and their application on the reference and downcore databasessemi; and to convert the prepared databases into the file formats of IKM and MAT software for estimation of paleohydrographic parameters.

Revised calibration of the MBT-CBT paleotemperature proxy

The MBT-CBT proxy for the reconstruction of paleotemperatures and past soil pH is based on the distribution of branched glycerol dialkyl glycerol tetraether (brGDGT) membrane lipids. The Methylation of Branched Tetraether (MBT) and the Cyclisation of Branched Tetraether (CBT) indices were developed to quantify these distributions, and significant empirical relations between these indices and annual mean air temperature (MAT) and/or soil pH were found in a large data set of soils. In this study, we extended this soil dataset to 278 globally distributed surface soils. Of these soils, 26% contains all nine brGDGTs, while in 63% of the soils the seven most common brGDGTs were detected, and the latter were selected for calibration purposes. This resulted in new transfer functions for the reconstruction of pH based on the CBT index: pH = 7.90-1.97 × CBT (r**2 = 0.70; RMSE = 0.8; n = 176), as well as for MAT based on the CBT index and methylation index based on the seven most abundant GDGTs (defined as MBT'): MAT = 0.81-5.67 × CBT + 31.0 × MBT' (r**2 = 0.59; RMSE = 5.0 °C; n = 176). The new transfer function for MAT has a substantially lower correlation coefficient than the original equation (r**2 = 0.77). To investigate possible improvement of the correlation, we used our extended global surface soil dataset to statistically derive the indices that best describe the relations of brGDGT composition with MAT and soil pH. These new indices, however, resulted in only a relatively minor increase in correlation coefficients, while they cannot be explained straightforwardly by physiological mechanisms. The large scatter in the calibration cannot be fully explained by local factors or by seasonality, but MAT for soils from arid regions are generally substantially (up to 20 °C) underestimated, suggesting that absolute brGDGT-based temperature records for these areas should be interpreted with caution. The applicability of the new MBT'-CBT calibration function was tested using previously published MBT-CBT-derived paleotemperature records covering the last deglaciation in Central Africa and East Asia, the Eocene-Oligocene boundary and the Paleocene-Eocene thermal maximum. The results show that trends remain similar in all records, but that absolute temperature estimates and the amplitude of temperature changes are lower for most records, and generally in better agreement with independent proxy data.