Late-Glacial and early Holocene marine 14C reservoir ages off North Iceland (16-9 cal kyr BP)

Tephrochronological age models and 48 14C age determinations on molluscs and foraminifera (planktonic and benthic) are applied for the calculation of marine 14C reservoir age variability during a time period covering the Heinrich event H1 to early Holocene (16–9 cal kyr BP). Our data source consists of four high-resolution marine sediment cores (HM107-04, HM107-05, MD99-2271, MD99-2275) from the North Icelandic shelf. The marine reservoir age (ΔR) is found to be extremely variable, ranging from 385 to 1065 14C years. Extreme ΔR values occur at the end of H1, with values around 1000 14C years (~15 cal kyr BP), probably due to reduced northward flow of well-ventilated subtropical surface waters and a southward expansion of polar waters, as well as an expansion of sea ice limiting air-sea gas exchange. With the onset of the Bølling-Allerød interstadial, the ΔR values decrease towards 0 14C years suggesting a more vigorous North Atlantic Current and an active meridional overturning circulation system. During the Younger Dryas stadial, ΔR values are consistently around 700 14C years suggesting e renewed expansion of polar waters and a weakened meridional overtuning circulation. Interestingly, ΔR values remain high (~200 14C years) at the onset of the Holocene suggesting continued high influence of polar waters. Subsequently, ΔR values rapidly decrease to ~¬ 250 14C years around 11 cal kyr BP, indicating increased air-sea CO2 exchange with the coeval atmosphere. The ΔR values average around 0 14C years from around 10.5 to 9.0 cal kyr BP.

Marine04 marine radiocarbon age calibration, 0-26 cal kyr BP

New radiocarbon calibration curves, IntCal04 and Marine04, have been constructed and internationally ratified to replace the terrestrial and marine components of IntCal98. The new calibration data sets extend an additional 2000 yr, from 0-26 cal kyr BP (Before Present, 0 cal. BP = AD 1950), and provide much higher resolution, greater precision, and more detailed structure than IntCal98. For the Marine04 curve, dendrochronologically-dated tree-ring samples, converted with a box diffusion model to marine mixed-layer ages, cover the period from 0-10.5 call kyr BR Beyond 10.5 cal kyr BP, high-resolution marine data become available from foraminifera in varved sediments and U/Th-dated corals. The marine records are corrected with site-specific C-14 reservoir age information to provide a single global marine mixed-layer calibration from 10.5-26.0 cal kyr BR A substantial enhancement relative to IntCal98 is the introduction of a random walk model, which takes into account the uncertainty in both the calendar age and the C-14 age to calculate the underlying calibration curve (Buck and Blackwell, this issue). The marine data sets and calibration curve for marine samples from the surface mixed layer (Marine04) are discussed here. The tree-ring data sets, sources of uncertainty, and regional offsets are presented in detail in a companion paper by Reimer et al. (this issue).

IntCal04 terrestrial radiocarbon age calibration, 0-26 cal kyr BP

A new calibration curve for the conversion of radiocarbon ages to calibrated (cal) ages has been constructed and internationally ratified to replace ImCal98, which extended from 0-24 cal kyr BP (Before Present, 0 cal BP = AD 1950). The new calibration data set for terrestrial samples extends from 0-26 cal kyr BP, but with much higher resolution beyond 11.4 cal kyr BP than ImCal98. Dendrochronologically-dated tree-ring samples cover the period from 0-12.4 cal kyr BP. Beyond the end of the tree rings, data from marine records (corals and foraminifera) are converted to the atmospheric equivalent with a site-specific marine reservoir correction to provide terrestrial calibration from 12.4-26.0 cal kyr BP. A substantial enhancement relative to ImCal98 is the introduction of a coherent statistical approach based on a random walk model, which takes into account the uncertainty in both the calendar age and the C-14 age to calculate the underlying calibration curve (Buck and Blackwell, this issue). The tree-ring data sets, sources of uncertainty, and regional offsets are discussed here. The marine data sets and calibration curve for marine samples from the surface mixed layer (Marine 04) are discussed in brief, but details are presented in Hughen et al. (this issue a). We do not make a recommendation for calibration beyond 26 cal kyr BP at this time; however, potential calibration data sets are compared in another paper (van der Plicht et al., this issue).

