Individual foraminifera variability from the western and eastern equatorial Pacific during the Late Holocene and Last Glacial Maximum

El Niño-Southern Oscillation (ENSO) is a major source of global interannual variability, but its response to climate change is uncertain. Paleoclimate records from the Last Glacial Maximum (LGM) provide insight into ENSO behavior when global boundary conditions (ice sheet extent, atmospheric partial pressure of CO2) were different from those today. In this work, we reconstruct LGM temperature variability at equatorial Pacific sites using measurements of individual planktonic foraminifera shells. A deep equatorial thermocline altered the dynamics in the eastern equatorial cold tongue, resulting in reduced ENSO variability during the LGM compared to the Late Holocene. These results suggest that ENSO was not tied directly to the east-west temperature gradient, as previously suggested. Rather, the thermocline of the eastern equatorial Pacific played a decisive role in the ENSO response to LGM climate.

LPJ-LMfire simulations for Europe for the Last Glacial Maximum and preindustrial control, link to netCDF files

This data sets contains LPJ-LMfire dynamic global vegetation model output covering Europe and the Mediterranean for the Last Glacial Maximum (LGM; 21 ka) and for a preindustrial control simulation (20th century detrended climate). The netCDF data files are time averages of the final 30 years of the model simulation. Each netCDF file contains four or five variables: fractional cover of 9 plant functional types (PFTs; cover), total fractional coverage of trees (treecover), population density of hunter-gatherers (foragerPD; only for the "people" simulations), fraction of the gridcell burned on 30-year average (burnedf), and vegetation net primary productivity (NPP). The model spatial resolution is 0.5-degrees For the LGM simulations, LPJ-LMfire was driven by the PMIP3 suite of eight GCMs for which LGM climate simulations were available. Also provided in this archive is the result of an LPJ-LMfire run that was forced by the average climate of all GCMs (the "GCM-mean" files), and the average of each of the individual LPJ-LMfire runs over the eight LGM scenarios individually (the "LPJ-mean" files). The model simulations are provided that include the influence of human presence on the landscape (the "people" files), and in a "world without humans" scenario (the "natural" files). Finally this archive contains the preindustrial reference simulation with and without human influence ("PI_reference_people" and "PI_reference_nat", respectively). There are therefore 22 netCDF files in this archive: 8 each of LGM simulations with and without people (total 16) and the "GCM mean" simulation (2 files) and the "LPJ mean" aggregate (2 files), and finally the two preindustrial "control" simulations ("PI"), with and without humans (2 files). In addition to the LPJ-LMfire model output (netCDF files), this archive also contains a table of arboreal pollen percent calculated from pollen samples dated to the LGM at sites throughout (lgmAP.txt), and a table containing the location of archaeological sites dated to the LGM (LGM_archaeological_site_locations.txt).

Last glacial maximum sediment record of the southern Weddell Sea shelf

Abstract: The history of grounded ice-sheet extent on the southern Weddell Sea shelf during the Last Glacial Maximum (LGM) and the timing of post-LGM ice-sheet retreat are poorly constrained. Several glaciological models reconstructed widespread grounding and major thickening of the Antarctic Ice Sheet in the Weddell Sea sector at the LGM. In contrast, recently published onshore data and modelling results concluded only very limited LGM-thickening of glaciers and ice streams feeding into the modern Filchner and Ronne ice shelves. These studies concluded that during the LGM ice shelves rather than grounded ice covered the Filchner and Ronne troughs, two deep palaeo-ice stream troughs eroded into the southern Weddell Sea shelf. Here we review previously published and unpublished marine geophysical and geological data from the southern Weddell Sea shelf. The stratigraphy and geometry of reflectors in acoustic sub-bottom profiles are similar to those from other West Antarctic palaeo-ice stream troughs, where grounded ice had advanced to the shelf break at the LGM. Numerous cores from the southern Weddell Sea shelf recovered sequences with properties typical for subglacially deposited tills or subglacially compacted sediments. These data sets give evidence that grounded ice had advanced across the shelf during the past, thereby grounding in even the deepest parts of the Filchner and Ronne troughs. Radiocarbon dates from glaciomarine sediments overlying the subglacial deposits are limited, but indicate that the ice grounding occurred at the LGM and that ice retreat started before ~15.1 corrected 14C kyrs before present (BP) on the outer shelf and before ~7.7 corrected 14C kyrs BP on the inner shelf, which is broadly synchronous with ice retreat in other Antarctic sectors. The apparent mismatch between the ice-sheet reconstructions from marine and terrestrial data can be attributed to ice streams with very low surface profiles (similar to those of "ice plains") that had advanced through Filchner Trough and Ronne Trough at the LGM. Considering the global sea-level lowstand of ~130 metres below present, a low surface slope of the expanded LGM-ice sheet in the southern Weddell Sea can reconcile grounding-line advance to the shelf break with limited thickening of glaciers and ice streams in the hinterland. This scenario implies that ice-sheet growth in the Weddell Sea sector during the LGM and ice-sheet drawdown throughout the last deglaciation could only have made minor contributions to the major global sea-level fluctuations during these times.

