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.

Abrupt changes in high-latitude nutrient supply to the Atlantic during the last glacial cycle

The supply of nutrients to the low-latitude thermocline is largely controlled by intermediate-depth waters formed at the surface in the high southern latitudes. Silicic acid is an essential macronutrient for diatoms, which are responsible for a significant portion of marine carbon export production. Changes in ocean circulation, such as those observed during the last deglaciation, would influence the nutrient composition of the thermocline and, therefore, the relative abundance of diatoms in the low latitudes. Here we present the first record of the silicic acid content of the Atlantic over the last glacial cycle. Our results show that at intermediate depths of the South Atlantic, the silicic acid concentration was the same at the Last Glacial Maximum (LGM) as it is today, overprinted by high silicic acid pulses that coincided with abrupt changes in ocean and atmospheric circulation during Heinrich Stadials and the Younger Dryas. We suggest these pulses were caused by changes in intermediate water formation resulting from shifts in the subpolar hydrological cycle, with fundamental implications for the nutrient supply to the Atlantic.

Peak Last Glacial weathering intensity on the North American continent recorded by the authigenic Hf isotope composition of North Atlantic deep-sea sediments

We have retrieved radiogenic hafnium (Hf) isotope compositions (ɛHf) from authigenic Fe–Mn oxyhydroxides of deep northwest Atlantic sediments deposited over the past 26 ka to investigate the oceanic evidence of changes in dissolved weathering inputs from NE America during the last deglaciation. The extraction of seawater-derived Hf isotopic compositions from Fe–Mn oxyhydroxides is not a standard procedure. Comparisons between the Al/Hf ratios and Hf isotopic compositions of the chemically extracted authigenic phase on the one hand, and those of the corresponding detrital fractions on the other, provide evidence that the composition of past seawater has been reliably obtained for most sampled depths with our leaching procedures. This is endorsed most strongly by data for a sediment core from 4250 m water depth at the deeper Blake Ridge, for which consistent replicates were produced throughout. The Hf isotopic composition of the most recent sample in this core also closely matches that of nearby present day central North Atlantic seawater. Comparison with previously published seawater Nd and Pb isotope compositions obtained on the same cores shows that both Hf and Pb were released incongruently during incipient chemical weathering, but responded differently to the deglacial retreat of the Laurentide Ice Sheet. Hafnium was released more congruently during peak glacial conditions of the Last Glacial Maximum (LGM) and changed to typical incongruent interglacial ɛHf signatures either during or shortly after the LGM. This indicates that some zircon-derived Hf was released to seawater during the LGM. Conversely, there is no clear evidence for an increase in the influence of weathering of Lu-rich mineral phases during deglaciation, possibly since relatively unradiogenic Hf contributions from feldspar weathering were superimposed. While the authigenic Pb isotope signal in the same marine sediment samples traced peak chemical weathering rates on continental North America during the transition to the Holocene a similar incongruent excursion is notably absent in the Hf isotope record. The early change towards more radiogenic ɛHf in relation to the LGM may provide direct evidence for the transition from a cold-based to a warm-based Laurentide Ice Sheet on the Atlantic sector of North America.

Deep thermohaline Circulation in the low-latitude Atlantic during the Last Glacial

Present-day low-latitude eastern and western Atlantic basins are geochemically distinct below the sill depth of the Mid-Atlantic Ridge. While Antarctic Bottom Water (AABW) circulates freely in the western Atlantic, flow into the eastern Atlantic is restricted below 4 km which results in filling the abyssal depths of this basin with water of geochemical similarity to nutrient depleted North Atlantic Deep Water. Using carbon isotopes and Cd/Ca ratios in benthic foraminifera we reconstruct the geochemistry of these basins during the last glacial maximum. Results indicate that deep eastern and western Atlantic basins became geochemically identical during the last glacial. This was achieved by shoaling of the upper surface of AABW above the sill depth of the Mid-Atlantic Ridge, which allowed bottom waters in both basins to be filled with the same water mass. Although AABW became the dominant water mass in the deep eastern Atlantic basin during the glacial, Holocene-glacial delta13C-PO4 shifts in this basin are in Redfield proportions, unlike the disproportionate Holocene-glacial delta13C-PO4 shifts observed in the Southern Ocean. By examining the composition of deep and intermediate waters throughout the Atlantic, we show that this effect was induced by a change in gradient of the delta13C-PO4 deepwater mixing line during glacial times. Evidence from high-latitude planktonic data suggests that the change in gradient of the deepwater mixing line was brought about through a significant reduction in the thermodynamic effect on Southern Ocean surface waters. By using coupled delta13C-PO4 data to constrain the composition of end member water masses in the glacial Atlantic, we conclude that deep waters in the low-latitude glacial Atlantic were composed of a mixture of northern and southern source waters in a ratio of 1:3.

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