Orbital and millenial-scale environmental changes between 64 and 25 ka BP recorded in Black Sea sediments

High-resolution pollen and dinoflagellate cyst records from sediment core M72/5-25-GC1 were used to reconstruct vegetation dynamics in northern Anatolia and surface conditions of the Black Sea between 64 and 20 ka BP. During this period, the dominance of Artemisia in the pollen record indicates a steppe landscape and arid climate conditions. However, the concomitant presence of temperate arboreal pollen suggests the existence of glacial refugia in northern Anatolia. Long-term glacial vegetation dynamics reveal two major arid phases ~64–55 and 40–32 ka BP, and two major humid phases ~54–45 and 28–20 ka BP, correlating with higher and lower summer insolation, respectively. Dansgaard–Oeschger (D–O) cycles are clearly indicated by the 25-GC1 pollen record. Greenland interstadials are characterized by a marked increase in temperate tree pollen, indicating a spread of forests due to warm/wet conditions in northern Anatolia, whereas Greenland stadials reveal cold and arid conditions as indicated by spread of xerophytic biomes. There is evidence for a phase lag of ~500 to 1500 yr between initial warming and forest expansion, possibly due to successive changes in atmospheric circulation in the North Atlantic sector. The dominance of Pyxidinopsis psilata and Spiniferites cruciformis in the dinocyst record indicates brackish Black Sea conditions during the entire glacial period. The decrease of marine indicators (marine dinocysts, acritarchs) at ~54 ka BP and increase of freshwater algae (Pediastrum, Botryococcus) from 32 to 25 ka BP reveals freshening of the Black Sea surface water. This freshening is possibly related to humid phases in the region, to connection between Caspian Sea and Black Sea, to seasonal freshening by floating ice, and/or to closer position of river mouths due to low sea level. In the southern Black Sea, Greenland interstadials are clearly indicated by high dinocyst concentrations and calcium carbonate content, as a result of an increase in primary productivity. Heinrich events show a similar impact on the environment in the northern Anatolia/Black Sea region as Greenland stadials.

Climate warming and vegetation response after Heinrich event 1 (16 700–16 000 cal yr BP) in Europe south of the Alps

Chironomids preserved in a sediment core from Lago di Origlio (416 m a.s.l.), a lake in the foreland of the Southern Swiss Alps, allowed quantitative reconstruction of Late Glacial and Early Holocene summer temperatures using a combined Swiss–Norwegian temperature inference model based on chironomid assemblages from 274 lakes. We reconstruct July air temperatures of ca. 10 °C between 17 300 and 16 000 cal yr BP, a rather abrupt warming to ca. 12.0 °C at ca. 16 500–16 000 cal yr BP, and a strong temperature increase at the transition to the Bølling/Allerød interstadial with average temperatures of about 14 °C. During the Younger Dryas and earliest Holocene similar temperatures are reconstructed as for the interstadial. The rather abrupt warming at 16 500–16 000 cal yr BP is consistent with sea-surface temperature as well as speleothem records, which indicate a warming after the end of Heinrich event 1 (sensu stricto) and before the Bølling/Allerød interstadial in southern Europe and the Mediterranean Sea. Pollen records from Origlio and other sites in southern Switzerland and northern Italy indicate an early reforestation of the lowlands 2000–1500 yr prior to the large-scale afforestation of Central Europe at the onset of the Bølling/Allerød period at ca. 14 700–14 600 cal yr BP. Our results suggest that these early afforestation processes in the formerly glaciated areas of northern Italy and southern Switzerland have been promoted by increasing temperatures.

Direct linking of Greenland and Antarctic ice cores at the Toba eruption (74 ka BP)

