Sedimentary evolution and environmental history of Lake Van (Turkey) over the past 600,000 years

The lithostratigraphic framework of Lake Van, eastern Turkey, has been systematically analysed to document the sedimentary evolution and the environmental history of the lake during the past ca 600,000 years. The lithostratigraphy and chemostratigraphy of a 219 m long drill core from Lake Van serves to separate global climate oscillations from local factors caused by tectonic and volcanic activity. An age model was established based on the climatostratigraphic alignment of chemical and lithological signatures, validated by 40Ar/39Ar ages. The drilled sequence consists of ca 76% lacustrine carbonaceous clayey silt, ca 2% fluvial deposits, ca 17% volcaniclastic deposits and 5% gaps. Six lacustrine lithotypes were separated from the fluvial and event deposits, such as volcaniclastics (ca 300 layers) and graded beds (ca 375 layers), and their depositional environments are documented. These lithotypes are: (i) graded beds frequently intercalated with varved clayey silts reflect rising lake-levels during the terminations; (ii) varved clayey silts reflect strong seasonality and an intralake oxic–anoxic boundary, for example, lake-level highstands during interglacials/interstadials; (iii) CaCO3-rich banded sediments are representative of a lowering of the oxic-anoxic boundary, for example, lake-level decreases during glacial inceptions; (iv) CaCO3-poor banded and mottled clayey silts reflect an oxic–anoxic boundary close to the sediment-water interface, for example, lake-level low-stands during glacials/stadials; (v) diatomaceous muds were deposited during the early beginning of the lake as a fresh water system; and (vi) fluvial sands and gravels indicate the initial flooding of the lake basin. The recurrence of lithologies (i) to (iv) follows the past five glacial/interglacial cycles. A 20 m thick disturbed unit reflects an interval of major tectonic activity in Lake Van at ca 414 ka BP.

A reconstruction of radiocarbon production and total solar irradiance from the Holocene ¹⁴C and CO₂ records: implications of data and model uncertainties

Radiocarbon production, solar activity, total solar irradiance (TSI) and solar-induced climate change are reconstructed for the Holocene (10 to 0 kyr BP), and TSI is predicted for the next centuries. The IntCal09/SHCal04 radiocarbon and ice core CO2 records, reconstructions of the geomagnetic dipole, and instrumental data of solar activity are applied in the Bern3D-LPJ, a fully featured Earth system model of intermediate complexity including a 3-D dynamic ocean, ocean sediments, and a dynamic vegetation model, and in formulations linking radiocarbon production, the solar modulation potential, and TSI. Uncertainties are assessed using Monte Carlo simulations and bounding scenarios. Transient climate simulations span the past 21 thousand years, thereby considering the time lags and uncertainties associated with the last glacial termination. Our carbon-cycle-based modern estimate of radiocarbon production of 1.7 atoms cm−2 s−1 is lower than previously reported for the cosmogenic nuclide production model by Masarik and Beer (2009) and is more in-line with Kovaltsov et al. (2012). In contrast to earlier studies, periods of high solar activity were quite common not only in recent millennia, but throughout the Holocene. Notable deviations compared to earlier reconstructions are also found on decadal to centennial timescales. We show that earlier Holocene reconstructions, not accounting for the interhemispheric gradients in radiocarbon, are biased low. Solar activity is during 28% of the time higher than the modern average (650 MeV), but the absolute values remain weakly constrained due to uncertainties in the normalisation of the solar modulation to instrumental data. A recently published solar activity–TSI relationship yields small changes in Holocene TSI of the order of 1 W m−2 with a Maunder Minimum irradiance reduction of 0.85 ± 0.16 W m−2. Related solar-induced variations in global mean surface air temperature are simulated to be within 0.1 K. Autoregressive modelling suggests a declining trend of solar activity in the 21st century towards average Holocene conditions.

Independent variations of CH₄ emissions and isotopic composition over the past 160,000 years