The reconstruction of the stable carbon isotope evolution in atmospheric CO2 (d13Catm ), 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 d13Catm before, during and after the Marine Isotope Stage 5.5 (155 000 to 105 000 years BP). 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 permil shift to heavier values between the mean d13Catm level in the Penultimate (~ 140 000 years BP) and Last Glacial Maximum (~ 22 000 years 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 years, but with different phasing and magnitudes. Furthermore, a 5000 years lag in the CO2 decline relative to EDC temperatures is confirmed during the glacial inception at the end of MIS 5.5 (120 000 years BP). Based on our isotopic data this lag can be explained by terrestrial carbon release and carbonate compensation.

Cape Lookout cold-water coral area, coral ages and environmental parameters for the last glacial cycle

Climatic and oceanographic changes, as occurring at a glacial-interglacial scale, may alter the environmental conditions needed for the development of prolific cold-water coral reefs and mounds. Studies constraining the temporal distribution of cold-water corals in the NE Atlantic suggested the cyclic changes of the Atlantic Meridional Overturning Circulation as the main driver for the development and dispersal of cold-water coral ecosystems. However, conclusions were hindered by lack of data from the NW Atlantic. Aiming to overcome this lack of data, the temporal occurrence of cold-water corals in the Cape Lookout area along the southeastern US margin was explored by U-series dating. Furthermore, the local influence of the regional water masses, namely the Gulf Stream, on cold-water coral proliferation and occurrence since the Last Glacial Maximum was examined. Results suggest that the occurrence of cold-water corals in the Cape Lookout area is restricted to interglacial periods, with corals being present during the last ~7 kyr and also during the Eemian (~125 ka). The reconstructed local environmental conditions suggest an offshore displacement of the Gulf Stream and increased influence from the Mid-Atlantic Bight shelf waters during the last glacial period. During the deglacial sea level rise, the Gulf Stream moved coastward providing present-day-like conditions to the surface waters. Nevertheless, present-day conditions at the ocean sea floor were not established before 7.5 cal ka BP once the ultimate demise of the Laurentide ice-sheet caused the final sea level rise and the displacement of the Gulf Stream to its present location. Occasional presence of the Gulf Stream over the site during the Mid- to Late Holocene coincides with enhanced bottom current strength and a slightly higher bottom water temperature, which are environmental conditions that are favorable for cold-water coral growth.

High-resolution multi-proxy record of the last glacial period from Lake Van, eastern Anatolian high plateau, Turkey

Detailed analyses of the Lake Van pollen, Ca/K ratio and stable oxygen isotope record allow the identification of millennial-scale vegetation and environmental changes in eastern Anatolia throughout the last glacial (~75-15 ka BP). The climate within the last glacial was cold and dry, with low arboreal pollen (AP) levels. The driest and coldest period corresponds to Marine Isotope Stage (MIS) 2 (~28-14.5 ka BP) dominated by the highest values of xerophytic steppe vegetation. Our high-resolution multi proxy record shows rapid expansions and contractions of tree populations that reflects variability in temperature and moisture availability. This rapid vegetation and environmental changes can be linked to the stadial-interstadial pattern of the Dansgaard-Oeschger (DO) events as recorded in the Greenland ice cores. Periods of reduced moisture availability were characterized by enhanced xerophytic species and high terrigenous input from the Lake Van catchment area. Furthermore, comparison with the marine realm reveals that the complex atmosphere-ocean interaction can be explained by the strength and position of the westerlies, which is responsible for the supply of humidity in eastern Anatolia. Influenced by diverse topography of the Lake Van catchment, larger DO interstadials (e.g. DO 19, 17-16, 14, 12 and 8) show the highest expansion of temperate species within the last glacial. However, Heinrich events (HE), characterized by highest concentrations of ice-rafted debris (IRD) in marine sediments, are identified in eastern Anatolia by AP values not lower and high steppe components not more abundant than during DO stadials. In addition, this work is a first attempt to establish a continuous microscopic charcoal record over the last glacial in the Near East, which documents an initial immediate response to millennial-scale climate and environmental variability and enables us to shed light on the history of fire activity during the last glacial.

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.