550 Geowissenschaften
Refine
Year of publication
Document Type
- Article (888)
- 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 (775)
- 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)
Ob wir bei der Erwärmung eher in Richtung plus 2 Grad oder eher in Richtung plus 5 Grad steuern, entscheidet sich zu einem nicht unbedeutenden Teil über den Meeren. Dort beschatten niedrige, flache, sehr großflächige Wolkenschichten die darunterliegende Wasseroberfläche und sorgen so für Abkühlung. Auch in Zukunft noch?
The current state of research about ancient settlements within the Nile Delta allows the hypothesizing of fluvial connections to ancient settlements all over the Nile Delta. Previous studies suggest a larger Nile branch close to Kom el-Gir, an ancient settlement hill in the northwestern Nile Delta. To contribute new knowledge to this little-known site and prove this hypothesis, this study aims at using small-scale paleogeographic investigations to reconstruct an ancient channel system in the surroundings of Kom el-Gir. The study pursues the following: (1) the identification of sedimentary environments via stratigraphic and portable X-ray fluorescence (pXRF) analyses of the sediments, (2) the detection of fluvial elements via electrical resistivity tomography (ERT), and (3) the synthesis of all results to provide a comprehensive reconstruction of a former fluvial network in the surroundings of Kom el-Gir. Therefore, auger core drillings, pXRF analyses, and ERT were conducted to examine the sediments within the study area. Based on the evaluation of the results, the study presents clear evidence of a former channel system in the surroundings of Kom el-Gir. Thereby, it is the combination of both methods, 1-D corings and 2-D ERT profiles, that derives a more detailed illustration of previous environmental conditions which other studies can adopt. Especially within the Nile Delta which comprises a large number of smaller and larger ancient settlement hills, this study's approach can contribute to paleogeographic investigations to improve the general understanding of the former fluvial landscape.
Previous investigation of seismic anisotropy indicates the presence of a simple mantle flow regime beneath the Turkish-Anatolian Plateau and Arabian Plate. Numerical modeling suggests that this simple flow is a component of a large-scale global mantle flow associated with the African superplume, which plays a key role in the geodynamic framework of the Arabia-Eurasia continental collision zone. However, the extent and impact of the flow pattern farther east beneath the Iranian Plateau and Zagros remains unclear. While the relatively smoothly varying lithospheric thickness beneath the Anatolian Plateau and Arabian Plate allows progress of the simple mantle flow, the variable lithospheric thickness across the Iranian Plateau is expected to impose additional boundary conditions on the mantle flow field. In this study, for the first time, we use an unprecedented data set of seismic waveforms from a network of 245 seismic stations to examine the mantle flow pattern and lithospheric deformation over the entire region of the Iranian Plateau and Zagros by investigation of seismic anisotropy. We also examine the correlation between the pattern of seismic anisotropy, plate motion using GPS velocities and surface strain fields. Our study reveals a complex pattern of seismic anisotropy that implies a similarly complex mantle flow field. The pattern of seismic anisotropy suggests that the regional simple mantle flow beneath the Arabian Platform and eastern Turkey deflects as a circular flow around the thick Zagros lithosphere. This circular flow merges into a toroidal component beneath the NW Zagros that is likely an indicator of a lateral discontinuity in the lithosphere. Our examination also suggests that the main lithospheric deformation in the Zagros occurs as an axial shortening across the belt, whereas in the eastern Alborz and Kopeh-Dagh a belt-parallel horizontal lithospheric deformation plays a major role.
Atmosphärische Schwerewellen spielen eine wichtige Rolle für die Zirkulation der mittleren Atmosphäre, die wiederum die Troposphäre auf saisonalen und längeren Zeitskalen beeinflusst, und stellen somit ein Schlüsselelement für das Wetter- und Klimageschehen dar. Eine adäquate Beschreibung des Lebenszyklus atmosphärischer Schwerewellen in den operationellen Modellen zur Wettervorhersage und Klimasimulation ist daher sehr wünschenswert. Um zu einer verbesserten mathematischen Darstellung der Schwerewellendynamik in den Modellen beizutragen, wurden in den vergangenen Jahren zahlreiche numerische Studien durchgeführt. Wenngleich auch viele der ablaufenden Prozesse gegenwärtig gut verstanden sind, stellt die Wechselwirkung zwischen den mesoskaligen Schwerewellen und den synoptischskaligen Prozessen aufgrund der hohen Komplexität der Strömung weiterhin eine besondere Herausforderung für die Erforschung der Schwerewellenaktivität dar und erfordert oftmals hochaufgelöste numerische Simulationen über große Modelldomänen.
Folglich ist es wichtig, dass die angewendeten numerischen Verfahren effizient sind und möglichst idealisierte, aber dennoch atmosphärenähnliche Szenarien simulieren. In dieser Arbeit wird ein effizientes numerisches Verfahren zur Modellierung der Dynamik interner Schwerewellen sowie deren Einfluss auf die Zirkulation der mittleren Atmosphäre entwickelt.
Dabei wird die Diskretisierung des pseudo-inkompressiblen Finite-Volumen-Modells auf einem versetzten Gitter von Rieper et al. (2013), welches der Einfachheit halber Schallwellen aus der Dynamik herausfiltert und zur Untersuchung adiabatischer Atmosphärenprozesse auf der f-Ebene entwickelt wurde, im wesentlichen durch zwei Komponenten erweitert: 1) die Anwendung eines semi-impliziten Zeitschrittverfahrens auf die Bewegungsgleichungen zur Integration der Auftriebs- und Corioliseffekte und 2) die Berücksichtigung einer Heizung durch einen thermischen Relaxationsansatz, welcher in der Troposphäre ein baroklin instabiles Strömungsprofil erzeugt und eine zeitabhängige Dynamik des Hintergrundzustands zulässt. Zur Überprüfung der korrekten Implementierung der Erweiterungen werden eine Reihe von atmosphärischen Standardteststudien durchgeführt, welche die Konvergenzeigenschaften sowie die Effizienz des Verfahrens validieren. Darüber hinaus zeigen die Testfälle, dass die Ergebnisse des Modells mit anderen veröffentlichten Arbeiten sehr gut übereinstimmen.
Schließlich wird als Anwendungstestfall eine mesoskalige Simulation barokliner Instabilität in der Troposphäre durchgeführt, welche ferner die darin enthaltene kleinskalige Wellenaktivität sowie deren Einfluss auf die mittlere Atmosphäre modelliert. Die abschließende Betrachung der zonal und zeitlich gemittelten Felder zeigt die erwartete Zonalwindumkehr in der Höhe.
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.
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.
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.
The formation of terrestrial planets was a complex process which begun in the very early stage of the Solar System in the protoplanetary disk (PPD). Chondrites are fragments of planet precursors, which have never experienced differentiation and can help to reconstruct the first processes leading to planet formation. The main components of chondrites are chondrules, calcium-aluminum-rich inclusions (CAIs), amoeboid olivine aggregates (AOAs), metals and fine-grained material. Each of these components formed by a complex mechanism involving aggregation and/or melting. Previous research has already provided an overall view of the formation of these objects, however, there are still open questions regarding the aggregation behavior of particles, the heating mechanism(s) and the thermal history of CAIs, AOAs and chondrules. For instance, the involvement of flash-heating events and electrostatics in the aggregation and melting of these objects has been a keen topic of discussion.
The aim of this doctoral thesis was to develop and carry out an experiment to study various early Solar System processes under long-term microgravity. In the project with the acronym EXCISS (Experimental Chondrule Formation aboard the ISS), free-floating, 126(23)µm-sized Mg2SiO4 dust particles were exposed to electric fields and electric discharges.
The experimental set-up was installed inside a 10x10x15 cm3-sized container and consisted of an arc generation unit connected to the sample chamber, a camera with an optical system, a power supply unit with lithium-ion batteries and the EXCISS mainboard with a Raspberry Pi Zero and mass storage devices. The sample chamber was manufactured from quartz glass and the experiments were filmed. The complete experiment container was subsequently returned to the Goethe University and the samples were analyzed with scanning electron microscopy, electron backscatter diffraction and synchrotron micro-CT.
Video analysis has shown that particles, which were agitated by electric discharges, align in chains within the electric field with their longest axis parallel to the electric field lines. Consequently, electric fields could have influenced the inner structure and porosity of particle aggregates in the PPD.
The discharge experiments produced fused aggregates and individual melt spherules.
The fused aggregates share many morphological characteristics with natural fluffy-type CAIs and some igneous CAIs found in chondrites. Consequently, CAIs could have formed by the aggregation of particles with various degrees of melting. Further, a small amount of melting could have supplied the required stability for such fractal structures to have survived transportation and aggregation to, and subsequent compaction within, developing planetesimals.
Some initial particles were completely melted by the arc discharges and formed melt spherules. The newly formed olivines crystallized with a preferred orientation of the [010] axis perpendicular to the surface of the spherule. Similar preferred orientations have been found in natural chondrules. However, the microstructure differs from the results of previous experiments on Earth, which show, for example, crystal settling on one side of the sample because of the influence of gravity. Furthermore, the melt spherules show evidence for an interaction of the melt with the surrounding hot gas. Therefore, microgravity experiments with more advanced experimental parameters bear great potential for future chondrule formation experiments.
Responses of southern ocean seafloor habitats and communities to global and local drivers of change
(2021)
Knowledge of life on the Southern Ocean seafloor has substantially grown since the beginning of this century with increasing ship-based surveys and regular monitoring sites, new technologies and greatly enhanced data sharing. However, seafloor habitats and their communities exhibit high spatial variability and heterogeneity that challenges the way in which we assess the state of the Southern Ocean benthos on larger scales. The Antarctic shelf is rich in diversity compared with deeper water areas, important for storing carbon (“blue carbon”) and provides habitat for commercial fish species. In this paper, we focus on the seafloor habitats of the Antarctic shelf, which are vulnerable to drivers of change including increasing ocean temperatures, iceberg scour, sea ice melt, ocean acidification, fishing pressures, pollution and non-indigenous species. Some of the most vulnerable areas include the West Antarctic Peninsula, which is experiencing rapid regional warming and increased iceberg-scouring, subantarctic islands and tourist destinations where human activities and environmental conditions increase the potential for the establishment of non-indigenous species and active fishing areas around South Georgia, Heard and MacDonald Islands. Vulnerable species include those in areas of regional warming with low thermal tolerance, calcifying species susceptible to increasing ocean acidity as well as slow-growing habitat-forming species that can be damaged by fishing gears e.g., sponges, bryozoan, and coral species. Management regimes can protect seafloor habitats and key species from fishing activities; some areas will need more protection than others, accounting for specific traits that make species vulnerable, slow growing and long-lived species, restricted locations with optimum physiological conditions and available food, and restricted distributions of rare species. Ecosystem-based management practices and long-term, highly protected areas may be the most effective tools in the preservation of vulnerable seafloor habitats. Here, we focus on outlining seafloor responses to drivers of change observed to date and projections for the future. We discuss the need for action to preserve seafloor habitats under climate change, fishing pressures and other anthropogenic impacts.
This paper investigates the global stratospheric Brewer–Dobson circulation (BDC) in the ERA5 meteorological reanalysis from the European Centre for Medium-Range Weather Forecasts (ECMWF). The analysis is based on simulations of stratospheric mean age of air, including the full age spectrum, with the Lagrangian transport model CLaMS (Chemical Lagrangian Model of the Stratosphere), driven by reanalysis winds and total diabatic heating rates. ERA5-based results are compared to results based on the preceding ERA-Interim reanalysis. Our results show a significantly slower BDC for ERA5 than for ERA-Interim, manifesting in weaker diabatic heating rates and higher age of air. In the tropical lower stratosphere, heating rates are 30 %–40 % weaker in ERA5, likely correcting a bias in ERA-Interim. At 20 km and in the Northern Hemisphere (NH) stratosphere, ERA5 age values are around the upper margin of the uncertainty range from historical tracer observations, indicating a somewhat slow–biased BDC. The age trend in ERA5 over the 1989–2018 period is negative throughout the stratosphere, as climate models predict in response to global warming. However, the age decrease is not linear but steplike, potentially caused by multi-annual variability or changes in the observations included in the assimilation. During the 2002–2012 period, the ERA5 age shows a similar hemispheric dipole trend pattern as ERA-Interim, with age increasing in the NH and decreasing in the Southern Hemisphere (SH). Shifts in the age spectrum peak and residual circulation transit times indicate that reanalysis differences in age are likely caused by differences in the residual circulation. In particular, the shallow BDC branch accelerates in both reanalyses, whereas the deep branch accelerates in ERA5 and decelerates in ERA-Interim.
Extreme convective precipitation is expected to increase with global warming. However, the rate of increase and the understanding of contributing processes remain highly uncertain. We investigated characteristics of convective rain cells like area, intensity, and lifetime as simulated by a convection-permitting climate model in the area of Germany under historical (1976–2005) and future (end-of-century, RCP8.5 scenario) conditions. To this end, a tracking algorithm was applied to 5-min precipitation output. While the number of convective cells is virtually similar under historical and future conditions, there are more intense and larger cells in the future. This yields an increase in hourly precipitation extremes, although mean precipitation decreases. The relative change in the frequency distributions of area, intensity, and precipitation sum per cell is highest for the most extreme percentiles, suggesting that extreme events intensify the most. Furthermore, we investigated the temperature and moisture scaling of cell characteristics. The temperature scaling drops off at high temperatures, with a shift in drop-off towards higher temperatures in the future, allowing for higher peak values. In contrast, dew point temperature scaling shows consistent rates across the whole dew point range. Cell characteristics scale at varying rates, either below (mean intensity), at about (maximum intensity and area), or above (precipitation sum) the Clausius–Clapeyron rate. Thus, the widely investigated extreme precipitation scaling at fixed locations is a complex product of the scaling of different cell characteristics. The dew point scaling rates and absolute values of the scaling curves in historical and future conditions are closest for the highest percentiles. Therefore, near-surface humidity provides a good predictor for the upper limit of for example, maximum intensity and total precipitation of individual convective cells. However, the frequency distribution of the number of cells depending on dew point temperature changes in the future, preventing statistical inference of extreme precipitation from near-surface humidity.
The reanalysis products and derived products, ERA5 (Copernicus Climate Change Service, 2018) and W5E5 (WATCH Forcing Data (WFD) methodology applied to ERA5) (LANGE ET AL., 2021) have been recently published initiating a new phase of scientific research utilizing these datasets. ERA5 and W5E5 offer the possibility to reduce insecurities in model results through their improved quality compared to previous climate reanalyses (CUCCHI ET AL., 2020). The suitability of either climate forcing as input for the hydrological model WaterGAP and the influence of the models specific calibration routine has been evaluated with four model experiments. The model was validated by analysing the models ability to produce reasonable values for global water balance components and to reproduce observed discharge in 1427 basins as well as total water storage anomalies in 143 basins using well established efficiency metrics. Bias correction of W5E5 was found to lead to more global realistic mean precipitation and consequently discharge and AET values. In an uncalibrated model setup ERA5 results in better performances across all efficiency metrics. Model results produced with W5E5 as climate input were strongly improved through calibration ultimately leading to the best performances out of all four model experiments. However, model performances considerably improved through calibration with both climate forcings hence calibration was found to have the strongest effect on model performance. Furthermore, spatial differences in performance of either forcing were identified. Snow-dominated regions show an overall better performance with ERA5, while wetter and warmer regions are better represented with W5E5. Finally, it can be concluded that W5E5 should be preferred as climate input for impact modelling; however, depending on the spatial scale and region ERA5 should at least be considered, in particular for snow-dominated regions.
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.
We present the results of a multi-disciplinary investigation on a deciduous human tooth (Pradis 1), recently recovered from the Epigravettian layers of the Grotte di Pradis archaeological site (Northeastern Italian Prealps). Pradis 1 is an exfoliated deciduous molar (Rdm2), lost during life by an 11–12-year-old child. A direct radiocarbon date provided an age of 13,088–12,897 cal BP (95% probability, IntCal20). Amelogenin peptides extracted from tooth enamel and analysed through LC–MS/MS indicate that Pradis 1 likely belonged to a male. Time-resolved 87Sr/86Sr analyses by laser ablation mass spectrometry (LA-MC-ICPMS), combined with dental histology, were able to resolve his movements during the first year of life (i.e. the enamel mineralization interval). Specifically, the Sr isotope ratio of the tooth enamel differs from the local baseline value, suggesting that the child likely spent his first year of life far from Grotte di Pradis. Sr isotopes are also suggestive of a cyclical/seasonal mobility pattern exploited by the Epigravettian human group. The exploitation of Grotte di Pradis on a seasonal, i.e. summer, basis is also indicated by the faunal spectra. Indeed, the nearly 100% occurrence of marmot remains in the entire archaeozoological collection indicates the use of Pradis as a specialized marmot hunting or butchering site. This work represents the first direct assessment of sub-annual movements observed in an Epigravettian hunter-gatherer group from Northern Italy.
Wildfire is the most common disturbance type in boreal forests and can trigger significant changes in forest composition. Waterlogging in peatlands determines the degree of tree cover and the depth of the burnt horizon associated with wildfires. However, interactions between peatland moisture, vegetation composition and flammability, and fire regime in forest and forested peatland in Eurasia remain largely unexplored, despite their huge extent in boreal regions. To address this knowledge gap, we reconstructed the Holocene fire regime, vegetation composition, and peatland hydrology at two sites located in predominantly light taiga (Pinus sylvestris Betula) with interspersed dark taiga communities (Pinus sibirica, Picea obovata, Abies sibirica) in western Siberia in the Tomsk Oblast, Russia. We found marked shifts in past water levels over the Holocene. The probability of fire occurrence and the intensification of fire frequency and severity increased at times of low water table (drier conditions), enhanced fuel dryness, and an intermediate dark-to-light taiga ratio. High water level, and thus wet peat surface conditions, prevented fires from spreading on peatland and surrounding forests. Deciduous trees (i.e. Betula) and Sphagnum were more abundant under wetter peatland conditions, and conifers and denser forests were more prevalent under drier peatland conditions. On a Holocene scale, severe fires were recorded between 7.5 and 4.5 ka with an increased proportion of dark taiga and fire avoiders (Pinus sibirica at Rybnaya and Abies sibirica at Ulukh–Chayakh) in a predominantly light taiga and fire-resister community characterised by Pinus sylvestris and lower local water level. Severe fires also occurred over the last 1.5 kyr and were associated with a declining abundance of dark taiga and fire avoiders, an expansion of fire invaders (Betula), and fluctuating water tables. These findings suggest that frequent, high-severity fires can lead to compositional and structural changes in forests when trees fail to reach reproductive maturity between fire events or where extensive forest gaps limit seed dispersal. This study also shows prolonged periods of synchronous fire activity across the sites, particularly during the early to mid-Holocene, suggesting a regional imprint of centennial- to millennial-scale Holocene climate variability on wildfire activity. Humans may have affected vegetation and fire from the Neolithic; however, increasing human presence in the region, particularly at the Ulukh–Chayakh Mire over the last 4 centuries, drastically enhanced ignitions compared to natural background levels. Frequent warm and dry spells predicted by climate change scenarios for Siberia in the future will enhance peatland drying and may convey a competitive advantage to conifer taxa. However, dry conditions will probably exacerbate the frequency and severity of wildfire, disrupt conifers' successional pathway, and accelerate shifts towards deciduous broadleaf tree cover. Furthermore, climate–disturbance–fire feedbacks will accelerate changes in the carbon balance of boreal peatlands and affect their overall future resilience to climate change.
Eisbildende Prozesse sind für die Wolkenbildung von großer Bedeutung und haben erhebliche Auswirkungen auf das Wetter und Klima der Erde, indem sie den Strahlungsantrieb und die Niederschlagsbildung beeinflussen. In den mittleren Breiten entsteht der meiste Niederschlag in sogenannten Mischphasenwolken (MPC), welche sowohl aus unterkühlten Wolkentröpfchen als auch aus Eiskristallen bestehen. Bei Temperaturen zwischen 0°C und -38°C erfolgt die Bildung von Eiskristallen in MPC in Gegenwart von Aerosolpartikeln, die als sogenannte Eiskeime (INP) die Fähigkeit besitzen, auf ihrer Oberfläche Eis zu nukleieren. Trotz der großen wissenschaftlichen Fortschritte in den letzten Jahrzehnten, weist der heterogene Eisbildungsprozess, als einer der wichtigsten in der Atmosphäre auftretenden Aerosol-Wolken-Wechselwirkungsprozesse, immer noch große Unsicherheiten auf. Um zukünftige Klimavorhersagen und -projektionen in Modellen besser abbilden zu können, ist es somit notwendig den Wissensgrad der räumlichen und zeitlichen Heterogenität von INP in Bezug auf Herkunft, Anzahl und Zusammensetzung zu erhöhen. Im Zentrum dieser Arbeit steht der Eiskeimzähler FINCH (Fast Ice Nucleus Chamber), der für Labor- und Feldexperimente von der Johann Wolfgang Goethe-Universität in Frankfurt am Main entwickelt wurde. Durch das Mischen des Probenstroms mit einem warm-feuchten und einem kalten-trockenen Luftstrom wird eine Übersättigung in der in-situ Eiskammer erreicht, die benötigt wird, eisbildende Partikel zu aktivieren. Die aktivierten Partikel können beim Durchströmen der Kammer zu Wassertropfen oder Eiskristallen anwachsen. Am Ausgang der Kammer wird die Anzahl und Größe der Partikel durch die FINCH-Optik erfasst. Als grundlegender Schritt und aufbauend auf den Charakterisierungsmessungen von Frank (2017) wurden in der vorliegenden Arbeit die Leistung, die Zuverlässigkeit sowie die Reproduzierbarkeit von FINCH in Validierungsexperimenten im Labor überprüft. Im Zuge dessen wurden heterogene Gefrierexperimente mit definierten Referenzaerosolproben (bspw. K-Feldspat) bei wasserübersättigten Bedingungen und verschiedenen Gefriertemperaturen durchgeführt. Für den Großteil der erzielten Resultate konnte eine zufriedenstellende Übereinstimmung mit Literaturwerten von anderen INP-Messinstrumenten aus der ganzen Welt erzielt werden. Es zeigte sich, dass die Leistungsfähigkeit von FINCH messtechnische Limitationen für Messexperimente bei Temperaturen >-10°C und <-30°C aufweist, was eine Einschränkung des Messbereichs bedeutet. Hinsichtlich der Quantifizierung des Unsicherheitsbereiches des Messgerätes in Bezug auf Temperatur und relativer Feuchte bedarf es im Nachgang an dieser Arbeit weiterer Charakterisierungsmessungen. Im Rahmen der Ice Nuclei Research Unit (INUIT) Forschergruppe wurde FINCH mit einem gepumpten Gegenstrom-Impaktor PCVI und dem online Einzelpartikel-Massenspektrometer ALABAMA gekoppelt. Diese spezielle Messmethodik dient zur chemischen und mikrophysikalischen Charakterisierung der INP und der Eispartikelresiduen (IPR). Der Fokus lag zunächst darauf die Funktionalität des gekoppelten Messsystems im Labor zu überprüfen. Ausführliche Charakterisierungsmessungen zeigten unter eisübersättigten und unterkühlten Bedingungen, dass das Prinzip der Trennung der INP von nicht-aktivierten Aerosolen und unterkühlten Tropfen hinter FINCH durch den PCVI funktioniert. Ebenso konnten erste quantitative Aussagen zur chemischen Zusammensetzung der IPR getroffen werden. Es zeigte sich, dass bei den Aktivierungsexperimenten ein geringer Anteil an Partikeltypen metallischer Art von ALABAMA detektiert wurden, der nicht dem untersuchten Aerosoltyp zugeordnet werden konnte. Der Ursprung dieser Kontamination konnte im Rahmen dieser Arbeit nicht abschließend geklärt werden und bedarf weiterer Validierungsmessungen im Labor. Atmosphärische Eiskeimkonzentrationen wurden im Rahmen von Feldmesskampagnen an der Hochalpinen Forschungsstation Jungfraujoch (JFJ) in den Schweizer Alpen und am Campus Riedberg der Johann Wolfgang Goethe-Universität in Frankfurt am Main untersucht. Hier konnten erste Erfahrungen mit Außenluftmessungen bezüglich der Leistungsfähigkeit und der Nachweisgrenze (LOD) des Messgerätes gesammelt werden. Durch den Einfluss der freien Troposphäre am JFJ waren die Messungen hauptsächlich von aerosolpartikelarmer Luft mit einer geringen Anzahl von Eiskeimen geprägt, so dass sich die gemessenen INP-Konzentrationen oftmals unter die Nachweisgrenze von FINCH fielen. Unter Einsatz eines Aerosolkonzentrators konnte die Detektionseffizienz verbessert und das LOD herabgesetzt werden. Am JFJ wurden die INP im Mittel bei einer Temperatur von -23°C und einem Wassersättigungsverhältnis von 107% beprobt. Die mediane (mittlere) INP-Konzentration inklusive LOD lag bei 2,1 (3,3) sL-1 und oberhalb des LOD bei 3,1 (4,5) sL-1. Ein Vergleich mit den Messungen am Campus Riedberg unter annähernd gleichen Bedingungen resultiert in ähnlichen Konzentrationen.
Following votes in the Coniacian Working Group, the Cretaceous Subcommission and the International Commission on Stratigraphy, on May 1st, 2021, the International Union of Geological Sciences voted unanimously to ratify the Global Stratotype Section and Point (GSSP) proposal for the base of the Coniacian Stage of the Upper Cretaceous Series and Cretaceous System. The lower boundary of the Coniacian Stage is placed at the base of Bed 46 of the Salzgitter-Salder section in northern Germany. The boundary is defined by the first appearance of the inoceramid bivalve species Cremnoceramus deformis erectus (Meek) and complemented by the Navigation carbon isotope event. Additional data include the bivalve genus Didymotis, foraminifera, ammonite, nannofossil and organic-walled dinoflagellate cyst events. Three auxiliary sections (Słupia Nadbrzeżna, central Poland; Střeleč, Czech Republic; El Rosario, NE Mexico) supplement the details of the boundary record in various facies, and in differing geographic and biogeographic contexts.
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.
Africa's protected areas (PAs) are the last stronghold of the continent's unique biodiversity, but they appear increasingly threatened by climate change, substantial human population growth, and land-use change. Conservation planning is challenged by uncertainty about how strongly and where these drivers will interact over the next few decades. We investigated the combined future impacts of climate-driven vegetation changes inside African PAs and human population densities and land use in their surroundings for 2 scenarios until the end of the 21st century. We used the following 2 combinations of the shared socioeconomic pathways (SSPs) and representative greenhouse gas concentration pathways (RCPs): the “middle-of-the-road” scenario SSP2–RCP4.5 and the resource-intensive “fossil-fueled development” scenario SSP5–RCP8.5. Climate change impacts on tree cover and biome type (i.e., desert, grassland, savanna, and forest) were simulated with the adaptive dynamic global vegetation model (aDGVM). Under both scenarios, most PAs were adversely affected by at least 1 of the drivers, but the co-occurrence of drivers was largely region and scenario specific. The aDGVM projections suggest considerable climate-driven tree cover increases in PAs in today's grasslands and savannas. For PAs in West Africa, the analyses revealed climate-driven vegetation changes combined with hotspots of high future population and land-use pressure. Except for many PAs in North Africa, future decreases in population and land-use pressures were rare. At the continental scale, SSP5–RCP8.5 led to higher climate-driven changes in tree cover and higher land-use pressure, whereas SSP2–RCP4.5 was characterized by higher future population pressure. Both SSP–RCP scenarios implied increasing challenges for conserving Africa's biodiversity in PAs. Our findings underline the importance of developing and implementing region-specific conservation responses. Strong mitigation of future climate change and equitable development scenarios would reduce ecosystem impacts and sustain the effectiveness of conservation in Africa.
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.
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).
Questions: Both species turnover and intraspecific trait variation can affect plant assemblage dynamics along environmental gradients. Here, we asked how community assemblage patterns in relation to species turnover and intraspecific variation differ between endemic and non-endemic species. We hypothesized that endemic species show lower intraspecific variation than non-endemic species because they tend to have high rates of in situ speciation, whereas non-endemic species are expected to have a larger gene pool and higher phenotypic plasticity.
Location: La Palma, Canary Islands.
Methods: We established 44 sampling sites along a directional gradient of precipitation, heat load, soil nitrogen, phosphorus and pH. Along this gradient, we estimated species abundances and measured three traits (plant height, leaf area and leaf thickness) on perennial endemic and non-endemic plant species. In total, we recorded traits for 1,223 plant individuals of 43 species. Subsequently, we calculated community-weighted mean traits to measure the relative contribution of species turnover, intraspecific variation and their covariation along the analysed gradient.
Results: The contribution of intraspecific variation to total variation was similar in endemic and non-endemic assemblages. For plant height, intraspecific variation explained roughly as much variation as species turnover. For leaf area and leaf thickness, intraspecific variation explained almost no variation. Species turnover effects mainly drove trait responses along the environmental gradient, but intraspecific variation was important for responses in leaf area to precipitation.
Conclusions: Despite their distinct evolutionary history, endemic and non-endemic plant assemblages show similar patterns in species turnover and intraspecific variation. Our results indicate that species turnover is the main component of trait variation in the underlying study system. However, intraspecific variation can increase individual species’ fitness in response to precipitation. Overall, our study challenges the theory that intraspecific trait variation is more important for the establishment of non-endemic species compared with endemic species.
Metasomatic reaction zones between mafic and ultramafic rocks exhumed from subduction zones provide a window into mass-transfer processes at high pressure. However, accurate interpretation of the rock record requires distinguishing high-pressure metasomatic processes from inherited oceanic signatures prior to subduction. We integrated constraints from bulk-rock geochemical compositions and petrophysical properties, mineral chemistry, and thermodynamic modeling to understand the formation of reaction zones between juxtaposed metagabbro and serpentinite as exemplified by the Voltri Massif (Ligurian Alps, Italy). Distinct zones of variably metasomatized metagabbro are dominated by chlorite, amphibole, clinopyroxene, epidote, rutile, ilmenite, and titanite between serpentinite and eclogitic metagabbro. Whereas the precursor serpentinite and oxide gabbro formed and were likely already in contact in an oceanic setting, the reaction zones formed by diffusional Mg-metasomatism between the two rocks from prograde to peak, to retrograde conditions in a subduction zone. Metasomatism of mafic rocks by Mg-rich fluids that previously equilibrated with serpentinite could be widespread along the subduction interface, within the subducted slab, and the mantle wedge. Furthermore, the models predict that talc formation by Si-metasomatism of serpentinite in subduction zones is limited by pressure-dependent increase in the silica activity buffered by the serpentine-talc equilibrium. Elevated activities of aqueous Ca and Al species would also favor the formation of chlorite and garnet. Accordingly, unusual conditions or processes would be required to stabilize abundant talc at high P-T conditions. Alternatively, a different set of mineral assemblages, such as serpentine- or chlorite-rich rocks, may be controlling the coupling-decoupling transition of the plate interface.
Talc formation via silica-metasomatism of ultramafic rocks is believed to play key roles in subduction zone processes. Yet, the conditions of talc formation remain poorly constrained. We used thermodynamic reaction-path models to assess the formation of talc at the slab-mantle interface and show that it is restricted to a limited set of pressure–temperature conditions, protolith, and fluid compositions. In contrast, our models predict that chlorite formation is ubiquitous at conditions relevant to the slab-mantle interface of subduction zones. The scarcity of talc and abundance of chlorite is evident in the rock record of exhumed subduction zone terranes. Talc formation during Si-metasomatism may thus play a more limited role in volatile cycling, strain localization, and in controlling the decoupling-coupling transition of the plate interface. Conversely, the observed and predicted ubiquity of chlorite corroborates its prominent role in slab-mantle interface processes that previous studies attributed to talc.
Key Points:
Limited talc formation by Si-metasomatism of ultramafic rocks in subduction zones
Chlorite formation is likely pervasive at the slab-mantle interface
Preferential formation of chlorite has wide-ranging chemical and physical implications for subduction zone processes
Plain Language Summary: In subduction zones, talc can form during chemical reactions of mantle rocks with silica-enriched fluids at the interface between descending oceanic plates and the overriding mantle. Its formation and distribution in subduction zones are believed to affect the volatile budget, rheological properties, and the down-dip limit of the decoupling of the slab-mantle interface. Therefore, illuminating the conditions that facilitate talc formation at high pressure-temperature conditions is key in assessing its roles in fundamental subduction zone processes. Using thermodynamic reaction-path models, we show that the formation of talc at the slab-mantle interface is restricted to a limited set of environmental conditions, because its formation is highly sensitive to the compositions of the mantle rocks and reactant fluids. Contrary to common belief, talc is unlikely to form in high abundance in ultramafic rocks metasomatized by Si-rich slab-derived fluids. Rather, our models predict the ubiquitous formation of chlorite along with other silicate minerals during Si-metasomatism due to the competing effects from other dissolved components that favor their formation over talc. This study calls into question the importance of talc during Si-metasomatism in subduction zones but highlights the more predominant role of chlorite.
Polarization of Λ and ¯Λ hyperons along the beam direction in Pb-Pb collisions at √sNN=5.02 TeV
(2022)
The polarization of the Λ and ¯Λ hyperons along the beam (z) direction, Pz, has been measured in Pb-Pb collisions at √sNN=5.02 TeV recorded with ALICE at the Large Hadron Collider (LHC). The main contribution to Pz comes from elliptic flow-induced vorticity and can be characterized by the second Fourier sine coefficient Pz,s2=⟨Pzsin(2φ−2Ψ2)⟩, where φ is thhyperon azimuthal emission angle and Ψ2 is the elliptic flow plane angle. We report the measurement of Pz,s2 for different collision centralities and in the 30%–50% centrality interval as a function of the hyperon transverse momentum and rapidity. The Pz,s2 is positive similarly as measured by the STAR Collaboration in Au-Au collisions at √sNN=200 GeV, with somewhat smaller amplitude in the semicentral collisions. This is the first experimental evidence of a nonzero hyperon Pz in Pb-Pb collisions at the LHC. The comparison of the measured Pz,s2 with the hydrodynamic model calculations shows sensitivity to the competing contributions from thermal and the recently found shear-induced vorticity, as well as to whether the polarization is acquired at the quark-gluon plasma or the hadronic phase.
Reliable identification of chondrules, calcium-aluminum-rich inclusions (CAIs), carbonate grains, and Ca-phosphate grains at depth within untouched, unprepared chondritic samples by a nondestructive analytical method, such as synchrotron X-ray fluorescence (SXRF) computed tomography (CT), is an essential first step before intrusive analytical and sample preparation methods are performed. The detection of a local Ca-enrichment could indicate the presence of such a component, all of which contain Ca as major element and/or Ca-bearing minerals, allowing it to be precisely located at depth within a sample. However, the depth limitation from which Ca-K fluorescence can travel through a chondrite sample (e.g., ∼115 µm through material of 1.5 g cm−3) to XRF detectors leaves many Ca-bearing components undetected at deeper depths. In comparison, Sr-K lines travel much greater distances (∼1700 µm) through the same sample density and are, thus, detected from much greater depths. Here, we demonstrate a clear, positive, and preferential correlation between Ca and Sr and conclude that Sr-detection can be used as proxy for the presence of Ca (and, thus, Ca-bearing components) throughout mm-sized samples of carbonaceous chondritic material. This has valuable implications, especially for sample return missions from carbonaceous C-type asteroids, such as Ryugu or Bennu. Reliable localization, identification, and targeted analysis by SXRF of Ca-bearing chondrules, CAIs, and carbonates at depth within untouched, unprepared samples in the initial stages of a multianalysis investigation insures the valuable information they hold of pre- and post-accretion processes in the early solar system is neither corrupted nor destroyed in subsequent processing and analyses.
As part of two drilling campaigns of the International Continental Scientific Drilling Program (ICDP), several geophysical borehole measurements were carried out by the Leibniz Institute for Applied Geophysics (LIAG) in two lakes. The acquired data was used to answer stratigraphic and paleoclimatic research questions, including the establishment of robust age-depth models and the construction of continuous lithological profiles.
Lake Towuti is located on Sulawesi (Indonesia), within the "Indo-Pacific Warm Pool" (IPWP), a globally important region for atmospheric heat and moisture budgets. The lake exists for approximately one million years, but its exact age is uncertain. We present the first agedepth model for the approximately 100 m continuous sediment sequence from the central part of the lake. The basis for this model is the magnetic susceptibility measured in the borehole and a tephra layer with an age of about 797 ka at 72 m depth. Our age-depth model is inferred from cyclostratigraphic analysis of borehole data and covers a period from 903 ± 11 to 131 ± 67 ka. We suggest that orbital eccentricity and/or changes between global cold and warm periods are responsible for hydroclimatic changes in the IPWP, that these changes affect sedimentation processes in Lake Towuti, and that we can measure and observe this effect in the sediment properties today. Additionally, we created a continuous artificial lithological profile from a series of different borehole data using cluster analysis. This provides information from parts of the borehole where no sediment is available due to core loss.
Lake Ohrid is 1.36 million years old and is located on the Balkan Peninsula on the border between Albania and North Macedonia. The primary hole 'DEEP' in the central part of the lake has been the subject of several investigations, but information about sediments of the marginal locations 'Pestani' and 'Cerava' have not been published yet. In our study, we use natural gamma radiation (GR) measured in the borehole to generate an age-depth model for DEEP. This is performed using the correlation of GR to the global LR04 reference record of Lisiecki and Raymo (2005).
The age information is then transferred via prominent seismic marker horizons to the other two sites, Pestani and Cerava, where it provides the first age-control points for the construction of age-depth models from correlation of GR to LR04. The generated age-depth models are tested using cyclostratigraphic methods, but the limits of this approach are revealed. At DEEP, sedimentation rates (SR) from the cyclostratigraphic method and the correlative approach differ by 2.8 %, at Pestani this difference is 16.7 %, and at Cerava the quality of the data does not allow a reliable evaluation of SR using the cyclostratigraphic approach. We used cluster analysis to construct artificial lithological profiles at all three sites and integrated them into the respective age-depth models. This enables us to determine which sediment types were deposited at what time, and we recognize the change between warm and cold periods in the sediment properties at all three locations. The analyses in this study were all performed on borehole and seismic data and thus do not involve sediment core data. Especially at Pestani and Cerava, new insights into the sedimentological history of Lake Ohrid could be obtained.
In the last part we discuss the occurrence of the half-precession (HP) signal in the European region during the last one million years. The focus is on Lake Ohrid, but a range of other proxies, from the eastern Mediterranean, across the European continent, up to Greenland are analyzed in regards to HP. Applying filters, we focus on the frequency range with a period of 13-8.5 ka and only HP remains in the records. We use correlative methods to determine the clarity of the HP signal in proxies distributed across the European realm. Additionally, we determined the development of HP over time. The HP signal is clearest in the southeast and decreases toward the north. It is further more pronounced in interglacial periods and in the younger part (<621 ka) of most proxies. We suggest that there are mechanisms that transmit the HP signal from its origin near the equator to higher latitudes via different processes. In this context, for instance, the African monsoon, the Nile River and the Mediterranean outflow via the Strait of Gibraltar can be important factors.
New particle formation in the upper free troposphere is a major global source of cloud condensation nuclei (CCN)1,2,3,4. However, the precursor vapours that drive the process are not well understood. With experiments performed under upper tropospheric conditions in the CERN CLOUD chamber, we show that nitric acid, sulfuric acid and ammonia form particles synergistically, at rates that are orders of magnitude faster than those from any two of the three components. The importance of this mechanism depends on the availability of ammonia, which was previously thought to be efficiently scavenged by cloud droplets during convection. However, surprisingly high concentrations of ammonia and ammonium nitrate have recently been observed in the upper troposphere over the Asian monsoon region5,6. Once particles have formed, co-condensation of ammonia and abundant nitric acid alone is sufficient to drive rapid growth to CCN sizes with only trace sulfate. Moreover, our measurements show that these CCN are also highly efficient ice nucleating particles—comparable to desert dust. Our model simulations confirm that ammonia is efficiently convected aloft during the Asian monsoon, driving rapid, multi-acid HNO3–H2SO4–NH3 nucleation in the upper troposphere and producing ice nucleating particles that spread across the mid-latitude Northern Hemisphere.
State of the simulation of mesoscale winds in the Mediterranean and opportunities for improvements
(2022)
The Mediterranean region is a densely populated and economically relevant area with complex orography including mountain ranges, islands, and straits. In combination with pressure gradients, this creates many mesoscale wind systems that cause, e.g., wind gusts and wildfire risk in the Mediterranean. This article reviews the recent state of the science of several mesoscale winds in the Mediterranean and associated processes. Previous work, including case studies on several time ranges and resolutions, as well as studies on these winds under future climate conditions, is discussed. Simulations with grid spacings of 25 to 50 km can reproduce winds driven by large-scale pressure patterns such as Mistral, Tramontane, and Etesians. However, these simulations struggle with the correct representation of winds channeled in straits and mountain gaps and around islands. Grid spacings of 1–3 km are certainly necessary to resolve these small-scale features. The smaller grid spacings are widely used in case studies, but not yet in simulations over large areas and long periods, which also could help to understand the interaction between small-scale phenomena in separate locations. Furthermore, by far not all Mediterranean straits, islands, and mountain gaps were studied in-depth and many interesting Mediterranean small-scale winds still need to be studied.
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.
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.
Oceanic islands only comprise a small amount of the Earth’s land area but harbour a disproportionate amount of global biodiversity. This vast diversity is not only reflected in the taxonomic uniqueness of island biota but also in the remarkable evolution of functional traits. Functional traits, i.e. measurable characteristics that strongly influence the fitness of species, determine how a species responds to its environment and can help to gain more insights into the biogeographical, ecological and evolutionary processes that have shaped island biodiversity. However, research in island biogeography has primarily focused on species richness, and knowledge of functional trait patterns on oceanic islands is scarce. Hence, in this dissertation, I have explored how trait-based approaches can increase our understanding of how biodiversity on oceanic islands assembles and how it is driven by the environment. The Canary Islands (Spain) are a particularly suitable model system to investigate patterns and drivers of biodiversity. The archipelago is characterised by a high variation in environmental heterogeneity and inhabits a unique and well-described native flora. Therefore, I have investigated five principal research questions using the flora (Spermatophytes) of the Canary Islands as a study object. First, I have analysed how climate and biogeography shape the assembly of the Canary Islands flora using a novel trait-based approach. Second, the question of whether rare climates link to functional trait distinctiveness in the native Canary Islands flora was addressed. Third, I have examined how intraspecific trait variation is represented in the native flora of oceanic islands focusing on the succulent scrub of La Palma (Canary Islands). Fourth, this dissertation investigated whether scientific floras can be reliable sources for trait data of plants native to oceanic islands. Finally, I have explored how climate change may impact the native Canary Islands flora by analysing possible climate change-induced shifts in plant species distribution and plant traits.
The results of my dissertation expand the understanding of the importance of biogeography and the environment in determining the functional composition of island floras. I have assessed that traits of endemic plant species did not expand the functional trait space of the Canary Islands but were packed with the ones of non-endemic species. This result hints at a trait convergence in endemic species, possibly driven by non-adaptive speciation processes. Moreover, I have evidenced that humidity is a critical driver of functional diversity in native plant assemblages and particularly leads to a high trait convergence in arid environments via environmental filtering. In contrast, alien species have expanded the Canary Islands flora’s functional trait space. I further have shown that in contrast to native species assemblages, alien species assemblages are characterised by an increasing functional diversity with increasing aridity. This contrasting pattern of functional diversity could pose a potential risk to the native flora of the Canary Islands as a low functional diversity is expected to reduce the resilience of species assemblages to the establishment of more functionally diverse alien plant species. However, in this dissertation, I also have revealed that endemic plant species on the Canary Islands show a high intraspecific variation in arid environments, possibly as an adaptation to environmental stress. Intraspecific variation could help endemic plant species have a competitive advantage over alien species and be more resilient to environmental changes. Furthermore, in this dissertation, I have shown that scientific floras and taxonomic monographs could be used to gain information on quantitative functional traits of plants native to oceanic islands. This finding is particularly relevant for advances in trait-based research, as coverage of trait data for oceanic island floras is extremely poor in global trait databases. Hence, for some of the studies included in this dissertation, trait data were retrieved from scientific floras and taxonomic monographs and used to answer novel scientific research questions. Thus, I have used trait data from the literature to analyse the effect of climate change on the range size of plants native to the Canary Islands. Identifying plant species of particular conservation concern is critical on oceanic islands as many island species have limited distributions and small population sizes, and their niche tracking is impeded by insularity. I have revealed that single-island endemic plants gain less and lose more climatically suitable areas than archipelago endemic and non-endemic native plants due to a climate change-induced decrease in precipitation until 2100...
Diamant hat besondere physikalische und optische Eigenschaften sowie eine starke Resistenz gegenüber Strahlenschädigung. Diese Eigenschaften ermöglichen eine vielfältige Anwendung von Diamant in Wissenschaft und Technik, wie zum Beispiel als Sensormaterial in Strahlungsdetektoren.
Kubisches Zirconiumdioxid (ZrO2) wird aufgrund seiner mechanisch und optisch ähnlichen Eigenschaften unter anderem an Stelle von Diamant eingesetzt. Es ist ebenfalls ein geeignetes Material für viele technische Anwendungen und wird durch seine Strahlenresistenz in Strahlungsumgebungen verwendet. Da beide Materialien in diesem Anwendungsbereich hoher energetischer Strahlung ausgesetzt sind, sind Reaktionen auf die Bestrahlung wie etwa strukturelle Veränderungen oder die Änderungen von Materialeigenschaften von großem Interesse.
In der vorliegenden Arbeit wurde die Morphologie, Struktur und physikalischen Eigenschaften von Diamant und Yttriumoxid-stabilisiertem kubischem ZrO2 nach der Bestrahlung mit 14 MeV Au-Ionen und 1.6 GeV Au-Ionen untersucht. Die durch die Bestrahlung verursachten Veränderungen der Oberflächen und der bestrahlten Volumina wurden mit diversen komplementären analytischen Methoden charakterisiert, bewertet und für die verschiedenen Materialien und Ionenenergien verglichen.
Mittels Röntgenfluoreszenzmessungen wurde die Verteilung und Menge an implantiertem Au semi-quantitativ ermittelt. Die Oberflächen der Proben wurden mit optischer Mikroskopie, Rasterkraftmikroskopie, Rasterelektronenmikroskopie, Röntgenreflektometrie und Elektronenrückstreubeugung untersucht. Strukturelle Veränderungen wurden mit Raman-Spektroskopie analysiert. Der elektrische Widerstand, die Dichte, die Härte sowie das Ätzverhalten der bestrahlten Proben wurden ermittelt und geben Auskunft über die Änderung physikalischer Eigenschaften der Materialien.
Diamant und kubisches ZrO2 reagieren sehr unterschiedlich auf die Bestrahlung mit Au-Ionen gleicher Energien und Fluenzen. Die Diamantproben zeigen nach der Bestrahlung mit 14 MeV Au-Ionen deutliche Veränderungen und Schädigungen der Oberfläche sowie des bestrahlten Volumens. Es wird eine Änderung der Struktur, der Dichte, der Härte, des elektrischen Widerstands sowie des Ätzverhaltens der Proben beobachtet, was auf die Amorphisierung von Diamant zurückgeführt wird. Kubisches ZrO2 ist deutlich strahlungsresistenter gegenüber der Bestrahlung mit 14 MeV Au-Ionen. Es werden keine signifikanten strukturellen Änderungen im getesteten Fluenzbereich beobachtet.
Die mit 1.6 GeV Au-Ionen bestrahlten Diamanten zeigen nur geringe Schädigungen und keine deutliche Änderung der Struktur oder der physikalischen Eigenschaften. Die kubischen ZrO2 Proben sind als Folge der Bestrahlung mit 1.6 GeV Au-Ionen zerbrochen, was auf hohe interne Spannung durch Defektbildung zurückgeführt wird.
Extreme convective precipitation events are among the most severe hazards in central Europe and are expected to intensify under global warming. However, the degree of intensification and the underlying processes are still uncertain. In this thesis, recent advances in continuous, radar-based precipitation monitoring and convection-permitting climate modeling are used to investigate Lagrangian properties of convective rain cells such as precipitation intensity, cell area, and precipitation sum and their relationship to large-scale, environmental conditions.
Firstly, convective precipitation objects are tracked in a gauge-adjusted radar-data set and the properties of these cells are related to large-scale environmental variables to investigate the observed super-Clausius-Clapeyron (CC) scaling of convective extreme precipitation. The Lagrangian precipitation sum of convective cells increases with dew point temperature at rates well above the CC-rate with increasing rates for higher dew point temperatures. These varying, high rates are caused by a covarying increase of CAPE with dew point temperature as well as the effect of high vertical wind shear causing an increase in cell area and thus precipitation sum. At the same time, cells move faster at high vertical wind shear so that Eulerian scaling rates are lower than Lagrangian but still above the CC-rate. The results show that wind shear and static instability need to be taken into account when transferring precipitation scaling under current climate conditions to future conditions. Secondly, the representation of convective cell properties in the convection-permitting climate model COSMO-CLM is evaluated. The model can simulate the observed frequency distributions of cell properties such as lifetime, area, mean and maximum intensity, and precipitation sum. The increase of area and intensity with lifetime is also well captured despite an underestimation of the intensity of the most severe cells. Furthermore, the model can represent the temperature scaling of intensity, area, and precipitation sum but fails to simulate the observed increase of lifetime. Thus, the model is suitable to study climatologies of convective storms in Germany. Thirdly, two COSMO-CLM projections at the end of the century under emission scenario RCP8.5 were investigated. While the number of convective cells and their lifetime remain approximately constant compared to present conditions, intensity and area increase strongly. The relative increase of intensity and area is largest for the highest percentiles meaning that extreme events intensify the most. The characteristic afternoon maximum of convective precipitation is damped, and shifted to later times of day which leads to an increase of nighttime precipitation in the future. Scaling rates of cell properties with dew point temperature are nearly identical in present and future in the simulation driven by the EC-Earth model which means that the upper limit of cell properties like intensity, area, and precipitation sum could be predicted from near-surface dew point temperature. However, this result could not be reproduced by the simulation driven by MIROC5 and needs further investigation.
Auf der Suche nach Erfahrungen in den Tropen setzte der Geografiestudent Jürgen Runge das erste Mal in Togo seinen Fuß auf den afrikanischen Kontinent. Aus einem etwas holprigen Start wurde eine große Zuneigung zu Zentral- und Westafrika. Heute ist Runge Direktor des Zentrums für Interdisziplinäre Afrikaforschung an der Goethe-Universität und forscht gemeinsam mit Partnern der Region vor allem zu Landschaftsentwicklung, Flusssedimenten und Klimawandel.
Cratonic eclogite is the product of oceanic crust subduction into the subcontinental lithospheric mantle, and it also is a fertile diamond source rock. In contrast to matrix minerals in magma-borne xenoliths, inclusions in diamond are shielded from external fluids, retaining more pristine information on the state of the eclogite source at the time of encapsulation. Vanadium is a multi-valent element and a widely used elemental redox proxy. Here, we show that that xenolithic garnet has lower average V abundances than garnet inclusions. This partly reflects crystal-chemical controls, whereby higher average temperatures recorded by inclusions, accompanied by enhanced Na2O and TiO2 partitioning into garnet, facilitate V incorporation at the expense of clinopyroxene. Unexpectedly, although diamond formation is strongly linked to metasomatism and xenoliths remained open systems, V concentrations are similar for bulk eclogites reconstructed from inclusions and from xenoliths. This suggests an oxygen-conserving mechanism for eclogitic diamond formation, and implies that eclogite is an efficient system to buffer fO2 over aeons of lithospheric mantle modification by subduction-derived and other fluids.
The Altenberg–Teplice Volcanic Complex (ATVC) is a large ~ NNW–SSE trending volcano-plutonic system in the southern part of the Eastern Erzgebirge (northern Bohemian Massif, south-eastern Germany and northern Czech Republic). This study presents high precision U–Pb CA-ID-TIMS zircon ages for the pre-caldera volcano-sedimentary Schönfeld–Altenberg Complex and various rocks of the caldera stage: the Teplice rhyolite, the microgranite ring dyke, and the Sayda-Berggießhübel dyke swarm. These data revealed a prolonged time gap of ca. 7–8 Myr between the pre-caldera stage (Schönfeld–Altenberg Complex) and the climactic caldera stage. The volcanic rocks of the Schönfeld–Altenberg Complex represent the earliest volcanic activity in the Erzgebirge and central Europe at ca. 322 Ma. The subsequent Teplice rhyolite was formed during a relatively short time interval of only 1–2 Myr (314–313 Ma). During the same time interval (314–313 Ma), the microgranite ring dyke intruded at the rim of the caldera structure. In addition, one dyke of the Sayda-Berggiesshübel dyke swarm was dated at ca. 314 Ma, while another yielded a younger age (ca. 311 Ma). These data confirm the close genetic and temporal relationship of the Teplice rhyolite, the microgranite ring dyke, and (at least part of) the Sayda-Berggießhübel dyke swarm. Remarkably, the caldera formation in the south of the Eastern Erzgebirge (caldera stage of ATVC: 314–313 Ma) and that in the north (Tharandt Forest caldera: 314–312 Ma) occurred during the same time. These data document a large ~ 60 km NNW–SSE trending magmatic system in the whole Eastern Erzgebirge. For the first time, Hf-O-isotope zircon data was acquired on the ring dyke from the ATVC rocks to better characterize its possible sources. The homogeneous Hf-O-isotope zircon data from the microgranite ring dyke require preceding homogenization of basement rocks. Some small-scale melts that were produced during Variscan amphibolite-facies metamorphism show similar Hf-O-isotope characteristics and can therefore be considered as the most probable source for the microgranite ring dyke melt. In addition, a second source with low oxygen isotope ratios (e.g. basic rocks) probably contributed to the melt and possibly triggered the climactic eruption of the Teplice rhyolite as well as the crystal-rich intrusion of the ring dyke.
The climate system is one of the classical examples of a complex dynamical system consisting of interacting sub-systems through mass, momentum, and energy exchange across various spatial and temporal scales. This thesis aims to detect and quantify sub-component interactions from an information exchange (IE) perspective. For this purpose, IE estimators derived from information theory are explored and applied to the available climate data obtained from observations, reanalysis, global and regional climate models. Specifically, this thesis investigates the usefulness of information theory methods for process-oriented climate model evaluation.
Firstly, methods derived from the concepts of information theory such as transfer entropy and information flow along with their linear and non-linear estimation techniques are initially tested and applied to idealized two-dimensional dynamical systems. The results revealed an expected direction and magnitude of IE providing insights into underlying dynamics. However, as expected the linear estimators are robust for linear systems but fail for non-linear systems. Though the non-linear estimators (kernel and kraskov) showed expected results for all the idealized systems, their free tuning parameters are to be tested for consistent results. Moreover, these methods are sensitive to the available time series length.
A real world example case study involving the dynamics between the Indian and Pacific oceans revealed a physically consistent bi-directional IE. However, unexpected IE was detected in the example of North Atlantic and European air temperatures indicating hidden drivers. Though IE provides insights into system dynamics, the availability of time series length and the system at hand must be carefully taken into account before inferring any possible interpretations of the results.
Quantifying the IE from El-Ni\~{n}o southern oscillation (ENSO) and Indian Ocean Dipole (IOD) to the Indian Summer Monsoon Rainfall (ISMR) with the observational and reanalysis data sets revealed that both ENSO and IOD are synergistic predictors for the inter-annual variability of the ISMR over central India i.e., the monsoon core region. Though the investigated three Global Climate Models (GCM) could not reveal the underlying IE dynamics of ENSO, IOD, and ISMR, a Regional Climate Model (RCM) simulation downscaling one of the GCMs with realistic large scale signals across the lateral boundaries showed good agreement with the observations.
Evaluating a coupled regional climate modeling system driven by two different global data sets with IE estimators revealed significant differences between the process chains linking the north-west Mediterranean sea surface temperatures, evaporation, wind speed, and the Vb-cyclone induced precipitation over Danube, Odra, and Elbe catchments in the historical period (1951-2005). Detailed investigation revealed that the north-west Mediterranean Sea in the coupled regional simulation driven by ERA-20C reanalysis corresponded to the Vb-cyclone precipitation over the three catchments while no such correspondence is noted in the EC-EARTH driven simulation. This discrepancy is attributed to the inheritance of the simulation biases from GCM into the RCM. In the future period (1965-2099), no significant changes in the processes are noted from the simulation.
Overall, this thesis used IE estimators in investigating the underlying dynamics of climate system and climate models. The estimators proved useful in providing insights into climate system dynamics assisting in a process based climate model evaluation.
Carbon is an element that controls planetary habitability, and is fundamental for life on Earth. Its behaviour has important consequences for the global climate system, the origin and evolution of life on Earth. While the biosphere and atmosphere’s carbon cycle only accounts for less than 1% of the global carbon budget, hidden reservoirs of deep carbon in the Earth’s interior comprise the predominant storage of carbon on the planet. At the Earth’s surface, 60-70 % of carbon is hosted by carbonate minerals, which are then transported to the Earth’s interior, mainly in the form of sediments, by subduction of the oceanic lithosphere. Subducting plates are subjected to decarbonation, dehydration, and melting with CO2 release via supra-subduction volcanism. Nevertheless, part of the subducted carbonates’ may survive and be further transported to the deep mantle. Direct evidence of the existence of carbonates in the Earth’s interior, possibly reaching down to the lower mantle, comes from the finding of syngenetic inclusions of carbonates in diamonds and mantle xenoliths. The presence of carbonates in the deep Earth has a critical effect on the physical properties of the mantle. Melting and chemical speciation of the mantle are strongly affected by the form of C and carbonate stability. Therefore, the study of the stability and physical properties of carbonates at high pressures and temperatures is fundamental, because understanding the processes involved in the deep carbon cycle helps to improve our picture of the whole mantle.
The systematic characterization of the elastic properties of carbonates as a function of their structure and chemical composition is of great importance because it may allow to identify their presence and distribution by seismology. Inverting seismic observations to successfully constrain the chemical composition and mineralogy of the Earth’s interior requires knowledge of the physical properties of all possible Earth’s materials at pressures and temperatures applicable to the Earth’s interior. Up to now, a multitude of studies has focused on the construction of phase diagrams and structural transitions by means of X-ray diffraction and vibrational spectroscopy experiments.
Few studies are available on the complete elastic tensor of carbonates, however most of the datasets are not accompanied by an accurate characterization of the samples, which are often solid solutions and the exact chemical composition, density or the details about the experimental methods used are not presented. The aim of this thesis is to study the effect of chemical composition on the elastic properties of carbonates, providing a reliable dataset on the elasticity of the main carbonates. In particular, the elastic properties of crystalline aragonite, CaCO3, and Fe-dolomite, (Ca, Mg, Fe)(CO3)2, with different compositions were studied by Brillouin spectroscopy at ambient conditions. Brillouin spectroscopy was also used to investigate the elastic behaviour of amorphous calcium carbonate samples with different water contents (up to 18 wt%) at high pressures, up to 20 GPa.
Furthermore, the importance of cationic substitution on the structure and high pressure behaviour of carbonates was investigated by studying a synthetic CaCO3-SrCO3 solid solution at ambient conditions and at high pressures, up to 10 GPa, by single crystal X-ray diffraction. Finally, the study of the effect of composition on the elastic properties of families of isostructural solids was also extended to a different class of materials, the metal guanidinium formates. The elasticity of a family of perovskite metal organic frameworks, metal guanidinium formates C(NH2)3MII(HCOO)3, with MII =Mn, Zn, Cu, Co, Cd and Ca was investigated by combining Brillouin spectroscopy, resonant ultrasound spectroscopy, density functional theory and thermal diffuse scattering analysis.
Reise ohne Wiederkehr
(2022)
Seasonal forecasting systems still have difficulties predicting temperature over continental regions, while their performance is better over some maritime regions. On the other hand, the land surface is a substantial source of (sub-)seasonal predictability. A crucial land surface component in focus here is the snow cover, which stores water and modulates the surface radiation balance. This paper’s goal is to attribute snow cover seasonal forecasting biases and lack of skill to either initialization or parameterization errors. For this purpose, we compare the snow representation in five seasonal forecasting systems (from DWD, ECMWF, Météo-France, CMCC, and ECCC) and their performances in predicting snow and 2-m temperature over a Siberian region against ERA5 reanalysis and station data. Although all systems use similar atmospheric and land initialization approaches and data, their snow and temperature biases differ in sign and amplitude. Too-large initial snow biases persist over the forecast period, delaying and prolonging the melting phase. The simplest snow scheme (used in DWD’s system) shows too-early and fast melting in spring. However, systems including multi-layer snow schemes (Météo-France and CMCC) do not necessarily perform better. Both initialization and parameterization are causes of snow biases, but, depending on the system, one can be more dominant.
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.
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.
Atmospheric particles play an important role in the radiative balance of the Earth, as well as they affect human health and air quality. Hence, the chemical characterization constitutes a crucial task to determinate their properties, sources and fate. Particularly, the analysis of nanoparticles (d<100 nm) represents an analytical challenge, since these particles are abundant in number but have very little mass.
This accumulative thesis focuses on the chemical characterization of nanoparticles, performed in both laboratory and field studies. Here, I present four manuscripts, two of which are my main project as a lead author.
The first manuscript (Caudillo et al., 2021) focuses on the gas and the particle phase originated from biogenic precursor gases (α-pinene and isoprene). The experiments were performed in the CLOUD chamber at CERN to simulate pure biogenic new particle formation. Both gas and particle phases are measured with a nitrate CI-APi-TOF mass spectrometer, while the TD-DMA is coupled to it for particle-phase measurements, this setup allows a direct comparison as both measurements use the identical chemical ionization and detector. This study demonstrates the suitability of the TD-DMA for measuring newly formed nanoparticles and it confirms that isoprene suppresses new particle formation but contributes to the growth of newly formed particles.
The second manuscript (Caudillo et al., 2022) presents an intercomparison of four different techniques (including the TD-DMA) for measuring the chemical composition of SOA nanoparticles. The measurements were conducted in the CLOUD chamber. The intercomparison was done by contrasting the observed chemical composition, the calculated volatility, and the thermal desorption behavior (for the thermal desorption techniques). The methods generally agreed on the most important compounds that are found in the nanoparticles. However, they did see different parts of the organic spectrum. Potential explanations for these differences are suggested.
The third manuscript (Ungeheuer al., 2022) presents both laboratory and ambient measurements to investigate the ability of lubricant oil to form new particles. These new particles are an important source of ultrafine particles in the areas nearby large airports. The ambient measurements were performed downwind of Frankfurt International Airport, and it was found that the fraction of lubricant oil is largest in the smallest particles. In the laboratory, the main finding was that evaporated lubricant oil nucleates and forms new particles rapidly. The results suggest that nucleation of lubricant oil and subsequent particle growth can occur in the cooling exhaust plumes of aircraft-turbofans.
The fourth manuscript (Wang et al., 2022) is a new particle formation study in the CLOUD chamber at CERN. This study shows that nitric acid, sulfuric acid, and ammonia interact synergistically and rapidly form particles under upper free tropospheric conditions. These particles can grow by condensation (driven by the availability of ammonia) up to CCN sizes and INP particles. The ability of these particles to act as a CCN and INP was also investigated and it was found to be as efficient as for desert dust. This mechanism constitutes an important finding and it can account for previous observations of high concentrations of ammonia and ammonium nitrate over the Asia monsoon region.
Nontarget screening exhibits a seasonal cycle of PM2.5 organic aerosol composition in Beijing
(2022)
The molecular composition of atmospheric particulate matter (PM) in the urban environment is complex, and it remains a challenge to identify its sources and formation pathways. Here, we report the seasonal variation of the molecular composition of organic aerosols (OA), based on 172 PM2.5 filter samples collected in Beijing, China, from February 2018 to March 2019. We applied a hierarchical cluster analysis (HCA) on a large nontarget-screening data set and found a strong seasonal difference in the OA chemical composition. Molecular fingerprints of the major compound clusters exhibit a unique molecular pattern in the Van Krevelen-space. We found that summer OA in Beijing features a higher degree of oxidation and a higher proportion of organosulfates (OSs) in comparison to OA during wintertime, which exhibits a high contribution from (nitro-)aromatic compounds. OSs appeared with a high intensity in summer-haze conditions, indicating the importance of anthropogenic enhancement of secondary OA in summer Beijing. Furthermore, we quantified the contribution of the four main compound clusters to total OA using surrogate standards. With this approach, we are able to explain a small fraction of the OA (∼11–14%) monitored by the Time-of-Flight Aerosol Chemical Speciation Monitor (ToF-ACSM). However, we observe a strong correlation between the sum of the quantified clusters and OA measured by the ToF-ACSM, indicating that the identified clusters represent the major variability of OA seasonal cycles. This study highlights the potential of using nontarget screening in combination with HCA for gaining a better understanding of the molecular composition and the origin of OA in the urban environment.
The modern precipitation balance in southeastern (SE) Brazil is regulated by the South American summer Monsoon and threatened by global climate change. On glacial-interglacial timescales, monsoon intensity was strongly controlled by precession-forced changes in insolation. To date, relatively little is known about the spatiotemporal distribution of tropical precipitation in SE Brazil and the resulting variability of fluvial discharge on glacial-interglacial timescales. Here, we present X-ray diffraction-derived mineralogical data for the 150–70 ka period (marine isotope stage (MIS) 6 to MIS 5) from the Doce River basin. This area was sensitive to changes in monsoonal precipitation intensity due to its proximity to the South Atlantic Convergence Zone. The data, obtained from a marine sediment core (M125-55–7) close to the Doce river mouth (20°S), show pronounced changes in the Doce River suspension load’s mineralogical composition on glacial-interglacial and precessional timescales. While the ratio of silicates to carbonates displays precession-paced changes, the mineralogical composition of the carbonate-free fraction discriminates between two assemblages which strongly vary between glacial and interglacial time scales, with precession-forced variability only visible in MIS 5. The first assemblage, dominated by high contents of kaolinite and gibbsite, indicates intensified lowland erosion of mature tropical soils. The second one, characterized by higher contents of the well-ordered illite, quartz and albite, points to intensified erosion of immature soils in the upper Doce Basin. High kaolinite contents in the silicate fraction prevailed in late MIS 6 and indicate pronounced lowland soil erosion along a steepened topographic gradient. The illite-rich mineral assemblage was more abundant in MIS 5, particularly during times of high austral summer insolation, indicating strong monsoonal rainfall and intense physical erosion in the upper catchment. When the summer monsoon weakened in times of lower insolation, the mineral assemblage was dominated by kaolinite again, indicative of lower precipitation and runoff in the upper catchment and dominant lowland erosion.
This thesis is focusing on the impact of Paratethys and Mediterranean water bodies over the Eurasian climate and the interplay between climate, tectonics and biosphere during the late Miocene. This target was the interval between 12.7 and 7.65 Ma for Paratethys, following the Eastern Paratethys restriction and isolation, and 7.2−6.5 Ma (the early Messinian) in Mediterranean, zooming on the effects of gateway restrictions over the eastern Mediterranean and the new born Aegean domain. In both cases restriction is overlapping with large scale climatic changes and tectonic reconfiguration, leading a sort of symbiotic relationship.
Paratethys was a giant epicontinental sea that covered a large part of Eurasia since Paleogene. Due to the Eurasia-Afro-Arabia collision and formation of the Alpine-Himalayan belt (Rögl, 1999; Popov et al., 2006), the Paratethys was divided during the late Miocene in smaller basins that in time were isolated of each other. The protracted isolation and intense continentalisation of paratethyan realm led to changes in humidity distribution, basin connectivity, sediment sources and salinity. These changes had in turn major consequences over water circulation, water availability, vegetation cover and biota. These changes are more intense after 11.6 Ma, when the Eastern Paratethys lost any sustained marine connection, evolving into an enclosed system with endemic fauna (Harzhauser and Piller, 2007).
Mediterranean Sea is a Mezozoic oceanic relic squeezed between Africa, Europe, Anatolia and Arabia, as Africa continued to subduct beneath the European plate. As opposed to Paratethys, it maintained the open connection with the ocean until Messinian, when the two Atlantic gateways (Betic and Rifian corridors) closed for a short time, isolating the basin. The cut off resulted in a dramatic drop down and onset of evaporitic precipitation in marginal basins, the event receiving the name of Messinian Salinity Crisis (5.97−5.55). The restriction affected all marine ecosystems, due to changes in salinity and stratification of water column.
The main objectives of this thesis were:
(1) build valid paleo-temperatures records for both basins based on biomarkers;
(2) reconstruct the hydrology for the late Miocene time interval;
(3) identify vegetation composition and changes;
(4) identify paleo-fires in the late Miocene sediment records;
(5) identify the biotic response to the overall climate and tectonic changes.
All the above objectives were attained with results published in specific journals (Chapters 5−7).
Based on Panagia section (Taman Peninsula, Russia) the longest Paratethys temperature record was completed (~5 Myr), covering the interval between 12.7 and 7.65 Ma. A comprehensive SST and MAT records was obtained, as well as soil pH and carbon (δ13C) and hydrogen (δ2H) stable isotopic compositions on n-alkanes and alkenones. The main findings are concentrated around three prolonged periods with severe droughts affecting the late Miocene circum-Paratethys region peaking at 9.65, 9.4 and 7.9 Ma, associated with a transition towards open land vegetation, intensification of fire activity and enhanced evaporation and aridity.
The time intervals with dryer conditions recorded in Panagia coincide with periods of mammal turnover and dispersal in Eurasia indicating that major environmental changes occurred in the circum-Paratethys region and Paratethys fragmentation had a great impact on the terrestrial ecosystems, when periods of prolonged droughts generated biotic crises and animal displacements across the Eurasian continent. The δ13CC29n-alkane values and charcoal morphologies from Panagia indicate an increased contribution of C4 plants adapted to drier conditions at 9.66 Ma. Similarly high δ13CC29n-alkane values continue until 9.4 Ma, when in Western Europe increased seasonality accelerated the demise of the evergreen subtropical woodlands and expansion of grasslands from Anatolia and Middle East to Europe.
As a result of basin fragmentation and climatic stress, the Eastern Paratethys sub-basins progressively lost their marine properties and turned into brackish-fresh water bodies fed primarily by riverine input. The shallower areas became in time emerged, obstructing connections and isolating the biota, inducing rapid adjusting or extinctions. Thus, the Paratethys harbored a highly endemic fauna (Rögl, 1999), such as dwarf whales, dolphins, seals (among mammals), as well as fish and other taxa (mollusks, ostracods, diatoms, foraminifera, algae, etc.).
Collectively the data structured and analyzed in chapter five support a model in which the Eastern Paratethys evolved as a largely (en)closed system, registering paleoenvironmental signals that are governed by interbasinal connectivity (or lack of it) and regional climate changes in the basin catchment. Acting as an important source of humidity for Western and Central Asia, the size and areal extent of the Paratethys water body is likely to have had a major impact on hydroclimate patterns in the Eurasian interior, with the cumulative fluctuations in both hydrology and surface temperature enhancing the aridity and seasonality, with different partition of moisture over the year. Our combined data suggests a decoupling of Paratethys from the global system as isolation advanced, dominated by regional tectonics and ultimately the Paratethys volume and areal extent reduction.
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.
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.
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.
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.
"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.
Sedimentary charcoal records are widely used to reconstruct regional changes in fire regimes through time in the geological past. Existing global compilations are not geographically comprehensive and do not provide consistent metadata for all sites. Furthermore, the age models provided for these records are not harmonised and many are based on older calibrations of the radiocarbon ages. These issues limit the use of existing compilations for research into past fire regimes. Here, we present an expanded database of charcoal records, accompanied by new age models based on recalibration of radiocarbon ages using IntCal20 and Bayesian age-modelling software. We document the structure and contents of the database, the construction of the age models, and the quality control measures applied. We also record the expansion of geographical coverage relative to previous charcoal compilations and the expansion of metadata that can be used to inform analyses. This first version of the Reading Palaeofire Database contains 1676 records (entities) from 1480 sites worldwide. The database (RPDv1b – Harrison et al., 2021) is available at https://doi.org/10.17864/1947.000345.
Melting inside earth is a common phenomenon and can be observed in many different regions where melt travels through the mantle and crust to eventually reach the surface where it crystallizes to build large volcanic provinces, whole stratigraphic layers of flood basalts, or even the oceanic crust. Often, melt reaching the surface is a good source of information. It can be used to achieve a better understanding about processes taking place in deeper regions inside the mantle and it is therefore essential to fundamentally understand melting and melt percolation processes. In order to achieve a deeper understanding, the aim of this thesis is to investigate processes that are connected to melting by using numerical models.
The physical model used is a so called two-phase flow model which describes the ability of melt to percolate through a viscously deforming, partially molten matrix. A famous feature of two-phase flow are solitary porosity waves, which are waves of locally higher porosity ascending through a partially molten background, keeping its shape constant, driven by decompaction and compaction of the solid matrix in front and behind the wave.
The viscosity law for shear- and volume viscosity was strongly simplified in most previous studies that modeled solitary waves. Often the porosity dependency is underestimated or its influence on the volume viscosity is even neglected, leading to too high viscosities. In this work more realistic laws are used that strongly decrease for small melt fractions. Those laws are incorporated into a 2D Finite Difference mantle convection code with two-phase flow to study the ascent of solitary porosity waves.
The model results show that an initial Gaussian-shaped wave rapidly evolves into a solitary wave with a certain amplitude, traveling upwards with constant velocity. Even though strongly weaker viscosities are used, the effect on dispersion curves and wave shape are only minor as long as the background porosity is rather small. The results are still in agreement to semi-analytical solutions which neglect shear stresses in the melt segregation equation. Higher background porosities and wave amplitudes lead to significant decrease in phase velocity and wave width, as the viscosity is strongly effected. However, the models show that solitary waves are still a possible mechanism for more realistic matrix viscosities.
While the ascending of porosity waves are mostly described by the movement of fluid melt, partially molten regions inside Earth trigger upwelling of both, solid and fluid phases, which can be called diapirism. While diapirs can have a wide variety of wavelengths, porosity waves are restricted to a few times the compaction length. The size of a melt perturbation in terms of compaction length therefore describes whether material is transported by diapirism or porosity waves. In this thesis we study the transition from diapiric rise to solitary porosity waves by systematically changing the size of a porosity perturbation from 1.8 to 120 times the compaction length. In case of a perturbation of the size of a few times the compaction length a single porosity wave will emerge, either with a positive or negative vertical matrix flux and if melt is not allowed to move relative to the matrix a diapir will emerge. In between these physical end members a regime can be observed where the partially molten perturbation will split up into numerous solitary waves, whose phase velocity is low compared to the Stokes velocity and the swarm of solitary waves will ascend jointly as a diapir, slowly elongating due to a higher amplitude main solitary wave.
Solitary waves will always emerge from a melt perturbation as long as two-phase flow is enabled, but the time for a solitary wave to emerge increases non-linearly with the perturbation radius in terms of compaction length. In nature, in many cases this time might be too long for solitary waves to emerge.
Another important feature when it comes to two-phase flow is the transport of trace elements in melt. Incompatible elements prefer to go into the melt, which eventually enriches the area where it crystallizes again. In order to model this redistribution, the code FDCON was extended to allow for fully consistent transport of elements in melt, including melting, freezing and re-equilibration with time. A 2D model, a simple representation of a volcanic back arc, is set up to investigate the behavior of trace elements. The influence of retention number and re-equilibration time is examined. Lava-lamp like convection can be observed in the lower part of the model, producing melt, that eventually leads to enrichment in trace elements in the upper high-viscous layer. The total enrichment in this layer approaches an asymptotic value and a 0D model is introduced to recreate this behavior.
Titanite is a potentially powerful U–Pb petrochronometer that may record metamorphism, metasomatism, and deformation. Titanite may also incorporate significant inherited Pb, which may lead to inaccurate and geologically ambiguous U–Pb dates if a proper correction is not or cannot be applied. Here, we present laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS)-derived titanite U–Pb dates and trace element concentrations for two banded calcsilicate gneisses from south-central Maine, USA (SSP18-1A and SSP18-1B). Single spot common Pb-corrected dates range from 400 to 280 Ma with ±12–20 Ma propagated 2SE. Titanite grains in sample SSP18-1B exhibit regular core-to-rim variations in texture, composition, and date. We identify four titanite populations: (1) 397 ± 5 Ma (95% CL) low Y + HREE cores and mottled grains, (2) 370 ± 7 Ma high Y + REE mantles and cores, (3) 342 ± 6 Ma cores with high Y + REE and no Eu anomaly, and (4) 295 ± 6 Ma LREE-depleted rims. We interpret the increase in titanite Y + HREE between ca. 397 and ca. 370 Ma to constrain the timing of diopside fracturing and recrystallization and amphibole breakdown. Apparent Zr-in-titanite temperatures (803 ± 36°C at 0.5 ± 0.2 GPa) and increased XDi suggest a thermal maximum at ca. 370 Ma. Population 3 domains dated to ca. 342 Ma exhibit no Eu anomaly and are observed only in compositional bands dominated by diopside (>80 vol%), suggesting limited equilibrium between titanite and plagioclase. Finally, low LREE and high U/Th in Population 4 titanite dates the formation of hydrous phases, such as allanite, during high XH2O fluid infiltration at ca. 295 Ma. In contrast to the well-defined date–composition–texture relationships observed for titanite from SSP18-1B, titanite grains from sample SSP18-1A exhibit complex zoning patterns and little correlation between texture, composition, and date. We hypothesize that the incorporation of variable amounts of radiogenic Pb from dissolved titanite into recrystallized domains resulted in mixed dates spanning 380–330 Ma. Although titanite may reliably record multiple phases of metamorphism, these data highlight the importance of considering U–Pb data along with chemical and textural data to screen for inherited radiogenic Pb.
Sulfur in the slab: a sulfur-isotopes and thermodynamic-modeling perspective from exhumed terranes
(2022)
Sulfur is a key element in the subduction zone-volcanic arc system; however, the mechanism(s) that recycle sulfur from the slab into the overlying volcanic arc are debated. Here we summarize recent advances in quantifying this component of the deep sulfur cycle. First, primary metamorphic or inherited sulfides in oceanic-type eclogites are only rarely observed as inclusions and are typically absent from the rock matrix. Additionally, sulfides are relatively common in rocks metasomatized at the slab-mantle interface by slab-derived fluids during exhumation. Combined, these two observations suggest that sulfur loss from subducted mafic crust is relatively efficient. Thermodynamic modeling in Perple_X using the Holland and Powell (2011) database combined with the Deep Earth Water model suggests that the efficiency and speciation of sulfur loss varies depending on the degree of seafloor alteration prior to subduction and the geothermal gradient of the slab. In relatively cold subduction zones, such as Honshu, slab-fluids derived from subducted mafic crust are predicted to exhibit elevated concentrations of HSO4-, SO42-, HSO3-, and CaSO4(aq), whereas hot subduction zones, such as Cascadia, are predicted to produce slab fluids enriched in HS- and H2S at lower pressures. The oxidation of sulfur expelled from subducted pyrite is balanced by the reduction of Fe3+ to Fe2+, consistent with the low Fe3+/SFe of exhumed eclogites relative to blueschists and altered oceanic crust. Where oxidized S-bearing fluids are produced, they are anticipated to interact with more reduced rocks at the slab-mantle interface and within the mantle wedge, resulting in sulfide precipitation and significant isotopic fractionation. The δ34S values of slab fluids are estimated to fall between -11 and +8 ‰. Rayleigh fractionation during progressive fluid-rock interaction results in fractionations of tens of per mil as oxidized species are depleted and sulfides are precipitated, resulting in δ34S values of sulfides that easily span the -21.7 to +13.9 ‰ range observed in metasomatic sulfides in exhumed high-pressure rocks. However, in subduction zones where reduced species prevail, the S isotopic signature of slab fluids is expected to reflect their source and will exhibit a narrower range in δ34S values. As a result, the δ34S values measured in arc magmas may not always be a reliable indicator of the contribution of different components of the slab, such as sediments vs. AOC. Additionally, the impact of S recycling on the oxygen fugacity of arc magmas is expected to vary both spatially and temporally throughout Earth history.
Ob Klimawandel oder Luftverschmutzung: Die chemischen und physikalischen Prozesse in der Atmosphäre haben wichtige Auswirkungen auf die menschliche Gesundheit und Ökosysteme. Dabei ist die Atmosphäre mehr als ein Gemisch aus Stickstoff, Sauerstoff, Wasserdampf, Helium und Kohlenstoffdioxid. Es gibt zahlreiche Spurengase, deren Gesamtanteil am Volumen weniger als 1 % ausmacht. In dieser Arbeit werden Stickstoffoxide, Schwefeldioxid, Kohlenstoffmonoxid und Schwefelsäure näher betrachtet, die im Rahmen der flugzeugbasierten Messkampagne Chemistry of the Atmosphere: field experiment in Europe (CAFE-EU)/BLUESKY gemessen wurden.
Die Stickstoffoxide NO und NO2, als NOx zusammengefasst, besitzen hauptsächlich anthropogene Quellen, allen voran fossile Verbrennung und industrielle Prozesse. Zwischen NO und NO2 besteht ein photochemisches Gleichgewicht, sodass in der Atmosphäre vor allem NO2 in relevanten Konzentrationen vorkommt; dies wirkt aufgrund der Bildung von Salpetersäure, HNO3, in wässriger Lösung beim Einatmen ätzend und ist entsprechend gesundheitsschädlich. Troposphärisches Ozon, O3, wesentlicher Bestandteil von Sommersmog, wird hauptsächlich durch die Reaktion von NO mit Peroxiden (HO2 und RO2) gebildet. In der Stratosphäre entstehen NOx hauptsächlich durch die Photodissoziation von Lachgas, N2O, das aufgrund seiner langen Lebenszeit von der Tropo- in die Stratosphäre transportiert werden kann und dort die wichtigste Stickstoffquelle darstellt. In der Stratosphäre tragen NOx zum katalytischen Abbaumechanismus des Ozons bei (Bliefert, 2002; Seinfeld and Pandis, 2016).
Schwefeldioxid, SO2, ist ein toxisches Gas, dessen atmosphärische Quellen hauptsächlich anthropogen sind, nämlich fossile Verbrennung und industrielle Prozesse; Senken sind trockene und feuchte Deposition, wobei letztere zu saurem Regen führen kann. Seit den 1980ern sinken die globalen SO2-Emissionen. SO2 kann in der Atmosphäre zu Sulfat und Schwefelsäure oxidiert werden, was Hauptbestandteil des Wintersmogs ist. Der wichtigste Mechanismus ist die Oxidation mit dem Hydroxylradikal, OH˙, unter Beteiligung von Wasserdampf. In der Stratosphäre ist Carbonylsulfid, OCS, die wichtigste Schwefelquelle, da es analog zum N2O dank seiner langen Lebenszeit von der Tropo- in die Stratosphäre transportiert werden kann (Bliefert, 2002; Seinfeld und Pandis, 2016). Typische Konzentrationen von Schwefelsäure sind 105 cm–3 nachts und 107 cm–3 tagsüber in der Troposphäre sowie 105 cm–3 tagsüber in der Stratosphäre (Clarke et al., 1999; Weber et al., 1999; Fiedler et al., 2005; Arnold, 2008; Kürten et al., 2016; Berresheim et al., 2000).
Kohlenstoffmonoxid, CO, ist ein toxisches Gas, das zu gleichen Teilen durch direkte Emissionen (v.a. Biomasseverbrennung und fossile Verbrennung) und In-situ-Oxidation (v.a. von Methan, Isopren und industriellen Kohlenwasserstoffen) in die Atmosphäre gelangt. Die Hauptsenke ist die Reaktion mit OH˙ in der Troposphäre. Seit 2000 sinkt die globale CO-Konzentration (Bliefert, 2002).
Doch neben Gasen sind auch Aerosolpartikel fester Bestandteil des Gemisches Luft, welche luftgetragene feste oder flüssige Teilchen sind. Primäre Aerosolpartikel werden direkt als solche in die Atmosphäre emittiert, während sekundäre Aerosolpartikel in der Atmosphäre gebildet werden, indem gasförmige Vorläufersubstanzen mit geringer Flüchtigkeit auf primären Partikeln kondensieren oder durch Zusammenclustern und Anwachsen komplett neue Partikel bilden. Aerosolpartikel ermöglichen als Wolkenkondensationskeime erst die Bildung von Wolken und wirken somit – neben ihrem direkten reflektierenden Effekt – durch Änderung der Wolkenbedeckung und -eigenschaften insgesamt kühlend aufs Klima und beeinflussen die lokalen und globalen Wasserkreisläufe. Doch sie haben auch negative Auswirkungen auf die menschliche Gesundheit und sind für eine Verkürzung der durchschnittlichen Lebensdauer in Regionen mit hohen Feinstaubbelastungen verantwortlich (Seinfeld und Pandis, 2016; Bellouin et al., 2020; World Health Organization, 2016).
Neben den bisher betrachteten neutralen, also ungeladenen Gasen und Partikeln sind Ionen in der Gasphase sowie geladene Partikel ebenfalls Bestandteil der Atmosphäre. Sie spielen bei vielen atmosphärischen Prozessen eine wichtige Rolle, wie etwa bei Gewittern, Radiowellenübertragung und ionen-induzierter Nukleation von Aerosolpartikeln. Die Hauptquellen für Ionisation in der Tropo- und Stratosphäre ist die galaktische kosmische Strahlung, die entgegen ihrem Namen hauptsächlich aus Protonen und α-Partikeln (primäre Partikel genannt) besteht und in der Erdatmosphäre durch Kollision mit Luftmolekülen Teilchenschauer von sekundären Partikeln (u.a. Myonen, Pionen und Neutrinos) hervorruft. Die primären und sekundären Partikel können die Luftmoleküle ionisieren unter Entstehung von N+, N2+, O+, O2+ und Elektronen. Sauerstoff reagiert rasch mit letzteren zu O– und O2–. Diese Kationen und Anionen reagieren weiter, bis Ionenclustern der Summenformeln (HNO3)n(H2O)mNO3– und H+(H2O)n(B)m gebildet werden, wobei B Basen wie Methanol, Aceton, Ammoniak oder Pyridin sind. Weitere Ionisationsquellen sind der Zerfall des Radioisotops 222Rn in Bodennähe und ionisierende Solarstrahlung oberhalb der Stratosphäre. Atmosphärische Ionen haben zwei wichtige Senken: die Wiedervereinigung, auch Rekombination genannt, bei der sich ein Kation und ein Anion gegenseitig neutralisieren sowie das Anhaften an Aerosolpartikeln. Letztere Senke ist vor allem in der Troposphäre aufgrund der relativ hohen Konzentration an Aerosolpartikeln relevant (Arnold, 2008; Viggiano und Arnold, 1995; Bazilevskaya et al., 2008; Hirsikko et al., 2011).
This thesis presents the experimental and numerical analysis of seismic waves that are produced by wind farms. With the aim to develop renewable energies rapidly, the number of wind turbines has been increased in recent years. Ground motions induced by their operation can be observed by seismometers several kilometers away. Hence, the seismic noise level can be significantly increased at the seismic station. Therefore, this study combines long-term experiments and numerical simulations to improve the understanding of the seismic wavefields emitted by complete wind farms and to advance the prediction of signal amplitudes.
Firstly, wind-turbine induced signals that are measured at a small wind farm close to Würzburg (Germany) are correlated with the operational data of the turbines. The frequency-dependent decay of signal amplitudes with distance from the wind farm is modeled using an analytical method including the complex effects of interferences of the wavefields produced by the multiple wind turbines. Specific interference patterns significantly affect the wave propagation and therefore the signal amplitude in the far field of a wind farm. Since measurements inside the wind turbines show that the assumption of in-phase vibrating wind turbines is inappropriate, an approach to calculate representative seismic radiation patterns from multiple wind turbines, which allows the prediction of amplitudes in the far field of a complete wind farm, is proposed.
In a second study, signals with a frequency of 1.15 Hz, produced by the Weilrod wind farm (north of Frankfurt, Germany) are observed at the seismological observatory TNS (Taunus), which is located at a distance of 11 km from the wind farm. The propagation of the wavefield emitted by the wind farm is numerically modeled in 3D, using the spectral element method. It is shown that topographic effects can cause local signal amplitude reductions, but also signal amplification along the travel path of the seismic wave. The comparison of simulations with and without topography reveals that the reduction and amplification are spatially linked to the shape of the topography, which could be an explanation for the relatively high signal amplitude observed at TNS.
Finally, the reduction of the impact of wind turbines on seismic measurements using borehole installations is studied using 2D numerical models. Possible effects of the seismic velocity, attenuation, and layering of the subsurface are demonstrated. Results show that a borehole can be very effective in reducing the observed high-frequency signals emitted by wind turbines. However, a borehole might not be beneficial if signals with frequencies of about 1 Hz (or lower) are of interest, due significant wavelength-dependent effects. The estimations of depth-dependent amplitudes with a layered subsurface are validated with existing data from wind-turbine-induced signals measured at the top and bottom of two boreholes.
The experimental analysis of measurements conducted at wind farms and the advances of modeling such signals improve the understanding of the propagation of wind-farm induced seismic wave fields. Furthermore, the methods developed in this work have a high potential of universal application to the prediction of signal amplitudes at seismometers close to wind farms with arbitrary layout and geographic location.
Garnet xenocrysts from kimberlites provide unique insights into the composition, structure and evolution of the subcontinental lithospheric mantle (SCLM). For example, different metasomatic events in the SCLM are reflected in compositional differences between garnet xenocrysts. As mantle metasomatism largely controls the physical and chemical properties of the SCLM, it exerts first order control over the genesis of kimberlitic magmas and diamond formation. However, dating mantle lithologies and processes is complicated by high ambient temperatures that allow the equilibration of most isotopic systems up to the time of kimberlite eruption. As a consequence, the temporal connection between metasomatic events in the mantle and kimberlite genesis is commonly ambiguous.
In this study, we applied LA-ICPMS U-Pb dating to 43 harzburgitic, lherzolithic and megacrystic garnet xenocrysts from the ~376 Ma diamondiferous V. Grib kimberlite, Russia, in order to investigate the link between different types of mantle metasomatism and kimberlite genesis.
Our results indicate that, with two possible exceptions, only harzburgitic garnet overlaps in age with the kimberlite eruption, whereas lherzolitic and megacrystic garnet crystals are ~20 to 130 million years older. Furthermore, garnet U-Pb ages and Ni-in-garnet temperatures of ~820 to 1200 °C do not correlate. This, and the high closure temperature of U-Pb in garnet (≥900 °C) suggests that the garnet U-Pb ages indeed reflect metasomatic events in the SCLM. However, the U-Pb ages could also reflect cooling ages. In this case, the metasomatic events recorded in the garnet crystals must still have occurred up to ~130 million years prior to the eruption of the V. Grib kimberlite.
These findings have far-reaching implications for the genesis of (diamondiferous) kimberlites, as they clearly show that the time lag between metasomatic events in the SCLM, as recorded in kimberlitic garnet xenocrysts, and kimberlite eruption may extend to tens of millions of years.
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.
Beim Clusterprojekt ELEMENTS von Goethe-Universität, TU Darmstadt, Universität Gießen und GSI Helmholtzzentrum für Schwerionenforschung arbeiten Theorie und Experiment Hand in Hand, um die Struktur der Materie unter extremen Bedingungen zu verstehen. So wird ersichtlich, warum etwa Kollisionen von Neutronensternen viele der schweren Elemente auf unserem Planeten geliefert haben.
Vor mehr als vier Milliarden Jahren prägten heiße Magma-Ozeane die Oberfläche der Erde. Als die Erde allmählich abkühlte, bildeten sich an manchen Stellen Krusten, später die ersten Kontinente. Die Geowissenschaftlerin Dr. Sonja Aulbach erforscht die damals ablaufenden Prozesse anhand spezieller Gesteinsproben und mit hochmoderner Analysetechnik.
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.
Vladimir Vernadsky's concept of living matter is central to his biogeochemistry, the science he founded. For several reasons, his original understanding of living matter is one of the most complex notions in the history of the life sciences. First, biogeochemistry is by definition an interdisciplinary enterprise that embraces biology, including evolutionary theory, geology, and chemistry, and combines them into a unique research program. Second, if understood in the original sense as used by Vernadsky, living matter is a concept built into idiosyncratic metaphysics constructed around the so-called principle of life's eternity. Third, the concept of living matter reflects the specificity of Vernadsky's sophisticated philosophy of science as he insisted that 'scientific thought' is a planetary phenomenon as well as a geological force. In our contribution, we will introduce Vernadsky's concept of living matter in its historical context. Accordingly, we will also give some chronology of Vernadsky's work related to the growth of his biosphere concept highlighting the 'Ukrainian' period as it is in this period that he intensively elaborated on the notion of living matter. This will be followed by his theory of living matter as it was formulated in his major works of the later period. We are going to locate the notion of living matter within Vernadsky's theoretical system and demonstrate that he regarded his theory of the living as an evolutionary theory complementary to that of Charles Darwin from the very beginning. Additionally, we will briefly present Vladimir Beklemishev's concept of 'geomerida' which he developed at approximately the same time as Vernadsky was elaborating on his 'living matter' to highlight the specificity of the latter's methodology.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
New U–Pb ages of detrital and igneous zircons of the Uppermost Unit of Crete shed light on its provenance and on Eohellenic to Eoalpine imprints in the eastern Mediterranean. The detrital zircons of all nappes show Variscan ages and are characterized by a Minoan-type age spectrum, which is typical for the NE margin of Gondwana. Parts of the metasedimentary rocks are unexpectedly young. Their detrital zircon ages continue via the Permian until the Late Triassic, Middle Jurassic and Early Cretaceous. The high-grade metamorphic rocks of the Asterousia crystalline complex are likely equivalents of the low-grade metamorphic trench and fore-arc deposits of the Vatos nappe pointing to Late Cretaceous slab roll back. The presence of both late Permian detrital zircons and Late Cretaceous arc-type granitoids suggest that the Uppermost Unit of Crete is derived from the late Permian/Late Cretaceous magmatic belt situated north of the Sava–Vardar–Izmir–Ankara Suture in the Strandja–Rhodope area. To achieve their recent position on Crete, the nappes had to travel more than 500 km. The traveling path is well tracked by rocks of the Upper Cycladic Unit, which are similar to those of the Uppermost Unit of Crete. The large displacement of the Cretan nappes was controlled not only by nappe transport, but probably also by dextral strike–slip along the North Anatolian Fault Zone and related counterclockwise rotation of the Anatolian block since the Eocene.
Chloritoid and kyanite coexist in metapelites from the high-pressure/low-temperature Massa Unit in the Alpi Apuane metamorphic complex (Northern Apennines, Italy). The composition of chloritoid is extremely variable throughout the Massa Unit. Fe-chloritoid occurs in association with hematite-free, graphite-bearing schists, whereas strongly zoned Fe-Mg chloritoid is found with hematite and kyanite. We investigated the effect of different bulk Fe2O3 contents in controlling chloritoid composition through phase equilibria modelling of four selected samples, representative of the different chloritoid-bearing parageneses found in the Massa Unit. The ferric iron content, measured through wet chemical titration, ranges from 0 (graphite-chloritoid schist) to 73% of the total iron (hematite-chloritoid schist). We show that Mg-rich chloritoid compositions and stability of kyanite at greenschist to blueschist facies conditions can be reproduced in the MnO–Na2O–K2O–FeO–MgO–Al2O3–SiO2–H2O–TiO2–O (MnNKFMASHTO) chemical system only considering the presence of significant amounts of ferric iron as part of the bulk composition. The stabilization of kyanite at lower grade is directly linked to the presence of Fe2O3, which renders the reactive bulk rock composition effectively enriched in Al2O3 with respect to Fe and Mg. We also document that high Fe2O3 contents exacerbate the effect of chloritoid fractionation, producing strongly zoned Fe-Mg-chloritoid grains. Finally, the P–T modelling of the Massa Units performed in this study allows, for the first time, the recognition of a two-stage evolution at peak conditions, with an earlier pressure peak (1.2–1.3 GPa at 350–400°C), and a later thermal peak (0.7–1.1 GPa at 440–480°C), compatible with subduction, underthrusting and exhumation of the Adria continental margin during growth of the Northern Apennine orogenic wedge.
The oxidation state of sulfur in slab fluids is controversial, with both dominantly oxidized and reduced species proposed. Here we use in situ X-ray absorption spectroscopy analysis of sulfur-in-apatite to monitor changes in the oxidation state of sulfur during high-P metasomatism by slab fluids in the subduction channel. Our samples include a 73 cm continuous transect of reaction zones between a metagabbroic eclogite block and serpentinite matrix from a mélange zone on the island of Syros, Greece. The block core consists of garnet, omphacite, phengite, paragonite, epidote-clinozoisite, and rutile. In this region, apatite is only observed as elongate inclusions in omphacite cores. From the core outwards micas are increasingly replaced by epidote-clinozoisite, garnets are smaller and more frequent, pyrite + bornite is observed as inclusions in recrystallized omphacite, and apatite is increasingly abundant in the matrix and inclusions in garnet. A major transition at 48 cm separates an assemblage of Ca-Na amphibole, omphacite, chlorite, pyrite, and apatite from the inner garnet-bearing eclogite assemblages. Omphacite disappears from the assemblage at ~56 cm and amphibole compositions sharply transition to tremolite at 59 cm. Finally, the assemblage tremolite + talc + pyrite is observed after ~70 cm.Apatites in the eclogite assemblages exclusively display S6+ peaks in their absorption spectra. This includes apatite inclusions in omphacite in the least altered lithology, as well as matrix apatite and isolated apatite inclusions in garnet in the outermost metasomatized eclogite zone. In the intermediate pyrite-rich (~1-5 vol %) amphibole + omphacite + chlorite zone, apatite displays a strong S1- absorption peak in most grains, with rare analyses showing mixed S1- and S6+. Finally, apatite in the outermost tremolite-bearing assemblages only displays a S6+ peak. The pyrite-rich zone at 48 cm occurs at the initial interface between the serpentinite matrix and eclogite block, characterized by a dramatic decrease in Na content and Mg#. Our data suggest that reduction of S6+ in infiltrating fluids to S1- in pyrite became focused as Fe diffused across the steep Mg# gradient, resulting in pyrite precipitation. In contrast, S reduction in the Mg-rich tremolite-dominant portions of the transect was limited by a lack of Fe, resulting in low modes of pyrite and fluid buffered S6+ in apatite. Finally, S6+-bearing apatite is also observed in reaction zone lithologies from elsewhere on Syros, suggesting our observations are not isolated.Two important conclusions are drawn from these data and observations: (1) In the case of Syros, slab fluids at eclogite-facies conditions carried oxidized S6+, and (2) The interaction of these fluids with eclogites composed of ferrous-Fe silicates resulted in extensive sulfide precipitation.
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.
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.
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.
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.
Non-matrix-matched calibration of laser ablation ICPMS (trace/major) element data is a common quantification strategy. However, LA sampling is associated with downhole elemental fractionation, potentially causing inaccuracies if the magnitude of fractionation between the sample and reference material (RM) differs. Here, we estimate fractionation factors (FFs) for different elements (El) in a range of RMs relative to NIST SRM610/612 (FFEl/Ca-NIST) and evaluate element-specific corrections for downhole fractionation using these measured FFEl/Ca-NIST. Significantly different mean El/Ca values were observed before and after correction, particularly for the alkali elements (all RMs), and B, Fe, and Zn (some RMs), notably improving accuracy, especially for the alkali elements. In cases where this methodology does not result in an accuracy improvement, this may help identify underlying issues in reported/reference values for RMs, given that this phenomenon should be accounted for. Overall, we recommend considering routine assessment of FFs and applying a FF correction to enhance data quality.
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.
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.
Carbonate archives record a brief snapshot of the ambient Earth’s surface conditions at their deposition. However, the geologically reasonable extraction and interpretation of geochemical proxy data from ancient, diagenetically altered rock archives is fraught with problems. Three issues stand out: the dichotomy between petrographic and geochemical alteration; the lack of quantitative age constraints for specific diagenetic phases resulting in a poorly constrained admixture of local, basin-wide and over-regional (far-field) features; and an often insufficient understanding of the temperatures and compositions of diagenetic fluids. Here, the archive of Devonian marine limestones exposed to multiple far-field diagenetic events is used as an example to explore the above-listed issues. Methods applied include petrography, micro XRF, fluid inclusion data, clumped isotopes, δ13C and δ18O isotopes, 87Sr/86Sr ratios and quartz trace element data. Devonian limestones studied here were overprinted by two cross-cutting regional fault zones (T ≈ 230 °C) by multiple events between the Variscan Orogeny and the late Paleogene. The following processes are recorded: (i) protolith deposition and partial dolomitisation during rapid burial in the Middle/Late Devonian (T ≈ 180 °C); (ii) deep burial to ca 6.5 km and tectonic/hydrothermal overprint during the Variscan Orogeny in the Carboniferous (T ≈ 90–230 °C); (iii) rapid uplift to 1–2 km burial depth at the end of the Variscan Orogeny and hypogene karstification (T ≈ 50 to 100 °C) initiated by regional geology in the Permian/Triassic; (iv) tectonic/hydrothermal overprint during the opening of the Proto-Atlantic Ocean between the Early Jurassic and the Early Cretaceous (T ≈ 50 to 130 °C); (v) tectonic/hydrothermal overprint including renewed hypogene karstification and hydrothermal calcite cement precipitation (T ≈ 50 to 180 °C) during Alpine Orogeny between the Late Cretaceous and late Paleogene. Despite this complex series of diagenetic events, the protolith limestones largely preserved their respective Middle/Late Devonian dissolved inorganic carbon (DIC) and 87Sr/86Sr signatures. This study documents that geochemical proxy data, placed into their petrographic, paleotemperature, and local to over-regional context, significantly increases the ability to extract quantitative information from ancient carbonate rock archives. Research shown here has wider relevance for carbonate archive research in general.
Highlights
• Constrictional structures range from dome-and-basin folds to coeval folds and boudins.
• Under bulk constriction, the competent layer rotates slower than a passive plane.
• Extension-parallel and –perpendicular folds grow simultaneously.
• Extension-perpendicular folds affect previous boudins.
Abstract
We conducted scaled analogue modelling to show the influence of varying single layer initial orientation on the geometry of folds and boudins in a bulk constrictional strain field. The initial angle between the plane of shortening and the competent layer (θZ(i)) was incrementally increased from 0° to 90° by multiples of 11.25°. While the amount of layer thickening decreased with increasing θZ(i), the deformation structures produced range from pure dome-and-basin folds to coeval folds and boudins. Based on the attitude of fold axes, there are extension-parallel (FEPR) and extension-perpendicular (FEPP) folds, with axes subparallel and subperpendicular to the principal stretching axis (X), respectively. Coeval growth of FEPR folds and boudins occurred when θZ(i) > ca. 25°. The FEPP folds can be subdivided into a first type which affect the entire layer (if θZ(i) ranges between 11.25 and 78.75°) and a second type, referred to as FBEPP folds, which are affecting pre-existing boudins if θZ(i) > 45°. The interlimb angle of all types of folds increases with increasing θZ(i). Folds and boudins similar to the ones produced in this study can be found in salt domes and in tectonites of subduction zones.
Raumkonstruktionen, die über Handlung und Bedeutungszuschreibungen im Kontext sozialer Medien entstehen, sind ein Fallbeispiel dafür, wie vor dem Hintergrund von Digitalität Fachinhalte re-innoviert werden müssen. Mit dem Ziel eines Beitrags zur Konzeption professionellen Lehrkräftewissens im Kontext von Digitalität werden, ausgehend von der Fragestellung, welche professionellen Fähigkeiten Lehrkräfte benötigen, um Raumkonstruktionen im Unterricht zu thematisieren, in dieser kumulativen Dissertation Raumkonstruktionen aus normativer und empirischer Perspektive als exemplarischer geographischer Fachinhalt adressiert.
Als theoretischer Rahmen dient dabei das TPACK Modell von Mishra & Koehler (2006), das professionelles Lehrkräftewissen in die Bereiche fachliches, pädagogisches und technologisches Wissen sowie deren Überschneidungsbereiche einteilt. Zunächst erfolgt eine Anwendung des Modells als Reflexionsperspektive auf Fachinhalte zur Erzeugung normativer Fähigkeitsbeschreibungen. Diese Fähigkeitsbeschreibungen fließen in die Entwicklung eines TPACK-Selbsteinschätzungsfragebogens für Lehramtsstudierende der Geographie ein. Im Rahmen der damit durchgeführten Studie (n= 364) zeigen sich auf deskriptiver Ebene vergleichsweise niedrige Selbsteinschätzungen der Bereiche inhaltlichen und fachdidaktischen Wissens. Durch die Anwendung einer konfirmatorischen Faktorenanalyse kann das TPACK Modell als zufriedenstellend für die Beschreibung der Daten identifiziert werden. Auffällig ist allerdings die niedrige Korrelation des Konstrukts technologischen Wissens mit den angrenzenden Wissensbereichen. In Bezug auf die Selbsteinschätzungen im phasen- und kontextübergreifenden Vergleich lässt sich, ausgehend von linearen Regressionsanalysen, eine tendenzielle Zunahme der selbsteingeschätzten Fähigkeiten entlang der Fachsemesterzahl ermitteln. Mittels Zweistichproben-t-Tests können außerdem höhere Selbsteinschätzungen der pädagogischen Wissensbereiche durch Studierende, die als Vertretungslehrkräfte tätig sind, festgestellt werden.
In Bezug auf die zur Thematisierung von Raumkonstruktionen im Unterricht benötigten Fähigkeiten ist zunächst die Relevanz der Förderung der Integration fachlichen, pädagogischen und fachdidaktischen Wissens hervorzuheben. Die Studie gibt darüber hinaus Hinweise auf eine niedrigere Bedeutung technologischen Wissens im Hinblick auf soziale Medien als Beispiele alltäglicher Technologien. Vor dem Hintergrund des positiven Effekts von Selbstwirksamkeitsprozessen bieten die Ergebnisse Implikationen für eine Diskussion von praktischen Erfahrungen als Aspekt professionellen Lehrkräftewissens im Sinne der Förderung einer kritisch-reflexiven Auseinandersetzung mit der Tätigkeit als Vertretungslehrkraft aus fachdidaktischer Perspektive. Insgesamt leistet die kumulative Dissertation einen Beitrag zum Diskurs um fachliche geographische Bildung im Kontext von Digitalität und zur Konzeption des professionellen Lehrkräftewissens vor diesem Hintergrund.
Highlights
• Solidification and cooling of an intruded dyke or sill within the middle or shallow crust generate stresses of order 200 MPa, which relax on time scales of a few years to million years.
• Stresses may exceed the brittle strength forming tensile fractures.
• Combined with the pressure gradient within the over-pressurized felsic melts, this explains the migration of felsic contact melt into shrinkage cracks (Sederholm-type veins).
Abstract
Rapid emplacement of a mafic dyke or sill at mid-crustal depth heats and possibly melts the felsic wall rock followed by solidification. Associated volume changes generate stresses, possibly enforcing brittle failure and melt migration. We model the evolution of melting, solidification, temperature, and stress including visco-elastic relaxation in 1D - dykes or -sills using realistic rock rheologies of the Weschnitz pluton (Odenwald). For deep emplacement (Case 1, 15.3 km) extensive contact melting of the wall rock occurs, for shallow emplacement (Case 2, 10 km) it is negligible. The stresses are zero at high melt fractions, but increase during solidification and cooling: The intrusion orthogonal stress is always zero. The intrusion parallel stress σ‖ within the intrusion is tensile (O(200 MPa)). It relaxes on a time scale between a few years (Case 1) and 0.6 m.y. (Case 2). Within the wall rock σ‖ is compressive during heating, but becomes tensile under solidification and cooling. Wall rock stresses relax on a time scale of months to 100 years. A Deborah number is defined based on viscous to thermal relaxation allowing generalization of our results. Adding lithostatic stresses, the total stresses of Case 1 remain below the brittle strength, while for Case 2 they may exceed it. Adding the lithostatic pressure to the melt pressure, the effective stresses exceed the brittle strength and intrusion orthogonal tensile fractures are predicted. Combined with the pressure gradient within the over-pressurized felsic melts generated in the wall rock, this explains the migration of felsic contact melt into shrinkage cracks of the mafic sill in the Weschnitz pluton.