During the last glacial cycle, greenhouse gas concentrations fluctuated on decadal and longer timescales. Concentrations of methane, as measured in polar ice cores, show a close connection with Northern Hemisphere temperature variability, but the contribution of the various methane sources and sinks to changes in concentration is still a matter of debate. Here we assess changes in methane cycling over the past 160,000 years by measurements of the carbon isotopic composition delta C-13 of methane in Antarctic ice cores from Dronning Maud Land and Vostok. We find that variations in the delta C-13 of methane are not generally correlated with changes in atmospheric methane concentration, but instead more closely correlated to atmospheric CO2 concentrations. We interpret this to reflect a climatic and CO2-related control on the isotopic signature of methane source material, such as ecosystem shifts in the seasonally inundated tropical wetlands that produce methane. In contrast, relatively stable delta C-13 values occurred during intervals of large changes in the atmospheric loading of methane. We suggest that most methane sources-most notably tropical wetlands-must have responded simultaneously to climate changes across these periods.

Dating groundwater in the Bohemian Cretaceous Basin: Understanding tracer variations in the subsurface

The northern section of the Bohemian Cretaceous Basin has been the site of intensive U exploitation with harmful impacts on groundwater quality. The understanding of groundwater flow and age distribution is crucial for the prediction of the future dispersion and impact of the contamination. State of the art tracer methods (3H, 3He, 4He, 85Kr, 39Ar and 14C) were, therefore, used to obtain insights to ageing and mixing processes of groundwater along a north–south flow line in the centre of the two most important aquifers of Cenomanian and middle Turonian age. Dating of groundwater is particularly complex in this area as: (i) groundwater in the Cenomanian aquifer is locally affected by fluxes of geogenic and biogenic gases (e.g. CO2, CH4, He) and by fossil brines in basement rocks rich in Cl and SO4; (ii) a thick unsaturated zone overlays the Turonian aquifer; (iii) a periglacial climate and permafrost conditions prevailed during the Last Glacial Maximum (LGM), and iv) the wells are mostly screened over large depth intervals. Large disagreements in 85Kr and 3H/3He ages indicate that processes other than ageing have affected the tracer data in the Turonian aquifer. Mixing with older waters (>50 a) was confirmed by 39Ar activities. An inverse modelling approach, which included time lags for tracer transport throughout the unsaturated zone and degassing of 3He, was used to estimate the age of groundwater. Best fits between model and field results were obtained for mean residence times varying from modern up to a few hundred years. The presence of modern water in this aquifer is correlated with the occurrence of elevated pollution (e.g. nitrates). An increase of reactive geochemical indicators (e.g. Na) and radiogenic 4He, and a decrease in 14C along the flow direction confirmed groundwater ageing in the deeper confined Cenomanian aquifer. Radiocarbon ages varied from a few hundred years to more than 20 ka. Initial 14C activity for radiocarbon dating was calibrated by means of 39Ar measurements. The 14C age of a sample recharged during the LGM was further confirmed by depleted stable isotope signatures and near freezing point noble gas temperature. Radiogenic 4He accumulated in groundwater with concentrations increasing linearly with 14C ages. This enabled the use of 4He to validate the dating range of 14C and extend it to other parts of this aquifer. In the proximity of faults, 39Ar in excess of modern concentrations and 14C dead CO2 sources, elevated 3He/4He ratios and volcanic activity in Oligocene to Quaternary demonstrate the influence of gas of deeper origin and impeded the application of 4He, 39Ar and 14C for groundwater dating.

Greenland accumulation and its connection to the large-scale atmospheric circulation in ERA-Interim and paleoclimate simulations

