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

Spatial heterogeneity and lake morphology affect diffusive greenhouse gas emission estimates of lakes

Most estimates of diffusive flux (F) of methane (CH4) and carbon dioxide (CO2) from lakes are based on single-point flux chamber measurements or on piston velocity (k) modeled from wind speed and single-point measurements of surface water gas concentrations (C-aq). We analyzed spatial variability of F of CH4 and CO2, as well as C-aq and k in 22 European lakes during late summer. F and k were higher in the lake centers, leading to considerable bias when extrapolating single-point chamber measurements to whole-lake estimates. The ratio of our empirical k estimates to wind speed-modeled k was related to lake size and shape, suggesting a lake morphology effect on the relationship between wind speed and k. This indicates that the error inherent to established wind speed models can be reduced by determining k and C-aq at multiple sites on lakes to calibrate wind speed-modeled k to the local system.

Lake surface water temperatures of European Alpine lakes (1989–2013) based on the Advanced Very High Resolution Radiometer (AVHRR) 1 km data set

Lake water temperature (LWT) is an important driver of lake ecosystems and it has been identified as an indicator of climate change. Consequently, the Global Climate Observing System (GCOS) lists LWT as an essential climate variable. Although for some European lakes long in situ time series of LWT do exist, many lakes are not observed or only on a non-regular basis making these observations insufficient for climate monitoring. Satellite data can provide the information needed. However, only few satellite sensors offer the possibility to analyse time series which cover 25 years or more. The Advanced Very High Resolution Radiometer (AVHRR) is among these and has been flown as a heritage instrument for almost 35 years. It will be carried on for at least ten more years, offering a unique opportunity for satellite-based climate studies. Herein we present a satellite-based lake surface water temperature (LSWT) data set for European water bodies in or near the Alps based on the extensive AVHRR 1 km data record (1989–2013) of the Remote Sensing Research Group at the University of Bern. It has been compiled out of AVHRR/2 (NOAA-07, -09, -11, -14) and AVHRR/3 (NOAA-16, -17, -18, -19 and MetOp-A) data. The high accuracy needed for climate related studies requires careful pre-processing and consideration of the atmospheric state. The LSWT retrieval is based on a simulation-based scheme making use of the Radiative Transfer for TOVS (RTTOV) Version 10 together with ERA-interim reanalysis data from the European Centre for Medium-range Weather Forecasts. The resulting LSWTs were extensively compared with in situ measurements from lakes with various sizes between 14 and 580 km2 and the resulting biases and RMSEs were found to be within the range of −0.5 to 0.6 K and 1.0 to 1.6 K, respectively. The upper limits of the reported errors could be rather attributed to uncertainties in the data comparison between in situ and satellite observations than inaccuracies of the satellite retrieval. An inter-comparison with the standard Moderate-resolution Imaging Spectroradiometer (MODIS) Land Surface Temperature product exhibits RMSEs and biases in the range of 0.6 to 0.9 and −0.5 to 0.2 K, respectively. The cross-platform consistency of the retrieval was found to be within ~ 0.3 K. For one lake, the satellite-derived trend was compared with the trend of in situ measurements and both were found to be similar. Thus, orbital drift is not causing artificial temperature trends in the data set. A comparison with LSWT derived through global sea surface temperature (SST) algorithms shows lower RMSEs and biases for the simulation-based approach. A running project will apply the developed method to retrieve LSWT for all of Europe to derive the climate signal of the last 30 years. The data are available at doi:10.1594/PANGAEA.831007.

The stable carbon isotopic composition of Daphnia ephippia in small, temperate lakes reflects in-lake methane availability

Daphnia can ingest methane-oxidizing bacteria and incorporate methanogenic carbon into their biomass, leading to low stable carbon isotope ratios (expressed as δ13C values) of their tissue. Therefore, δ13C analysis of Daphnia resting eggs (ephippia) in lake sediment records can potentially be used to reconstruct past in-lake availability of methane (CH4). However, detailed multilake studies demonstrating that δ13C values of recently deposited Daphnia ephippia (δ13Cephippia) are systematically related to in-lake CH4 concentrations (CH4aq) are still missing. We measured δ13Cephippia from surface sediments of 15 small lakes in Europe, and compared these values with late-summer CH4aq. δ13Cephippia ranged from −51.6‰ to −25.9‰, and was strongly correlated with CH4aq in the surface water and above the sediment (r −0.73 and −0.77, respectively), whereas a negative rather than the expected positive correlation was found with δ13C values of carbon dioxide (CO2) (r −0.54), and no correlation was observed with CO2aq. At eight sites, offsets between δ13 CCO2 and δ13Cephippia exceeded offsets between δ13 CCO2 and δ13Calgae reported in literature. δ13Cephippia was positively correlated with δ13C values of sedimentary organic matter (r 0.54), but up to 20.7‰ lower in all except one of the lakes (average −6.1‰). We conclude that incorporation of methanogenic carbon prior to ephippia formation must have been widespread by Daphnia in our study lakes, especially those with high CH4aq. Our results suggest a systematic relationship between δ13Cephippia values and CH4aq in small temperate lakes, and that δ13Cephippia analysis on sediment records may provide insights into past changes in in-lake CH4aq.