Sedimentary record of the western Amundsen Sea Embayment after the last glacial maximum

Reliable dating of glaciomarine sediments deposited on the Antarctic shelf since the Last Glacial Maximum (LGM) is very challenging because of the general absence of calcareous (micro-) fossils and the recycling of fossil organic matter. As a consequence, radiocarbon (14C) ages of the acid-insoluble organic fraction (AIO) of the sediments bear uncertainties that are very difficult to quantify. In this paper we present the results of three different chronostratigraphic methods to date a sedimentary unit consisting of diatomaceous ooze and diatomaceous mud that was deposited following the last deglaciation at five core sites on the inner shelf in the western Amundsen Sea (West Antarctica). In three cores conventional 14C dating of the AIO in bulk sediment samples yielded age reversals down-core, but at all sites the AIO 14C ages obtained from diatomaceous ooze within the diatom-rich unit yielded similar uncorrected 14C ages ranging from 13,517±56 to 11,543±47 years before present (yr BP). Correction of these ages by subtracting the core-top ages, which are assumed to reflect present-day deposition (as indicated by 21044 Pb dating of the sediment surface at one core site), yielded ages between ca. 10,500 and 8,400 calibrated years before present (cal yr BP). Correction of the AIO ages of the diatomaceous ooze by only subtracting the marine reservoir effect (MRE) of 1,300 years indicated deposition of the diatom-rich sediments between 14,100 and 11,900 cal yr BP. Most of these ages are consistent with age constraints between 13.0 and 8.0 ka BP for the diatom-rich unit, which we obtained by correlating the relative palaeomagnetic intensity (RPI) records of three of the sediment cores with global and regional reference curves for palaeomagnetic intensity. As a third dating technique we applied conventional 53 radiocarbon dating of the AIO included in acid-cleaned diatom hard parts that were extracted from the diatomaceous ooze. This method yielded uncorrected 14C ages of only 5,111±38 and 5,106±38 yr BP, respectively. We reject these young ages, because they are likely to be overprinted by the adsorption of modern atmospheric carbon dioxide onto the surfaces of the extracted diatom hard parts prior to sample graphitisation and combustion for 14C dating. The deposition of the diatom-rich unit in the western Amundsen Sea suggests deglaciation of the inner shelf before ca. 13 ka BP. The deposition of diatomaceous oozes on other parts of the Antarctic shelf around the same time, however, seems to be coincidental rather than directly related.

Sea-surface temperature and sea ice distribution of the Southern Ocean at the EPILOG Last Glacial Maximum