Changes in Greenland accumulation and the stability in the relationship between accumulation variability and large-scale circulation are assessed by performing time-slice simulations for the present day, the preindustrial era, the early Holocene, and the Last Glacial Maximum (LGM) with a comprehensive climate model. The stability issue is an important prerequisite for reconstructions of Northern Hemisphere atmospheric circulation variability based on accumulation or precipitation proxy records from Greenland ice cores. The analysis reveals that the relationship between accumulation variability and large-scale circulation undergoes a significant seasonal cycle. As the contributions of the individual seasons to the annual signal change, annual mean accumulation variability is not necessarily related to the same atmospheric circulation patterns during the different climate states. Interestingly, within a season, local Greenland accumulation variability is indeed linked to a consistent circulation pattern, which is observed for all studied climate periods, even for the LGM. Hence, it would be possible to deduce a reliable reconstruction of seasonal atmospheric variability (e.g., for North Atlantic winters) if an accumulation or precipitation proxy were available that resolves single seasons. We further show that the simulated impacts of orbital forcing and changes in the ice sheet topography on Greenland accumulation exhibit strong spatial differences, emphasizing that accumulation records from different ice core sites regarding both interannual and long-term (centennial to millennial) variability cannot be expected to look alike since they include a distinct local signature. The only uniform signal to external forcing is the strong decrease in Greenland accumulation during glacial (LGM) conditions and an increase associated with the recent rise in greenhouse gas concentrations.

The effect of abrupt climatic warming on biogeochemical cycling and N2O emissions in a terrestrial ecosystem

The large, rapid increase in atmospheric N2O concentrations that occurred concurrent with the abrupt warming at the end of the Last Glacial period might have been the result of a reorganization in global biogeochemical cycles. To explore the sensitivity of nitrogen cycling in terrestrial ecosystems to abrupt warming, we combined a scenario of climate and vegetation composition change based on multiproxy data for the Oldest Dryas–Bølling abrupt warming event at Gerzensee, Switzerland, with a biogeochemical model that simulates terrestrial N uptake and release, including N2O emissions. As for many central European sites, the pollen record at the Gerzensee is remarkable for the abundant presence of the symbiotic nitrogen fixer Hippophaë rhamnoides (L.) during the abrupt warming that also marks the beginning of primary succession on immature glacial soils. Here we show that without additional nitrogen fixation, climate change results in a significant increase of N2O emissions of approximately factor 3.4 (from 6.4 ± 1.9 to 21.6 ± 5.9 mg N2O–N m− 2 yr− 1). Each additional 1000 mg m− 2 yr− 1 of nitrogen added to the ecosystem through N-fixation results in additional N2O emissions of 1.6 mg N2O–N m− 2 yr− 1 for the time with maximum H. rhamnoides coverage. Our results suggest that local reactions of emissions to abrupt climate change could have been considerably faster than the overall atmospheric concentration changes observed in polar ice. Nitrogen enrichment of soils due to the presence of symbiotic N-fixers during early primary succession not only facilitates the establishment of vegetation on soils in their initial stage of development, but can also have considerable influence on biogeochemical cycles and the release of reactive nitrogen trace gases to the atmosphere.

Temperature reconstruction from 10 to 120 kyr b2k from the NGRIP ice core

In order to reconstruct the temperature of the North Greenland Ice Core Project (NGRIP) site, new measurements of δ15N have been performed covering the time period from the beginning of the Holocene to Dansgaard–Oeschger (DO) event 8. Together with previously measured and mostly published δ15N data, we present for the first time a NGRIP temperature reconstruction for the whole last glacial period from 10 to 120 kyr b2k (thousand years before 2000 AD) including every DO event based on δ15N isotope measurements combined with a firn densification and heat diffusion model. The detected temperature rises at the onset of DO events range from 5 °C (DO 25) up to 16.5 °C (DO 11) with an uncertainty of ±3 °C. To bring measured and modelled data into agreement, we had to reduce the accumulation rate given by the NGRIP ss09sea06bm timescale in some periods by 30 to 35%, especially during the last glacial maximum. A comparison between reconstructed temperature and δ18Oice data confirms that the isotopic composition of the stadial was strongly influenced by seasonality. We evidence an anticorrelation between the variations of the δ18Oice sensitivity to temperature (referred to as α) and obliquity in agreement with a simple Rayleigh distillation model. Finally, we suggest that α might be influenced by the Northern Hemisphere ice sheet volume.

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.

