Sea-surface temperature reconstruction of sediments from the North Pacific and North Atlantic

Holocene climate variability is investigated in the North Pacific and North Atlantic realms, using alkenone-derived sea-surface temperature (SST) records as well as a millennial scale simulation with a coupled atmosphere-ocean general circulation model (AOGCM). The alkenone SST data indicate a temperature increase over almost the entire North Pacific from 7 cal kyr BP to the present. A dipole pattern with a continuous cooling in the northeastern Atlantic and a warming in the eastern Mediterranean Sea and the northern Red Sea is detected in the North Atlantic realm. Similarly, SST variations are opposite in sign between the northeastern Pacific and the northeastern Atlantic. A 2300 year long AOGCM climate simulation reveals a similar SST seesaw between the northeastern Pacific and the northeastern Atlantic on centennial time scales. Our analysis of the alkenone SST data and the model results suggests fundamental inter-oceanic teleconnections during the Holocene.

Chemical composition of sediment core GeoB12309-5 from the northern Arabian Sea, its pore water profile and age depth relationship as well as CTD data obtained during cruise M74/3

The Arabian Sea off the Pakistan continental margin is characterized by one of the world's largest oxygen minimum zones (OMZ). The lithology and geochemistry of a 5.3 m long gravity core retrieved from the lower boundary of the modern OMZ (956 m water depth) were used to identify late Holocene changes in oceanographic conditions and the vertical extent of the OMZ. While the lower part of the core (535 - 465 cm, 5.04 - 4.45 cal kyr BP, Unit 3) is strongly bioturbated indicating oxic bottom water conditions, the upper part of the core (284 - 0 cm, 2.87 cal kyr BP to present, Unit 1) shows distinct and well-preserved lamination, suggesting anoxic bottom waters. The transitional interval from 465 to 284 cm (4.45 - 2.87 cal kyr BP, Unit 2) contains relicts of lamination which are in part intensely bioturbated. These fluctuations in bioturbation intensity suggest repetitive changes between anoxic and oxic/suboxic bottom-water conditions between 4.45 - 2.87 cal kyr BP. Barium excess (Baex) and total organic carbon (TOC) contents do not explain whether the increased TOC contents found in Unit 1 are the result of better preservation due to low BWO concentrations or if the decreased BWO concentration is a result of increased productivity. Changes in salinity and temperature of the outflowing water from the Red Sea during the Holocene influenced the water column stratification and probably affected the depth of the lower boundary of the OMZ in the northern Arabian Sea. Even if we cannot prove certain scenarios, we propose that the observed downward shift of the lower boundary of the OMZ was also impacted by a weakened Somali Current and a reduced transport of oxygen-rich Indian Central Water into the Arabian Sea, both as a response to decreased summer insolation and the continuous southward shift of the Intertropical Convergence Zone during the late Holocene.

High resolution isotope and Mg/Ca record of sediment core GeoB12615-4

The importance of intermediate water masses in climate change and ocean circulation has been emphasized recently. In particular, Southern Ocean Intermediate Waters (SOIW), such as Antarctic Intermediate Water and Subantarctic Mode Water, are thought to have acted as active interhemispheric transmitter of climate anomalies. Here we reconstruct changes in SOIW signature and spatial and temporal evolution based on a 40 kyr time series of oxygen and carbon isotopes as well as planktic Mg/Ca based thermometry from Site GeoB12615-4 in the western Indian Ocean. Our data suggest that SOIW transmitted Antarctic temperature trends to the equatorial Indian Ocean via the "oceanic tunnel" mechanism. Moreover, our results reveal that deglacial SOIW carried a signature of aged Southern Ocean deep water. We find no evidence of increased formation of intermediate waters during the deglaciation.

Dinoflagellate cyst of sediment core M35003-4 southeast of Grenada

