Age modeling of young non-varved lake sediments: methods and limits. Examples from two lakes in Central Chile

High-resolution and highly precise age models for recent lake sediments (last 100–150 years) are essential for quantitative paleoclimate research. These are particularly important for sedimentological and geochemical proxies, where transfer functions cannot be established and calibration must be based upon the relation of sedimentary records to instrumental data. High-precision dating for the calibration period is most critical as it determines directly the quality of the calibration statistics. Here, as an example, we compare radionuclide age models obtained on two high-elevation glacial lakes in the Central Chilean Andes (Laguna Negra: 33°38′S/70°08′W, 2,680 m a.s.l. and Laguna El Ocho: 34°02′S/70°19′W, 3,250 m a.s.l.). We show the different numerical models that produce accurate age-depth chronologies based on 210Pb profiles, and we explain how to obtain reduced age-error bars at the bottom part of the profiles, i.e., typically around the end of the 19th century. In order to constrain the age models, we propose a method with five steps: (i) sampling at irregularly-spaced intervals for 226Ra, 210Pb and 137Cs depending on the stratigraphy and microfacies, (ii) a systematic comparison of numerical models for the calculation of 210Pb-based age models: constant flux constant sedimentation (CFCS), constant initial concentration (CIC), constant rate of supply (CRS) and sediment isotope tomography (SIT), (iii) numerical constraining of the CRS and SIT models with the 137Cs chronomarker of AD 1964 and, (iv) step-wise cross-validation with independent diagnostic environmental stratigraphic markers of known age (e.g., volcanic ash layer, historical flood and earthquakes). In both examples, we also use airborne pollutants such as spheroidal carbonaceous particles (reflecting the history of fossil fuel emissions), excess atmospheric Cu deposition (reflecting the production history of a large local Cu mine), and turbidites related to historical earthquakes. Our results show that the SIT model constrained with the 137Cs AD 1964 peak performs best over the entire chronological profile (last 100–150 years) and yields the smallest standard deviations for the sediment ages. Such precision is critical for the calibration statistics, and ultimately, for the quality of the quantitative paleoclimate reconstruction. The systematic comparison of CRS and SIT models also helps to validate the robustness of the chronologies in different sections of the profile. Although surprisingly poorly known and under-explored in paleolimnological research, the SIT model has a great potential in paleoclimatological reconstructions based on lake sediments

Late Holocene summer temperatures in the central Andes reconstructed from the sediments of high-elevatin Laguna Chepical, Chile (32° S)

High-resolution reconstructions of climate variability that cover the past millennia are necessary to improve the understanding of natural and anthropogenic climate change across the globe. Although numerous records are available for the mid- and high-latitudes of the Northern Hemisphere, global assessments are still compromised by the scarcity of data from the Southern Hemisphere. This is particularly the case for the tropical and subtropical areas. In addition, high elevation sites in the South American Andes may provide insight into the vertical structure of climate change in the mid-troposphere. This study presents a 3000 yr-long austral summer (November to February) temperature reconstruction derived from the 210Pb- and 14C-dated organic sediments of Laguna Chepical (32°16' S, 70°30' W, 3050 m a.s.l.), a high-elevation glacial lake in the subtropical Andes of central Chile. Scanning reflectance spectroscopy in the visible light range provided the spectral index R570/R630, which reflects the clay mineral content in lake sediments. For the calibration period (AD 1901–2006), the R570/R630 data were regressed against monthly meteorological reanalysis data, showing that this proxy was strongly and significantly correlated with mean summer (NDJF) temperatures (R3 yr = −0.63, padj = 0.01). This calibration model was used to make a quantitative temperature reconstruction back to 1000 BC. The reconstruction (with a model error RMSEPboot of 0.33 °C) shows that the warmest decades of the past 3000 yr occurred during the calibration period. The 19th century (end of the Little Ice Age (LIA)) was cool. The prominent warmth reconstructed for the 18th century, which was also observed in other records from this area, seems systematic for subtropical and southern South America but remains difficult to explain. Except for this warm period, the LIA was generally characterized by cool summers. Back to AD 1400, the results from this study compare remarkably well to low altitude records from the Chilean Central Valley and southern South America. However, the reconstruction from Laguna Chepical does not show a warm Medieval Climate Anomaly during the 12–13th century, which is consistent with records from tropical South America. The Chepical record also indicates substantial cooling prior to 800 BC. This coincides with well-known regional as well as global glacier advances which have been attributed to a grand solar minimum. This study thus provides insight into the climatic drivers and temperature patterns in a region for which currently very few data are available. It also shows that since ca. AD 1400, long-term temperature patterns were generally similar at low and high altitudes in central Chile.