Estimation and calibration of the water isotope differential diffusion length in ice core records

Palaeoclimatic information can be retrieved from the diffusion of the stable water isotope signal during firnification of snow. The diffusion length, a measure for the amount of diffusion a layer has experienced, depends on the firn temperature and the accumulation rate. We show that the estimation of the diffusion length using power spectral densities (PSDs) of the record of a single isotope species can be biased by uncertainties in spectral properties of the isotope signal prior to diffusion. By using a second water isotope and calculating the difference in diffusion lengths between the two isotopes, this problem is circumvented. We study the PSD method applied to two isotopes in detail and additionally present a new forward diffusion method for retrieving the differential diffusion length based on the Pearson correlation between the two isotope signals. The two methods are discussed and extensively tested on synthetic data which are generated in a Monte Carlo manner. We show that calibration of the PSD method with this synthetic data is necessary to be able to objectively determine the differential diffusion length. The correlation-based method proves to be a good alternative for the PSD method as it yields precision equal to or somewhat higher than the PSD method. The use of synthetic data also allows us to estimate the accuracy and precision of the two methods and to choose the best sampling strategy to obtain past temperatures with the required precision. In addition to application to synthetic data the two methods are tested on stable-isotope records from the EPICA (European Project for Ice Coring in Antarctica) ice core drilled in Dronning Maud Land, Antarctica, showing that reliable firn temperatures can be reconstructed with a typical uncertainty of 1.5 and 2 °C for the Holocene period and 2 and 2.5 °C for the last glacial period for the correlation and PSD method, respectively.

Late Quaternary history of North Eurasian Norway spruce ( Picea abies ) and Siberian spruce ( Picea obovata ) inferred from macrofossils, pollen and cytoplasmic DNA variation

Aim We used combined palaeobotanical and genetic data to assess whether Norway spruce (Picea abies) and Siberian spruce (Picea obovata), two major components of the Eurasian boreal forests, occupied separate glacial refugia, and to test previous hypotheses on their distinction, geographical delimitation and introgression. Location The range of Norway spruce in northern Europe and Siberian spruce in northern Asia. Methods Pollen data and recently compiled macrofossil records were summarized for the Last Glacial Maximum (LGM), late glacial and Holocene. Genetic variation was assessed in 50 populations using one maternally (mitochondrial nad1) and one paternally (chloroplast trnT–trnL) inherited marker and analysed using spatial analyses of molecular variance (SAMOVA). Results Macrofossils showed that spruce was present in both northern Europe and Siberia at the LGM. Congruent macrofossil and pollen data from the late glacial suggested widespread expansions of spruce in the East European Plain, West Siberian Plain, southern Siberian mountains and the Baikal region. Colonization was largely completed during the early Holocene, except in the formerly glaciated area of northern Europe. Both DNA markers distinguished two highly differentiated groups that correspond to Norway spruce and Siberian spruce and coincide spatially with separate LGM spruce occurrences. The division of the mtDNA variation was geographically well defined and occurred to the east of the Ural Mountains along the Ob River, whereas the cpDNA variation showed widespread admixture. Genetic diversity of both DNA markers was higher in western than in eastern populations. Main conclusions North Eurasian Norway spruce and Siberian spruce are genetically distinct and occupied separate LGM refugia, Norway spruce on the East European Plain and Siberian spruce in southern Siberia, where they were already widespread during the late glacial. They came into contact in the basin of the Ob River and probably hybridized. The lower genetic diversity in the eastern populations may indicate that Siberian spruce suffered more from past climatic fluctuations than Norway spruce.

