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)
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.