High-resolution, well-dated calcareous dinoflagellate cyst and organic carbon records from a 58 kyr sediment core (M35003-4) located southeast of the island of Grenada show that rapid and pronounced changes in cyst association and accumulation and organic carbon deposition occurred, controlled by (1) a significant southward shift in the position of the North Equatorial Current during the last glacial period and the Younger Dryas cold interval and (2) rapid changes in local productivity in marine isotopic stage 3 that are associated with variations in Orinoco River nutrient discharge and coastal upwelling strength. Prominent cyst accumulation peaks representing extremely oligotrophic and stratified thermocline conditions mimic the Greenland ice core and northern Atlantic Dansgaard/Oeschger stadials and Heinrich events. We provide new evidence for a coupled tropical/high-latitude Atlantic climate system during the last glacial period and suggest that changes in the zonality of the low-latitude winds may play an important role in modulating rapid interhemispheric climate variability.

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.

Seawater carbonate chemistry and biological processes of foraminifera, Globigerinoides ruber and Globigerinella siphonifera during experiments, 2011

Specimens of two species of planktic foraminifera, Globigerinoides ruber and Globigerinella siphonifera, were grown under controlled laboratory conditions at a range of temperatures (18-31 °C), salinities (32-44 psu) and pH levels (7.9-8.4). The shells were examined for their calcium isotope compositions (d44/40Ca) and strontium to calcium ratios (Sr/Ca) using Thermal Ionization Mass Spectrometry and Inductively Coupled Plasma Mass Spectrometry. Although the total variation in d44/40Ca (~0.3 per mill) in the studied species is on the same order as the external reproducibility, the data set reveals some apparent trends that are controlled by more than one environmental parameter. There is a well-defined inverse linear relationship between d44/40Ca and Sr/Ca in all experiments, suggesting similar controls on these proxies in foraminiferal calcite independent of species. Analogous to recent results from inorganically precipitated calcite, we suggest that Ca isotope fractionation and Sr partitioning in planktic foraminifera are mainly controlled by precipitation kinetics. This postulation provides us with a unique tool to calculate precipitation rates and draws support from the observation that Sr/Ca ratios are positively correlated with average growth rates. At 25 °C water temperature, precipitation rates in G. siphonifera and G. ruber are calculated to be on the order of 2000 and 3000 µmol/m**2/h, respectively. The lower d44/40Ca observed at 29 °C in both species is consistent with increased precipitation rates at high water temperatures. Salinity response of d44/40Ca (and Sr/Ca) in G. siphonifera implies that this species has the highest precipitation rates at the salinity of its natural habitat, whereas increasing salinities appear to trigger higher precipitation rates in G. ruber. Isotope effects that cannot be explained by precipitation rate in planktic foraminifera can be explained by a biological control, related to a vacuolar pathway for supply of ions during biomineralization and a pH regulation mechanism in these vacuoles. In case of an additional pathway via cross-membrane transport, supplying light Ca for calcification, the d44/40Ca of the reservoir is constrained as -0.2 per mill relative to seawater. Using a Rayleigh distillation model, we calculate that calcification occurs in a semi-open system, where less than half of the Ca supplied by vacuolization is utilized for calcite precipitation. Our findings are relevant for interpreting paleo-proxy data on d44/40Ca and Sr/Ca in foraminifera as well as understanding their biomineralization processes.

Age determination of sediment cores from the Panama Basin

