Refine
Year of publication
Document Type
- Article (2857)
- Doctoral Thesis (505)
- Postprint (135)
- Other (64)
- Review (52)
- Monograph/Edited Volume (34)
- Preprint (17)
- Habilitation Thesis (12)
- Conference Proceeding (7)
- Master's Thesis (6)
- Part of a Book (3)
- Bachelor Thesis (2)
- Report (2)
- Moving Images (1)
Is part of the Bibliography
- yes (3697) (remove)
Keywords
- climate change (51)
- Holocene (44)
- erosion (28)
- permafrost (28)
- Himalaya (27)
- remote sensing (26)
- Climate change (23)
- Tibetan Plateau (22)
- Earthquake source observations (21)
- Pollen (21)
Institute
- Institut für Geowissenschaften (3697) (remove)
The microstructure of permafrost ground contains clues to its formation and hence its preconditioning to future change.
We applied X-ray computed microtomography (CT) to obtain high-resolution data (Delta x = 50 mu m) of the composition of a 164 cm long permafrost core drilled in a Yedoma upland in north-eastern Siberia.
The CT analysis allowed the microstructures to be directly mapped and volumetric contents of excess ice, gas inclusions, and two distinct sediment types to be quantified. Using laboratory measurements of coarsely resolved core samples, we statistically estimated the composition of the sediment types and used it to indirectly quantify volumetric contents of pore ice, organic matter, and mineral material along the core.
We conclude that CT is a promising method for obtaining physical properties of permafrost cores which opens novel research potentials.
Terrestrial landscape dynamics are dominated by the production, mobilisation, transfer and deposition of sediment. Numerous chemical elements are carried by sediments, making them a key component for ecological processes, as soil constitution and thus plants and animal ecosystems depends on them, and by extension the human species. They are also essential for climate evolution and regulation as marine sedimentation acts as a carbon sink. However, the processes at the origin of their production, mobilisation, transfer can occur suddenly with high-energy content – such as volcanic eruption, mass wasting or flooding events and wildfires – shaking ecosystems and shaping landforms. Besides, in the last era, the human species has shown its ability to disturb landscape dynamics and change sediments cycles. Thus, there is a need for predictive understanding of the processes involved. This relies on un-derstanding of the mechanisms of key processes and their controls, and knowledge of the state and evolution of the Earth’s surface. Classic approaches to these challenges include empirical observations and numerical modeling of geochemical fluxes and surface processes, as well as the study of terrestrial sedimentary archives to better understand the parameters at stake in landscape dynamics and climate changes and the different actions and retroactions between the production, mobilisation, transfer and deposition of sediments which ultimately shape landscapes. Environmental seismology complements these approaches.
Environmental seismology is the science topic investigating the source functions and propagation properties of seismic vibrations triggered by processes happening at or near the Earth’s surface, below and above it - cryosphere, hydrosphere, atmosphere, human habitat, biosphere, etc., to obtain insight into all these physical processes. Indeed, from mass wasting event to rivers, from wild species to hu-man, all these processes are generating seismic waves. Environmental seismology is a rather recent field, with new branches rapidly expanding and at various stages of scientific progress. This thesis is motivated by the goal of learning more on two major natural process hazards (river bedload transport and mass wasting) as well as on human-generated acoustic hazard, while covering the axis of funda-mental research progression, from data exploration and method and theory development to proof of concept, with the twin aims of developing a better understanding of the operation of these specific processes, and of advancing the methods we have at our disposal to study them.
First, I provide a benchmark for assessment of the reliability of existing seismic bedload model inver-sion to retrieve bedload flux from seismic data. Bedload flux measurements are essential to better understand river dynamics, and this can be achieved with environmental seismology. However, due to a lack of well-constrained validation data, the accuracy of the resulting inversions is unknown. I address this gap in Chapter 2.2, reporting a seismic field experiment, and comparing the results to high-quality independent bedload measurements to constrain a seismic bedload model. The study shows that the quality of bedload flux estimates from seismic data strongly depends on the quality of the input data for the model. Direct measurements of relevant parameters, chiefly seismic ground proper-ties needed for the Green's function and the grain size distribution of the moving bedload, considerably improve the model quality over generic approaches using empirical or theoretical functions. I also pro-vide a numerical tool to facilitate the use of water turbulence and bedload seismic inversion models: seismic models for bedload flux and water turbulence require painstaking work to constrain parame-ters describing the ground properties by active seismic study or analysis of passive seismic data, and the grain size distribution via independent measurements. Reasonable predictions can be achieved by using a Monte Carlo approach to optimize the free parameters with respect to the target parameters. The validation of the tool, in Chapter 2.3, with independent measurements of water depth and bedload flux at a study site on the Eshtemoa river in Israël makes it available for reliable use at other sites. The work reported in this chapter has been published in Lagarde et al. 2021 and Dietze et al. 2019b.
In a second study, reported in Chapter 3, I investigate the formation of a failure plane prior to a rock-slide. A better understanding of the dynamics of the preparation phase is essential to determine the timing, volume and mobilization mechanism of a rock slope failure, and this can be achieved with envi-ronmental seismology. I take advantage of a network of seismic stations close to an instable slope recording cracking signals prior to the slope failure, and use a machine learning technique based on hidden Markov models to isolate these signals from the seismic data, retrieving the cumulative number of cracking events in a period of 20 days prior to a large rockslide and 10 days after. The trajectory of the cumulative number of cracks shifts from a rather linear shape in the two weeks prior to the rock-slide to an S-shaped development in the last 27 h before failure. I interpret this change as a switch from initial distributed cracking to localised damage accumulation in the hours prior to the failure. I develop a simple physical model to explain the temporal evolution of crack activity during the S-shape phase, revealing the importance of an internal parameter, the total crack boundary length as the dominant control parameter of failure plane evolution. This study has been published as Lagarde et al. 2023.
Third, I develop a model converting acoustic signals to seismic signals. A part of the acoustic vibrations generated on the Earth’s surface is converted to seismic signals at the ground interface. Consequently, noise pollution may be translated to slope fatigue and rock micro (or macro) fracturing resulting in a degrading effect on landforms. Moreover, this pollution can have negative impacts, such as physical, physiological as well as psychological effects on animal species. At present, the impact of seismic pol-lution generated by acoustic sources is difficult to evaluate. In Chapter 4, I improve and implement a model converting the acoustic pressure generated by a source in the atmosphere to the corresponding seismic signal for a receiver within the ground. The ground is considered to be a porous elastic medium in which wave behaviour can be approximated by the Biot-Stoll model. The model is extended for use of a temporal pressure pulse as an input, and to produce output on a plan 2D map, where the wind effect on the acoustic to seismic coupling can be reproduced. I invest extensive effort in making the model user-friendly, as the project aims at reaching a large audience, comprising, for example, geo-morphologists, biologists and sociologists. Finally, the model is subjected to synthetic testing as well as a qualitative comparison of the predicted ground particle velocity and the seismic signal of a real heli-copter flight as a source of acoustic input.
These studies advance understanding of the operation of specific natural processes in channels and on hillslopes, and bring us closer to designing functioning early warning systems for mass-wasting and flood events. This thesis also raises questions that have not been considered before, such as the con-tribution of human acoustic pollution to the seismic hum and its impact on the natural environment, or the importance of cracks in the self-development of the failure plane prior to slope. Together, these studies question general assumptions usually made regarding the triggering of mass wasting or the hillslope-channel connectivity. Beyond this, the thesis covers the axis of fundamental research progres-sion, from data exploration and method and theory development to proof of concept, and shows how in the rapidly developing field of environmental seismology, an active awareness of progress can help strengthen and accelerate general advances.
Clusty is a new open source toolbox dedicated to earthquake clustering based on waveforms recorded across a network of seismic stations. Its main application is the study of active faults and the detection and characterization of faults and fault networks. By using a density-based clustering approach, earthquakes pertaining to a common fault can be recognized even over long fault segments, and the first-order geometry and extent of active faults can be inferred. Clusty implements multiple techniques to compute a waveform based network similarity from maximum cross-correlation coefficients at multiple stations. The clustering procedure is designed to be transparent and parameters can be easily tuned. It is supported by a number of analysis visualization tools which help to assess the homogeneity within each cluster and the differences among distinct clusters. The toolbox returns graphical representations of the results. A list of representative events and stacked waveforms facilitate further analyses like moment tensor inversion. Results obtained in various frequency bands can be combined to account for large magnitude ranges. Thanks to the simple configuration, the toolbox is easily adaptable to new data sets and to large magnitude ranges. To show the potential of our new toolbox, we apply Clusty to the aftershock sequence of the M-w 6.9 25 October 2018 Zakynthos (Greece) Earthquake. Thanks to the complex tectonic setting at the western termination of the Hellenic Subduction System where multiple faults and faulting styles operate simultaneously, the Zakynthos data set provides an ideal case-study for our clustering analysis toolbox. Our results support the activation of several faults and provide insight into the geometry of faults or fault segments. We identify two large thrust faulting clusters in the vicinity of the main shock and multiple strike-slip clusters to the east, west and south of these clusters. Despite its location within the largest thrust cluster, the main shock does not show a high waveform similarity to any of the clusters. This is consistent with the results of other studies suggesting a complex failure mechanism for the main shock. We propose the existence of conjugated strike-slip faults in the south of the study area. Our waveform similarity based clustering toolbox is able to reveal distinct event clusters which cannot be discriminated based on locations and/or timing only. Additionally, the clustering results allows distinction between fault and auxiliary planes of focal mechanisms and to associate them to known active faults.
Impact of Late Pleistocene climate variability on paleo-erosion rates in the western Himalaya
(2022)
It has been proposed that at short timescales of 10(2)-10(5) yr, climatic variability can explain variations in sediment flux, but in orogens with pronounced climatic gradients rate changes caused by the oscillating efficiency in rainfall, runoff, and/or sediment transport and deposition are still not well-constrained.
To explore landscape responses under variable climatic forcing, we evaluate time windows of prevailing sediment aggradation and related paleo-erosion rates from the southern flanks of the Dhauladhar Range in the western Himalaya.
We compare past and present Be-10-derived erosion rates of well-dated Late Pleistocene fluvial landforms and modern river sediments and reconstruct the sediment aggradation and incision history based on new luminescence data.
Our results document significant variations in erosion rates ranging from 0.1 to 3.4 mm/yr over the Late Pleistocene.
We find that, during times of weak monsoon intensity, the moderately steep areas (hillslope angles of 27 +/- 13 degrees) erode at lower rates of 0.1-0.4 mm/yr compared to steeper (>40 degrees) crestal regions of the Dhauladhar Range that erode at 0.8-1.3 mm/yr.
In contrast, during several millennia of stronger monsoon intensity, both the moderately steep and high slope areas record higher erosion rates (>1-3.4 mm/yr). Lithological clast-count analysis shows that this increase of erosion is focused in the moderately steep areas, where Lesser Himalayan rocks are exposed.
Our data thus highlight the highly non-linear response of climatic forcing on landscape evolution and suggest complex depositional processes and sedimentary signals in downstream areas. (C) 2021 Elsevier B.V. All rights reserved.
Bastnsite [REE(CO3)F] is the main mineral of REE ore deposits in carbonatites. Synthetic bastnsite-like compounds were precipitated from aqueous solutions by many different methods, but previous attempts to model magmatic crystallization of bastnsite from hydrous calciocarbonatite melts were unsuccessful.
Here we present the first experimental evidence that bastnsite and two other REE carbonates, burbankite, and lukechangite, can crystallize from carbonatite melt in the synthetic system La(CO3)F-CaCO3-Na2CO3 at temperatures between 580 and 850 degrees C and a pressure 100 MPa.
The experiments on starting mixtures of reagent-grade CaCO3, Na2CO3, La-2(CO3)(3), and LaF3 were carried out in cold-seal rapid-quench pressure vessels.
The studied system is an isobaric pseudoternary join of a quinary system where CO2 and fluorides act as independent components. Liquidus phases in the run products are calcite, nyerereite, Na carbonate, bastnsite, burbankite solid solution (Na,Ca)(3)(Ca,La)(3)(CO3)(5), and lukechangite Na3La2(CO3)(4)F. Calcite and bastnite form a eutectic in the boundary join La(CO3)F-CaCO3 at 780 +/- 20 degrees C and 58 wt% La(CO3)F. Phase equilibria in the boundary join La(CO3)F-Na2CO3 are complicated by peritectic reaction between Ca-free end-member of burbankite solid solution petersenite (Pet) and lukechangite (Luk) with liquid (L): Na4La2(CO3)(5) (Pet) + NaF (L) = Na3La2(CO3)(4)F (Luk) + Na2CO3 (Nc).
The right-hand-side assemblage becomes stable below 600 +/- 20 degrees C. In ternary mixtures, bastnsite (Bst), burbankite (Bur), and calcite (Cc) are involved in another peritectic reaction: 2La(CO3)F (Bst) + CaCO3 (Cc) + 2Na(2)CO(3) (L) = Na2CaLa2(CO3)(5) (Bur) + 2NaF (L). Burbankite in equilibrium with calcite replaces bastnsite below 730 +/- 20 degrees C. Stable solidus assemblages in the pseudoternary system are: basnsite-burbankite-fluorite-calcite, basnasite-burbankite-fluorite-lukechangite, bastnsite-burbankite-lukechangite, burbankite-lukechangite-nyerereitecalcite, and burbankite-lukechangite-nyerereite-natrite. Addition of 10 wt% Ca-3(PO4)(2) to one of the ternary mixtures resulted in massive crystallization of La-bearing apatite and monazite and complete disappearance of bastnsite and burbankite.
Our results confirm that REE-bearing phosphates are much more stable than carbonates and fluorocarbonates.
Therefore, primary crystallization of the latter from common carbonatite magmas is unlikely. Possible exceptions are carbonatites at Mountain Pass that are characterized by very low P2O5 concentrations (usually at or below 0.5 wt%) and extremely high REE contents in the order of a few weight percent or more. In other carbonatites, bastnsite and burbankite probably crystallized from highly concentrated alkaline carbonate-chloride brines that were found in melt inclusions and are thought to be responsible for widespread fenitization around carbonatite bodies.
The chaotic spatio-temporal electrical activity during life-threatening cardiac arrhythmias like ventricular fibrillation is governed by the dynamics of vortex-like spiral or scroll waves. The organizing centers of these waves are called wave tips (2D) or filaments (3D) and they play a key role in understanding and controlling the complex and chaotic electrical dynamics. Therefore, in many experimental and numerical setups it is required to detect the tips of the observed spiral waves. Most of the currently used methods significantly suffer from the influence of noise and are often adjusted to a specific situation (e.g. a specific numerical cardiac cell model). In this study, we use a specific type of deep neural networks (UNet), for detecting spiral wave tips and show that this approach is robust against the influence of intermediate noise levels. Furthermore, we demonstrate that if the UNet is trained with a pool of numerical cell models, spiral wave tips in unknown cell models can also be detected reliably, suggesting that the UNet can in some sense learn the concept of spiral wave tips in a general way, and thus could also be used in experimental situations in the future (ex-vivo, cell-culture or optogenetic experiments).
The European Alps are amongst the regions with highest glacier mass loss rates over the last decades. Under the threat of ongoing climate change, the ability to predict glacier mass balance changes for water and risk management purposes has become imperative. This raises an urgent need for reliable glacier models. The European Alps do not only host glaciers, but also numerous caves containing carbonate formations, called speleothems. Previous studies have shown that those speleothems also grew during times when the cave was covered by a warm-based glacier. In this thesis, I utilise speleothems from the European Alps as archives of local, environmental conditions related to mountain glacier evolution.
Previous studies have shown that speleothem isotope data from the Alps can be strongly affected by in-cave processes. Therefore, part of this thesis focusses on developing an isotope evolution model, which successfully reproduces differences between contemporaneous growing speleothems. The model is used to propose correction approaches for prior calcite precipitation effects on speleothem oxygen isotopes (δ18O). Applications on speleothem records from caves outside of the Alps demonstrate that corrected δ18O agrees better with other records and climate model simulations.
Existing speleothem growth histories and carbon isotope (δ13C) records from Alpine caves located at different elevations are used to infer soil vs. glacier cover and the thermal regime of the glacier over the last glacial cycle. The compatibility with glacier evolution models is statistically assessed. A general agreement between speleothem δ13C-derived information on soil vs. glacier presence and modelled glacier coverage is found. However, glacier retreat during Marine Isotope Stage (MIS) 3 seems to be underestimated by the model. Furthermore, speleothem data provides evidence of surface temperature above the freezing point which is, however, not fully reproduced by the simulations.
History of glacier cover and their thermal regime is explored for the high-elevation cave system Melchsee-Frutt in the Swiss Alps. Based on new (MIS 9b – MIS 7b, MIS 2) and available speleothem δ13C (MIS 7a – 5d) data, warm-based glacier cover is inferred for MIS 8, 7d, 6, and 2. Also a short period of cold-based ice coverage is found for early MIS 6. In a detailed multi-proxy analysis (δ18O, δ13C, Mg/Ca and Sr/Ca), millennial-scale changes in the glacier-related source of the water infiltrating in the karst during MIS 8 and 7d are found and linked to Northern Hemisphere climate variability.
While speleothem records from high-elevation cave sites in the Alps exhibit huge potential for glacier reconstruction, several limitations remain, which are discussed throughout this thesis. Ultimately, recommendations are given to further leverage subglacial speleothems as an archive of glacier dynamics.
From June to August 2021, we deployed a dense seismic nodal network across the Hengill geothermal area in southwest Iceland to image and characterize faults and high-temperature zones at high resolution.
The nodal network comprised 498 geophone nodes spread across the northern Nesjavellir and southern Hverahlio geothermal fields and was complemented by an existing permanent and temporary backbone seismic network of a total of 44 short-period and broadband stations.
In addition, we recorded distributed acoustic sensing data along two fiber optic telecommunication cables near the Nesjavellir geothermal power plant with commercial interrogators.
During the time of deployment, a vibroseis survey took place around the Nesjavellir power plant.
Here, we describe the network and the recorded datasets.
Furthermore, we showsome initial results that indicate a high data quality and highlight the potential of the seismic records for various follow up studies, such as high-resolution event location to delineate faults and body- and surface-wave tomographies to image the subsurface velocity structure in great detail.
Lakes act as important sinks for inorganic and organic sediment components. However, investigations of sedimentary carbon budgets within glacial lakes are currently absent from Arctic Siberia. The aim of this paper is to provide the first reconstruction of accumulation rates, sediment and carbon budgets from a lacustrine sediment core from Lake Rauchuagytgyn, Chukotka (Arctic Siberia). We combined multiple sediment biogeochemical and sedimentological parameters from a radiocarbon-dated 6.5m sediment core with lake basin hydroacoustic data to derive sediment stratigraphy, sediment volumes and infill budgets. Our results distinguished three principal sediment and carbon accumulation regimes that could be identified across all measured environmental proxies including early Marine Isotope Stage 2 (MIS2) (ca. 29-23.4 ka cal BP), mid-MIS2-early MIS1 (ca. 23.4-11.69 ka cal BP) and the Holocene (ca. 11.69-present). Estimated organic carbon accumulation rates (OCARs) were higher within Holocene sediments (average 3.53 gOCm(-2) a(-1)) than Pleistocene sediments (average 1.08 gOCm(-2) a(-1)) and are similar to those calculated for boreal lakes from Quebec and Finland and Lake Baikal but significantly lower than Siberian thermokarst lakes and Alberta glacial lakes. Using a bootstrapping approach, we estimated the total organic carbon pool to be 0.26 +/- 0.02 Mt and a total sediment pool of 25.7 +/- 1.71 Mt within a hydroacoustically derived sediment volume of ca. 32 990 557m(3). The total organic carbon pool is substantially smaller than Alaskan yedoma, thermokarst lake sediments and Alberta glacial lakes but shares similarities with Finnish boreal lakes. Temporal variability in sediment and carbon accumulation dynamics at Lake Rauchuagytgyn is controlled predominantly by palaeoclimate variation that regulates lake ice-cover dynamics and catchment glacial, fluvial and permafrost processes through time. These processes, in turn, affect catchment and within-lake primary productivity as well as catchment soil development. Spatial differences compared to other lake systems at a trans-regional scale likely relate to the high-latitude, mountainous location of Lake Rauchuagytgyn.
Martian atmospheric spectral end-members retrieval from ExoMars Thermal Infrared (TIRVIM) data
(2022)
Key knowledge about planetary composition can be recovered from the study of thermal infrared spectral range datasets.
This range has a huge diagnostic potential because it contains diagnostic absorptions from a planetary surface and atmosphere. The main goal of this study is to process and interpret the dataset from the Thermal Infrared channel (TIRVIM) which is part of the Atmospheric Chemistry Suite of the ExoMars2016 Trace Gas Orbiter mission to find and characterize dust and water ice clouds in the atmosphere.
The method employed here is based on the application of principal component analysis and target transformation techniques to extract the independent variable components present in the analyzed dataset. Spectral shapes of both atmospheric dust and water ice aerosols have been recovered from the analysis of TIRVIM data.
The comparison between our results with those previously obtained on Thermal Emission Spectrometer (TES) data and with previous analysis on TIRVIM data, validates the methodology here applied, showing that it allows to correctly recover the atmospheric spectral endmembers present in the TIRVIM data.
Moreover, comparison with atmospheric retrievals on PFS, TES and IRIS data, allowed us to assess the temporal stability and homogeneity of dust and water ice components in the Martian atmosphere over a time period of almost 50 years.
The early exhumation history of the Tauern Window in the European Eastern Alps and its surface expression is poorly dated and quantified, partly because thermochronological and provenance information are sparse from the Upper Austrian Northern Alpine Foreland Basin. For the first time, we combine a single-grain double-dating approach (Apatite Fission Track and U-Pb dating) with trace-element geochemistry analysis on the same apatites to reconstruct the provenance and exhumation history of the late Oligocene/early Miocene Eastern Alps. The results from 22 samples from the Chattian to Burdigalian sedimentary infill of the Upper Austrian Northern Alpine Foreland Basin were integrated with a 3D seismic-reflection data set and published stratigraphic reports. Our highly discriminative data set indicates an increasing proportion of apatites (from 6% to 23%) with Sr/Y values <0.1 up-section and an increasing amount of apatites (from 24% to 38%) containing >1,000 ppm light rare-earth elements from Chattian to Burdigalian time. The number of U-Pb ages with acceptable uncertainties increases from 40% to 59% up-section, with mostly late Variscan/Permian ages, while an increasing number of grains (10%-27%) have Eocene or younger apatite fission track cooling ages. The changes in the apatite trace-element geochemistry and U-Pb data mirror increased sediment input from an >= upper amphibolite-facies metamorphic source of late Variscan/Permian age - probably the otztal-Bundschuh nappe system - accompanied by increasing exhumation rates indicated by decreasing apatite fission track lag times. We attribute these changes to the surface response to upright folding and doming in the Penninic units of the future Tauern Window starting at 29-27 Ma. This early period of exhumation (0.3-0.6 mm/a) is triggered by early Adriatic indentation along the Giudicarie Fault System.
Natural gas can be temporarily stored in a variety of underground facilities, such as depleted gas and oil fields, natural aquifers and caverns in salt rocks. Being extensively monitored during operations, these systems provide a favourable opportunity to investigate how pressure varies in time and space and possibly induces/triggers earthquakes on nearby faults. Elaborate and detailed numerical modelling techniques are often applied to study gas reservoirs. Here we show the possibilities and discuss the limitations of a flexible and easily formulated tool that can be straightforwardly applied to simulate temporal pore-pressure variations and study the relation with recorded microseismic events. We use the software POEL (POroELastic diffusion and deformation) which computes the poroelastic response to fluid injection/extraction in a horizontally layered poroelastic structure. We further develop its application to address the presence of vertical impermeable faults bounding the reservoir and of multiple injection/extraction sources. Exploiting available information on the reservoir geometry and physical parameters, and records of injection/extraction rates for a gas reservoir in southern Europe, we perform an extensive parametric study considering different model configurations. Comparing modelled spatiotemporal pore-pressure variations with in situ measurements, we show that the inclusion of vertical impermeable faults provides an improvement in reproducing the observations and results in pore-pressure accumulation near the faults and in a variation of the temporal pore-pressure diffusion pattern. To study the relation between gas storage activity and recorded local microseismicity, we applied different seismicity models based on the estimated porepressure distribution. This analysis helps to understand the spatial distribution of seismicity and its temporal modulation. The results show that the observed microseismicity could be partly linked to the storage activity, but the contribution of tectonic background seismicity cannot be excluded.
Reservoir-triggered seismicity has been observed near dams during construction, impoundment, and cyclic filling in many parts of the earth. In Turkey, the number of dams has increased substantially over the last decade, with Ataturk Dam being the largest dam in Turkey with a total water capacity of 48.7 billion m(3). After the construction of the dam, the monitoring network has improved. Considering earthquakes above the long-term completeness magnitude of M-C = 3.5, the local seismicity rate has substantially increased after the filling of the reservoir. Recently, two damaging earthquakes of M-w 5.5 and M-w 5.1 occurred in the town of Samsat near the Ataturk Reservoir in 2017 and 2018, respectively. In this study, we analyze the spatio-temporal evolution of seismicity and its source properties in relation to the temporal water-level variations and the stresses resulting from surface loading and pore-pressure diffusion. We find that water-level and seismicity rate are anti-correlated, which is explained by the stabilization effect of the gravitational induced stress imposed by water loading on the local faults. On the other hand, we find that the overall effective stress in the seismogenic zone increased over decades due to pore-pressure diffusion, explaining the enhanced background seismicity during recent years. Additionally, we observe a progressive decrease of the Gutenberg-Richter b-value. Our results indicate that the stressing rate finally focused on the region where the two damaging earthquakes occurred in 2017 and 2018.
The accelerating climatic changes and new infrastructure development across the Arctic require more robust risk and environmental assessment, but thus far there is no consistent record of human impact. We provide a first panarctic satellite-based record of expanding infrastructure and anthropogenic impacts along all permafrost affected coasts (100 km buffer, approximate to 6.2 Mio km(2)), named the Sentinel-1/2 derived Arctic Coastal Human Impact (SACHI) dataset. The completeness and thematic content goes beyond traditional satellite based approaches as well as other publicly accessible data sources. Three classes are considered: linear transport infrastructure (roads and railways), buildings, and other impacted area. C-band synthetic aperture radar and multi-spectral information (2016-2020) is exploited within a machine learning framework (gradient boosting machines and deep learning) and combined for retrieval with 10 m nominal resolution. In total, an area of 1243 km(2) constitutes human-built infrastructure as of 2016-2020. Depending on region, SACHI contains 8%-48% more information (human presence) than in OpenStreetMap. 221 (78%) more settlements are identified than in a recently published dataset for this region. 47% is not covered in a global night-time light dataset from 2016. At least 15% (180 km(2)) correspond to new or increased detectable human impact since 2000 according to a Landsat-based normalized difference vegetation index trend comparison within the analysis extent. Most of the expanded presence occurred in Russia, but also some in Canada and US. 31% and 5% of impacted area associated predominantly with oil/gas and mining industry respectively has appeared after 2000. 55% of the identified human impacted area will be shifting to above 0 C-circle ground temperature at two meter depth by 2050 if current permafrost warming trends continue at the pace of the last two decades, highlighting the critical importance to better understand how much and where Arctic infrastructure may become threatened by permafrost thaw.
We revise the conceptual model of calcite varves and present, for the first time, a dual lake monitoring study in two alkaline lakes providing new insights into the seasonal sedimentation processes forming these varves. The study lakes, Tiefer See in NE Germany and Czechowskie in N Poland, have distinct morphology and bathymetry, and therefore, they are ideal to decipher local effects on seasonal deposition. The monitoring setup in both lakes is largely identical and includes instrumental observation of (i) meteorological parameters, (ii) chemical profiling of the lake water column including water sampling, and (iii) sediment trapping at both bi-weekly and monthly intervals. We then compare our monitoring data with varve micro-facies in the sediment record. One main finding is that calcite varves form complex laminae triplets rather than simple couplets as commonly thought. Sedimentation of varve sub-layers in both lakes is largely dependent on the lake mixing dynamics and results from the same seasonality, commencing with diatom blooms in spring turning into a pulse of calcite precipitation in summer and terminating with a re-suspension layer in autumn and winter, composed of calcite patches, plant fragments and benthic diatoms. Despite the common seasonal cycle, the share of each of these depositional phases in the total annual sediment yield is different between the lakes. In Lake Tiefer See calcite sedimentation has the highest yields, whereas in Lake Czechowskie, the so far underestimated re-suspension sub-layer dominates the sediment accumulation. Even in undisturbed varved sediments, re-suspended material becomes integrated in the sediment fabric and makes up an important share of calcite varves. Thus, while the biogeochemical lake cycle defines the varves' autochthonous components and micro-facies, the physical setting plays an important role in determining the varve sub-layers' proportion.
Large earthquakes can increase the amount of water feeding stream flows, raise groundwater levels, and thus grant plant roots more access to water in water-limited environments. We examine growth and photosynthetic responses of Pine plantations to the Maule M-w 8.8 earthquake in headwater catchments of Chile's Coastal Range. We combine high-resolution wood anatomic (lumen area) and biogeochemical (delta 13C of wood cellulose) proxies of daily to weekly tree growth sampled from trees on floodplains and close to ridge lines. We find that, immediately after the earthquake, at least two out of six tree trees on valley floors had increased lumen area and decreased delta 13C, while trees on hillslopes had a reverse trend. Our results indicate a control of soil water on this response, largely consistent with models that predict how enhanced postseismic vertical soil permeability causes groundwater levels to rise on valley floors, but fall along the ridges. Statistical analysis with boosted regression trees indicates that streamflow discharge gained predictive importance for photosynthetic activity on the ridges, but lost importance on the valley floor after the earthquake. We infer that earthquakes may stimulate ecohydrological conditions favoring tree growth over days to weeks by triggering stomatal opening. The weak and short-lived signals that we identified, however, show that such responses are only valid under water-limited, rather than energy-limited tree, growth. Hence, dendrochronological studies targeted at annual resolution may overlook some earthquake effects on tree vitality.
Seafloor spreading at slow rates can be accommodated on large-offset oceanic detachment faults (ODFs), that exhume lower crustal and mantle rocks in footwall domes termed oceanic core complexes (OCCs). Footwall rocks experience large rotation during exhumation, yet important aspects of the kinematics-particularly the relative roles of solid-block rotation and flexure-are not clearly understood. Using a high-resolution numerical model, we explore the exhumation kinematics in the footwall beneath an emergent ODF/OCC. A key feature of the models is that footwall motion is dominated by solid-block rotation, accommodated by the nonplanar, concave-down fault interface. A consequence is that curvature measured along the ODF is representative of a neutral stress configuration, rather than a "bent" one. Instead, it is in the subsequent process of "apparent unbending" that significant flexural stresses are developed in the model footwall. The brittle strain associated with apparent unbending is produced dominantly in extension, beneath the OCC, consistent with earthquake clustering observed in the Trans-Atlantic Geotraverse at the Mid-Atlantic Ridge.
The wetland cover is defined as the spatially homogenous region of a wetland attributed to the underlying biophysical conditions such as vegetation, turbidity, hydric soil, and the amount of water.
Here, we present a novel method to derive the wetland-cover types (WCTs) combining three commonly used multispectral indices, NDVI, MNDWI, and NDTI, in three large Ramsar wetlands located in different geomorphic and climatic settings across India. These wetlands include the Kaabar Tal, a floodplain wetland in east Ganga Plains, Chilika Lagoon, a coastal wetland in eastern India, and Nal Sarovar in semi-arid western India.
The novelty of our approach is that the derived WCTs are stable in space and time, and therefore, a given WCT across different wetlands or within different zones of a large wetland will imply similar underlying biophysical attributes.
The WCTs can therefore provide a novel tool for monitoring and change detection of wetland cover types.
We have automated the proposed WCT algorithm using the Google Earth Engine (GEE) environment and by developing ArcGIS tools. The method can be implemented on any wetland and using any multispectral imagery dataset with visible and NIR bands.
The proposed methodology is simple yet robust and easy to implement and, therefore, holds significant importance in wetland monitoring and management.
The Alpine mountains in central Europe are characterized by a heterogeneous crust accumulating different tectonic units and blocks in close proximity to sedimentary foreland basins. Centroid moment tensor inversion provides insight into the faulting mechanisms of earthquakes and related tectonic processes but is significantly aggravated in such an environment. Thanks to the dense AlpArray seismic network and our flexible bootstrap-based inversion tool Grond, we are able to test different setups with respect to the uncertainties of the obtained moment tensors and centroid locations. We evaluate the influence of frequency bands, azimuthal gaps, input data types, and distance ranges and study the occurrence and reliability of non-double-couple (DC) components. We infer that for most earthquakes (M-w >= 3.3) a combination of time domain full waveforms and frequency domain amplitude spectra in a frequency band of 0.02-0.07 Hz is suitable. Relying on the results of our methodological tests, we perform deviatoric moment tensor (MT) inversions for events with M-w > 3.0. Here, we present 75 solutions for earthquakes between January 2016 and December 2019 and analyze our results in the seismotectonic context of historical earthquakes, seismic activity of the last 3 decades, and GNSS deformation data. We study regions of comparably high seismic activity during the last decades, namely the Western Alps, the region around Lake Garda, and the eastern Southern Alps, as well as clusters further from the study region, i.e., in the northern Dinarides and the Apennines. Seismicity is particularly low in the Eastern Alps and in parts of the Central Alps. We apply a clustering algorithm to focal mechanisms, considering additional mechanisms from existing catalogs. Related to the N-S compressional regime, E-W-to-ENE-WSW-striking thrust faulting is mainly observed in the Friuli area in the eastern Southern Alps. Strike-slip faulting with a similarly oriented pressure axis is observed along the northern margin of the Central Alps and in the northern Dinarides. NW-SE-striking normal faulting is observed in the NW Alps, showing a similar strike direction to normal faulting earthquakes in the Apennines. Both our centroid depths and hypocentral depths in existing catalogs indicate that Alpine seismicity is predominantly very shallow; about 80% of the studied events have depths shallower than 10 km.
The magnitude of earthquakes on continental normal faults rarely exceeds 7.0 Mw. However, because of their vicinity to large population centers they can be highly destructive.
Long recurrence time, relatively small deformations, and limited observations hinder our understanding of the deformation patterns and mechanisms controlling the magnitude of events.
Here, this problem is addressed with 2D thermomechanical modeling of normal fault seismic cycles.
The 2020 Samos, Greece Mw7.0 earthquake is used as an example as it is one of the largest and most studied continental normal fault earthquakes. The modeling approach employs visco-elasto-plastic rheology, compressibility, free surface, and a rate-and-state friction law for the fault.
Modeling of the Samos earthquake suggests the pore fluid pressure ratio on the fault ranges from 0 to 0.7. The model demonstrates that most of the deformation during interseismic and coseismic periods, besides on the fault, occurs in the hanging wall and footwall below the seismogenic part of the fault. The largest vertical surface displacement during the earthquake is the subsidence of the hanging wall in the vicinity of the fault, while the uplift of the footwall and remote part of the hanging wall is significantly smaller.
Modeling of the seismic cycles on normal faults with different setups shows the dependency of the magnitude on the thermal profile and dipping angle of the fault; low heat flow and low dipping angle are favorable conditions for the largest events, while steep normal faults in the areas of high heat flow tend to have the smallest magnitudes.
In a warming Arctic, permafrost-related disturbances, such as retrogressive thaw slumps (RTS), are becoming more abundant and dynamic, with serious implications for permafrost stability and bio-geochemical cycles on local to regional scales. Despite recent advances in the field of earth observation, many of these have remained undetected as RTS are highly dynamic, small, and scattered across the remote permafrost region. Here, we assessed the potential strengths and limitations of using deep learning for the automatic segmentation of RTS using PlanetScope satellite imagery, ArcticDEM and auxiliary datasets. We analyzed the transferability and potential for pan-Arctic upscaling and regional cross-validation, with independent training and validation regions, in six different thaw slump-affected regions in Canada and Russia. We further tested state-of-the-art model architectures (UNet, UNet++, DeepLabv3) and encoder networks to find optimal model configurations for potential upscaling to continental scales. The best deep learning models achieved mixed results from good to very good agreement in four of the six regions (maxIoU: 0.39 to 0.58; Lena River, Horton Delta, Herschel Island, Kolguev Island), while they failed in two regions (Banks Island, Tuktoyaktuk). Of the tested architectures, UNet++ performed the best. The large variance in regional performance highlights the requirement for a sufficient quantity, quality and spatial variability in the training data used for segmenting RTS across diverse permafrost landscapes, in varying environmental conditions. With our highly automated and configurable workflow, we see great potential for the transfer to active RTS clusters (e.g., Peel Plateau) and upscaling to much larger regions.
The formation of the Central Andes dates back to similar to 50 Ma, but its most pronounced episode, including the growth of the Altiplano-Puna Plateau and pulsatile tectonic shortening phases, occurred within the last 25 Ma.
The reason for this evolution remains unexplained. Using geodynamic numerical modeling we infer that the primary cause of the pulses of tectonic shortening and growth of the Central Andes is the changing geometry of the subducted Nazca plate, and particularly the steepening of the mid-mantle slab segment which results in a slowing down of the trench retreat and subsequent increase in shortening of the advancing South America plate.
This steepening first happens after the end of the flat slab episode at similar to 25 Ma, and later during the buckling and stagnation of the slab in the mantle transition zone. Processes that mechanically weaken the lithosphere of the South America plate, as suggested in previous studies, enhance the intensity of the shortening events.
These processes include delamination of the mantle lithosphere and weakening of foreland sediments.
Our new modeling results are consistent with the timing and amplitude of the deformation from geological data in the Central Andes at the Altiplano latitude.
Plain Language Summary
The Central Andes is a subduction-type orogeny that formed as a result of the interaction between the Nazca oceanic plate and the South American continental plate over the last 50 million years. Growth of the Andes is primarily the result of crustal shortening. Nevertheless, "geological" data compiled from previous studies have shown that phases of drastic pulsatile shortening occur at 15 and 5 Ma.
In this study, we used high-resolution 2D numerical geodynamic simulations to investigate the link between oceanic and continental plate dynamics and their interaction. We find that when the oceanic plate steepens in the mantle transition zone, the trench retreat is hindered. Coupled with the weakening of the continental plate through the slab flattening and subsequent delamination of the lithospheric mantle, this leads to pulsatile shortening phases of a magnitude equivalent to that suggested by the data.
Magma-filled dikes may feed erupting fissures that lead to alignments of craters developing at the surface, yet the details of activity and migrating eruptions at the crater row are difficult to monitor and are hardly understood.
The 2021 Tajogaite eruption at the Cumbre Vieja, La Palma (Spain), lasted 85 days and developed a pronounced alignment of craters that may be related to changes within the volcano edifice.
Here, we use COSMO-SkyMed satellite radar data and ground-based time-lapse photographs, offering a high-resolution dataset to explore the locations and characteristics of evolving craters.
Our results show that the craters evolve both gradually and suddenly and can be divided into three main phases. Phase 1, lasting the first 6 weeks of the eruption, was characterized by a NW-SE linear evolution of up to seven craters emerging on the growing cone.
Following two partial collapses of the cone to the northwest and a seismicity increase at depth, Phase 2 started and caused a propagation of the main activity toward the southeastern side, together with the presence of up to 11 craters along this main NW-SE trend. Associated with strong deep and shallow earthquakes, Phase 3 was initiated and continued for the final 2 weeks of the eruption, expressed by the development of up to 18 craters, which became dominant and clustered in the southeastern sector in early December 2021. In Phase 3, a second and oblique alignment and surface fracture was identified.
Our findings that crater and eruption changes coincide together with an increase in seismic activity at depth point to a deep driver leading to crater and morphology changes at the surface.
These also suggest that crater distributions might allow for improved monitoring of changes occurring at depth, and vice versa, such that strong seismicity changes at depth may herald the migration and new formation of craters, which have major implications for the assessment of tephra and lava flow hazards on volcanoes.
The Sea of Marmara is a tectonically active basin that straddles the North Anatolian Fault Zone (NAFZ), a major strike-slip fault that separates the Eurasian and Anatolian tectonic plates. The Main Marmara Fault (MMF), which is part of the NAFZ, contains an approximately 150 km long seismotectonic segment that has not ruptured since 1766. A key question for seismic hazard and risk assessment is whether or not the next rupture along this segment is likely to produce one major earthquake or a series of smaller earthquakes. Geomechanical characteristics such as along-strike variations in rock strength may provide an important control on seismotectonic segmentation. We find that variations in lithospheric strength throughout the Marmara region control the mechanical segmentation of the MMF and help explain its long-term seismotectonic segmentation. In particular, a strong crust that is mechanically coupled to the upper mantle spatially correlates with aseismic patches, where the MMF bends and changes its strike in response to the presence of high-density lower crustal bodies. Between the bends, mechanically weaker crustal domains that are decoupled from the mantle indicate a predominance of creeping. These results are highly relevant for the ongoing debate regarding the characteristics of the Marmara seismic gap, especially in view of the seismic hazard (Mw > 7) in the densely populated Marmara region.
On 12 August 2021, a > 220 s lasting complex earthquake with Mw > 8.2 hit the South Sandwich Trench. Due to its remote location and short interevent times, reported earthquake parameters varied significantly between different international agencies. We studied the complex rupture by combining different seismic source characterization techniques sensitive to different frequency ranges based on teleseismic broadband recordings from 0.001 to 2 Hz, including point and finite fault inversions and the back-projection of high-frequency signals.
We also determined moment tensor solutions for 88 aftershocks. The rupture initiated simultaneously with a rupture equivalent to a M-w 7.6 thrust earthquake in the deep part of the seismogenic zone in the central subduction interface and a shallow megathrust rupture, which propagated unilaterally to the south with a very slow rupture velocity of 1.2 km/s and varying strike following the curvature of the trench.
The slow rupture covered nearly two-thirds of the entire subduction zone length, and with M-w 8.2 released the bulk of the total moment of the whole earthquake.
Tsunami modeling indicates the inferred shallow rupture can explain the tsunami records. The southern segment of the shallow rupture overlaps with another activation of the deeper part of the megathrust equivalent to M-w 7.6. The aftershock distribution confirms the extent and curvature of the rupture. Some mechanisms are consistent with the mainshocks, but many indicate also activation of secondary faults. Rupture velocities and radiated frequencies varied strongly between different stages of the rupture, which might explain the variability of published source parameters.
Plain Language Summary
The earthquake of 12 August 2021 along the deep-sea trench of the South Sandwich Islands in the South Atlantic reached a magnitude of 8.2 and triggered a tsunami. The automatic earthquake parameter determination of different agencies showed very different results shortly after the earthquake and partially underestimated the tsunami potential of the earthquake.
A possible reason was the complex rupture process and that the tsunami was generated by a long and shallow slow slip rupture sandwiched between more conventional fast slip subevents at its northern and southern ends. In addition, the fault surface, which extended over 450 km, was highly curved striking 150 degrees-220 degrees.
We investigated the different components of the seismic wavefields in different frequency ranges and with different methods.
The analysis shows how even complex earthquakes can be deciphered by combining analyzing methods. The comparison with aftershocks and the triggered tsunami waves confirms our model that explains the South Sandwich rupture by four subevents in the plate boundary along the curved deep-sea trench. Here, the depth, rupture velocities, and slip on each segment of the rupture vary considerably.
The method can also be applied to other megathrust earthquakes and help to further improve tsunami warnings in the future.
New middle Miocene to Pliocene (~14–3 Ma) apatite fission track (AFT) cooling ages combined with published K–Ar/Ar–Ar and zircon fission track (ZFT) ages from the Hazara and Swat regions of Pakistan are used to explain the Oligocene to Pliocene structural evolution in the Western Himalaya. The structural model explains the distribution of K–Ar/Ar–Ar ages in three distinct age groups (Proterozoic, Paleozoic-Mesozoic, and Eocene to Oligocene). The Proterozoic to Mesozoic sequence of northern Hazara and Swat experienced elevated temperature and pressure conditions, evident by reset Eocene to Oligocene K–Ar/Ar–Ar hornblende and Eocene to Miocene muscovite ages, caused by Kohistan overthrusting the Indian margin during and after the India–Asia collision. Samples from the Indus syntaxis with Paleo to Mesoproterozoic K–Ar/Ar–Ar hornblende ages and Eocene to Oligocene Ar–Ar muscovite ages show no signs of Cenozoic metamorphism; these samples were thermally imprinted up to the Ar–Ar muscovite closure temperature. Neoproterozoic to Lower Paleozoic rocks from the southern parts of Hazara and Swat show Mesozoic to Oligocene partially reset Ar–Ar muscovite ages and preservation of Ordovician metamorphism. The combined analysis of published K–Ar/Ar–Ar (muscovite), ZFT, and new AFT ages (~14–12 Ma) suggests that the Main Central thrust/Panjal thrust was active from Oligocene to early Miocene (~30–18 Ma), and the Nathia-Gali and Main Boundary thrusts were active from the middle to late Miocene (~14–9 Ma) in the Hazara area. New and published AFT ages (~6–3 Ma) from the Indus syntaxis suggest that early Pliocene tectonic thickening in the hinterland formed the N–S trending Indus anticline, creating an erosional half window in the Main Mantle thrust, forming the Indus syntaxis, and dividing the Main Central thrust sheet into the Hazara and Swat segments.
We propose lacunarity as a novel recurrence quantification measure and illustrate its efficacy to detect dynamical regime transitions which are exhibited by many complex real-world systems. We carry out a recurrence plot-based analysis for different paradigmatic systems and nonlinear empirical data in order to demonstrate the ability of our method to detect dynamical transitions ranging across different temporal scales. It succeeds to distinguish states of varying dynamical complexity in the presence of noise and non-stationarity, even when the time series is of short length. In contrast to traditional recurrence quantifiers, no specification of minimal line lengths is required and geometric features beyond linear structures in the recurrence plot can be accounted for. This makes lacunarity more broadly applicable as a recurrence quantification measure. Lacunarity is usually interpreted as a measure of heterogeneity or translational invariance of an arbitrary spatial pattern. In application to recurrence plots, it quantifies the degree of heterogeneity in the temporal recurrence patterns at all relevant time scales. We demonstrate the potential of the proposed method when applied to empirical data, namely time series of acoustic pressure fluctuations from a turbulent combustor. Recurrence lacunarity captures both the rich variability in dynamical complexity of acoustic pressure fluctuations and shifting time scales encoded in the recurrence plots. Furthermore, it contributes to a better distinction between stable operation and near blowout states of combustors.
The present-day structure of the Eastern Cordillera in NW Argentina is governed by structural and lithological heterogeneities inherited from preceding deformation phases, which influence the localization of newly-formed faults and the inversion of pre-existing structures.
The Salta Rift Basin formed during a Late Jurassic-Cretaceous extensional phase and created a dominant structural and stratigraphic imprint in NW Argentina that is partic-ularly evident within the Eastern Cordillera, where uplift and exhumation have exposed the Salta Group syn-rift succession.
Although in general, the Salta Group rests upon Paleozoic rocks, locally the Tacuru Group forms an intermediate succession, consisting of interfingering eolian sandstones and proximal fault-related conglomerates with a Jurassic maximum depositional age. This succession might be the key to unraveling the Mesozoic history of NW Argentina, prior to the deposition of the Salta Group.
The conglomerates represent the earliest deposits related to extension in the western Lomas de Olmedo sub-basin, which is also documented in predominantly Jurassic zircon (U-Th-Sm)/He cooling ages of the rift shoulders. The detrital zircon U-Pb age signature and sandstone provenance of the Tacuru Group conglomerates differs strongly from the Salta Group syn-rift strata, which show a more regional signal.
These variations and the angularity of the unconformity may be connected to a rotation of the extension direction in the western Lomas de Olmedo sub-basin.
In the Gasht-Masuleh area in the Alborz Mountains, gabbroic magma intruded Palaeozoic metasediments and Mesozoic sediments and crystallised as isotropic and cumulate gabbros. LREE enrichment points to relatively low degrees of mantle melting and depletion of Ti, Nb and Ta relative to primitive mantle points to an arc related component in the magma. Clinopyroxene compositions indicate MORB to arc signatures. U–Pb zircon crystallisation ages of 99.5 ± 0.6 Ma and 99.4 ± 0.6 Ma and phlogopite 40Ar/39Ar ages of 97.1 ± 0.4 Ma, 97.5 ± 0.4 Ma, 97.1 ± 0.1 Ma, within 2σ error, indicate that gabbro intrusion occurred in the (Albian-)Cenomanian (mid-Cretaceous). As active subduction did not take place in the Cretaceous in North Iran, the small volume mafic magmatism in the Gasht-Masuleh area must be due to local, extension-related mantle melting. Melting was most likely caused by far field effects triggered by roll-back of the Neo-Tethys subducting slab. As subduction took place at a distance of ~ 400 km (present distance) from the Alborz Mountains, the observed arc geochemical signatures must be inherited from a previous subduction event and concomitant mantle metasomatism, possibly in combination with contamination of the magma by crustal material.
In recurrence analysis, the tau-recurrence rate encodes the periods of the cycles of the underlying high-dimensional time series. It, thus, plays a similar role to the autocorrelation for scalar time-series in encoding temporal correlations.
However, its Fourier decomposition does not have a clean interpretation. Thus, there is no satisfactory analogue to the power spectrum in recurrence analysis.
We introduce a novel method to decompose the tau-recurrence rate using an over-complete basis of Dirac combs together with sparsity regularization.
We show that this decomposition, the inter-spike spectrum, naturally provides an analogue to the power spectrum for recurrence analysis in the sense that it reveals the dominant periodicities of the underlying time series.
We show that the inter-spike spectrum correctly identifies patterns and transitions in the underlying system in a wide variety of examples and is robust to measurement noise.
Late Miocene-Pliocene onset of fluvial incision of the Cauca River Canyon in the Northern Andes
(2022)
The incision of kilometer-scale canyons into high-standing topography is often used to constrain the surface uplift history of mountain ranges, controlled by tectonic and geodynamic processes.
However, changes in climate may also be responsible for canyon incision. This study deciphers the timing of incision of the similar to 2.5-km-deep Cauca River Canyon in the Central Cordillera of the Northern Andes using the cooling (exhumation) history of rocks from the canyon walls and a regional analysis of channel steepness in rivers.
Ten bedrock samples and one detrital sample were collected on the eastern border of the canyon between 300 m and 2300 m of elevation.
Bedrock and detrital AFT data yield ages from 50 to 38 Ma, while two bed-rock AHe ages from the valley bottom yield ages of 7-6 Ma.
The AHe ages and inverse thermal history models reveal a previously unidentified late Miocene (ca. 7-6 Ma) pulse of exhumation that we interpret as the age of a single incision event that formed the Cauca River Canyon.
We conclude that the Cauca River Canyon was carved as a response to rock uplift in the northern Central Cordillera and propagation of an erosion wave into the mountain range starting in the latest Miocene.
Strain NGK35T is a motile, Gram-stain-negative, rod-shaped (1.0-2.1 mu m long and 0.6-0.8 mu m wide), aerobic bacterium that was isolated from plastic-polluted landfill soil. The strain grew at temperatures between 6 and 37 degrees C (optimum, 28 degrees C), in 0-10 % NaCl (optimum, 1 %) and at pH 6.0-9.5 (optimum, pH 7.5-8.5).
It was positive for cytochrome c oxidase, catalase as well as H2S production, and hydrolysed casein and urea. It used a variety of different carbon sources including citrate, lactate and pyruvate.
The predominant membrane fatty acids were C-16:1 cis9 and C-16:0, followed by C-17:0 cyclo and C-18:1 cis11. The major polar lipids were phosphatidylglycerol and phosphatidylethanolamine, followed by diphosphatidyglycerol. The only quinone was ubiquinone Q-8. Phylogenetic analysis based on 16S rRNA gene sequences indicated that strain NGK35(T) belongs to the genus Paenalcaligenes (family Alcaligenaceae), appearing most closely related to Paenalcaligenes hominis CCUG 53761A(T) (96.90 %) and Paenalcaligenes suwonensis ABC02-12(T) (96.94 %).
The genomic DNA G+C content of strain NGK35(T) was 52.1 mol%. Genome-based calculations (genome-to-genome distance, average nucleotide identity and DNA G+C content) clearly indicated that the isolate represents a novel species within the genus Paenalcaligenes.
Based on phenotypic and molecular characterization, strain NGK35(T) can clearly be differentiated from its phylogenetic neighbours establishing a novel species, for which the name Paenalcaligenes niemegkensis sp. nov. is proposed.
The type strain is NGK35T (=DSM 113270(T)=NCCB 100854(T)).
A marine sediment record from the central Bering Sea, spanning the last 20 thousand years (ka), was studied to unravel the depositional history with regard to terrigenous sediment supply and biogenic sedimentation. Methodic approaches comprised the inference of accumulation rates of siliciclastic and biogenic components, grain-size analysis, and (clay) mineralogy, as well as paleoclimatic modelling. Changes in the depositional history provides insight into land-ocean linkages of paleoenvironmental changes. During the finale of the Last Glacial Maximum, the depositional environment was characterized by hemipelagic background sedimentation. A marked change in the terrigenous sediment provenance during the late Heinrich 1 Stadial (15.7-14.5 ka), indicated by increases in kaolinite and a high glaciofluvial influx of clay, gives evidence of the deglaciation of the Brooks Range in the hinterland of Alaska. This meltwater pulse also stimulated the postglacial onset of biological productivity. Glacial melt implies regional climate warming during a time of widespread cooling on the northern hemisphere. Our simulation experiment with a coupled climate model suggests atmospheric teleconnections to the North Atlantic, with impacts on the dynamics of the Aleutian Low system that gave rise to warmer winters and an early onset of spring during that time. The late deglacial period between 14.5 and 11.0 ka was characterized by enhanced fluvial runoff and biological productivity in the course of climate amelioration, sea-level rise, seasonal sea-ice retreat, and permafrost thaw in the hinterland. The latter processes temporarily stalled during the Younger Dryas stadial (12.9-11.7 ka) and commenced again during the Preboreal (earliest Holocene), after 11.7 ka. High river runoff might have fertilized the Bering Sea and contributed to enhanced upper ocean stratification. Since 11.0 ka, advanced transgression has shifted the coast line and fluvial influence of the Yukon River away from the study site. The opening of the Bering Strait strengthened contour currents along the continental slope, leaving behind winnowed sand-rich sediments through the early to mid-Holocene, with non-deposition occurring since about 6.0 ka.
Rapid humidity changes across the Northern South China Sea during the last similar to 40 kyrs
(2021)
A key aspect of East Asian climate is its summer monsoonal system which influences nearly one-third of the world's population. Recent results indicate that the primary response of the East Asian summer monsoon (EASM) to anthropogenic forced climate warming may be a shift in geographical range instead of an intensity change, which would lead to spatial coexistence of floods and droughts over southeastern Asia. The predicted EASM variability in the future has made it paramount to study its past changes and the associated tempo-spatial pattern of aridity and humidity in its purview. In order to decipher past changes in EASM, we applied a multi-proxy geochemical approach to the sediment core ORI-891-16-P1 located in the northern South China Sea. The position of this sediment core on top of a seamount makes it uniquely sensitive to changes in the terrigenous input into northern South China Sea unbiased by sea level-induced downslope transport processes. Utilizing the ln(Ti/Ca) ratio throughout the sediment sequence we trace terrigenous influx changes reflecting EASM prevalence during the last similar to 40 kyrs. Based on the comparison of our results to previous studies we infer that the Last Glacial Maximum (LGM; similar to 20 ka BP) was characterized by a steep N-S humidity gradient. This spatial pattern was in line with a southward shift or contraction of the summer monsoonal trough of 10-15 degrees from its current position toward the centre of the South China Sea. Superimposed on orbital time scale fluctuations we also find strong indication of millennial-scale variability related to Heinrich Stadials. The impact of Heinrich Stadials on the EASM seems amplified during insolation minima, while high summer insolation seems to buffer the monsoonal system to such perturbations. We infer that (i) the humidity-aridity distribution during the LGM mimics predictions of the proposed future EASM configuration, and (ii) that the sensitivity of the EASM to weakening in the Atlantic Meridional Overturning Circulation is the strongest since the last glacial.
With the present study, we introduce a fast and robust method to calculate the source displacement spectra of small earthquakes on a local to regional scale. The work is based on the publicly available Qopen method of full envelope inversion, which is further tuned for the given purpose. Important source parameters-seismic moment, moment magnitude, corner frequency, and high-frequency fall off-are determined from the source spectra by fitting a simple earthquake source model. The method is demonstrated by means of a data set comprising the 2018 West Bohemia earthquake swarm. We report moment magnitudes, corner frequencies, and centroid moment tensors inverted from short-period body waves with the Grond package for all earthquakes with a local magnitude larger than 1.8. Moment magnitudes calculated by envelope inversion show a very good agreement to moment magnitudes resulting from the probabilisitc moment tensor inversion. Furthermore, source displacement spectra from envelope inversion show a good agreement with spectra obtained by multiple taper analysis of the direct onsets of body waves but are not affected by the large scatter of the second. The seismic moments obtained with the envelope inversion scale with corner frequencies according to M-0 proportional to f(c)(-4.7). Earthquakes of the present data set result in a smaller stress drop for smaller magnitudes. Self-similarity of earthquake rupture is not observed. In addition, we report frequency-dependent site amplification at the used stations.
The Shanderman lamprophyre dykes crop out in the western part of the Alborz Mountains (Talesh).
These rocks are classified as camptonites, composed of primary olivine, Ti-rich diopside, kaersutite, biotite, plagioclase, K-feldspar, and minor Ti-rich spinels, magnetite, pentlandite-pyrrhotite/chalcopyrite, and powellite-scheelite. Secondary analcime-wairakite, serpentines, and prehnite are common minor minerals within the studied rocks.
Olivine, Ti-rich diopside, spinel, and amphibole show distinct chemical zoning. Spinels display a core-to-rim decrease in Cr2O3, MgO, and Al2O3 concentrations and an increase in TiO2 and FeOT (total Fe as FeO), reflecting the oxidation state increase due to hydrothermal fluid influx. Low SiO2 contents (< 42 wt%), high MgO (12.44 to 13.98 wt%), and Fe2O3T (12.76 to 13.43 wt%), Cr (318-537 mu g/g) and Ni (231-327 mu g/g) contents indicate the ultrabasic nature of the rocks.
The samples show potassic character (2.1-2.8 wt% K2O), along with elevated LREE and LILE, and also exhibit minor positive Eu anomalies (Eu/Eu* = 1.09 to 1.20).
Olivine-spinel geothermometry indicates a maximum crystallization temperature of 1227 degrees C (ave. 988 degrees C +/- 65 degrees C).
Exsolution of pentlandite-pyrrhotite/chalcopyrite solid solutions occurred during magma cooling and crystallization. At lower temperatures, analcime-wairakite and prehnite partially replaced plagioclases.
The geochemical modeling of the rocks indicates the Shanderman lamprophyre magmas were derived from low-grade melting (< 5%) of amphibole-bearing garnet lherzolite source without or with very few phlogopites.
The primary magma of Shanderman lamprophyres was derived from a depth of similar to 135 km by partial melting of a metasomatized mantle source in a post-collisional environment.
Ground subsidence caused by natural or anthropogenic processes affects major urban areas worldwide. Sinkhole formation and infrastructure fractures have intensified in the federal capital of Maceio (Alagoas, Brazil) since early 2018, forcing authorities to relocate affected residents and place buildings under demolition. In this study, we present a 16-year history (2004-2020) of surface displacement, which shows precursory deformations in 2004-2005, reaching a maximum cumulative subsidence of approximately 200 cm near the Mundau Lagoon coast in November 2020. By integrating the displacement observations with numerical source modelling, we suggest that extensive subsidence can be primarily associated with the removal of localized, deep-seated material at the location and depth where salt is mined. We discuss the accelerating subsidence rates, influence of severe precipitation events on the aforementioned geological instability, and related hazards. This study suggests that feedback destabilization mechanisms may arise in evaporite systems due to anthropogenic activities, fostering enhanced and complex superficial ground deformation.
Geochemical homogeneity in shale is often assumed when tracing subsurface fluids and characterizing sedimentary basins. This study presents measurements of the bulk gas composition, stable isotopes, and noble gas volume fraction and isotopes for shale gas samples collected from gas wells in the Wufeng-Longmaxi Shale, the southern Sichuan Basin, China. The dryness [C-1 /(C-2 + C-3)] ranging from 166.3 to 251.2, combined with delta C-13(1) and delta DC1 that vary from -28.8 to -27.3 parts per thousand and - 153 to -145 parts per thousand, respectively, point to a late mature thermogenic origin of hydrocarbon gas. He-3/He-4 ratios of gas samples are around 0.01 times the air value suggesting dominantly crust-derived He. Ne-21/Ne-22 and Ar-40/Ar-36 ratios of many gas samples are higher than the corresponding air values indicating the mixing of crustal and atmospheric noble gases. Multiple dichotomous patterns are observed in noble gas signatures of forelimb and backlimb samples, and depression and crest samples. Ne-20/Ne-22 ratios of some crest samples are higher than that of depression samples in the backlimb, pointing to the presence of diffusion-driven fractionation that is likely caused by the long-distance migration from depression to crest. Elemental ratios of air-derived noble gas isotopes - Ne-22/Ar-36, Kr-84/Ar-36, and Xe-132/Ar-36 are compared to the recharge water values, suggesting the interactions of oil, gas, and water phases in the shale over geologic time. Forelimb samples generally display older ages than backlimb samples, indicating a larger flux of external radiogenic He-4 due to the higher density of deep faults in the forelimb area caused by the basementinvolved deformation. The basement-involved deformation also causes pore collapse especially in the forelimb leading to a lower porosity that results in a more pristine noble gas signature in the forelimb due to the reduced impact of younger recharge water.
Equilibrium mass-dependent ("stable") isotopic fractionation of an element during magmatic processes is driven by a contrast in bonding environment between minerals and silicate melt, which is expressed as an isotopic fractionation factor.
A quantitative understanding of such isotopic fractionation factors is vital to interpret observed isotopic variations in magmatic rocks.
It is well known that the local environment and the bond strength of an element dictate the sign and magnitude of isotopic fractionation between minerals, but it is uncertain how the structure and chemical composition of a silicate melt can affect mineral-melt isotopic fractionation factors.
To explore this, we studied the coordination environment of nickel (Ni) in different silicate glasses using extended X-ray absorption fine structure (EXAFS) measurements at the German synchrotron X-ray source (DESY).
We determined -Ni-O bond lengths in a suite of synthetic but near-natural silicate glasses using EXAFS and found that the former vary systematically with melt alkalinity, which is best described by the parameter ln[1 + (Na + K)/Ca]. With increasing melt alkalinity, Ni occupies more IV-fold coordinated sites, which are associated with a shorter -Ni-O bond length. Next, we use the ionic model, which allows to predict isotopic fractionation factors based on the difference in bond length between two phases.
We find that more alkaline melts have a stronger preference for the heavier isotopes of Ni than less alkaline melts. This implies that the magnitude of mineral-melt Ni isotope fractionation factors, for instance between olivine and melt, will depend on the alkalinity of the melt.
At magmatic temperatures, however, the variation in fractionation factors caused by melt alkalinity will rarely exceed 0.05 parts per thousand and is thus mostly negligible, in particular in the realm of basaltic melt compositions. Nevertheless, the relationship between melt alkalinity and fractionation factor reported here can be used to extrapolate empirical data for mineral-melt Ni isotope fractionation factors, once such data become available, to the full range of magma compositions on Earth and other Solar System bodies.
Origin and migration of fluoride in the area of the Aluto Volcanic Complex (Main Ethiopian Rift)
(2022)
Fluoride-enriched ground and surface waters represent a major health risk for the local population in many areas along the East African Rift.
The present study investigates the origin of fluoride and the reason for its accumulation in the rift waters, following two hypotheses: (i) fluid-rock-interactions release fluoride from minerals into the water and (ii) magmatic-derived fluoride-containing liquids and gases migrate along permeable fault zones until they dissolve in ground-and surface water or be released to the atmosphere.
Rock-, gas, water-, soil-, and plant samples were collected from the area within and close by the Aluto Volcanic Complex, which is part of the Main Ethiopian Rift. Most analyzed waters showed fluoride concentrations above the drinking water limit (> 1.5 mg/L) with the highest values in hot springs (up to 70 mg/L) and the geothermal well (76 mg/L).
In the solid phase, a high fluoride content was found in acid volcanic rocks (ignimbrite: 4391 ppm; rhyolite: 3248 ppm) as well as in pumice (up to 1955 ppm). The fluoride content of soil samples collected within the volcanic complex varied between 82 and 1036 ppm, whereas former lake sediments from outside the Aluto Volcanic Complex showed higher fluoride contents ranging from 674 to 8171 ppm. Identified fluoride-rich minerals are various amphiboles (about 3 wt.-% F-), fluor-apatite (about 3 wt % F-), minerals of the fluorcaphite group (about 5 wt.-% F-), parisite (up to 9 wt.-% F-), and fluorite (CaF2).
Elevated fluoride concentrations were also measured in some gas samples from fumaroles (up to 50 ppm) and in plant samples collected next to the fumaroles (up to 65 ppm). Leaching experiments of solid samples with deionized water demonstrated that fluoride can easily be mobilized from pumice and sediments but hardly from consolidated rocks.
This fluoride release increased with temperature (up to 150 C) and correlated roughly with dis-solved silica indicating the binding of some fluoride to the amorphous or weakly crystalline silica fraction. Based on these results it was concluded that fluoride migrates via different pathways through the environment: At the depth of the magma chamber during magmatic differentiation fluoride enriches initially in the magmatic melt and accumulates in some late-crystallizing minerals of igneous rocks such as fluorite or mica. Upon volcanic eruption fluoride is predominantly incorporated in the glass (ignimbrite) and ash phase (pumice).
On the surface, these rocks are exposed to weathering and fluoride leaches partly out into the aqueous phase. Soft and porous rocks such as pumice release fluoride first whereas hard extrusive/igneous rocks are less prone to weathering and retain the fluoride. Pumice and (lake) sediments might act both, as source and as sink for fluoride in the area.
Although some fluoride might drain from the surface into the deeper subsoil, we conclude that magmatic fluids (liquid and gas) contribute more substantially to the overall fluoride ground water content because (i) of the much higher fluoride content in deep geothermal waters and hot springs as compared to surface-near waters; (ii) active geothermal surface manifestations located along fault zones indicate that fluids migrate from deep magmatic intrusions (as gas and liquid) towards the surface, where fluoride dissolves in groundwater; (iii) and the good correlation between bicarbonate (deriving from dissolution of magmatic CO2) and fluoride content in all analyzed water samples.
Forward stratigraphic modelling is a fast-developing modelling approach, used to test conceptual models, and predict stratigraphic architecture and depositional facies from basin to reservoir scales.
Published subsurface applications demonstrate its added value by integrating multidisciplinary data as well as geological concepts into its constraints.
When applied to carbonate depositional systems, composed of multiple sediment factories, the cooperating and interdependent production mechanisms remain poorly studied.
By applying the technique to a well-studied section of the Maldives carbonate platform, a specific model design-adapted to the geological age and setting, and constrained by available data-sheds light on the interaction of its carbonate producers.
The results yield a naturalistic depositional facies distribution and offer insight in the changing relationship between biotic communities during the platform evolution.
After calibration, the reference model unequivocally links the formerly proposed genetic model to the seismostratigraphic architecture.
Furthermore, the results show how environmental changes (seemingly of secondary impact compared to changes in physical accommodation in the stratigraphic record) can induce substantial fluctuations in carbonate production rates of biotic communities, affect the ecological accommodation, and thus impact the platform architecture. Therefore, it is crucial to treat carbonate production rates during periods of environmental change as variables with associated uncertainties in a forward stratigraphic model setup.
Floodplain wetlands are critical for sustaining various ecological and hydrological functions in a riverine environment. Severe anthropogenic alterations and human occupation of floodplains have threatened these wetlands in several parts of the world. A major handicap in designing sustainable restoration and monitoring strategies for these wetlands is the lack of scientific process-based understanding and information on the basin-scale controls of their degradation. Here, we offer a novel approach to integrate the connectivity of the wetlands with the surrounding landscape along with other attributes such as stream density, hydrometeorological parameters, and groundwater dynamics to explain their degradation and then to prioritise them for restoration and monitoring.
We hypothesise that the best possible connectivity scenario for the existence of a wetland would be if (a) the wetland has a high connectivity with its upslope area, and (b) the wetland has a low connectivity with its downslope region.
The first condition ensures the flow of water into the wetland and the second condition allows longer water residence time in the wetland. Accordingly, we define four connectivity-based wetland health scenarios-good, no impact, bad, and worst.
We have implemented the proposed method in 3226 wetlands in the Ramganga Basin in north India. Further, we have applied specific selection criteria, such as distance from the nearest stream and stream density, to prioritise the wetlands for restoration and monitoring.
We conclude that the connectivity analysis offers a quick process-based assessment of wetlands' health status and serves as an important criterion to prioritise the wetlands for developing appropriate management strategies.
This research demonstrated the application of hydrochemical data and stable water isotopes of delta O-18 and delta D (or delta 2H) for evaluating the relationship between surface water in Lake Bosumtwi and the underlying groundwater system.
It aimed at determining the presence or absence of a hydraulic relationship and for evaluating the possible direction of flow at the interface between the two reservoirs. The study also estimated evaporative losses of infiltrating rainwater as it transits the unsaturated zone and provided important information on the hydrological processes prevailing in the area.
The results of Q-Mode hierarchical cluster analysis (HCA) clearly differentiate the lake water from the groundwater based on their spatial relationship.
These results indicated that groundwater recharge occurs on the hilltops of the crater, where it is slightly acidic with low levels of dissolved ions, characterised by short residence time and rapid unrestricted vertical infiltration and recharge.
The groundwater becomes more mineralized with longer contact times and deeper circulation with the host rock, while it flows from the recharge areas towards the lake at lower elevations.
Analyses of delta O-18 and delta D showed a high evaporation rate on the lake surface (90%) with a significant evaporative enrichment, whereas groundwater showed no significant isotopic variations.
Thus suggesting that the aquifers have been recharged by recent meteoric water that has undergone some evaporative enrichment since the study established an evaporation rate of water infiltrating the unsaturated zone ranging from 54 to 60%. Both reservoirs do not appear to be hydraulically connected, and where such a connection exists, it is expected to favour the lake.
A circular, single-contig Methanobacterium sp. metagenome-assembled genome (MAG) was recovered from high-CO2 enrichments inoculated with drill core material from the tectonic Eger Rift terrestrial subsurface.
Annotation of the recovered MAG highlighted putative methanogenesis genes, providing valuable information on archaeal activity in the deep biosphere.
Refractive index provides fundamental insights into the electronic structure of materials. At high pressure, however, the determination of refractive index and its wavelength dispersion is challenging, which limits our understanding of how physical properties of even simple materials, such as MgO, evolve with pressure.
Here, we report on the measurement of room-temperature refractive index of MgO up to similar to 140 GPa.
The refractive index of MgO at 600 nm decreases by similar to 2.4% from similar to 1.737 at 1 atm to similar to 1.696 (+/- 0.017) at similar to 140 GPa.
Despite the index at 600 nm is essentially pressure independent, the absolute wavelength dispersion of the refractive index at 550-870 nm decreases by similar to 28% from similar to 0.015 at 1 atm to similar to 0.011 (+/- 8.04 x 10(-4)) at similar to 103 GPa.
Single-effective-oscillator analysis of our refractive index data suggests that the bandgap of MgO increases by similar to 1.1 eV from 7.4 eV at 1 atm to similar to 8.5 (+/- 0.6) eV at similar to 103 GPa.
The Gram-negative bacterium Paenalcaligenes niemegkensis NGK35(T) was isolated from plastic debris in an abandoned landfill.
It has the ability to grow on polyethylene and hexadecane as the sole carbon sources. Here, we report the corresponding draft genome, which contains 3.66 Mbp and is characterized by a G+C content of 52.1%.
International migration patterns, at the global level, can to a large extent be explained through economic factors in origin and destination countries.
On the other hand, it has been shown that global climate change is likely to affect economic development over the coming decades.
Here, we demonstrate how these future climate impacts on national income levels could alter the global migration landscape.
Using an empirically calibrated global migration model, we investigate two separate mechanisms. The first is through destination-country income, which has been shown consistently to have a positive effect on immigration.
As countries' income levels relative to each other are projected to change in the future both due to different rates of economic growth and due to different levels of climate change impacts, the relative distribution of immigration across destination countries also changes as a result, all else being equal.
Second, emigration rates have been found to have a complex, inverted U-shaped dependence on origin-country income.
Given the available migration flow data, it is unclear whether this dependence-found in spatio-temporal panel data-also pertains to changes in a given migration flow over time. If it does, then climate change will additionally affect migration patterns through origin countries' emigration rates, as the relative and absolute positions of countries on the migration "hump" change.
We illustrate these different possibilities, and the corresponding effects of 3 degrees C global warming (above pre-industrial) on global migration patterns, using climate model projections and two different methods for estimating climate change effects on macroeconomic development.
Seismicity models are probabilistic forecasts of earthquake rates to support seismic hazard assessment.
Physics-based models allow extrapolating previously unsampled parameter ranges and enable conclusions on underlying tectonic or human-induced processes.
The Coulomb Failure (CF) and the rate-and-state (RS) models are two widely used physics-based seismicity models both assuming pre-existing populations of faults responding to Coulomb stress changes.
The CF model depends on the absolute Coulomb stress and assumes instantaneous triggering if stress exceeds a threshold, while the RS model only depends on stress changes.
Both models can predict background earthquake rates and time-dependent stress effects, but the RS model with its three independent parameters can additionally explain delayed aftershock triggering.
This study introduces a modified CF model where the instantaneous triggering is replaced by a mean time-to-failure depending on the absolute stress value.
For the specific choice of an exponential dependence on stress and a stationary initial seismicity rate, we show that the model leads to identical results as the RS model and reproduces the Omori-Utsu relation for aftershock decays as well stress-shadowing effects.
Thus, both CF and RS models can be seen as special cases of the new model. However, the new stress response model can also account for subcritical initial stress conditions and alternative functions of the mean time-to-failure depending on the problem and fracture mode.
The most profound consequences of the presence of Ca-Mg carbonates (CaCO3-MgCO3) in the Earth's upper mantle may be to lower the melting temperatures of the mantle and control the melt composition.
Low-degree partial melting of a carbonate-bearing mantle produces CO2-rich, silica-poor melts compositionally imposed by the melting relations of carbonates.
Thus, understanding the melting relations in the CaCO3-MgCO3 system facilitates the interpretation of natural carbonate-bearing silicate systems.
We report the melting relations of the CaCO3-MgCO3 system and the partition coefficient of trace elements between carbonates and carbonate melt from experiments at high pressure (6 and 9 GPa) and temperature (1300-1800 degrees C) using a rocking multi-anvil press. In the absence of water, Ca-Mg carbonates are stable along geothermal gradients typical of subducting slabs.
Ca-Mg carbonates ( similar to Mg0.1-0.9Ca0.9-0.1CO3) partially melt beneath mid-ocean ridges and in plume settings. Ca-Mg carbonates melt incongruently, forming periclase crystals and carbonate melt between 4 and 9 GPa.
Furthermore, we show that the rare earth element (REE) signature of Group-I kimberlites, namely strong REE fractionation and depletion of heavy REE relative to the primitive mantle, is resembled by carbonate melt in equilibrium with Ca-bearing magnesite and periclase at 6 and 9 GPa.
This suggests that the dolomite-magnesite join of the CaCO3-MgCO3 system might be useful to approximate the REE signature of carbonate-rich melts parental to kimberlites.
A metagenome-assembled genome (MAG), named Methanosarcina sp. strain ERenArc_MAG2, was obtained from a 3-month-old H-2/CO2 atmosphere enrichment culture, originally inoculated with 60-m deep drill core sediment collected from the tectonic Eger Rift terrestrial subsurface.
Annotation of the recovered draft genome revealed putative archaeal methanogenesis genes in the deep biosphere.
A metagenome-assembled genome (MAG), named Methanosarcina sp. strain ERenArc_MAG2, was obtained from a 3-month-old H-2/CO2 atmosphere enrichment culture, originally inoculated with 60-m deep drill core sediment collected from the tectonic Eger Rift terrestrial subsurface.
Annotation of the recovered draft genome revealed putative archaeal methanogenesis genes in the deep biosphere.