Plant-wax, alkenone and major element analyses of ODP Site 659 for the Last Glacial cycle and associated coretop data

We present a hydrologic reconstruction of the Sahara-Sahel transition, covering the complete last glacial cycle (130 ka), based on a combination of plant-wax-specific hydrogen (dD) and carbon isotopes (d13C). The dD and d13C signatures of long-chain n-alkanes from ODP Site 659 off NW Africa reveal a significant anti-correlation. Complementary to published pollen data, we infer that this plant-wax signal reflects sensitive responses of the vegetation cover to precipitation changes in the Sahel region, as well as varying contributions from biomes north of the Sahara (C3 domain) by North-East Trade Winds (NETW). During arid phases, especially the northern parts of the Sahel likely experienced crucial water stress, which resulted in a pronounced contraction of the vegetation cover, thus reducing the amount of C4 plant waxes from the region. The increase in NETW strength during dry periods further promoted a more pronounced C3-plant-wax signal derived from the North African C3 plant domain. During humid periods, the C4-dominated Sahelian environments spread northward into the Saharan realm, in association with lower NETW inputs of C3 plant waxes. Arid-humid cycles deduced from plant-wax dD are in accordance with concomitant changes in weathering intensity reflected in varying major element distributions. Environmental shifts are generally linked to periods with large fluctuations in Northern Hemisphere summer insolation. During Marine Isotope Stages 2 and 3, when insolation variability was low, coupling of the hydrologic regime to alkenone-based estimates of NE Atlantic sea-surface temperatures becomes apparent.

Modeling the sensitivity of precipitation oxygen isotopes to the Last Glacial Maximum boundary conditions: A study using Community Atmosphere Model (CAM3.0)

