Timing and climatic imprint of Greenland interstadials recorded in stalagmites from Northern Turkey

A 50 kyr-long exceptionally well-dated and highly resolved stalagmite oxygen (δ 18O) and carbon (δ 13C) isotope record from Sofular Cave in northwestern Turkey helps to further improve the dating of Greenland Interstadials (GI) 1, and 3–12. Timing of most GI in the Sofular record is consistent within ±10 to 300 years with the “iconic” Hulu Cave record. Larger divergences (>500 years) between Sofular and Hulu are only observed for GI 4 and 7. The Sofular record differs from the most recent NGRIP chronology by up to several centuries, whereas age offsets do not increase systematically with depth. The Sofular record also reveals a rapid and sensitive climate and ecosystem response in the eastern Mediterranean to GI, whereas a phase lag of ∼100 years between climate and full ecosystem response is evident. Finally, results of spectral analyses of the Sofular isotope records do not support a 1,470-year pacing of GI.

The speliothem record of climate variability in Southern Arabia.

Uranium-series dated stalagmites from Oman indicate that pluvial conditions prevailed from 6.3 to 10.5, 78 to 82, 120 to 130, 180 to 200 and 300 to 330 kyr B.P.; all of these periods coincide with peak interglacials. Oxygen (δ18O) and hydrogen (δD) isotope ratios of speleothem calcite and fluid inclusions reveal the source of moisture and provide information on the amount of precipitation, respectively. δ18O and δD values of stalagmites deposited during peak interglacials vary between −8 and −4 ‰ (VPDB) and −53 and −20‰ (Vienna Standard Mean Ocean Water [VSMOW]) respectively, whereas modern stalagmites range from −2.6 to −1.1‰ in δ18O (VPDB) and −7.6 and −3.3‰ in δD (VSMOW), respectively. The growth and isotopic records indicate that during peak interglacial periods, the limit of the monsoon rainfall was shifted far north of its present location and each pluvial period was coinciding with an interglacial stage of the marine oxygen isotope record.

Pediastrum species as potential indicators of lake-level change in tropical South America

We explored the potential for using Pediastrum (Meyen), a genus of green alga commonly found in palaeoecological studies, as a proxy for lake-level change in tropical South America. The study site, Laguna La Gaiba (LLG) (17°45′S, 57°40′W), is a broad, shallow lake located along the course of the Paraguay River in the Pantanal, a 135,000-km2 tropical wetland located mostly in western Brazil, but extending into eastern Bolivia. Fourteen surface sediment samples were taken from LLG across a range of lake depths (2-5.2 m) and analyzed for Pediastrum. We found seven species, of which P. musteri (Tell et Mataloni), P. argentiniense (Bourr. et Tell), and P. cf. angulosum (Ehrenb.) ex Menegh. were identified as potential indicators of lake level. Results of the modern dataset were applied to 31 fossil Pediastrum assemblages spanning the early Holocene (12.0 kyr BP) to present to infer past lake level changes qualitatively. Early Holocene (12.0-9.8 kyr BP) assemblages do not show a clear signal, though abundance of P. simplex (Meyen) suggests relatively high lake levels. Absence of P. musteri, characteristic of deep, open water, and abundance of macrophyte-associated taxa indicate lake levels were lowest from 9.8 to 3.0 kyr BP. A shift to wetter conditions began at 4.4 kyr BP, indicated by the appearance of P. musteri, though inferred lake levels did not reach modern values until 1.4 kyr BP. The Pediastrum-inferred mid-Holocene lowstand is consistent with lower precipitation, previously inferred using pollen from this site, and is also in agreement with evidence for widespread drought in the South American tropics during the middle Holocene. An inference for steadily increasing lake level from 4.4 kyr BP to present is consistent with diatom-inferred water level rise at Lake Titicaca, and demonstrates coherence with the broad pattern of increasing monsoon strength from the late Holocene until present in tropical South America.

High-latitude climate sensitivity to ice-sheet forcing over the last 120kyr

