Refine
Year of publication
Document Type
- Article (593) (remove)
Has Fulltext
- yes (593)
Is part of the Bibliography
- no (593)
Keywords
- climate change (11)
- Climate change (5)
- Atmospheric chemistry (4)
- Geochemistry (4)
- loess (4)
- Biogeochemistry (3)
- COSMO-CLM (3)
- Klimawandel (3)
- Palaeoclimate (3)
- Salinity (3)
Institute
- Geowissenschaften (593) (remove)
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.
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.
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.