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.

A tephra lattice for Greenland and a reconstruction of volcanic events spanning 25–45 ka b2k

Tephra layers preserved within the Greenland ice-cores are crucial for the independent synchronisation of these high-resolution records to other palaeoclimatic archives. Here we present a new and detailed tephrochronological framework for the time period 25,000 e 45,000 a b2k that brings together results from 4 deep Greenland ice-cores. In total, 99 tephra deposits, the majority of which are preserved as cryptotephra, are described from the NGRIP, NEEM, GRIP and DYE-3 records. The major element signatures of single glass shards within these deposits indicate that 93 are basaltic in composition all originating from Iceland. Specifically, 43 originate from Grimsv € otn, 20 are thought to be sourced from the Katla volcanic system and 17 show affinity to the Kverkfj € oll system. Robust geochemical characterisations, independent ages derived from the GICC05 ice-core chronology, and the stratigraphic positions of these deposits relative to the Dansgaard-Oeschger climate events represent a key framework that provides new information on the frequency and nature of volcanic events in the North Atlantic region between GS-3 and GI-12. Of particular importance are 19 tephra deposits that lie on the rapid climatic transitions that punctuate the last glacial period. This framework of well-constrained, time-synchronous tie-lines represents an important step towards the independent synchronisation of marine, terrestrial and ice-core records from the North Atlantic region, in order to assess the phasing of rapid climatic changes during the last glacial period.

Doppler Orbit Determination of Deep Space Probes by the Bernese GNSS Software: First Results of the Combined Orbit Determination from DSN and Inter-Satellite Ka-Band Data from the Grail Mission

Navigation of deep space probes is most commonly operated using the spacecraft Doppler tracking technique. Orbital parameters are determined from a series of repeated measurements of the frequency shift of a microwave carrier over a given integration time. Currently, both ESA and NASA operate antennas at several sites around the world to ensure the tracking of deep space probes. Just a small number of software packages are nowadays used to process Doppler observations. The Astronomical Institute of the University of Bern (AIUB) has recently started the development of Doppler data processing capabilities within the Bernese GNSS Software. This software has been extensively used for Precise Orbit Determination of Earth orbiting satellites using GPS data collected by on-board receivers and for subsequent determination of the Earth gravity field. In this paper, we present the currently achieved status of the Doppler data modeling and orbit determination capabilities in the Bernese GNSS Software using GRAIL data. In particular we will focus on the implemented orbit determination procedure used for the combined analysis of Doppler and intersatellite Ka-band data. We show that even at this earlier stage of the development we can achieve an accuracy of few mHz on two-way S-band Doppler observation and of 2 µm/s on KBRR data from the GRAIL primary mission phase.

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).

Central European vegetation response to abrupt climate change at 8.2 ka

Oxygen isotope records show a major climatic reversal at 8.2 ka in Greenland and Europe. Annually laminated sediments from two lakes in Switzerland and Germany were sampled contiguously to assess the response of European vegetation to climate change ca. 8.2 ka with time resolution and precision comparable to those of the Greenland ice cores. The pollen assemblages show pronounced and immediate responses (0–20 yr) of terrestrial vegetation to the climatic change at 8.2 ka. A sudden collapse of Corylus avellana (hazel) was accompanied by the rapid expansion of Pinus (pine), Betula (birch), and Tilia (linden), and by the invasion of Fagus silvatica (beech) and Abies alba (fir). Vegetational changes suggest that climatic cooling reduced drought stress, allowing more drought-sensitive and taller growing species to out-compete Corylus avellana by forming denser forest canopies. Climate cooling at 8.2 ka and the immediate reorganization of terrestrial ecosystems has gone unrecognized by previous pollen studies. On the basis of our data we conclude that the early Holocene high abundance of C. avellana in Europe was climatically caused, and we question the conventional opinion that postglacial expansions of F. silvatica and A. alba were controlled by low migration rates rather than by climate. The close connection between climatic change and vegetational response at a subcontinental scale implies that forecasted global warming may trigger rapid collapses, expansions, and invasions of tree species.

Sedimentological processes and environmental variability at Lake Ohrid (Macedonia, Albania) between 637 ka and the present