Ignoring heterozygous sites biases phylogenomic estimates of divergence times: implications for the evolutionary history of Microtus voles

Phylogenetic reconstruction of the evolutionary history of closely related organisms may be difficult because of the presence of unsorted lineages and of a relatively high proportion of heterozygous sites that are usually not handled well by phylogenetic programs. Genomic data may provide enough fixed polymorphisms to resolve phylogenetic trees, but the diploid nature of sequence data remains analytically challenging. Here, we performed a phylogenomic reconstruction of the evolutionary history of the common vole (Microtus arvalis) with a focus on the influence of heterozygosity on the estimation of intraspecific divergence times. We used genome-wide sequence information from 15 voles distributed across the European range. We provide a novel approach to integrate heterozygous information in existing phylogenetic programs by repeated random haplotype sampling from sequences with multiple unphased heterozygous sites. We evaluated the impact of the use of full, partial, or no heterozygous information for tree reconstructions on divergence time estimates. All results consistently showed four deep and strongly supported evolutionary lineages in the vole data. These lineages undergoing divergence processes split only at the end or after the last glacial maximum based on calibration with radiocarbon-dated paleontological material. However, the incorporation of information from heterozygous sites had a significant impact on absolute and relative branch length estimations. Ignoring heterozygous information led to an overestimation of divergence times between the evolutionary lineages of M. arvalis. We conclude that the exclusion of heterozygous sites from evolutionary analyses may cause biased and misleading divergence time estimates in closely related taxa.

Depositional modes and lake-level variability at Lake Towuti, Indonesia, during the past ~29 kyr BP

Lake Towuti (2.5°S, 121.5°E) is a long-lived, tectonic lake located on the Island of Sulawesi, Indonesia, and in the center of the Indo-Pacific warm pool (IPWP). Lake Towuti is connected with upstream lakes Matano and Mahalona through the Mahalona River, which constitutes the largest inlet to the lake. The Mahalona River Delta is prograding into Lake Towuti’s deep northern basin thus exerting significant control on depositional processes in the basin. We combine high-resolution seismic reflection and sedimentological datasets from a 19.8-m-long sediment piston core from the distal edge of this delta to characterize fluctuations in deltaic sedimentation during the past ~29 kyr BP and their relation to climatic change. Our datasets reveal that, in the present, sedimentation is strongly influenced by deposition of laterally transported sediments sourced from the Mahalona River Delta. Variations in the amount of laterally transported sediments, as expressed by coarse fraction amounts in pelagic muds and turbidite recurrence rates and cumulative thicknesses, are primarily a function of lake-level induced delta slope instability and delta progradation into the basin. We infer lowest lake-levels between ~29 and 16, a gradual lake level rise between ~16 and 11, and high lake-levels between ~11 and 0 kyr BP. Periods of highest turbidite deposition, ~26 to 24 and ~18 to 16 kyr BP coincide with Heinrich events 2 and 1, respectively. Our lake-level reconstruction therefore supports previous observations based on geochemical hydroclimate proxies of a very dry last glacial and a wet Holocene in the region, and provides new evidence of millennial-scale variations in moisture balance in the IPWP.

Untangling a Holocene pollen record with forest model simulations and independent climate data