Quantitative high-resolution warm season rainfall recorded in varved sediments of Lake Oeschinen, northern Swiss Alps: calibration and validation AD 1901–2008

High-resolution, well-calibrated records of lake sediments are critically important for quantitative climate reconstructions, but they remain a methodological and analytical challenge. While several comprehensive paleotemperature reconstructions have been developed across Europe, only a few quantitative high-resolution studies exist for precipitation. Here we present a calibration and verification study of lithoclastic sediment proxies from proglacial Lake Oeschinen (46°30′N, 7°44′E, 1,580 m a.s.l., north–west Swiss Alps) that are sensitive to rainfall for the period AD 1901–2008. We collected two sediment cores, one in 2007 and another in 2011. The sediments are characterized by two facies: (A) mm-laminated clastic varves and (B) turbidites. The annual character of the laminae couplets was confirmed by radiometric dating (210Pb, 137Cs) and independent flood-layer chronomarkers. Individual varves consist of a dark sand-size spring-summer layer enriched in siliciclastic minerals and a lighter clay-size calcite-rich winter layer. Three subtypes of varves are distinguished: Type I with a 1–1.5 mm fining upward sequence; Type II with a distinct fine-sand base up to 3 mm thick; and Type III containing multiple internal microlaminae caused by individual summer rainstorm deposits. Delta-fan surface samples and sediment trap data fingerprint different sediment source areas and transport processes from the watershed and confirm the instant response of sediment flux to rainfall and erosion. Based on a highly accurate, precise and reproducible chronology, we demonstrate that sediment accumulation (varve thickness) is a quantitative predictor for cumulative boreal alpine spring (May–June) and spring/summer (May–August) rainfall (rMJ = 0.71, rMJJA = 0.60, p < 0.01). Bootstrap-based verification of the calibration model reveals a root mean squared error of prediction (RMSEPMJ = 32.7 mm, RMSEPMJJA = 57.8 mm) which is on the order of 10–13 % of mean MJ and MJJA cumulative precipitation, respectively. These results highlight the potential of the Lake Oeschinen sediments for high-resolution reconstructions of past rainfall conditions in the northern Swiss Alps, central and eastern France and south-west Germany.

Chronology of Lake El'gygytgyn sediments – a combined magnetostratigraphic, palaeoclimatic and orbital tuning study based on multi-parameter analyses

A 318-metre-long sedimentary profile drilled by the International Continental Scientific Drilling Program (ICDP) at Site 5011-1 in Lake El’gygytgyn, Far East Russian Arctic, has been analysed for its sedimentologic response to global climate modes by chronostratigraphic methods. The 12 km wide lake is sited off-centre in an 18 km large crater that was created by the impact of a meteorite 3.58 Ma ago. Since then sediments have been continuously deposited. For establishing their chronology, major reversals of the earth’s magnetic field provided initial tie points for the age model, confirming that the impact occurred in the earliest geomagnetic Gauss chron. Various stratigraphic parameters, reflecting redox conditions at the lake floor and climatic conditions in the catchment were tuned synchronously to Northern Hemisphere insolation variations and the marine oxygen isotope stack, respectively. Thus, a robust age model comprising more than 600 tie points could be defined. It could be shown that deposition of sediments in Lake El’gygytgyn occurred in concert with global climatic cycles. The upper �160m of sediments represent the past 3.3 Ma, equivalent to sedimentation rates of 4 to 5 cm ka−1, whereas the lower 160m represent just the first 0.3 Ma after the impact, equivalent to sedimentation rates in the order of 45 cm ka−1. This study also provides orbitally tuned ages for a total of 8 tephras deposited in Lake El’gygytgyn.

Environmental change within the Balkan region during the past ca. 50 ka recorded in the sediments from lakes Prespa and Ohrid