High-resolution percent Corg and delta18Oforam records obtained from Panama Basin core Atlantis II 54-25PC and additional data from nearby core P7 show that enhanced burial of organic carbon has characterized every major glacial period for the last 500 kyr in that area. Both Corg concentration and mass accumulation rate profiles exhibit a sawtooth pattern with maxima occurring typically in the later stages of glacial periods. Comparison with dust records suggests that the carbon accumulation rate profile reflects both the upwelling history and a variable rate of iron input during the late Quaternary. The sawtooth character may derive from increased wind velocities and rates of upwelling during glacials which are indirectly related to ice volume (Sarnthein et al., 1988). The rapid decline in export production at the end of glacials in the equatorial Pacific may be attributed to the retreat of ice sheets (thus reduced wind velocities and upwelling) coupled with a coincident decline in atmospheric dust load and/or delivery rate. The Corg accumulation rate profiles do not correlate well with atmospheric CO2 records. For example, atmospheric CO2 was already at a minimum 40 kyr ago when production in the Panama Basin began increasing dramatically, commensurate with an increase in global dust levels. Using the relationship between the degree of photosynthetic fractionation and the concentration of free CO2 in the surface ocean postulated by Popp et al. (1989), delta13Corg measurements made on core P7 show that Panama Basin surface waters have been supplying CO2 to the atmosphere continually for at least the last 50 kyr. There is no evidence for a flux of CO2 into the surface ocean in this area at any time during this period despite the higher production. If the Panama Basin cores are representative of the eastern and central equatorial Pacific, then these observations weaken the influence on CO2 drawdown postulated for increased glacial productivity at low latitudes.

Stable isotope record, Mg/Ca ratios and sea surface reconstruction for sediment core MD06-3067

We present sea surface, upper thermocline, and benthic d18O data, as well as temperature and paleoproductivity proxy data, from the International Marine Global Change Study Program (IMAGES) Core MD06-3067 (6°31'N, 126°30'E, 1575 m water depth), located in the western equatorial Pacific Ocean within the flow path of the Mindanao Current. Our records reveal considerable glacial-interglacial and suborbital variability in the Mindanao Dome upwelling over the last 160 kyr. Dome activity generally intensified during glacial intervals resulting in cooler thermocline waters, whereas it substantially declined during interglacials, in particular in the early Holocene and early marine oxygen isotope stage (MIS) 5e, when upwelling waters did not reach the thermocline. During MIS 3 and MIS 2, enhanced surface productivity together with remarkably low SST and low upper ocean thermal contrast provide evidence for episodic glacial upwelling to the surface, whereas transient surface warming marks periodic collapses of the Mindanao Dome upwelling during Heinrich events. We attribute the high variability during MIS 3 and MIS 2 to changes in the El Niño Southern Oscillation state that affected boreal winter monsoonal winds and upper ocean circulation. Glacial upwelling intensified when a strong cyclonic gyre became established, whereas El Niño-like conditions during Heinrich events tended to suppress the cyclonic circulation, reducing Ekman transport. Thus, our findings demonstrate that variations in the Mindanao Dome upwelling are closely linked to the position and intensity of the tropical convection and also reflect far-field influences from the high latitudes.

Results of Heinrich event 1 simulation from the University of Victoria Earth System-Climate Model (UVic ESCM) and the Community Climate System Model (CCSM3)

We investigated changes in tropical climate and vegetation cover associated with abrupt climate change during Heinrich Event 1 (HE1, ca. 17.5 ka BP) using two different global climate models: the University of Victoria Earth System-Climate Model (UVic ESCM) and the Community Climate System Model version 3 (CCSM3). Tropical South American and African pollen records suggest that the cooling of the North Atlantic Ocean during HE1 influenced the tropics through a southward shift of the rain belt. In this study, we simulated the HE1 by applying a freshwater perturbation to the North Atlantic Ocean. The resulting slowdown of the Atlantic Meridional Overturning Circulation was followed by a temperature seesaw between the Northern and Southern Hemispheres, as well as a southward shift of the tropical rain belt. The shift and the response pattern of the tropical vegetation around the Atlantic Ocean were more pronounced in the CCSM3 than in the UVic ESCM simulation. For tropical South America, opposite changes in tree and grass cover were modeled around 10° S in the CCSM3 but not in the UVic ESCM. In tropical Africa, the grass cover increased and the tree cover decreased around 15° N in the UVic ESCM and around 10° N in the CCSM3. In the CCSM3 model, the tree and grass cover in tropical Southeast Asia responded to the abrupt climate change during the HE1, which could not be found in the UVic ESCM. The biome distributions derived from both models corroborate findings from pollen records in southwestern and equatorial western Africa as well as northeastern Brazil.