Interpretation of ice-core records is currently limited by paucity of modelling at adequate temporal and spatial resolutions. Several key questions relate to mechanisms of polar amplification and inter-hemispheric coupling on glacial/interglacial timescales. Here, we present the first results from a large set of global ocean–atmosphere climate model ‘snap-shot’ simulations covering the last 120 000 years using the Hadley Centre climate model (HadCM3) at up to 1 kyr temporal resolution. Two sets of simulations were performed in order to examine the roles of orbit and greenhouse gases versus ice-sheet forcing of orbital-scale climate change. A series of idealised Heinrich events were also simulated, but no changes to aerosols or vegetation were prescribed. This paper focuses on high latitudes and inter-hemispheric linkages. The simulations reproduce polar temperature trends well compared to ice-core reconstructions, although the magnitude is underestimated. Polar amplification varies with obliquity, but this variability is dampened by including variations in land ice coverage, while the overall amplification factor increases. The relatively constant amplification of Antarctic temperatures (with ice-sheet forcing included) suggests it is possible to use Antarctic temperature reconstructions to estimate global changes (which are roughly half the magnitude). Atlantic Ocean overturning circulation varies considerably only with the introduction of Northern Hemisphere ice sheets, but only weakens in the North Atlantic in the deep glacial, when ocean–sea-ice feedbacks result in the movement of the region of deep convection to lower latitudes and with the introduction of freshwater to the surface North Atlantic in order to simulate Heinrich events.

Interhemispheric coupling, the West Antarctic Ice Sheet and warm Antarctic interglacials

Ice core evidence indicates that even though atmospheric CO2 concentrations did not exceed ~300 ppm at any point during the last 800 000 years, East Antarctica was at least ~3–4 °C warmer than preindustrial (CO2~280 ppm) in each of the last four interglacials. During the previous three interglacials, this anomalous warming was short lived (~3000 years) and apparently occurred before the completion of Northern Hemisphere deglaciation. Hereafter, we refer to these periods as "Warmer than Present Transients" (WPTs). We present a series of experiments to investigate the impact of deglacial meltwater on the Atlantic Meridional Overturning Circulation (AMOC) and Antarctic temperature. It is well known that a slowed AMOC would increase southern sea surface temperature (SST) through the bipolar seesaw and observational data suggests that the AMOC remained weak throughout the terminations preceding WPTs, strengthening rapidly at a time which coincides closely with peak Antarctic temperature. We present two 800 kyr transient simulations using the Intermediate Complexity model GENIE-1 which demonstrate that meltwater forcing generates transient southern warming that is consistent with the timing of WPTs, but is not sufficient (in this single parameterisation) to reproduce the magnitude of observed warmth. In order to investigate model and boundary condition uncertainty, we present three ensembles of transient GENIE-1 simulations across Termination II (135 000 to 124 000 BP) and three snapshot HadCM3 simulations at 130 000 BP. Only with consideration of the possible feedback of West Antarctic Ice Sheet (WAIS) retreat does it become possible to simulate the magnitude of observed warming.

An oceanic origin for the increase of atmospheric radiocarbon during the Younger Dryas

Variations in carbon-14 to carbon-12 ratio in the atmosphere (Δ14Catm) provide a powerful diagnostic for elucidating the timing and nature of geophysical and anthropological change. The (Atlantic) marine archive suggests a rapid Δ14Catm increase of 50‰ at the onset of the Younger Dryas (YD) cold reversal (12.9–11.7 kyr BP), which has not yet been satisfactorily explained in terms of magnitude or causal mechanism, as either a change in ocean ventilation or production rate. Using Earth-system model simulations and comparison of marine-based radiocarbon records from different ocean basins, we demonstrate that the YD Δ14Catm increase is smaller than suggested by the marine archive. This is due to changes in reservoir age, predominantly caused by reduced ocean ventilation.

Glacial-interglacial changes in dust deposition on the Chinese Loess Plateau

