Refine
Year of publication
- 2021 (148) (remove)
Document Type
- Article (102)
- Doctoral Thesis (26)
- Postprint (16)
- Monograph/Edited Volume (1)
- Conference Proceeding (1)
- Habilitation Thesis (1)
- Master's Thesis (1)
Is part of the Bibliography
- yes (148) (remove)
Keywords
- PHREEQC (5)
- Earthquake source observations (4)
- reactive transport (4)
- ALOS World 3D (2)
- ASTER GDEM (2)
- Alpen (2)
- Alps (2)
- Arctic Ocean (2)
- Argentina (2)
- Argentinien (2)
Institute
- Institut für Geowissenschaften (148) (remove)
Geodetic studies of crustal deformation using Global Navigation Satellite System (GNSS, earlier commonly referred to as Global Positioning System, GPS) measurements at CSIR-NGRI started in 1995 with the installation of a permanent GNSS station at CSIR-NGRI Hyderabad which later became an International GNSS Service (IGS) site. The CSIR-NGRI started expanding its GNSS networks after 2003 with more focussed studies through installation in the NE India, Himalayan arc, Andaman subduction zone, stable and failed rift regions of India plate. In each instance, these measurements helped in unravelling the geodynamics of the region and seismic hazard assessment, e.g., the discovery of a plate boundary fault in the Indo-Burmese wedge, rate and mode of strain accumulation and its spatial variation in the Garhwal-Kumaun and Kashmir region of the Himalayan arc, the influence of non-tectonic deformation on tectonic deformation in the Himalayan arc, nature of crustal deformation through earthquake cycle in the Andaman Sumatra subduction zone, and localised deformation in the intraplate region and across the paleo rift regions. Besides these, GNSS measurements initiated in the Antarctica region have helped in understanding the plate motion and influence of seasonal variations on deformation. Another important by-product of the GNSS observations is the capabilities of these observations in understanding the ionospheric variations due to earthquake processes and also due to solar eclipse. We summarize these outcomes in this article.
Raman spectroscopic quantification of tetrahedral boron in synthetic aluminum-rich tourmaline
(2021)
The Raman spectra of five B-[4]-bearing tourmalines of different composition synthesized at 700 degrees C/4.0 GPa (including first-time synthesis of Na-Li-B-[4]-tourmaline, Ca-Li-B-[4]-tourmaline, and Ca-bearing square-B-[4]-tourmaline) reveal a strong correlation between the tetrahedral boron content and the summed relative intensity of all OH-stretching bands between 3300-3430 cm(-1). The band shift to low wavenumbers is explained by strong O3-H center dot center dot center dot O5 hydrogen bridge bonding. Applying the regression equation to natural B-[4]-bearing tourmaline from the Koralpe (Austria) reproduces the EMPA-derived value perfectly [EMPA: 0.67(12) B-[4] pfu vs. Raman: 0.66(13) B-[4] pfu]. This demonstrates that Raman spectroscopy provides a fast and easy-to-use tool for the quantification of tetrahedral boron in tourmaline. The knowledge of the amount of tetrahedral boron in tourmaline has important implications for the better understanding and modeling of B-isotope fractionation between tourmaline and fluid/melt, widely used as a tracer of mass transfer processes.
In this study, we investigate numerically the hydro-mechanical behavior of fractured crystalline rock due to one of the five hydraulic stimulations at the Pohang Enhanced Geothermal site in South Korea. We use the commercial code FracMan (Golder Associates) that enables studying hydro-mechanical coupled processes in fractured media in three dimensions combining the finite element method with a discrete fracture network. The software is used to simulate fluid pressure perturbation at fractures during hydraulic stimulation. Our numerical simulation shows that pressure history matching can be obtained by partitioning the treatment into separate phases. This results in adjusted stress-aperture relationships. The evolution of aperture adjustment implies that the stimulation mechanism could be a combination of hydraulic fracturing and shearing. The simulated extent of the 0.01 MPa overpressure contour at the end of the treatment equals to similar to 180 m around the injection point.
Large earthquakes are usually modeled with simple planar fault surfaces or a combination of several planar fault segments. However, in general, earthquakes occur on faults that are non-planar and exhibit significant geometrical variations in both the along-strike and down-dip directions at all spatial scales. Mapping of surface fault ruptures and high-resolution geodetic observations are increasingly revealing complex fault geometries near the surface and accurate locations of aftershocks often indicate geometrical complexities at depth. With better geodetic data and observations of fault ruptures, more details of complex fault geometries can be estimated resulting in more realistic fault models of large earthquakes. To address this topic, we here parametrize non-planar fault geometries with a set of polynomial parameters that allow for both along-strike and down-dip variations in the fault geometry. Our methodology uses Bayesian inference to estimate the non-planar fault parameters from geodetic data, yielding an ensemble of plausible models that characterize the uncertainties of the non-planar fault geometry and the fault slip. The method is demonstrated using synthetic tests considering slip spatially distributed on a single continuous finite non-planar fault surface with varying dip and strike angles both in the down-dip and along-strike directions. The results show that fault-slip estimations can be biased when a simple planar fault geometry is assumed in presence of significant non-planar geometrical variations. Our method can help to model earthquake fault sources in a more realistic way and may be extended to include multiple non-planar fault segments or other geometrical fault complexities.
Earthquake source parameters such as seismic stress drop and corner frequency are observed to vary widely, leading to persistent discussion on potential scaling of stress drop and event size. Physical mechanisms that govern stress drop variations arc difficult to evaluate in nature and are more readily studied in controlled laboratory experiments. We perform two stick-slip experiments on fractured (rough) and cut (smooth) Westerly granite samples to explore fault roughness effects on acoustic emission (AE) source parameters. We separate large stick-slip events that generally saturate the seismic recording system from populations of smaller AE events which are sensitive to fault stresses prior to slip. AE event populations show many similarities to natural seismicity and may be interpreted as laboratory equivalent of natural microseismic events. We then compare the temporal evolution of mechanical data such as measured stress release during slip to temporal changes in stress drops derived from Alis using the spectral ratio technique. We report on two primary observations: (1) In contrast to most case studies for natural earthquakes, we observe a strong increase in seismic stress drop with AE size. (2) The scaling of stress drop with magnitude is governed by fault roughness, whereby the rough fault shows a more rapid increase of the stress drop magnitude relation with progressing large stick-slip events than the smooth fault. The overall range of AE sizes on the rough surface is influenced by both the average grain size and the width of the fault core. The magnitudes of the smallest AE events on smooth faults may also be governed by grain size. However, AEs significantly grow beyond peak roughness and the width of the fault core. Our laboratory tests highlight that source parameters vary substantially in the presence of fault zone heterogeneity (i.e. roughness and narrow grain size distribution), which may affect seismic energy partitioning and static stress drops of small and large AE events.
Fast Holocene slip and localized strain along the Liquiñe-Ofqui strike-slip fault system, Chile
(2021)
In active tectonic settings dominated by strike-slip kinematics, slip partitioning across subparallel faults is a common feature; therefore, assessing the degree of partitioning and strain localization is paramount for seismic hazard assessments. Here, we estimate a slip rate of 18.8 +/- 2.0 mm/year over the past 9.0 +/- 0.1 ka for a single strand of the Liquirie-Ofqui Fault System, which straddles the Main Cordillera in Southern Chile. This Holocene rate accounts for similar to 82% of the trench-parallel component of oblique plate convergence and is similar to million-year estimates integrated over the entire fault system. Our results imply that strain localizes on a single fault at millennial time scale but over longer time scales strain localization is not sustained. The fast millennial slip rate in the absence of historical Mw> 6.5 earthquakes along the Liquine-Ofqui Fault System implies either a component of aseismic slip or Mw similar to 7 earthquakes involving multi-trace ruptures and > 150-year repeat times. Our results have implications for the understanding of strike-slip fault system dynamics within volcanic arcs and seismic hazard assessments.
Rupture directivity, implying a predominant earthquake rupture propagation direction, is typically inferred upon the identification of 2D azimuthal patterns of seismic observations for weak to large earthquakes using surface-monitoring networks. However, the recent increase of 3D monitoring networks deployed in the shallow subsurface and underground laboratories toward the monitoring of microseismicity allows to extend the directivity analysis to 3D modeling, beyond the usual range of magnitudes. The high-quality full waveforms recorded for the largest, decimeter-scale acoustic emission (AE) events during a meter-scale hydraulic fracturing experiment in granites at similar to 410 m depth allow us to resolve the apparent durations observed at each AE sensor to analyze 3D-directivity effects. Unilateral and (asymmetric) bilateral ruptures are then characterized by the introduction of a parameter kappa, representing the angle between the directivity vector and the station vector. While the cloud of AE activity indicates the planes of the hydrofractures, the resolved directivity vectors show off-plane orientations, indicating that rupture planes of microfractures on a scale of centimeters have different geometries. Our results reveal a general alignment of the rupture directivity with the orientation of the minimum horizontal stress, implying that not only the slip direction but also the fracture growth produced by the fluid injections is controlled by the local stress conditions.
Despite advanced seismological techniques, automatic source characterization for microseismic earthquakes remains difficult and challenging since current inversion and modelling of high-frequency signals are complex and time consuming. For real-time applications such as induced seismicity monitoring, the application of standard methods is often not fast enough for true complete real-time information on seismic sources. In this paper, we present an alternative approach based on recent advances in deep learning for rapid source-parameter estimation of microseismic earthquakes. The seismic inversion is represented in compact form by two convolutional neural networks, with individual feature extraction, and a fully connected neural network, for feature aggregation, to simultaneously obtain full moment tensor and spatial location of microseismic sources. Specifically, a multibranch neural network algorithm is trained to encapsulate the information about the relationship between seismic waveforms and underlying point-source mechanisms and locations. The learning-based model allows rapid inversion (within a fraction of second) once input data are available. A key advantage of the algorithm is that it can be trained using synthetic seismic data only, so it is directly applicable to scenarios where there are insufficient real data for training. Moreover, we find that the method is robust with respect to perturbations such as observational noise and data incompleteness (missing stations). We apply the new approach on synthesized and example recorded small magnitude (M <= 1.6) earthquakes at the Hellisheioi geothermal field in the Hengill area, Iceland. For the examined events, the model achieves excellent performance and shows very good agreement with the inverted solutions determined through standard methodology. In this study, we seek to demonstrate that this approach is viable for microseismicity real-time estimation of source parameters and can be integrated into advanced decision-support tools for controlling induced seismicity.
Data driven high resolution modeling and spatial analyses of the COVID-19 pandemic in Germany
(2021)
The SARS-CoV-2 virus has spread around the world with over 100 million infections to date, and currently many countries are fighting the second wave of infections. With neither sufficient vaccination capacity nor effective medication, non-pharmaceutical interventions (NPIs) remain the measure of choice.
However, NPIs place a great burden on society, the mental health of individuals, and economics. Therefore the cost/benefit ratio must be carefully balanced and a target-oriented small-scale implementation of these NPIs could help achieve this balance.
To this end, we introduce a modified SEIRD-class compartment model and parametrize it locally for all 412 districts of Germany. The NPIs are modeled at district level by time varying contact rates. This high spatial resolution makes it possible to apply geostatistical methods to analyse the spatial patterns of the pandemic in Germany and to compare the results of different spatial resolutions.
We find that the modified SEIRD model can successfully be fitted to the COVID-19 cases in German districts, states, and also nationwide. We propose the correlation length as a further measure, besides the weekly incidence rates, to describe the current situation of the epidemic.
Woody plants are expanding into the Arctic in response to the warming climate. The impact on arctic plant communities is not well understood due to the limited knowledge about plant assembly rules.
Records of past plant diversity over long time series are rare. Here, we applied sedimentary ancient DNA metabarcoding targeting the P6 loop of the chloroplast trnL gene to a sediment record from Lake Ilirney (central Chukotka, Far Eastern Russia) covering the last 28 thousand years.
Our results show that forb-rich steppe-tundra and dwarf-shrub tundra dominated during the cold climate before 14 ka, while deciduous erect-shrub tundra was abundant during the warm period since 14 ka. Larix invasion during the late Holocene substantially lagged behind the likely warmest period between 10 and 6 ka, where the vegetation biomass could be highest.
We reveal highest richness during 28-23 ka and a second richness peak during 13-9 ka, with both periods being accompanied by low relative abundance of shrubs. During the cold period before 14 ka, rich plant assemblages were phylogenetically clustered, suggesting low genetic divergence in the assemblages despite the great number of species. This probably originates from environmental filtering along with niche differentiation due to limited resources under harsh environmental conditions. In contrast, during the warmer period after 14 ka, rich plant assemblages were phylogenetically overdispersed.
This results from a high number of species which were found to harbor high genetic divergence, likely originating from an erratic recruitment process in the course of warming. Some of our evidence may be of relevance for inferring future arctic plant assembly rules and diversity changes. By analogy to the past, we expect a lagged response of tree invasion. Plant richness might overshoot in the short term; in the long-term, however, the ongoing expansion of deciduous shrubs will eventually result in a phylogenetically more diverse community.
Semi-distributed hydrological and water quality models are increasingly used as innovative and scientific-based management tools.
However, their application is usually restricted to the gauging stations where they are originally calibrated, limiting their spatial capability.
In this study, the semi-distributed hydrological water quality model HYPE (HYdrological Predictions for the Environment) was tested spatially to represent nitrate-N (NO3- N) and total phosphorus (TP) concentrations and loads of the nested and heterogeneous Selke catchment (463 km(2)) in central Germany.
First, an automatic calibration procedure and uncertainty analysis were conducted using the DiffeRential Evolution Adaptive Metropolis (DREAM) tool to simulate discharge, NO3--N and TP concentrations. A multi-site and multi-objective calibration approach was applied using three main gauging stations, covering the most important hydro-meteorological and physiographical characteristics of the whole catchment. Second, the model's capability was tested to represent further internal stations, which were not initially considered for calibration. Results showed that discharge was well represented by the model at all three main stations during both calibration (1994-1998) and validation (1999-2014) periods with lowest Nash-Sutcliffe Efficiency (NSE) of 0.71 and maximum Percentage BIAS (PBIAS) of 18.0%.
The model was able to reproduce the seasonal dynamics of NO3--N and TP concentrations with low predictive uncertainty at the three main stations, reflected by PBIAS values in the ranges from 16.1% to 6.4% and from 20.0% to 11.5% for NO3--N and TP load simulations, respectively.
At internal stations, the model could represent reasonably well the seasonal variation of nutrient concentrations with PBIAS values in the ranges from 9.0% to 14.2% for NO3--N and from 25.3% to 34.3% for TP concentration simulations.
Overall, results suggested that the spatial validation of a nutrient transport model can be better ensured when a multi-site and multi-objective calibration approach using archetypical gauging stations is implemented.
Further, results revealed that the delineation of sub-catchments should put more focus on hydro-meteorological conditions than on land-use features.
Seeing beyond the outcrop
(2021)
Paleokarst breccias are a common feature of sedimentary rift basins. The Billefjorden Trough in the High Arctic archipelago of Svalbard is an example of such a rift. Here the Carboniferous stratigraphy exhibits intervals of paleokarst breccias formed by gypsum dissolution. In this study we integrate digital outcrop models (DOMs) with a 2D ground penetrating radar (GPR) survey to extrapolate external irregular paleokarst geometries beyond the 2D outcrops. DOMs are obtained through combining a series of overlapping photographs with structure-frommotion photogrammetry, to create mmto dm-resolution georeferenced DOMs. GPR is typically used for surveying the shallow subsurface and relies on detecting the contrasts in electro-magnetic permittivity. We defined three geophysical facies based on their appearance in GPR. By integrating subsurface geophysical data with DOMs we were able to correlate reflection patterns in GPR with outcrop features. The chaotic nature of paleokarst breccias is seen both in outcrop and GPR. Key horizons in outcrop and the GPR profiles allow tying together observations between these methods. Furthermore, we show that this technique expands the twodimensional outcrop surface into a three-dimensional domain, thus complementing, strengthening and extending outcrop interpretations.
Seismic scattering and absorption of oceanic lithospheric S waves in the Eastern North Atlantic
(2021)
The scattering and absorption of high-frequency seismic waves in the oceanic lithosphere is to date only poorly constrained by observations. Such estimates would not only improve our understanding of the propagation of seismic waves, but also unravel the small-scale nature of the lithosphere and its variability. Our study benefits from two exceptional situations: (1) we deployed over 10 months a mid-aperture seismological array in the central part of the Eastern North Atlantic in 5 km water depth and (2) we could observe in total 340 high-frequency (up to 30 Hz) Po and So arrivals with tens to hundreds of seconds long seismic coda from local and regional earthquakes in a wide range of backazimuths and epicentral distances up to 850 km with a travel path in the oceanic lithosphere. Moreover, the array was located about 100 km north of the Gloria fault, defining the plate boundary between the Eurasian and African plates at this location which also allows an investigation of the influence of an abrupt change in lithospheric age (20 Ma in this case) on seismic waves. The waves travel with velocities indicating upper-mantle material. We use So waves and their coda of pre-selected earthquakes to estimate frequency-dependent seismic scattering and intrinsic attenuation parameters. The estimated scattering attenuation coefficients are between 10(-4) and 4 x 10(-5) m(-1) and are typical for the lithosphere or the upper mantle. Furthermore, the total quality factors for So waves below 5 Hz are between 20 and 500 and are well below estimates from previous modelling for observations in the Pacific Ocean. This implies that the Atlantic Ocean is more attenuative for So waves compared to the Pacific Ocean, which is inline with the expected behaviour for the lithospheric structures resulting from the slower spreading rates in the Atlantic Ocean. The results for the analysed events indicate that for frequencies above 3 Hz, intrinsic attenuation is equal to or slightly stronger than scattering attenuation and that the So-wave coda is weakly influenced by the oceanic crust. Both observations are in agreement with the proposed propagation mechanism of scattering in the oceanic mantle lithosphere. Furthermore, we observe an age dependence which shows that an increase in lithospheric age is associated with a decrease in attenuation. However, we also observe a trade-off of this age-dependent effect with either a change in lithospheric thickness or thermal variations, for example due to small-scale upwellings in the upper mantle in the southeast close to Madeira and the Canaries. Moreover, the influence of the nearby Gloria fault is visible in a reduction of the intrinsic attenuation below 3 Hz for estimates across the fault. This is the first study to estimate seismic scattering and absorption parameters of So waves for an area with several hundreds of kilometres radius centred in the Eastern North Atlantic and using them to characterize the nature of the oceanic lithosphere.
The investigation of stresses, faults, structure and seismic hazards requires a good understanding and mapping of earthquake rupture and slip. Constraining the finite source of earthquakes from seismic and geodetic waveforms is challenging because the directional effects of the rupture itself are small and dynamic numerical solutions often include a large number of free parameters. The computational effort is large and therefore difficult to use in an exploratory forward modelling or inversion approach. Here, we use a simplified self-similar fracture model with only a few parameters, where the propagation of the fracture front is decoupled from the calculation of the slip. The approximative method is flexible and computationally efficient. We discuss the strengths and limitations of the model with real-case examples of well-studied earthquakes. These include the M-w 8.3 2015 Illapel, Chile, megathrust earthquake at the plate interface of a subduction zone and examples of continental intraplate strike-slip earthquakes like the M-w 7.1 2016 Kumamoto, Japan, multisegment variable slip event or the M-w 7.5 2018 Palu, Indonesia, supershear earthquake. Despite the simplicity of the model, a large number of observational features ranging from different rupture-front isochrones and slip distributions to directional waveform effects or high slip patches are easy to model. The temporal evolution of slip rate and rise time are derived from the incremental growth of the rupture and the stress drop without imposing other constraints. The new model is fast and implemented in the open-source Python seismology toolbox Pyrocko, ready to study the physics of rupture and to be used in finite source inversions.
Pressure induced structural changes in silicate melts have a great impact on their physico-chemical properties and hence on their behaviour in the deep Earth's interior. In order to gain a deeper understanding we have studied the densification mechanism in multicomponent aluminosilicate glasses (albitic and albit-diopside composition) by means of extended X-ray absorption fine structure spectroscopy coupled to a diamond anvil cell up to 164 GPa. We have monitored the structural modifications from the network-former Ge as well as the network-modifier Sr. Notably, we tracked the evolution of Ge-O and Sr-O bond lengths (RGe-O, RSr-O) and their coordination number with pressure. We show that RGe-O increases strongly up to about 32 GPa, whereas RSr-O increases only slightly up to similar to 26 GPa. We assign these extensions to the increase of the coordination number from 4 to 6 (Ge) and from similar to 6 to at least 9 (Sr). Upon further compression RGe-O and RSr-O exhibit a continuous decrease to the highest probed pressure. These bond contractions, notably of RGe-O, that are continuous and exceed the one observed in pure SiO2 and GeO2, reflect a higher structural flexibility of multi-component glasses compared to those simple systems. Particularly, the high fraction of non-bridging oxygen atoms due to the presence of Na, Sr, Ca, Mg in the studied glasses, favours the simple compression of the highly-coordinated polyhedra of Si and Ge at pressure greater than 30 GPa. This is in strong contrast to pure oxides where cation polyhedral distortions govern the densification mechanism of the glass. The results of this study demonstrate that low field-strength alkali and alkaline earth cations, ubiquitous in deep Earth's melts, have a profound influence on the densification mechanism of glasses. Our results provide important constrains for interpreting the observed low velocity anomalies at the Earth's core-mantle boundary that have been, beyond others, referred to the presence of high-density melts. The hypothesis that non-buoyant melts at the Earth's core-mantle boundary can be formed by peculiar structural transformations in melts leading to higher coordination numbers compared to their crystalline equivalents is not supported from the present observations. The present results rather suggest that if velocity anomalies are to be explained by melts, these likely have considerable differences in chemical composition to the surrounding crystalline phase assemblage.
The within-site variability in site response is the randomness in site response at a given site from different earthquakes and is treated as aleatory variability in current seismic hazard/risk analyses.
In this study, we investigate the single-station variability in linear site response at K-NET and KiK-net stations in Japan using a large number of earthquake recordings.
We found that the standard deviation of the horizontal-to-vertical Fourier spectral ratio at individual sites, that is single-station horizontal-to-vertical spectral ratio (HVSR) sigma sigma(HV,s), approximates the within-site variability in site response quantified using surface-to-borehole spectral ratios (for oscillator frequencies higher than the site fundamental frequency) or empirical ground-motion models.
Based on this finding, we then utilize the single-station HVSR sigma as a convenient tool to study the site-response variability at 697 KiK-net and 1169 K-NET sites.
Our results show that at certain frequencies, stiff, rough and shallow sites, as well as small and local events tend to have a higher sigma(HV,s).
However, when being averaged over different sites, the single-station HVSR sigma, that is sigma(HV), increases gradually with decreasing frequency. In the frequency range of 0.25-25 Hz, sigma(HV) is centred at 0.23-0.43 in ln scales (a linear scale factor of 1.26-1.54) with one standard deviation of less than 0.1. sigma(HV) is quite stable across different tectonic regions, and we present a constant, as well as earthquake magnitude- and distance-dependent sigma(HV) models.
Carbonate minerals are common in both marine and lacustrine records, and are frequently used for paleoenvironmental reconstructions. The sedimentary sequence of the endorheic Dead Sea and its precursors contain aragonite laminae that provide a detailed sedimentary archive of climatic, hydrologic, limnologic and environmental conditions since the Pleistocene. However, the interpretation of these archives requires a detailed understanding of the constraints and mechanisms affecting CaCO3 precipitation, which are still debated. The implications of aragonite precipitation in the Dead Sea and in its late Pleistocene predecessor (Lake Lisan) were investigated in this study by mixing natural and synthetic brines with a synthetic bicarbonate solution that mimics flash-floods composition, with and without the addition of extracellular polymeric substances (EPS). Aragonite precipitation was monitored, and precipitation rates and carbonate yields were calculated and are discussed with respect to modern aquatic environments. The experimental insights on aragonite precipitation are then integrated with microfacies analyses in order to reconstruct and constrain prevailing limnogeological processes and their hydroclimatic drivers under low (interglacial) and high (glacial) lake level stands. Aragonite precipitation took place within days to several weeks after the mixing of the brines with a synthetic bicarbonate solution. Incubation time was proportional to bicarbonate concentration, and precipitation rates were partially influenced by ionic strength. Additionally, extracellular polymeric substances inhibited aragonite precipitation for several months. As for the lake's water budget, our calculations suggest that the precipitation of a typical aragonite lamina (0.5 mm thick) during high lake stand requires unreasonable freshwater inflow from either surface or subsurface sources. This discrepancy can be resolved by considering one or a combination of the following scenarios; (1) discontinuous aragonite deposition over parts of the lake floor; (2) supply of additional carbonate flux (or fluxes) to the lake from aeolian dust and the remobilization and dissolution of dust deposits at the watershed; (3) carbonate production via oxidation of organic carbon by sulfate-reducing bacteria. Altogether, it is suggested that aragonite laminae thickness cannot be directly interpreted for quantitatively reconstructing the hydrological balance for the entire lake, they may still prove valuable for identifying inherent hydroclimatic periodicities at a single site.
Sedimentary ancient DNA-based studies have been used to probe centuries of climate and environmental changes and how they affected cyanobacterial assemblages in temperate lakes. Due to cyanobacteria containing potential bloom-forming and toxin-producing taxa, their approximate reconstruction from sediments is crucial, especially in lakes lacking long-term monitoring data. To extend the resolution of sediment record interpretation, we used high-throughput sequencing, amplicon sequence variant (ASV) analysis, and quantitative PCR to compare pelagic cyanobacterial composition to that in sediment traps (collected monthly) and surface sediments in Lake Tiefer See. Cyanobacterial composition, species richness, and evenness was not significantly different among the pelagic depths, sediment traps and surface sediments (p > 0.05), indicating that the cyanobacteria in the sediments reflected the cyanobacterial assemblage in the water column. However, total cyanobacterial abundances (qPCR) decreased from the metalimnion down the water column. The aggregate-forming (Aphanizomenon) and colony-forming taxa (Snowella) showed pronounced sedimentation. In contrast, Planktothrix was only very poorly represented in sediment traps (meta- and hypolimnion) and surface sediments, despite its highest relative abundance at the thermocline (10 m water depth) during periods of lake stratification (May-October). We conclude that this skewed representation in taxonomic abundances reflects taphonomic processes, which should be considered in future DNA-based paleolimnological investigations.
Since the beginning of the Anthropocene, lacustrine biodiversity has been influenced by climate change and human activities. These factors advance the spread of harmful cyanobacteria in lakes around the world, which affects water quality and impairs the aquatic food chain. In this study, we assessed changes in cyanobacterial community dynamics via sedimentary DNA (sedaDNA) from well-dated lake sediments of Lake Tiefer See, which is part of the Klocksin Lake Chain spanning the last 350 years. Our diversity and community analysis revealed that cyanobacterial communities form clusters according to the presence or absence of varves. Based on distance-based redundancy and variation partitioning analyses (dbRDA and VPA) we identified that intensified lake circulation inferred from vegetation openness reconstructions, delta C-13 data (a proxy for varve preservation) and total nitrogen content were abiotic factors that significantly explained the variation in the reconstructed cyanobacterial community from Lake Tiefer See sediments. Operational taxonomic units (OTUs) assigned to Microcystis sp. and Aphanizomenon sp. were identified as potential eutrophication-driven taxa of growing importance since circa common era (ca. CE) 1920 till present. This result is corroborated by a cyanobacteria lipid biomarker analysis. Furthermore, we suggest that stronger lake circulation as indicated by non-varved sediments favoured the deposition of the non-photosynthetic cyanobacteria sister clade Sericytochromatia, whereas lake bottom anoxia as indicated by subrecent- and recent varves favoured the Melainabacteria in sediments. Our findings highlight the potential of high-resolution amplicon sequencing in investigating the dynamics of past cyanobacterial communities in lake sediments and show that lake circulation, anoxic conditions, and human-induced eutrophication are main factors explaining variations in the cyanobacteria community in Lake Tiefer See during the last 350 years.