Refine
Year of publication
Document Type
- Doctoral Thesis (34) (remove)
Language
- English (34) (remove)
Has Fulltext
- yes (34)
Is part of the Bibliography
- no (34) (remove)
Keywords
- Adaptation (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)
- Dynamics (1)
Institute
- Geowissenschaften / Geographie (34) (remove)
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).
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.
During the last decades mammalian intracranial structures like the ethmoidal region have rarely been a focus of morphological studies, as they required invasive techniques. Contrary, the ontogeny of the fetal nasal capsule could easily be investigated based on histological material. Since the early 21st century modern imaging techniques like high-resolution computed tomography (μCT) reveal non-destructive insights into the mammalian skull. Furthermore, visualization software enables the virtual reconstruction of the tissues and additionally their morphometric analyses. However, the use of morphometric approaches on the nasal cavity is still scarce. Moreover, the turbinal skeleton is generally regarded as a unit, or the rostral respiratory part is compared to the caudal olfactory part; but the distinct olfactory turbinals have been considered only in a few studies.
The present study focuses on the highly diverse facial shape of the dog (Canis lupus familiaris) that evolved during domestication. Due to human-controlled breeding and care the natural selective pressure in prehistoric dogs has been replaced continually by artificial selection. As a consequence, harmful mutations on gene loci which e.g., control facial length growth got fixed within an extremely short time. According to veterinarian studies the turbinals of short snouted breeds continue their growth after the elongation of the facial bones has stopped prematurely. However, such investigations are based on low-resolution CT or MRT data and the morphological descriptions are vague. Referring to the elongation of the face in dolichocephalic breeds no former study has dealt with the detailed morphology of their turbinal skeleton so far.
The current study is based on comparative anatomical, morphometric, morphofunctional, and ontogenetic patterns of the dog’s turbinal skeleton. The 32 macerated skulls and four histological serial sections represent eleven breeds which cover different snout lengths (brachycephalic, mesaticephalic, dolichocephalic; according to two length indices), functional groups (scent hound, sighthound, companion/toy), and breeding histories (ancient pure-breeding associated with an unchanged appearance, modern time fashion breeding). The nasal cavity of the selected skulls was μCT-scanned and virtual 3D models of the turbinal skeleton were reconstructed. The breeds have been compared with each other in their number of olfactory turbinals, in the morphology of all turbinals and the lamina semicircularis as well as in their morphometrics and ontogeny. Based on morphological and ontogenetic patterns a new terminology of the interturbinals was established. The morphometric data covers the measurement of the relative turbinal surface area (IAT) and the calculation of the surface density (SDEN) and the turbinal complexity (TC). For the latter parameter a new morphometric approach was developed. For the ontogenetic comparison histological serial sections of perinatal dog stages have been consulted. As the dog’s ancestor macerated skulls of three adult Eurasian wolves (Canis lupus lupus) function for outgroup comparison and represent the grundplan with which the breeds are compared.
The results support former studies concerning a species-specific number of the fronto- and ethmoturbinals: in the Eurasian wolf and all postnatal dogs under study three ethmoturbinals and three frontoturbinals are observed. Additionally, two types of interturbinals are distinguished, namely four prominent interturbinals which are present in nearly all individuals and show a homologous pattern, and a variable number of additional interturbinals which differ in their shape among the dogs. Generally, longer snouted breeds have more additional interturbinals, so the total number of olfactory turbinals is increased to a maximum of 16 in the borzoi, whereas several short snouted breeds have only nine olfactory turbinals due to the loss of additional interturbinals and one prominent interturbinal. Regarding ontogeny the growth of the respiratory and the olfactory turbinals and the lamina semicircularis is highly associated with the growth of the facial bones after birth. As the viscerocranium of brachycephalic breeds is subjected to a postnatal growth inhibition the ethmoidal region stops growing prematurely, too. The turbinals of both functional parts develop less accessory lamellae that results in the reduction of the three morphometric parameters IAT, SDEN, and TC. The increase of all these three parameters with increasing snout length proves a correlation between both variables in the maxilloturbinal, all olfactory turbinals, and the lamina semicircularis in the dog. With the help of the perinatal dog stages plesiomorphic patterns which are present in all adult specimens (e.g., separation of ethmoturbinal I into two laminae, the presence of the uncinate process) were distinguished from less established morphological traits which get preferably reduced in association with brachycephaly (e.g., the anterior process of the posterior lamina of ethmoturbinal I, the caudal processes of frontoturbinal 1 and 2 within the frontal sinus due to the latter’s reduction). Obviously, the driving mechanism behind these and further variations are mutations on gene loci which control ontogenetic processes: the in other studies already described postnatal growth inhibition in the dermal bones of the midface of brachycephalic breeds seems to have a similar effect on the ethmoidal region. The results of the present study serve as basis for the evaluation how far the bony turbinals’ morphology, morphometrics, and ontogeny might be associated with physiological, genetic, neurological, and phylogenetic patterns. Additionally, the growth patterns of the hard tissues need to be compared to those of the soft tissues (i.e. the nasal epithelium).
As part of two drilling campaigns of the International Continental Scientific Drilling Program (ICDP), several geophysical borehole measurements were carried out by the Leibniz Institute for Applied Geophysics (LIAG) in two lakes. The acquired data was used to answer stratigraphic and paleoclimatic research questions, including the establishment of robust age-depth models and the construction of continuous lithological profiles.
Lake Towuti is located on Sulawesi (Indonesia), within the "Indo-Pacific Warm Pool" (IPWP), a globally important region for atmospheric heat and moisture budgets. The lake exists for approximately one million years, but its exact age is uncertain. We present the first agedepth model for the approximately 100 m continuous sediment sequence from the central part of the lake. The basis for this model is the magnetic susceptibility measured in the borehole and a tephra layer with an age of about 797 ka at 72 m depth. Our age-depth model is inferred from cyclostratigraphic analysis of borehole data and covers a period from 903 ± 11 to 131 ± 67 ka. We suggest that orbital eccentricity and/or changes between global cold and warm periods are responsible for hydroclimatic changes in the IPWP, that these changes affect sedimentation processes in Lake Towuti, and that we can measure and observe this effect in the sediment properties today. Additionally, we created a continuous artificial lithological profile from a series of different borehole data using cluster analysis. This provides information from parts of the borehole where no sediment is available due to core loss.
Lake Ohrid is 1.36 million years old and is located on the Balkan Peninsula on the border between Albania and North Macedonia. The primary hole 'DEEP' in the central part of the lake has been the subject of several investigations, but information about sediments of the marginal locations 'Pestani' and 'Cerava' have not been published yet. In our study, we use natural gamma radiation (GR) measured in the borehole to generate an age-depth model for DEEP. This is performed using the correlation of GR to the global LR04 reference record of Lisiecki and Raymo (2005).
The age information is then transferred via prominent seismic marker horizons to the other two sites, Pestani and Cerava, where it provides the first age-control points for the construction of age-depth models from correlation of GR to LR04. The generated age-depth models are tested using cyclostratigraphic methods, but the limits of this approach are revealed. At DEEP, sedimentation rates (SR) from the cyclostratigraphic method and the correlative approach differ by 2.8 %, at Pestani this difference is 16.7 %, and at Cerava the quality of the data does not allow a reliable evaluation of SR using the cyclostratigraphic approach. We used cluster analysis to construct artificial lithological profiles at all three sites and integrated them into the respective age-depth models. This enables us to determine which sediment types were deposited at what time, and we recognize the change between warm and cold periods in the sediment properties at all three locations. The analyses in this study were all performed on borehole and seismic data and thus do not involve sediment core data. Especially at Pestani and Cerava, new insights into the sedimentological history of Lake Ohrid could be obtained.
In the last part we discuss the occurrence of the half-precession (HP) signal in the European region during the last one million years. The focus is on Lake Ohrid, but a range of other proxies, from the eastern Mediterranean, across the European continent, up to Greenland are analyzed in regards to HP. Applying filters, we focus on the frequency range with a period of 13-8.5 ka and only HP remains in the records. We use correlative methods to determine the clarity of the HP signal in proxies distributed across the European realm. Additionally, we determined the development of HP over time. The HP signal is clearest in the southeast and decreases toward the north. It is further more pronounced in interglacial periods and in the younger part (<621 ka) of most proxies. We suggest that there are mechanisms that transmit the HP signal from its origin near the equator to higher latitudes via different processes. In this context, for instance, the African monsoon, the Nile River and the Mediterranean outflow via the Strait of Gibraltar can be important factors.
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.
The weather of the atmospheric boundary layer significantly affects our life on Earth. Thus, a realistic modelling of the atmospheric boundary layer is crucial. Hereby, the processes of the atmospheric boundary layer depend on an accurate representation of the land-atmosphere coupling in the model. In this context the land surface temperature (LST) plays an important role. In this thesis, it is examined if the assimilation of LST can lead to improved estimates of the boundary layer and its processes.
To properly assimilate the LST retrievals, a suitable model equivalent in the weather prediction model is necessary. In the weather forecast model of the German Weather Service used here, the LST is modelled without a vegetation temperature. To compensate for this deficit, two different vegetation parameterizations were investigated and the better one, a conductivity scheme, was implemented. In order to make optimal use of the influence of the assimilation of the LST observation on the model system, it is useful to pass on the information of the observation to land and atmosphere already in the assimilation step. For that reason, a fully coupled land-atmosphere prediction model was used. Therefore, the existing control vector of the assimilation system, a local ensemble transform Kalman filter, was extended by the soil temperature and moisture. In two-day case studies in March and August 2017, different configurations of the augmented assimilation system were evaluated based on observing system simulation experiments (OSSE).
LST was assimilated hourly over two days in the weakly and strongly coupled assimilation system. In addition, every six hours a free 24-hour forecast was simulated. The experiments were validated with the simulated truth (a high-resolution model run) and compared against an experiment without assimilation. It was shown that the prediction of the boundary layer temperature, especially during the day, and the prediction of the soil temperature, during the whole day and night, could be improved.
The best impact of LST assimilation was achieved with the fully coupled system. The humidity variables of the model benefited only partially from the LST assimilation. For this reason, covariances in the model ensemble were investigated in more detail. To check their compatibility with the high-resolution model run the ensemble consistency score was introduced. It was found that the covariances between the LST and the temperatures of the high-resolution model run were better represented in the ensemble than those between the LST and the humidity variables.
Geochemical investigations on biogenic carbonates are commonly conducted to reconstruct the environmental conditions of the past. However, different carbonate producers incorporate elements to varying degrees, due to biological vital effects. Detecting and quantifying these effects is crucial to produce reliable reconstructions. These paleoreconstructions are of great importance to evaluate the consequences of our recent climate change and identify control mechanisms on the distribution of endangered species such as Desmophyllum pertusum. In chapter three we tested Mg/Ca, Sr/Ca and Na/Ca ratios on this species, among other coldwater scleractinians, to test if they provide reliable proxy information. The results reveal no apparent control of Mg/Ca or Sr/Ca ratios through seawater temperature, salinity or pH. Na/Ca ratios appear to be partly controlled by the seawater temperature, which is also true for other aragonitic organisms such as warm-water corals and the bivalve Mytilus edulis. However, a large variability complicates possible reconstructions by means of Na/Ca. In addition, we explore different models to explain the apparent temperature effect on Na/Ca ratios based on temperature sensitive Na and Ca pumping enzymes.
The bivalve Acesta excavata is commonly found in cold-water coral reefs among the North Atlantic, together with D. pertusum. Multiple linear regression analysis, presented in chapter four, indicates that up to 79% of the elemental variability in Mg/Ca, Sr/Ca and Na/Ca is explainable with temperature and salinity as independent predictor variables. Vital effects, for instance growth rate effects, are evident and make paleoreconstructions not feasible. Furthermore, organic material embedded in the shell, as well as possible stress effects can drastically change the elemental composition. Removal of these organic matrices from bulk samples for LA-ICP-MS (laser ablation inductively coupled mass spectrometer) measurements by means of oxidative cleaning is not possible, but Na/Ca ratios decrease after this cleaning. This is presumably an effect of leaching and not caused by the removal of organic matrices.
Interesting biogeochemical relations were found in the parasitic foraminifera H. sarcophaga. We report Mg/Ca, Sr/Ca, Na/Ca and Mn/Ca ratios measured in H. sarcophaga from two different host species (A. excavata and D. pertusum) in chapter five. Sr/Ca ratios are significantly higher in foraminifera that lived on D. pertusum. This could indicate that dissolved host material is utilized in shell calcification of H. sarcophaga, given the naturally higher strontium concentration in the aragonite of D. pertusum. Mn/Ca ratios are highest in foraminifera that lived on A. excavata but did not fully penetrate the host’s shell. Most likely, this represents a juvenile stadium of the foraminifera during which it feeds on the organic
periostracum of the bivalve, which is enriched in Mn and Fe. The isotopic compositions are similarly affected, both δ18O and δ13C values are significantly lower in foraminifera that lived 23on D. pertusum compared to specimen that lived on A. excavata. Again, this might represent the uptake of dissolved host material or different pH regimes in the calcifying fluid of the hosts (bivalve < 8, coral > 8) that control the extent of hydration/hydroxylation reactions. Temperature reconstructions are possible using stable oxygen isotopes on this foraminifera species; however, the results are only reliable if the foraminifera lived on A. excavata. Samples of H. sarcophaga from D. pertusum would lead to overestimations of the seawater temperature due to the lower δ18O values.
Apart from biological vital effects, storage and preservation methods can significantly change the geochemical composition of different marine biogenic carbonates. In chapter six this is presented on the example of ethanol preservation, a common technique to allow extended storage of biogenic samples. The investigation reveals a significant decrease of Mg/Ca and Na/Ca ratios even after only 45 days storage in ultrapure ethanol. Sr/Ca ratios on the other hand are not influenced.
Besides temperature, salinity and pH further environmental parameters are important such as nutrient availability, especially for the distribution of cold-water corals. In chapter seven we extend the investigations on A. excavata by including the elemental ratios Ba/Ca, Mn/Ca and P/Ca. We expected P/Ca to be helpful in the otherwise difficult process of dentifying growth increments. Based on our observations we had to refute this theory. P/Ca ratios are not systematically enriched in the vicinity of growth lines. Instead, we found a regular sequence of peaks of Ba/Ca, P/Ca and Mn/Ca. This sequence as well as the peaks in general are potentially caused by equential blooms of different algae, diatoms and other planktonic organisms ...
Extreme convective precipitation events are among the most severe hazards in central Europe and are expected to intensify under global warming. However, the degree of intensification and the underlying processes are still uncertain. In this thesis, recent advances in continuous, radar-based precipitation monitoring and convection-permitting climate modeling are used to investigate Lagrangian properties of convective rain cells such as precipitation intensity, cell area, and precipitation sum and their relationship to large-scale, environmental conditions.
Firstly, convective precipitation objects are tracked in a gauge-adjusted radar-data set and the properties of these cells are related to large-scale environmental variables to investigate the observed super-Clausius-Clapeyron (CC) scaling of convective extreme precipitation. The Lagrangian precipitation sum of convective cells increases with dew point temperature at rates well above the CC-rate with increasing rates for higher dew point temperatures. These varying, high rates are caused by a covarying increase of CAPE with dew point temperature as well as the effect of high vertical wind shear causing an increase in cell area and thus precipitation sum. At the same time, cells move faster at high vertical wind shear so that Eulerian scaling rates are lower than Lagrangian but still above the CC-rate. The results show that wind shear and static instability need to be taken into account when transferring precipitation scaling under current climate conditions to future conditions. Secondly, the representation of convective cell properties in the convection-permitting climate model COSMO-CLM is evaluated. The model can simulate the observed frequency distributions of cell properties such as lifetime, area, mean and maximum intensity, and precipitation sum. The increase of area and intensity with lifetime is also well captured despite an underestimation of the intensity of the most severe cells. Furthermore, the model can represent the temperature scaling of intensity, area, and precipitation sum but fails to simulate the observed increase of lifetime. Thus, the model is suitable to study climatologies of convective storms in Germany. Thirdly, two COSMO-CLM projections at the end of the century under emission scenario RCP8.5 were investigated. While the number of convective cells and their lifetime remain approximately constant compared to present conditions, intensity and area increase strongly. The relative increase of intensity and area is largest for the highest percentiles meaning that extreme events intensify the most. The characteristic afternoon maximum of convective precipitation is damped, and shifted to later times of day which leads to an increase of nighttime precipitation in the future. Scaling rates of cell properties with dew point temperature are nearly identical in present and future in the simulation driven by the EC-Earth model which means that the upper limit of cell properties like intensity, area, and precipitation sum could be predicted from near-surface dew point temperature. However, this result could not be reproduced by the simulation driven by MIROC5 and needs further investigation.
The climate system is one of the classical examples of a complex dynamical system consisting of interacting sub-systems through mass, momentum, and energy exchange across various spatial and temporal scales. This thesis aims to detect and quantify sub-component interactions from an information exchange (IE) perspective. For this purpose, IE estimators derived from information theory are explored and applied to the available climate data obtained from observations, reanalysis, global and regional climate models. Specifically, this thesis investigates the usefulness of information theory methods for process-oriented climate model evaluation.
Firstly, methods derived from the concepts of information theory such as transfer entropy and information flow along with their linear and non-linear estimation techniques are initially tested and applied to idealized two-dimensional dynamical systems. The results revealed an expected direction and magnitude of IE providing insights into underlying dynamics. However, as expected the linear estimators are robust for linear systems but fail for non-linear systems. Though the non-linear estimators (kernel and kraskov) showed expected results for all the idealized systems, their free tuning parameters are to be tested for consistent results. Moreover, these methods are sensitive to the available time series length.
A real world example case study involving the dynamics between the Indian and Pacific oceans revealed a physically consistent bi-directional IE. However, unexpected IE was detected in the example of North Atlantic and European air temperatures indicating hidden drivers. Though IE provides insights into system dynamics, the availability of time series length and the system at hand must be carefully taken into account before inferring any possible interpretations of the results.
Quantifying the IE from El-Ni\~{n}o southern oscillation (ENSO) and Indian Ocean Dipole (IOD) to the Indian Summer Monsoon Rainfall (ISMR) with the observational and reanalysis data sets revealed that both ENSO and IOD are synergistic predictors for the inter-annual variability of the ISMR over central India i.e., the monsoon core region. Though the investigated three Global Climate Models (GCM) could not reveal the underlying IE dynamics of ENSO, IOD, and ISMR, a Regional Climate Model (RCM) simulation downscaling one of the GCMs with realistic large scale signals across the lateral boundaries showed good agreement with the observations.
Evaluating a coupled regional climate modeling system driven by two different global data sets with IE estimators revealed significant differences between the process chains linking the north-west Mediterranean sea surface temperatures, evaporation, wind speed, and the Vb-cyclone induced precipitation over Danube, Odra, and Elbe catchments in the historical period (1951-2005). Detailed investigation revealed that the north-west Mediterranean Sea in the coupled regional simulation driven by ERA-20C reanalysis corresponded to the Vb-cyclone precipitation over the three catchments while no such correspondence is noted in the EC-EARTH driven simulation. This discrepancy is attributed to the inheritance of the simulation biases from GCM into the RCM. In the future period (1965-2099), no significant changes in the processes are noted from the simulation.
Overall, this thesis used IE estimators in investigating the underlying dynamics of climate system and climate models. The estimators proved useful in providing insights into climate system dynamics assisting in a process based climate model evaluation.
State-of-the-art climate models contain, to a significant degree, empirical components. In particular, subgrid-scale (SGS) parameterizations are usually highly tuned against observations or high-resolution model data. While this enables the models to minimize the error during hindcasts, it is not guaranteed that it yields a benefit for climate projections because of climate change. In this thesis the Fluctuation-Dissipation theorem (FDT) is used to update the statistics of the system in the presence of an external forcing. If the empirical parameters are tuned objectively to the data (i.e., they depend on the statistics of the data), then they might be updated with the FDT. This ansatz is tested within a framework of a semi-empirical model (SEM) based on the leading variance patterns of a quasigeostrophic three-layer model (QG3LM) and supplemented by a purely data-driven parameterization. We show that the FDT is able to successfully update the tuning parameters of the data-driven SGS closure, resulting in a systematic improvement in model performance in comparison to an untreated SEM. Ideally, SGS parameterizations should contain little to no tuning parameters. Thus, complementary to the FDT approach we investigate a stochastic SGS closure constrained by first principles that is calculated using the stochastic mode reduction (SMR). The SMR allows for an analytic derivation of the SGS closure from the model equations while requiring only minimal tuning. We successfully apply the SMR to the QG3LM and construct the reduced stochastic model (RSM). Furthermore, we show that the RSM is more robust against an external forcing than the SEM. Additionally, we find that, under appropriate conditions, the FDT is able to update the empirical parts of the RSM. Yet, only for the response in mean streamfunction the RSM provides useful results, while the response in covariance of the streamfunction is incorrect for most cases. Nevertheless, we obtain a remarkably accurate response in both moments for the RSM in an idealized setting. In combination with the results of the FDT study this indicates that the considered RSM is too low dimensional and encourages us to investigate the response of larger RSMs in the future.
Carbon is an element that controls planetary habitability, and is fundamental for life on Earth. Its behaviour has important consequences for the global climate system, the origin and evolution of life on Earth. While the biosphere and atmosphere’s carbon cycle only accounts for less than 1% of the global carbon budget, hidden reservoirs of deep carbon in the Earth’s interior comprise the predominant storage of carbon on the planet. At the Earth’s surface, 60-70 % of carbon is hosted by carbonate minerals, which are then transported to the Earth’s interior, mainly in the form of sediments, by subduction of the oceanic lithosphere. Subducting plates are subjected to decarbonation, dehydration, and melting with CO2 release via supra-subduction volcanism. Nevertheless, part of the subducted carbonates’ may survive and be further transported to the deep mantle. Direct evidence of the existence of carbonates in the Earth’s interior, possibly reaching down to the lower mantle, comes from the finding of syngenetic inclusions of carbonates in diamonds and mantle xenoliths. The presence of carbonates in the deep Earth has a critical effect on the physical properties of the mantle. Melting and chemical speciation of the mantle are strongly affected by the form of C and carbonate stability. Therefore, the study of the stability and physical properties of carbonates at high pressures and temperatures is fundamental, because understanding the processes involved in the deep carbon cycle helps to improve our picture of the whole mantle.
The systematic characterization of the elastic properties of carbonates as a function of their structure and chemical composition is of great importance because it may allow to identify their presence and distribution by seismology. Inverting seismic observations to successfully constrain the chemical composition and mineralogy of the Earth’s interior requires knowledge of the physical properties of all possible Earth’s materials at pressures and temperatures applicable to the Earth’s interior. Up to now, a multitude of studies has focused on the construction of phase diagrams and structural transitions by means of X-ray diffraction and vibrational spectroscopy experiments.
Few studies are available on the complete elastic tensor of carbonates, however most of the datasets are not accompanied by an accurate characterization of the samples, which are often solid solutions and the exact chemical composition, density or the details about the experimental methods used are not presented. The aim of this thesis is to study the effect of chemical composition on the elastic properties of carbonates, providing a reliable dataset on the elasticity of the main carbonates. In particular, the elastic properties of crystalline aragonite, CaCO3, and Fe-dolomite, (Ca, Mg, Fe)(CO3)2, with different compositions were studied by Brillouin spectroscopy at ambient conditions. Brillouin spectroscopy was also used to investigate the elastic behaviour of amorphous calcium carbonate samples with different water contents (up to 18 wt%) at high pressures, up to 20 GPa.
Furthermore, the importance of cationic substitution on the structure and high pressure behaviour of carbonates was investigated by studying a synthetic CaCO3-SrCO3 solid solution at ambient conditions and at high pressures, up to 10 GPa, by single crystal X-ray diffraction. Finally, the study of the effect of composition on the elastic properties of families of isostructural solids was also extended to a different class of materials, the metal guanidinium formates. The elasticity of a family of perovskite metal organic frameworks, metal guanidinium formates C(NH2)3MII(HCOO)3, with MII =Mn, Zn, Cu, Co, Cd and Ca was investigated by combining Brillouin spectroscopy, resonant ultrasound spectroscopy, density functional theory and thermal diffuse scattering analysis.
Mistral and Tramontane are wind systems in southern France and the western Mediterranean Sea. Both are caused by similar synoptic situations and channeled in valleys. Their relevance for the climate of the western Mediterranean region motivated this work. The representation of Mistral and Tramontane in regional climate simulations was surveyed with the models ALADIN, WRF, PROMES, COSMO-CLM, RegCM, and LMDZ. ERA-Interim and global CMIP5 simulations (MPI-ESM, CMCC-CM, HadGEM2-ES, and CNRM-CM5) provided the lateral boundary data for the regional simulations regarding the 20th century and two representative concentration pathways for the 21st century (RCP4.5 and RCP8.5).
A Mistral and Tramontane time series, a principal component analysis of pressure fields, and a Bayesian network were combined to develop a classification algorithm to identify pressure patterns in favor of Mistral and Tramontane. The regional climate models were able to reproduce the observed climatology of Mistral and Tramontane. Compared to observational data (SAFRAN and QuikSCAT), the simulations underestimate the wind speed over the Mediterranean Sea, mainly at the borders of the main flow. Simulations with smaller grid spacing showed better agreement with the observations.
A sensitivity study tested the influence of the Charnock parameter on the Mistral wind field. Its value impacted both wind speed and wind direction. Decreasing the orographic resolution in idealized simulations using COSMO-CLM caused a reduction in wind speed and a broader flow area. Including a parameterization for subgrid scale orography improved the simulation. However, an accurate simulation of Mistral and Tramontane still requires a high-resolution orography.
The classification algorithm also was applied to pressure fields from regional climate simulations driven by global simulation data. At the end of the 21st century, only small, non-significant changes in the number of Mistral days per year occur in the projection simulations. The number of Tramontane days per year decreased significantly.
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.
Mongolia covers a huge area in Asia and provides excellent Palaeozoic successions although large regions still lack detailed information on fauna and flora in Palaeozoic rocks. Of special interest is the Central Asian Orogenic Belt (CAOB), one of the largest collisional complexes on Earth. The CAOB is composed of a large number of terranes, continental margins, island arcs, backarc/forearc basins and accretionary wedges (Badarch et al. 2002; Safonova et al. 2017). Many terranes and regions underwent strong metamorphism, such as areas north of the Main Mongolian Lineament. To the south, Palaeozoic rocks exhibit low-grade metamorphism, thrusting and folding. Two promising long successions of Palaeozoic rocks were studied, namely the Hushoot Shiveetiin gol section and the Bayankhoshuu Ruins section. The intention of the PhD was to study marine facies settings in Palaeozoic rocks of southern and southwestern Mongolia, in an area little is known in terms of biostratigraphy and events. In order to get a better understanding on events and what might have been the driving forces I studied sections in Mongolia which have not been in the focus of research in the last decades. In order to complement studies on Late Devonian events elsewhere, I decided to study sections in open ocean environments (CAOB), far away from mainly studied epicontinental areas.As stated in many publications, events are the driving force for evolution. They exhibit dramatical changes in the palaeontological record of organisms and they are often associated with dramatic extinctions (Walliser 1996) and anoxic sediments (but not always). Mass extinctions are episodes in which a large number of plant and animal species became extinct within a few thousand to a hundred thousand years. Most events (first and second order events) are traceable worldwide. For instance, in the Late Devonian, mass extinction events recognized at the Frasnian–Famennian (F/F) stage boundary and at the Devonian–Carboniferous (D/C) boundary (McLaren and Goodfellow 1990; Sepkoski 1996; Walliser 1996). Overall, 19% of all families, 50% of all genera and at least 70% of all species became extinct (Raup and Sepkoski 1982; Sepkoski 1996; McGhee et al. 2013), but extinctions also concerned palaeoecosystems and due to Late Devonian events whole coral/stromatoporoid reef ecosystems became extinct. What might be the reason(s) for these dramatical changes? Why did the carbonate factory brake down? There are several reasons which have been discussed. For instance, extraterrestrial bolide impacts (McLaren 1970), anoxia within the water column due to climate changes (House 1985; Becker and House 1994; Caplan and Bustin 1999; Bond and Wignall 2005), transgressions and regressions (Newell 1967; Hallam and Wignall 1999; Purdy 2008; Ruban 2010, 2013; Smith and Benson 2013), eutrophication and increased sediment transport (Joachimski et al. 1993; Schobben et al. 2016), explosive volcanism (Paschall et al. 2019) and/ or large igneous provinces (LIP’s, Ernst et al. 2019; Racki et al. 2020, among others) and much more has been considered. The main problem is that previous studies preferably have been done along former epicontinental margins between Laurassia and Gondwana. Less information is available in deep open oceanic successions and shallow-water areas around island arcs. Generally, events are characterized by bituminous rocks, such as black limestones and shales within marine realms, but they do not necessarily occur everywhere (see Carmichael et al. 2016). So, the question was: “Can we recognize Palaeozoic event layers or equivalents in the CAOB?” If so, is this comparable to already published data? What might be the main trigger, as there are many reasons still in discussion? And finally, what has happened in the aftermath of those events as it seems that the CAOB may have acted as a refugium for at least some groups, such as for crinoids (Waters and Webster 2009; Tolokonnikova and Ernst 2010). In this thesis, I will provide data from isolated ecosystems from a rather unstudied region, which will complement studies from other parts of the World.
This Ph.D. thesis demonstrates i) the highly precise performance of refined and new analytical setups for clumped isotope analysis (Δ47 and Δ48) and ii) the applicability of clumped isotope analyses to biogenic and abiogenic carbonated apatite (Δ47) and abiogenic carbonates (Δ47 and Δ48) for research related to paleothermophysiology and paleoclimatology, whereas the overall analytical precision has been increased.
A comprehensive Δ47 dataset with 122 replicate analyses is provided from which the temperature dependence of Δ47 for (bio)apatite (Δ47-1/T2) is calculated between 1 °C and 80 °C. The temperature dependence of oxygen isotope equilibrium fractionation between carbonated synthetic apatite and water (1,000ln(αCHAP-H2O)) is experimentally determined. When applied to tooth enameloid from a modern Greenland shark (Somniosus microcephalus), a Late Miocene megatooth shark (Carcharodon megalodon), and an Upper Cretaceous Tyrannosaurus rex, reconstructed Δ47-based temperatures and δ18OH2O are in line with previously published data.
An analytical setup for highly precise clumped isotope analysis is described that allows for the simultaneous measurement of ∆47 and ∆48 in CO2 with external reproducibilities close to the respective shot-noise limits. The analyte gases originate from pure carbonates that were digested in hypersaturated orthophosphoric acid and purified using a fully automated device. Δ47 data sets with 117 replicate analyses in total on 22 pedogenic carbonate nodules from two Spanish Middle Miocene sections reveal the continental Southern European thermal structure during the end of the Middle Miocene Climatic Optimum (MCO) and the complete Middle Miocene Climatic Transition (MMCT; from 15.33 to 12.98 Ma).
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.
The overarching goal of the thesis was to create a holistic predictive framework, a vegetation model, by improving the representations of and interactions between the biosphere, hydrosphere, atmosphere and pedosphere. Vegetation models rep- resent a crucial component of Earth system model since the properties of the land surface, via interactions with the atmosphere, can have extremely large climatic effect. Yet, there remains great uncertainty associated with the dynamics of the vegetated land surface. Various vegetation models have been critiqued for numerous reasons including overly simplistic representations of vegetation, prescribed vegetation, poor representations of diversity, inaccurate representations of competition, non-transparent model calibration, and poor responses to drought. The purpose of the creation of this "next generation" model was to address deficiencies common to current vegetation modelling paradigms.
The representation of the biosphere within this framework was improved via two separate development axes. First, ecological realism was improved by integrating concepts from community assembly theory, co-existence theory, and evolutionary theory. Explicitly, rather than defining teleonomic rules to define plant behaviour the process of natural selection is modelled. By modelling the pro- cess of natural selection and its affect on relative fitness, myriad "rules" which continually adapt to biotic and abiotic conditions "come out" as a consequence of the modelled dynamics rather than being "put in". In aDGVM2 (adaptive Dynamic Global Vegetation model 2) communities of plants and their trait values evolve through time, this evolution is constrained by trade-offs between traits. Poorly performing individuals are more likely to die and produce fewer copies of themselves, this results in a filtering of trait values. Further, the community and species’ trait values can evolve through successive generations via reproduction, mutation and crossover which we approximate by using a genetic optimisation algorithm. Thus, a plant community consisting of individuals and species with potentially novel and diverse trait values is assembled iteratively through time.
We tested the assertion that improved integration of concepts from community assembly, evolutionary, and co-existence theory could address limitations of DGVMs in Chapter 2. We demonstrated that such an approach does indeed allow diverse communities of plants to emerge from the modelling framework. We showed that the position of the emergent communities in trait space differed along abiotic gradients and that, in simulations where reproductive isolation was simulated, communities emerged which were composed of multiple co-existing clusters in trait-space. Simulated trait values of co-existing strategies emerging from aDGVM2 were often multimodal, indicative of the emergence of multiple life- history strategies.
Second, to successfully model how natural selection forms a community requires accurate representation of how resource availability affects fitness. In the majority of dynamic global vegetation models (DGVMs) there is no real representation of plant hydraulics with plant water availability being calculated as a simple function of relative soil moisture content and root fractions across a number of soil layers. Worryingly, a number vegetation models appear to under represent the magnitude of these observed responses to drought. This was deficiency was ad- dressed in Chapter 3 by designing a simplified version of the cohesion tension theory of sapwood ascent where elements determining plant conductances are considered in series and implementing a set of trait trade-offs which influence a plant’s hydraulic strategy whereby hydraulic safety trades-off against xylem and leaf conductivity.
Interactions between the biosphere, pedosphere, and hydrosphere can also potentially mediate water resource availability and thus fitness. In the majority of DGVMs the volume of soil explored and explorable by plant roots in fixed glob- ally and usually constrained to a depth not greater than 3m. However, we know that soils can have a strong effect on vegetation distributions, that soil depth is not constant globally, and that plants root to variable depths.
In Chapter 4 I explored interactions between soil depth, plant rooting and the emergent properties of communities and highlighted the importance of considering interactions between the biosphere, hydrosphere, pedosphere, and fire. Here I demonstrated that, in addition to fire and precipitation, edaphic constraints on the volume of soil explorable by plant roots (e.g. by shallow soils, lateritic layers, anoxic conditions due to water logging, toxicity resulting from heavy metal concentrations) can affect the process of plant community assembly, alter the mean values of multiple traits in communities, and the trait diversity of communities.
...
Climate controls the broad-scale distribution of vegetation and change in climate will alter the vegetation distribution, biome boundaries, biodiversity, phenology and supply of ecosystem services. A better understanding of the consequences of climate change is required, particularly in under-investigated regions such as tropical Asia, i.e., South and South-east Asia, which is a host to 7 of the 36 global biodiversity hotspots. Conservation strategies would also require an in-depth understanding of the response of vegetation to climate change. Therefore, the main objective of this thesis was to investigate the impact of climate change and rising CO2 vegetation in tropical Asia. Dynamic global vegetation model (DGVMs) are the well-known tools to investigate vegetation-climate interactions and climate change impacts on ecosystems. In this thesis, I used a complex trait-based DGVM called adaptive dynamic vegetation model version 2 (aDGVM2).
In Chapter 1, I presented a brief background of the phytogeography and discussed the exiting knowledge gap on vegetation-climate interactions in the region. One major disadvantage for available DGVMs studies for the tropical Asia is that most of them have used fixed plant functional types (PFTs) and do not explicitly represent the distinct varieties of vegetation type of the region such as Asian savannas. In Chapter 2, I discussed at great length to improve DGVMs for South Asia and discussed ways to include them in the model for better representation of region vegetation-climate interaction.
I upgraded the current version of aDGVM2 and added a new vegetation type i.e., C3 grasses, and modified the sub-module to simulate photosynthesis for each individual plants to aDGVM2. In chapter 3, I used this updated version of aDGVM2 to simulate the current and future vegetation distribution in South Asia under RCP4.5 and RCP8.5 (RCP: representative concentration pathway). The model predicted an increase in biomass, canopy cover, and tree height under the presence of CO2 fertilization, which triggered transitions towards tree-dominated biomes by the end of the 21st century under both RCPs. I found that vegetation along the Western Ghats and the Himalayas are more susceptible to change due to climate change and open biomes such as grassland and savanna are prone to woody encroachment.
In Chapter 4, the study domain was extended to include South-east Asia to verify if the model configuration used in Chapter 3 can also simulate vegetation patterns in tropical Asia. The aDGVM2 simulations showed a robust trend of increasing vegetation biomass and transitions from small deciduous vegetation to taller evergreen vegetation across most of tropical Asia. Shifts in plant phenology also affect ecosystem carbon cycles and ecosystem feedback to climate, yet the quantification of such impacts remains challenging. The study showed increased biomass due to CO2 fertilization, indicates that the region can remain a carbon sink given there is no other resource limitation. However, nutrient limitations on CO2 fertilization effects were not included in the study, and carbon sink potential has to be seen with caution.
In Chapter 5, I focused on Asian savannas, which have been mismanaged since the colonial era due to misinterpretation as a degraded forest. I proposed a biome classification scheme to distinguish between degraded forest or woodland and savanna based on the abundance of grass biomass and canopy cover. I found that considering vegetation systems as woodland or degraded forest could easily be mistaken as a potential for forest restoration within a tree-centric perspective. This would put approximately 35% to 40% of a unique savanna biome at risk. Although projected woody encroachments may imply a transition toward the forest that benefits climate mitigation. This raises potential conflicts of interest between biodiversity conservation in open ecosystems, i.e., savanna and active afforestation, to enhance carbon sequestration. Proper management strategies should be taken into account to maintain a balance for both objective
In conclusion, the model predicted that vegetation in South and South-East Asia would significantly shift towards tree-dominated biomes due to CO2-induced fertilization of C3-photosynthesis. The simulation under fixed CO2 and rising CO2 scenarios clearly showed that rising level of atmospheric CO2 is responsible for most of the predicted change in biome properties. This study is an important step towards understanding ecosystems of South and Southeast Asia, specifically savannas. The aDGVM2 can serve as tools to inform decision making for climate adaptation and mitigation for savanna. The thesis, thus contributes to our ability to improve conservation strategies to mitigate the consequences of climate change.
The formation of terrestrial planets was a complex process which begun in the very early stage of the Solar System in the protoplanetary disk (PPD). Chondrites are fragments of planet precursors, which have never experienced differentiation and can help to reconstruct the first processes leading to planet formation. The main components of chondrites are chondrules, calcium-aluminum-rich inclusions (CAIs), amoeboid olivine aggregates (AOAs), metals and fine-grained material. Each of these components formed by a complex mechanism involving aggregation and/or melting. Previous research has already provided an overall view of the formation of these objects, however, there are still open questions regarding the aggregation behavior of particles, the heating mechanism(s) and the thermal history of CAIs, AOAs and chondrules. For instance, the involvement of flash-heating events and electrostatics in the aggregation and melting of these objects has been a keen topic of discussion.
The aim of this doctoral thesis was to develop and carry out an experiment to study various early Solar System processes under long-term microgravity. In the project with the acronym EXCISS (Experimental Chondrule Formation aboard the ISS), free-floating, 126(23)µm-sized Mg2SiO4 dust particles were exposed to electric fields and electric discharges.
The experimental set-up was installed inside a 10x10x15 cm3-sized container and consisted of an arc generation unit connected to the sample chamber, a camera with an optical system, a power supply unit with lithium-ion batteries and the EXCISS mainboard with a Raspberry Pi Zero and mass storage devices. The sample chamber was manufactured from quartz glass and the experiments were filmed. The complete experiment container was subsequently returned to the Goethe University and the samples were analyzed with scanning electron microscopy, electron backscatter diffraction and synchrotron micro-CT.
Video analysis has shown that particles, which were agitated by electric discharges, align in chains within the electric field with their longest axis parallel to the electric field lines. Consequently, electric fields could have influenced the inner structure and porosity of particle aggregates in the PPD.
The discharge experiments produced fused aggregates and individual melt spherules.
The fused aggregates share many morphological characteristics with natural fluffy-type CAIs and some igneous CAIs found in chondrites. Consequently, CAIs could have formed by the aggregation of particles with various degrees of melting. Further, a small amount of melting could have supplied the required stability for such fractal structures to have survived transportation and aggregation to, and subsequent compaction within, developing planetesimals.
Some initial particles were completely melted by the arc discharges and formed melt spherules. The newly formed olivines crystallized with a preferred orientation of the [010] axis perpendicular to the surface of the spherule. Similar preferred orientations have been found in natural chondrules. However, the microstructure differs from the results of previous experiments on Earth, which show, for example, crystal settling on one side of the sample because of the influence of gravity. Furthermore, the melt spherules show evidence for an interaction of the melt with the surrounding hot gas. Therefore, microgravity experiments with more advanced experimental parameters bear great potential for future chondrule formation experiments.
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.
In light of the global sea-level rise and climate change of the 21th century, it is important to look back into the recent past in order to understand what the future might hold. A multi-proxy data set was compiled to evaluate the influence of geomorphological and environmental factors, such as antecedent topography, subsidence, sea level and climate, on reef, sand apron and lagoon development in modern carbonate platforms through the Holocene. Therefore, a combination of remote sensing and morphological data from 122 modern carbonate platforms and atolls in the Atlantic, Indian and Pacific Oceans were conducted, along with a case study from the oceanic (Darwinian) barrier-reef system of Bora Bora, French Polynesia, South Pacific.
The influence of antecedent topography and platform size as factors controlling Holocene sand apron development and extension in modern atolls and carbonate platforms is hypothesized. Antecedent topography describes the elevation and relief of the underlying Pleistocene topography (karst) and determines the distance from the sea floor to the rising postglacial sea level. Maximum lagoon depth and marginal reef thickness, when available in literature, were used as proxies for antecedent topography. Sand apron proportions of 122 atolls and carbonate platforms from the Atlantic, Indian and Pacific Oceans were quantified and correlated to maximum lagoon depth, total platform area and marginal reef thickness. This study shows that sand apron proportions increase with decreasing lagoon depths. Sand apron proportions also increase with decreasing platform area. The interaction of antecedent topography and Holocene sea-level rise is responsible for variations in accommodation space and at least determines the extension of the lateral expansion of sand aprons. In general, sand apron formation started when marginal reefs approached relative sea level. Spatial and regional variations in sea-level history let sand apron formation start earlier in the Indo-Pacific region (transgressive-regressive) than in the Western Atlantic Ocean (transgressive).
The influence of sea level, antecedent topography and subsidence of a volcanic island on late Quaternary reef development was evaluated based on six rotary core transects on the barrier and fringing reefs of Bora Bora. This study was designed to revalue the Darwinian model, the subsidence theory of reef development, which genetically connects fringing reef, barrier reef and atoll development by continuous subsidence of the volcanic basement. Postglacial sea-level rise, and to a minor degree subsidence, were identified as major factors controlling Holocene reef development in that they have created accommodation space and controlled reef architecture. Antecedent topography was also an important factor because the Holocene barrier reef is located on a Pleistocene barrier reef forming a topographic high. Pleistocene soil and basalt formed the pedestal of the fringing reef. Uranium-Thorium dating shows that barrier and fringing reefs developed contemporaneously during the Holocene.
In the barrier–reef lagoon of Bora Bora, the influence of environmental factors, such as sea level and climate, tsunamis and tropical cyclones controlling Holocene sediment dynamics was evaluated based on sedimentological, paleontological, geochronological and geochemical data. The lagoonal succession comprises mixed carbonate-siliciclastic sediments overlying peat and Pleistocene soil. The multi-proxy data set shows variations in grain-size, total organic carbon (proxy for primary productivity), Ca and Cl element intensities (proxies for carbonate availability and lagoonal salinity) during the mid-late Holocene. These patterns could result from event sedimentation during storms and correlate to event deposits found in nearby Tahaa, probably induced by elevated cyclone activity. Accordingly, elevated erosion and runoff from the volcanic island and lower lagoonal salinity would be a result of rainfall during repeated cyclone landfall. However, Ti/Ca and Fe/Ca ratios as proxies for terrigenous sediment delivery peaked out in the early Holocene and declined since the mid-Holocene. Benthic foraminifera assemblages do not indicate reef-to-lagoon transport. Alternatively, higher and sustained hydrodynamic energy is probably induced by stronger trade winds and a higher-than-present sea level during the mid-late Holocene. The increase in mid-late Holocene sediment dynamics within the back-reef lagoon is supposed to display sediment-load shedding of sand aprons due to the oversteepening of slopes at sand apron/lagoon edges during their progradation rather than an increase in tropical storm activity during that time.
The influence of sea-level and climate changes on sediment import, composition and distribution in the Bora Bora lagoon during the Holocene is validated. Lagoonal facies succession comprises siderite-rich marly wackestones, foraminifera-siderite wackestones, mollusk-foraminifera marly packstones and mollusk-rich wackestones during the early-mid Holocene, and mudstones since the mid-late Holocene. During the early Holocene, enhanced weathering and iron input from the volcanic island due to wetter climate conditions led to the formation of siderite within the lagoonal sediments. The geochemical composition of these siderites shows that precipitation was driven by microbial activity and iron reduction in the presence of dissolved bicarbonate. Chemical substitutions at grain margins illustrate changes in the oxidation state and probably reflect changes in pore water chemistry due to sea-level rise and climate change (rainfall). In the late Holocene, sediment transport into the lagoon is hampered by motus on the windward side of the lagoon, which led to early submarine lithification within the lagoon.
Hydro-climatic causes of widespread floods in central Europe : on rain-on-snow and Vb-cyclone events
(2021)
The presented work investigates the hydro-meteorological and hydro-climatological drivers of widespread floods in Central Europe during the past century. Due to the strong seasonality of the detected flood drivers, the thesis is divided into two parts: the first part focuses on widespread winter floods and the second one on extreme summer floods. For analysing past flood events, we profited from the dynamically downscaled centennial ERA-20C reanalysis (continuously from 1901—2010). The downscaling was performed over Europe with a coupled regional atmosphere-ocean model (COSMO-CLM+NEMO) to represent the water cycle more realistic. These high resolution atmospheric data allowed us to study the four-dimensional atmospheric state during selected floods during the early decades of the 20th century for the first time with such a high temporal and spatial resolution.
During the winter half-year, the observed floods were particularly widespread. High peak discharges were recorded simultaneously in the Rhine, Elbe, and Danube catchments. Most of these trans-basin floods were compound events caused by rainfall during extensive snowmelt (i.e., rain-on-snow events). Interestingly, the winter flood time series exhibited a remarkable high flood frequency during the 1940s and 1980s, while other decades were flood-poor. We detected a synchronization of the inter-annual flood frequency with the superposition of the North Atlantic Oscillation (NAO) and the Scandinavian pattern (SCA). The negative NAO phase is often associated with large snowfall and cyclone tracks over southern Europe, while the negative SCA pattern correlates with total precipitation in the affected river catchments.
During the summer half-year, most extreme floods in Central Europe were caused by so-called Vb-cyclones propagating from the Mediterranean Sea north-eastward to Central Europe. So far in the literature, only a few Vb-events, which occurred during the past two decades, have been analysed. We extended the previous case studies by several past Vb-cyclone floods since 1900. We investigated the processes that intensify Vb-cyclone precipitation with Lagrangian moisture-source diagnostics and the parametric transfer entropy measure TE-linear. Overall, an enhanced and dynamically driven moisture uptake over the Mediterranean Sea was found to be characteristic for Vb-events with heavy precipitation. This is supported by high information exchange from evaporation over the western basin of the Mediterranean Sea towards heavy precipitation in the Odra catchment. The dominating moisture uptake regions during the investigated events were, however, the European continent and the North Sea. A possible cause could be the pre-moistening of non-saturated continental moisture sources upstream of the affected river catchments as indicated by significant information exchange from land surface evaporation and soil moisture content along the Vb-cyclone pathway. Besides, evaporation over the Mediterranean Sea might contribute to Vb-cyclone intensification in the early stages of their development through latent heat release. On the catchment scale, orographic rainfall and convective precipitation further enhance the flood triggering rainfall. As expected, the Vb-cyclones mainly trigger precipitation along west-east orientated mountain ranges such as the Alps or Ore mountains due to their meridional pathway. Remarkably, during summer, we detected a convective fraction of up to 90% during the afternoons of individual days and up to 23% on average (based on convective cell tracking and convection-permitting simulations of selected flood events since 1900).
The presented analyses deepened the knowledge on atmospheric and hydroclimatic drivers of widespread floods in Central Europe. This will serve as a basis for future studies on the predictability of floods induced by rain-on-snow and Vb-cyclone precipitation events in the context of a changing climate.
The analysis of the global stratospheric meridional circulation, known as the Brewer-Dobson circulation, is an essential part of both experimental and theoretical atmospheric sciences. This large-scale circulation has a crucial influence on the global burden of greenhouse gases and ozone depleting substances throughout the complete atmosphere. This makes it an important factor for the Earth’s radiative budget, which is perceptible at the Earth’s surface despite the remote location of the stratosphere. In the course of climate change it is generally expected that also the Brewer-Dobson circulation undergoes significant changes in structure and strength, although the exact repercussions are still uncertain and thus remain an open scientific question. A general problem for the observational investigation of the dynamical processes in the stratosphere is that residual mean transport cannot be measured directly and hence requires the use of sophisticated proxies. Many studies in the past consider the so-called mean age of air, which is a measure of the average time an air parcel has spent in the stratosphere since passing a certain reference point. While changes in the strength and structure can be detected and visualized using mean age of air, a more thorough distinction between the different involved transport mechanisms of the circulation (residual circulation, mixing) cannot be made. For that, consideration of a full distribution of all relevant transit times through the stratosphere, an age spectrum, is favorable and a powerful tool to analyze the spatial structure as well as possible future changes in detail. Mean age of air and age spectra can be readily derived in atmospheric modeling studies, but an observationally based retrieval is challenging. Mean age of air is usually approximated from measurements of very long-lived trace gas species that act as a dynamical tracer for the stratosphere. The retrieval of age spectra from observations, however, remains an open task for which different methods have been proposed in the past, that often require a combination of strong assumptions and model data explicitly. This is a major issue for a precise and independent investigation of stratospheric dynamics based on measurements. The focus of this cumulative dissertation is on the development process and application of an inversion method to derive stratospheric age spectra from mixing ratios of chemically active substances that combines an applicable and precise ansatz with a minimized amount of necessary model data. Chemically active species have the important benefit that chemistry and transport in the stratosphere are strongly correlated so that the state of depletion of a trace gas can give some information on certain parts of the age spectrum. Considering a sufficient number of distinct trace gases simultaneously, a full approximation of the age spectrum should be possible. The main section of this thesis is split into three parts, which follow the main aspects and key results of the three publications involved (Hauck et al., 2019, 2020; Keber et al., 2020). The newly developed inverse method is based upon the previously established ansatz by Schoeberl et al. (2005), but constrains the shape of the age spectrum by a single parameter inverse Gaussian function. This keeps the balance between applicability and accuracy with a limited amount of measurement data. Additionally, the method introduces a seasonal scaling factor that imposes higher order maxima and minima onto the intrinsically monomodal spectrum based on the seasonal cycle of the tropical upward mass flux to incorporate phases of weaker and stronger transport. A proof of concept of the inverse method is provided using an idealized simulation of the ECHAM/MESSy Atmospheric Chemistry (EMAC) model, where the method is applied to a set of artificial radioactive trace gases with known chemical lifetime. The results imply that the method works properly and retrieves age spectra that match the EMAC reference spectra significantly well on the global and seasonal scale. Only in the lower stratosphere, the performance of the inverse method on the seasonal scale decreases as entrainment into the stratosphere is considered only across the tropical tropopause. Transport across the local extratropical tropopause, however, is a key feature for trace gases in the extratropical lowermost stratosphere so that this entrainment must be included explicitly.
In the second part, the discovered problems are approached to make the inverse method applicable to observations. The formulation of the method is extended to incorporate transport explicitly across the tropical (30° S – 30° N), northern extratropical (30° N – 90° N), and southern extratropical tropopause (30° S – 90° S) each with a single age spectrum that can be inverted independently.
...
Historic amphibian settlements in the northwestern Nile delta - a geoarchaeological perspective
(2020)
No concise picture of the archaeological and palaeoecological evolution can be drawn for the northwestern Nile delta, and archaeological records show significant population dynamics that still need explanation and spur the need for further palaeoenvironmental research. This study delivers a set of new methods especially in the fields of remote sensing and data analytics that can be regarded as important milestones and foundations for further palaeoenvironmental research in the area. Additionally, it shows new insights for individual time slices.
This geoarchaeological project is a cooperation with the archaeological excavations of the German Archaeological Institute (DAI) in Buto and Kom el’Gir. It expands the work of Wunderlich (1989) which laid important foundations in understanding the origin of the initial landscape that was later colonized in different cultural stages showing different dynamics, settlement intensities and even long phases of abandonment or breaks in between. This forms the starting point for relating the population dynamics of the different cultural phases reaching from Predynastic (prior to 3150 before Christ) up to the Greco-Roman era (~anno Domini 650) to the environmental history and events that occurred in the area. It is very likely that environmental changes such as the shifting of major water routes, inundation or paludification of larger areas or other environmental events affected settlements and human life in the area.
In the fields of remote sensing new methods are presented to complete information on the location of ancient settlements, and complex workflows are developed that allow the tracing of subsurface structures via indirect analysis of vegetation growth in larger time series data. It was verified that a relationship exists between vegetation performance, the appearance of archaeologic material in the topsoil, and the location of former Nile river branches.
Together with a new high resolution digital elevation model (DEM) based on TanDEM-X data, new interpretations with a high spatial significance are possible. For individual time slices, namely the Late Dynastic and Greco-Roman era, this work delivers a detailed landscape description suggesting a finely ramified subdelta, with all settlements placed on alluvial levees. This explains the massive increase in settlements in the Ptolemaic, Roman and in particular late Roman periods (4th century before Christ – anno Domini 7th century).
We sampled the Nile delta clays together with the channels and the material of the archaeologic excavations in vibracores and profile walls. This geologic inspection of the subsurface together with geochemical results from a handheld portable X-ray fluorescence device (pXRF) allowed new interpretations of the landscape and environmental history. For example, we used geochemical data to distinguish between artificial and natural channels as a measure for the anthropogenic influence, a proxy for past environmental characteristics and lastly as a basis for a new dating method. Many of the channels, for instance, were dated by our own 14C datings, comparisons with the previous work ofWunderlich (1989) and application of new dating approach based on machine learning with artificial neural networks. Additionally, we run a full methodological approach, and examine the applicability of pXRF methods in general, and test the quality of the data to detect distinct geochemical differences between the main settlement phases with advanced methods in data analytics. The dating is based, for example, on the training of artificial neural networks with pXRF data from archaeological material of well-dated context to date test data of cultural layers within the vibracores. With this method the homogeneous Nile alluvium, cultural layers and channels can be dated roughly and, as a result, fundamental changes in the landscape can be linked with the settlement history of Buto and neighboring tells.
Shrubs are a characteristic component of savannas, where they coexist with trees and grasses. They are often part of woody encroachment phenomena, which have been observed globally, and the determinant of shrub encroachment cases, which are particularly of concern in African savannas. In response to climate change and land use change, African savannas are vulnerable to biome shifts and shrub encroachment is a process driving and explaining this risk.
We contribute to furthering the understanding of shrubs biogeography and ecology by considering the number of stems of woody plants to characterise shrubs phenotype and strategy. We postulate that shrubs are multi-stemmed, compared to single-stemmed trees and integrate this assumption in aDGVM2 (adaptive Dynamic Global Vegetation Model 2). Modelling a trait representing the number of stems of a woody plant implies a trade-off between single-stemmed plants having higher height growth potential and multi-stemmed plants having higher hydraulic capacity but limited height growth. Multi-stemmed individuals, being shorter, are more likely to suffer severe damage from fires than tall single-stemmed trees managing to grow their crown out of the flame zone.
We simulate potential vegetation over sub-Saharan Africa at 1° spatial resolution, with aDGVM2 and compare it to simulations without our shrub model turned on. We also test the impact of fire by including or excluding it from our simulations. To assess the accuracy and relevance of our approach, we benchmark our overall model’s performance against multiple satellite derived products of above ground biomass (AGBM), and against specific field measurements of AGBM. We further benchmark our results against vegetation cover type derived from satellite data.
We demonstrate that shrubs can be modelled as multi-stemmed woody plants in African savannas based on whole-plant trait trade-off without being predefined as static functional types. Indeed, the addition of our shrub model to aDGVM2 allows for shrubs to emerge dynamically through community assembly processes without a priori categorisation. Our shrub model also improves the simulated vegetation patterns simulated by aDGVM2 in sub-Saharan Africa, particularly in savannas. The simulated pattern of stem number per woody individual broadly follows our assumptions about biogeographic patterns as it is lowest in equatorial African forests and increases in savannas and grasslands as precipitation decreases. Shrubs are more abundant in more water-stressed regions where they have a competitive advantage over trees due to their increased relative water transport potential. However, in arid and hyper-arid regions, further investigations are required. Simulated shrub prevalence is higher in more open and fire prone landscapes, where woody cover and biomass are reduced.
Adding shrubs to aDGVM2, while increasing complexity allows for greater simulated diversity. As resilience and resistance of ecosystems have been shown to be influenced by diversity, such model development is necessary to improve our ability to forecast ecosystems responses to changes. However, there are challenges to fully tap this benefit. Assessing the accuracy and relevance of our approach is challenging. Data and simulations are conceptually different which limit the possibility to conclude based on comparison. Benchmarking challenge is exacerbated by the variability existing among satellite derived products and site studies observations. In areas of extremely low biomass and vegetation cover, such as deserts and semi-deserts, the accuracy of our model is more concerning as small differences in absolute values are relatively more important.
Categorisation of life-forms shapes our understanding of their ecology and biogeography, thus, consensus about their definition is direly needed. To contribute to this debate, we investigate how vegetation distribution patterns arising from our shrub model inform our understanding of shrub biogeography. First, shrub distribution in trait space (considering stem number), relatively to environmental drivers, concurs with our assumptions. Second, shrub spatial distribution is consistent with our characterisation assumptions. Third, the role of simulated shrubs in an ecosystem supports realistic ecological dynamics. Our model allows for, shrubs to exhibit a specific phenotype, but also a specific life-strategy, which we characterise in terms of persistence strategy (shrubs are mainly resprouters, in contrast to trees, which can be either resprouters or reseeders) and in terms of resource acquisition (rooting strategy) and allocation (carbon investment). Adding stem count as a trait to aDGVM2 increase the range of simulated functional diversity.
Our shrub model allows for aDGVM2 to simulate realistic ratio of grass to woody vegetation across sub-saharan Africa. Similarly, it simulates ratio of shrubs to trees consistent with our hypotheses.
...
As a cognitively-mediated response, autonomous adaptation at farm-gate levels constitutes reactionary actions by farmers against climate impacts. These actions are shaped by interacting factors such as household characteristics, livelihood scope and resources. It is driven by the goal of adapting cultivated farmlands to climate and for sustaining crop yields. Thus, interest in balancing adaptation goals with protection of vegetation conditions is less of a priority. Lack of research interest in understanding the gap between objectives of reactionary adaptation and protection of surface conditions (vegetation canopies) is a gap in research. In many studies, farm-gate level adaptation is described as a set of zero-feedback actions in response to climate impacts. This perception conceals the stress and impact-engendering attribute of reactionary adaptation. Inspired towards addressing this conceptual gap; this study investigates impact of farmers’ reactionary adaptation on vegetation cover in Keffi, Nasarawa, Nigeria. A twenty-year time-series NDVI and rainfall datasets are linearly regressed to examine the extent of NDVI-rainfall sensitivity. A weak linear relationship between NDVI and rainfall in Keffi for the period, 1999-2018 is observed. At a regression slope of 0.001, R squared, R2=0.129 (implying that only about 13% of the variability in NDVI in Keffi are explained by rainfall amount) and a bivariate regression coefficient, r=0.359; statistical evidence shows that rainfall amount are not significant predictors of NDVI in Keffi. In investigating the possible interference of non-rainfall factors on vegetation productivity (NDVI) in Keffi; a residual trend (RESTREND) analysis was carried out. Regression of residuals from NDVI-Rainfall linear regression produced a R=0.192 with a negative and downwards slope. The downward character of the RESTREND slope is suggestive of non-rainfall factors contained in the residuals. In validating the RESTREND analysis, a comparative analysis between observed and predicted NDVI derived from a reference NDVI value of 0.46 was carried out. The NDVI value of 0.46, is empirically assumed to be average NDVI value expected at a minimum rainfall amount of 850mm/year reported in tropical Savanna ecosystems. Using this empirical relationship, NDVI values were predicted for Keffi. Even at higher rainfall amounts≈1340mm/year, amounts were unable to produce corresponding higher NDVI values; rather a more plausible correlation between reference-derived predicted NDVI values and rainfall was obtained. A further analysis with predicted NDVI values, based on 1999 NDVI value in Keffi returned higher NDVI units than observed NDVI values. This strengthens the attribution of the possible interference of rainfall-NDVI sensitivity by non-rainfall factors like human activities on vegetation productivity. Surface soil analysis to exclude potential impacts of soil nutrients and moisture deficiency on vegetation productivity, showed that soil had insignificant effect on vegetation dynamics. Further inferential analysis, using the inter-annual NDVI and the reclassified bi-decadal NDVI maps showed that spatial vegetation distribution in Keffi were driven by farmers inter-annual rotational cultivation footprints than rainfall variability. With a three-class categorization, “gain, loss and significant loss”, the spatial distribution of vegetation in Keffi between (1999-2008) and (2009-2018) was assessed. Temporal condition (stressed and healthy) across the three classes supports the attribution of farmers’ reactionary adaptation and cultivation practices on the dynamic spatial vegetation distribution. Between 1999 -2018, an increase in areas with significant vegetation loss (42%), so with a decrease of -25% in areas with healthy vegetation was observed. The character of vegetation cover across the two decadal time slices, reflects landuse intensity and unsustainable farming practices. Preferences for modification of cultivation practices and changes in seed by farmers exerts positive feedbacks on vegetation cover. Higher statistical measures, 38.4% (yearly cropping) and 44% (shifting cultivation with less fallow periods) were observed in the chi-square analysis. These measures were higher than 2.0% relating to shifting cultivation with more fallow periods. While 11.6% farmers noted cultural practices as reasons for preferred cultivation methods, 48.4% farmers attributed climate as reason behind cultivation modification. This was higher than 24.4% who linked issues of tenure rights to cultivation practices. With preferences for yield- breaching strategies, the non-receding cultivation and shorter fallow practices in Keffi triggers feedback on vegetation dynamics. Evidence from this study shows that the NDVI-rainfall functional sensitivity in Keffi is plausibly dampened by effects of reactionary farm-gate level adaptation practices.
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.
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.
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 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.
During my PhD, I was applying the clumped isotope technique to modern brachiopods and fossil belemnites, and I conducted methodological work. Carbonate clumped isotope thermometry is a tool to reconstruct carbonate precipitation temperatures. In contrast to oxygen isotope thermometry, i.e., the δ18O-thermometer, the carbonate clumped isotope thermometer does not require an estimate for the oxygen isotope composition of the seawater, as it considers the fractionation of isotopes exclusively amongst carbonate isotopologues. The ∆47 value of a carbonate expresses the abundance of the 13C–18O bond bearing carbonate isotopologue, within the carbonate, relative to its random distribution. In thermodynamic equilibrium, the ∆47 value of a given carbonate is solely a function of the carbonate precipitation temperature. However, kinetic isotope fractionations, i.e., vital effects, driven by diffusion, pH or incomplete oxygen isotope exchange between water and dissolved inorganic carbonate species can cause the carbonate to be precipitated with isotopic compositions that are offset from those predicted for thermodynamic equilibrium.
Brachiopods serve as important geochemical archives of past climate conditions. To investigate the nature and significance of kinetic controls on brachiopod shell δ18O and ∆47 values, in collaboration with the BASE-LiNE Earth ITN, I analysed the bulk and clumped isotope compositions of eighteen modern brachiopod shells, collected from different geographic locations and water depths that cover a substantial range of growth temperatures. Growth temperatures and seawater δ18O values for each brachiopod were independently determined. Most of the analysed brachiopods exhibit combined offsets from clumped and oxygen isotope equilibrium, and there is a significant negative correlation between the offset values. The observed correlation slope between offset ∆47 and offset δ18O point to the importance of kinetic effects associated with Knudsen diffusion and incomplete hydration and hydroxylation of CO2 (aq), occurring during biomineralisation. The correlations between the growth rates of the analysed brachiopods and both the offset ∆47 and the offset δ18O values provide further arguments for the presence of kinetic effects. In conclusion, the oxygen and clumped isotope composition of modern brachiopod shells are affected by growth rate-induced kinetic effects that hinder their use for palaeoceanography.
One of the most important events in human history occurred during the Early Pleistocene: the dispersal of early hominins out of Africa and into Europe and Asia. In Western Europe, the earliest evidences of the genus Homo have been found in the Baza Basin, at the sites of Orce in the SE of the Iberian Peninsula. These sites contain fossils and lithic industry dated approximately as 1.4–1.3 Ma.While hominin remains and artifacts at Orce, as well as the accompanying fauna, have been extensively studied, the properties and evolution of the Early Pleistocene vegetation in the basin remain unknown. The general effect of climate change on the expansion of early hominins from Africa into Eurasia still remains unclear. It is not known if the Early Pleistocene climate changes and the development of glacials periods led to the extirpation of European communities, or if those communities were able to endure and persist through such adverse climatic periods. This open question highlights the need for climate and environmental analyses for the time before, during and after the first presence of Homo in Europe. This PhD thesis contributes to that need by the presentation of the first long pollen record of the Baza Basin, where the oldest hominin sites in Western Europe are found.
Air-sea feedbacks between the Mediterranean Sea and the atmosphere on various temporal and spatial scales play a major role in the Mediterranean regional climate system and beyond. The Mediterranean Sea is a source of moisture due to excess evaporation and, on a long-term average, is associated with a warming of the lower atmosphere in contact with the sea surface due to heat loss at the air-sea interface. The complex air-sea interactions and feedbacks in the Mediterranean basin strongly modulate the sea surface fluxes and favor several cyclogenetic activities under certain meteorological conditions. Examples of such cyclonic activities are medicanes (Mediterranean hurricanes) and Vb-cyclones. Medicanes are mesoscale, marine, and warm-core Mediterranean cyclones that exhibit some similarities to tropical cyclones, while Vb-cyclones are extra-tropical cyclones, that propagate from the Western Mediterranean Sea and travel across the Eastern European Alps into the Central European region. Extremely strong winds and heavy precipitation associated with these cyclones can lead to severe destruction and flooding. Changes in the intensity and frequency of these cyclones are also projected under changing future climate conditions, where the Mediterranean region has been identified as a hotspot in terms of rising temperatures.
The development of high-resolution regional climate models (RCMs) has progressed our understanding of the processes characterizing the Mediterranean climate. However, large uncertainties still exist regarding the estimates of air-sea fluxes, which, in turn, affect the simulation of the Mediterranean climate. Several factors can be attributed to such discrepancies, such as data quality, temporal and spatial resolution, and the misrepresentation of physical processes. To overcome some of these inconsistencies and deficiencies of the existing climate simulations, a new high-resolution atmosphere-ocean regional coupled model (AORCM) has been developed to simulate the air-sea feedback mechanisms. This coupled model incorporates the coupling of RCM COSMO-CLM (CCLM) and the regional ocean model NEMO-MED12 for the Mediterranean Sea (MED) as well as NEMO-NORDIC for the North- and Baltic Sea (NORDIC). Several experiments were performed using both the coupled and uncoupled models to investigate the impact of air-sea interactions and feedbacks on sea surface heat fluxes, wind speed, and on the formation of Mediterranean cyclones (i.e., medicanes and Vb-cyclones). These experiments were performed using different horizontal atmospheric grid resolutions to analyze the effect of resolution on sea surface heat fluxes, wind speed, and the development of medicanes.
The results of the present study indicate that a finer atmospheric grid resolution ([is as appreciated as]9 vs. [is as appreciated as]50 km) improved the wind speed simulations (particularly near coastal areas) and subsequently improved the simulations of the turbulent heat fluxes. Both parameters were better simulated in the coupled simulations than in the uncoupled simulations, but coupling introduced a warm SST bias in winter. Radiation fluxes were slightly better represented in coarse-grid simulations than in fine-grid simulations. However, the higher-resolution coupled model could reproduce the observed net outgoing total surface heat flux over the Mediterranean Sea. In addition to that sub diurnal SST variations have a strong effect on sub-daily heat fluxes and wind speed but minor effects at longer timescales. Regarding the impact of atmospheric grid resolution ([is as appreciated as]50, 25, and [is as appreciated as]9 km) and ocean coupling on medicanes, it was detected that the coupled model with a finer atmospheric grid ([is as appreciated as]9 km) was able to not only reproduce most medicane events, but also improved the track length, warm core, and wind speed compared to the uncoupled model. The coupled model with the coarse-grid ([is as appreciated as]50 and [is as appreciated as]25 km) did not show any improvement in simulating medicanes compared to the uncoupled model. The spectral nudging technique, applied on the wind components above 850 hPa in the interior domain to keep large-scale circulation close to the driving data (i.e., ERAInterim reanalysis), improved the accuracy of the times and locations of generated medicanes, but no improvement was found in the track length and intensity.
Concerning the role of the Mediterranean Sea coupling on Vb cyclones, the investigation showed that atmosphere-ocean coupling had an overall positive impact, although with a strong case-by-case variation, on the trajectories and intensity of Vb-cyclones as a result of the variation in moisture source for each event. In general, all model configurations could replicate Vbcyclones, their trajectories, and associated precipitation fields. The average structure of the precipitation field was best represented in the coupled simulations. Coupling of the North- and Baltic Seas also showed an improvement in some of the simulated Vb-cyclones.
The atmosphere-ocean coupling showed an overall positive impact on the simulation of sea surface heat fluxes and Mediterranean cyclones (medicanes and Vb-cyclones). Moreover, the representation of sea surface heat fluxes, wind speed, and medicane features was more realistic when using a finer atmospheric grid resolution (less than 10 km). The present study suggests that the combination of a finer atmospheric grid resolution together with atmosphere-ocean coupling is advantageous in simulating the Mediterranean climate system.
Floodplains and other wetlands depend on seasonal river flooding and play an important role in the terrestrial water cycle. They influence evapotranspiration, water storage and river discharge dynamics, and they are the habitat of a large number of animals and plants. Thus, to assess the Earth’s system and its changes, a robust understanding of the dynamics of floodplain wetlands including inundated areas, water storages, and water flows is required.
This PhD thesis aims at improving the modeling of large floodplains and wetlands within the global-scale hydrological model WaterGAP, in order to better estimate water flows and water storage variations in different storage compartments. Within the scope of this thesis, I have developed a new approach to simulate dynamic floodplain inundation on a global-scale. This approach introduces an algorithm into WaterGAP, which has a spatial resolution of 0.5 degree (longitude and latitude) globally. The new approach uses subgrid-scale topography, based on high-resolution digital elevation models, to describe the floodplain elevation profile within each grid cell by applying a hypsographic curve. The approach comprises the modeling of a two-way river-floodplain interaction, the separate downstream water transport within the river and the floodplain – both with temporally and spatially different variable flow velocities – and the floodplain-groundwater interactions. The WaterGAP version that includes the floodplain algorithm, WaterGAP 2.2b_fpl, estimates floodplain and river water storage, inundated area and water table elevation, and also simulates backwater effects.
WaterGAP 2.2b_fpl was applied to model river discharge, river flow velocity, water storages, water heights and surface water extent on a global-scale. Model results were comprehensively validated against ground observations and remote sensing data. Overall, the modeled and observed data are in agreement. In comparison to the former version WaterGAP 2.2b, the model performance has improved significantly. The improvements are most remarkable in the Amazon River basin. However, the seasonal variation of surface water extent and total water storage anomalies are still too low in many regions on the globe when compared to observations. A detailed analysis of the simulated results suggests that in the Amazon River basin the introduction of backwater effects is important for realistically simulating water storages and surface water extent. Future efforts should focus on the simulation of water levels in order to better model the flow routing according to water slope. To further improve the model performance in specific regions, I recommend that the globally constant model parameters that affect inundation initiation, river-floodplain interaction, DEM correction for vegetation, and backwater amount at basin or subbasin-scale be adjusted.