550 Geowissenschaften
Refine
Year of publication
Document Type
- Article (889)
- Doctoral Thesis (194)
- Contribution to a Periodical (32)
- Book (26)
- Working Paper (22)
- Part of Periodical (20)
- Conference Proceeding (18)
- Part of a Book (9)
- Diploma Thesis (8)
- diplomthesis (7)
- Other (5)
- Report (5)
- Master's Thesis (4)
- Review (4)
- Periodical (3)
- Bachelor Thesis (2)
- Habilitation (2)
- Preprint (2)
Language
Keywords
- climate change (11)
- Climate change (7)
- Klima (7)
- Klimaänderung (7)
- Modellierung (7)
- COSMO-CLM (6)
- Klimawandel (6)
- precipitation (6)
- Atmospheric chemistry (5)
- Deutschland (5)
Institute
- Geowissenschaften (776)
- Geowissenschaften / Geographie (135)
- Biodiversität und Klima Forschungszentrum (BiK-F) (61)
- Geographie (61)
- Senckenbergische Naturforschende Gesellschaft (54)
- Präsidium (44)
- Extern (28)
- Biowissenschaften (22)
- Institut für Ökologie, Evolution und Diversität (8)
- Physik (8)
"High-aluminous coal" is an important coal kind and widely distributed in North China in age of Permo-Carboniferous period. To explore their occurrence state, a total of 15 harmful elements (Li, Ga, In, Cd, Cr, Pb, Be, Mn, Zn, Ag, Co, Ni, Cu, Ba and U) in the No.9 coal and No.11 coal collected from Pingshuo mining district were determined by inductively coupled plasma mass spectrometry (ICP-MS) and scanning electron microscope with energy spectrum (SEM-EDX). The results showed that the content of Li, Ga, In, Pb, Ag and U were all exceed the world hard coal. In view of the result of clustering analysis within trace elements, it was found that Co, Ni, Zn, Cu, Ag and Cr were mainly associated with sulfide minerals due to their common sulfophilic property. Manganese was mainly occurred in carbonate minerals, while Ba, Cd and U were mainly associated with total minerals. In addition, Pb was related to sulfides and Be is mainly distributed in clay minerals. The enrichment of such harmful elements in Pingshuo coal was caused by the combined effect of transgression and input of terrestrial materials in the peat accumulation stage. Li, Ga, In and Ag have reached the harmful grade.
Partial melting of crustal and mantle rocks under pressure from impedance spectroscopy measurements
(2004)
The purpose of this work is to achieve a better understanding of the physical properties of rocks during partial melting processes. The electrical conductivity of some crustal and upper mantle rocks was measured prior and above the melting under pressure. The variations of the electrical conductivity were compared with the distribution of melt in partially molten rock samples. The electrical conductivity was estimated from the impedance spectroscopy at temperatures between 800 and 1450˚C and at pressures between 0.3 and 2 GPa. These measurements were performed in a piston cylinder apparatus. At temperatures above the melting, samples were equilibrated during a long time and subsequently quenched. Thin sections were prepared and topology, volume fraction and chemical composition of melt was analyzed by using a microprobe. Above the solidus temperature, the electrical conductivity increases for about 1 to 2 orders of magnitude in comparison with non-melted rocks. The "melt effect" seems to reflect the formation of an interconnected network of melt. When a complete melt connectivity is established, the charge transport follows the network of the formed melt films at grain boundaries. Usually, it takes a long time in order to reach a steady state of the electrical resistance in partially molten rocks. Only when a steady state of the electrical resistance is achieved, the bulk conductivity of a sample can be measured properly. The time-independent electrical conductivity were found only after 200 h of annealing time at a desired temperature.
Usually, the measurements of a dihedral angle on grain-liquid interfaces in rocks show that the wetting of grain faces start to develop at temperatures slightly above the solidus temperature. The development of these faces should lead to a continuous melt network even at small melt fractions of few wt.%. This result is not confirmed by our electrical conductivity measurements. The complete interconnection of the melt phase, which was mark by an increase of the electrical conductivity, corresponds to a temperature significantly above the solidus temperature, for at least 30-50˚C. The development of these faces stimulate a significant increase of the electrical conductivity, and corresponds to the occurence of at least 5 wt.% of a melt fraction. This result could be explained by deviations from the textural equilibrium of a melt phase topology in partially molten samples due to heterogeneous grain size distribution, misorientation of grains and anisotropy of the superficial energy of adjacent grain boundaries.
Some mixing models that allow to calculate the electrical conductivity of a composite as a function of a melt fraction were examined and the results of these calculations are discussed.
The experimental results were compared to the conductivity data obtained from magnetotelluric and electromagnetic measurements in the Northern part of mid-Atlantic ridge where a series of magma chambers are presumably located. There is a good agreement between our conductivity values for a melt fraction of 10-13 the conductivity estimated in the Reykjanes ridge zone.
Tectonics and geothermal gradients from subduction to collision in the NW Variscan Iberian Massif
(2022)
The earliest tectonometamorphic record of tectonic slices incorporated to the base of an orogen holds the key to understand how an orogen is built. The tectonic pile of the NW Iberian section of the Variscan Orogen includes tectonic slices separated by crustal-scale thrusts. The earliest tectonometamorphic record in the uppermost parautochthon is calculated at 11–14 kbar and 450–500°C (P-T gradient about 13°C/km), suggesting a subduction-related metamorphic recrystallization at lower pressure than the overlying Lower Allochthon. Early conditions calculated in the autochthon (9–10 kbar and 425–450°C; 16°C/km) point to a relatively ‘cold’ collisional setting. Higher thermal gradients obtained from some sections of the autochthon (11–12 kbar and 700–725°C; 21°C/km) and the Lower Parautochthon (7.5 kbar and 550–700°C; 24–31°C/km), correspond to more advanced and ‘hot’ stages of collision. New U–Pb monazite geochronology indicates a 318–311 Ma age for the final formation of HT domes in the region. We propose the rapid decrease in P-T gradient (from <10 to 16°C/km) documents a fail to sustain further burial along a regular subduction zone. We consider the subsequent increase in the geothermal gradient (from 16 to 31°C/km) as the culmination of previous crustal accretion and the onset of crustal underthrusting and later processes in a collisional stage. We propose these switches in the early tectonometamorphic record of individual tectonic slices as potential markers to track the transition from subduction to collision in collisional orogens.
While high-quality climate reconstructions of some past warm periods in the Cenozoic era now exist, the geological processes responsible for driving the observed longterm changes in atmospheric CO2 are not sufficiently well understood. The long-term change in atmospheric CO2 across the Cenozoic has been proposed to be driven by processes such as terrestrial weathering, organic carbon production and burial, reverse weathering, and volcanic degassing. One way of constraining the relative importance of the various driving forces proposed so far is to better understand the degree to which ocean chemistry has changed because the chemistry of seawater responds to geologic processes that drive atmospheric CO2. In addition, knowledge of the concentration of the major elements in seawater is crucial for accurately applying proxies such as those based on the boron isotopic composition and Mg/Ca of marine carbonates (a proxy for palaeo pH/CO2 and palaeotemperature, respectively). Previously reported records of seawater composition are primarily derived from fluid inclusions in marine evaporites; however, the results are sparse due to the limited availability of such deposits. In this thesis, changes in the Eocene seawater chemistry were reconstructed using trace element (elements/Ca) and isotopic (δ26Mg) proxies in a Larger Benthic Foraminifera (LBFs), i.e., Nummulites sp., to constrain the driving processes of long-term changes in seawater chemistry.
To achieve the objective of this thesis, first, a measurement protocol was established using LA-ICPMS to measure the K/Ca ratio simultaneously with other element/calcium ratios, which is challenging due to the interference of ArH+ on K+. Utilising this newly established measurement protocol, laboratory-cultured Operculina ammonoides grown at different seawater calcium concentrations ([Ca2+]), repeated at different temperatures, as well as modern O. ammonoides collected from different regions exhibiting a range of seawater parameters, were investigated. A significant correlation was observed between K/Casw and K/CaLBF, allowing K/CaLBF to potentially be used as a proxy for seawater major ion reconstructions. In addition, modern O. ammonoides demonstrated no significant influence of most seawater parameters (temperature, salinity, pH, or [CO32-]) on K/CaLBF. Modern O.
ammonoides were also assessed for their Mg isotopic composition (δ26Mg), revealing no significant effect of temperature or salinity on δ26MgLBF. Furthermore, the Mg isotopic fractionation in O. ammonoides was found to be close to that of inorganic calcite, indicating minimal vital effects in these large benthic foraminifera.
Operculina ammonoides is the nearest living relative of the abundant Eocene genus Nummulites, enabling the reconstruction of seawater chemistry using the calibration based on O. ammonoides. The trace elemental/calcium proxies, including Na/Ca, K/Ca, and Mg/Ca, as well as the δ26Mg proxy, were investigated in Eocene Nummulites. The result showed that during the Eocene, [Ca2+]sw was 1.6-2 times higher, while [K+]sw was ~2 times lower than the modern seawater composition. Furthermore, [Mg2+]sw decreased from the early Eocene (54.3− +9 7..69 mmol kg-1 at ~55 Ma) to Late Eocene (37.8− +4 4..3 4 mmol kg-1 at ~31 Ma), followed by
an increase toward modern seawater [Mg]. In contrast, the variability in δ26Mgsw values remained within a narrow range of ~0.3 ‰ throughout the Cenozoic. The reconstructed [Ca2+]sw agrees with the suggestion that Cenozoic seawater chemistry changes can be explained via a change in the seafloor spreading rate. When combined with existing records, the observed minimal change in δ26Mgsw with an increase in [Mg2+]sw suggests an additional possible role of a decrease in the formation of authigenic clay minerals coincident with the Cenozoic decline in deep ocean temperature, which is also supported by the increase in the [K+]sw reconstructed here for the first time. This finding highlights that the reduction in seafloor-spreading rate and decline in reverse weathering during the Cenozoic era has played a significant role in the evolution of seawater chemistry, emphasizing the importance of these processes in driving long-term changes in the carbon cycle.
Highlights
• We find DBrfluid/melt = 1.19 to 3.92 for experimental Br degassing from basalt magma into aqueous fluids.
• D <1 under almost dry conditions propose only little Br degassing for dry intra-plate volcanism relative to volcanic arcs.
• An annual global Br flux of 23.5–72.9 × 109 g/y into the atmosphere was calculated.
Abstract
We present the first in-situ partitioning data for bromine between a natural basaltic melt and a coexisting fluid. For this study hydrothermal diamond anvil cell experiments at pressures up to 1.7 GPa were conducted. We combined laser heating to melt the basalt glass with external heating to lower the temperature gradient in the cell and to initiate circulation for the aqueous fluid. Bromine concentrations were measured in-situ with X-ray fluorescence in the basaltic melts, glasses, and in the fluid. From the results we calculated partition coefficients of DBrfluid/melt = 1.19 to 3.92 in the range of 0.4 to 1 GPa for aqueous fluids. Experiments with neon as the surrounding fluid (DBrfluid/melt = 0.38 ± 0.01 at 1.1 GPa) suggest that Br-release from a basalt into volatiles that have no bonding affinity with Br is weak. This should be the case for dry intra-plate volcanic eruptions. From the experimentally gained partition coefficients and from global Br concentration values in melt inclusions of arc magmas, we calculated an annual global Br flux of 23.5–72.9 × 109 g/y.
We present a new experimental dataset on the impact of the heavy halogens chlorine, bromine and iodine on the Raman water bands concerning pressure and their concentration at room temperature. These experiments were conducted at ambient temperature, with variations in halogen concentration and pressure ranging from 0 to 1.4 GPa.
The strength of the Raman water band shift change increases with the ionic size from chlorine, over bromine, to iodine. Our experiments further demonstrate that increased pressure diminishes the impact of the halogen shift change to a varying extent for each of the three halogens. This finding can have significant implications for the salinity calculation of fluid inclusions in minerals such as quartz or olivine. Particularly in the low salinity range, the concentration can be markedly underestimated if the pressure effect is neglected. For experiments in diamond anvil cells involving halogens dissolved in water, the change in Raman water band shifts can serve either as a new tool to monitor pressure, or to monitor the salinity.
Upper mantle shear zones are complex systems where deformation is commonly closely interacting with metamorphic (solid-solid) and/or melt/fluid-rock reactions. Here, feedback processes between deformation, reactions, grain size reduction and phase mixing result in strain weakening and the localization of deformation. The expression of these interlinked processes is portrayed by the microfabrics of strained peridotites and pyroxenites. The present thesis is focusing on these processes and their impact on the deformation in three upper mantle shear zones situated in the peridotite massifs of Lanzo (Italian Alps), Erro-Tobbio (Italian Alps) and Ronda (Betic Cordillera, Spain). In all three shear zones, the presence of melt led to phase mixing either by interstitial crystallization of pyroxenes from a Si-saturated and partially also highly evolved melt or by melt-rock reactions of pyroxene porphyroclasts with a Si-undersaturated melt. The effect of melt on the localization of strain is twofold and variable. Enhanced deformation by melt-wetted boundaries is assumed for all shear zones. Additionally, phase mixing by crystallization of interstitial pyroxenes or melt-rock reactions reduce or maintain the grain size by the formation of fine grained neoblasts and secondary phase boundary pinning. In this regard, pre- to early syn-kinematic, map-scale percolation of OH-bearing, evolved melts in the NW Ronda peridotite massif and the associated crystallization of interstitial pyroxenes result in the activation of grain size sensitive deformation mechanisms in the entire melt-effected area. In the rocks collected at Erro-Tobbio, syn-kinematic melt-rock reactions of pyroxene porphyroclasts and Si-undersaturated melt led to the formation of ultramylonitic neoblast tails (grain size ~10 μm). Compared to the adjacent coarser-grained olivine-dominated matrix, the activation of diffusion creep led to an increase in the strain rate by an order of magnitude within interconnected ultramylonitic layers. Strain localization and softening in ultramylonitic layers are also documented in the Lanzo samples. Neoblast tails of pyroxene porphyroclasts were likewise identified as their precursor. The phase assemblage of the tails, including ortho- and clinopyroxene, olivine, plagioclase, and spinel (± amphibole), and their geochemical trends suggest, unlike in Erro-Tobbio, a formation by continuous net-transfer reactions enhanced by the spinel lherzolite to plagioclase lherzolite transition.
The new results obtained from the three studied shear zones underscore the importance of reactions for the interlinked processes of grain size reduction, phase mixing, strain localization and strain softening in upper mantle shear zones. Concerning strain localization, the nature of the reaction (solid-solid, melt/fluid-rock) seems to play a subordinate role compared to its timing. Pre- to early syn-kinematic melt-triggered reactions result in strain localization along map-scale shear zones. Late stage syn-kinematic melt-rock or metamorphic reactions under high stress conditions are capable of localizing the deformation along discrete, sub-centimeter thick ultramylonites.
Highlights
• We present the first results of a deep learning model based on a convolutional neural network for earthquake magnitude estimation, using HR-GNSS displacement time series.
• The influence of different dataset configurations, such as station numbers, epicentral distances, signal duration, and earthquake size, were analyzed to figure out how the model can be adapted to various scenarios.
• The model was tested using real data from different regions and magnitudes, resulting in the best cases with 0.09 ≤ RMS ≤ 0.33.
Abstract
High-rate Global Navigation Satellite System (HR-GNSS) data can be highly useful for earthquake analysis as it provides continuous high-frequency measurements of ground motion. This data can be used to analyze diverse parameters related to the seismic source and to assess the potential of an earthquake to prompt strong motions at certain distances and even generate tsunamis. In this work, we present the first results of a deep learning model based on a convolutional neural network for earthquake magnitude estimation, using HR-GNSS displacement time series. The influence of different dataset configurations, such as station numbers, epicentral distances, signal duration, and earthquake size, were analyzed to figure out how the model can be adapted to various scenarios. We explored the potential of the model for global application and compared its performance using both synthetic and real data from different seismogenic regions. The performance of our model at this stage was satisfactory in estimating earthquake magnitude from synthetic data with 0.07 ≤ RMS ≤ 0.11. Comparable results were observed in tests using synthetic data from a different region than the training data, with RMS ≤ 0.15. Furthermore, the model was tested using real data from different regions and magnitudes, resulting in the best cases with 0.09 ≤ RMS ≤ 0.33, provided that the data from a particular group of stations had similar epicentral distance constraints to those used during the model training. The robustness of the DL model can be improved to work independently from the window size of the time series and the number of stations, enabling faster estimation by the model using only near-field data. Overall, this study provides insights for the development of future DL approaches for earthquake magnitude estimation with HR-GNSS data, emphasizing the importance of proper handling and careful data selection for further model improvements.
PolarCAP – A deep learning approach for first motion polarity classification of earthquake waveforms
(2022)
Highlights
• We present PolarCAP, a deep learning model that can classify the polarity of a waveform with a 98% accuracy.
• The first-motion polarity of seismograms is a useful parameter, but its manual determination can be laborious and imprecise.
• We demonstrate that in several cases the model can assign trace polar-ity more accurately than a human analyst.
Abstract
The polarity of first P-wave arrivals plays a significant role in the effective determination of focal mechanisms specially for smaller earthquakes. Manual estimation of polarities is not only time-consuming but also prone to human errors. This warrants a need for an automated algorithm for first motion polarity determination. We present a deep learning model - PolarCAP that uses an autoencoder architecture to identify first-motion polarities of earth-quake waveforms. PolarCAP is trained in a supervised fashion using more than 130,000 labelled traces from the Italian seismic dataset (INSTANCE) and is cross-validated on 22,000 traces to choose the most optimal set of hyperparameters. We obtain an accuracy of 0.98 on a completely unseen test dataset of almost 33,000 traces. Furthermore, we check the model generalizability by testing it on the datasets provided by previous works and show that our model achieves a higher recall on both positive and negative polarities.
Highlights
• Germany plans more long-distances water transfers to secure drinking water supply.
• Long-distance water transfers can unfold lock-ins that limit adaptive water governance.
• Our interdisciplinary case study shows how lock-ins emerge over different spaces and times.
• Commercialisation of water but also local protests contributed to various lock-ins.
• We therefore call for context-specific assessments of potentials and risks of LDWT.
Abstract
Germany plans to expand water transfers over long distances in the light of numerous and pressing challenges for drinking water supply. Research on inter- and intrabasin water transfers warns, however, that major investments in large-scale infrastructure systems accompanied by institutional logics and political interests often lead to a so-called lock-in. As a consequence, long-distance water transfers can limit the potential for adaptive water governance in the involved supply areas over decades with negative impacts for people and the environment. By using a case study in Germany as an example, we researched when, where and how such lock-ins around long-distance water transfers emerge. In the infrastructural development of the Elbaue-Ostharz transfer system we found various lock-ins that overlap in space and time. Some are located at the centre others at the margins of the infrastructure and commercialization of the water sector as well as hydraulic and hygienic concerns interlock with local protests in a way that the expansion of the long-distance water transfer infrastructure is presented continuously as imperative. Our findings contribute to a relational understanding of lock-ins of long-distance water transfers as contingent and diverse processes. Given the widespread occurrence of lock-ins, we argue for a context-specific assessment of potentials and risks of long-distance water transfers in times of multiple crises.
Highlights
• New fumarole and thermal water data for Askja and Kverkfjöll volcanoes, Iceland.
• Data compared to modelled compositions and fluxes of magmatic gas.
• Fumarole compositions compatible with origin of CO2 and S from degassing intrusions.
• Intrusive magmatic fluxes sufficient to sustain hydrothermal fluxes of CO2 and S in Iceland
• Magma degassing insignificant/minor source of H2O and Cl to Icelandic hydrothermal fluids
Abstract
Mantle volatiles are transported to Earth's crust and surface by basaltic volcanism. During subaerial eruptions, vast amounts of carbon, sulfur and halogens can be released to the atmosphere during a short time-interval, with impacts ranging in scale from the local environment to the global climate. By contrast, passive volatile release at the surface originating from magmatic intrusions is characterized by much lower flux, yet may outsize eruptive volatile quantities over long timescales. Volcanic hydrothermal systems (VHSs) act as conduits for such volatile release from degassing intrusions and can be used to gauge the contribution of intrusive magmatism to global volatile cycles. Here, we present new compositional and isotopic (δD and δ18O-H2O, 3He/4He, δ13C-CO2, Δ33S-δ34S-H2S and SO4) data for thermal waters and fumarole gases from the Askja and Kverkfjöll volcanoes in central Iceland. We use the data together with magma degassing modelling and mass balance calculations to constrain the sources of volatiles in VHSs and to assess the role of intrusive magmatism to the volcanic volatile emission budgets in Iceland.
The CO2/ΣS (10−30), 3He/4He (8.3–10.5 RA; 3He/4He relative to air), δ13C-CO2 (−4.1 to −0.2 ‰) and Δ33S-δ34S-H2S (−0.031 to 0.003 ‰ and −1.5 to +3.6‰) values in high-gas flux fumaroles (CO2 > 10 mmol/mol) are consistent with an intrusive magmatic origin for CO2 and S at Askja and Kverkfjöll. We demonstrate that deep (0.5–5 kbar, equivalent to ∼2–18 km crustal depth) decompression degassing of basaltic intrusions in Iceland results in CO2 and S fluxes of 330–5060 and 6–210 kt/yr, respectively, which is sufficient to account for the estimated CO2 flux of Icelandic VHSs (3365–6730 kt/yr), but not the VHS S flux (220–440 kt/yr). Secondary, crystallization-driven degassing from maturing intrusions and leaching of crustal rocks are suggested as additional sources of S. Only a minor proportion of the mantle flux of Cl is channeled via VHSs whereas the H2O flux remains poorly constrained, because magmatic signals in Icelandic VHSs are masked by a dominant shallow groundwater component of meteoric water origin. These results suggest that the bulk of the mantle CO2 and S flux to the atmosphere in Iceland is supplied by intrusive, not eruptive magmatism, and is largely vented via hydrothermal fields.
Highlights
• Subcrustal earthquakes detected beneath Fogo volcano, Cape Verde.
• At the focal depth of 40 km temperatures are likely too high for brittle failure.
• The earthquakes may originate from magma injection into a deep subcrustal reservoir.
• This observation indicates a distinct magma supply system of Fogo volcano.
Abstract
Fogo volcano belongs to the Cape Verde hotspot and its most recent eruption occurred from November 2014 to February 2015. From January to December 2016 we operated a temporary seismic network and array on Fogo and were able to locate 289 earthquakes in total. Array analysis shows that most of the events occur within the crust at distances >25 km near the neighboring island of Brava. However, on 15th August 2016 the network recorded an isolated cluster of >20 earthquakes, 13 of which could be located beneath the southern part of Fogo. The differences between S- and P-wave arrival times at steep incidence clearly indicate focal depths between approximately 38 and 44 km whereas receiver-function analyses place the Moho discontinuity at depths between 11 and 14 km. Thus, the earthquakes are located well within the upper mantle directly beneath Fogo. In view of the elevated upper-mantle temperatures within a hotspot regime, we propose that fracturing induced by magmatic injection is the most likely cause for the observed deep earthquakes.
Highlights
• Full automatized analysis of teleseismic XKS shear wave splitting.
• Rapid analysis of large seismological data sets.
• Automated window selection and quality classification.
• Application to the USArray Transportable Array including expansion to Alaska.
• Improved statistical evidence and objectivity of derived effective splitting.
Abstract
Recent technological advances have led to community wide use of large-scale seismic experiments which produce seismic data on previously impossible scales. Standard processing procedures thus require automatization to facilitate a fast and objective analysis of the data. Among these, XKS-splitting is an important tool to derive first insights into the Earth's deformation regimes at depth by studying seismic anisotropy. Most often, shear-wave splitting is interpreted to represent crystallographic preferred orientation (CPO) of mantle minerals like olivine as dominating feature and can thus be used as a proxy of mantle flow processes. Here, we introduce an addition to the MATLAB®-based SplitRacer tool box (Reiss and Rümpker 2017) which automatizes the entire XKS-splitting procedure. This is achieved by the automatization of 1) choosing a time window based on spectral analyses and 2) categorization of results based on three different XKS-splitting methods (energy minimization, rotation correlation and splitting intensity). This provides effective and objective results for splitting as well as null-measurement results. This extension allows to use SplitRacer without a graphical interface and introduces a bootstrapping statistics as error estimate of the single layer joint splitting method. The procedures are designed to allow a fast and more objective analysis of a vast amount of data, as produced by recent seismic deployments (e.g. USArray, AlpArray). We test this automatization by applying the analysis to the USArray data set, which has approximately 1900 stations with between two to fifteen years of data. We can reproduce the general pattern of the results from former studies with the more objective automatic analysis. Based on a joint-splitting approach, we approximate the splitting effect at individual stations by a single anisotropic layer. As we include null-measurements as well as a larger data set as previous studies, we can provide improved statistical evidence for these effective splitting parameters.
Highlights
• We show the first observations of seismo-acoustic tremor at Oldoinyo Lengai, the world's only active carbonatite volcano.
• We observe significant changes in seismic and acoustic tremor properties and their correlation in one year of data collection.
• Using satellite-based thermal data, we identify different volcanic processes (degassing, lava pond dynamics and spattering).
Abstract
We analyze volcanic tremor from Oldoinyo Lengai, Tanzania, which is currently the only active volcano on Earth producing carbonatitic lavas. Here, we use data from the recent SEISVOL deployment and focus on a co-located seismic and infrasound station about 200 m below the summit. We show the very first observations of seismo-acoustic tremor caused by carbonatitic eruptions. This seismo-acoustic tremor is highly variable throughout the ∼one year of data which we characterize by analyzing its seismic amplitude, duration, recurrence, dominant seismic frequency and harmonics. Frequency gliding occurs frequently and over short (minutes to hours) to long time scales (hours to days) and likely reflects different time-dependent mechanisms, such as evenly-spaced repeating events with a change in inter-event times, crater dynamics that alter resonators, and dike intrusions. Seismic and acoustic wavefields correlate well for stronger eruptive sequences but are only partially coherent which suggests that high-frequency seismic tremor (up to 25 Hz) may be caused by the low viscosity of the carbonatitic melt and not by ground-coupled airwaves. In addition, the comparison between seismic-acoustic and satellite InfraRed thermal data allows us to infer different volcanic activity styles which partially alternate throughout the year: intrusive activity and the construction of hornitos, degassing, activity from a lava pond, and varying styles of extrusive activity, in particular spattering. Our study provides important insights into the eruption dynamics of this peculiar volcano which suggests shallow melt storage within the crater floor.
The metasomatised continental mantle may play a key role in the generation of some ore deposits, in particular mineral systems enriched in platinum-group elements (PGE) and Au. The cratonic lithosphere is the longest-lived potential source for these elements, but the processes that facilitate their pre-concentration in the mantle and their later remobilisation to the crust are not yet well-established. Here, we report new results on the petrography, major-element, and siderophile- and chalcophile-element composition of native Ni, base metal sulphides (BMS), and spinels in a suite of well-characterised, highly metasomatised and weakly serpentinised peridotite xenoliths from the Bultfontein kimberlite in the Kaapvaal Craton, and integrate these data with published analyses. Pentlandite in polymict breccias (failed kimberlite intrusions at mantle depth) has lower trace-element contents (e.g., median total PGE 0.72 ppm) than pentlandite in phlogopite peridotites and Mica-Amphibole-Rutile-Ilmenite-Diopside (MARID) rocks (median 1.6 ppm). Spinel is an insignificant host for all elements except Zn, and BMS and native Ni account for typically <25% of the bulk-rock PGE and Au. High bulk-rock Te/S suggest a role for PGE-bearing tellurides, which, along with other compounds of metasomatic origin, may host the missing As, Ag, Cd, Sb, Te and, in part, Bi that are unaccounted for by the main assemblage.
The close spatial relationship between BMS and metasomatic minerals (e.g., phlogopite, ilmenite) indicates that the lithospheric mantle beneath Bultfontein was resulphidised by metasomatism after initial melt depletion during stabilisation of the cratonic lithosphere. Newly-formed BMS are markedly PGE-poor, as total PGE contents are <4.2 ppm in pentlandite from seven samples, compared to >26 ppm in BMS in other peridotite xenoliths from the Kaapvaal craton. This represents a strong dilution of the original PGE abundances at the mineral scale, perhaps starting from precursor PGE alloy and small volumes of residual BMS. The latter may have been the precursor to native Ni, which occurs in an unusual Ni-enriched zone in a harzburgite and displays strongly variable, but overall high PGE abundances (up to 81 ppm). In strongly metasomatised peridotites, Au is enriched relative to Pd, and was probably added along with S. A combination of net introduction of S, Au +/− PGE from the asthenosphere and intra-lithospheric redistribution, in part sourced from subducted materials, during metasomatic events may have led to sulphide precipitation at ~80–120 km beneath Bultfontein. This process locally enhanced the metallogenic fertility of this lithospheric reservoir. Further mobilisation of the metal budget stored in these S-rich domains and upwards transport into the crust may require interaction with sulphide-undersaturated melts that can dissolve sulphides along with the metals they store.
We investigate the applicability of the well-known multilevel Monte Carlo (MLMC) method to the class of density-driven flow problems, in particular the problem of salinisation of coastal aquifers. As a test case, we solve the uncertain Henry saltwater intrusion problem. Unknown porosity, permeability and recharge parameters are modelled by using random fields. The classical deterministic Henry problem is non-linear and time-dependent, and can easily take several hours of computing time. Uncertain settings require the solution of multiple realisations of the deterministic problem, and the total computational cost increases drastically. Instead of computing of hundreds random realisations, typically the mean value and the variance are computed. The standard methods such as the Monte Carlo or surrogate-based methods are a good choice, but they compute all stochastic realisations on the same, often, very fine mesh. They also do not balance the stochastic and discretisation errors. These facts motivated us to apply the MLMC method. We demonstrate that by solving the Henry problem on multi-level spatial and temporal meshes, the MLMC method reduces the overall computational and storage costs. To reduce the computing cost further, parallelization is performed in both physical and stochastic spaces. To solve each deterministic scenario, we run the parallel multigrid solver ug4 in a black-box fashion.
Die Welt im Wasserstress
(2024)
Wie haben sich die Wasserresourcen in den letzten 120 Jahren verändert? Und was passiert, wenn es bis Ende des 21. Jahrhunderts noch einmal zwei Grad wärmer wird als heute? Fragen wie diese beantwortet das globale Wasser-Modell WaterGAP, das maßgeblich vom Institut für Physische Geographie der Goethe-Universität und von der Ruhr-Universität Bochum entwickelt wird. Bislang ließen sich die damit erzeugten Daten nur von Expertinnen und Experten nutzen. Eine neue Web-App ändert das nun. Entwickelt wurde sie von dem französischen Geodaten-Unternehmen Ageoce, das dafür mit der Goethe-Universität kooperierte.
Reconstructing Oligocene-Miocene paleoelevation contributes to our understanding of the evolutionary history of the European Alps and sheds light on geodynamic and Earth’s surface processes involved in the development of Alpine topography. Despite being one of the most intensively explored mountain ranges worldwide, constraints on the elevation history of the European Alps, however, remain scarce. Here we present stable and clumped isotope geochemistry 15 measurements to provide a new paleoelevation estimate for the mid-Miocene (~14.5 Ma) European Central Alps. We apply stable isotope δ-δ paleoaltimetry on near sea level pedogenic carbonate oxygen isotope (δ18O) records from the Northern Alpine Foreland Basin (Swiss Molasse Basin) and high-Alpine phyllosilicate hydrogen isotope (δD) records from the Simplon Fault Zone (Swiss Alps). We further explore Miocene paleoclimate and paleoenvironmental conditions in the Swiss Molasse Basin through carbonate stable (δ18O, δ13C) and clumped (Δ47) isotope data from three foreland basin sections in different 20 alluvial megafan settings (proximal, mid-fan, and distal). Combined pedogenic carbonate δ18O values and Δ47 temperatures (30 ± 5°C) yield a near sea level precipitation δ18Ow value of -5.8 ± 0.2‰ and in conjunction with the high-Alpine phyllosilicate δD record suggest that the region surrounding the SFZ attained surface elevations of >4000 m no later than the mid-Miocene. Our near sea level δ18Ow estimate is supported by paleoclimate (iGCM Echam5-wiso) modeled δ18O values, which vary between -4.2 and -7.6‰ for the Northern Alpine Foreland Basin.
Stratospheric inorganic chlorine (Cly) is predominantly released from long-lived chlorinated source gases and, to a small extent, very short-lived chlorinated substances. Cly includes the reservoir species (HCl and ClONO2) and active chlorine species (i.e., ClOx). The active chlorine species drive catalytic cycles that deplete ozone in the polar winter stratosphere. This work presents calculations of inorganic chlorine (Cly) derived from chlorinated source gas measurements on board the High Altitude and Long Range Research Aircraft (HALO) during the Southern Hemisphere Transport, Dynamic and Chemistry (SouthTRAC) campaign in austral late winter and early spring 2019. Results are compared to Cly in the Northern Hemisphere derived from measurements of the POLSTRACC-GW-LCYCLE-SALSA (PGS) campaign in the Arctic winter of 2015/2016. A scaled correlation was used for PGS data, since not all source gases were measured. Using the SouthTRAC data, Cly from a scaled correlation was compared to directly determined Cly and agreed well. An air mass classification based on in situ N2O measurements allocates the measurements to the vortex, the vortex boundary region, and midlatitudes. Although the Antarctic vortex was weakened in 2019 compared to previous years, Cly reached 1687±19 ppt at 385 K; therefore, up to around 50 % of total chlorine was found in inorganic form inside the Antarctic vortex, whereas only 15 % of total chlorine was found in inorganic form in the southern midlatitudes. In contrast, only 40 % of total chlorine was found in inorganic form in the Arctic vortex during PGS, and roughly 20 % was found in inorganic form in the northern midlatitudes. Differences inside the two vortices reach as much as 540 ppt, with more Cly in the Antarctic vortex in 2019 than in the Arctic vortex in 2016 (at comparable distance to the local tropopause). To our knowledge, this is the first comparison of inorganic chlorine within the Antarctic and Arctic polar vortices. Based on the results of these two campaigns, the differences in Cly inside the two vortices are substantial and larger than the inter-annual variations previously reported for the Antarctic.
Biogenic organic precursors play an important role in atmospheric new particle formation (NPF). One of the major precursor species is α-pinene, which upon oxidation can form a suite of products covering a wide range of volatilities. Highly oxygenated organic molecules (HOMs) comprise a fraction of the oxidation products formed. While it is known that HOMs contribute to secondary organic aerosol (SOA) formation, including NPF, they have not been well studied in newly formed particles due to their very low mass concentrations. Here we present gas- and particle-phase chemical composition data from experimental studies of α-pinene oxidation, including in the presence of isoprene, at temperatures (−50 and −30 ∘C) and relative humidities (20 % and 60 %) relevant in the upper free troposphere. The measurements took place at the CERN Cosmics Leaving Outdoor Droplets (CLOUD) chamber. The particle chemical composition was analyzed by a thermal desorption differential mobility analyzer (TD-DMA) coupled to a nitrate chemical ionization–atmospheric pressure interface–time-of-flight (CI-APi-TOF) mass spectrometer. CI-APi-TOF was used for particle- and gas-phase measurements, applying the same ionization and detection scheme. Our measurements revealed the presence of C8−10 monomers and C18−20 dimers as the major compounds in the particles (diameter up to ∼ 100 nm). Particularly, for the system with isoprene added, C5 (C5H10O5−7) and C15 compounds (C15H24O5−10) were detected. This observation is consistent with the previously observed formation of such compounds in the gas phase. However, although the C5 and C15 compounds do not easily nucleate, our measurements indicate that they can still contribute to the particle growth at free tropospheric conditions. For the experiments reported here, most likely isoprene oxidation products enhance the growth of particles larger than 15 nm. Additionally, we report on the nucleation rates measured at 1.7 nm (J1.7 nm) and compared with previous studies, we found lower J1.7 nm values, very likely due to the higher α-pinene and ozone mixing ratios used in the present study.