550 Geowissenschaften
Refine
Year of publication
Document Type
- Article (889)
- Doctoral Thesis (194)
- Contribution to a Periodical (32)
- Book (26)
- Working Paper (22)
- Part of Periodical (20)
- Conference Proceeding (18)
- Part of a Book (9)
- Diploma Thesis (8)
- diplomthesis (7)
- Other (5)
- Report (5)
- Master's Thesis (4)
- Review (4)
- Periodical (3)
- Bachelor Thesis (2)
- Habilitation (2)
- Preprint (2)
Language
Keywords
- climate change (11)
- Climate change (7)
- Klima (7)
- Klimaänderung (7)
- Modellierung (7)
- COSMO-CLM (6)
- Klimawandel (6)
- precipitation (6)
- Atmospheric chemistry (5)
- Deutschland (5)
Institute
- Geowissenschaften (776)
- Geowissenschaften / Geographie (135)
- Biodiversität und Klima Forschungszentrum (BiK-F) (61)
- Geographie (61)
- Senckenbergische Naturforschende Gesellschaft (54)
- Präsidium (44)
- Extern (28)
- Biowissenschaften (22)
- Institut für Ökologie, Evolution und Diversität (8)
- Physik (8)
"High-aluminous coal" is an important coal kind and widely distributed in North China in age of Permo-Carboniferous period. To explore their occurrence state, a total of 15 harmful elements (Li, Ga, In, Cd, Cr, Pb, Be, Mn, Zn, Ag, Co, Ni, Cu, Ba and U) in the No.9 coal and No.11 coal collected from Pingshuo mining district were determined by inductively coupled plasma mass spectrometry (ICP-MS) and scanning electron microscope with energy spectrum (SEM-EDX). The results showed that the content of Li, Ga, In, Pb, Ag and U were all exceed the world hard coal. In view of the result of clustering analysis within trace elements, it was found that Co, Ni, Zn, Cu, Ag and Cr were mainly associated with sulfide minerals due to their common sulfophilic property. Manganese was mainly occurred in carbonate minerals, while Ba, Cd and U were mainly associated with total minerals. In addition, Pb was related to sulfides and Be is mainly distributed in clay minerals. The enrichment of such harmful elements in Pingshuo coal was caused by the combined effect of transgression and input of terrestrial materials in the peat accumulation stage. Li, Ga, In and Ag have reached the harmful grade.
Partial melting of crustal and mantle rocks under pressure from impedance spectroscopy measurements
(2004)
The purpose of this work is to achieve a better understanding of the physical properties of rocks during partial melting processes. The electrical conductivity of some crustal and upper mantle rocks was measured prior and above the melting under pressure. The variations of the electrical conductivity were compared with the distribution of melt in partially molten rock samples. The electrical conductivity was estimated from the impedance spectroscopy at temperatures between 800 and 1450˚C and at pressures between 0.3 and 2 GPa. These measurements were performed in a piston cylinder apparatus. At temperatures above the melting, samples were equilibrated during a long time and subsequently quenched. Thin sections were prepared and topology, volume fraction and chemical composition of melt was analyzed by using a microprobe. Above the solidus temperature, the electrical conductivity increases for about 1 to 2 orders of magnitude in comparison with non-melted rocks. The "melt effect" seems to reflect the formation of an interconnected network of melt. When a complete melt connectivity is established, the charge transport follows the network of the formed melt films at grain boundaries. Usually, it takes a long time in order to reach a steady state of the electrical resistance in partially molten rocks. Only when a steady state of the electrical resistance is achieved, the bulk conductivity of a sample can be measured properly. The time-independent electrical conductivity were found only after 200 h of annealing time at a desired temperature.
Usually, the measurements of a dihedral angle on grain-liquid interfaces in rocks show that the wetting of grain faces start to develop at temperatures slightly above the solidus temperature. The development of these faces should lead to a continuous melt network even at small melt fractions of few wt.%. This result is not confirmed by our electrical conductivity measurements. The complete interconnection of the melt phase, which was mark by an increase of the electrical conductivity, corresponds to a temperature significantly above the solidus temperature, for at least 30-50˚C. The development of these faces stimulate a significant increase of the electrical conductivity, and corresponds to the occurence of at least 5 wt.% of a melt fraction. This result could be explained by deviations from the textural equilibrium of a melt phase topology in partially molten samples due to heterogeneous grain size distribution, misorientation of grains and anisotropy of the superficial energy of adjacent grain boundaries.
Some mixing models that allow to calculate the electrical conductivity of a composite as a function of a melt fraction were examined and the results of these calculations are discussed.
The experimental results were compared to the conductivity data obtained from magnetotelluric and electromagnetic measurements in the Northern part of mid-Atlantic ridge where a series of magma chambers are presumably located. There is a good agreement between our conductivity values for a melt fraction of 10-13 the conductivity estimated in the Reykjanes ridge zone.
Tectonics and geothermal gradients from subduction to collision in the NW Variscan Iberian Massif
(2022)
The earliest tectonometamorphic record of tectonic slices incorporated to the base of an orogen holds the key to understand how an orogen is built. The tectonic pile of the NW Iberian section of the Variscan Orogen includes tectonic slices separated by crustal-scale thrusts. The earliest tectonometamorphic record in the uppermost parautochthon is calculated at 11–14 kbar and 450–500°C (P-T gradient about 13°C/km), suggesting a subduction-related metamorphic recrystallization at lower pressure than the overlying Lower Allochthon. Early conditions calculated in the autochthon (9–10 kbar and 425–450°C; 16°C/km) point to a relatively ‘cold’ collisional setting. Higher thermal gradients obtained from some sections of the autochthon (11–12 kbar and 700–725°C; 21°C/km) and the Lower Parautochthon (7.5 kbar and 550–700°C; 24–31°C/km), correspond to more advanced and ‘hot’ stages of collision. New U–Pb monazite geochronology indicates a 318–311 Ma age for the final formation of HT domes in the region. We propose the rapid decrease in P-T gradient (from <10 to 16°C/km) documents a fail to sustain further burial along a regular subduction zone. We consider the subsequent increase in the geothermal gradient (from 16 to 31°C/km) as the culmination of previous crustal accretion and the onset of crustal underthrusting and later processes in a collisional stage. We propose these switches in the early tectonometamorphic record of individual tectonic slices as potential markers to track the transition from subduction to collision in collisional orogens.
While high-quality climate reconstructions of some past warm periods in the Cenozoic era now exist, the geological processes responsible for driving the observed longterm changes in atmospheric CO2 are not sufficiently well understood. The long-term change in atmospheric CO2 across the Cenozoic has been proposed to be driven by processes such as terrestrial weathering, organic carbon production and burial, reverse weathering, and volcanic degassing. One way of constraining the relative importance of the various driving forces proposed so far is to better understand the degree to which ocean chemistry has changed because the chemistry of seawater responds to geologic processes that drive atmospheric CO2. In addition, knowledge of the concentration of the major elements in seawater is crucial for accurately applying proxies such as those based on the boron isotopic composition and Mg/Ca of marine carbonates (a proxy for palaeo pH/CO2 and palaeotemperature, respectively). Previously reported records of seawater composition are primarily derived from fluid inclusions in marine evaporites; however, the results are sparse due to the limited availability of such deposits. In this thesis, changes in the Eocene seawater chemistry were reconstructed using trace element (elements/Ca) and isotopic (δ26Mg) proxies in a Larger Benthic Foraminifera (LBFs), i.e., Nummulites sp., to constrain the driving processes of long-term changes in seawater chemistry.
To achieve the objective of this thesis, first, a measurement protocol was established using LA-ICPMS to measure the K/Ca ratio simultaneously with other element/calcium ratios, which is challenging due to the interference of ArH+ on K+. Utilising this newly established measurement protocol, laboratory-cultured Operculina ammonoides grown at different seawater calcium concentrations ([Ca2+]), repeated at different temperatures, as well as modern O. ammonoides collected from different regions exhibiting a range of seawater parameters, were investigated. A significant correlation was observed between K/Casw and K/CaLBF, allowing K/CaLBF to potentially be used as a proxy for seawater major ion reconstructions. In addition, modern O. ammonoides demonstrated no significant influence of most seawater parameters (temperature, salinity, pH, or [CO32-]) on K/CaLBF. Modern O.
ammonoides were also assessed for their Mg isotopic composition (δ26Mg), revealing no significant effect of temperature or salinity on δ26MgLBF. Furthermore, the Mg isotopic fractionation in O. ammonoides was found to be close to that of inorganic calcite, indicating minimal vital effects in these large benthic foraminifera.
Operculina ammonoides is the nearest living relative of the abundant Eocene genus Nummulites, enabling the reconstruction of seawater chemistry using the calibration based on O. ammonoides. The trace elemental/calcium proxies, including Na/Ca, K/Ca, and Mg/Ca, as well as the δ26Mg proxy, were investigated in Eocene Nummulites. The result showed that during the Eocene, [Ca2+]sw was 1.6-2 times higher, while [K+]sw was ~2 times lower than the modern seawater composition. Furthermore, [Mg2+]sw decreased from the early Eocene (54.3− +9 7..69 mmol kg-1 at ~55 Ma) to Late Eocene (37.8− +4 4..3 4 mmol kg-1 at ~31 Ma), followed by
an increase toward modern seawater [Mg]. In contrast, the variability in δ26Mgsw values remained within a narrow range of ~0.3 ‰ throughout the Cenozoic. The reconstructed [Ca2+]sw agrees with the suggestion that Cenozoic seawater chemistry changes can be explained via a change in the seafloor spreading rate. When combined with existing records, the observed minimal change in δ26Mgsw with an increase in [Mg2+]sw suggests an additional possible role of a decrease in the formation of authigenic clay minerals coincident with the Cenozoic decline in deep ocean temperature, which is also supported by the increase in the [K+]sw reconstructed here for the first time. This finding highlights that the reduction in seafloor-spreading rate and decline in reverse weathering during the Cenozoic era has played a significant role in the evolution of seawater chemistry, emphasizing the importance of these processes in driving long-term changes in the carbon cycle.
Highlights
• We find DBrfluid/melt = 1.19 to 3.92 for experimental Br degassing from basalt magma into aqueous fluids.
• D <1 under almost dry conditions propose only little Br degassing for dry intra-plate volcanism relative to volcanic arcs.
• An annual global Br flux of 23.5–72.9 × 109 g/y into the atmosphere was calculated.
Abstract
We present the first in-situ partitioning data for bromine between a natural basaltic melt and a coexisting fluid. For this study hydrothermal diamond anvil cell experiments at pressures up to 1.7 GPa were conducted. We combined laser heating to melt the basalt glass with external heating to lower the temperature gradient in the cell and to initiate circulation for the aqueous fluid. Bromine concentrations were measured in-situ with X-ray fluorescence in the basaltic melts, glasses, and in the fluid. From the results we calculated partition coefficients of DBrfluid/melt = 1.19 to 3.92 in the range of 0.4 to 1 GPa for aqueous fluids. Experiments with neon as the surrounding fluid (DBrfluid/melt = 0.38 ± 0.01 at 1.1 GPa) suggest that Br-release from a basalt into volatiles that have no bonding affinity with Br is weak. This should be the case for dry intra-plate volcanic eruptions. From the experimentally gained partition coefficients and from global Br concentration values in melt inclusions of arc magmas, we calculated an annual global Br flux of 23.5–72.9 × 109 g/y.
We present a new experimental dataset on the impact of the heavy halogens chlorine, bromine and iodine on the Raman water bands concerning pressure and their concentration at room temperature. These experiments were conducted at ambient temperature, with variations in halogen concentration and pressure ranging from 0 to 1.4 GPa.
The strength of the Raman water band shift change increases with the ionic size from chlorine, over bromine, to iodine. Our experiments further demonstrate that increased pressure diminishes the impact of the halogen shift change to a varying extent for each of the three halogens. This finding can have significant implications for the salinity calculation of fluid inclusions in minerals such as quartz or olivine. Particularly in the low salinity range, the concentration can be markedly underestimated if the pressure effect is neglected. For experiments in diamond anvil cells involving halogens dissolved in water, the change in Raman water band shifts can serve either as a new tool to monitor pressure, or to monitor the salinity.
Upper mantle shear zones are complex systems where deformation is commonly closely interacting with metamorphic (solid-solid) and/or melt/fluid-rock reactions. Here, feedback processes between deformation, reactions, grain size reduction and phase mixing result in strain weakening and the localization of deformation. The expression of these interlinked processes is portrayed by the microfabrics of strained peridotites and pyroxenites. The present thesis is focusing on these processes and their impact on the deformation in three upper mantle shear zones situated in the peridotite massifs of Lanzo (Italian Alps), Erro-Tobbio (Italian Alps) and Ronda (Betic Cordillera, Spain). In all three shear zones, the presence of melt led to phase mixing either by interstitial crystallization of pyroxenes from a Si-saturated and partially also highly evolved melt or by melt-rock reactions of pyroxene porphyroclasts with a Si-undersaturated melt. The effect of melt on the localization of strain is twofold and variable. Enhanced deformation by melt-wetted boundaries is assumed for all shear zones. Additionally, phase mixing by crystallization of interstitial pyroxenes or melt-rock reactions reduce or maintain the grain size by the formation of fine grained neoblasts and secondary phase boundary pinning. In this regard, pre- to early syn-kinematic, map-scale percolation of OH-bearing, evolved melts in the NW Ronda peridotite massif and the associated crystallization of interstitial pyroxenes result in the activation of grain size sensitive deformation mechanisms in the entire melt-effected area. In the rocks collected at Erro-Tobbio, syn-kinematic melt-rock reactions of pyroxene porphyroclasts and Si-undersaturated melt led to the formation of ultramylonitic neoblast tails (grain size ~10 μm). Compared to the adjacent coarser-grained olivine-dominated matrix, the activation of diffusion creep led to an increase in the strain rate by an order of magnitude within interconnected ultramylonitic layers. Strain localization and softening in ultramylonitic layers are also documented in the Lanzo samples. Neoblast tails of pyroxene porphyroclasts were likewise identified as their precursor. The phase assemblage of the tails, including ortho- and clinopyroxene, olivine, plagioclase, and spinel (± amphibole), and their geochemical trends suggest, unlike in Erro-Tobbio, a formation by continuous net-transfer reactions enhanced by the spinel lherzolite to plagioclase lherzolite transition.
The new results obtained from the three studied shear zones underscore the importance of reactions for the interlinked processes of grain size reduction, phase mixing, strain localization and strain softening in upper mantle shear zones. Concerning strain localization, the nature of the reaction (solid-solid, melt/fluid-rock) seems to play a subordinate role compared to its timing. Pre- to early syn-kinematic melt-triggered reactions result in strain localization along map-scale shear zones. Late stage syn-kinematic melt-rock or metamorphic reactions under high stress conditions are capable of localizing the deformation along discrete, sub-centimeter thick ultramylonites.
Highlights
• We present the first results of a deep learning model based on a convolutional neural network for earthquake magnitude estimation, using HR-GNSS displacement time series.
• The influence of different dataset configurations, such as station numbers, epicentral distances, signal duration, and earthquake size, were analyzed to figure out how the model can be adapted to various scenarios.
• The model was tested using real data from different regions and magnitudes, resulting in the best cases with 0.09 ≤ RMS ≤ 0.33.
Abstract
High-rate Global Navigation Satellite System (HR-GNSS) data can be highly useful for earthquake analysis as it provides continuous high-frequency measurements of ground motion. This data can be used to analyze diverse parameters related to the seismic source and to assess the potential of an earthquake to prompt strong motions at certain distances and even generate tsunamis. In this work, we present the first results of a deep learning model based on a convolutional neural network for earthquake magnitude estimation, using HR-GNSS displacement time series. The influence of different dataset configurations, such as station numbers, epicentral distances, signal duration, and earthquake size, were analyzed to figure out how the model can be adapted to various scenarios. We explored the potential of the model for global application and compared its performance using both synthetic and real data from different seismogenic regions. The performance of our model at this stage was satisfactory in estimating earthquake magnitude from synthetic data with 0.07 ≤ RMS ≤ 0.11. Comparable results were observed in tests using synthetic data from a different region than the training data, with RMS ≤ 0.15. Furthermore, the model was tested using real data from different regions and magnitudes, resulting in the best cases with 0.09 ≤ RMS ≤ 0.33, provided that the data from a particular group of stations had similar epicentral distance constraints to those used during the model training. The robustness of the DL model can be improved to work independently from the window size of the time series and the number of stations, enabling faster estimation by the model using only near-field data. Overall, this study provides insights for the development of future DL approaches for earthquake magnitude estimation with HR-GNSS data, emphasizing the importance of proper handling and careful data selection for further model improvements.
PolarCAP – A deep learning approach for first motion polarity classification of earthquake waveforms
(2022)
Highlights
• We present PolarCAP, a deep learning model that can classify the polarity of a waveform with a 98% accuracy.
• The first-motion polarity of seismograms is a useful parameter, but its manual determination can be laborious and imprecise.
• We demonstrate that in several cases the model can assign trace polar-ity more accurately than a human analyst.
Abstract
The polarity of first P-wave arrivals plays a significant role in the effective determination of focal mechanisms specially for smaller earthquakes. Manual estimation of polarities is not only time-consuming but also prone to human errors. This warrants a need for an automated algorithm for first motion polarity determination. We present a deep learning model - PolarCAP that uses an autoencoder architecture to identify first-motion polarities of earth-quake waveforms. PolarCAP is trained in a supervised fashion using more than 130,000 labelled traces from the Italian seismic dataset (INSTANCE) and is cross-validated on 22,000 traces to choose the most optimal set of hyperparameters. We obtain an accuracy of 0.98 on a completely unseen test dataset of almost 33,000 traces. Furthermore, we check the model generalizability by testing it on the datasets provided by previous works and show that our model achieves a higher recall on both positive and negative polarities.
Highlights
• Germany plans more long-distances water transfers to secure drinking water supply.
• Long-distance water transfers can unfold lock-ins that limit adaptive water governance.
• Our interdisciplinary case study shows how lock-ins emerge over different spaces and times.
• Commercialisation of water but also local protests contributed to various lock-ins.
• We therefore call for context-specific assessments of potentials and risks of LDWT.
Abstract
Germany plans to expand water transfers over long distances in the light of numerous and pressing challenges for drinking water supply. Research on inter- and intrabasin water transfers warns, however, that major investments in large-scale infrastructure systems accompanied by institutional logics and political interests often lead to a so-called lock-in. As a consequence, long-distance water transfers can limit the potential for adaptive water governance in the involved supply areas over decades with negative impacts for people and the environment. By using a case study in Germany as an example, we researched when, where and how such lock-ins around long-distance water transfers emerge. In the infrastructural development of the Elbaue-Ostharz transfer system we found various lock-ins that overlap in space and time. Some are located at the centre others at the margins of the infrastructure and commercialization of the water sector as well as hydraulic and hygienic concerns interlock with local protests in a way that the expansion of the long-distance water transfer infrastructure is presented continuously as imperative. Our findings contribute to a relational understanding of lock-ins of long-distance water transfers as contingent and diverse processes. Given the widespread occurrence of lock-ins, we argue for a context-specific assessment of potentials and risks of long-distance water transfers in times of multiple crises.
Highlights
• New fumarole and thermal water data for Askja and Kverkfjöll volcanoes, Iceland.
• Data compared to modelled compositions and fluxes of magmatic gas.
• Fumarole compositions compatible with origin of CO2 and S from degassing intrusions.
• Intrusive magmatic fluxes sufficient to sustain hydrothermal fluxes of CO2 and S in Iceland
• Magma degassing insignificant/minor source of H2O and Cl to Icelandic hydrothermal fluids
Abstract
Mantle volatiles are transported to Earth's crust and surface by basaltic volcanism. During subaerial eruptions, vast amounts of carbon, sulfur and halogens can be released to the atmosphere during a short time-interval, with impacts ranging in scale from the local environment to the global climate. By contrast, passive volatile release at the surface originating from magmatic intrusions is characterized by much lower flux, yet may outsize eruptive volatile quantities over long timescales. Volcanic hydrothermal systems (VHSs) act as conduits for such volatile release from degassing intrusions and can be used to gauge the contribution of intrusive magmatism to global volatile cycles. Here, we present new compositional and isotopic (δD and δ18O-H2O, 3He/4He, δ13C-CO2, Δ33S-δ34S-H2S and SO4) data for thermal waters and fumarole gases from the Askja and Kverkfjöll volcanoes in central Iceland. We use the data together with magma degassing modelling and mass balance calculations to constrain the sources of volatiles in VHSs and to assess the role of intrusive magmatism to the volcanic volatile emission budgets in Iceland.
The CO2/ΣS (10−30), 3He/4He (8.3–10.5 RA; 3He/4He relative to air), δ13C-CO2 (−4.1 to −0.2 ‰) and Δ33S-δ34S-H2S (−0.031 to 0.003 ‰ and −1.5 to +3.6‰) values in high-gas flux fumaroles (CO2 > 10 mmol/mol) are consistent with an intrusive magmatic origin for CO2 and S at Askja and Kverkfjöll. We demonstrate that deep (0.5–5 kbar, equivalent to ∼2–18 km crustal depth) decompression degassing of basaltic intrusions in Iceland results in CO2 and S fluxes of 330–5060 and 6–210 kt/yr, respectively, which is sufficient to account for the estimated CO2 flux of Icelandic VHSs (3365–6730 kt/yr), but not the VHS S flux (220–440 kt/yr). Secondary, crystallization-driven degassing from maturing intrusions and leaching of crustal rocks are suggested as additional sources of S. Only a minor proportion of the mantle flux of Cl is channeled via VHSs whereas the H2O flux remains poorly constrained, because magmatic signals in Icelandic VHSs are masked by a dominant shallow groundwater component of meteoric water origin. These results suggest that the bulk of the mantle CO2 and S flux to the atmosphere in Iceland is supplied by intrusive, not eruptive magmatism, and is largely vented via hydrothermal fields.
Highlights
• Subcrustal earthquakes detected beneath Fogo volcano, Cape Verde.
• At the focal depth of 40 km temperatures are likely too high for brittle failure.
• The earthquakes may originate from magma injection into a deep subcrustal reservoir.
• This observation indicates a distinct magma supply system of Fogo volcano.
Abstract
Fogo volcano belongs to the Cape Verde hotspot and its most recent eruption occurred from November 2014 to February 2015. From January to December 2016 we operated a temporary seismic network and array on Fogo and were able to locate 289 earthquakes in total. Array analysis shows that most of the events occur within the crust at distances >25 km near the neighboring island of Brava. However, on 15th August 2016 the network recorded an isolated cluster of >20 earthquakes, 13 of which could be located beneath the southern part of Fogo. The differences between S- and P-wave arrival times at steep incidence clearly indicate focal depths between approximately 38 and 44 km whereas receiver-function analyses place the Moho discontinuity at depths between 11 and 14 km. Thus, the earthquakes are located well within the upper mantle directly beneath Fogo. In view of the elevated upper-mantle temperatures within a hotspot regime, we propose that fracturing induced by magmatic injection is the most likely cause for the observed deep earthquakes.
Highlights
• Full automatized analysis of teleseismic XKS shear wave splitting.
• Rapid analysis of large seismological data sets.
• Automated window selection and quality classification.
• Application to the USArray Transportable Array including expansion to Alaska.
• Improved statistical evidence and objectivity of derived effective splitting.
Abstract
Recent technological advances have led to community wide use of large-scale seismic experiments which produce seismic data on previously impossible scales. Standard processing procedures thus require automatization to facilitate a fast and objective analysis of the data. Among these, XKS-splitting is an important tool to derive first insights into the Earth's deformation regimes at depth by studying seismic anisotropy. Most often, shear-wave splitting is interpreted to represent crystallographic preferred orientation (CPO) of mantle minerals like olivine as dominating feature and can thus be used as a proxy of mantle flow processes. Here, we introduce an addition to the MATLAB®-based SplitRacer tool box (Reiss and Rümpker 2017) which automatizes the entire XKS-splitting procedure. This is achieved by the automatization of 1) choosing a time window based on spectral analyses and 2) categorization of results based on three different XKS-splitting methods (energy minimization, rotation correlation and splitting intensity). This provides effective and objective results for splitting as well as null-measurement results. This extension allows to use SplitRacer without a graphical interface and introduces a bootstrapping statistics as error estimate of the single layer joint splitting method. The procedures are designed to allow a fast and more objective analysis of a vast amount of data, as produced by recent seismic deployments (e.g. USArray, AlpArray). We test this automatization by applying the analysis to the USArray data set, which has approximately 1900 stations with between two to fifteen years of data. We can reproduce the general pattern of the results from former studies with the more objective automatic analysis. Based on a joint-splitting approach, we approximate the splitting effect at individual stations by a single anisotropic layer. As we include null-measurements as well as a larger data set as previous studies, we can provide improved statistical evidence for these effective splitting parameters.
Highlights
• We show the first observations of seismo-acoustic tremor at Oldoinyo Lengai, the world's only active carbonatite volcano.
• We observe significant changes in seismic and acoustic tremor properties and their correlation in one year of data collection.
• Using satellite-based thermal data, we identify different volcanic processes (degassing, lava pond dynamics and spattering).
Abstract
We analyze volcanic tremor from Oldoinyo Lengai, Tanzania, which is currently the only active volcano on Earth producing carbonatitic lavas. Here, we use data from the recent SEISVOL deployment and focus on a co-located seismic and infrasound station about 200 m below the summit. We show the very first observations of seismo-acoustic tremor caused by carbonatitic eruptions. This seismo-acoustic tremor is highly variable throughout the ∼one year of data which we characterize by analyzing its seismic amplitude, duration, recurrence, dominant seismic frequency and harmonics. Frequency gliding occurs frequently and over short (minutes to hours) to long time scales (hours to days) and likely reflects different time-dependent mechanisms, such as evenly-spaced repeating events with a change in inter-event times, crater dynamics that alter resonators, and dike intrusions. Seismic and acoustic wavefields correlate well for stronger eruptive sequences but are only partially coherent which suggests that high-frequency seismic tremor (up to 25 Hz) may be caused by the low viscosity of the carbonatitic melt and not by ground-coupled airwaves. In addition, the comparison between seismic-acoustic and satellite InfraRed thermal data allows us to infer different volcanic activity styles which partially alternate throughout the year: intrusive activity and the construction of hornitos, degassing, activity from a lava pond, and varying styles of extrusive activity, in particular spattering. Our study provides important insights into the eruption dynamics of this peculiar volcano which suggests shallow melt storage within the crater floor.
The metasomatised continental mantle may play a key role in the generation of some ore deposits, in particular mineral systems enriched in platinum-group elements (PGE) and Au. The cratonic lithosphere is the longest-lived potential source for these elements, but the processes that facilitate their pre-concentration in the mantle and their later remobilisation to the crust are not yet well-established. Here, we report new results on the petrography, major-element, and siderophile- and chalcophile-element composition of native Ni, base metal sulphides (BMS), and spinels in a suite of well-characterised, highly metasomatised and weakly serpentinised peridotite xenoliths from the Bultfontein kimberlite in the Kaapvaal Craton, and integrate these data with published analyses. Pentlandite in polymict breccias (failed kimberlite intrusions at mantle depth) has lower trace-element contents (e.g., median total PGE 0.72 ppm) than pentlandite in phlogopite peridotites and Mica-Amphibole-Rutile-Ilmenite-Diopside (MARID) rocks (median 1.6 ppm). Spinel is an insignificant host for all elements except Zn, and BMS and native Ni account for typically <25% of the bulk-rock PGE and Au. High bulk-rock Te/S suggest a role for PGE-bearing tellurides, which, along with other compounds of metasomatic origin, may host the missing As, Ag, Cd, Sb, Te and, in part, Bi that are unaccounted for by the main assemblage.
The close spatial relationship between BMS and metasomatic minerals (e.g., phlogopite, ilmenite) indicates that the lithospheric mantle beneath Bultfontein was resulphidised by metasomatism after initial melt depletion during stabilisation of the cratonic lithosphere. Newly-formed BMS are markedly PGE-poor, as total PGE contents are <4.2 ppm in pentlandite from seven samples, compared to >26 ppm in BMS in other peridotite xenoliths from the Kaapvaal craton. This represents a strong dilution of the original PGE abundances at the mineral scale, perhaps starting from precursor PGE alloy and small volumes of residual BMS. The latter may have been the precursor to native Ni, which occurs in an unusual Ni-enriched zone in a harzburgite and displays strongly variable, but overall high PGE abundances (up to 81 ppm). In strongly metasomatised peridotites, Au is enriched relative to Pd, and was probably added along with S. A combination of net introduction of S, Au +/− PGE from the asthenosphere and intra-lithospheric redistribution, in part sourced from subducted materials, during metasomatic events may have led to sulphide precipitation at ~80–120 km beneath Bultfontein. This process locally enhanced the metallogenic fertility of this lithospheric reservoir. Further mobilisation of the metal budget stored in these S-rich domains and upwards transport into the crust may require interaction with sulphide-undersaturated melts that can dissolve sulphides along with the metals they store.
We investigate the applicability of the well-known multilevel Monte Carlo (MLMC) method to the class of density-driven flow problems, in particular the problem of salinisation of coastal aquifers. As a test case, we solve the uncertain Henry saltwater intrusion problem. Unknown porosity, permeability and recharge parameters are modelled by using random fields. The classical deterministic Henry problem is non-linear and time-dependent, and can easily take several hours of computing time. Uncertain settings require the solution of multiple realisations of the deterministic problem, and the total computational cost increases drastically. Instead of computing of hundreds random realisations, typically the mean value and the variance are computed. The standard methods such as the Monte Carlo or surrogate-based methods are a good choice, but they compute all stochastic realisations on the same, often, very fine mesh. They also do not balance the stochastic and discretisation errors. These facts motivated us to apply the MLMC method. We demonstrate that by solving the Henry problem on multi-level spatial and temporal meshes, the MLMC method reduces the overall computational and storage costs. To reduce the computing cost further, parallelization is performed in both physical and stochastic spaces. To solve each deterministic scenario, we run the parallel multigrid solver ug4 in a black-box fashion.
Die Welt im Wasserstress
(2024)
Wie haben sich die Wasserresourcen in den letzten 120 Jahren verändert? Und was passiert, wenn es bis Ende des 21. Jahrhunderts noch einmal zwei Grad wärmer wird als heute? Fragen wie diese beantwortet das globale Wasser-Modell WaterGAP, das maßgeblich vom Institut für Physische Geographie der Goethe-Universität und von der Ruhr-Universität Bochum entwickelt wird. Bislang ließen sich die damit erzeugten Daten nur von Expertinnen und Experten nutzen. Eine neue Web-App ändert das nun. Entwickelt wurde sie von dem französischen Geodaten-Unternehmen Ageoce, das dafür mit der Goethe-Universität kooperierte.
Reconstructing Oligocene-Miocene paleoelevation contributes to our understanding of the evolutionary history of the European Alps and sheds light on geodynamic and Earth’s surface processes involved in the development of Alpine topography. Despite being one of the most intensively explored mountain ranges worldwide, constraints on the elevation history of the European Alps, however, remain scarce. Here we present stable and clumped isotope geochemistry 15 measurements to provide a new paleoelevation estimate for the mid-Miocene (~14.5 Ma) European Central Alps. We apply stable isotope δ-δ paleoaltimetry on near sea level pedogenic carbonate oxygen isotope (δ18O) records from the Northern Alpine Foreland Basin (Swiss Molasse Basin) and high-Alpine phyllosilicate hydrogen isotope (δD) records from the Simplon Fault Zone (Swiss Alps). We further explore Miocene paleoclimate and paleoenvironmental conditions in the Swiss Molasse Basin through carbonate stable (δ18O, δ13C) and clumped (Δ47) isotope data from three foreland basin sections in different 20 alluvial megafan settings (proximal, mid-fan, and distal). Combined pedogenic carbonate δ18O values and Δ47 temperatures (30 ± 5°C) yield a near sea level precipitation δ18Ow value of -5.8 ± 0.2‰ and in conjunction with the high-Alpine phyllosilicate δD record suggest that the region surrounding the SFZ attained surface elevations of >4000 m no later than the mid-Miocene. Our near sea level δ18Ow estimate is supported by paleoclimate (iGCM Echam5-wiso) modeled δ18O values, which vary between -4.2 and -7.6‰ for the Northern Alpine Foreland Basin.
Stratospheric inorganic chlorine (Cly) is predominantly released from long-lived chlorinated source gases and, to a small extent, very short-lived chlorinated substances. Cly includes the reservoir species (HCl and ClONO2) and active chlorine species (i.e., ClOx). The active chlorine species drive catalytic cycles that deplete ozone in the polar winter stratosphere. This work presents calculations of inorganic chlorine (Cly) derived from chlorinated source gas measurements on board the High Altitude and Long Range Research Aircraft (HALO) during the Southern Hemisphere Transport, Dynamic and Chemistry (SouthTRAC) campaign in austral late winter and early spring 2019. Results are compared to Cly in the Northern Hemisphere derived from measurements of the POLSTRACC-GW-LCYCLE-SALSA (PGS) campaign in the Arctic winter of 2015/2016. A scaled correlation was used for PGS data, since not all source gases were measured. Using the SouthTRAC data, Cly from a scaled correlation was compared to directly determined Cly and agreed well. An air mass classification based on in situ N2O measurements allocates the measurements to the vortex, the vortex boundary region, and midlatitudes. Although the Antarctic vortex was weakened in 2019 compared to previous years, Cly reached 1687±19 ppt at 385 K; therefore, up to around 50 % of total chlorine was found in inorganic form inside the Antarctic vortex, whereas only 15 % of total chlorine was found in inorganic form in the southern midlatitudes. In contrast, only 40 % of total chlorine was found in inorganic form in the Arctic vortex during PGS, and roughly 20 % was found in inorganic form in the northern midlatitudes. Differences inside the two vortices reach as much as 540 ppt, with more Cly in the Antarctic vortex in 2019 than in the Arctic vortex in 2016 (at comparable distance to the local tropopause). To our knowledge, this is the first comparison of inorganic chlorine within the Antarctic and Arctic polar vortices. Based on the results of these two campaigns, the differences in Cly inside the two vortices are substantial and larger than the inter-annual variations previously reported for the Antarctic.
Biogenic organic precursors play an important role in atmospheric new particle formation (NPF). One of the major precursor species is α-pinene, which upon oxidation can form a suite of products covering a wide range of volatilities. Highly oxygenated organic molecules (HOMs) comprise a fraction of the oxidation products formed. While it is known that HOMs contribute to secondary organic aerosol (SOA) formation, including NPF, they have not been well studied in newly formed particles due to their very low mass concentrations. Here we present gas- and particle-phase chemical composition data from experimental studies of α-pinene oxidation, including in the presence of isoprene, at temperatures (−50 and −30 ∘C) and relative humidities (20 % and 60 %) relevant in the upper free troposphere. The measurements took place at the CERN Cosmics Leaving Outdoor Droplets (CLOUD) chamber. The particle chemical composition was analyzed by a thermal desorption differential mobility analyzer (TD-DMA) coupled to a nitrate chemical ionization–atmospheric pressure interface–time-of-flight (CI-APi-TOF) mass spectrometer. CI-APi-TOF was used for particle- and gas-phase measurements, applying the same ionization and detection scheme. Our measurements revealed the presence of C8−10 monomers and C18−20 dimers as the major compounds in the particles (diameter up to ∼ 100 nm). Particularly, for the system with isoprene added, C5 (C5H10O5−7) and C15 compounds (C15H24O5−10) were detected. This observation is consistent with the previously observed formation of such compounds in the gas phase. However, although the C5 and C15 compounds do not easily nucleate, our measurements indicate that they can still contribute to the particle growth at free tropospheric conditions. For the experiments reported here, most likely isoprene oxidation products enhance the growth of particles larger than 15 nm. Additionally, we report on the nucleation rates measured at 1.7 nm (J1.7 nm) and compared with previous studies, we found lower J1.7 nm values, very likely due to the higher α-pinene and ozone mixing ratios used in the present study.
Reconstructing Oligocene–Miocene paleoelevation contributes to our understanding of the evolutionary history of the European Alps and sheds light on geodynamic and Earth surface processes involved in the development of Alpine topography. Despite being one of the most intensively explored mountain ranges worldwide, constraints on the elevation history of the European Alps remain scarce. Here we present stable and clumped isotope measurements to provide a new paleoelevation estimate for the mid-Miocene (∼14.5 Ma) European Central Alps. We apply stable isotope δ–δ paleoaltimetry to near-sea-level pedogenic carbonate oxygen isotope (δ18O) records from the Northern Alpine Foreland Basin (Swiss Molasse Basin) and high-Alpine phyllosilicate hydrogen isotope (δD) records from the Simplon Fault Zone (Swiss Alps). We further explore Miocene paleoclimate and paleoenvironmental conditions in the Swiss Molasse Basin through carbonate stable (δ18O, δ13C) and clumped (Δ47) isotope data from three foreland basin sections in different alluvial megafan settings (proximal, mid-fan, and distal). Combined pedogenic carbonate δ18O values and Δ47 temperatures (30±5 ∘C) yield a near-sea-level precipitation δ18Ow value of ‰ and, in conjunction with the high-Alpine phyllosilicate δD value of ‰, suggest that the region surrounding the Simplon Fault Zone attained surface elevations of >4000 m no later than the mid-Miocene. Our near-sea-level δ18Ow estimate is supported by paleoclimate (iGCM ECHAM5-wiso) modeled δ18O values, which vary between −4.2 ‰ and −7.6 ‰ for the Northern Alpine Foreland Basin.
Stratospheric inorganic chlorine (Cly) is predominantly released from long-lived chlorinated source gases and, to a small extent, very short-lived chlorinated substances. Cly includes the reservoir species (HCl and ClONO2) and active chlorine species (i.e., ClOx). The active chlorine species drive catalytic cycles that deplete ozone in the polar winter stratosphere. This work presents calculations of inorganic chlorine (Cly) derived from chlorinated source gas measurements on board the High Altitude and Long Range Research Aircraft (HALO) during the Southern Hemisphere Transport, Dynamic and Chemistry (SouthTRAC) campaign in austral late winter and early spring 2019. Results are compared to Cly in the Northern Hemisphere derived from measurements of the POLSTRACC-GW-LCYCLE-SALSA (PGS) campaign in the Arctic winter of 2015/2016. A scaled correlation was used for PGS data, since not all source gases were measured. Using the SouthTRAC data, Cly from a scaled correlation was compared to directly determined Cly and agreed well. An air mass classification based on in situ N2O measurements allocates the measurements to the vortex, the vortex boundary region, and midlatitudes. Although the Antarctic vortex was weakened in 2019 compared to previous years, Cly reached 1687±19 ppt at 385 K; therefore, up to around 50 % of total chlorine was found in inorganic form inside the Antarctic vortex, whereas only 15 % of total chlorine was found in inorganic form in the southern midlatitudes. In contrast, only 40 % of total chlorine was found in inorganic form in the Arctic vortex during PGS, and roughly 20 % was found in inorganic form in the northern midlatitudes. Differences inside the two vortices reach as much as 540 ppt, with more Cly in the Antarctic vortex in 2019 than in the Arctic vortex in 2016 (at comparable distance to the local tropopause). To our knowledge, this is the first comparison of inorganic chlorine within the Antarctic and Arctic polar vortices. Based on the results of these two campaigns, the differences in Cly inside the two vortices are substantial and larger than the inter-annual variations previously reported for the Antarctic.
Biogenic organic precursors play an important role in atmospheric new particle formation (NPF). One of the major precursor species is α-pinene, which upon oxidation can form a suite of products covering a wide range of volatilities. Highly oxygenated organic molecules (HOMs) comprise a fraction of the oxidation products formed. While it is known that HOMs contribute to secondary organic aerosol (SOA) formation, including NPF, they have not been well studied in newly formed particles due to their very low mass concentrations. Here we present gas- and particle-phase chemical composition data from experimental studies of α-pinene oxidation, including in the presence of isoprene, at temperatures (−50 and −30 ∘C) and relative humidities (20 % and 60 %) relevant in the upper free troposphere. The measurements took place at the CERN Cosmics Leaving Outdoor Droplets (CLOUD) chamber. The particle chemical composition was analyzed by a thermal desorption differential mobility analyzer (TD-DMA) coupled to a nitrate chemical ionization–atmospheric pressure interface–time-of-flight (CI-APi-TOF) mass spectrometer. CI-APi-TOF was used for particle- and gas-phase measurements, applying the same ionization and detection scheme. Our measurements revealed the presence of C8−10 monomers and C18−20 dimers as the major compounds in the particles (diameter up to ∼ 100 nm). Particularly, for the system with isoprene added, C5 (C5H10O5−7) and C15 compounds (C15H24O5−10) were detected. This observation is consistent with the previously observed formation of such compounds in the gas phase. However, although the C5 and C15 compounds do not easily nucleate, our measurements indicate that they can still contribute to the particle growth at free tropospheric conditions. For the experiments reported here, most likely isoprene oxidation products enhance the growth of particles larger than 15 nm. Additionally, we report on the nucleation rates measured at 1.7 nm (J1.7 nm) and compared with previous studies, we found lower J1.7 nm values, very likely due to the higher α-pinene and ozone mixing ratios used in the present study.
Although global- and catchment-scale hydrological models are often shown to accurately simulate long-term runoff time-series, far less is known about their suitability for capturing hydrological extremes, such as droughts. Here we evaluated simulations of hydrological droughts from nine catchment scale hydrological models (CHMs) and eight global scale hydrological models (GHMs) for eight large catchments: Upper Amazon, Lena, Upper Mississippi, Upper Niger, Rhine, Tagus, Upper Yangtze and Upper Yellow. The simulations were conducted within the framework of phase 2a of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP2a). We evaluated the ability of the CHMs, GHMs and their respective ensemble means (Ens-CHM and Ens-GHM) to simulate observed hydrological droughts of at least one month duration, over 31 years (1971–2001). Hydrological drought events were identified from runoff-deficits and the Standardised Runoff Index (SRI). In all catchments, the CHMs performed relatively better than the GHMs, for simulating monthly runoff-deficits. The number of drought events identified under different drought categories (i.e. SRI values of -1 to -1.49, -1.5 to -1.99, and ≤-2) varied significantly between models. All the models, as well as the two ensemble means, have limited abilities to accurately simulate drought events in all eight catchments, in terms of their occurrence and magnitude. Overall, there are opportunities to improve both CHMs and GHMs for better characterisation of hydrological droughts.
Lightning climate change projections show large uncertainties caused by limited empirical knowledge and strong assumptions inherent to coarse-grid climate modeling. This study addresses the latter issue by implementing and applying the lightning potential index parameterization (LPI) into a fine-grid convection-permitting regional climate model (CPM). This setup takes advantage of the explicit representation of deep convection in CPMs and allows for process-oriented LPI inputs such as vertical velocity within convective cells and coexistence of microphysical hydrometeor types, which are known to contribute to charge separation mechanisms. The LPI output is compared to output from a simpler flash rate parameterization, namely the CAPE × PREC parameterization, applied in a non-CPM on a coarser grid. The LPI’s implementation into the regional climate model COSMO-CLM successfully reproduces the observed lightning climatology, including its latitudinal gradient, its daily and hourly probability distributions, and its diurnal and annual cycles. Besides, the simulated temperature dependence of lightning reflects the observed dependency. The LPI outperforms the CAPE × PREC parameterization in all applied diagnostics. Based on this satisfactory evaluation, we used the LPI to a climate change projection under the RCP8.5 scenario. For the domain under investigation centered over Germany, the LPI projects a decrease of 4.8% in flash rate by the end of the century, in opposition to a projected increase of 17.4% as projected using the CAPE × PREC parameterization. The future decrease of LPI occurs mostly during the summer afternoons and is related to (i) a change in convection occurrence and (ii) changes in the microphysical mixing. The two parameterizations differ because of different convection occurrences in the CPM and non-CPM and because of changes in the microphysical mixing, which is only represented in the LPI lightning parameterization.
Global analysis of halogenated trace gases in the UTLS: from long-lived to short-lived substances
(2023)
In this dissertation, the distribution of chlorinated and brominated substances in the upper troposphere and lower stratosphere is investigated. These substances contribute significantly to the catalytic decomposition of ozone and are involved in the recurrent formation of the polar ozone hole in the Antarctic winter and spring. The Montreal Protocol, a multilateral environmental treaty to protect the ozone layer, has successfully reduced emissions of long-lived chlorine- and bromine-containing substances. Short-lived chlorinated and brominated substances, some of which are natural and anthropogenic in origin, are not regulated by the Montreal Protocol and it can be assumed that their relative contribution to the stratospheric halogen budget will increase, while the contribution of long-lived compounds will steadily decrease. The distribution of long- and short-lived halogenated substances are part of current research. For the upper troposphere and lower stratosphere, the very short-lived substances are particularly important. The lower stratosphere needs special investigation in this respect, since its composition is influenced by different transport processes. The influences on ozone trends in the lower stratosphere are subject to great uncertainties. Especially in the Southern Hemisphere, the number of observations is very limited.
In this work, the GhOST (Gas chromatograph for Observational Studies using Tracers) instrument was used during the SouthTRAC measurement campaign on the German HALO (High Altitude and LOng range) research aircraft, providing observations of halogenated hydrocarbons in Antarctic late winter to early spring 2019, a generally poorly sampled region. The polar vortex was, compared to previous years, significantly weaker and shifted towards the eastern South Pacific and South America. From the airborne measurements of chlorinated source gases, inorganic chlorine (the sum of active chlorine and reservoir gases; Cly) could be inferred with the result that Cly within the vortex increased up to 1687 ± 19 ppt at 385 K potential temperature, accounting for about 50 % of the total chlorine within the vortex and only 15 % of the total chlorine in the southern mid-latitudes. A comparison with the Northern Hemisphere could be made using the PGS measurement campaign in the Arctic winter 2015/2016. Under comparable conditions (season and distance from the tropopause), only 40 % of the total chlorine was in the inorganic form within the Arctic polar vortex and about 20 % was found in the mid-latitudes of the Northern Hemisphere. In addition, about 540 ppt more Cly was present in the Antarctic vortex than in the Arctic vortex, exceeding the annual variations previously reported for Antarctica.
The mean age of air plays an essential role in the derivation of Cly via the organic source gases, as was done in this work. A new method for determining the mean age of air from observational data has been introduced that accounts for extra-tropical input to the stratosphere in addition to tropical input. This new method was compared with the previously used method, which considered only the tropical input. The new method shows more realistic values especially near the tropopause. On average, the air of the lower stratosphere in the Northern Hemisphere was older than in the Southern Hemisphere by about 0.5 ± 0.3 years. About 65 K above the tropopause, the pattern changed with older air in the mid-latitudes of the Northern Hemisphere, but older air in high latitudes of the Southern Hemisphere, which implies differences in the strength and isolation of the respective polar vortex as well as the wave forcing in the shallow branch of the Brewer-Dobson circulation of the respective hemisphere. This is in good agreement with the distribution of Cly. The difference in the lower stratosphere was not clearly evident with the old method and it can be assumed that investigations of the differences in Cly of Northern and Southern Hemisphere will benefit from the new method.
Finally, the global and seasonal distribution of the two most important representatives of the short-lived brominated substances, CH2Br2 and CHBr3, was investigated. For this purpose, two additional HALO measurement campaigns have been used, the 2012 TACTS measurement campaign and the 2017 WISE measurement campaign, as well as the HIAPER Pole-to-Pole Observations (HIPPO) and Atmospheric Tomography (ATom) measurement campaigns. Observations of CH2Br2 show a pronounced seasonality in the free and upper troposphere of both hemispheres with slightly larger values in the Northern Hemisphere. CHBr3, on the other hand, shows a generally higher variability and lower seasonality with larger mixing ratios at mid and high latitudes in the northern hemispheric winter and autumn. A comparison of the lower stratosphere is limited to autumn and spring of both hemispheres due to the limited data basis of the observations. The distributions in each spring are similar (less than 0.1 ppt differences for e.g., CH2Br2). In hemispheric autumn, larger differences are evident with substantially smaller mixing ratios in the southern hemispheric lower stratosphere. This suggests that the transport processes of the two hemispheres may be different and implies that the input of tropospheric air (flushing) to the Northern Hemisphere lowest stratosphere is more efficient than in the Southern Hemisphere. Vertical profiles of CH2Br2 and CHBr3 in the mid-latitudes of both hemispheres and resulting vertical gradients support this conjecture. However, the Southern Hemisphere data set is insufficient to quantify this difference and further measurements are needed.
Semi-arid African ecosystems influence trends and variability in global terrestrial carbon dynamics. However, there are uncertainties in potential effects of future climates for semi-arid ecosystems, especially for niche ecosystems. At the same time, African ecosystems provide the livelihoods and ecosystem services for around 1.4 billion people. Future population growth and associated changes in land use pose a challenge for the protection of African biodiversity. Therefore, this work focussed on future impacts of climate change on African ecosystems and carbon dynamics and also for African protected areas (PAs), where they may cooccur with other global change factors. Another focus was on uncertainties associated with future projections and with modelling the Nama Karoo, as an example of a semi-arid niche ecosystem. Dynamic vegetation models (DVMs) were the main research tool.
In Chapter 2, we analysed climate change impacts on African ecosystems and carbon pools until the end of the 21st century and associated uncertainties based on an ensemble of vegetation simulations with the DVM adaptive dynamic vegetation model (aDGVM). We investigated the impact of increased atmospheric CO2 concentrations and two climate change scenarios (medium (RCP4.5) and high emissions (RCP8.5); RCP - representative concentration pathway) on vegetation changes. Differences in the simulated vegetation were primarily driven by assumptions about the influence of CO2 on plants. Elevated CO2 concentrations led to increased total aboveground vegetation biomass and shrub encroachment into grasslands and savannas for both climate scenarios. In simulations without the direct influence of CO2 on plants, there was hardly any shrub encroachment and vegetation biomass decreased or varied between a slight decrease in some cases and a slight increase in others. Based on these results, biome changes due to climate change are likely in Africa in the future. Due to the large uncertainties in future projections, strategies to adapt to climate change must be flexible.
The simulated vegetation in Chapter 2 represented potential, natural vegetation and is particularly suitable to investigate PAs. However, PAs do not exist isolated from their environment and social developments. In Chapter 3, the vegetation projections with CO2 effect from Chapter 2 were combined with projections for population density and land use. Except for many PAs in North Africa, most PAs were adversely affected by at least one of the three drivers by the end of the 21st century in both investigated scenarios ("middle-of-the-road" and "fossil-fuelled development"). Cooccurrence of the drivers varied by region and scenario for PAs. Both scenarios implied increasing challenges for the conservation of African biodiversity in PAs. The impact of climate change on vegetation is likely to be exacerbated by socio-economic change for most African PAs. Strong mitigation of future climate change together with equitable societal development may facilitate successful ecosystem conservation.
The simulations in Chapters 2 and 3 showed large-scale patterns of vegetation change, but their low resolution makes them unsuitable for local analyses. In Chapter 4, the challenges of simulating smaller scale, semi-arid ecosystems and their carbon cycle were analysed for the Nama Karoo with the aDGVM2 and its shrub module. The aDGVM2 is based on the aDGVM, but represents plants more flexibly. In all tested aDGVM2 configurations, the carbon fluxes improved compared to initial simulations but still overestimated them. The measured morphology of the dwarf shrubs and soil water dynamics were not reproduced in aDGVM2. Semi-arid soil water dynamics and coping strategies of semi-arid dwarf shrubs under drought stress are not adequately implemented in the aDGVM2. Further field research on semi-arid water and carbon dynamics of vegetation is necessary to parameterise the aDGVM2 for dwarf shrubs. If these challenges are overcome, DVMs can be a powerful tool for much-needed research on the impacts of climate change on the Nama Karoo.
The analyses have shown that climate change under medium to high emission scenarios is likely to lead to large-scale changes in ecosystems and the carbon balance in Africa. Because lower emissions scenarios come with less uncertainty, climate change adaptation strategies likely need to be less complex or extensive if climate change is minimised. For African PAs, the challenges of climate change may be exacerbated by socio-economic factors to a regionally varying extent. This research suggests that successful ecosystem conservation depends on climate change mitigation measures and ensuring equitable, sustainable development. The shown uncertainties, e.g., in the implementation of the CO2 effect on plants or vegetation dynamics in more niche ecosystems, help to focus future research efforts and increase our understanding of the range of plausible futures we may need to adapt to.
First-principles modeling techniques offer the ability to simulate a wide range of systems under different physical conditions, such as temperature, pressure, and composition, without relying on empirical knowledge. Density functional theory (DFT), a quantum mechanical method, has become an exceptionally successful framework for materials science modeling. Employing DFT makes it possible to gain valuable insights into the fundamental state of a system, enabling the reliable determination of equilibrium crystal structures. Over time, DFT has become an essential tool that can be incorporated into various schemes for predicting the properties of a material related to its structure, insulating/metallic behavior, magnetism, and optics. DFT is regularly applied in numerous fields, spanning from fundamental subjects in condensed matter physics to the study of large-scale phenomena in geosciences. In the latter, the effectiveness of DFT stems from its ability to simulate the properties found on the Earth, other planets, and meteorites, which may pose challenges for their direct study or laboratory investigation.
In this thesis, a comprehensive examination of a family of monosulfides and a perovskite heterostructure was conducted. These materials are relevant for their potential applications in technology, energy harvesting, and in the case of monosulfides, their speculated abundance on the planet Mercury.
Firstly, a DFT approach was used to analyze two non-magnetic monosulfides, CaS and MgS. We determined their structural properties and then focused on the modeling of their reflectivity in the infrared region. The calculation of the reflectivity considered both harmonic and anharmonic contributions. In the harmonic limit, the non-analytic correction was employed to accurately determine the LO/TO splitting, which is necessary to delimit the retstrahlend band, that is, the maximum of the reflectivity. The anharmonic effects given by up to three-phonon and isotopic scatterings, which were included using perturbation theory, primarily smeared the reflectivity spectra edges in the high-wave region.
Secondly, four polymorphs of MnS were studied using a combination of first-principles methods to simulate their antiferromagnetic (AFM) and paramagnetic (PM) states. The integration of DFT+$U$ with special quasirandom structures (SQS) supercells, and occupation matrix control techniques was crucial for achieving convergence, structural optimization accuracy, and obtaining finite energy band gaps and local magnetic moments in the PM phases. The addition of the Hubbard $U$ correction was necessary to treat the highly-correlated Mn $d$-electrons. The success of our approach was clear based on our electronic structure predictions for the PM rock-salt B1-MnS polymorph. Experimentally this phase has been observed to be an insulator, but multiple \emph{ab initio} works resulted previously in metallic behavior. Our computations, on the other hand, predicted insulating and magnetic properties that compare well with available measurements. Additionally, the pressure-field stability of the four MnS polymorphs was studied. In the case of the PM phases, B1-MnS was identified to be the most stable up to about 21 GPa, then transforming into the B31-MnS polymorph. This finding was in close agreement with high-pressure experiments reporting a similar phase transformation. The optical properties of B1-, B4-, and B31-MnS were also simulated. The SQS technique was used to obtain soft-mode-free phonon band structures within the harmonic approximation. Then, the anharmonic effects were included, and the reflectivity was calculated for B1-MnS and B4-MnS. In both cases, a good agreement for the LO/TO splitting was achieved in comparison to experimental results.
Lastly, the oxygen-deficient heterostructure of LaAlO$_{3-\delta}$ /SrTiO$_{3-\delta}$ was investigated also employing DFT+$U$, with a particular emphasis on the potential impact of vacancy clustering at the interface. Six distinct configurations of pairs of vacancies were studied and their energies were compared to find the most stable one. The orbital reconstruction of Ti orbitals was also examined based on their location with respect to the vacancies and the local magnetic moments were calculated. The final results showed that linearly arranged vacancies located opposite to Ti ions give the most energetically stable configuration.
Climatology of morphology and cloud-radiative properties of marine low-level mixed-phase clouds
(2023)
Marine stratocumuli cover about 40 - 60% of the ocean surface. They self-organize into different morphological regimes. The two organized cellular regimes are called open and closed mesoscale-cellular convective (MCC) clouds. In mid-to-high latitudes, open and closed cells are the two most frequent types of MCC clouds. In particular, many MCC clouds consist of a mixture of vapor, liquid droplets, and ice particles, referred to as mixed-phase clouds (MPCs). Even for the same cloud fraction, the albedo of open cells is, on average, lower than that of closed MCC clouds. Cloud phase and morphology individually influence the cloud radiative effect. Thus, this thesis investigates the relationships between the cloud phase, MCC organization, cell size, and differences regarding the cloud-radiative effect.
This thesis focuses on space-borne retrievals to achieve extensive temporal and spatial coverage. The liDAR-raDAR (DARDAR) version 2 product collocates two active and one passive satellite: CloudSat, Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO), and Moderate Resolution Imaging Spectroradiometer (MODIS). The cloud phase of DARDAR is vertically integrated to establish a single cloud phase at each data point. The MCC classification data set based on the liquid water path (LWP) of MODIS scenes is collocated with the DARDAR product to determine the MCC organization. Cell-size statistics of both MCC clouds are obtained using a marker-based image segmentation method on MODIS reflectance scenes. In addition, based on MODIS reflectance scenes, a convolutional neural network (CNN) is developed to classify open and closed MCC scenes to avoid missing mature MPCs with a low LWP.
The first part of this thesis explores the relationships between cloud phase, morphology, and cloud albedo in the Southern Ocean (SO). At a given cloud-top temperature (CTT), seasonal changes in the mixed-phase fraction, defined as the number of MPCs divided by the sum of MPC and supercooled liquid cloud (SLC) pixels, are stronger than the morphological changes. Therefore, external factors seem to influence these changes instead of morphology. The dependence of cloud phase on cloud-top height (CTH) is more substantial than on CTT in clouds with CTHs below 2.5 km. The previously observed acceleration of closed-to-open transition in MPCs, known as preconditioning, is not the primary driver of climatological cloud morphology statistics in the SO. The morphological differences in cloud albedo are more pronounced in SLCs than in MPCs. This change in albedo alters the cloud radiative effect in the SO by 21Wm−2 to 39Wm−2 depending onseason and cloud phase.
Open and closed MCC clouds exhibit larger equivalent cell diameters in the MPCs than in SLCs in austral summer, whereas, in austral winter, the SLCs are larger. The cell’s aspect ratio accounts for varying CTHs. Closed cells have smaller aspect ratios than open cells, so their cell diameter is smaller, independent of CTH. While the seasonal differences in closed cells are due to changes in CTH, the seasonal aspect ratio differences in open cells are mainly caused by MPCs. With increasing aspect ratios, the cloud albedo decreases in both open and closed MCC clouds, with the most substantial decrease in open MPCs clouds. This leads to cloud-radiative changes of 60 - 75Wm−2 in the SO, depending on cloud phase and aspect ratio.
The established CNN exhibits a good accuracy of 80.6%, with even higher accuracies in the Open (85.5%) and Closed (87.3%) categories. The global MCC climatology based on the CNN generally agrees well with previous MCC distributions. The most notable difference occurs in the Northern Hemisphere (NH) in boreal winter, with a higher occurrence frequency of closed and open MCC clouds. This might indicate missing MPCs in previous studies based on the LWP and some restricted to warm cloud scenes. Thus, the developed CNN seems to better represent the different morphologies in MPCs than in previous classifications.
In conclusion, this thesis shows that understanding the dependencies of cloud phase, cloud morphology, and cell size is important to enhance predictions of the cloud-radiative effect and thus, it is important to evaluate how cloud phase, cloud morphology, and cellsize change in a warming climate.
Due to massive energetic investments in woody support structures, trees are subject to unique physiological, mechanical, and ecological pressures not experienced by herbaceous plants. Despite a wealth of studies exploring trait relationships across the entire plant kingdom, the dominant traits underpinning these unique aspects of tree form and function remain unclear. Here, by considering 18 functional traits, encompassing leaf, seed, bark, wood, crown, and root characteristics, we quantify the multidimensional relationships in tree trait expression. We find that nearly half of trait variation is captured by two axes: one reflecting leaf economics, the other reflecting tree size and competition for light. Yet these orthogonal axes reveal strong environmental convergence, exhibiting correlated responses to temperature, moisture, and elevation. By subsequently exploring multidimensional trait relationships, we show that the full dimensionality of trait space is captured by eight distinct clusters, each reflecting a unique aspect of tree form and function. Collectively, this work identifies a core set of traits needed to quantify global patterns in functional biodiversity, and it contributes to our fundamental understanding of the functioning of forests worldwide.
Die nachfolgende Dissertation wurde an der Goethe-Universität Frankfurt am Institut für Geowissenschaften (FB 11) in der Arbeitsgruppe Kristallographie und Mineralogie (AG Winkler) verfasst. Die Betreuung der hier durchgeführten Arbeiten erfolgte hauptsächlich durch Prof. B. Winkler in Zusammenarbeit mit Dr. L. Bayarjargal, PD Dr. E. Haussühl und PD Dr. V. Vinograd. Bei dem vorliegenden Manuskript handelt es sich um eine kumulative bzw. publikationsbasierte Dissertation, welche die Forschungsergebnisse verschiedener bereits veröffentlichter wissenschaftlicher Fachartikel zusammenfasst.
Die Arbeit beschreibt verschiedene Synthesen und Untersuchungen an Carbonaten und teilt sich im Wesentlichen in zwei Abschnitte. Zum einen wurden Experimente mit Carbonaten bei Extrembedingungen bzw. unter hohen Drücken und hohen Temperaturen durchgeführt, wie sie auch im Inneren der Erde zu finden sind. Im zweiten Teil wurden Carbonate bei Raumbedingungen synthetisiert und der Einbau von Seltenerdelementen untersucht. Grundsätzlich werden jedoch in beiden Teilen dieser Arbeit die Strukturen und Eigenschaften verschiedener Carbonate und eine mögliche Kationensubstitution bzw. die Synthese isostruktureller Verbindungen erforscht.
The Izu–Bonin–Mariana volcanic arc is situated at a convergent plate margin where subduction initiation triggered the formation of MORB-like forearc basalts as a result of decompression melting and near-trench spreading. International Ocean Discovery Program (IODP) Expedition 352 recovered samples within the forearc basalt stratigraphy that contained unusual macroscopic globular textures hosted in andesitic glass (Unit 6, Hole 1440B). It is unclear how these andesites, which are unique in a stratigraphic sequence dominated by forearc basalts, and the globular textures therein may have formed. Here, we present detailed textural evidence, major and trace element analysis, as well as B and Sr isotope compositions, to investigate the genesis of these globular andesites. Samples consist of K2O-rich basaltic globules set in a glassy groundmass of andesitic composition. Between these two textural domains a likely hydrated interface of devitrified glass occurs, which, based on textural evidence, seems to be genetically linked to the formation of the globules. The andesitic groundmass is Cl rich (ca. 3000 µg/g ), whereas globules and the interface are Cl poor (ca. 300 µg/g ). Concentrations of fluid-mobile trace elements also appear to be fractionated in that globules and show enrichments in B, K, Rb, Cs, and Tl, but not in Ba and W relative to the andesitic groundmass, whereas the interface shows depletions in the latter, but is enriched in the former. Interestingly, globules and andesitic groundmass have identical Sr isotopic composition within analytical uncertainty ( 87Sr∕86Sr of 0.70580 ± 10 ), indicating that they likely formed from the same source. However, globules show high δ11 B (ca. + 7 ‰ ), whereas their host andesites are isotopically lighter (ca. – 1 ‰ ), potentially indicating that whatever process led to their formation either introduced heavier B isotopes to the globules, or induced stable isotope fractionation of B between globules and their groundmass. Based on the bulk of the textural information and geochemical data obtained from these samples, we conclude that these andesites likely formed as a result of the assimilation of shallowly altered oceanic crust (AOC) during forearc basaltic magmatism. Assimilation likely introduced radiogenic Sr, as well as heavier B isotopes to comparatively unradiogenic and low δ11B forearc basalt parental magmas (average 87Sr∕86Sr of 0.703284). Moreover, the globular textures are consistent with their formation being the result of fluid-melt immiscibility that was potentially induced by the rapid release of water from assimilated AOC whose escape likely formed the interface. If the globular textures present in these samples are indeed the result of fluid-melt immiscibility, then this process led to significant trace element and stable isotope fractionation. The textures and chemical compositions of the globules highlight the need for future experimental studies aimed at investigating the exsolution process with respect to potential trace element and isotopic fractionation in arc magmas that have perhaps not been previously considered.
Knowledge about the initial tectonic and depositional dynamics, as well as the influence of early rifting on climate and environmental evolution remains speculative to a large extent, because sediments are usually deeply buried. Within the East African Rift System, inversion tectonics uplifted a few of these successions to the surface hence presenting rare windows into the pre-rift depositional history. One such example, an exceptional 700 m long and up to 60 m high fresh road cut provided the opportunity to study in detail initial rift successions of the southern Albertine Rift (Western Uganda). This focusses on the basal and poorly known Middle to Late Miocene in order to unravel the climatic, environmental, hydrological and tectonic evolution of the initial Albertine Rift. A large and robust multi-proxy dataset was gathered comprising 169 m of stratigraphic thickness, which spans from 14.5 to 4.9 Ma according to a revised lithostratigraphic model. Fieldwork comprised logging of the sedimentary record, spectral gamma ray, magnetic susceptibility and 2D wall mapping with photomosaics. Additionally, the sections were sampled for bulk mineral and clay mineral analysis. The succession exposes a suite of lithofacies and architectural elements detailing the evolution of a fluvio-lacustrine system. Five depositional environments were identified which show an overall back-stepping trend from an alluvial plain to a delta plain and finally palustrine/shallow lacustrine conditions. Mesoscale base-level cycles, preservation potential of architectural elements, and stacking pattern exhibit limited accommodation space. However, it increases over time. This overall trend indicates increasing tectonic subsidence, which can be explained by flexural downwarp within the pre-rift phase and in the upper part grading into fault-controlled crustal extension of the syn-rift phase, which more and more disrupted a large-scale river system. From the Middle Miocene up to the early Pliocene, this study revealed that palaeoclimate trends become marked by increasing and more fluctuating Th concentrations, loss of feldspar, intercalated lenses of hydroxosulphate minerals, and a shift from smectite-dominated to kaolinite-dominated clays. These signals are all interpreted as detrital except for the hydroxosulphates, and they mirror the increasing intensity of chemical weathering and stripping of soils in the catchment. A trend towards increasing humidity is supported by an increase in lacustrine sediment facies and a lake-level rise. Nevertheless, intercalation of hydroxosulphate, ferricretes and pedogenised horizons prove ongoing seasonality and dry intervals. Finally, based on a revised stratigraphic model a sequence stratigraphic correlation of the outcrop's depositional cycles with basin-scale cycles is presented. According to these cycles, transition from the pre-rift to the syn-rift stage is marked by an unconformity and a tectonic pulse in the latest Miocene. However, the response of fluvial supply, the depositional system as well as climate conditions are less punctuated and characterised by gradual trends and temporal delays. The long pre-rift phase (ca 10 Myr) and the gradual transition to the syn-rift phase is in accordance with the active rifting model, which is based on thermal thinning of the lithosphere by asthenospheric upwelling.
Herein, the high-pressure/high-temperature synthesis (11 GPa, 650 °C) of Tb3B10O17(OH)5 in a modified Walker-type multianvil device is presented. The structure of this rare-earth borate was determined by single-crystal X-ray diffraction methods and was found to crystallize orthorhombically in the space group Pmn21 (no. 31) with the unit cell parameters a = 16.2527(4), b = 4.4373(1), and c = 8.8174(2) Å. The new compound was further characterized using infrared spectroscopy, energy-dispersive X-ray spectroscopy, second harmonic generation (SHG) measurements, and temperature-dependent X-ray powder diffraction. Tb3B10O17(OH)5 decomposes to β-Tb(BO2)3 at temperatures higher than 460 °C. With increasing temperatures, the formation of μ-TbBO3 was observed, which transforms to π-TbBO3 upon cooling.
High-pressure/high-temperature synthesis of the new boron-rich terbium hydroxyborate Tb3B12O19(OH)7
(2023)
Monoclinic Tb3B12O19(OH)7 was obtained by multianvil high-pressure/high-temperature syntheses at 6 GPa and 650 °C. The crystal structure was investigated by single-crystal X-ray diffraction methods and space group C2 (no. 5) with the unit cell parameters a = 24.2299(5) Å, b = 4.4667(1) Å, c = 7.0964(2) Å, β = 94.58(1)°, and two formula units per cell were revealed. Powder X-ray diffraction, infrared spectroscopy and the investigation of its second harmonic generation properties support the proposed structural model.
In this dissertation, different aspects of turbulent transport and thermally driven flows over complex terrain are investigated. Two publications concentrate on the vertical heat and moisture exchange in the convective boundary layer over mountainous terrain. To study this, Large-Eddy Simulation (LES) is used. Both turbulent and advective transport mechanisms are evaluated over the simple orography of a quasi-two-dimensional, periodic valley with prescribed surface fluxes. Here, terrain elevation varies along only one of the horizontal coordinate axes. Even a relatively shallow orography, possibly unresolved in existing numerical weather prediction models, modifies the domain-averaged moisture and temperature profiles. For the analysis, the flow is decomposed into a local turbulent part, a local mean circulation, and a large-scale part. An analysis of the turbulent kinetic energy and turbulent heat and moisture flux budgets shows that the thermal circulation significantly contributes to the vertical transport. It is found that thermal upslope winds are important for the moisture transport from the valley to the mountain tops. In total, moisture export out of the valley is mostly accomplished by the mean circulation. On the temperature distribution, which is horizontally relatively homogeneous, the thermal circulation has a weaker impact. If an upper-level wind is present, it interacts with the thermal circulation. This weakens the vertical transport of moisture and thus reduces its export out of the valley. The heat transport is less affected by the upper-level wind because of its weaker dependence on the thermal circulation. These findings were corroborated in a more realistic experiment simulating the full diurnal cycle using radiation forcing and an interactive land surface model.
Based on these results, coherent turbulent structures in the convective boundary layer over non-flat terrain are studied in further detail. A conditional sampling method based on the concentration of a decaying passive tracer is implemented in order to identify the boundary-layer plumes objectively. Conditional sampling allows to quantify the contribution of plume structures to the vertical transport of heat and moisture. In case of the idealized valley, vertical transport by coherent structures is the dominant contribution to the turbulent components of both heat and moisture flux. It is comparable in magnitude to the advective transport by the mean slope-wind circulation, although it is more important for heat than for moisture transport. A set of less idealized simulations considers the flow over three-dimensional terrain. In this case, conditional sampling is carried out by using a simple domain-decomposition approach. We demonstrate that thermal updrafts are generally more frequent on hill tops than over the surroundings, but they are less persistent on the windward sides when large-scale winds are present in the free atmosphere.
The tools for flow decomposition and budget analysis are also applied in another idealized case with a quasi-two-dimensional valley featuring the stable boundary layer. Here, the formation of a low stratus cloud is investigated. The main driver for the cloud formation is radiative cooling due to outgoing longwave radiation. Despite a purely horizontal flow, the advection terms in the prognostic equations for heat and moisture produce vertical mixing across the upper cloud edge leading to a loss of cloud water content. However, this behavior is not due to any kind of thermally-driven circulation. Instead, this spurious mixing is caused by the diffusive error of the advection scheme in regions where the sloping surfaces of the terrain-following vertical coordinate intersect the cloud top. It is shown that the intensity of the (spurious) numerical diffusion strongly depends on the horizontal resolution, the order of advection, and the choice of the scalar advection scheme. A LES with 4 m horizontal resolution serves as a reference. For horizontal resolutions of a few hundred meters, carried out with a model setup as it is used in Numerical Weather Prediction, a strong reduction of the simulated liquid-water path is observed. In order to keep the (spurious) numerical diffusion at coarser resolutions small, at least a fifth-order advection scheme should be used. In the present case, a WENO scalar advection scheme turns out to increase the numerical diffusion along a sharp cloud edge compared to an upwind scheme. Furthermore, the choice of the vertical coordinate has a strong impact on the simulated liquid-water path over orography. With a modified definition of the terrain-following sigma coordinate, it is possible to produce cloud water where the classical sigma coordinate does not allow any cloud formation.
Local climate change risk assessments (LCCRAs) are best supported by a quantitative integration of physical hazards, exposures and vulnerabilities that includes the characterization of uncertainties. We propose to use Bayesian Networks (BNs) for this task and show how to integrate freely-available output of multiple global hydrological models (GHMs) into BNs, in order to probabilistically assess risks for water supply. Projected relative changes in hydrological variables computed by three GHMs driven by the output of four global climate models were processed using MATLAB, taking into account local information on water availability and use. A roadmap to set up BNs and apply probability distributions of risk levels under historic and future climate and water use was co-developed with experts from the Maghreb (Tunisia, Algeria, Morocco) who positively evaluated the BN application for LCCRAs. We conclude that the presented approach is suitable for application in the many LCCRAs necessary for successful adaptation to climate change world-wide.
At high pressures, autoionization – along with polymerization and metallization – is one of the responses of simple molecular systems to a rise in electron density. Nitrosonium nitrate (NO+NO3−), known for this property, has attracted a large interest in recent decades and was reported to be synthesized at high pressure and high temperature from a variety of nitrogen–oxygen precursors, such as N2O4, N2O and N2–O2 mixtures. However, its structure has not been determined unambiguously. Here, we present the first structure solution and refinement for nitrosonium nitrate on the basis of single-crystal X-ray diffraction at 7.0 and 37.0 GPa. The structure model (P21/m space group) contains the triple-bonded NO+ cation and the NO3− sp2-trigonal planar anion. Remarkably, crystal-chemical considerations and accompanying density-functional-theory calculations show that the oxygen atom of the NO+ unit is positively charged – a rare occurrence when in the presence of a less-electronegative element.
An update of the two-energy turbulence scheme is presented, the 2TE + APDF scheme. The original version of the two-energy scheme is able to successfully model shallow convection without the need of an additional parameterization for non-local fluxes. However, the performance of the two-energy scheme is worse in stratocumulus cases, where it tends to overestimate the erosion of the stable layers. We have identified the causes: the non-local stability parameter does not consider local stratification, the scheme lacks an internal parameter that could distinguish between a shallow convection regime and a stratocumulus regime, and it uses an inflexible turbulence length scale formulation. To alleviate this problem, we propose several modifications: an update of the stability parameter, a modified computation of the turbulence length scale, and the introduction of the entropy potential temperature to distinguish between a shallow convection and a stratocumulus regime. In addition, the two-energy scheme is coupled to a simplified assumed probability density function method in order to achieve a more universal representation of the cloudy regimes. The updated turbulence scheme is evaluated for several idealized cases and one selected real case in the ICOsahedral Nonhydrostatic (ICON) modeling framework. The results show that the updated scheme corrects the overmixing problem in the stratocumulus cases. The performance of the updated scheme is comparable to the operational setup, and can be thus used instead of the operational turbulence and shallow convection scheme in ICON. Additionally, the updated scheme improves the coupling with dynamics, which is beneficial for the modeling of coherent flow structures in the atmospheric boundary layer.
The accessories perovskite, pyrochlore, zirconolite, calzirtite and melanite from carbonatites and carbonate-rich foidites from the Kaiserstuhl are variously suited for the in situ determination of their U–Pb ages and Sr, Nd- and Hf-isotope ratios by LA-ICP-MS. The 143Nd/144Nd ratios may be determined precisely in all five phases, the 176Hf/177Hf ratios only in calzirtite and the 87Sr/86Sr ratios in perovskites and pyrochlores. The carbonatites and carbonate-rich foidites belong to one of the three magmatic groups that Schleicher et al. (1990) distinguished in the Kaiserstuhl on the basis of their Sr, Nd and Pb isotope ratios. Tephrites, phonolites and essexites (nepheline monzogabbros) form the second and limburgites (nepheline basanites) and olivine nephelinites the third. Our 87Sr/86Sr isotope data from the accessories overlap with the carbonatite and olivine nephelinite fields defined by Schleicher et al. (1990) but exhibit a much narrower range. These and the εNd and εHf values plot along the mantle array in the field of oceanic island basalts relatively close to mid-ocean ridge basalts. Previously reported K–Ar, Ar–Ar and fission track ages for the Kaiserstuhl lie between 16.2 and 17.8 Ma. They stem entirely from the geologically older tephrites, phonolites and essexites. No ages existed so far for the geologically younger carbonatites and carbonate-rich foidites except for one apatite fission track age (15.8 Ma). We obtained precise U–Pb ages for zirconolites and calzirtites of 15.66, respectively 15.5 Ma (± 0.1 2σ) and for pyrochlores of 15.35 ± 0.24 Ma. Only the perovskites from the Badberg soevite yielded a U–P concordia age of 14.56 ± 0.86 Ma while the perovskites from bergalites (haüyne melilitites) only gave 206Pb/238U and 208Pb/232Th ages of 15.26 ± 0.21, respectively, 15.28 ± 0.48 Ma. The main Kaiserstuhl rock types were emplaced over a time span of 1.6 Ma almost 1 million years before the carbonatites and carbonate-rich foidites. These were emplaced within only 0.32 Ma.
Living on the edge: environmental variability of a shallow late Holocene cold-water coral mound
(2022)
Similar to their tropical counterparts, cold-water corals (CWCs) are able to build large three-dimensional reef structures. These unique ecosystems are at risk due to ongoing climate change. In particular, ocean warming, ocean acidification and changes in the hydrological cycle may jeopardize the existence of CWCs. In order to predict how CWCs and their reefs or mounds will develop in the near future one important strategy is to study past fossil CWC mounds and especially shallow CWC ecosystems as they experience a greater environmental variability compared to other deep-water CWC ecosystems. We present results from a CWC mound off southern Norway. A sediment core drilled from this relatively shallow (~ 100 m) CWC mound exposes in full detail hydrographical changes during the late Holocene, which were crucial for mound build-up. We applied computed tomography, 230Th/U dating, and foraminiferal geochemical proxy reconstructions of bottom-water-temperature (Mg/Ca-based BWT), δ18O for seawater density, and the combination of both to infer salinity changes. Our results demonstrate that the CWC mound formed in the late Holocene between 4 kiloannum (ka) and 1.5 ka with an average aggradation rate of 104 cm/kiloyears (kyr), which is significantly lower than other Holocene Norwegian mounds. The reconstructed BWTMg/Ca and seawater density exhibit large variations throughout the entire period of mound formation, but are strikingly similar to modern in situ observations in the nearby Tisler Reef. We argue that BWT does not exert a primary control on CWC mound formation. Instead, strong salinity and seawater density variation throughout the entire mound sequence appears to be controlled by the interplay between the Atlantic Water (AW) inflow and the overlying, outflowing Baltic-Sea water. CWC growth and mound formation in the NE Skagerrak was supported by strong current flow, oxygen replenishment, the presence of a strong boundary layer and larval dispersal through the AW, but possibly inhibited by the influence of fresh Baltic Water during the late Holocene. Our study therefore highlights that modern shallow Norwegian CWC reefs may be particularly endangered due to changes in water-column stratification associated with increasing net precipitation caused by climate change.
We performed an experiment under long-term microgravity conditions aboard the International Space Station (ISS) to obtain information on the energetics and experimental constraints required for the formation of chondrules in the solar nebula by ’nebular lightning’. As a simplified model system, we exposed porous forsterite (Mg2 SiO4) dust particles to high-energetic arc discharges. The characterization of the samples after their return by synchrotron microtomography and scanning electron microscopy revealed that aggregates had formed, consisting of several fused Mg2SiO4 particles. The partial melting and fusing of Mg2SiO4 dust particles under microgravity conditions leads to a strong reduction of their porosity. The experimental outcomes vary strongly in their appearance from small spherical melt-droplets (∅≈90 µm) to bigger and irregularly shaped aggregates (∅≈350 µm). Our results provided new constraints with respect to energetic aspects of chondrule formation and a roadmap for future and more complex experiments on Earth and in microgravity conditions.
Several past summer floods in Central Europe were associated with so-called Vb‑cyclones propagating from the Mediterranean Sea north-eastward to Central Europe. This study illustrates the usefulness of the parametric transfer entropy measure TE‑linear in investigating heavy Vb‑cyclone precipitation events in the Odra catchment (Poland). With the application of the TE‑linear approach, we confirm the impact of the Mediterranean Sea on precipitation intensification. Moreover, we also detect significant information exchange to Vb‑cyclone precipitation from evaporation over the European continent along the typical Vb‑cyclone pathway. Thus, the Mediterranean Sea could enhance the Vb‑cyclone precipitation by pre-moistening continental moisture source regions that contribute to precipitation downstream in the investigated catchments. Overall, the transfer entropy approach with the measure TE‑linear proved to be computationally effective and complementary to traditional methods such as Lagrangian and Eulerian diagnostics.
This work describes the development and characterization of two instruments and their data evaluation, which contributes to a better understanding of new particle formation and growth, as well as their interactions with clouds. Both instruments were characterized at the Cosmics Leaving Outdoor Droplets (CLOUD) experiment at the European Center for Nuclear Research (CERN).
Climatic and biogeographical drivers of functional diversity in the flora of the Canary Islands
(2022)
Aim: Functional traits can help us to elucidate biogeographical and ecological processes driving assemblage structure. We analysed the functional diversity of plant species of different evolutionary origins across an island archipelago, along environmental gradients and across geological age, to assess functional aspects of island biogeographical theory.
Location: Canary Islands, Spain.
Major taxa studied: Spermatophytes.
Time period: Present day.
Methods: We collected data for four traits (plant height, leaf length, flower length and fruit length) associated with resource acquisition, competitive ability, reproduction and dispersal ability of 893 endemic, non-endemic native and alien plant species (c. 43% of the Canary Island flora) from the literature. Linking these traits to species occurrences and composition across a 500 m × 500 m grid, we calculated functional diversity for endemic, non-endemic native and alien assemblages using multidimensional functional hypervolumes and related the resulting patterns to climatic (humidity) and island biogeographical (geographical isolation, topographic complexity and geological age) gradients.
Results: Trait space of endemic and non-endemic native species overlapped considerably, and alien species added novel trait combinations, expanding the overall functional space of the Canary Islands. We found that functional diversity of endemic plant assemblages was highest in geographically isolated and humid grid cells. Functional diversity of non-endemic native assemblages was highest in less isolated and humid grid cells. In contrast, functional diversity of alien assemblages was highest in arid ecosystems. Topographic complexity and geological age had only a subordinate effect on functional diversity across floristic groups.
Main conclusions: We found that endemic and non-endemic native island species possess similar traits, whereas alien species tend to expand functional space in ecosystems where they have been introduced. The spatial distribution of the functional diversity of floristic groups is very distinct across environmental gradients, indicating that species assemblages of different evolutionary origins thrive functionally in dissimilar habitats.
In order to gain further insights into early solar system aggregation processes, we carried out an experiment on board the International Space Station, which allowed us to study the behavior of dust particles exposed to electric arc discharges under long-term microgravity. The experiment led to the formation of robust, elongated, fluffy aggregates, which were studied by scanning electron microscopy, electron backscatter diffraction, and synchrotron micro-computed tomography. The morphologies of these aggregates strongly resemble the typical shapes of fractal fluffy-type calcium-aluminum-rich inclusions (CAIs). We conclude that a small amount of melting could have supplied the required stability for such fractal structures to have survived transportation and aggregation to and compaction within planetesimals. Other aggregates produced in our experiment have a massy morphology and contain relict grains, likely resulting from the collision of grains with different degrees of melting, also observed in some natural CAIs. Some particles are surrounded by igneous rims, which remind in thickness and crystal orientation of Wark–Lovering rims; another aggregate shows similarities to disk-shaped CAIs. These results imply that a (flash-)heating event with subsequent aggregation could have been involved in the formation of different morphological CAI characteristics.
Convective rain cell properties and the resulting precipitation scaling in a warm-temperate climate
(2022)
Convective precipitation events have been shown to intensify at rates exceeding the Clausius–Clapeyron rate (CC rate) of ca. 7% K−1 under current climate conditions. In this study, we relate atmospheric variables (low-level dew point temperature, convective available potential energy, and vertical wind shear), which are regarded as ingredients for severe deep convection, to properties of convective rain cells (cell area, maximum precipitation intensity, lifetime, precipitation sum, and cell speed). The rain cell properties are obtained from a rain gauge-adjusted radar dataset in a mid-latitude region, which is characterized by a temperate climate with warm summers (Germany). Different Lagrangian cell properties scale with dew point temperature at varying rates. While the maximum precipitation intensity of cells scales consistently at the CC rate, the area and precipitation sum per cell scale at varying rates above the CC rate. We show that this super-CC scaling is caused by a covarying increase of convective available potential energy with dew point temperature. Wind shear increases the precipitation sum per cell mainly by increasing the spatial cell extent. From a Eulerian point of view, this increase is partly compensated by a higher cell velocity, which leads to Eulerian precipitation scaling rates close to and slightly above the CC rate. Thus, Eulerian scaling rates of convective precipitation are modulated by convective available potential energy and vertical wind shear, making it unlikely that present scaling rates can be applied to future climate conditions. Furthermore, we show that cells that cause heavy precipitation at fixed locations occur at low vertical wind shear and, thus, move relatively slowly compared to typical cells.
Moisture sources of heavy precipitation in Central Europe in synoptic situations with Vb-cyclones
(2022)
During the past century, several extreme summer floods in Central Europe were associated with so-called Vb-cyclones propagating from the Mediterranean Sea north-eastward to Central Europe. The processes intensifying the precipitation in synoptic situations with Vb-cyclones in the Danube, Elbe, and Odra catchments are only partially understood. Our study aims to investigate these processes with Lagrangian moisture-source diagnostics for 16 selected Vb-events. Moreover, we analyse the characteristics of typical moisture source regions during 1107 Vb-events from 1901 to 2010 based on ERA-20C reanalysis dynamically downscaled with COSMO-CLM+NEMO. We observe moisture contributions by various source regions highlighting the complex dynamical interplay of different air masses leading to moisture convergence in synoptic situations with Vb-cyclones. Overall, up to 80% of the precipitation originates from the European continent, indicating the importance of continental moisture recycling, especially within the respective river catchment. Other major moisture uptake regions are the North Sea, the Baltic Sea, the North Atlantic, and for a few events the Black Sea. Remarkably, anomalies in these oceanic source regions show no connection to precipitation amounts in synoptic situations with Vb-cyclones. In contrast, the Vb-cyclones with the highest precipitation are associated with anomalously high evaporation in the Mediterranean Sea, even though the Mediterranean Sea is only a minor moisture source region on average. Interestingly, the evaporation anomalies are not connected with sea-surface temperature but with wind-speed anomalies (Spearman’s rank correlation coefficient R≈0.7, significant with p<0.01) indicating mainly dynamically driven evaporation. The particular role of the Mediterranean Sea hints towards possible importance of Mediterranean moisture for the early-stage intensification of Vb-cyclones and the pre-moistening of the continental uptake regions upstream of the target catchments.
The Yacoraite Formation (Salta rift, Argentina) consists of Maastrichtian–Danian lacustrine carbonate and siliciclastic deposits with interbedded volcanic ash layers, organized in four third-order stratigraphic sequences. It offers the exceptional opportunity to jointly apply in situ zircon and carbonate U-Pb geochronology that resulted in two distinct depositional age depth models. Ages of the youngest zircon population from ash layers were linearly interpolated to derive a zircon depositional age depth model. A carbonate depositional age depth model was instead obtained from dated carbonate phases including microbialites, ooids, oncoids of calcitic and dolomitic mineralogy as well as early lacustrine calcite cements. Mean ages were defined from different carbonate phases belonging to the same layer and then linearly interpolated. Sedimentation rates were calculated from both depth models between pairs of dated samples and used to estimate the age of sequence boundaries, as well as the duration of the four stratigraphic sequences. The zircon and carbonate depositional age depth models agree with biostratigraphic constraints and exhibit excellent consistency. The onset and end of sedimentation were estimated at 68.2 ± 0.9 Ma and 62.3 ± 0.6 Ma (duration ca 5.7 Ma) via zircon geochronology and at 67.9 ± 1.7 Ma and 61.9 ± 1.3 Ma (duration ca 6.0 Ma) via carbonate geochronology. Results from this study show that with suitable samples and a newly implemented working strategy, in situ U-Pb dating of depositional and early diagenetic carbonates represent a valuable chronostratigraphic tool for estimating sedimentation rate and duration in poorly time-framed depositional systems.