The Chinese Loess Plateau (CLP) contains an extensive record of aeolian deposition through multiple glacial–interglacial cycles. Independent chronologies based on pedostratigraphy, magnetic susceptibility, radiocarbon and luminescence dating were developed for 79 sites and used to estimate aeolian mass accumulation rates (MARs) for marine isotope stages 1–5. The regional median value of MAR for Stage 2 is 310 g/m2/yr compared to an estimate of 65 g/m2/yr for Stage 5. Estimated MARs from individual sites for Stage 2 are approximately 4.3 times greater than MARs for Stage 5 and 2.1 times greater than for Stage 1. MAR values at individual sites are consistently highest in the northwest and lowest in the southwest of the CLP during all marine isotope stages. MARs estimated on sections through loess terraces are consistently higher than MAR estimates at other sites, indicating that local recycling of loess material from exposed river valley deposits has been significant throughout the last 130 kyr. Although the spatial and temporal patterns in MAR are robust, there are uncertainties about the magnitude of these changes due to (a) lack of bulk density measurements and uncertainties in the chronologies for individual sites, (b) site and chronological biases in the sampling used to derive regional estimates, and (c) the unquantified nature of human impact on accumulation rates during the late Holocene. Nevertheless, the records from the CLP pose a number of challenges which could be addressed by numerical models of the palaeo-dust cycle.

Antarctic Peninsula ice sheet evolution during the Cenozoic Era

The Antarctic Peninsula region is currently undergoing rapid environmental change, resulting in the thinning, acceleration and recession of glaciers and the sequential collapse of ice shelves. It is important to view these changes in the context of long-term palaeoenvironmental complexity and to understand the key processes controlling ice sheet growth and recession. In addition, numerical ice sheet models require detailed geological data for tuning and testing. Therefore, this paper systematically and holistically reviews published geological evidence for Antarctic Peninsula Ice Sheet variability for each key locality throughout the Cenozoic, and brings together the prevailing consensus of the extent, character and behaviour of the glaciations of the Antarctic Peninsula region. Major contributions include a downloadable database of 186 terrestrial and marine calibrated dates; an original reconstruction of the LGM ice sheet; and a new series of isochrones detailing ice sheet retreat following the LGM. Glaciation of Antarctica was initiated around the Eocene/Oligocene transition in East Antarctica. Palaeogene records of Antarctic Peninsula glaciation are primarily restricted to King George Island, where glacigenic sediments provide a record of early East Antarctic glaciations, but with modification of far-travelled erratics by local South Shetland Island ice caps. Evidence for Neogene glaciation is derived primarily from King George Island and James Ross Island, where glaciovolcanic strata indicate that ice thicknesses reached 500–850 m during glacials. This suggests that the Antarctic Peninsula Ice Sheet draped, rather than drowned, the topography. Marine geophysical investigations indicate multiple ice sheet advances during this time. Seismic profiling of continental shelf-slope deposits indicates up to ten large advances of the Antarctic Peninsula Ice Sheet during the Early Pleistocene, when the ice sheet was dominated by 40 kyr cycles. Glacials became more pronounced, reaching the continental shelf edge, and of longer duration during the Middle Pleistocene. During the Late Pleistocene, repeated glacials reached the shelf edge, but ice shelves inhibited iceberg rafting. The Last Glacial Maximum (LGM) occurred at 18 ka BP, after which transitional glaciomarine sediments on the continental shelf indicate ice-sheet retreat. The continental shelf contains large bathymetric troughs, which were repeatedly occupied by large ice streams during Pleistocene glaciations. Retreat after the LGM was episodic in the Weddell Sea, with multiple readvances and changes in ice-flow direction, but rapid in the Bellingshausen Sea. The late Holocene Epoch was characterised by repeated fluctuations in palaeoenvironmental conditions, with associated glacial readvances. However, this has been subsumed by rapid warming and ice-shelf collapse during the twentieth century.

How well can we simulate past climates? Evaluating the models using global palaeoenvironmental datasets.

