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

Combining sedimentological, trace metal (Mn, Mo) and molecular evidence for reconstructing past water-column redox conditions: The example of meromictic Lake Cadagno (Swiss Alps)

Here, we present sedimentological, trace metal, and molecular evidence for tracking bottom water redox-state conditions during the past 12,500 years in nowadays sulfidic and meromictic Lake Cadagno (Switzerland). A 10.5 m long sediment core from the lake covering the Holocene period was investigated for concentration variations of the trace metals Mn and Mo (XRF core scanning and ICP-MS measurements), and for the presence of anoxygenic phototrophic sulfur bacteria (carotenoid pigment analysis and 16S rDNA real time PCR). Our trace metal analysis documents an oxic-intermediate-sulfidic redox-transition period beginning shortly after the lake formation similar to 12.5 kyr ago. The oxic period is characterized by low sedimentary Mn and Mo concentrations, as well as by the absence of any remnants of anoxygenic phototrophic sulfur bacteria. Enhanced accumulation/preservation of Mn (up to 5.6 wt%) in the sediments indicates an intermediate, Mn-enriched oxygenation state with fluctuating redox conditions during a similar to 2300-year long transition interval between similar to 12.1 and 9.8 kyr BP. We propose that the high Mn concentrations are the result of enhanced Mn2+ leaching from the sediments during reducing conditions and subsequent rapid precipitation of Mn-(oxyhydr) oxide minerals during episodic and short-term water-column mixing events mainly due to flood-induced underflows. At 9800 +/- 130 cal yr BP, a rapid transition to fully sulfidic conditions is indicated by the marked enrichment of Mo in the sediments (up to 490 ppm), accompanied by an abrupt drop in Mn concentrations and the increase of molecular biomarkers that indicate the presence of anoxygenic photosynthetic bacteria in the water column. Persistently high Mo concentrations >80 ppm provide evidence that sulfidic conditions prevailed thereafter until modern times, without any lasting hypolimnetic ventilation and reoxygenation. Hence, Lake Cadagno with its persistently stable chemocline offers a framework to study in great temporal detail over similar to 12 kyr the development of phototrophic sulfur bacteria communities and redox processes in a sulfidic environment, possibly depicting analogous conditions in an ancient ocean. Our study underscores the value of combining sedimentological, geochemical, and microbiological approaches to characterize paleo-environmental and -redox conditions in lacustrine and marine settings.

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.