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.

Late Pleistocene oxygen isotope record of biognic silica

We determined the isotopic composition of oxygen in marine diatoms in eight deep-sea cores recovered from the Atlantic sector of the Southern Ocean. The analytical reproducibility and core-to-core consistency of the isotopic signal suggests that diatom delta18O can be used as a new paleocenographic tool to reconstruct past variations in surface water characteristics and to generate 18O -isotope-based stratigraphy for the Southern Ocean. The data indicate that diatom delta18O reflects sea surface temperature and seawater isotopic composition and that diatoms retain their isotopic signal on timescales of a least 430 ka. The delta18O analyses of different diatom assemblages reveal that the isotopic signal is free of species effects and that the common Antarctic species have the same water-opal fractionation. The transition from the last glacial maximum (LGM) to the Holocene is fully recorded in high sedimentation rate cores. An 18O enrichment during the LGM, a post-LGM meltwater spike and an input of meltwater during the late Holocene are the main isotopic features observed in down core records. The origin of this meltwater was very likely melting icebergs and/or continental ice or by melting sea ice that had accumulated snow. The most pronounced meltwater effects are recorded in cores that are associated with the Weddel gyre. Our results provide the basis for extending isotope studies to oceanic regions devoid of carbonate; further, isotopic stratigraphies may be constructed for records and regions where they were previously not possible.

Micropaleontological and geochemical proxy studies of sediment core GeoB7926-2

[1] The low-latitude upwelling regime off the Mauritanian coast in the subtropical NE Atlantic accounts for a significant part of global export production. Although productivity variations in coastal upwelling areas are usually attributed to changes in wind stress and upwelling intensity, productivity dynamics off Mauritania are less straightforward because of the complex atmospheric and hydrographic setting. Here we integrate micropaleontological (diatoms) and geochemical (bulk biogenic sediment components, X-ray fluorescence, and alkenones) proxies to examine on submillennial-to-millennial changes in diatom production that occurred off Mauritania, NW Africa, for the last 25 ka. During the Last Glacial Maximum (LGM, 19.0-23.0 ka B.P.), moderate silicate content of upwelled waters coupled with weakened NE trade winds determined moderate diatom productivity. No significant cooling is observed during the LGM, suggesting that our alkenone-based SST reconstruction represents a local, upwelling-related signal rather than a global insolation related one. Extraordinary increases in diatom and opal concentrations during Heinrich event 1 (H1, 15.5-18.0 ka B.P.) and the Younger Dryas (YD, 13.5-11.5 ka B.P.) are attributed to enhanced upwelling of silica-rich waters and an enlarged upwelling filament, due to more intense NE trade winds. The synchronous increase of CaCO3 and K intensity and the decreased opal and diatoms values mark the occurrence of the Bølling/Allerød (BA, 13.5-15.5 ka B.P.) due to weakened eolian input and more humid conditions on land. Although the high export of diatoms is inextricably linked to upwelling intensity off Mauritania, variability in the nutrient content of the thermocline also plays a decisive role.

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.

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.

Averages and ranges of paleoenvironmental proxy records of ODP Hole 108-658C (Table 1)