Global syntheses of palaeoenvironmental data are required to test climate models under conditions different from the present. Data sets for this purpose contain data from spatially extensive networks of sites. The data are either directly comparable to model output or readily interpretable in terms of modelled climate variables. Data sets must contain sufficient documentation to distinguish between raw (primary) and interpreted (secondary, tertiary) data, to evaluate the assumptions involved in interpretation of the data, to exercise quality control, and to select data appropriate for specific goals. Four data bases for the Late Quaternary, documenting changes in lake levels since 30 kyr BP (the Global Lake Status Data Base), vegetation distribution at 18 kyr and 6 kyr BP (BIOME 6000), aeolian accumulation rates during the last glacial-interglacial cycle (DIRTMAP), and tropical terrestrial climates at the Last Glacial Maximum (the LGM Tropical Terrestrial Data Synthesis) are summarised. Each has been used to evaluate simulations of Last Glacial Maximum (LGM: 21 calendar kyr BP) and/or mid-Holocene (6 cal. kyr BP) environments. Comparisons have demonstrated that changes in radiative forcing and orography due to orbital and ice-sheet variations explain the first-order, broad-scale (in space and time) features of global climate change since the LGM. However, atmospheric models forced by 6 cal. kyr BP orbital changes with unchanged surface conditions fail to capture quantitative aspects of the observed climate, including the greatly increased magnitude and northward shift of the African monsoon during the early to mid-Holocene. Similarly, comparisons with palaeoenvironmental datasets show that atmospheric models have underestimated the magnitude of cooling and drying of much of the land surface at the LGM. The inclusion of feedbacks due to changes in ocean- and land-surface conditions at both times, and atmospheric dust loading at the LGM, appears to be required in order to produce a better simulation of these past climates. The development of Earth system models incorporating the dynamic interactions among ocean, atmosphere, and vegetation is therefore mandated by Quaternary science results as well as climatological principles. For greatest scientific benefit, this development must be paralleled by continued advances in palaeodata analysis and synthesis, which in turn will help to define questions that call for new focused data collection efforts.

Potential role of vegetation feedback in the climate sensitivity of high-latitude regions: a case study at 6000 years B.P.

Previous climate model simulations have shown that the configuration of the Earth's orbit during the early to mid-Holocene (approximately 10–5 kyr) can account for the generally warmer-than-present conditions experienced by the high latitudes of the northern hemisphere. New simulations for 6 kyr with two atmospheric/mixed-layer ocean models (Community Climate Model, version 1, CCMl, and Global ENvironmental and Ecological Simulation of Interactive Systems, version 2, GENESIS 2) are presented here and compared with results from two previous simulations with GENESIS 1 that were obtained with and without the albedo feedback due to climate-induced poleward expansion of the boreal forest. The climate model results are summarized in the form of potential vegetation maps obtained with the global BIOME model, which facilitates visual comparisons both among models and with pollen and plant macrofossil data recording shifts of the forest-tundra boundary. A preliminary synthesis shows that the forest limit was shifted 100–200 km north in most sectors. Both CCMl and GENESIS 2 produced a shift of this magnitude. GENESIS 1 however produced too small a shift, except when the boreal forest albedo feedback was included. The feedback in this case was estimated to have amplified forest expansion by approximately 50%. The forest limit changes also show meridional patterns (greatest expansion in central Siberia and little or none in Alaska and Labrador) which have yet to be reproduced by models. Further progress in understanding of the processes involved in the response of climate and vegetation to orbital forcing will require both the deployment of coupled atmosphere-biosphere-ocean models and the development of more comprehensive observational data sets

Eemian and Holocene sea-surface conditions in the southern Black Sea: organic-walled dinoflagellate cyst record from core 22-GC3/3