To understand the validity of d18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (d18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and d18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and d18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in d18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the d18Oprecip distribution among the simulations.

Mean sortable silt values from Holocene and glacial Scotia Sea sediments

The Antarctic Circumpolar Current is key to the mixing and ventilation of the world's oceans. This current flows from west to east between about 45° and 70° S connecting the Atlantic, Pacific and Indian oceans, and is driven by westerly winds and buoyancy forcing. High levels of productivity in the current regulate atmospheric CO2 concentrations. Reconstructions of the current during the last glacial period suggest that flow speeds were faster or similar to present, and it is uncertain whether the strength and position of the westerly winds changed. Here we reconstruct Antarctic Circumpolar Current bottom speeds through the constricting Drake Passage and Scotia Sea during the Last Glacial Maximum and Holocene based on the mean grain size of sortable silt from a suite of sediment cores. We find essentially no change in bottom flow speeds through the region, and, given that the momentum imparted by winds, and modulated by sea-ice cover, is balanced by the interaction of these flows with the seabed, this argues against substantial changes in wind stress. However, glacial flow speeds in the sea-ice zone south of 56° S were significantly slower than present, whereas flow in the north was faster, but not significantly so. We suggest that slower flow over the rough topography south of 56° S may have reduced diapycnal mixing in this region during the last glacial period, possibly reducing the diapycnal contribution to the Southern Ocean overturning circulation.

Geochemistry of last glacial sediments of Isla de los Estados, southeastern Tierra del Fuego

The island of Isla de los Estados is situated at 54.5°S, 64°W, east of Argentinian Tierra del Fuego, and is located in a sensitive geographic position in relation to the zonal circulation between Antarctica and South America. Its terrestrial records of the last deglaciation, recording atmospheric conditions but within an oceanic setting, can help to clarify changes of regional circulation patterns, both atmospheric and marine. Here, we present geochemical analyses from 16-10 ka cal BP of a peat core from Lago Galvarne Bog at the northern coast of the island, and a lake sediment core from Laguna Cascada 3 km further south. The data comprise TC, TN, loss on ignition analyses and continuous XRF scanning on both cores as well as age-depth modeling based on AMS-14C dating. Deglaciation and onset of peat formation in the coastal areas began before 16 ka cal BP followed by a rapid glacial retreat and the start of lacustrine sedimentation further inland. Data suggest initially windy conditions with permafrost succeeded by gradually warmer and wetter conditions until ca 14.5 ka cal BP. The warming trend slows down until ca 13.5 ka cal BP, followed by arid conditions culminating around 12.8 ka cal BP. Our data suggest fairly warm conditions and the establishment of denser peat and forest vegetation ca 10.6 ka cal BP, contemporaneous with the onset of the Antarctic thermal optimum. This indicates large-scale shifts in the placement of zonal flow and the Westerlies at the beginning of the Holocene.

Last Glacial Maximum carbon isotopic composition

We present Holocene and last glacial maximum (LGM) oxygen and carbon isotope measurements on Planulina wuellerstorfi in six southeast Pacific cores. Sedimentation rates are low in this part of the ocean, and measurements were made on individual foraminiferal shells in order to identify the Holocene and glacial individuals on the basis of their extreme d18O. The new d13C data were combined with previous P. wuellerstorfi data for interpretation of global thermohaline circulation. Data from the Southern Ocean were examined closely for regional coherency and a few anomalous d13C values suspected of having productivity overprint were removed. The resulting global d13C distributions and gradients indicate that the deep water circulation was similar during the Holocene and LGM. This interpretation brings d13C data to a better agreement with Cd/Ca data and marks a sharp contrast with a widely held view based on d13C measurements that the glacial Southern Ocean was the terminus of the thermohaline circulation. The proposed presence of glacial North Atlantic Deep Water does not necessarily contradict the postulated presence of Glacial North Atlantic Intermediate Water.

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.

Biomarker and XRD results for the last 3.4 Ma from IODP Site 306-U1313

Various types of abrupt/millennial-scale climate variability such as Dansgaard/Oeschger and Heinrich Events characterized the last glacial period. Over the last decade, a number of studies demonstrated that such millennial-scale climate variability was not limited to the last glacial but inherent to Quaternary climate. Here we review the occurrence and origin of millennial ice-rafting events in the North Atlantic during the late Pliocene and Pleistocene (last 3.4 Ma) with a special focus on North Atlantic Hudson Strait (HS) Heinrich(-like) Events. Besides a clear biomarker signature, we show that Heinrich Layers 5, 4, 2, and 1 in marine sediment cores from across the North Atlantic all bear the organic geochemical fingerprint of the Hudson area. Using this framework and combining previously published results, detailed investigations into the organic and inorganic chemistry of ice-rafted debris (IRD) found across the North Atlantic demonstrate that prior to MIS 16 (~ 650 ka) IRD in the North Atlantic did not originate from the Hudson area of northern Canada. The signature of this early IRD is distinctly different compared to that of HS Heinrich Layers. Rather ice-rafting events during the late Pliocene and early Pleistocene predominantly emanated from the calving of the Greenland and Fennoscandian ice sheets and possibly minor contributions from local ice streams from the North American and British ice sheets. Compared to North Atlantic HS Heinrich Events, these early Pleistocene IRD-events had a limited impact on surface water characteristics in the North Atlantic. North Atlantic HS Heinrich(-like) Events first occurred during MIS 16. At the same time, the dominant frequency in silicate-rich IRD accumulation shifted from the obliquity (41-ka) to a 100-ka frequency across the North Atlantic. Iceberg survivability or a change in iceberg trajectory likely did not control this change in IRD-regime. These results lend further support for the existing hypothesis that an increase in size (thickness) of the Laurentide ice sheet controls the occurrence of North Atlantic HS Heinrich Events, favoring an internal dynamic mechanism for their occurrence.

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.

Dust fluxes and iron fertilization in Holocene and Last Glacial Maximum climates

Mineral dust aerosols play a major role in present and past climates. To date, we rely on climate models for estimates of dust fluxes to calculate the impact of airborne micronutrients on biogeochemical cycles. Here we provide a new global dust flux data set for Holocene and Last Glacial Maximum (LGM) conditions based on observational data. A comparison with dust flux simulations highlights regional differences between observations and models. By forcing a biogeochemical model with our new data set and using this model's results to guide a millennial-scale Earth System Model simulation, we calculate the impact of enhanced glacial oceanic iron deposition on the LGM-Holocene carbon cycle. On centennial timescales, the higher LGM dust deposition results in a weak reduction of <10?ppm in atmospheric CO2 due to enhanced efficiency of the biological pump. This is followed by a further ~10?ppm reduction over millennial timescales due to greater carbon burial and carbonate compensation.