@article{HannemannEulenfeldKruegeretal.2021, author = {Hannemann, Katrin and Eulenfeld, Tom and Kr{\"u}ger, Frank and Dahm, Torsten}, title = {Seismic scattering and absorption of oceanic lithospheric S waves in the Eastern North Atlantic}, series = {Geophysical journal international}, volume = {229}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggab493}, pages = {948 -- 961}, year = {2021}, abstract = {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.}, language = {en} } @article{DahmHeimannMetzetal.2021, author = {Dahm, Torsten and Heimann, Sebastian and Metz, Malte and Isken, Marius Paul}, title = {A self-similar dynamic rupture model based on the simplified wave-rupture analogy}, series = {Geophysical journal international / the Royal Astronomical Society, the Deutsche Geophysikalische Gesellschaft and the European Geophysical Society}, volume = {225}, journal = {Geophysical journal international / the Royal Astronomical Society, the Deutsche Geophysikalische Gesellschaft and the European Geophysical Society}, number = {3}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggab045}, pages = {1586 -- 1604}, year = {2021}, abstract = {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.}, language = {en} } @article{KrstulovićRosaBiedermannetal.2021, author = {Krstulović, Marija and Rosa, Angelika D. and Biedermann, Nicole and Irifune, Tetsuo and Wilke, Max}, title = {Structural changes in aluminosilicate glasses up to 164 GPa and the role of alkali, alkaline earth cations and alumina in the densification mechanism}, series = {Chemical geology : official journal of the European Association for Geochemistry}, volume = {560}, journal = {Chemical geology : official journal of the European Association for Geochemistry}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0009-2541}, doi = {10.1016/j.chemgeo.2020.119980}, pages = {14}, year = {2021}, abstract = {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.}, language = {en} } @article{ZhuCottonKwaketal.2021, author = {Zhu, Chuanbin and Cotton, Fabrice and Kwak, Dong-Youp and Ji, Kun and Kawase, Hiroshi and Pilz, Marco}, title = {Within-site variability in earthquake site response}, series = {Geophysical journal international}, volume = {229}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggab481}, pages = {1268 -- 1281}, year = {2021}, abstract = {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.}, language = {en} } @article{SimonPapadopoulosGuillemoteauetal.2021, author = {Simon, Fran{\c{c}}ois-Xavier and Papadopoulos, Nikos and Guillemoteau, Julien and Oikonomou, Dimitris and Simirdanis, Kleanthis}, title = {Multi-frequency loop electromagnetic system measurement on shallow offshore archaeological site of Oulos}, series = {ArcheoSciences : revue d'arch{\´e}om{\´e}trie / Groupe des M{\´e}thodes Pluridisciplinaires Contribuant {\`a} l'Arch{\´e}ologie (GMPCA)}, volume = {45}, journal = {ArcheoSciences : revue d'arch{\´e}om{\´e}trie / Groupe des M{\´e}thodes Pluridisciplinaires Contribuant {\`a} l'Arch{\´e}ologie (GMPCA)}, number = {1}, publisher = {Presses Universitaires de Rennes}, address = {Rennes}, issn = {1960-1360}, doi = {10.4000/archeosciences.9690}, pages = {215 -- 218}, year = {2021}, language = {en} } @article{BenDorFlaxLevitanetal.2021, author = {Ben Dor, Yoav and Flax, Tomer and Levitan, Itamar and Enzel, Yehouda and Brauer, Achim and Erel, Yigal}, title = {The paleohydrological implications of aragonite precipitation under contrasting climates in the endorheic Dead Sea and its precursors revealed by experimental investigations}, series = {Chemical geology : official journal of the European Association for Geochemistry}, volume = {576}, journal = {Chemical geology : official journal of the European Association for Geochemistry}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0009-2541}, doi = {10.1016/j.chemgeo.2021.120261}, pages = {18}, year = {2021}, abstract = {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.}, language = {en} } @article{NwosuRoeserYangetal.2021, author = {Nwosu, Ebuka Canisius and Roeser, Patricia Angelika and Yang, Sizhong and Ganzert, Lars and Dellwig, Olaf and Pinkerneil, Sylvia and Brauer, Achim and Dittmann, Elke and Wagner, Dirk and Liebner, Susanne}, title = {From water into sediment-tracing freshwater cyanobacteria via DNA analyses}, series = {Microorganisms : open access journal}, volume = {9}, journal = {Microorganisms : open access journal}, number = {8}, publisher = {MDPI}, address = {Basel}, issn = {2076-2607}, doi = {10.3390/microorganisms9081778}, pages = {20}, year = {2021}, abstract = {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.}, language = {en} } @article{NwosuBrauerKaiseretal.2021, author = {Nwosu, Ebuka Canisius and Brauer, Achim and Kaiser, J{\´e}r{\^o}me and Horn, Fabian and Wagner, Dirk and Liebner, Susanne}, title = {Evaluating sedimentary DNA for tracing changes in cyanobacteria dynamics from sediments spanning the last 350 years of Lake Tiefer See, NE Germany}, series = {Journal of paleolimnology}, volume = {66}, journal = {Journal of paleolimnology}, number = {3}, publisher = {Springer}, address = {Dordrecht}, issn = {0921-2728}, doi = {10.1007/s10933-021-00206-9}, pages = {279 -- 296}, year = {2021}, abstract = {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.}, language = {en} } @article{ZhuTianKempkaetal.2021, author = {Zhu, Zhennan and Tian, Hong and Kempka, Thomas and Jiang, Guosheng and Dou, Bin and Mei, Gang}, title = {Mechanical behaviors of granite after thermal treatment under loading and unloading conditions}, series = {Natural resources research / sponsored by the International Association for Mathematical Geology}, volume = {30}, journal = {Natural resources research / sponsored by the International Association for Mathematical Geology}, number = {3}, publisher = {Springer Science + Business Media B.V.}, address = {New York, NY [u.a.]}, issn = {1520-7439}, doi = {10.1007/s11053-021-09815-7}, pages = {2733 -- 2752}, year = {2021}, abstract = {Understanding the mechanical behaviors of granite after thermal treatment under loading and unloading conditions is of utmost relevance to deep geothermal energy recovery. In the present study, a series of loading and unloading triaxial compression tests (20, 40 and 60 MPa) on granite specimens after exposure to different temperatures (20, 200, 300, 400, 500 and 600 degrees C) was carried out to quantify the combined effects of thermal treatment and loading/unloading stress conditions on granite strength and deformation. Changes in the microstructure of granite exposed to high temperatures were revealed by optical microscopy. The experimental results indicate that both, thermal treatment and loading/unloading stress conditions, degrade the mechanical behaviors and further decrease the carrying capacity of granite. The gradual degradation of the mechanical characteristics of granite after thermal treatment is mainly associated with the evolution of thermal micro-cracks based on optical microscopy observations. The unloading stress state induces the extension of tension cracks parallel to the axial direction, and thus, the mechanical properties are degraded. Temperatures above 400 degrees C have a more significant influence on the mechanical characteristics of granite than the unloading treatment, whereby 400 degrees C can be treated as a threshold temperature for the delineation of significant deterioration. This study is expected to support feasibility and risk assessments by means of providing data for analytical calculations and numerical simulations on granite exposed to high temperatures during geothermal energy extraction.}, language = {en} } @article{JoziNajafabadiHaberlandRybergetal.2021, author = {Jozi Najafabadi, Azam and Haberland, Christian and Ryberg, Trond and Verwater, Vincent F. and Breton, Eline le and Handy, Mark R. and Weber, Michael}, title = {Relocation of earthquakes in the southern and eastern Alps (Austria, Italy) recorded by the dense, temporary SWATH-D network using a Markov chain Monte Carlo inversion}, series = {Solid earth : SE ; an interaktive open access journal of the European Geosciences Union}, volume = {12}, journal = {Solid earth : SE ; an interaktive open access journal of the European Geosciences Union}, number = {5}, publisher = {Copernicus}, address = {G{\"o}ttingen}, organization = {AlpArray Working Grp, AlpArray SWATH-D Working Grp}, issn = {1869-9529}, doi = {10.5194/se-12-1087-2021}, pages = {1087 -- 1109}, year = {2021}, abstract = {In this study, we analyzed a large seismological dataset from temporary and permanent networks in the southern and eastern Alps to establish high-precision hypocenters and 1-D V-P and V-P/V-S models. The waveform data of a subset of local earthquakes with magnitudes in the range of 1-4.2 M-L were recorded by the dense, temporary SWATH-D network and selected stations of the AlpArray network between September 2017 and the end of 2018. The first arrival times of P and S waves of earthquakes are determined by a semi-automatic procedure. We applied a Markov chain Monte Carlo inversion method to simultaneously calculate robust hypocenters, a 1-D velocity model, and station corrections without prior assumptions, such as initial velocity models or earthquake locations. A further advantage of this method is the derivation of the model parameter uncertainties and noise levels of the data. The precision estimates of the localization procedure is checked by inverting a synthetic travel time dataset from a complex 3-D velocity model and by using the real stations and earthquakes geometry. The location accuracy is further investigated by a quarry blast test. The average uncertainties of the locations of the earthquakes are below 500m in their epicenter and similar to 1.7 km in depth. The earthquake distribution reveals seismicity in the upper crust (0-20 km), which is characterized by pronounced clusters along the Alpine frontal thrust, e.g., the Friuli-Venetia (FV) region, the Giudicarie-Lessini (GL) and Schio-Vicenza domains, the Austroalpine nappes, and the Inntal area. Some seismicity also occurs along the Periadriatic Fault. The general pattern of seismicity reflects head-on convergence of the Adriatic indenter with the Alpine orogenic crust. The seismicity in the FV and GL regions is deeper than the modeled frontal thrusts, which we interpret as indication for southward propagation of the southern Alpine deformation front (blind thrusts).}, language = {en} }