The Toba eruption that occurred some 74 ka ago in Sumatra, Indonesia, is among the largest volcanic events on Earth over the last 2 million years. Tephra from this eruption has been spread over vast areas in Asia, where it constitutes a major time marker close to the Marine Isotope Stage 4/5 boundary. As yet, no tephra associated with Toba has been identified in Greenland or Antarctic ice cores. Based on new accurate dating of Toba tephra and on accurately dated European stalagmites, the Toba event is known to occur between the onsets of Greenland interstadials (GI) 19 and 20. Furthermore, the existing linking of Greenland and Antarctic ice cores by gas records and by the bipolar seesaw hypothesis suggests that the Antarctic counterpart is situated between Antarctic Isotope Maxima (AIM) 19 and 20. In this work we suggest a direct synchronization of Greenland (NGRIP) and Antarctic (EDML) ice cores at the Toba eruption based on matching of a pattern of bipolar volcanic spikes. Annual layer counting between volcanic spikes in both cores allows for a unique match. We first demonstrate this bipolar matching technique at the already synchronized Laschamp geomagnetic excursion (41 ka BP) before we apply it to the suggested Toba interval. The Toba synchronization pattern covers some 2000 yr in GI-20 and AIM-19/20 and includes nine acidity peaks that are recognized in both ice cores. The suggested bipolar Toba synchronization has decadal precision. It thus allows a determination of the exact phasing of inter-hemispheric climate in a time interval of poorly constrained ice core records, and it allows for a discussion of the climatic impact of the Toba eruption in a global perspective. The bipolar linking gives no support for a long-term global cooling caused by the Toba eruption as Antarctica experiences a major warming shortly after the event. Furthermore, our bipolar match provides a way to place palaeo-environmental records other than ice cores into a precise climatic context.

A compilation of Western European terrestrial records 60–8 ka BP: towards an understanding of latitudinal climatic gradients

Terrestrial records of past climatic conditions, such as lake sediments and speleothems, provide data of great importance for understanding environmental changes. However, unlike marine and ice core records, terrestrial palaeodata are often not available in databases or in a format that is easily accessible to the non-specialist. As a consequence, many excellent terrestrial records are unknown to the broader palaeoclimate community and are not included in compilations, comparisons, or modelling exercises. Here we present a compilation of Western European terrestrial palaeo-records covering, entirely or partially, the 60–8-ka INTIMATE time period. The compilation contains 56 natural archives, including lake records, speleothems, ice cores, and terrestrial proxies in marine records. The compilation is limited to include records of high temporal resolution and/or records that provide climate proxies or quantitative reconstructions of environmental parameters, such as temperature or precipitation, and that are of relevance and interest to a broader community. We briefly review the different types of terrestrial archives, their respective proxies, their interpretation and their application for palaeoclimatic reconstructions. We also discuss the importance of independent chronologies and the issue of record synchronization. The aim of this exercise is to provide the wider palaeo-community with a consistent compilation of high-quality terrestrial records, to facilitate model-data comparisons, and to identify key areas of interest for future investigations. We use the compilation to investigate Western European latitudinal climate gradients during the deglacial period and, despite of poorly constrained chronologies for the older records, we summarize the main results obtained from NW and SW European terrestrial records before the LGM.

Quantitative summer and winter temperature reconstructions from pollen and chironomid data between 15 and 8 ka BP in the Baltic–Belarus area

New pollen based reconstructions of summer (May-to-August) and winter (December-to-February) temperatures between 15 and 8 ka BP along a S-N transect in the Baltic-Belarus (BB) area display trends in temporal and spatial changes in climate variability. These results are completed by two chironomid-based July mean temperature reconstructions. The magnitude of change compared with modern temperatures was more prominent in the northern part of BB area. The 4 C degrees winter and 2 C degrees summer warming at the start of GI-1 was delayed in the BB area and Lateglacial maximum temperatures were reached at ca 13.6 ka BP, being 4 C degrees colder than the modern mean. The Younger Dryas cooling in the area was 5 C degrees colder than present, as inferred by all proxies. In addition, our analyses show an early Holocene divergence in winter temperature trends with modern values reaching 1 ka earlier (10 ka BP) in southern BB compared to the northern part of the region (9 ka BP).

Results of modelled Arctic river discharge trends since 7 ka BP