Lake ice phenology from AVHRR data for European lakes: An automated two-step extraction method

Several lake ice phenology studies from satellite data have been undertaken. However, the availability of long-term lake freeze-thaw-cycles, required to understand this proxy for climate variability and change, is scarce for European lakes. Long time series from space observations are limited to few satellite sensors. Data of the Advanced Very High Resolution Radiometer (AVHRR) are used in account of their unique potential as they offer each day global coverage from the early 1980s expectedly until 2022. An automatic two-step extraction was developed, which makes use of near-infrared reflectance values and thermal infrared derived lake surface water temperatures to extract lake ice phenology dates. In contrast to other studies utilizing thermal infrared, the thresholds are derived from the data itself, making it unnecessary to define arbitrary or lake specific thresholds. Two lakes in the Baltic region and a steppe lake on the Austrian–Hungarian border were selected. The later one was used to test the applicability of the approach to another climatic region for the time period 1990 to 2012. A comparison of the extracted event dates with in situ data provided good agreements of about 10 d mean absolute error. The two-step extraction was found to be applicable for European lakes in different climate regions and could fill existing data gaps in future applications. The extension of the time series to the full AVHRR record length (early 1980 until today) with adequate length for trend estimations would be of interest to assess climate variability and change. Furthermore, the two-step extraction itself is not sensor-specific and could be applied to other sensors with equivalent near- and thermal infrared spectral bands.

Ecosystem regimes and responses in a coupled ancient lake system from MIS 5b to present: the diatom record of lakes Ohrid and Prespa

We reconstruct the aquatic ecosystem interactions since the last interglacial period in the oldest, most diverse, hydrologically connected European lake system, by using palaeolimnological diatom and selected geochemistry data from Lake Ohrid “DEEP site” core and equivalent data from Lake Prespa core, Co1215. Driven by climate forcing, the lakes experienced two adaptive cycles during the last 92 ka: "interglacial and interstadial" and "glacial" cycle. The short-term ecosystems reorganizations, e.g. regime shifts within these cycles substantially differ between the lakes, as evident from the inferred amplitudes of variation. The deeper Lake Ohrid shifted between ultra oligo- and oligotrophic regimes in contrast to the much shallower Lake Prespa, which shifted from a deeper, (oligo-) mesotrophic to a shallower, eutrophic lake and vice versa. Due to the high level of ecosystem stability (e.g. trophic state, lake level), Lake Ohrid appears relatively resistant to external forcing, such as climate and environmental change. Recovering in a relatively short time from major climate change, Lake Prespa is a resilient ecosystem. At the DEEP site, the decoupling between the lakes' response to climate change is marked in the prolonged and gradual changes during the MIS 5/4 and 2/1 transitions. These response differences and the lakes' different physical and chemical properties may limit the influence of Lake Prespa on Lake Ohrid. Regime shifts of Lake Ohrid due to potential hydrological change in Lake Prespa are not evident in the data presented here. Moreover, a complete collapse of the ecosystems functionality and loss of their diatom communities did not happen in either lake for the period presented in the study.

Patterns of morphological changes and hybridization between sympatric whitefish morphs (Coregonus spp.) in a Swiss lake: a role for eutrophication?

Whitefish, genus Coregonus, show exceptional levels of phenotypic diversity with sympatric morphs occurring in numerous postglacial lakes in the northern hemisphere. Here, we studied the effects of human-induced eutrophication on sympatric whitefish morphs in the Swiss lake, Lake Thun. In particular, we addressed the questions whether eutrophication (i) induced hybridization between two ecologically divergent summer-spawning morphs through a loss of environmental heterogeneity, and (ii) induced rapid adaptive morphological changes through changes in the food web structure. Genetic analysis based on 11 microsatellite loci of 282 spawners revealed that the pelagic and the benthic morph represent highly distinct gene pools occurring at different relative proportions on all seven known spawning sites. Gill raker counts, a highly heritable trait, showed nearly discrete distributions for the two morphs. Multilocus genotypes characteristic of the pelagic morph had more gill rakers than genotypes characteristic of benthic morph. Using Bayesian methods, we found indications of recent but limited introgressive hybridization. Comparisons with historical gill raker data yielded median evolutionary rates of 0.24 haldanes and median selection intensities of 0.27 for this trait in both morphs for 1948-2004 suggesting rapid evolution through directional selection at this trait. However, phenotypic plasticity as an alternative explanation for this phenotypic change cannot be discarded. We hypothesize that both the temporal shifts in mean gill raker counts and the recent hybridization reflect responses to changes in the trophic state of the lake induced by pollution in the 1960s, which created novel selection pressures with respect to feeding niches and spawning site preferences.