Abstract. Lakes Prespa and Ohrid, in the Balkan region, are considered to be amongst the oldest lakes in Europe. Both lakes are hydraulically connected via karst aquifers. From Lake Ohrid, several sediment cores up to 15m long have been studied over the last few years. Here, we document the first long sediment record from nearby Lake Prespa to clarify the influence of Lake Prespa on Lake Ohrid and the environmental history of the region. Radiocarbon dating and dated tephra layers provide robust age control and indicate that the 10.5m long sediment record from Lake Prespa reaches back to 48 ka. Glacial sedimentation is characterized by low organic matter content and absence of carbonates in the sediments, which indicate oligotrophic conditions in both lakes. Holocene sedimentation is characterized by particularly high carbonate content in Lake Ohrid and by particularly high organic matter content in Lake Prespa, which indicates a shift towards more mesotrophic conditions in the latter. Long-term environmental change and short-term events, such as related to the Heinrich events during the Pleistocene or the 8.2 ka cooling event during the Holocene, are well recorded in both lakes, but are only evident in certain proxies. The comparison of the sediment cores from both lakes indicates that environmental change affects particularly the trophic state of Lake Prespa due to its lower volume and water depth.

Fourier transform infrared spectroscopy, a new method for rapid determination of total organic and inorganic carbon and biogenic silica concentration in lake sediments

Abstract We demonstrate the use of Fourier transform infrared spectroscopy (FTIRS) to make quantitative measures of total organic carbon (TOC), total inorganic carbon (TIC) and biogenic silica (BSi) concentrations in sediment. FTIRS is a fast and costeffective technique and only small sediment samples are needed (0.01 g). Statistically significant models were developed using sediment samples from northern Sweden and were applied to sediment records from Sweden, northeast Siberia and Macedonia. The correlation between FTIRS-inferred values and amounts of biogeochemical constituents assessed conventionally varied between r = 0.84–0.99 for TOC, r = 0.85– 0.99 for TIC, and r = 0.68–0.94 for BSi. Because FTIR spectra contain information on a large number of both inorganic and organic components, there is great potential for FTIRS to become an important tool in paleolimnology.

Lipid biomarkers in Holocene and glacial sediments from ancient Lake Ohrid (Macedonia, Albania)

Abstract. Organic matter preserved in Lake Ohrid sediments originates from aquatic and terrestrial sources. Its variable composition reflects climate-controlled changes in the lake basin’s hydrology and related organic matter export, i.e. changes in primary productivity, terrestrial plant matter input and soil erosion. Here, we present first results from lipid biomarker investigations of Lake Ohrid sediments from two near-shore settings: site Lz1120 near the southern shore, with low-lying lands nearby and probably influenced by river discharge, and site Co1202 which is close to the steep eastern slopes. Variable proportions of terrestrial n-alkanoic acids and n-alkanols as well as compositional changes of !- hydroxy acids document differences in soil organic matter supply between the sites and during different climate stages (glacial, Holocene, 8.2 ka cooling event). Changes in the vegetation cover are suggested by changes in the dominant chain length of terrestrial n-alkanols. Effective microbial degradation of labile organic matter and in situ contribution of organic matter derived from the microbes themselves are both evident in the sediments. We found evidence for anoxic conditions within the photic zone by detecting epicholestanol and tetrahymanol from sulphur-oxidising phototrophic bacteria and bacterivorous ciliates and for the influence of a settled human community from the occurrence of coprostanol, a biomarker for human and animal faeces (pigs, sheep, goats), in an early Holocene sample. This study illustrates the potential of lipid biomarkers for future environmental reconstructions using one of Europe’s oldest continental climate archives, Lake Ohrid.

Small aliquot and single grain IRSL and post-IR IRSL dating of fluvial and alluvial sediments from the Pativilca valley, Peru