Adaptation potential of forests to rapid climatic changes can be assessed from vegetation dynamics during past climatic changes as preserved in fossil pollen data. However, pollen data reflect the integrated effects of climate and biotic processes, such as establishment, survival, competition, and migration. To disentangle these processes, we compared an annually laminated late Würm and Holocene pollen record from the Central Swiss Plateau with simulations of a dynamic forest patch model. All input data used in the simulations were largely independent from pollen data; i.e. the presented analysis is non-circular. Temperature and precipitation scenarios were based on reconstructions from pollen-independent sources. The earliest arrival times of the species at the study site after the last glacial were inferred from pollen maps. We ran a series of simulations under different combinations of climate and immigration scenarios. In addition, the sensitivity of the simulated presence/absence of four major species to changes in the climate scenario was examined. The pattern of the pollen record could partly be explained by the used climate scenario, mostly by temperature. However, some features, in particular the absence of most species during the late Würm could only be simulated if the winter temperature anomalies of the used scenario were decreased considerably. Consequently, we had to assume in the simulations, that most species immigrated during or after the Younger Dryas (12 000 years BP), Abies and Fagus even later. Given the timing of tree species immigration, the vegetation was in equilibrium with climate during long periods, but responded with lags at the time-scale of centuries to millennia caused by a secondary succession after rapid climatic changes such as at the end of Younger Dryas, or immigration of dominant taxa. Climate influenced the tree taxa both directly and indirectly by changing inter-specific competition. We concluded, that also during the present fast climatic change, species migration might be an important process, particularly if geographic barriers, such as the Alps are in the migrational path.

Long distance dispersal shaped patterns of human genetic diversity in Eurasia

Most previous attempts at reconstructing the past history of human populations did not explicitly take geography into account, or considered very simple scenarios of migration and ignored environmental information. However, it is likely that the Last Glacial Maximum (LGM) affected the demography and the range of many species, including our own. Moreover, long-distance dispersal (LDD) may have been an important component of human migrations, allowing fast colonization of new territories and preserving high levels of genetic diversity. Here, we use a high-quality microsatellite dataset genotyped in 22 populations to estimate the posterior probabilities of several scenarios for the settlement of the Old World by modern humans. We considered models ranging from a simple spatial expansion to others including LDD and a LGM-induced range contraction, as well as Neolithic demographic expansions. We find that scenarios with LDD are much better supported by data than models without LDD. Nevertheless, we show evidence that LDD events to empty habitats were strongly prevented during the settlement of Eurasia. This unexpected absence of LDD ahead of the colonization wave front could have been caused by an Allee effect, either due to intrinsic causes such as an inbreeding depression built during the expansion, or to extrinsic causes such as direct competition with archaic humans. Overall, our results suggest only a relatively limited effect of the LGM-contraction on current patterns of human diversity. This is in clear contrast with the major role of LDD migrations, which have potentially contributed to the intermingled genetic structure of Eurasian populations.

A 22,000 14C year BP sediment and pollen record of climate change from Laguna Miscanti (23°S), northern Chile

Lake sediments and pollen, spores and algae from the high-elevation endorheic Laguna Miscanti (22°45′S, 67°45′W, 4140 m a.s.l., 13.5 km2 water surface, 10 m deep) in the Atacama Desert of northern Chile provide information about abrupt and high amplitude changes in effective moisture. Although the lack of terrestrial organic macrofossils and the presence of a significant 14C reservoir effect make radiocarbon dating of lake sediments very difficult, we propose the following palaeoenvironmental history. An initial shallow freshwater lake (ca. 22,000 14C years BP) disappeared during the extremely dry conditions of the Last Glacial Maximum (LGM; 18,000 14C years BP). That section is devoid of pollen. The late-glacial lake transgression started around 12,000 14C years BP, peaked in two phases between ca. 11,000 and <9000 14C years BP, and terminated around 8000 14C years BP. Effective moisture increased more than three times compared to modern conditions (∼200 mm precipitation), and a relatively dense terrestrial vegetation was established. Very shallow hypersaline lacustrine conditions prevailed during the mid-Holocene until ca. 3600 14C years BP. However, numerous drying and wetting cycles suggest frequent changes in moisture, maybe even individual storms during the mid-Holocene. After several humid spells, modern conditions were reached at ca. 3000 14C years BP. Comparison between limnogeological data and pollen of terrestrial plants suggest century-scale response lags. Relatively constant concentrations of long-distance transported pollen from lowlands east of the Andes suggest similar atmospheric circulation patterns (mainly tropical summer rainfall) throughout the entire period of time. These findings compare favorably with other regional paleoenvironmental data.