Geowissenschaften
Refine
Year of publication
Document Type
- Article (593)
- Doctoral Thesis (139)
- Book (25)
- Contribution to a Periodical (21)
- Conference Proceeding (19)
- Working Paper (19)
- Part of Periodical (10)
- Diploma Thesis (8)
- diplomthesis (7)
- Preprint (6)
Is part of the Bibliography
- no (872)
Keywords
- climate change (13)
- Climate change (5)
- Geochemistry (5)
- Klimaänderung (5)
- Atmospheric chemistry (4)
- Boden (4)
- Deutschland (4)
- Klima (4)
- Modellierung (4)
- Stratosphäre (4)
Institute
- Geowissenschaften (872)
- Senckenbergische Naturforschende Gesellschaft (59)
- Präsidium (49)
- Biodiversität und Klima Forschungszentrum (BiK-F) (45)
- Geographie (26)
- Biowissenschaften (15)
- Medizin (11)
- Physik (7)
- Institut für Ökologie, Evolution und Diversität (6)
- Biochemie und Chemie (5)
Non-forest ecosystems, dominated by shrubs, grasses and herbaceous plants, provide ecosystem services including carbon sequestration and forage for grazing, yet are highly sensitive to climatic changes. Yet these ecosystems are poorly represented in remotely-sensed biomass products and are undersampled by in-situ monitoring. Current global change threats emphasise the need for new tools to capture biomass change in non-forest ecosystems at appropriate scales. Here we assess whether canopy height inferred from drone photogrammetry allows the estimation of aboveground biomass (AGB) across low-stature plant species sampled through a global site network. We found mean canopy height is strongly predictive of AGB across species, demonstrating standardised photogrammetric approaches are generalisable across growth forms and environmental settings. Biomass per-unit-of-height was similar within, but different among, plant functional types. We find drone-based photogrammetry allows for monitoring of AGB across large spatial extents and can advance understanding of understudied and vulnerable non-forested ecosystems across the globe.
Highlights
• Northern and eastern grassland-savanna boundary defined by minimum temperature.
• Dynamics of fire, frost and growing season temperatures combine to produce this limit.
• Western limit is related to moisture availability.
• Modern, high-resolution climate data enables refinement of bioclimatic limits.
• Reparameterisation improves global model performance at regional scale.
Abstract
Understanding the controls of biome distributions is crucial for assessing terrestrial ecosystem functioning and its response to climate change. We analysed to what extent differences in climate factors (minimum temperatures, water availability, and growing season temperatures (degree days above 5 °C (GDD5)) might explain the poorly understood borders between grasslands, savannas and shrublands in eastern South Africa. The results were used to improve bioclimatic limits in the dynamic global vegetation model (DGVM) LPJ-GUESS. The vegetation model was also used to explore the role of fire in the biome borders. Results show no clear differences between the adjacent biomes in water availability. Treeless grasslands primarily occur in areas with minimum temperatures and GDD5 values below that of savannas. The standard fire module in LPJ-GUESS is not able to reproduce observed burned area patterns in the study region, but simulations with prescribed fire return intervals show that a combination of low temperatures and fire can explain the treeless state of the grassland biome. These results confirm earlier hypotheses that a combination of low winter temperatures, causing frost damage to trees, and low growing season temperatures that impede tree sapling growth and recruitment, particularly under re-occurring fires, drive the grassland-savanna border. With these insights implemented, the LPJ-GUESS simulation results substantially improved grass distribution in the grassland biome, but challenges remain concerning the grassland-shrubland boundary, tree-grass competition and prognostic fire modelling.
Measurements of vertical velocity from vertically pointing Doppler lidars are used to derive the profiles of vertical velocity variance. Observations were taken during the FESSTVaL (Field Experiment on Submesoscale Spatio-Temporal Variability in Lindenberg) campaign during the warm seasons of 2020 and 2021. Normalized by the square of convective velocity scale, the average vertical velocity variance profile follows the universal profile of Lenschow et al. (1980), however, daily profiles still show a high day-to-day variability. We found that moisture transport and the content of moisture in the boundary layer could explain the remaining variability of the normalized vertical velocity variance. The magnitude of the normalized vertical velocity variance is highest on clear-sky days, and decreases as the relative humidity increase and surface latent heat flux decrease in cloud-topped and rainy days. This suggests that moisture content and moisture transport are limiting factors for the intensity of turbulence in the convective boundary layer. We also found that the intensity of turbulence decreases with an increase in boundary layer cloud fraction during FESSTVaL, while the latent heating in the cloud layer was not a relevant source of turbulence in this case. We conclude that a new vertical velocity scale has to be defined that would take into account the moist processes in the convective boundary layer.
Measurements of vertical velocity from vertically pointing Doppler lidars are used to derive the profiles of normalized vertical velocity variance. Observations were taken during the FESSTVaL (Field Experiment on Submesoscale Spatio-Temporal Variability in Lindenberg) campaign during the warm seasons of 2020 and 2021. Normalized by the square of the convective velocity scale, the average vertical velocity variance profile follows the universal profile of Lenschow et al. (1980). However, daily profiles still show a high day-to-day variability. We found that moisture transport and the content of moisture in the boundary layer could explain the remaining variability of the normalized vertical velocity variance. The magnitude of the normalized vertical velocity variance is highest on clear-sky days and decreases as the absolute humidity increases and surface latent heat flux decreases on cloud-topped days. This suggests that moisture content and moisture transport are limiting factors for the intensity of turbulence in the convective boundary layer. We also found that the intensity of turbulence decreases with an increase in the boundary layer cloud fraction during FESSTVaL, while the latent heating in the cloud layer was not a relevant source of turbulence in this case. We conclude that a new vertical velocity scale has to be defined that would take into account the moist processes in the convective boundary layer.
A comprehensive study of sillenite Bi12SiO20 single-crystal properties, including elastic stiffness and piezoelectric coefficients, dielectric permittivity, thermal expansion and molar heat capacity, is presented. Brillouin-interferometry measurements (up to 27 GPa), which were performed at high pressures for the first time, and ab initio calculations based on density functional theory (up to 50 GPa) show the stability of the sillenite structure in the investigated pressure range, in agreement with previous studies. Elastic stiffness coefficients c11 and c12 are found to increase continuously with pressure while c44 increases slightly for lower pressures and remains nearly constant above 15 GPa. Heat-capacity measurements were performed with a quasi-adiabatic calorimeter employing the relaxation method between 2 K and 395 K. No phase transition could be observed in this temperature interval. Standard molar entropy, enthalpy change and Debye temperature are extracted from the data. The results are found to be roughly half of the previous values reported in the literature. The discrepancy is attributed to the overestimation of the Debye temperature which was extracted from high-temperature data. Additionally, Debye temperatures obtained from mean sound velocities derived by Voigt-Reuss averaging are in agreement with our heat-capacity results. Finally, a complete set of electromechanical coefficients was deduced from the application of resonant ultrasound spectroscopy between 103 K and 733 K. No discontinuities in the temperature dependence of the coefficients are observed. High-temperature (up to 1100 K) resonant ultrasound spectra recorded for Bi12MO20 crystals revealed strong and reversible acoustic dissipation effects at 870 K, 960 K and 550 K for M = Si, Ge and Ti, respectively. Resonances with small contributions from the elastic shear stiffness c44 and the piezoelectric stress coefficient e123 are almost unaffected by this dissipation.
New U–Pb ages of detrital and igneous zircons of the Uppermost Unit of Crete shed light on its provenance and on Eohellenic to Eoalpine imprints in the eastern Mediterranean. The detrital zircons of all nappes show Variscan ages and are characterized by a Minoan-type age spectrum, which is typical for the NE margin of Gondwana. Parts of the metasedimentary rocks are unexpectedly young. Their detrital zircon ages continue via the Permian until the Late Triassic, Middle Jurassic and Early Cretaceous. The high-grade metamorphic rocks of the Asterousia crystalline complex are likely equivalents of the low-grade metamorphic trench and fore-arc deposits of the Vatos nappe pointing to Late Cretaceous slab roll back. The presence of both late Permian detrital zircons and Late Cretaceous arc-type granitoids suggest that the Uppermost Unit of Crete is derived from the late Permian/Late Cretaceous magmatic belt situated north of the Sava–Vardar–Izmir–Ankara Suture in the Strandja–Rhodope area. To achieve their recent position on Crete, the nappes had to travel more than 500 km. The traveling path is well tracked by rocks of the Upper Cycladic Unit, which are similar to those of the Uppermost Unit of Crete. The large displacement of the Cretan nappes was controlled not only by nappe transport, but probably also by dextral strike–slip along the North Anatolian Fault Zone and related counterclockwise rotation of the Anatolian block since the Eocene.
Highlights
• Solidification and cooling of an intruded dyke or sill within the middle or shallow crust generate stresses of order 200 MPa, which relax on time scales of a few years to million years.
• Stresses may exceed the brittle strength forming tensile fractures.
• Combined with the pressure gradient within the over-pressurized felsic melts, this explains the migration of felsic contact melt into shrinkage cracks (Sederholm-type veins).
Abstract
Rapid emplacement of a mafic dyke or sill at mid-crustal depth heats and possibly melts the felsic wall rock followed by solidification. Associated volume changes generate stresses, possibly enforcing brittle failure and melt migration. We model the evolution of melting, solidification, temperature, and stress including visco-elastic relaxation in 1D - dykes or -sills using realistic rock rheologies of the Weschnitz pluton (Odenwald). For deep emplacement (Case 1, 15.3 km) extensive contact melting of the wall rock occurs, for shallow emplacement (Case 2, 10 km) it is negligible. The stresses are zero at high melt fractions, but increase during solidification and cooling: The intrusion orthogonal stress is always zero. The intrusion parallel stress σ‖ within the intrusion is tensile (O(200 MPa)). It relaxes on a time scale between a few years (Case 1) and 0.6 m.y. (Case 2). Within the wall rock σ‖ is compressive during heating, but becomes tensile under solidification and cooling. Wall rock stresses relax on a time scale of months to 100 years. A Deborah number is defined based on viscous to thermal relaxation allowing generalization of our results. Adding lithostatic stresses, the total stresses of Case 1 remain below the brittle strength, while for Case 2 they may exceed it. Adding the lithostatic pressure to the melt pressure, the effective stresses exceed the brittle strength and intrusion orthogonal tensile fractures are predicted. Combined with the pressure gradient within the over-pressurized felsic melts generated in the wall rock, this explains the migration of felsic contact melt into shrinkage cracks of the mafic sill in the Weschnitz pluton.
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.
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.
Rationale: Potassium (K) is a major component of several silicate minerals and seawater, and, therefore, constraining past changes in the potassium cycle is a promising way of tracing large-scale geological processes on Earth. However, [K] measurement using inductively coupled plasma mass spectrometry (ICP-MS) is challenging due to an ArH+ interference, which may be of a similar magnitude to the K+ ion beam in samples with <0.1% m/m [K].
Methods: In this work, we investigated the effect of the ArH+ interference on K/Ca data quality by comparing results from laser-ablation (LA)-ICP-MS measured in medium and high mass resolution modes and validating our LA results via solution ICP-optical emission spectroscopy (OES) and solution ICP-MS measurements. To do so, we used a wide range of geological reference materials, with a particular focus on marine carbonates, which are potential archives of past changes in the K cycle but are typically characterised by [K] < 200 μg/g. In addition, we examine the degree to which trace-element data quality is driven by downhole fractionation during LA-ICP-MS measurements.
Results: Our results show that medium mass resolution (MR) mode is sufficiently capable of minimising the effect of the ArH+ interference on K+. However, the rate of downhole fractionation for Na and K varies between different samples as a result of their differing bulk composition, resulting in matrix-specific inaccuracy. We show how this can be accounted for via downhole fractionation corrections, resulting in an accuracy of better than 1% and a long-term reproducibility (intermediate precision) of <6% (relative standard deviation) in JCp-1NP using LA-ICP-MS in MR mode.
Conclusion: Our [K] measurement protocol is demonstrably precise and accurate and applicable to a wide range of materials. The measurement of K/Ca in relatively low-[K] marine carbonates is presented here as a key example of a new application opened up by these advances.
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.
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.
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.