Lake Ohrid (Macedonia, Albania) is thought to be more than 1.2 million years old and host more than 300 endemic species. As a target of the International Continental scientific Drilling Program (ICDP), a successful deep drilling campaign was carried out within the scope of the Scientific Collaboration on Past Speciation Conditions in Lake Ohrid (SCOPSCO) project in 2013. Here, we present lithological, sedimentological, and (bio-)geochemical data from the upper 247.8 m composite depth of the overall 569 m long DEEP site sediment succession from the central part of the lake. According to an age model, which is based on 11 tephra layers (first-order tie points) and on tuning of bio-geochemical proxy data to orbital parameters (second-order tie points), the analyzed sediment sequence covers the last 637 kyr. The DEEP site sediment succession consists of hemipelagic sediments, which are interspersed by several tephra layers and infrequent, thin (< 5 cm) mass wasting deposits. The hemipelagic sediments can be classified into three different lithotypes. Lithotype 1 and 2 deposits comprise calcareous and slightly calcareous silty clay and are predominantly attributed to interglacial periods with high primary productivity in the lake during summer and reduced mixing during winter. The data suggest that high ion and nutrient concentrations in the lake water promoted calcite precipitation and diatom growth in the epilimnion during MIS15, 13, and 5. Following a strong primary productivity, highest interglacial temperatures can be reported for marine isotope stages (MIS) 11 and 5, whereas MIS15, 13, 9, and 7 were comparably cooler. Lithotype 3 deposits consist of clastic, silty clayey material and predominantly represent glacial periods with low primary productivity during summer and longer and intensified mixing during winter. The data imply that the most severe glacial conditions at Lake Ohrid persisted during MIS16, 12, 10, and 6, whereas somewhat warmer temperatures can be inferred for MIS14, 8, 4, and 2. Interglacial-like conditions occurred during parts of MIS14 and 8.

Northern Mediterranean climate since the Middle Pleistocene: a 637 ka stable isotope record from Lake Ohrid (Albania/Macedonia)

Lake Ohrid (Macedonia/Albania) is an ancient lake with unique biodiversity and a site of global significance for investigating the influence of climate, geological, and tectonic events on the generation of endemic populations. Here, we present oxygen (δ18O) and carbon (δ13C) isotope data from carbonate over the upper 243 m of a composite core profile recovered as part of the Scientific Collaboration on Past Speciation Conditions in Lake Ohrid (SCOPSCO) project. The investigated sediment succession covers the past ca. 637 ka. Previous studies on short cores from the lake (up to 15 m, < 140 ka) have indicated the total inorganic carbon (TIC) content of sediments to be highly sensitive to climate change over the last glacial–interglacial cycle. Sediments corresponding to warmer periods contain abundant endogenic calcite; however, an overall low TIC content in glacial sediments is punctuated by discrete bands of early diagenetic authigenic siderite. Isotope measurements on endogenic calcite (δ18Oc and δ13Cc) reveal variations both between and within interglacials that suggest the lake has been subject to palaeoenvironmental change on orbital and millennial timescales. We also measured isotope ratios from authigenic siderite (δ18Os and δ13Cs) and, with the oxygen isotope composition of calcite and siderite, reconstruct δ18O of lake water (δ18Olw) over the last 637 ka. Interglacials have higher δ18Olw values when compared to glacial periods most likely due to changes in evaporation, summer temperature, the proportion of winter precipitation (snowfall), and inflow from adjacent Lake Prespa. The isotope stratigraphy suggests Lake Ohrid experienced a period of general stability from marine isotope stage (MIS) 15 to MIS 13, highlighting MIS 14 as a particularly warm glacial. Climate conditions became progressively wetter during MIS 11 and MIS 9. Interglacial periods after MIS 9 are characterised by increasingly evaporated and drier conditions through MIS 7, MIS 5, and the Holocene. Our results provide new evidence for long-term climate change in the northern Mediterranean region, which will form the basis to better understand the influence of major environmental events on biological evolution within Lake Ohrid.

Environmental control on the occurrence of high-coercivity magnetic minerals and formation of iron sulfides in a 640 ka sediment sequence from Lake Ohrid (Balkans)

The bulk magnetic mineral record from Lake Ohrid, spanning the past 637 kyr, reflects large-scale shifts in hydrological conditions, and, superimposed, a strong signal of environmental conditions on glacial–interglacial and millennial timescales. A shift in the formation of early diagenetic ferrimagnetic iron sulfides to siderites is observed around 320 ka. This change is probably associated with variable availability of sulfide in the pore water. We propose that sulfate concentrations were significantly higher before  ∼  320 ka, due to either a higher sulfate flux or lower dilution of lake sulfate due to a smaller water volume. Diagenetic iron minerals appear more abundant during glacials, which are generally characterized by higher Fe / Ca ratios in the sediments. While in the lower part of the core the ferrimagnetic sulfide signal overprints the primary detrital magnetic signal, the upper part of the core is dominated by variable proportions of high- to low-coercivity iron oxides. Glacial sediments are characterized by high concentration of high-coercivity magnetic minerals (hematite, goethite), which relate to enhanced erosion of soils that had formed during preceding interglacials. Superimposed on the glacial–interglacial behavior are millennial-scale oscillations in the magnetic mineral composition that parallel variations in summer insolation. Like the processes on glacial–interglacial timescales, low summer insolation and a retreat in vegetation resulted in enhanced erosion of soil material. Our study highlights that rock-magnetic studies, in concert with geochemical and sedimentological investigations, provide a multi-level contribution to environmental reconstructions, since the magnetic properties can mirror both environmental conditions on land and intra-lake processes.