Several high-resolution proxy environmental records have been obtained for the last 35 kyr from ODP Hole 658C, a well-studied site ca. 200 km off Cap Blanc, NW Africa. The collective assessment based on the marine proxies (UK'37 SST, contents of TOC and chlorins, Upwelling Radiolarian Index and the percentage of Florisphaera profunda), surprisingly indicates that the last glacial maximum (LGM) was characterized by warmer sea surface temperature (SST), weaker upwelling, and lower marine productivity, compared with the preceding older glacial and subsequent deglaciation periods. Of the terrigenous proxies, the mean grain size of the non-carbonate fraction and the terrigenous alkane content indicate that wind strength and aridity were high. The weaker upwelling at the 658 site during the LGM may have resulted from changes in the strength and direction of the wind systems and/or shifts in the position and geometry of the upwelling cell.

Coccolithophore taxa and PCA factor analysis of ODP Site 1233

We provide high-resolution sea surface temperature (SST) and paleoproductivity data focusing on Termination 1. We describe a new method for estimating SSTs based on multivariate statistical analyses performed on modern coccolithophore census data, and we present the first downcore reconstructions derived from coccolithophore assemblages at Ocean Drilling Project (ODP) Site 1233 located offshore Chile. We compare our coccolithophore SST record to alkenone-based SSTs as well as SST reconstructions based on dinoflagellates and radiolaria. All reconstructions generally show a remarkable concordance. As in the alkenone SST record, the Last Glacial Maximum (LGM, 19-23 kyr B.P.) is not clearly defined in our SST reconstruction. After the onset of deglaciation, three major warming steps are recorded: from 18.6 to 18 kyr B.P. (~2.6°C), from 15.7 to 15.3 kyr B.P. (~2.5°C), and from 13 to 11.4 kyr B.P. (~3.4°C). Consistent with the other records from Site 1233 and Antarctic ice core records, we observed a clear Holocene Climatic Optimum (HCO) from ~8-12 kyr B.P. Combining the SST reconstruction with coccolith absolute abundances and accumulation rates, we show that colder temperatures during the LGM are linked to higher coccolithophore productivity offshore Chile and warmer SSTs during the HCO to lower coccolithophore productivity, with indications of weak coastal upwelling. We interpret our data in terms of latitudinal displacements of the Southern Westerlies and the northern margin of the Antarctic Circumpolar Current system over the deglaciation and the Holocene.

(Appendix table DR1) Foraminiferal distribution data for ODP Hole 178-1098B

Benthic foraminiferal data from Ocean Drilling Program Site 1098 indicate significant changes in deep-water conditions of the Palmer Deep, western Antarctic Peninsula margin, throughout the Holocene (13 ka to present). The earliest Holocene represents a period of transition from the Last Glacial Maximum (LGM). Cold bottom waters, similar to saline shelf water (SSW), dominated the middle Holocene. The late Holocene in the Palmer Deep has been characterized by alternating dominance of circumpolar deep water (CDW) and saline shelf water. These changes have global oceanographic and climatic implications. We suggest that the middle Holocene bottom-water record, in the absence of circumpolar deep water on the western Antarctic Peninsula shelf, indicates high saline shelf water production and/or weakened circumpolar deep water production during the middle Holocene climatic optimum. The late Holocene benthic foraminiferal record indicates rapidly fluctuating sea-ice conditions and may indicate a teleconnection between the South Pacific and Southern Ocean, thus having implications related to the Southern Oscillation Index.

A Late Pleistocene sea level stack

Late Pleistocene sea level has been reconstructed from ocean sediment core data using a wide variety of proxies and models. However, the accuracy of individual reconstructions is limited by measurement error, local variations in salinity and temperature, and assumptions particular to each technique. Here we present a sea level stack (average) which increases the signal-to-noise ratio of individual reconstructions. Specifically, we perform principal component analysis (PCA) on seven records from 0-430 ka and five records from 0-798 ka. The first principal component, which we use as the stack, describes ~80 % of the variance in the data and is similar using either five or seven records. After scaling the stack based on Holocene and Last Glacial Maximum (LGM) sea level estimates, the stack agrees to within 5 m with isostatically adjusted coral sea level estimates for Marine Isotope Stages 5e and 11 (125 and 400 ka, respectively). When we compare the sea level stack with the d18O of benthic foraminifera, we find that sea level change accounts for about ~40 % of the total orbital-band variance in benthic d18O, compared to a 65 % contribution during the LGM-to-Holocene transition. Additionally, the second and third principal components of our analyses reflect differences between proxy records associated with spatial variations in the d18O of seawater.

Late Weichselian fluvial evolution of the southern Kara Sea

Glaciations had a profound impact on the global sea-level and particularly on the Arctic environments. One of the key questions related to this topic is, how did the discharge of the Siberian Ob and Yenisei rivers interact with a proximal ice sheet? In order to answer this question high-resolution (1-12 kHz), shallow-penetration seismic profiles were collected on the passive continental margin of the Kara Sea Shelf to study the paleo-drainage pattern of the Ob and Yenisei rivers. Both rivers incised into the recent shelf, leaving filled and unfilled river channels and river canyons/valleys connecting to a complex paleo-drainage network. These channels have been subaerially formed during a regressive phase of the global sea-level during the Last Glacial Maximum. Beyond recent shelf depths of 120 m particle transport is manifested in submarine channel-levee complexes acting as conveyor for fluvial-derived fines. In the NE area, uniform draping sediments are observed. Major morphology determining factors are (1) sea-level fluctuations and (2) LGM ice sheet influence. Most individual channels show geometries typical for meandering rivers and appear to be an order of magnitude larger than recent channel profiles of gauge stations on land. The Yenisei paleo-channels have larger dimensions than the Ob examples and could be originated by additional water release during the melt of LGM Putoran ice masses. Asymmetrical submarine channel-levee complexes with channel depths of 60 m and more developed, in some places bordered by glacially dominated morphology, implying deflection by the LGM ice masses. A total of more than 12,000 km of acoustic profiles reveal no evidence for an ice-dammed lake of greater areal extent postulated by several workers. Furthermore, the existence of the channel-levee complexes is indicative of unhindered sediment flow to the north. Channels situated on the shelf above 120-m water depth exhibit no phases of ponding and or infill during sea-level lowstand. These findings denote the non-existence of an ice sheet on large areas of the Kara Sea shelf.

Glacial-interglacial sediment record from the North Atlantic

Studies from the subtropical western and eastern Atlantic Ocean, using the 231Pa/230Th ratio as a kinematic proxy for deep water circulation, provided compelling evidence for a strong link between climate and the rate of meridional overturning circulation (MOC) over the last deglaciation. In this study, we present a compilation of existing and new sedimentary 231Pa/230Th records from North Atlantic cores between 1710 and 4550 m water depth. Comparing sedimentary 231Pa/230Th from different depths provides new insights into the evolution of the geometry and rate of deep water formation in the North Atlantic during the last 20,000 years. The 231Pa/230Th ratio measured in upper Holocene sediments indicates slow water renewal above ?2500 m and rapid flushing below, consistent with our understanding of modern circulation. In contrast, during the Last Glacial Maximum (LGM), Glacial North Atlantic Intermediate Water (GNAIW) drove a rapid overturning circulation to a depth of at least ?3000 m depth. Below ~4000 m, water renewal was much slower than today. At the onset of Heinrich event 1, transport by the overturning circulation declined at all depths. GNAIW shoaled above 3000 m and significantly weakened but did not totally shut down. During the Bølling-Allerød (BA) that followed, water renewal rates further decreased above 2000 m but increased below. Our results suggest for the first time that ocean circulation during that period was quite distinct from the modern circulation mode, with a comparatively higher renewal rate above 3000 m and a lower renewal rate below in a pattern similar to the LGM but less accentuated. MOC during the Younger Dryas appears very similar to BA down to 2000 m and slightly slower below.

Sedimentological investigations and age constraints of sediment cores from the outer Abbot trough, Amundsen Sea shelf

Knowing the extent of the West Antarctic Ice Sheet (WAIS) at the Last Glacial Maximum (LGM) is crucial for initiating and calibrating numerical ice sheet models that can predict future ice-sheet change and contributions to sea level. However, empirical data are lacking for key areas of outer continental shelves, where the LGM-WAIS must have terminated. We present detailed marine geophysical and geological data documenting an up to ~12 m-thick sequence of glaciomarine sediments within a relict glacial trough in the outer parts of the Amundsen Sea Embayment. Continuous deposition must have persisted here since at least >40 ka BP, pre-dating the established LGM by >13,000 years. Observations constrain the LGM grounding line to a distinct grounding-zone wedge ~100 km inland from the continental shelf edge. Thus, a substantial shelf area (~6000 km**2) remained ice free through the last glacial cycle.