550 Geowissenschaften
Refine
Year of publication
Document Type
- Article (884)
- Doctoral Thesis (193)
- Contribution to a Periodical (32)
- Book (26)
- Working Paper (22)
- Part of Periodical (20)
- Conference Proceeding (18)
- Part of a Book (9)
- diplomthesis (8)
- Diploma Thesis (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 (772)
- Geowissenschaften / Geographie (133)
- Biodiversität und Klima Forschungszentrum (BiK-F) (61)
- Geographie (61)
- Senckenbergische Naturforschende Gesellschaft (53)
- Präsidium (44)
- Extern (28)
- Biowissenschaften (22)
- Institut für Ökologie, Evolution und Diversität (8)
- Physik (8)
Carbonate archives record a brief snapshot of the ambient Earth’s surface conditions at their deposition. However, the geologically reasonable extraction and interpretation of geochemical proxy data from ancient, diagenetically altered rock archives is fraught with problems. Three issues stand out: the dichotomy between petrographic and geochemical alteration; the lack of quantitative age constraints for specific diagenetic phases resulting in a poorly constrained admixture of local, basin-wide and over-regional (far-field) features; and an often insufficient understanding of the temperatures and compositions of diagenetic fluids. Here, the archive of Devonian marine limestones exposed to multiple far-field diagenetic events is used as an example to explore the above-listed issues. Methods applied include petrography, micro XRF, fluid inclusion data, clumped isotopes, δ13C and δ18O isotopes, 87Sr/86Sr ratios and quartz trace element data. Devonian limestones studied here were overprinted by two cross-cutting regional fault zones (T ≈ 230 °C) by multiple events between the Variscan Orogeny and the late Paleogene. The following processes are recorded: (i) protolith deposition and partial dolomitisation during rapid burial in the Middle/Late Devonian (T ≈ 180 °C); (ii) deep burial to ca 6.5 km and tectonic/hydrothermal overprint during the Variscan Orogeny in the Carboniferous (T ≈ 90–230 °C); (iii) rapid uplift to 1–2 km burial depth at the end of the Variscan Orogeny and hypogene karstification (T ≈ 50 to 100 °C) initiated by regional geology in the Permian/Triassic; (iv) tectonic/hydrothermal overprint during the opening of the Proto-Atlantic Ocean between the Early Jurassic and the Early Cretaceous (T ≈ 50 to 130 °C); (v) tectonic/hydrothermal overprint including renewed hypogene karstification and hydrothermal calcite cement precipitation (T ≈ 50 to 180 °C) during Alpine Orogeny between the Late Cretaceous and late Paleogene. Despite this complex series of diagenetic events, the protolith limestones largely preserved their respective Middle/Late Devonian dissolved inorganic carbon (DIC) and 87Sr/86Sr signatures. This study documents that geochemical proxy data, placed into their petrographic, paleotemperature, and local to over-regional context, significantly increases the ability to extract quantitative information from ancient carbonate rock archives. Research shown here has wider relevance for carbonate archive research in general.
Ecophysiological studies on Antarctic cryptophytes to assess whether climatic changes such as ocean acidification and enhanced stratification affect their growth in Antarctic coastal waters in the future are lacking so far. This is the first study that investigated the combined effects of increasing availability of pCO2 (400 and 1000 µatm) and irradiance (20, 200 and 500 μmol photons m−2 s −1) on growth, elemental composition and photophysiology of the Antarctic cryptophyte Geminigera cryophila. Under ambient pCO2, this species was characterized by a pronounced sensitivity to increasing irradiance with complete growth inhibition at the highest light intensity. Interestingly, when grown under high pCO2 this negative light effect vanished and it reached highest rates of growth and particulate organic carbon production at the highest irradiance compared to the other tested experimental conditions. Our results for G. cryophila reveal beneficial effects of ocean acidification in conjunction with enhanced irradiance on growth and photosynthesis. Hence, cryptophytes such as G. cryophila may be potential winners of climate change, potentially thriving better in more stratified and acidic coastal waters and contributing in higher abundance to future phytoplankton assemblages of coastal Antarctic waters.
Highlights
• Constrictional structures range from dome-and-basin folds to coeval folds and boudins.
• Under bulk constriction, the competent layer rotates slower than a passive plane.
• Extension-parallel and –perpendicular folds grow simultaneously.
• Extension-perpendicular folds affect previous boudins.
Abstract
We conducted scaled analogue modelling to show the influence of varying single layer initial orientation on the geometry of folds and boudins in a bulk constrictional strain field. The initial angle between the plane of shortening and the competent layer (θZ(i)) was incrementally increased from 0° to 90° by multiples of 11.25°. While the amount of layer thickening decreased with increasing θZ(i), the deformation structures produced range from pure dome-and-basin folds to coeval folds and boudins. Based on the attitude of fold axes, there are extension-parallel (FEPR) and extension-perpendicular (FEPP) folds, with axes subparallel and subperpendicular to the principal stretching axis (X), respectively. Coeval growth of FEPR folds and boudins occurred when θZ(i) > ca. 25°. The FEPP folds can be subdivided into a first type which affect the entire layer (if θZ(i) ranges between 11.25 and 78.75°) and a second type, referred to as FBEPP folds, which are affecting pre-existing boudins if θZ(i) > 45°. The interlimb angle of all types of folds increases with increasing θZ(i). Folds and boudins similar to the ones produced in this study can be found in salt domes and in tectonites of subduction zones.
Non-matrix-matched calibration of laser ablation ICPMS (trace/major) element data is a common quantification strategy. However, LA sampling is associated with downhole elemental fractionation, potentially causing inaccuracies if the magnitude of fractionation between the sample and reference material (RM) differs. Here, we estimate fractionation factors (FFs) for different elements (El) in a range of RMs relative to NIST SRM610/612 (FFEl/Ca-NIST) and evaluate element-specific corrections for downhole fractionation using these measured FFEl/Ca-NIST. Significantly different mean El/Ca values were observed before and after correction, particularly for the alkali elements (all RMs), and B, Fe, and Zn (some RMs), notably improving accuracy, especially for the alkali elements. In cases where this methodology does not result in an accuracy improvement, this may help identify underlying issues in reported/reference values for RMs, given that this phenomenon should be accounted for. Overall, we recommend considering routine assessment of FFs and applying a FF correction to enhance data quality.
"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.