A first chironomid-based summer temperature reconstruction (13–5 ka BP) around 49°N in inland Europe compared with local lake development

Temperature reconstructions for the end of the Pleistocene and the first half of the Holocene based on biotic proxies are rare for inland Europe around 49°N. We analysed a 7 m long sequence of lake deposits in the Vihorlat Mts in eastern Slovakia (820 m a.s.l.). Chironomid head capsules were used to reconstruct mean July temperature (TJuly), other proxies (diatoms, green algae, pollen, geochemistry) were used to reconstruct local environmental changes that might have affected the climate reconstruction, such as epilimnetic total phosphorus concentrations (TP), lake level changes and development of surrounding vegetation. During the Younger Dryas (YD), temperature fluctuated between 7 and 11 °C, with distinct, decadal to centennial scale variations, that agree with other palaeoclimate records in Europe such as δ18O content in stalagmites or Greenland ice cores. The results indicate that the site was somewhat colder than expected from the general south-to-north YD temperature gradient within Europe, possibly because of north-facing exposition. The warmer phases of the YD were characterised by low water level or even complete desiccation of the lake (12,200-12,400 cal yr BP). At the Late-Glacial/Holocene transition TJuly steeply increased from from 11 to 15.5 °C (11,700-11,400 cal yr BP) - the highest TJuly for entire sequence. This rapid climate change was reflected by all proxies as a compositional change and increasing species diversity. The open woodlands of Pinus, Betula, Larix and Picea were replaced by broad-leaved temperate forests dominated by Betula, later by Ulmus and finally by Corylus (ca 9700 cal yr BP). At the same time, input of eroded coarse-grained material into the lake decreased and organic matter (LOI) and biogenic silica increased. The Early-Holocene climate was rather stable till 8700 cal yr BP, with temporary decrease in TJuly around 11,200 cal yr BP. The lake was productive with a well-developed littoral, as indicated by both diatoms and chironomids. A distinct decline of TJuly to 10 °C between 8700 and 8000 cal yr BP was associated with decreasing chironomid diversity and increasing climate moistening indicated by pollen. Tychoplanktonic and phosphorus-demanding diatoms increased which might be explained by hydrological and land-cover changes. Later, a gradual warming started after 7000 cal yr BP and representation of macrophytes, periphytic diatoms and littoral chironomids increased. Our results suggest that the Holocene thermal maximum was taking place unusually early in the Holocene at our study site, but its timing might be affected by topography and mesoclimate. We further demonstrated that temperature changes had coincided with variations in local hydrology

Long term low latitude and high elevation cosmogenic 3He production rate inferred from a 107 ka-old lava flow in northern Chile; 22°S-3400 m a.s.l

Available geological calibration sites used to estimate the rate at which cosmogenic 3He is produced at the Earth’s surface are mostly clustered in medium to high latitudes. Moreover, most of them have exposure histories shorter than tens of thousands of years. This lack of sites prevents a qualitative assessment of available production models used to convert cosmogenic 3He concentrations into exposure ages and/or denudation rates. It thus limits our ability to take into account the atmospheric, geomagnetic and solar modulation conditions that might have affected the production of cosmogenic nuclides in the past for longer exposure histories and in low latitude regions. We present the cosmogenic 3He production rate inferred from a new geological calibration site located in northern Chile. Five samples were collected on the surface of the largest and best-preserved lava flow of the San Pedro volcano (21.934°S-68.510°W- 3390 m a.s.l), which displays pristine crease-structure features. 40Ar/39Ar dating yield a reliable plateau age of 107±12 ka for the eruption of this lava flow. Eight pyroxene aliquots separated from the surface samples yield a weighted average cosmogenic 3He concentration of 99.3±1.2 Mat.g-1 from which a local cosmogenic 3He production rate of 928±101 at.g-1.yr-1 is calculated. The local production rate is then scaled to a sea level high latitude (SLHL) reference position using different combinations of geographic spatialization schemes, atmosphere models and geomagnetic field reconstructions, yielding SLHL production rates between 103±11 and 130±14 at.g-1.yr-1 consistent with the most recent estimates available from the literature. Finally, we use the same scaling frameworks to re-evaluate the mean global-scale cosmogenic 3He production rate in olivine and pyroxene minerals at 120±16 at.g-1.yr-1 from the compilation of previously published calibration datasets.