Population genomic tests of models of adaptive radiation in Lake Victoria Region cichlid fish

Adaptive radiation is usually thought to be associated with speciation, but the evolution of intraspecific polymorphisms without speciation is also possible. The radiation of cichlid fish in Lake Victoria (LV) is perhaps the most impressive example of a recent rapid adaptive radiation, with 600+ very young species. Key questions about its origin remain poorly characterized, such as the importance of speciation versus polymorphism, whether species persist on evolutionary time scales, and if speciation happens more commonly in small isolated or in large connected populations. We used 320 individuals from 105 putative species from Lakes Victoria, Edward, Kivu, Albert, Nabugabo and Saka, in a radiation-wide amplified fragment length polymorphism (AFLP) genome scan to address some of these questions. We demonstrate pervasive signatures of speciation supporting the classical model of adaptive radiation associated with speciation. A positive relationship between the age of lakes and the average genomic differentiation of their species, and a significant fraction of molecular variance explained by above-species level taxonomy suggest the persistence of species on evolutionary time scales, with radiation through sequential speciation rather than a single starburst. Finally the large gene diversity retained from colonization to individual species in every radiation suggests large effective population sizes and makes speciation in small geographical isolates unlikely.

Segregation of species-specific male attractiveness in f(2) hybrid lake Malawi cichlid fish

Among the huge radiations of haplochromine cichlid fish in Lakes Malawi and Victoria, closely related species are often reproductively isolated via female mate choice although viable fertile hybrids can be produced when females are confined only with heterospecific males. We generated F(2) hybrid males from a cross between a pair of closely related sympatric cichlid fish from Lake Malawi. Laboratory mate choice experiments using microsatellite paternity analysis demonstrated that F(2) hybrid males differed significantly in their attractiveness to females of the two parental species, indicating heritable variation in traits involved in mate choice that may contribute to reproductive isolation between these species. We found no significant correlation between male mating success and any measurement of male colour pattern. A simple quantitative genetic model of reproductive isolation suggests that there may be as few as two chromosomal regions controlling species-specific attractiveness. We propose that adaptive radiation of Lake Malawi cichlids could be facilitated by the presence of genes with major effects on mate choice and reproductive isolation.

Age of cichlids: new dates for ancient lake fish radiations

Timing divergence events allow us to infer the conditions under which biodiversity has evolved and gain important insights into the mechanisms driving evolution. Cichlid fishes are a model system for studying speciation and adaptive radiation, yet, we have lacked reliable timescales for their evolution. Phylogenetic reconstructions are consistent with cichlid origins prior to Gondwanan landmass fragmentation 121-165 MYA, considerably earlier than the first known fossil cichlids (Eocene). We examined the timing of cichlid evolution using a relaxed molecular clock calibrated with geological estimates for the ages of 1) Gondwanan fragmentation and 2) cichlid fossils. Timescales of cichlid evolution derived from fossil-dated phylogenies of other bony fishes most closely matched those suggested by Gondwanan breakup calibrations, suggesting the Eocene origins and marine dispersal implied by the cichlid fossil record may be due to its incompleteness. Using Gondwanan calibrations, we found accumulation of genetic diversity within the radiating lineages of the African Lakes Malawi, Victoria and Barombi Mbo, and Palaeolake Makgadikgadi began around or after the time of lake basin formation. These calibrations also suggest Lake Tanganyika was colonized independently by the major radiating cichlid tribes that then began to accumulate genetic diversity thereafter. These results contrast with the widely accepted theory that diversification into major lineages took place within the Tanganyika basin. Together, this evidence suggests that ancient lake habitats have played a key role in generating and maintaining diversity within radiating lineages and also that lakes may have captured preexisting cichlid diversity from multiple sources from which adaptive radiations have evolved.

Unexplained gonad alterations in whitefish (Coregonus spp.) from Lake Thun, Switzerland: levels of persistent organic pollutants in different morphs