Based on the quantitative study of diatoms and radiolarians, summer sea-surface temperature (SSST) and sea ice distribution were estimated from 122 sediment core localities in the Atlantic, Indian and Pacific sectors of the Southern Ocean to reconstruct the last glacial environment at the EPILOG (19.5-16.0 ka or 23 000-19 000 cal yr. B.P.) time-slice. The statistical methods applied include the Imbrie and Kipp Method, the Modern Analog Technique and the General Additive Model. Summer SSTs reveal greater surface-water cooling than reconstructed by CLIMAP (Geol. Soc. Am. Map Chart. Ser. MC-36 (1981) 1), reaching a maximum (4-5 °C) in the present Subantarctic Zone of the Atlantic and Indian sector. The reconstruction of maximum winter sea ice (WSI) extent is in accordance with CLIMAP, showing an expansion of the WSI field by around 100% compared to the present. Although only limited information is available, the data clearly show that CLIMAP strongly overestimated the glacial summer sea ice extent. As a result of the northward expansion of Antarctic cold waters by 5-10° in latitude and a relatively small displacement of the Subtropical Front, thermal gradients were steepened during the last glacial in the northern zone of the Southern Ocean. Such reconstruction may, however, be inapposite for the Pacific sector. The few data available indicate reduced cooling in the southern Pacific and give suggestion for a non-uniform cooling of the glacial Southern Ocean.

Dissolution index of Globigerina bulloides in recent and Last Glacial Maximum sediments

The modern Atlantic Ocean, dominated by the interactions of North Atlantic Deep Water (NADW) and Antarctic Bottom Water (AABW), plays a key role in redistributing heat from the Southern to the Northern Hemisphere. In order to reconstruct the evolution of the relative importance of these two water masses, the NADW/AABW transition, reflected by the calcite lysocline, was investigated by the Globigerina bulloides dissolution index (BDX?). The depth level of the Late Glacial Maximum (LGM) calcite lysocline was elevated by several hundred metres, indicating a more corrosive water mass present at modern NADW level. Overall, the small range of BDX? data and the gradual decrease in preservation below the calcite lysocline point to a less stratified Atlantic Ocean during the LGM. Similar preservation patterns in the West and East Atlantic demonstrate that the modern west-east asymmetry did not exist due to an expansion of southern deep waters compensating for the decrease in NADW formation.

Holocene and Last Glacial Maximum Atlantic neodymium isotopes

Neodymium isotopes measured on chemically uncleaned planktic foraminifera from cores throughout the Atlantic Ocean. Samples are Holocene and Last Glacial Maximum in age.

North Atlantic Deep Water Production during the Last Glacial Maximum

Changes in deep ocean ventilation are commonly invoked as the primary cause of lower glacial atmospheric CO2. The water mass structure of the glacial deep Atlantic Ocean and the mechanism by which it may have sequestered carbon remain elusive. Here we present neodymium isotope measurements from cores throughout the Atlantic that reveal glacial-interglacial changes in water mass distributions. These results demonstrate the sustained production of North Atlantic Deep Water under glacial conditions, indicating that southern-sourced waters were not as spatially extensive during the Last Glacial Maximum as previously believed. We demonstrate that the depleted glacial delta C-13 values in the deep Atlantic Ocean cannot be explained solely by water mass source changes. A greater amount of respired carbon, therefore, must have been stored in the abyssal Atlantic during the Last Glacial Maximum. We infer that this was achieved by a sluggish deep overturning cell, comprised of well-mixed northern-and southern-sourced waters.

A reconstruction of atmospheric carbon dioxide and its stable carbon isotopic composition from the penultimate glacial maximum to the last glacial inception

The reconstruction of the stable carbon isotope evolution in atmospheric CO2 (δ13Catm), as archived in Antarctic ice cores, bears the potential to disentangle the contributions of the different carbon cycle fluxes causing past CO2 variations. Here we present a new record of δ13Catm before, during and after the Marine Isotope Stage 5.5 (155 000 to 105 000 yr BP). The dataset is archived on the data repository PANGEA® (www.pangea.de) under 10.1594/PANGAEA.817041. The record was derived with a well established sublimation method using ice from the EPICA Dome C (EDC) and the Talos Dome ice cores in East Antarctica. We find a 0.4‰ shift to heavier values between the mean δ13Catm level in the Penultimate (~ 140 000 yr BP) and Last Glacial Maximum (~ 22 000 yr BP), which can be explained by either (i) changes in the isotopic composition or (ii) intensity of the carbon input fluxes to the combined ocean/atmosphere carbon reservoir or (iii) by long-term peat buildup. Our isotopic data suggest that the carbon cycle evolution along Termination II and the subsequent interglacial was controlled by essentially the same processes as during the last 24 000 yr, but with different phasing and magnitudes. Furthermore, a 5000 yr lag in the CO2 decline relative to EDC temperatures is confirmed during the glacial inception at the end of MIS5.5 (120 000 yr BP). Based on our isotopic data this lag can be explained by terrestrial carbon release and carbonate compensation.