In order to compare the sea-surface conditions in the Black Sea during the Holocene and Eemian, sapropelic parts of marine core 22-GC3 (42°13.53′N/36°29.55′E, 838 m water depth) were studied for organic-walled dinoflagellate cyst content. The record shows a change from freshwater/brackish assemblages (Pyxidinopsis psilata, Spiniferites cruciformis, and Caspidinium rugosum) to more marine assemblages (Lingulodinium machaerophorum and Spiniferites ramosus complex) during each interglacial, due to the inflow of saline Mediterranean water. The lacustrine–marine transitions in 22-GC3 occurred at ~ 8.3 cal kyr BP during the early Holocene and ~ 128 kyr BP during the early Eemian, slightly later compared to the onset of interglacial conditions on the adjacent land. Dinoflagellate cyst assemblages reveal higher sea-surface salinity (~ 28–30) (e.g. Spiniferites pachydermus, Bitectatodinium tepikiense, and Spiniferites mirabilis) around ~ 126.5–121 kyr BP in comparison to the Holocene (~ 15–20) as well as relatively high sea-surface temperature (e.g. Tuberculodinium vancampoae, S. pachydermus, and S. mirabilis) especially at ~ 127.6–125.3 kyr BP. Establishment of high sea-surface salinity during the Eemian correlates very well with reconstructed relatively high global sea-level and is explained as a combined effect of increased Mediterranean supply and high temperatures at the beginning of the last interglacial. The observed changes in the dinocyst record highlight the importance of nutrients for the composition of the Eemian and Holocene dinocyst assemblages.

Interhemispheric dynamics of the African rainbelt during the late Quaternary

The spatial pattern of precipitation variability in tropical and subtropical Africa over the late Quaternary has long been debated. Prevailing hypotheses variously infer (1) insolation-controlled asymmetry of wet phases between hemispheres, (2) symmetric contraction and expansion of the tropical rainbelt, and (3) independent control on moisture available in Southern Africa via sea surface temperatures in the Indian Ocean. In this study we use climate-model simulations covering the last glacial cycle (120 kyr) with HadCM3 and the multi-model ensembles from PMIP3 (the Palaeoclimate Model Intercomparison Project) to investigate the long-term behaviour of the African rainbelt, and test these simulations against existing empirical palaeohydrological records. Through regional model-data comparisons we find evidence for the validity of several hypotheses, with various proposed processes occurring concurrently but with different regional emphasis (e.g. asymmetric shifts at the seasonal extremes and symmetric expansions/ contractions towards West equatorial regions). Crucially, variations in rainfall are associated with multiple forcing mechanisms that vary in their dominance both spatially and temporally over the glacial cycle; an important consideration when interpreting and extrapolating from often relatively short palaeoenvironmental records.

Timing and nature of the penultimate deglaciation in a high alpine stalagmite from Switzerland

The timing and nature of the penultimate deglaciation, also known as Termination II (T-II), is subject of controversial discussions due to the scarcity of precisely-dated palaeoclimate records. Here we present a new precisely-dated and highly-resolved multi-proxy stalagmite record covering T-II from the high alpine Schafsloch Cave in Switzerland, an area where climate is governed by the North Atlantic. The inception of stalagmite growth at 137.4 ± 1.4 kyr before present (BP) indicates the presence of drip water and cave air temperatures of above 0°C, and is related to a climate-induced change in the thermal state (from cold- to warm-based) of the glacier above the cave. The cessation of stalagmite growth between 133.1 ± 0.7 and 131.9 ± 0.6 kyr BP is most likely related to distinct drop in temperature associated with Heinrich stadial 11. The resumption of stalagmite growth at 131.9 ± 0.6 kyr BP is accompanied by an abrupt increase in temperature and precipitation as indicated by distinct shifts in the oxygen and carbon isotopic composition as well as in trace element concentrations. The mid-point of T-II is around 131.8 ± 0.6 kyr BP in the Schafsloch Cave record is significantly earlier compared to the age of 129.1 ± 0.1 kyr BP in the Sanbao Cave record from China. The different ages between both records can be best explained by the competing effects of insolation and glacial boundary forcing on seasonality and snow cover extent in Eurasia.