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.

Calcium isotope ratios of pore fluids and solid phase of organic rich sediments

Pore fluid calcium isotope, calcium concentration and strontium concentration data are used to measure the rates of diagenetic dissolution and precipitation of calcite in deep-sea sediments containing abundant clay and organic material. This type of study of deep-sea sediment diagenesis provides unique information about the ultra-slow chemical reactions that occur in natural marine sediments that affect global geochemical cycles and the preservation of paleo-environmental information in carbonate fossils. For this study, calcium isotope ratios (d44/40Ca) of pore fluid calcium from Ocean Drilling Program (ODP) Sites 984 (North Atlantic) and 1082 (off the coast of West Africa) were measured to augment available pore fluid measurements of calcium and strontium concentration. Both study sites have high sedimentation rates and support quantitative sulfate reduction, methanogenesis and anaerobic methane oxidation. The pattern of change of d44/40Ca of pore fluid calcium versus depth at Sites 984 and 1082 differs markedly from that of previously studied deep-sea Sites like 590B and 807, which are composed of nearly pure carbonate sediment. In the 984 and 1082 pore fluids, d44/40Ca remains elevated near seawater values deep in the sediments, rather than shifting rapidly toward the d44/40Ca of carbonate solids. This observation indicates that the rate of calcite dissolution is far lower than at previously studied carbonate-rich sites. The data are fit using a numerical model, as well as more approximate analytical models, to estimate the rates of carbonate dissolution and precipitation and the relationship of these rates to the abundance of clay and organic material. Our models give mutually consistent results and indicate that calcite dissolution rates at Sites 984 and 1082 are roughly two orders of magnitude lower than at previously studied carbonate-rich sites, and the rate correlates with the abundance of clay. Our calculated rates are conservative for these sites (the actual rates could be significantly slower) because other processes that impact the calcium isotope composition of sedimentary pore fluid have not been included. The results provide direct geochemical evidence for the anecdotal observation that the best-preserved carbonate fossils are often found in clay or organic-rich sedimentary horizons. The results also suggest that the presence of clay minerals has a strong passivating effect on the surfaces of biogenic carbonate minerals, slowing dissolution dramatically even in relation to the already-slow rates typical of carbonate-rich sediments.

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.

Stable isotope record and sea surface reconstruction of sediment core MD88-770

Hydrographical changes of the southern Indian Ocean over the last 230 kyr, is reconstructed using a 17-m-long sediment core (MD 88 770; 46°01'S 96°28'E, 3290m). The oxygen and carbon isotopic composition of planktonic (N. pachyderma sinistra and G. bulloides) and benthic (Cibicidoides wuellerstorfi, Epistominella exigua, and Melonis barleeanum) foraminifera have been analysed. Changes in sea surface temperatures (SST) are calculated using diatom and foraminiferal transfer functions. A new core top calibration for the Southern Ocean allows an extension of the method developed in the North Atlantic to estimate paleosalinities (Duplessy et al., 1991). The age scale is built using accelerator mass spectrometry (AMS) 14C dating of N. pachyderma s. for the last 35 kyr, and an astronomical age scale beyond. Changes in surface temperature and salinity clearly lead (by 3 to 7 kyr) deep water variations. Thus changes in deep water circulation are not the cause of the early response of the surface Southern Ocean to climatic changes. We suggest that the early warming and cooling of the Southern Ocean result from at least two processes acting in different orbital bands and latitudes: (1) seasonality modulated by obliquity affects the high-latitude ocean surface albedo (sea ice coverage) and heat transfer to and from the atmosphere; (2) low-latitude insolation modulated by precession influences directly the atmosphere dynamic and related precipitation/ evaporation changes, which may significantly change heat transfer to the high southern latitudes, through their control on latitudinal distribution of the major frontal zones and on the conditions of intermediate and deep water formation.