Infrared stimulated luminescence (IRSL) and post-IR IRSL are applied to small aliquots and single grains to determine the equivalent dose (De) of eleven alluvial and fluvial sediment samples collected in the Pativilca valley, Central Peru at ca. 10°S latitude. Small aliquot De distributions are rather symmetric and display over-dispersion values between 15 and 46%. Small aliquot g-values range between 4 and 8% per decade for the IRSL and 1 and 2% per decade for the post-IR IRSL signal. The single grain De distributions are highly over-dispersed with some of them skewed to higher doses, implying partial bleaching; this is especially true for the post-IR IRSL. Measurements of a modern analog reveal that residuals due to partial bleaching are present in both the IRSL as well as the post-IR IRSL signal. The g-values of individual grains exhibit a wide range with high individual uncertainties and might contribute significantly to the spread of the single grain De values, at least for the IRSL data. Electron Microprobe Analysis performed on single grains reveal that a varying K-content can be excluded as the origin of over-dispersion. Final ages for the different approaches are calculated using the Central Age Model and the Minimum Age Model (MAM). The samples are grouped into well-beached, potentially well-bleached and partially bleached according to the evaluation of the single grain distributions and the agreement of age estimates between methods. The application of the MAM to the single grain data resulted in consistent age estimates for both the fading corrected IRSL and the post-IR IRSL ages, and suggests that both approaches are suitable for dating these samples. Keywords

Lake sediments as natural seismographs: A compiled record of Late Quaternary earthquakes in Central Switzerland and its implication for Alpine deformation

Central Switzerland lies tectonically in an intraplate area and recurrence rates of strong earthquakes exceed the time span covered by historic chronicles. However, many lakes are present in the area that act as natural seismographs: their continuous, datable and high-resolution sediment succession allows extension of the earthquake catalogue to pre-historic times. This study reviews and compiles available data sets and results from more than 10 years of lacustrine palaeoseismological research in lakes of northern and Central Switzerland. The concept of using lacustrine mass-movement event stratigraphy to identify palaeo-earthquakes is showcased by presenting new data and results from Lake Zurich. The Late Glacial to Holocene mass-movement units in this lake document a complex history of varying tectonic and environmental impacts. Results include sedimentary evidence of three major and three minor, simultaneously triggered basin-wide lateral slope failure events interpreted as the fingerprints of palaeoseismic activity. A refined earthquake catalogue, which includes results from previous lake studies, reveals a non-uniform temporal distribution of earthquakes in northern and Central Switzerland. A higher frequency of earthquakes in the Late Glacial and Late Holocene period documents two different phases of neotectonic activity; they are interpreted to be related to isostatic post-glacial rebound and relatively recent (re-)activation of seismogenic zones, respectively. Magnitudes and epicentre reconstructions for the largest identified earthquakes provide evidence for two possible earthquake sources: (i) a source area in the region of the Alpine or Sub-Alpine Front due to release of accumulated north-west/south-east compressional stress related to an active basal thrust beneath the Aar massif; and (ii) a source area beneath the Alpine foreland due to reactivation of deep-seated strike-slip faults. Such activity has been repeatedly observed instrumentally, for example, during the most recent magnitude 4.2 and 3.5 earthquakes of February 2012, near Zug. The combined lacustrine record from northern and Central Switzerland indicates that at least one of these potential sources has been capable of producing magnitude 6.2 to 6.7 events in the past.

Lake sediments as archives of recurrence rates and intensities of past flood events

Palaeoflood hydrology is an expanding field as the damage potential of flood and flood-related processes are increasing with the population density and the value of the infrastructure. Assessing the risk of these hazards in mountainous terrain requires knowledge about the frequency and severness of such events in the past. A wide range of methods is employed using diverse biologic, geomorphic or geologic evidences to track past flood events. Impact of floods are studied and dated on alluvial fans and cones using for example the growth disturbance of trees (Stoffel and Bollschweiler 2008; Schneuwly-Bollschweiler and Stoffel 2012: this volume) or stratigraphic layers deposited by debris flows, allowing to reconstruct past flood frequencies (Bardou et~al. 2003). Further downstream, the classical approach of palaeoflood hydrology (Kochel and Baker 1982) utilizes geomorphic indicators such as overbank sediments, silt lines and erosion features of floods along a river (e.g. Benito and Thorndycraft 2005). Fine-grained sediment settles out of the river suspension in eddies or backwater areas, where the flow velocity of the river is reduced. Records of these deposits at different elevations across a river’s profile can be used to assess the discharge of the past floods. This approach of palaeoflood hydrology studies was successfully applied in several river catchments (e.g. Ely et al. 1993; Macklin and Lewin 2003; O’Connor et al. 1994; Sheffer et al. 2003; Thorndycraft et al. 2005; Thorndycraft and Benito 2006). All these different reconstruction methods have their own advantages and disadvantages, but often these studies have a limited time coverage and the records are potentially incomplete due to lateral limits of depositional areas and due to the erosional power of fluvial processes that remove previously deposited flood witnesses. Here, we present a method that follows the sediment particle transported by a flood event to its final sink: the lacustrine basin.