Sea surface temperature reconstruction for the middle Okinawa Trough during the last glacial-interglacial cycle using C-37 unsaturated alkenones

C-37 unsaturated alkenones were analyzed on a core retrieved from the middle Okinawa Trough. The calculated U-37(K') displays a trend generally parallel with those of the oxygen isotopic compositions of two planktonic foraminiferal species, Neogloboquadrina dutertrei and Globigerinoides sacculifer, suggesting that in this region, SST has varied in phase with global ice volume change since the last glacial -interglacial cycle. The U-37(K')-derived SST ranged from ca. 24.0 to 27.5 degrees C, with the highest value 27.5 degrees C occurring in marine isotope stage 5 and the lowest similar to 24.0 degrees C in marine isotope stage 2. This trend is consistent with the continental records from the East Asian monsoon domain and the marine records from the Equatorial Pacific. The deglacial increase of the U-37(K')-derived SST is similar to 2.4 degrees C from the Last Glacial Maximum to the Holocene. (c) 2007 Elsevier B.V. All rights reserved.

The last glacial cycle: transient simulations with an AOGCM

A number of transient climate runs simulating the last 120kyr have been carried out using FAMOUS, a fast atmosphere-ocean general circulation model (AOGCM). This is the first time such experiments have been done with a full AOGCM, providing a three-dimensional simulation of both atmosphere and ocean over this period. Our simulation thus includes internally generated temporal variability over periods from days to millennia, and physical, detailed representations of important processes such as clouds and precipitation. Although the model is fast, computational restrictions mean that the rate of change of the forcings has been increased by a factor of 10, making each experiment 12kyr long. Atmospheric greenhouse gases (GHGs), northern hemisphere ice sheets and variations in solar radiation arising from changes in the Earth's orbit are treated as forcing factors, and are applied either separately or combined in different experiments. The long-term temperature changes on Antarctica match well with reconstructions derived from ice-core data, as does variability on timescales longer than 10 kyr. Last Glacial Maximum (LGM) cooling on Greenland is reasonably well simulated, although our simulations, which lack ice-sheet meltwater forcing, do not reproduce the abrupt, millennial scale climate shifts seen in northern hemisphere climate proxies or their slower southern hemisphere counterparts. The spatial pattern of sea surface cooling at the LGM matches proxy reconstructions reasonably well. There is significant anti-correlated variability in the strengths of the Atlantic Meridional Overturning Circulation (AMOC) and the Antarctic Circumpolar Current (ACC) on timescales greater than 10kyr in our experiments. We find that GHG forcing weakens the AMOC and strengthens the ACC, whilst the presence of northern hemisphere ice-sheets strengthens the AMOC and weakens the ACC. The structure of the AMOC at the LGM is found to be sensitive to the details of the ice-sheet reconstruction used. The precessional component of the orbital forcing induces ~20kyr oscillations in the AMOC and ACC, whose amplitude is mediated by changes in the eccentricity of the Earth's orbit. These forcing influences combine, to first order, in a linear fashion to produce the mean climate and ocean variability seen in the run with all forcings.

Last glacial vegetation of northern Eurasia

In order to investigate the potential role of vegetation changes in megafaunal extinctions during the later part of the last glacial stage and early Holocene (42–10 ka BP), the palaeovegetation of northern Eurasia and Alaska was simulated using the LPJ-GUESS dynamic vegetation model. Palaeoclimatic driving data were derived from simulations made for 22 time slices using the Hadley Centre Unified Model. Modelled annual net primary productivity (aNPP) of a series of plant functional types (PFTs) is mapped for selected time slices and summarised for major geographical regions for all time slices. Strong canonical correlations are demonstrated between model outputs and pollen data compiled for the same period and region. Simulated aNPP values, especially for tree PFTs and for a mesophilous herb PFT, provide evidence of the structure and productivity of last glacial vegetation. The mesophilous herb PFT aNPP is higher in many areas during the glacial than at present or during the early Holocene. Glacial stage vegetation, whilst open and largely treeless in much of Europe, thus had a higher capacity to support large vertebrate herbivore populations than did early Holocene vegetation. A marked and rapid decrease in aNPP of mesophilous herbs began shortly after the Last Glacial Maximum, especially in western Eurasia. This is likely implicated in extinction of several large herbivorous mammals during the latter part of the glacial stage and the transition to the Holocene.

Mid-Holocene and glacial-maximum vegetation geography of the northern continents and Africa.

BIOME 6000 is an international project to map vegetation globally at mid-Holocene (6000 14C yr bp) and last glacial maximum (LGM, 18,000 14C yr bp), with a view to evaluating coupled climate-biosphere model results. Primary palaeoecological data are assigned to biomes using an explicit algorithm based on plant functional types. This paper introduces the second Special Feature on BIOME 6000. Site-based global biome maps are shown with data from North America, Eurasia (except South and Southeast Asia) and Africa at both time periods. A map based on surface samples shows the method’s skill in reconstructing present-day biomes. Cold and dry conditions at LGM favoured extensive tundra and steppe. These biomes intergraded in northern Eurasia. Northern hemisphere forest biomes were displaced southward. Boreal evergreen forests (taiga) and temperate deciduous forests were fragmented, while European and East Asian steppes were greatly extended. Tropical moist forests (i.e. tropical rain forest and tropical seasonal forest) in Africa were reduced. In south-western North America, desert and steppe were replaced by open conifer woodland, opposite to the general arid trend but consistent with modelled southward displacement of the jet stream. The Arctic forest limit was shifted slighly north at 6000 14C yr bp in some sectors, but not in all. Northern temperate forest zones were generally shifted greater distances north. Warmer winters as well as summers in several regions are required to explain these shifts. Temperate deciduous forests in Europe were greatly extended, into the Mediterranean region as well as to the north. Steppe encroached on forest biomes in interior North America, but not in central Asia. Enhanced monsoons extended forest biomes in China inland and Sahelian vegetation into the Sahara while the African tropical rain forest was also reduced, consistent with a modelled northward shift of the ITCZ and a more seasonal climate in the equatorial zone. Palaeobiome maps show the outcome of separate, independent migrations of plant taxa in response to climate change. The average composition of biomes at LGM was often markedly different from today. Refugia for the temperate deciduous and tropical rain forest biomes may have existed offshore at LGM, but their characteristic taxa also persisted as components of other biomes. Examples include temperate deciduous trees that survived in cool mixed forest in eastern Europe, and tropical evergreen trees that survived in tropical seasonal forest in Africa. The sequence of biome shifts during a glacial-interglacial cycle may help account for some disjunct distributions of plant taxa. For example, the now-arid Saharan mountains may have linked Mediterranean and African tropical montane floras during enhanced monsoon regimes. Major changes in physical land-surface conditions, shown by the palaeobiome data, have implications for the global climate. The data can be used directly to evaluate the output of coupled atmosphere-biosphere models. The data could also be objectively generalized to yield realistic gridded land-surface maps, for use in sensitivity experiments with atmospheric models. Recent analyses of vegetation-climate feedbacks have focused on the hypothesized positive feedback effects of climate-induced vegetation changes in the Sahara/Sahel region and the Arctic during the mid-Holocene. However, a far wider spectrum of interactions potentially exists and could be investigated, using these data, both for 6000 14C yr bp and for the LGM.