The Arctic hydrological cycle throughout the Holocene is analyzed based on the results of transient simulations with the coupled atmosphere-ocean circulation model ECHO-G. The results suggest a ~ 2 % increase of mid-Holocene to preindustrial Arctic river discharges for the Eurasian continent. However, rivers of the North America Arctic realm show a moderate runoff decline of approximately 4 to 5 % for the same period. The total river discharge into the Arctic Ocean has remained at an approximately constant preindustrial level since the mid Holocene. The positive discharge trend within Eurasia is caused by a more rapid decrease in local net evaporation compared to a smaller decline in advected moisture and hence precipitation. This effect is neither recognized within the North American Arctic domain nor in the far eastern part of the Eurasian Arctic realm. A detailed comparison of these model findings with a variety of proxy studies is conducted. The collected proxy records show trends of continental surface temperatures and precipitation rates that are consistent with the simulations. A continuation of the transient Holocene runs for the 19th and 20th century with increased greenhouse gases indicates an increase of the total river influx into the Arctic Ocean of up to 7.6 %. The Eurasian river discharges increase by 7.5 %, the North American discharges by up to 8.4 %. The most rapid increases have been detected since the beginning of the 20th century. These results are corroborated by the observed rising of Arctic river discharges during the last century which is attributed to anthropogenic warming. The acceleration of the Arctic hydrological cycle in the 20th century is without precedence in the Holocene.

Planktonic foraminifera oxygen isotope and trace metal data ODP Hole 172-1056A, 46-54 ka BP

Salinity increase in the subtropical gyre system may have pre-conditioned the North Atlantic Ocean for a rapid return to stronger overturning circulation and high-latitude warming following meltwater events during the Last Glacial period. Here we investigate the Gulf Stream - subtropical gyre system properties over Dansgaard-Oeschger (DO) cycles 14 to 12, including Heinrich ice-rafting event 5. During the Holocene and Last Glacial Maximum a positive gradient in surface dwelling planktonic foraminifera d18O (Globigerinoides ruber) can be observed between the Gulf Stream and subtropical gyre, due to decreasing temperature, increasing salinity, and a change from summer to year-round occurrence of G. ruber. We assess whether this gradient was a common feature during stadial-interstadial climate oscillations of Marine Isotope Stage 3, by comparing existing G. ruber d18O from ODP Site 1060 (subtropical gyre location) and new data from ODP Site 1056 (Gulf Stream location) between 54 and 46 ka. Our results suggest that this gradient was largely absent during the period studied. During the major warm DO interstadials 14 and 12 we infer a more zonal and wider Gulf Stream, influencing both ODP Sites 1056 and 1060. A Gulf Stream presence during these major interstadials is also suggested by the large vertical d18O gradient between shallow dwelling planktonic foraminifera species, especially G. ruber, and the deep dwelling species Globorotalia inflata at site 1056, which we associate with strong summer stratification and Gulf Stream presence. A major reduction in this vertical d18O gradient from 51 ka until the end of Heinrich event 5 at 48.5 ka suggests site 1056 was situated within the subtropical gyre in this mainly cold period, from which we infer a migration of the Gulf Stream to a position nearer to the continental shelf, indicative of a narrower Gulf Stream with possibly reduced transport.

Reference datasets and pollen record from sediment core YAK, Yakumo, SW Hokkaido

This study reconstructs middle and late Holocene vegetation and climate dynamics in the Oshima Peninsula, SW Hokkaido, using the published method of biome reconstruction and modern analogue technique applied to the Yakumo pollen record (42°17'03''N, 140°15'34''E) spanning the last 5500 years. Two previously published matrices assigning Japanese plant/pollen taxa to the major vegetation types (biomes) are tested using a newly compiled dataset of 78 surface pollen spectra from Hokkaido. With both matrices showing strengths and weaknesses in reconstructing cool mixed and temperate deciduous forests of Hokkaido, the results suggest the necessity to consider the whole list of identified terrestrial pollen taxa for generating robust vegetation reconstructions for northern Japan. Applied to the fossil pollen data, both biome-reconstruction approaches demonstrate consistently that oak-dominated cool mixed forest spread in the study region between 5.5 and 3.6 cal ka BP and was subsequently replaced by beech-dominated temperate deciduous forest. The pollen-based climate reconstruction suggests this change in the vegetation composition was caused by a shift from cooler and drier than present climate to warmer and wetter, similar to modern conditions about 3.6 cal ka BP. Comparing the pollen-based reconstruction results with the published marine records from the NW Pacific, the reconstructed vegetation and climate dynamics can be satisfactorily explained by the greater role played by the warm Tsushima Current in the Sea of Japan and in the Tsugaru Strait during the middle and late Holocene. An increase in sea surface temperatures west and south of the study site would favour air temperature rise and moisture uptake and cause an increase in precipitation and snow accumulation in the western part of Hokkaido during the late Holocene.