Atypical Mo isotope signatures in eastern Mediterranean sediments

Suboxic sapropel S1 sediments of the eastern Mediterranean carry an untypically light δ98/95Mo signal, being even lighter than oxic Mn-enriched sediments and Mn-crusts/nodules, which were previously thought to carry the lightest Mo isotope values. The evaluation and comparison of oxic S1, diagenetically overprinted oxic S1, and suboxic S1 sediments indicates that this light isotope signal is primarily an oxic signal that has been overprinted by secondary diagenetic processes occurring after sediment burial. Such secondary processes bear the potential of additional Mo isotope fractionation in particular in non-steady state diagenetic environments that involve the discontinuous re-location of the redox boundary.

The history of subaquatic volcanism recorded in the sediments of Lake Kivu

Subaquatic volcanic activity has been ongoing in Lake Kivu since the early Holocene and has a dynamic effect on the biological productivity in the surface water, and the preservation of carbonate in the deep anoxic water. Groundwater discharge into the lake’s deepwater propels the upward advection of the water column that ultimately supplies nutrients to the surface water for biological production. The amount of nutrients supplied from the deepwater can be increased suddenly by (1) a cold meteorological event that drives deep seasonal mixing resulting in increased nutrients from below and oxygen from above, or (2) subaquatic volcanic activity that induces a buoyant hydrothermal plume, which entrains nutrients from the deepwater and results in anoxia or suboxic conditions in the surface water. Previous sedimentological studies in Lake Kivu have hypothesized that regional climatic changes are responsible for sudden changes in the preservation of carbonates in the Main Basin. Here we reveal that sublacustrine volcanic events most likely induce the abrupt changes to the geochemistry in the sediment in Lake Kivu. An unprecedented look into the sediment stratigraphy and geochemistry from high-resolution seismic-reflection, and 15N-isotope analyses was conducted in the Main Basin. The results reveal that buoyant hydrothermal plumes caused by subaquatic volcanic activity are a possible trigger for increased biological productivity and organic matter preservation, and that ongoing hydrothermal activity increases the alkalinity in the deepwater, leading to carbonate preservation. The onset of carbonate preservation since the 1970s that is currently observed in the sediment could indicate that hydrothermal discharge has recently increased in the lake.

Vegetation and environmental changes in Northern Anatolia between 134 and 119 ka recorded in Black Sea sediments

This multiproxy study on SE Black Sea sediments provides the first detailed reconstruction of vegetation and environmental history of Northern Anatolia between 134 and 119 ka. Here, the glacial–interglacial transition is characterized by several short-lived alternating cold and warm events preceding a meltwater pulse (~ 130.4–131.7 ka). The latter is reconstructed as a cold arid period correlated to Heinrich event 11. The initial warming is evidenced at ~ 130.4 ka by increased primary productivity in the Black Sea, disappearance of ice-rafted detritus, and spreading of oaks in Anatolia. A Younger Dryas-type event is not identifiable. The Eemian vegetation succession corresponds to the main climatic phases in Europe: i) the Quercus–Juniperus phase (128.7–126.4 ka) indicates a dry continental climate; ii) the Ostrya–Corylus–Quercus–Carpinus phase (126.4–122.9 ka) suggests warm summers, mild winters, and high year-round precipitation; iii) the Fagus–Carpinus phase (122.9–119.5 ka) indicates cooling and high precipitation; and iv) increasing Pinus at ~ 121 ka marks the onset of cooler/drier conditions. Generally, pollen reconstructions suggest altitudinal/latitudinal migrations of vegetation belts in Northern Anatolia during the Eemian caused by increased transport of moisture. The evidence for the wide distribution of Fagus around the Black Sea contrasts with the European records and is likely related to climatic and genetic factors.