Institutes
Refine
Year of publication
Document Type
- Doctoral Thesis (44) (remove)
Language
- English (34)
- German (9)
- Multiple languages (1)
Has Fulltext
- yes (44)
Is part of the Bibliography
- no (44)
Keywords
- Adaptation (1)
- Armenien (1)
- Biomineralization (1)
- COVID-19 lockdown (1)
- Carbonate (1)
- Carbonated apatite (1)
- Chemistry (1)
- Climate Change (1)
- Climate System, Climate Models, Information theory (1)
- Clumped isotopes (1)
Institute
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.
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.
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.
First-principles modeling techniques offer the ability to simulate a wide range of systems under different physical conditions, such as temperature, pressure, and composition, without relying on empirical knowledge. Density functional theory (DFT), a quantum mechanical method, has become an exceptionally successful framework for materials science modeling. Employing DFT makes it possible to gain valuable insights into the fundamental state of a system, enabling the reliable determination of equilibrium crystal structures. Over time, DFT has become an essential tool that can be incorporated into various schemes for predicting the properties of a material related to its structure, insulating/metallic behavior, magnetism, and optics. DFT is regularly applied in numerous fields, spanning from fundamental subjects in condensed matter physics to the study of large-scale phenomena in geosciences. In the latter, the effectiveness of DFT stems from its ability to simulate the properties found on the Earth, other planets, and meteorites, which may pose challenges for their direct study or laboratory investigation.
In this thesis, a comprehensive examination of a family of monosulfides and a perovskite heterostructure was conducted. These materials are relevant for their potential applications in technology, energy harvesting, and in the case of monosulfides, their speculated abundance on the planet Mercury.
Firstly, a DFT approach was used to analyze two non-magnetic monosulfides, CaS and MgS. We determined their structural properties and then focused on the modeling of their reflectivity in the infrared region. The calculation of the reflectivity considered both harmonic and anharmonic contributions. In the harmonic limit, the non-analytic correction was employed to accurately determine the LO/TO splitting, which is necessary to delimit the retstrahlend band, that is, the maximum of the reflectivity. The anharmonic effects given by up to three-phonon and isotopic scatterings, which were included using perturbation theory, primarily smeared the reflectivity spectra edges in the high-wave region.
Secondly, four polymorphs of MnS were studied using a combination of first-principles methods to simulate their antiferromagnetic (AFM) and paramagnetic (PM) states. The integration of DFT+$U$ with special quasirandom structures (SQS) supercells, and occupation matrix control techniques was crucial for achieving convergence, structural optimization accuracy, and obtaining finite energy band gaps and local magnetic moments in the PM phases. The addition of the Hubbard $U$ correction was necessary to treat the highly-correlated Mn $d$-electrons. The success of our approach was clear based on our electronic structure predictions for the PM rock-salt B1-MnS polymorph. Experimentally this phase has been observed to be an insulator, but multiple \emph{ab initio} works resulted previously in metallic behavior. Our computations, on the other hand, predicted insulating and magnetic properties that compare well with available measurements. Additionally, the pressure-field stability of the four MnS polymorphs was studied. In the case of the PM phases, B1-MnS was identified to be the most stable up to about 21 GPa, then transforming into the B31-MnS polymorph. This finding was in close agreement with high-pressure experiments reporting a similar phase transformation. The optical properties of B1-, B4-, and B31-MnS were also simulated. The SQS technique was used to obtain soft-mode-free phonon band structures within the harmonic approximation. Then, the anharmonic effects were included, and the reflectivity was calculated for B1-MnS and B4-MnS. In both cases, a good agreement for the LO/TO splitting was achieved in comparison to experimental results.
Lastly, the oxygen-deficient heterostructure of LaAlO$_{3-\delta}$ /SrTiO$_{3-\delta}$ was investigated also employing DFT+$U$, with a particular emphasis on the potential impact of vacancy clustering at the interface. Six distinct configurations of pairs of vacancies were studied and their energies were compared to find the most stable one. The orbital reconstruction of Ti orbitals was also examined based on their location with respect to the vacancies and the local magnetic moments were calculated. The final results showed that linearly arranged vacancies located opposite to Ti ions give the most energetically stable configuration.
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.
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.
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.
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.
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).
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.