Since 2000, a surprisingly high number of macroscopical gonad alterations has been reported in whitefish (Coregonus spp.) from Lake Thun, Switzerland. This unique phenomenon is still unexplained and has received much public attention. As one possible trigger for these effects, the presence of persistent, bioaccumulative and toxic compounds acting as endocrine disruptors in the lake has been discussed. In this study, concentrations of selected persistent organic pollutants were examined in two morphs of whitefish from Lake Thun and their link to the observed abnormalities was investigated. Analyzed compound classes included polychlorinated biphenyls, polychlorinated dibenzo-p-dioxins and dibenzofurans, polychlorinated naphthalenes, polybrominated diphenyl ethers and hexabromocyclododecanes. The target substances were identified in all samples and concentrations of the analyzed compounds were highly correlated among each other. These correlations show that the analyzed substances have the same distribution pattern throughout the lake and that uptake, accumulation and elimination processes are similar. Significant differences in contaminant levels within the samples existed between the two analyzed morphs of whitefish, most likely due to different age, food patterns and growth rate. No difference in contaminant levels was observed between fish with abnormal gonads and fish with normal gonads, suggesting no causal link between the investigated lipophilic organohalogen compounds present in fish and the observed gonad abnormalities in whitefish from Lake Thun. A comparison to existing data shows that concentrations in Lake Thun whitefish are at the lower bound of contaminant levels in whitefish from Swiss lakes or from European waters.

Global change revealed by palaeolimnological records from remote lakes: a review

Over recent decades, palaeolimnological records from remote sites have provided convincing evidence for the onset and development of several facets of global environmental change. Remote lakes, defined here as those occurring in high latitude or high altitude regions, have the advantage of not being overprinted by local anthropogenic processes. As such, many of these sites record broad-scale environmental changes, frequently driven by regime shifts in the Earth system. Here, we review a selection of studies from North America and Europe and discuss their broader implications. The history of investigation has evolved synchronously with the scope and awareness of environmental problems. An initial focus on acid deposition switched to metal and other types of pollutants, then climate change and eventually to atmospheric deposition-fertilising effects. However, none of these topics is independent of the other, and all of them affect ecosystem function and biodiversity in profound ways. Currently, remote lake palaeolimnology is developing unique datasets for each region investigated that benchmark current trends with respect to past, purely natural variability in lake systems. Fostering conceptual and methodological bridges with other environmental disciplines will upturn contribution of remote lake palaeolimnology in solving existing and emerging questions in global change science and planetary stewardship.

Evidence for past variations in methane availability in a Siberian thermokarst lake based on δ13C of chitinous invertebrate remains

Understanding past methane dynamics in arctic wetlands and lakes is crucial for estimating future methane release. Methane fluxes from lake ecosystems have increasingly been studied, yet only few reconstructions of past methane emissions from lakes are available. In this study, we develop an approach to assess changes in methane availability in lakes based on δ13C of chitinous invertebrate remains and apply this to a sediment record from a Siberian thermokarst lake. Diffusive methane fluxes from the surface of ten newly sampled Siberian lakes and seven previously studied Swedish lakes were compared to taxon-specific δ13C values of invertebrate remains from lake surface sediments to investigate whether these invertebrates assimilated 13C-depleted carbon typical for methane. Remains of chironomid larvae of the tribe Orthocladiinae that, in the study lakes, mainly assimilate plant-derived carbon had higher δ13C than other invertebrate groups. δ13C of other invertebrates such as several chironomid groups (Chironomus, Chironomini, Tanytarsini, and Tanypodinae), cladocerans (Daphnia), and ostracods were generally lower. δ13C of Chironomini and Daphnia, and to a lesser extent Tanytarsini was variable in the lakes and lower at sites with higher diffusive methane fluxes. δ13C of Chironomini, Tanytarsini, and Daphnia were correlated significantly with diffusive methane flux in the combined Siberian and Swedish dataset (r = −0.72, p = 0.001, r = −0.53, p = 0.03, and r = −0.81, p < 0.001, respectively), suggesting that δ13C in these invertebrates was affected by methane availability. In a second step, we measured δ13C of invertebrate remains from a sediment record of Lake S1, a shallow thermokarst lake in northeast Siberia. In this record, covering the past ca 1000 years, δ13C of taxa most sensitive to methane availability (Chironomini, Tanytarsini, and Daphnia) was lowest in sediments deposited from ca AD 1250 to ca AD 1500, and after AD 1970, coinciding with warmer climate as indicated by an independent local temperature record. As a consequence the offset in δ13C between methane-sensitive taxa and bulk organic matter was higher in these sections than in other parts of the core. In contrast, δ13C of other invertebrate taxa did not show this trend. Our results suggest higher methane availability in the study lake during warmer periods and that thermokarst lakes can respond dynamically in their methane output to changing environmental conditions.