Refine
Year of publication
Document Type
- Doctoral Thesis (87) (remove)
Language
- English (87) (remove)
Has Fulltext
- yes (87)
Is part of the Bibliography
- no (87)
Keywords
- Boden (2)
- Caribbean (2)
- China (2)
- Roman (2)
- bank filtration (2)
- degradation (2)
- groundwater (2)
- 1,4-dioxane (1)
- 3D-Modellierung (1)
- Abbau (1)
Institute
- Geowissenschaften (87) (remove)
Paläoklimarekonstruktionen, die es sich zum Ziel gesetzt haben, Klima-Mensch Interaktionen auf lange Zeitreihen betrachtet zu erforschen, nehmen begünstigt durch die aktuell intensiv geführte Klimadebatte, einen immer größer werdenden Stellenwert in der öffentlichen und wissenschaftlichen Wahrnehmung ein. Denn trotz aller wissenschaftlicher Fortschritte, die in den vergangenen Jahrzehnten im Bereich der modernen Klimaforschung gemacht wurden, bleibt die zuverlässige Vorhersage und Modellierung von zukünftigen Klimaveränderungen noch immer eine der größten Herausforderungen unser heutigen Zeit. Betrachtet man die Karibik exemplarisch in diesem Rahmen, dann prognostizieren viele Modellrechnungen, infolge steigender Ozeantemperaturen, ein deutlich häufigeres Auftreten von tropischen Stürmen und Hurrikanen sowie eine Verschiebung hin zu höheren Sturmstärken. Dieser Trend stellt für die Karibik und viele daran angrenzende Staaten eine der größten Gefahren des modernen Klimawandels dar, den es wissenschaftlich über einen langen Zeitrahmen zu erforschen gilt.
Klimaprognosen stützen sich meist vollständig auf hoch-aufgelöste instrumentelle Datensätze. Diese sind aber alle durch einen wesentlichen Aspekt limitiert. Aufgrund ihrer eingeschränkten Verfügbarkeit (~150 Jahre) fehlt ihnen die erforderliche Tiefe, um die auf langen Zeitskalen operierenden Prozesse der globalen Klimadynamik adäquat abbilden zu können. Betrachtet man das Holozän in seiner Gesamtheit, so wurde die globale Klimadynamik über die vergangenen ~11,700 Jahre von periodisch auftretenden Prozessen und Abläufen gesteuert. Diese wirken grundsätzlich über Zeiträume von mehreren Jahrzehnten, teilweise Jahrhunderten und in einigen Fällen sogar Jahrtausenden. Viele dieser natürlichen Prozesse, können in der kurzen Instrumentellen Ära nicht gänzlich identifiziert und angemessen in Klimamodellen berücksichtig werden. Die alleinige Berücksichtigung der Instrumentellen Ära bietet daher nur eine eingeschränkte Perspektive, um die Ursachen und Abläufe von vergangenen sowie mögliche Folgen von zukünftigen Klimaveränderungen zu verstehen. Um diese Einschränkung zu überwinden, ist es somit erforderlich, dass die geowissenschaftliche Forschung mit Proxymethoden ein zusammenfassendes und mechanistisches Verständnis über alle Holozänen Klimaveränderungen erlangt.
Wenn man sich diese Limitierung, die ansteigenden Ozeantemperaturen und das in der Karibik in den vergangen 20 Jahren vermehrte Auftreten von starken tropischen Zyklonen ins Gedächtnis ruft, ist es nachvollziehbar, dass im Rahmen dieser Doktorarbeit ein zwei Jahrtausende langer und jährlich aufgelöster Klimadatensatz erarbeitet werden soll, der spät Holozäne Variationen von Ozeanoberflächenwasser-temperaturen (SST) und daraus resultierende lang-zeitliche Veränderungen in der Häufigkeit tropischer Zyklone widerspiegelt. In Zentralamerika wird das Ende der Maya Hochkultur (900-1100 n.Chr.) mit drastischen Umweltveränderungen (z.B. Dürren) assoziiert, die während der Mittelalterlichen Warmzeit (MWP; 900-1400 n.Chr.) durch eine globale Klimaveränderung hervorgerufen wurde. Die aus einem „Blue Hole“ abgeleiteten Informationen über Klimavariationen der Vergangenheit können als Referenz für die gegenwärtige Klimakriese verwendet werden.
Als „Blue Hole“ wird eine Karsthöhle bezeichnet, die sich subaerisch während vergangener Meeresspiegeltiefstände im karbonatischen Gerüst eines Riffsystems gebildet hat und in Folge eines Meeresspiegelanstiegs vollständig überflutet wurde. In einigen wenigen marinen „Blue Holes“ treten anoxische Bodenwasserbedingungen auf. Die in diesen anoxischen Karsthöhlen abgelagerten Abfolgen mariner Sedimente können als einzigartiges Klimaarchiv verwendet werden, da sie aufgrund des Fehlens von Bioturbation eine jährliche Schichtung (Warvierung) aufweisen.
In dieser kumulativen Dissertation über das „Great Blue Hole“ werden die Ergebnisse eines 3-jährigen Forschungsprojekts vorgestellt, dass das Ziel verfolgte einen wissenschaftlich herausragenden spät Holozänen Klimadatensatz für die süd-westliche Karibik zu erzeugen. Beim „Great Blue Hole“ handelt es sich um ein weltweit einzigartiges marines Sedimentarchiv für diverse spät Holozäne Klima-veränderungen, das im Zuge dieser Dissertation sowohl nach paläoklimatischen als auch nach sedimentologischen Fragestellungen untersucht wurde. Die vorliegende Doktorarbeit befasst sich im Einzelnen mit (1) der Ausarbeitung eines jährlich aufgelösten Archives für tropische Zyklone, (2) der Entwicklung eines jährlich aufgelösten SST Datensatzes und (3) einer kompositionellen Quantifizierung der sedimentären Abfolgen sowie einer faziell-stratigraphischen Charakterisierung von Schönwetter-Sedimenten und Sturmlagen. Zu jedem dieser drei Aspekte, wurde jeweils ein Fachartikel bei einer anerkannten wissenschaftlichen Fachzeitschrift mit „peer-review“ Verfahren veröffentlicht.
Der insgesamt 8.55 m lange Sedimentbohrkern („BH6“), der für diese Dissertation untersucht wurde, stammt vom Boden des 125 m tiefen und 320 m breiten „Great Blue Holes“, das sich in der flachen östlichen Lagune des 80 km vor der Küste von Belize (Zentralamerika) gelegenen „Lighthouse Reef“ Atolls befindet. Durch seine besondere Geomorphologie wirkt das, innerhalb des atlantischen „Hurrikan Gürtels“ positionierte, „Great Blue Hole“ wie eine gigantische Sedimentfalle. Die unter Schönwetter-Bedingungen kontinuierlich abgelagerten Abfolgen feinkörniger karbonatischer Sedimente, werden von groben Sturmlagen unterbrochen, die auf „over-wash“ Prozesse von tropischen Zyklonen zurückzuführen sind.
...
A graph theoretical approach to the analysis, comparison, and enumeration of crystal structures
(2008)
As an alternative approach to lattices and space groups, this work explores graph theory as a means to model crystal structures. The approach uses quotient graphs and nets - the graph theoretical equivalent of cells and lattices - to represent crystal structures. After a short review of related work, new classes of cycles in nets are introduced and their ability to distinguish between non-isomorphic nets and their computational complexity are evaluated. Then, two methods to estimate a structure’s density from the corresponding net are proposed. The first uses coordination sequences to estimate the number of nodes in a sphere, whereas the second method determines the maximal volume of a unit cell. Based on the quotient graph only, methods are proposed to determine whether nets consist of islands, chains, planes, or penetrating, disconnected sub-nets. An algorithm for the enumeration of crystal structures is revised and extended to a search for structures possessing certain properties. Particular attention is given to the exclusion of redundant nets and those, which, by the nature of their connectivity, cannot correspond to a crystal structure. Nets with four four-coordinated nodes, corresponding to sp3 hybridised carbon polymorphs with four atoms per unit cell, are completely enumerated in order to demonstrate the approach. In order to render quotient graphs and nets independent from crystal structures, they are reintroduced in a purely graph-theoretical way. Based on this, the issue of iso- and automorphism of nets is reexamined. It is shown that the topology of a net (that is the bonds in a crystal) constrains severely the symmetry of the embedding (that is the crystal), and in the case of connected nets the space group except for the setting. Several examples are studied and conclusions on phases are drawn (pseudo-cubic FeS2 versus pyrite; α- versus β- quartz; marcasite- versus rutile-like phases). As the automorphisms of certain quotient graphs stipulate a translational symmetry higher than an arbitrary embedding of the corresponding net would show, they are examined in more detail and a method to reduce the size of such quotient graphs is proposed. Besides two instructional examples with 2-dimensional graphs, the halite, calcite, magnesite, barytocalcite, and a strontium feldspar structures are discussed. For some of the structures it is shown that the quotient graph which is equivalent to a centred cell is reduced to a quotient graph equivalent to the primitive cell. For the partially disordered strontium feldspar, it is shown that even if it could be annealed to an ordered structure, the unit cell would likely remain unchanged. For the calcite and barytocalcite structures it is shown that the equivalent nets are not isomorphic.
The Earth’s surface condition we find today is a result of long exposure to metabolism of life forms. Particularly, molecular oxygen in the atmosphere is a feature which developed over time. The first substantial and lasting rise of atmospheric oxygen level happened ≈ 2.5 Ga ago, but localities are reported where transiently elevated oxygen levels appeared before this time-point. To trace the timing and circumstances of the earliest availability of free oxygen in the atmosphere is important to understand the habitats of early microbial life forms on Earth.
This thesis focuses to obtain information of oxygen levels and the related atmospheric cycling of metals in sediments of the 3.5 to 3.2 Ga Barberton Greenstone Belt. First, as iron was a ubiquitous constituent of Archean seawater, I investigated its isotopic composition in minerals of chemical sediments. Hereby, I tried to resolve the changes within the water basin on small scale sedimentary sequence cycles. Second, I focused on the minor constituents of Archean seawater. The Re-Os geochronologic system and the abundance patterns of the platinum-group elements were chosen to integrate information of oxygen promoted weathering of a large source area. To integrate information of a large time interval, the isotopes of uranium were investigated over a large stratigraphic section.
The two key findings of this thesis are:
• Quantitative oxidation of ferrous iron in surface layers of Paleoarchean seawater occurred during the onset and termination of hydrothermal FeIIaq delivery into shallow waters.
• Paleoarchean sedimentary successions of the Barberton Greenstone Belt lack any evidence of transient basin-scale oxygenation.
The Manzimnyama Iron Formation (IF, Fig Tree Group, Barberton Greenstone Belt, South Africa) has been deciphered to exist of cyclic stacks of lithostratigraphic units with varying amounts of iron oxide and carbonate minerals. In-situ femtosecond-Laser-Ablation ICP-MS iron isotope measurements showed that the majority of siderite (γ56Fe ≈ −0.5 ‰) precipitated directly from seawater of γ56Fe ≈ 0 ‰. Ferric iron from the surface layers is preserved in ≤ 1μ m hematite and in magnetite that has been grown within the consolidated sediment. During FeIIaq events, fine-grained hematite (γ56Fe ≈ 2.2 ‰) and magnetite (γ56Fe 0.5 to 0.8 ‰) indicate oxygen levels in surface waters of lower than 0.0002 μM. Upon onset and termination of iron oxide abundance, magnetite with γ56Fe ≈ 0 ‰ indicates that low concentrations of FeIIaq in surface waters were oxidized quantitatively. These observations demonstrate the existence of iron oxidation in Paleoarchean surface waters independent of FeIIaq concentration. This is the first investigation of Paleoarchean IF showing that lithostratigraphic cyclicity can be traced in iron isotopic composition of oxide minerals.
ID-ICP-MS measurement of Re, Ir, Ru, Pt and Pd, trace element (SF-ICP-MS) and ID-MCICP- MS uranium isotope determination have been applied to carbonaceous shale of the Mapepe Fm. (Fig Tree Group) after inverse Aqua Regia leaching and bulk digestion. The sediments reveal a silicified fraction which exhibits a seawater REE signature and a mixture of detrital and meteoritic PGE. Neither enrichment of the redox-sensitive elements Re or Mo nor fractionated uranium isotopes have been found on a stratigraphic interval of several hundred meters. The non-silica fraction shows no depletion of Re which indicates that the detrital material had no contact to oxidizing fluids. ID-TIMS measurements of Re and Os after the CrO3-SO4 Carius Tube method of two sample intervals showed that the Re-Os isotopic systems of the non-silica fractions are identical to two komatiite occurrences. Weltevreden Fm. and Komati Fm. rocks were uplifted, eroded and transported to the deep part of the sedimentary basin without any change to the Re-Os system. Negative fractionated uranium isotopes (γ238U = −0.41 ± 0.01 ‰) associated with detrital Ba-Cr-U occurrences suggest the existence of distal redox-processes that involve uranium species. This study demonstrates that over the time of exposure and deposition of the Mapepe Fm. sedimentation, free oxygen was not available for weathering in the catchment area.
The African continent is regularly portrayed as an indolent space with a well-known reputation as a chaotic continent. Viewed as lacking vision, means and capacities, Africa is perceived at best as a place that is marked by a permanent status quo, stagnation, or in worst case scenarios, as a declining continent. Various references to the continent are synonymous with famine, poverty, war, etc. Such portrayals are all the more intriguing given that the continent is known for its abundant natural resources, such as timber, oil, natural gas, minerals, etc., whose reserves are, moreover, not well known both by the African people and their leaders. As a result, there is still much progress to be made in tapping into the resources in order to improve the daily lives of African citizens.
In such a context dominated by infantile carelessness throughout the continent, the interventions of actors from outside the continent are the only hopes of bringing some vitality to this continent which is cloaked in "la grande nuit – the great darkness" (Mbembé 2013). Thus during the main sequences of recent history, representing different forms of Western penetration and activity on the African continent (slavery, imperialism, colonization), all the Western world’s contributions have obviously not sufficed to boost Africa and take it out of its never ending childhood. It has remained just as passive and apathetic today as it was yesterday.
The attraction of Asian actors to the continent is even more recent. And consistent with its abovementioned indolence, Africa is seen as an easy and defenceless prey for the Korean, Japanese, Indian, Malaysian, or Chinese conquerors. In the latter case, the insatiable appetite for natural resources whose reserves are being rapidly depleted is the cornerstone of their foreign aid policy. This led China to colonize the continent, showing a preference for Pariah Regimes which held no appeal for the West, by sending an army of workers to extract those resources (Lum et al. 2009), in defiance of all national and international regulations and based on completely opaque contracts.
Although the concept of African Agency was rapidly developed in several African countries, the aim of this study was more specific to Cameroon’s mining sector in which different entrepreneurs from abroad got involved over time. The thesis investigates whether indigenous citizens took part in any way in the development of mining projects in the country. Thus, the work assesses and analyses actions and reactions initiated and undertaken by local people in the context of China’s presence within Cameroon’s mining sector to promote and advance their interests over those of foreign investors. In addition, the author has no knowledge of any other study investigating African Agency in the mining sector as a whole in Cameroon.
In conducting this study, a multi-method research framework was developed including a series of methods used to collect data and analyse concepts of African Agency associated Political Ecology as they developed within Cameroon’s mining sector. Specifically, those methods comprised quantitative research when it came to collecting data using a positivist and empirical approach constructed by deducing evidence from statistical data collected by means of the 167 questionnaire surveys administered to local inhabitants and workers randomly selected on mining sites and in riparian communities. The questionnaires helped to capture Cameroonians' perceptions of the recent phenomenon of the gradual but significant influx of international actors and precisely Chinese players in the mining sector on the one hand, and on the other hand, observational data was collected across the GVC as developed in the Betare-Oya region. As a complement to the former technique, qualitative methods helped to study and deepen understanding of human behaviour and the social world in a holistic perspective through individual interviews, focus groups, and direct observations on the ground. In addition, the spatial analysis method based on the land use classification technique served to detect changes to land use/land cover that have been brought on by mechanised mining activities undertaken in this region. The sequencing of data collected and their processing from a ground theory perspective led to the formulation and specification of Cameroon’s Ecological Agency theory.
One of the earliest steps of this work consisted in a literature review and in placing the African Agency concept in a broader context. It then led to the state of the art, specifications about research content of the work and the main theories undergirding this thesis. Before examining developments that emerged during the last decade, a historical perspective was provided to the topic in order to show how African societies started mining operations and how they dealt with foreign partners interested in their mining resources. The aim was to show that while Western imperialism presented a challenge for the sector, it did not erase local participation, even despite the constraints associated with such involvement.
...
This work describes the development and characterization of two instruments and their data evaluation, which contributes to a better understanding of new particle formation and growth, as well as their interactions with clouds. Both instruments were characterized at the Cosmics Leaving Outdoor Droplets (CLOUD) experiment at the European Center for Nuclear Research (CERN).
The crude oil constituents benzene, toluene, ethylbenzene, and the three xylene isomers (BTEX) are the dominating groundwater contaminants originating from surface spill accidents by oil production facilities and with gasoline and jet fuel. Thereby BTEX posing a threat to the world´s scarce drinking water resources due to their water solubility and toxicity. An active remediation cleanup involving a BTEX event proves not only to be very expensive but almost impossible when it comes to the complete removal of contaminants from the subsurface. A favoured and common practice is combining an active remediation process focussing on the source of contamination coupled together with the monitoring of the residual contamination in the subsurface (monitored natural attenuation; MNA). MNA include all naturally occuring biological, chemical and physical processes in the subsurface. The general goal of this work was to improve the knowledge of biodegradation of aromatic hydrocarbons under anaerobic conditions in groundwater. For this groundwater and soil at the former military underground storage tank (UST) site Schäferhof – Süd near Nienburg/Weser (Niedersachsen, Germany) were sampled and analysed. The investigations were done in collaboration of the Umweltbundesamt, the universitys of Frankfurt and Bremen and the alphacon GmbH Ganderkesee. To investigate the extent of groundwater contamination, the terminal electron acceptor processes (TEAPs) and the metabolites of BTEX degradation in groundwater, six observation wells were sampled at regular intervals between January 2002 and September 2004. The wells were positioned in order to cover the upstream, the source area and the downstream of the presumed contamination source. Additionally, vertical sediment profiles were sampled and investigated with respect to spreading and concentration of BTEX in the subsurface. A large residual contamination involving BTEX is present in soil and groundwater at the studied locality. Maximum BTEX concentration values of 17 mg/kg were recorded in analysing sediment in the unsaturated zone. In the capillary fringe, values of 450 mg/kg were recorded (October 2004) and in the saturated zone maximum values of 6.7 mg/kg BTEX were detected. The groundwater samples indicate increasing BTEX concentrations in the groundwater flow direction (from 532 µg/l up to 3300 µg/l (mean values)). Biodegradation of aromatic hydrocarbons under anaerobic conditions in the sub surface at contaminated sites is characterised by generation of metabolites. From the monoaromatic hydrocarbons BTEX metabolites such as benzoic acid (BA) and the methylated homologs and C1-and C2-benzyl-succinic acids (BSA) are generated as intermediates. A solid-phase extraction method based on octadecyl-bonded silica sorbent has been developed to concentrate such metabolite compounds from water samples followed by derivatization and gas chromatography/mass spectrometry (GC/MS) of the extracts. The recovery rate range between 75 and 97%. The method detection limit was 0.8 µg/l. Organic acids were identified as metabolic by-products of biodegradation. Benzoic acid, C1-, C2- and C3-benzoic acid were determined in all contaminated wells with considerable concentrations. Furthermore, the depletion of the dominant terminal electron acceptors (TEAs) oxygen, nitrate, and sulphate and the production of dissolved ferrous iron and methane in groundwater indicate biological mediated processes in the plume evidently proving the occurrence of NA. A large overlap of different redox zones at the studied part of the plume has been observed. A important finding in this study is the strong influence of groundwater level fluctuations on the BTEX concentration in groundwater. A very dry summer in 2003 was recorded during the monitoring period, resulting on site in a drop of the groundwater level to 1.7 m and a concomitant increase of BTEX concentrations from 240 µg/l to 1300 µg/l. The groundwater level fluctuations, natural degradation and retention processes essentially influence BTEX concentrations in the groundwater. Groundwater level fluctuations have by far a stronger influence than the influence of biological degradation. Increasing BTEX concentrations are hence not a consequence of limited biological degradation. Another part of the study was to observe the isotopic fractionation of the electron acceptor Fe(III), due to biologically mediated reduction of Fe(III) to the watersoluble Fe(II) at the site and first field data are presented. Both groundwater and sediment samples were analysed with respect to their Fe isotopic compositions using high mass resolution Multi Collector-Inductively Coupled Plasma-Mass Spectrometry (MC-ICP-MS). The delta56Fe -values of groundwater samples taken from observation wells located downstream of the source area were isotopically lighter than delta56Fe -values obtained from groundwater in the uncontaminated well. The Fe isotopic composition of most parts of the sediment profile was similar to the Fe isotopic composition of uncontaminated groundwater. Thus, a significant iron isotope fractionation can be observed between sediment and groundwater downstream of the BTEX contamination.
Ziel der Arbeit war es, die Flugzeitmassenspektrometrie als neue Analysemethode für die instrumentelle Analytik halogenierter Spurengase in der Luft zu etablieren. Die grundle-gende Motivation dafür ist, dass anthropogene Emissionen vieler Vertreter dieser Sub-stanzklasse einen negativen Einfluss auf die Umwelt zeigen: in der Atmosphäre agieren die Substanzen bzw. ihre Abbauprodukte als Katalysatoren für den stratosphärischen Ozonab-bau und verstärken den Strahlungsantrieb der Erde durch Absorption elektromagnetischer Strahlung im sogenannten atmosphärischen Fenster. Um diese Effekte und deren Auswir-kung quantifizieren zu können, ist es notwendig, Konzentrationen und Trends der Substan-zen in der Atmosphäre zu überwachen. Nur so können Gegenmaßnahmen wie Produktions-reglementierungen geplant und bewertet werden. In Kombination mit inverser Modellie-rung können zudem Rückschlüsse über tatsächlich emittierten Mengen gezogen werden. Dies stellt den Anspruch an die Analytik, sehr geringe Mengen dieser Gase sehr präzise quantifizieren zu können, um auch schwache Trends zu erkennen. Zudem muss die Analy-semethode die Möglichkeit zu bieten, mit der wachsenden Anzahl bekannter und zu über-wachender Substanzen Schritt zu halten. Besonders für letzteren Aspekt bietet die Flug-zeitmassenspektrometrie einen entscheidenden Vorteil gegenüber der „konventionellen“ Methode, der Quadrupolmassenspektrometrie: sie zeichnet das gesamten Massenspektrum auf ohne dadurch an Empfindlichkeit einzubüßen. Um das atmosphärische Mischungsver-hältnis von Substanzen im Bereich von pmol mol−1 bis fmol mol−1 bestimmen zu können, muss das Quadrupolmassenspektrometer im Single Ion Monitoring Modus betrieben wer-den – so wird zwar eine hohe Sensitivität erreicht, es wird aber auch nur die Intensität eines bestimmten Masse zu Ladungsverhältnisses (kurz: Masse) zu einem Zeitpunkt aufgezeich-net. Ein Flugzeitmassenspektrometer hingegen extrahiert Ionen mit einer Frequenz im Ki-loherzbereich und zeichnet für jede Extraktion das vollständige Flugzeitspektrum und da-mit Massenspektrum auf.
Aufgabe dieser Arbeit war es, ein Flugzeitmassenspektrometer mit vorgeschalteter Pro-benanreicherungseinheit sowie Gaschromatograph zur Trennung des Subtanzgemisches vor der Detektion aufzubauen und Werkzeuge zur Datenauswertung zu entwickeln. Um einen zukünftigen Feldeinsatz vorzubereiten, sollte der Aufbau möglichst kompakt, mobil und vollständig automatisiert sein. Anschließend sollte Empfindlichkeit, Präzision und dynami-scher Messbereich geprüft, optimiert und die Anwendbarkeit zur Analyse halogenierter Spurengase gezeigt werden. Die Ergebnisse aus der in der vorliegenden Arbeit präsentier-ten Geräteentwicklung finden sich in drei Publikationen wieder, welche in thematischer Reihenfolge die Probenanreicherung (Obersteiner et al., 2016b), den Vergleich von Quadrupol- und Flugzeitmassenspektrometrie (Hoker et al., 2015) sowie Eigenschaften und Anwendung des neuen Aufbaus (Obersteiner et al., 2016a) behandeln. Mit den genannten Aufsätzen ist die Arbeitsgruppe Engel weltweit die erste, welche hochpräzise Analytik ha-logenierter Spurengase routinemäßig mittels Flugzeitmassenspektrometrie durchführt. Der nächste Schritt ist der Übergang von der Laboranwendung zur Feldmessung, z.B. in Form von bodenbasierter in situ Analyse troposphärischer Luftmassen am Taunus Observatorium auf dem Kleinen Feldberg. Da es bisher keine Messstation für die hier beschriebene analy-tische Fragestellung in Deutschland gibt, könnte eine deutliche Verbesserung der Überwa-chung halogenierter Treibhausgase und ozonzerstörender Substanzen in Europa erzielt wer-den. Weiterhin wäre eine Flugzeugapplikation in Zukunft denkbar, welche neben der durch das Flugzeitmassenspektrometer abgedeckten Substanzbandbreite auch von dessen hoher möglicher Spektrenrate profitieren könnte. In Kombination mit Hochgeschwindigkeitsgas-chromatographie könnte eine bisher unerreichte Zeitauflösung der Beprobung der Atmo-sphäre mittels Gaschromatographie-Massenspektrometrie erzielt werden.
The East African Rift System (EARS) was initiated in the Eocene epoch between 50 and 21 Ma probably due to the influence of mantle plumes that caused volcanism, flood basalts and rifting extensions in Ethiopa and the Afar region. As a result of magmatic intrusions and adiabatic decompression melting within the lithosphere caused by the impact of the Kenya plume, there was a southward propagation of the EARS of about 30 – 15 Ma from Ethiopia to Kenya, which coincide with the occurrence of volcanism. The EARS developed towards the south along the margins of the Tanzania Craton between 15 and 8 Ma. Previous findings of low-velocity anomalies within the upper mantle and the mantle transition zone indicate an upwelling of hot mantle material in the vicinity of the Afar region and the East African Rift. This study includes the analysis of P- and S-receiver functions in order to determine further impacts on the lithosphere from below. The aim was to determine the topographic undulations of further boundary layers and to identify their variability owing to the rifting processes and the formation of the EARS. The study area included the Tanzania Craton and the surrounding rift branches of the East African Rift System.
The region of the Rwenzori Mountains can be analysed in detail because of the large dataset of the RiftLink project. The use of the P-receiver function technique and the H-K stacking method enabled to determine different vP /vS ratios depending on the tectonic setting in the Rwenzori region: Rift shoulders (vP /vS =1.74), Albert Rift segment (vP /vS =1.80), Edward Rift segment (vP /vS =1.87) and Rwenzori Mountains (vP /vS =1.86). To determine the topography of the Moho, it is necessary to take into account the thickness of the sedimentary layer, the surface topography, the azimuthal variations in crustal thickness and the impact of local anomalies. After correcting these effects on the Moho depths, significant variations in Moho topography could be determined. The Moho depths range from 29 to 39 km beneath the rift shoulders of the Albertine Rift. Within the rift valley, the crustal thickness varies between 25 – 31 km in the Edward Rift segment and 22 – 30 km in the Albert Rift segment. An averaged crustal thickness of about 26 km within the rift valley indicates the lack of the crustal root beneath the Rwenzoris. Similar variations in crustal thickness were determined by using an automatic procedure for analysing S-receiver functions that was developed in this study.
The S-receiver functions are created by applying a rotation criterion in order to rotate the Z, N and E components into the L, Q and T components. It is necessary to perform trial rotations using different incident and azimuth angles to determine the correct rotation angles. The latter are identified by the use of the rotation criterion, including the amplitude ratio of the converted Moho signal to the direct S/SKS-wave signal. The L component is rotated correctly in the direction of the incident shear wave in the case of the maximum amplitude ratio. After analysing the frequency content of the receiver functions in order to sort out harmonic and long-periodic traces, the individual Moho signals are checked for consistency in order to remove atypic signals. To increase the signal-to-noise ratios on the traces, the S-receiver functions are stacked. For this purpose, the signals of the direct shear waves must originate from similar epicenters. On the basis of similar ray paths, the receiver functions show comparable waveforms and converted signals. To perform the stacking procedure, it is necessary to merge the datasets of the adjacent stations in order to obtain a sufficient number of receiver functions. This analysis is based on the assumption that the incident seismic waves arriving at the adjacent stations penetrate to some extent the same underground structures in the case of similar wave propagation paths. This approach accounts for the fact that the converted signals do not result exclusively from the piercing points at the boundary layers. Further signals originate from the conversions at the boundary layer within the Fresnel Zone. The piercing points are derived from the significant signals in the receiver functions. Depending on the order of arrival of the converted phases on the traces, the signals are attributed to the theoretical discontinuities DIS1, DIS2, DIS3 and DIS4. However, partly due to the low signal-to-noise ratios on the traces, it is difficult to identify the real conversions on the traces and to ensure that the converted signals are attributed to the correct boundary layers. For this reason, it is necessary to check the consistency of the conversion depths among each other. In the case of inconsistent conversion depths, the corresponding signals are either adjusted to another seismic boundary layer or removed from the dataset. To verify the functionality of the automatic procedure and to determine the resolvability with respect to two boundary layers, several models are tested including horizontal and dipping discontinuities. To resolve distinct discontinuities, their depths must differ by at least 60 km, otherwise, due to similar depth ranges of the different boundary layers, the converted signals cannot be separated from each other. As a consequence, the converted signals that originate from different discontinuities are attributed to a single one. Further tests including break-off edges of seismic discontinuities are performed to check the attributions of the converted signals to the discontinuities. Owing to the varying number of boundary layers, the converted signals cannot be attributed to the discontinuities according to the order of their arrivals on the traces. It is necessary to correct their attributions to the seismic discontinuities in order to resolve the boundary layers.
The crust-mantle boundary and further discontinuities within the lithospheric mantle are investigated by applying this automatic procedure. Depending on the tectonic setting, the conversion depths of the Moho range from about 30 – 45 km beneath the western rift shoulder to 20 – 35 km within the rift valley up to 30 – 40 km beneath the eastern rift shoulder. The long wavelengths of the shear waves hamper the correct identification of the converted phases in the S-receiver functions. With respect to the relative differences in conversion depth, the topographic undulations of the crust-mantle boundary are consistent with the Moho depths derived from P-receiver functions. In contrast to the Rwenzori region, it is difficult to resolve completely the trend of the Moho in the remaining area of the East African Rift due to the small dataset provided by IRIS. The results exibit an increase in crustal thickness to up to 45 km in the region of the Cenozoic volcanics such as Virunga, Kivu, Rungwe and Kenya. The greatest Moho depths of more than 50 km are located near Mount Kilimanjaro. In addition to the Moho, the analysis of the S-receiver functions revealed two further boundary layers at depths of 60 – 140 km and 110 – 260 km, which are associated with a mid-lithospheric discontinuity and the lithosphere-asthenosphere boundary, respectively. The shallowest conversion depths of the LAB are focussed to small-scale regions within the rift branches, namely the northern Albertine Rift, the Chyulu Hills and the Mozambique Belt, which are located around the Tanzania Craton. The larger thickness of the lithosphere beneath the cratonic terrain indicates that the Tanzania Craton is not significantly eroded. However, there are indications that the lithosphere beneath the craton and the rift branches is penetrated by ascending asthenospheric melts to depths of up to 140 and 60 km, respectively. The top of the ascending melts is associated with the occurrence of the mid-lithospheric discontinuity. The shallowest conversion depths of this boundary layer (60 – 90 km) are related to the rifted areas of the EARS and the Cenozoic volcanic provinces, which are located along the Albertine Rift, the Kenya Rift and the Rukwa-Malawi rift zones. The deepest conversion depths of up to 140 km are related to the Rwenzori Belt, the Ugandan Basement Complex and the interior of the Tanzania Craton.
In the past sixty years, excessive water consumption and dam construction have significantly influenced natural flow regimes and surface freshwater ecosystems throughout China, and thus resulted in serious environmental problems. In order to balance the competing water demands between human and environment and provide knowledge on sustainable water management, assessments on anthropogenic flow alterations and their impacts on aquatic and riparian ecosystems in China are needed.
In this study, the first evaluation on quantitative relationships between anthropogenic flow alterations and ecological responses in eleven river basins and watersheds in China was performed based on the data that could be obtained from published case studies. Quantitative relationships between changes in average annual discharge, seasonal low flow and seasonal high flow and changes in ecological indicators (fish diversity, fish catch and vegetation cover, etc.) were analyzed. The results showed that changes in riparian vegetation cover as well as changes in fish diversity and fish catch were strongly correlated with the changes in flow magnitude (r = 0.77, 0.66), especially with changes in average annual river discharge. In addition, more than half of the variations in vegetation cover could be explained by changes in average annual river discharge (r² = 0.63) and roughly 50 % changes in fish catch in arid and semi-arid region and 60% changes of fish catch in humid region could be related to alterations in average annual river discharge (r² = 0.53, 0.58).
In a supplementary analysis of this study, the first estimation on quantitative relationships between decreases in native fish species richness and anthropogenic flow alterations in 34 river basins and sub-basins in China was conducted. Linear relationships between losses of native fish species and five ecologically relevant flow indicators were analyzed by single and multiple regression models. For the single regression analysis, significant linear relationships were detected for the indicators of long-term average annual discharge (ILTA) and statistical low flow Q90 (IQ90). For the multiple regressions, no indicator other than ILTA has significant relationships with changes in number of fish species mainly due to collinearity. Two conclusions emerged from the analysis: 1) losses of fish species were positively correlated with changes in ILTA in China and 2) indicator of ILTA was dominant over other flow indicators included in this research for the given dataset. These results provide a guideline for the sustainable water resources management in rivers with high risk of fish extinction in China.
Atmospheric nanoaerosols have extensive effects on the Earth’s climate and human health. This cumulative work focuses on the development and characterization of instrumentation for measuring various parameters of atmospheric nanoaerosols, and its use to understand new particle formation from organic precursors. The principal research question is, how the chemical composition of nanoaerosol particles can be measured and how atmospheric chemistry influences aerosol processes, especially new particle formation and growth. Therefore, nanoaerosols are investigated under various aspects. More specifically, an instrument is developed to analyze nanoparticles, and field as well as chamber studies are conducted.
The main project is the instrument development of the Thermal Desorption Differential Mobility Analyzer (TD-DMA, project 1, Wagner et al. (2018)). This instrument analyzes the chemical composition of small aerosol particles. By characterization and testing in chamber experiments, it is proven to be suitable for the analysis of freshly nucleated particles.
The second project (Wagner et al. (2017)) applies a broad spectrum of aerosol measurement instruments for the characterization of aerosol particles produced by a skyscraper blasting. A comprehensive picture of the particle population emitted by the demolition is obtained.
Project 3 (K¨urten et al. (2016)) is also an ambient aerosol measurement, focusing of new particle formation in a rural area in central Germany, and the ability of a negative nitrate CI-APi-TOF to detect various substances in atmosphere. Project 4 (Heinritzi et al. (2016)) is a characterization of the negative nitrate CI-APi-TOF used in projects 1, 3, 5, 6, 7 and 8. The following projects focus on understanding new particle formation from atmospherically abundant organic precursors. Key instruments comprise the negative nitrate CI-APiTOF for gas-phase measurements of the nucleating species, and various sizing and counting instruments for quantifying the particle formation and growth. Project 5 (Kirkby et al. (2016)) shows that biogenic organic compounds formed from alpha-pinene can nucleate on their own without the influence of e.g. sulfuric acid. Project 6 (Tr¨ostl et al. (2016)) describe the subsequent growth of these particles. Project 7 (Stolzenburg et al. (2018)) covers the temperature dependence of this growth and in project 8 (Heinritzi et al. (2018)), the suppressing influence of isoprene on the new particle formation is assessed.
Groundwater is the largest source of accessible freshwater with its dynamics having significantly changed due to human withdrawals, and being projected to continue to as a result of climate change. The pumping of groundwater has led to lowered water tables, decreased base flow, and depletion.
Global hydrological models (GHMs) are used to simulate the global freshwater cycle, assessing impacts of changes in climate and human freshwater use. Currently, groundwater is commonly represented by a bucket-like linear storage component in these models. Bucket models, however, cannot provide information on the location of the groundwater table. Due to this limitation, they can only simulate groundwater discharge to surface water bodies but not recharge from surface water to groundwater and calculate no lateral and vertical groundwater flow whatsoever among grid cells. For instance this may lead to an underestimation of groundwater resources in semiarid areas, where groundwater is often replenished by surface water. In order to overcome these limitations it is necessary to replace the linear groundwater model in GHMs with a hydraulic head gradient-based groundwater flow model
This thesis presents the newly developed global groundwater model G3M and its coupling to the GHM WaterGAP spanning over 70,000 lines of newly developed code. Development and validation of the modeling software are discussed along with numerical challenges. Based on the newly developed software, a global natural equilibrium groundwater model is presented showing better agreements with observations than previous models. Groundwater discharge to rivers is found to be the most dominant flow component globally, compared to flows to other surface water bodies and lateral flows. Furthermore, first global maps of the distribution of gaining and losing surface water bodies are displayed.
For the purpose of determining the uncertainty in model outcomes a sensitivity study is conducted with an innovative approach through applying a global sensitivity analysis for a computationally complex model. First global maps of spatially distributed parameter sensitivities are presented. The results at hand indicate that globally simulated hydraulic heads are equally sensitive to hydraulic conductivity, groundwater recharge and surface water body elevation, even though parameter sensitivities do vary regionally.
A high resolution model of New Zealand is developed to further understand the involved uncertainties connected to the spatial resolution of the global model. This thesis finds that a new understanding is necessary how these models can be evaluated and that a simple increase in spatial resolution is not improving the model performance when compared to observations.
Alongside the assessment of the natural equilibrium, the concept of a fully coupled transient model as integrated storage component replacing the former model in the hydrological model WaterGAP is discussed. First results reveal that the model shows reasonable response to seasonal variability although it contains persistent head trends leading to global overestimates of water table depth due to an incomplete coupling. Nonetheless, WaterGAP-G3M is already able to show plausible long term storage trends for areas that are known to be affected by groundwater depletion. In comparison with two established regional models in the Central Valley the coupled model shows a highly promising simulation of storage declines.
The biomarker record in two different lakes in central Europe, Lake Albano and Lake Constance, is used to reflect environmental changes and lake system response during the Late Glacial and Holocene. Extractable organic compounds in lake sediments, which can be assigned to their biological source (biomarkers) function as fingerprints of past aquatic or land plant organisms. Using gas chromatography coupled with mass spectrometry, 21 different biomarkers (predominantly steroids and triterpenoids) as well as a variety of n-alkanes, nalkanols, and n-alkanoic acids could be identified in the sediment records of Lake Albano and Lake Constance. In the Holocene sediments of Lake Albano, the distribution of biomarkers such as dinosterol (dinoflagellates), isoarborinol, and diplopterol (aquatic organisms) indicate three biomarker zones: The period between 0-3,800 years BP (zone 3) is characterized by high concentrations of these biomarkers and others such as tetrahymanol and diploptene. Conversely, zone 2 (3,800-6,500 years BP) shows very low concentrations of all autochthonous biomarkers. In zone 1 (6,500–11,480 years BP), dinosterol, isoarborinol, and diplopterol range on a relatively high level, whereas diploptene and tetrahymanol display comparatively low concentrations. The results suggest at least two distinct changes in the predominance of primary producers during the Holocene, which are related to changes in the lake system such as lake mixing and water column stratification. This interpretation is consistent with previous investigations of Lake Albano sediments including pigment and hydrogen index data (Ariztegui et al., 1996b; Guilizzoni et al., 2002). Allochthonous biomarkers such as long-chain n-alkanes, amyrenones and friedelin indicate a development from forest to a more open landscape from 6,000 and 5.000 years BP, respectively. After a period of high concentrations during the first half of the Holocene, all biomarkers derived from deciduous trees exhibit relatively low values until around 1,000 years BP. Again, this is consistent with results from previous pollen investigations (Ariztegui et al., 2000). The sediment core from Upper Lake Constance comprises the Late Glacial and Holocene. It was analysed for biomarkers and inorganic tracers in order to compare the biomarker results with other proxy data from the same core. Magnetic susceptibility (MS) was measured to get a high-resolution stratigraphic framework of the core and to obtain further information about changes of the proportions of allochthonous and autochthonous input. Enhanced concentrations and accumulation rates of dinosterol (biomarker for dinoflagellates) and biogenic calcite give evidence of increasing lake productivity at the beginning of the Holocene followed by a decrease in bioproductivity after around 7,000 years BP. Younger Dryas sediments are characterized by low amounts of both dinosterol and biogenic calcite indicating a low productivity. The comparison of the concentrations and accumulation rates of b-sitosterol and stigmastanol with parameters reflecting lake productivity suggests that both steroids in Lake Constance sediments are mainly derived from terrigenous sources. Biomarkers as well as concentrations and accumulation rates of allochthonous inorganic compounds such as titanium, magnesium and strontium indicate a slightly enhanced allochthonous input after 8,500 years BP. Significant increase of erosive matter input from enhanced soil erosion is not observed before 4,000 years BP. This can be attributed to the combined effects of precipitation increase as a result of climatic deterioration and anthropogenic deforestation which is consistent with observations from other lakes in Central Europe. The MS record of Lake Constance confirms these results by tracing the climatically induced shifts of more intense bioproduction (low MS caused by increased calcite deposition) during the ‘climatic optimum’. This is followed by increasing input of terrigenous sediment compounds during colder and wetter periods which lead to higher MS values in the lake sediments. The occurrence of tetrahymanol in Lake Constance sediments questions the unambiguous use of tetrahymanol as an indicator for water column stratification. Anaerobic organic macroaggregates within the oxygenated, photic zone of the water column have to be considered as a possible living space for anaerobic microorganisms containing tetrahymanol. The direct comparison of two very different lakes Albano and Constance with respect to biomarkers indicating climate or environmental change provides a contribution to the recent biomarker research for a better understanding of biomarkers in lacustrine sediments.
In this thesis, laboratory investigations have been conducted to investigate several processes occurring during the melt segregation (crystal settling and compaction processes), as well as during emplacement of plutons. With the help of three different sets of centrifuge experiments rates of these three magmatic processes have been evaluated. In the first series of the centrifuge experiments, the diapiric ascent of buoyant material from two source layers at different depths was studied. Through five models, the hypothesis of ascending diapirs was tested and it was demonstrated whether a rising diapir ascends straight upward or if its ascent might be deviated by another buoyant, softer – and consequently easier to travel through – layer which is located within the overburden strata. We were interested under which conditions they can be formed. For this purpose we placed perturbations on top of both the buoyant layers; either with a set-off of both the protrusions (for three of these experiments), or with both protrusion sitting directly on top of each other (for one of the experiments). In the first experiment, we omitted the perturbations, to test which pathways diapirs take which grow from natural Rayleigh-Taylor instabilities. Three others experiments differed in the viscosity contrast between the overburden and the buoyant material. Through the experimental runs, the effects of different overburden viscosities and perturbation positions on the number of the diapirs were observed. The modeling results show that two diapirs rising from the offset perturbations do not take the same pathway through the overburden layer. Rather, each diapir takes a different pathway, with the deeper diapir piercing through its overburden while rising, regardless if it was a buoyant layer or denser overburden layers. However, when the two perturbations were situated directly above each other in the different PDMS layers, this resulted in the formation of one big diapir rather than several smaller ones, and the overburden layer was less deformed than with offset perturbations. Diapiric structures as those derived from the models without perturbation and where the perturbation are offset occur within Great Kavir Basin (Iran), where numerous salt diapirs grew from several salt horizons, which show a similar spatial distribution. The resulting structure observed in the model where the two perturbations situated directly above each other, is close to what is observed in composite batholiths such as the Flasergranitoid Zone within the Bergsträßer Odenwald Crystalline Complex (Germany). The second series of models were aimed to study crystal settling within a magma. For this purpose experiments with an artificial magma of 30 vol% olivine in 70 vol% basaltic melt were conducted to elucidate the formation mechanisms and time scales of gravitational cumulates. Through the experiments, two physical processes have been observed: (i) purely mechanical compaction, and (ii) chemical compaction induced by dissolution and re-precipitation of settled crystals. The results reveals that the mechanical settling of the dense olivine suspension occurs at about 1/6 the speed of simple Stokes settling, and a sedimentation exponent n of 4.1 is found. Evidences of chemical compaction induced by dissolution and re-precipitation of settled crystals have been highlighted by a detailed analysis of the fine structure of olivine grain boundaries. This last has revealed (1) the presence of Ca, which is characteristic only for MORB-melt, at the interface of two adjacent Ol-grains even when no melt is present; (2) a not fully crystallized boundary layer between two adjacent olivine grains. The crystal size distribution curves and the grain size growth exponent n ~3.6 indicate that diffusion controlled Ostwald ripening is the dominant crystal growth mechanism in concentrated magmatic suspensions. Finally, the formation times in natural olivine adcumulates have been calculated. The last series of centrifuge experiments deals with the crystal-melt settling-floating mechanism in a system composed of natural two pyroxene gabbro. The results have revealed a vertical evolution of the major and trace elements in the melt phase. Then, a numerical modelling of the sedimentation process of the crystals has been made in order to describe the compaction evolution with time. In comparing the numerical simulation with the centrifuge modelling, the stratification of the compacted layer in the runs is reproduced in numerical models. Moreover, on the base of the numerical and centrifuge modelling, a sedimentation exponent describing a deviation of settling in concentrated suspensions from Stokes sedimentation has been evaluated. Finally, the numerical simulation is applied to the Muskox intrusion to estimate the formation time and the melt fraction evolution in using the hindered sedimentation model calculations.
In dieser Arbeit wurde der chemische Ozonverlust in der arktischen Stratosphäre über elf Jahre hinweg, zwischen 1991 und 2002, mit Hilfe der so genannten "Ozon-Tracer Korrelationstechnik" (TRAC), untersucht. Bei dieser Methode werden Korrelationen zwischen Ozon und langlebigen Spurenstoffen im Verlauf des Winters im Polarwirbels beobachtet und so der jährliche akkumulierte Ozonverlust berechnet. Die Ergebnisse dieser Arbeit basieren im wesentlichen auf Messdaten der Satelliteninstrumente: HALOE (Halogen Occultation Experiment) auf UARS (Upper Atmosphere Research Satellite) und ILAS (Improved Limb Atmospheric Spectrometer) Instrument auf ADEOS (Advanced Earth Observing Satellite). Das HALOE Instrument misst seit Oktober 1991 kontinuierlich alle zwei bis drei Monate für einige Tage in höheren nördlichen Breiten. ILAS lieferte ausschließlich für den Winter 1996-97 Messungen, die über sieben Monate hinweg in hohen Breiten aufgenommen wurden. Aufgrund der eingeführten Erweiterungen und Verbesserungen der Methode in dieser Arbeit, konnte die Methode anhand einer detaillierten Studie für den Winter 1996-97 validiert werden. Die ILAS Messreihe wurde dazu verwendet, erstmals die Untersuchung der zeitlichen Entwicklung von Ozon-Tracer Korrelationen kontinuierlich für die gesamte Lebensdauer des Polarwirbels durchzuführen. Dabei wurden auch Korrelationen während der Bildung des Wirbels untersucht und im Besonderen mögliche Mischungsvorgänge zwischen Wirbelluft und Luftmassen außerhalb des Wirbels. Ausserdem wurde ein Vergleich der Ergebnisse von ILAS und HALOE Messdaten durchgeführt und Unterschiede in den Ergebnissen tiefgreifend analysiert. Basierend auf HALOE Messungen konnte die erweiterte TRAC Methode über elf Jahren hinweg angewendet werden. Damit war erstmals eine konsistente Analyse von Ozonverlust und Chloraktivierung über diesen Zeitraum möglich. Die Erweiterungen führten zu einer Verringerung und genauen Quantifizierung von Unsicherheiten der Ergebnisse. Ein deutlicher Zusammenhang zwischen meteorologischen Bedingungen, Chloraktivierung und dem chemischen Ozonverlust wurde deutlich. Weiterhin zeigte sich eine Abhängigkeit zwischen den meteorologischen Bedingungen und der Homogenität des Ozonverlustes innerhalb eines Winters, sowie der mögliche Einfluss von horizontaler Mischung auf Luftmassen in einem schwach ausgeprägten Polarwirbel. In dieser Arbeit wurde eine positive Korrelation zwischen den über die gesamte Lebensdauer des Wirbels auftretenden möglichen PSC-Flächen und den akkumulierten Ozonverlusten für die elf untersuchten Jahre deutlich. Es konnte darüber hinaus gezeigt werden, dass der Ozonverlust von deutlich mehr Einflüssen als nur von der Fläche möglichen PSC Auftretens bestimmt wird, sondern zum Beispiel von der Stärke der Sonneneinstrahlung abhängt. Außerdem lassen sich Auswirkungen von Vulkanausbrüchen, wie zum Beispiel im Jahr 1991 der des Mount Pinatubo, identifizieren.
Climatology of morphology and cloud-radiative properties of marine low-level mixed-phase clouds
(2023)
Marine stratocumuli cover about 40 - 60% of the ocean surface. They self-organize into different morphological regimes. The two organized cellular regimes are called open and closed mesoscale-cellular convective (MCC) clouds. In mid-to-high latitudes, open and closed cells are the two most frequent types of MCC clouds. In particular, many MCC clouds consist of a mixture of vapor, liquid droplets, and ice particles, referred to as mixed-phase clouds (MPCs). Even for the same cloud fraction, the albedo of open cells is, on average, lower than that of closed MCC clouds. Cloud phase and morphology individually influence the cloud radiative effect. Thus, this thesis investigates the relationships between the cloud phase, MCC organization, cell size, and differences regarding the cloud-radiative effect.
This thesis focuses on space-borne retrievals to achieve extensive temporal and spatial coverage. The liDAR-raDAR (DARDAR) version 2 product collocates two active and one passive satellite: CloudSat, Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO), and Moderate Resolution Imaging Spectroradiometer (MODIS). The cloud phase of DARDAR is vertically integrated to establish a single cloud phase at each data point. The MCC classification data set based on the liquid water path (LWP) of MODIS scenes is collocated with the DARDAR product to determine the MCC organization. Cell-size statistics of both MCC clouds are obtained using a marker-based image segmentation method on MODIS reflectance scenes. In addition, based on MODIS reflectance scenes, a convolutional neural network (CNN) is developed to classify open and closed MCC scenes to avoid missing mature MPCs with a low LWP.
The first part of this thesis explores the relationships between cloud phase, morphology, and cloud albedo in the Southern Ocean (SO). At a given cloud-top temperature (CTT), seasonal changes in the mixed-phase fraction, defined as the number of MPCs divided by the sum of MPC and supercooled liquid cloud (SLC) pixels, are stronger than the morphological changes. Therefore, external factors seem to influence these changes instead of morphology. The dependence of cloud phase on cloud-top height (CTH) is more substantial than on CTT in clouds with CTHs below 2.5 km. The previously observed acceleration of closed-to-open transition in MPCs, known as preconditioning, is not the primary driver of climatological cloud morphology statistics in the SO. The morphological differences in cloud albedo are more pronounced in SLCs than in MPCs. This change in albedo alters the cloud radiative effect in the SO by 21Wm−2 to 39Wm−2 depending onseason and cloud phase.
Open and closed MCC clouds exhibit larger equivalent cell diameters in the MPCs than in SLCs in austral summer, whereas, in austral winter, the SLCs are larger. The cell’s aspect ratio accounts for varying CTHs. Closed cells have smaller aspect ratios than open cells, so their cell diameter is smaller, independent of CTH. While the seasonal differences in closed cells are due to changes in CTH, the seasonal aspect ratio differences in open cells are mainly caused by MPCs. With increasing aspect ratios, the cloud albedo decreases in both open and closed MCC clouds, with the most substantial decrease in open MPCs clouds. This leads to cloud-radiative changes of 60 - 75Wm−2 in the SO, depending on cloud phase and aspect ratio.
The established CNN exhibits a good accuracy of 80.6%, with even higher accuracies in the Open (85.5%) and Closed (87.3%) categories. The global MCC climatology based on the CNN generally agrees well with previous MCC distributions. The most notable difference occurs in the Northern Hemisphere (NH) in boreal winter, with a higher occurrence frequency of closed and open MCC clouds. This might indicate missing MPCs in previous studies based on the LWP and some restricted to warm cloud scenes. Thus, the developed CNN seems to better represent the different morphologies in MPCs than in previous classifications.
In conclusion, this thesis shows that understanding the dependencies of cloud phase, cloud morphology, and cell size is important to enhance predictions of the cloud-radiative effect and thus, it is important to evaluate how cloud phase, cloud morphology, and cellsize change in a warming climate.
During this study clumped isotope analysis of carbonates was established at the Goethe University of Frankfurt, Germany. Therefore, preparation protocols and analytical parameters were elaborated to obtain precise and accurate Δ47 data. Briefly, analyte CO2 was cleaned cryogenically using glass extraction lines to remove traces of water that enable re-equilibration of C–O bonds in the gases. Furthermore, analyte CO2 was passed through a gas chromatograph (GC) to clean it from contaminants that produce isobaric interferences with m/z 47. Initially, phosphoric acid digestions of carbonates was conducted at 25 °C in McCrea-type reaction vessels. Afterwards samples were reacted at 90 °C using a common acid bath. Mass spectrometric analyses were performed using a MAT 253 equipped with a dual inlet system. Δ47 values were directly projected to the absolute scale using CO2 gases equilibrated at distinct temperatures.
In cooperation with Stefano Bernasconi and his research group at ETH Zurich we studied the non-linearity that occurs for the measurement of m/z 47. This effect results from secondary electrons created by the m/z 44 beam. These electrons cause a negative background on the m/z 47 collector. A correction procedure was proposed that relies on the determination of the negative background on the m/z 47 Faraday cup. This approach might reduce time-consuming analyses of heated gases which were used so far to account for the observed non-linearity. However, the suggested correction of the negative background on the m/z 47 cup is only applicable if the slit width of the m/z 44 beam is significantly wider than that of the m/z 47 beam.
This thesis, furthermore, presents a comparison of the different phosphoric acid digestion techniques which are commonly used for carbonate clumped isotope analysis. For calcitic and aragonitic material digested at 25 °C in McCrea-type vessels we observed that the sample size has an effect on Δ47 data: higher mean Δ47 values and a larger scatter of data were received for samples <7 mg than for larger aliquots. For carbonate samples digested at 90 °C in a common acid bath no sample size effect was determined. We assume that secondary re-equilibration of CO2 with water preferentially occurs at 25 °C producing the observed differences. However, a sample size effect can be avoided if reaction temperature is increased to 90 °C.
In order to make carbonate Δ47 data obtained from acid digestions at 90 °C comparable to Δ47 data received from reactions at 25 °C the difference of the acid fractionation factores (Δ47*25-90) between both temperatures has to be known. For the determination of the Δ47*25-90 value we have considered Δ47 data made at 25 °C from samples >7 mg only. For calicte and aragonite we obtained differences in fractionation factores of 0.075‰ and 0.066‰, respectively. These Δ47*25-90 values are coincident with the theoretical prediction of 0.069‰ proposed for calcite (Guo et al., 2009).
Moreover, this dissertation comprises a calibration study of the clumped isotope thermometer based on various natural calcites that grew between 9 and 38 °C. The samples include a brachiopod shell, a bivalve shell, an eggshell of an ostrich and foraminifera tests which formed from distinct biomineralizing processes. Furthermore we included an authigenic carbonate crystallized from biological-induced precipitation. The following linear relationship between 1/T2 and Δ47 was determined (with Δ47 in ‰ and T in K):
Δ47 = 0.0327 (± 0.0026) x 106 / T2 + 0.3030 (± 0.0308) (R2 = 0.9915)
This equation differs from the pioneering Ghosh et al. (2006a) calibration. However, our regression line is statistically indistinguishable from that of Henkes et al. (2013) which is based on aragonitic mollusks and calcitic brachiopod shells. Both studies have in common that calibration data were, at first, directly referenced to the absolute scale. In addition, both datasets rely on similar digestion techniques. Furthermore, the two calibrations are conform with the theoretical prediction of Guo et al. (2009).
The calcite calibration of the clumped isotope paleothermometer received in this study was applied to Δ47 data measured for Silurian brachiopods shells from Gotland/Sweden. Prior to isotopic analysis the fossils were intensively investigated for their preservation state (CL, SEM, trace elements). The lowest T(Δ47) values of ca. 28 to 33 °C were estimated from ultrastructurally well-preserved regions of some shells. For these samples also the lowest δ18Ow values of Silurian seawater were determined. These estimates of ca. −1‰ confirm the assumption that the δ18O value of the Silurian ocean was buffered to (0 ± 1)‰.
Nevertheless, most studied shells were characterized by a patchwork of pristine and altered shell portions resulting in elevated T(Δ47) values which plot mostly between 40 and 60 °C. Our results indicate that the clumped isotopic composition of the shells were altered at low water-rock ratios, not affecting the δ18O values. Δ47 and δ18O data of associated diagenetic phases (sparitic and micritic phases of the inner fillings of the fossils) provide evidence that the sparitic cements grew during several diagenetic events which occurred at different temperatures in fluid-buffered systems. We, furthermore, conclude that the micritic phases lithified at a very early diagenetic stage with the δ18O values being most probably close to a Silurian seawater composition
In dieser Arbeit werden Schmelz- und Anreicherungsprozesse des Erdmantels, sowie Kristallisationsereignisse der Erdkruste zweier ausgewählter Gebiete in Namibia und Spanien mithilfe geochemischer Methoden rekonstruiert und in einen zeitlichen Zusammenhang gebracht. Ein Vergleich der gewonnenen Ergebnisse beider Kompartimente soll dabei weitere Informationen liefern inwieweit Prozesse des Erdmantels und der Erdkruste miteinander verknüpft waren. Insbesondere soll ein weitere Beitrag zur aktuellen Diskussion geliefert werden, bei der sich das sogenannte „pulsed growth“ und „steady accumulation“ Modell gegenüberstehen (siehe Zusammenstellung Pearson et al., 2007). Zudem tragen die neu gewonnenen Daten dazu bei, die regionalen geologischen Gegebenheiten im besonderen Hinblick auf die geotektonische Geschichte besser zu verstehen.
Das Gibeon Kimberlit Feld befindet sich in der tektonischen Einheit des Rehoboth Terranes in Namibia und ist gekennzeichnet von Vulkanismus vor etwa 72.5 Ma (Davies et al., 2001), der Granat Peridotite und krustale Xenolithe mit an die Oberfläche beförderte. Eine klare Einordnung des Rehoboth Terranes in die Gesamtheit des Süd Afrikanischen Plattenverbunds ist noch nicht vollständig geklärt.
Die Südöstliche vulkanische Provinz in Spanien (SEVP) mit besonderem Hinblick auf die Region um Casas de Tallante stellt das zweite Probengebiet für diese Arbeit dar. Vor etwa 2.6 Ma (Bellon et al., 1983) kam es zur Extrusion von alkali-basaltischen Schmelzen, die zahlreiche Spinell / Plagioklas Peridotite mit sich brachten. Tufflagen, sowie die Matrix der Basalte ermöglichen einen Einblick in die untere Kruste der Region.
Untersuchungen der Erdmantelproben aus Namibia auf ihre Haupt- und Spurenelementchemie, sowie Lu-Hf und Sm-Nd Isotopie zeigten, dass zwei verschiedene Manteltypen vorliegen („N“ und „σ“ Typ), die zu einem Zeitpunkt um etwa 850 Ma („N“) und 1.9 Ga („σ“) angereichert wurden. Eine letzte Anreicherung beider Typen fand vermutlich während der Pan–Afrikanischen Orogenese um etwa 450 Ma statt. Die Reinterpretation eines zuvor publizierten Datensatzes (Pearson et al., 2004), suggeriert, dass es zu einer ersten Verarmung der σ Peridotite um etwa 2.9 Ga kam.
Untersuchungen der U-Pb und Hf Isotopie an Zirkonen aus der unteren Kruste des Probengebiets in Namibia ergaben, dass es zur Bildung von juvenilem Krustenmaterial vermutlich bereits im Archaikum kam (wie bereits vorgeschlagen durch z.B. Hoal et al., 1995; Franz et al., 1996), sowie in den Zeiträumen von 2.3 bis 2.7 und 1.5 bis 1.6 Ga, mit jeweils anschließendem krustalem Recycling und Krustenmischung. Eine Übereinstimmung von Mantel- und Krustenevents konnte für die Zeiträume von etwa 1.8, 0.8 - 0.9 Ga und 0.4 – 0.5 Ga gefunden werden. Eine mögliche erste Verarmung des σ Mantels wird bestätigt durch Zirkonalter im Bereich von 2.7 bis 2.9 Ga.
Die Analyse ausgewählter Spinell / Plagioklas Peridotite aus der SEVP, ergaben, dass ein heterogener Mantel mit mindestens 3 verschiedenen Typen vorliegt. Eine Korrelation der Lu-Hf Isotopie von 3 Proben dieses Probensatzes, sowie den Hf Isotopien einer weiteren Probe von Bianchini et al. (2011) suggerieren, dass es eventuell zu einem Verarmungsereignis zu einem Zeitpunkt von etwa 550 Ma kam. Sr Isotopien von Klinopyroxenen und Plagioklasen im Vergleich ergaben, dass die Sr Isotopie der Plagioklase, im Gegensatz zu den Klinopyroxenen, von denen der Alkali Basalte überprägt wurden.
Zirkonanalysen aus Lokalitäten innerhalb der SEVP (U-Pb, Hf) ergaben ein weitreichendes Altersspektrum, beginnend bei etwa 2-3 Ma bis hin ins Archaikum (2.7 bis 2.9 Ga) mit Provenance Ursprung aus Gondwana und dem Arabisch-Nubischen Schild. Die Kombination der U-Pb Altersinformationen mit den entsprechenden Hf Isotopien, zeigten, dass es vermutlich bereits im Archaikum zu juveniler Krustenbildung kam. Zirkone > 100 µm datieren den Zeitpunkt der Eruption der Alkali Basalte mit Altern um etwa 2.6 Ma und Hf Isotopien, die einem leicht verarmten Mantel entsprechen. Ein mögliches Verarmungsereignis im Erdmantel zu einem Zeitpunkt von etwa 550 Ma, ist im Einklang mit Krustenrecycling zu selbigem Zeitpunkt.
Die neugewonnenen Daten dieser Arbeit unterstützten das „pulsed growth“ Modell.
Literatur
Bellon, H., Bordet, P. and Montenat, C., 1983. Chronology of the Neogene Magmatism from Betic Ranges (Southern Spain). Bulletin De La Societe Geologique De France, 25(2): 205-217.
Bianchini, G., Beccaluva, L., Nowell, G.M., Pearson, D.G. and Siena, F., 2011. Mantle xenoliths from Tallante (Betic Cordillera): Insights into the multi-stage evolution of the south Iberian lithosphere. Lithos, 124(3-4): 308-318.
Davies, G.R., Spriggs, A.J. and Nixon, P.H., 2001. A non-cognate origin for the Gibeon kimberlite megacryst suite, Namibia: Implications for the origin of Namibian kimberlites. Journal of Petrology, 42(1): 159-172.
Franz, L., Brey, G.P. and Okrusch, M., 1996b. Steady state geotherm, thermal disturbances, and tectonic development of the lower lithosphere underneath the Gibeon Kimberlite Province, Namibia. Contributions to Mineralogy and Petrology, 126(1-2): 181-198.
Hoal, B.G., Hoal, K.E.O., Boyd, F.R. and Pearson, D.G., 1995. Age constraints on crustal and mantle lithosphere beneath the Gibeon kimberlite field, Namibia. South African Journal of Geology, 98(2): 112-118.
Pearson, D.G., Irvine, G.J., Ionov, D.A., Boyd, F.R. and Dreibus, G.E., 2004. Re-Os isotope systematics and platinum group element fractionation during mantle melt extraction: a study of massif and xenolith peridotite suites. Chemical Geology, 208(1-4): 29-59.
Pearson, D.G., Parman, S.W. and Nowell, G.M., 2007. A link between large mantle melting events and continent growth seen in osmium isotopes. Nature, 449(7159): 202-205.
Derivation and characterization of a new filter for nonlinear high-dimensional data assimilation
(2015)
Data assimilation (DA) combines model forecasts with real-world observations to achieve an optimal estimate of the state of a dynamical system. The quality of predictions in nonlinear and chaotic systems such as atmospheric or oceanic circulation is strongly sensitive to the initial conditions. Therefore, beyond the consistent reconstruction of past states, a primary relevance of advanced DA methods concerns the proper model initialization. The ensemble Kalman filter (EnKF) and its deterministic variants, mostly square root filters such as the ensemble transform Kalman filter (ETKF), represent a popular alternative to variational DA schemes. They are applied in a wide range of research and operations. Their forecast step employs an ensemble integration that fully respects the nonlinear nature of the analyzed system. In the analysis step, they implicitly assume the prior state and observation errors to be Gaussian. Consequently, in nonlinear systems, the mean and covariance of the analysis ensemble are biased and these filters remain suboptimal. In contrast, the fully nonlinear, non-Gaussian particle filter (PF) relies on Bayes' theorem without further assumptions, which guarantees an exact asymptotic behavior. However, it is exposed to weight collapse, particularly in higher-dimensional settings, known as the curse of dimensionality.
This work presents a new method to obtain an analysis ensemble with mean and covariance that exactly match the corresponding Bayesian estimates. This is achieved by a deterministic matrix square root transformation of the forecast ensemble, and subsequently a suitable random rotation that significantly contributes to filter stability while preserving the required second-order statistics. The forecast step remains as in the ETKF. The algorithm, which is fairly easy to implement and computationally efficient, is referred to as the nonlinear ensemble transform filter (NETF). The limitation with respect to fully-nonlinear filtering is that the NETF only considers the mean and covariance of the Bayesian analysis density, neglecting higher-order moments.
The properties and performance of the proposed algorithm are investigated via a set of experiments. The results indicate that such a filter formulation can increase the analysis quality, even for relatively small ensemble sizes, compared to other ensemble filters in nonlinear, non-Gaussian scenarios. They also confirm that localization enhances the applicability of this PF-inspired scheme in larger-dimensional systems. Finally, the novel filter is coupled to a large-scale ocean general circulation model with a realistic observation scenario. The NETF remains stable with a small ensemble size and shows a consistent behavior. Additionally, its analyses exhibit low estimation errors, as revealed by a comparison with a free ensemble integration and the ETKF. The results confirm that, in principle, the filter can be applied successfully and as simple as the ETKF in high-dimensional problems. No further modifications are needed, even though the algorithm is only based on the particle weights. Thus, it is able to overcome the curse of dimensionality, even in deterministic systems. This proves that the NETF constitutes a promising and user-friendly method for nonlinear high-dimensional DA.
Development of the flash-heating method for measuring melting temperatures in the diamond anvil cell
(2016)
A new ‘laser flash-heating’ method has been developed for measuring melting temperatures above 2000 K in a diamond anvil cell at gigapascals of pressure. It overcomes the general difficulties in detecting an onset of melting in a diamond anvil cell. It also circumvents the notorious experimental difficulties associated with the long heating durations of the CW laser-heating and the short timescales in the pulsed laser-heating and shock-compression experiments.
In this method, the duration of heating a sample is tuned to avoid chemical reactions of the sample with the diamond anvils and the surrounding pressure medium, while maintaining the accuracy of the temperature measurements. The absence of chemical reactions is confirmed by the EDS technique. Melt detection is now unambiguous from the analysis of textures on the surface and in depth of the recovered samples using the SEM and FIBM techniques, respectively. Using this method, the following has been achieved.
1. The melting curve of hcp-Re has been measured to 48 GPa, 4200 K for the first time. It has a significantly steeper slope than those observed for other transition metals like W and Mo with bcc structures. Above 20 GPa, Re becomes the most refractory metal surpassing W.
2. The melting curve of bcc-Mo has been measured to 45 GPa, 3100 K. It agrees with previous melt-slopes approaching zero value with pressure as reported in the LHDAC experiments using ADXRD and visual observation techniques for inferring the onset of melting. Flash-heating experiments at pressures higher than 50 GPa are required to further corroborate the flat melt-slope and resolve the long standing controversy about melting of Mo.
3. The melting curve of bcc-Ta has been measured to 85 GPa, 4300 K. Unlike in previous experiments using ADXRD and visual observation as probes, it has been tightly bracketed with an unambiguous detection of the onset of melting, without any chemical reaction. The present melting curve cannot be reconciled with shock measurements and theoretical predictions, and the precision of measurements calls for a reevaluation of theoretical, shock compression, and other DAC approaches to determine melting at high pressures. A further analysis with TEM technique for investigating the structure of the heated portion below and above melting temperatures of Ta may benefit in resolving various phase transitions predicted to explain the vast discrepancies in the reported melt-slopes.
When extrapolated to one atmosphere pressure, all the measured flashmelting curves agree with the known melting points.
The development of benthic foraminiferal assemblages during the past 6,000 yrs was investigated in Holocene sediment cores from three carbonate platforms (Turneffe Islands, Lighthouse Reef, and Glovers Reef) of Belize, Central America. Foraminiferal assemblages and their diversity were determined in different time periods to identify their dependence on environmental factors, such as lagoonal age, lagoonal depth, water circulation, substrate, bottom-water temperature, and salinity. Geochemical proxies (δ18O and δ13C), obtained from the common larger foraminifer Archaias angulatus were used to estimate Holocene seasonal BW-temperatures and climate variabilities. A total of 51 samples were taken from 12 vibracores for taxonomic determination and 10 to 15 subsamples of 32 tests of Archaias angulatus were used for stable oxygen and carbon isotope analyses. Based on cluster analyses, seven benthic foraminiferal assemblages are distinguished during the Holocene. The three platforms exhibit characteristic differences in benthic foraminiferal fauna and diversity, which are controlled by their respective environments during the last 6,000 yrs. Turneffe Islands has four benthic foraminiferal assemblages, which are typical for restricted lagoons with fluctuating salinity. Lighthouse Reef is inhabited by two benthic foraminifera associations, which are characteristic of high water exchange with the surrounding ocean and clear waters. Glovers Reef is characterized by two benthic foraminiferal assemblages, which occur in deeper lagoons with slow water circulation. In general, during the Holocene, the highest mean diversity, evenness, and richness of benthic foraminifera were found in the Turneffe Islands and the lowest occurred at Glovers Reef. The foraminiferal faunas of the Lighthouse and Glovers Reefs had been in a “Diversification Stage” since 6,000 yrs, whereas the foraminiferal fauna of the Turneffe Islands reflects the development from a “Colonisation” (~4,000 yrs BP) to a “Diversification Stage” (~2,000 yrs to present time). Lagoonal depth, water circulation, substrate, and BW-temperature have higher influence on foraminiferal diversity as compared to lagoonal size and age. The negative correlation between diversity and lagoonal depth is based on differences in light intensity and substrate. In contrast to Lighthouse Reef, the Turneffe Islands and Glovers Reef show decreasing diversity of benthic foraminifera with increasing lagoon depth, due to finer sediment, turbid waters and/or dense mangrove growth, which reduce the light intensity and the number of species. Water Circulation also affected the benthic foraminifera modes of living and their diversity during the last 6,000 yrs. Increasing abundances of infaunal taxa refer to restricted circulation and/or lower oxygen conditions, as assumed for the Turneffe Islands and Glovers Reef. Increasing abundances of epifaunal foraminifera, as observed in the Lighthouse Reef indicate better circulation and/or higher oxygen conditions. Holocene BW-temperature reconstructions based on δ18O of single Archaias angulatus tests do not correspond to typical Holocene climate models of the Caribbean. In the Belize area, mean BW-temperature trends indicate local climate variations. A decrease of δ13C values during the last 1,000 yrs could be related to the “Suess Effect”. The seasonal BW-temperature variations within single large benthic foraminifera tests correspond to present-day temperature fluctuations in the lagoons, and indicate higher temperatures in Summer and Autumn and lower temperatures in Winter and Spring.
Mantle convection is the process by which heat from the Earth’s core is transferred upwards to the surface and it is accepted to explain the dynamics of the Earth’s interior. On geological time-scales, mantle material flows like a viscous fluid as a consequence of the buoyancy forces arising from thermal expansion. Indeed, mantel convection provides a framework which links together the major disciplines, such as seismology, mineral physics, geochemistry tectonic and geology. The numerical model has been applied to understand the dynamic, structure and evaluation of the Earth, and other terrestrial planets and the investigations continue to explore, different aspects of the mantle convection.
In fact, to model this phenomenon, two complementary approaches are possible. On the one hand, one can solve self-consistently the equations of thermal convection, including parameters and employing physical relationships derived from mineral physics. Our understanding of mantle convection depends ultimately upon the success of such fully self-consistent dynamic models in explaining observable features of the flow. Although, these models presently unable to predict the actual convection pattern of the Earth, they are extremely useful to investigate general characteristics of given physical systems. On the other hand, to permit comparison with specific observables associated with the flow, one can consider a more restricted problem. Instead of focusing on the time evolution of mantle flow, if we know a priori the temperature - and hence presumably the density - anomalies that drive the convection, we can try to build a snapshot of the present-day flow pattern, consistent with those anomalies, that can successfully predict the observables. As matter of fact, the aim of this study is to investigate both approaches in comparison with the main geophysical constraints on mantle structure. These constraints include the geoid anomalies, the dynamic surface and core-mantle boundary topography and tectonic plate motions.
The most appropriate mathematical basis functions for describing a bounded and continuous function on a spherical surface are spherical harmonics. We may therefore expand the geodynamic observables in terms of spherical harmonics. We have investigated two methods of the global spherical harmonic analysis by specific attention to the dynamic geoid computation of the geodynamic models. The first method is the quadrature method in which the loss of the orthogonality of the Legendre functions in transition from continues to discrete case is the major drawback to the method. Particularly, we showed that in the absence of the tesseral harmonics, quadrature formulation leads to obtain inaccurate results. The second method is the least-squares which can be considered as the best linear unbiased estimator that provides the exact results. We showed that even with a low resolution grid data it is possible to reconstruct the data and achieve an accurate result by using this method, which is extremely remarkable in three-dimensional global convection studies. However, special care has to be taken since there is some source of errors that might influence the efficiency of this method.
In general, to better understanding of the properties of the mantle, it is useful to assess observable characteristics of plumes in the mantle, including geoid, topography and heat flow anomalies. However, only few studies exist on geoid and topography for axi-symmetric convection and their models were restricted to isoviscous (or stratified) mantle and low Rayleigh numbers. We studied fully coupled depth and temperature dependent Arrhenius type of viscosity in axi-symmetric spherical shell geometry in order to investigate the shape of geoid anomalies and dynamic topography above a plume. Indeed, the topography and geoid anomalies produced from plumes are sensitive to rheology of the mantle and rheology of the plume; both have effects on shape and amplitude of the geoid anomalies. As results we are able to define different classes of plumes by their geoid signals.
Mainly depth-dependent viscosity models show a geoid with negative sign above the plume which can turn to the positive sign by decrease the viscosity contrast. This can be considered as a transition between the strongly depth dependent and the constant viscosity case. Our results basically support the idea by Morgan [1965] and McKenzie [1977]. They have shown the magnitude and even the sign of the total gravity anomaly depend on the spatial variation in effective viscosity. In addition, Hager [1984] has concluded that the total gravity field is depend on the radial distribution of effective viscosity, and a small change in viscosity contrast leads to varying sign of the response function.
In the case of temperature-dependent viscosity, the formation of an immobile lithosphere is a natural result, and the flow as well as the total geoid becomes strongly time dependent. When we increase the activation energy, all geoids associated with the first arriving plumes look like bell shaped whereas for typical plumes, after reaching a statistical steady state, bell-shaped geoids with decreasing amplitude as well as linear flank shaped geoids are observed. It is surprising that in spite of large differences in lateral and depth varying viscosities, the shapes of the geoid anomalies remained rather similar. We also identified different behaviors in the combined model with temperature-and pressure-dependent viscosity. In fact, in spite of the strongly different rheology, the geoid anomalies in all cases were surprisingly similar. Furthermore, we proposed a scaling law for the geoid which makes our results directly applicable to other planets. Moreover, we can apply the results of our calculation to find relations between different rheology and sub-lid temperature, since we know that the mantle temperature can change significantly with variation in pressure-temperature dependent viscosity. It is also possible to define a range of stagnant lid thickness related to the amplitude of the geoid which can be reasonable for study of the lid thickness in Venus or Mars.
Nevertheless, in these series of models, we simplified a number of complexities within the Earth. One of the most important of such simplification is the Boussinesq approximation. This approximation is valid if the temperature scale height (i.e. the depth over which temperature increases by a factor of “ ” due to adiabatic compression) is much greater than the convection depth. However, a temperature scale height in the Earth’s mantle is at best only slightly greater than the mantle depth. Hence, the Boussinesq approximation could mask some very important stratification and compressibility effects that influence both the spatial and temporal structure of the convection. Therefore, in more advance models we considered compressibility in our mantle convection models, assuming that density vary both radially and laterally, being determined as a function of pressure and temperature through an appropriate equation of the state. Moreover, thermodynamic properties assumed to be a function of depth.
We examined the details of the structure of the spherical axi-symmetric Anelastic Liquid Approximation model (ALA) with special attention to the Arrhenius rheology, and compare it to the cases of compressible convection without depth dependent thermodynamical properties, and to cases of the extended Boussinesq approximation. At the same time, the effects of the interaction between temperature and pressure-dependent viscosity and thermodynamic parameters in the compressible mantle convection on the geoid and topography have been studied. We showed that assuming compressible convection with depth-dependent thermodynamic properties strongly influence the geoid undulations. Using compressible convection with constant thermodynamic properties is physically inconsistent and may lead to spurious results for the geoid and convection pattern. Indeed, by a systematic study of different approaches of compressibility in the spherical shell convection for different Arrhenius viscosity laws we proved that only in the unrealistic case of zero activation energy the different compressibility modes result in comparable convection and geoid patterns. In all other rheological cases, large differences have been obtained, that stressing the important role of consistent compressible thermodynamic properties for mantle convection.
In addition, we examine the impact of compressibility as well as different rheologies on the power law relation that connects the Nusselt number to the Rayleigh number. We have discovered that the power law index of the relationship is controlled by the rheology, independent of which approximation is used. Instead, the bound of this relation is controlled by a combination of different approximation and rheology.
Next, instead of focusing on the time evolution of mantle flow, we have carried out three-dimensional spherical shell models of mantle circulation to investigate the effects of joint radial and lateral viscosity variations on the Earth’s non-hydrostatic geoid, surface and core-mantle boundary topographies. These models include realistic lateral viscosity variations (LVV) in the lithosphere, upper mantle and lower mantle in combination with different stratified viscosity structures. We have demonstrated that the contradictory results concerning the effects of LVV can be clarified by the most straight-forward problem in geoid modeling; namely, rather poorly known stratified viscosity structure. We explored three classes of dynamic geoid models due to lateral viscosity variations. In the first class, the LVV strongly improved the fit to the observed geoid. Indeed, when the viscosity contrast between lower and upper mantles is not large enough to produce a good fit to geoid the LVVs are able to perform this action by adjusting amplitudes, so that it becomes comparable with observation. In the second class, inducing the LVV moderately improved the fit. Actually, when the geoid induced by a stratified viscosity structure already has a good correlation with observation, then the LVV causes its amplitude to further improve. In the last class, if the viscosity contrast between upper and lower mantle would be high enough, inducing LVV deteriorate the fit to the observed geoid.. Indeed, depending on the stratified viscosity, inducing the LVV may take place in one of these categories.
We also quantified the effects of LVV in the mantle and lithosphere individually. We found that the presence of LVV in the mantle (upper and lower) improves the fit to the observed geoid regardless of stratified viscosity. While LVV in the lithosphere is a crucial parameter, and dependent of the stratified viscosity, may increase or decrease the geoid fit. In fact, when the lower mantle considers being viscous enough, it would support the negative buoyancy of subducting slabs. Thus, it transmits some of the stress back to the top boundary and causes a weak coupling between slab and surface. Therefore, by including the low viscous plate boundaries in this model, the slabs and overriding plates decouples and the fit to the observed geoid degrades. In contrast, when the lower mantle viscosity is not sufficiently stiff, the presence of the low viscous plate boundaries assists to weaken the strong mechanical coupling between slab and surface. Hence, a better fit achieved.
Water is scarce in semi-arid and arid regions. Using alternative water sources (i.e. non-conventional water sources), such as municipal reuse water and harvested rain, contributes to using existing water resources more efficiently and productively. The aim of this study is to evaluate the two alternative water sources reuse water and harvested rain for the irrigation of small-holder agriculture from a system perspective. This helps decision and policy makers to have proper information about which system and technology to adopt under local conditions. For this, the evaluation included ecologic, societal, economic, institutional and political as well as technical aspects. For the evaluation, the study area in central-northern Namibia was chosen in the frame of the research and development project CuveWaters. The main methods used include a mathematical material flow analysis, the computation and modelling of crop requirements, a multi-criteria decision analysis using the Analytical Hierarchy Process (AHP) method and a financial cost-benefit analysis. From a systemic perspective, the proposed novel systems were compared to the exciting conventional infrastructure. The results showed that both water reuse and rainwater harvesting systems for the irrigation of small-holder horticulture offer numerous technological, ecologic, economic, societal, institutional and political benefits. Rainwater harvesting based gardens have a positive benefit-cost ratio under favorable conditions. Government programs could fund the infrastructure investment costs, while the micro-entrepreneur can assume a micro-credit to finance operation and maintenance costs. Installing sanitation in informal settlements and reusing municipal water for irrigation reduces the overall water demand of households and agriculture by 39%, compared to improving sanitation facilities in informal settlements without reusing the water for agriculture. Given that water is the limiting factor for crop fertigation, the generated nutrient-rich reuse water is sufficient to annually irrigate about 10 m2 to 13 m2 per sanitation user. Compared to crop nutrient requirements, there are too many nutrients in the reuse water. Thus when using nutrient-rich reuse water, no use of fertilizers and a careful salt management is necessary. When comparing this novel system with improved sanitation, advanced wastewater treatment and nutrient-rich water reuse to the conventional and to two adapted systems, results showed that the novel CuveWaters system is the best option for the given context in a semi-arid developing country. Therefore, the results of this study suggest a further roll-out of the novel CuveWaters system. The methodology developed and the results of this study demonstrated that taking sanitation users into consideration plays a major role for the planning of an integrated water reuse infrastructure because they are the determinant factor for the amount of available nutrient-rich reuse water. In addition, it could be shown that water reuse and rainwater harvesting systems for the irrigation of small-scale gardens provide a wide range of benefits and can be key to using scarce water resources more efficiently and to contributing to the Sustainable Development Goals.
The thesis is devoted to the study of the Antarctic polar vortex, mainly by analyzing data collected during APE-GAIA (1999) and ASHOE (1994) campaigns and recorded by the ADEOS satellite (1996-1997), and to improvement of the chromato-graphic processing schemes. A general introduction and overview of the campaigns and instruments relevant to the present work are given in Chapters 1 and 2. A relatively large part of the thesis (Chapters 3-5) is on improvement of the analysis of raw chromatographic data recorded during in-flight measurements of the trace gases. A Gaussian non-straight-base-line method, i.e. the Gaussian processing scheme (Chapter 3), is developed for better evaluation of the chromatographic peak size. Furthermore, a statistical cross-correlation method (Chapter 5) based on statistical behaviour of the whole chromatogram series fNchrg recorded, e.g., during a research flight or laboratory calibration, is developed and applied to measure the low-concentration trace gases. As demonstrated for HAGAR's chromatograms (HAGAR - High Altitude Gas Analyzer), the combination of the Gaussian fitting scheme for individual chromatograms and the statistical cross-correlation method for a series of subsequent chromatograms considerably improves and stabilizes quantitative analysis of in-flight chromatographic data. In this case, the detection accuracy of weak and noisy chromatographic signals can be improved by up to 40 %. A particular attention is paid to the in-flight two-standard calibration method. For this method, a special procedure, that allows to evaluate and effectively remove a weak background chromatographic signal associated with residual molecules in the carrier gas N2, is proposed and coded (Chapter 4). The developed approaches and methods are completely automized and, therefore, can be used for processing of in-flight chromatograms of recent and future field campaigns. The main part of the thesis (Chapters 6-8) deals with a two-dimensional quasi-Lagrangian coordinate system ... , based on a long-lived stratospheric trace gas i, and its systematic use for i = N2O in order to describe the structure of a well-developed Antarctic polar vortex, linearization and compactization of the tracer-tracer correlations in the polar vortex core (i.e. the stratospheric dynamics in this area), and the differential ozone losses in the Antarctic polar vortex area. In the coordinate system ... (...-method, Chapter 6), which refers to a well-developed polar vortex, the mixing ratio Âi is the vertical coordinate and ... = .... i is the reference profile in the vortex core) is the meridional coordinate. The quasi-Lagrangian coordinates ... are much more long-lived comparing with the standard quasi-isentropic coordinates, potential temperature ... and equivalent latitude ..e, do not require explicit reference to geographic space, and therefore well-suited for studying the dynamics of the Antarctic polar vortex and the relevant ozone loss processes. By using the introduced coordinate system ... to analyze the well-developed Antarctic vortex investigated in the APE-GAIA campaign, it is shown, in concurrence with the conclusion of A. M. Lee et al. (2001), that the Antarctic vortex area can be described in terms of the well-mixed and well-isolated vortex core, relatively wide vortex boundary region and adjoining surf zone. In this case, the reference profile ... i , which is compact in a well-developed and isolated polar vortex core [J. B. Greenblatt et al. (2002)], can be found by combining airborne (and/or balloon) data with high-altitude satellite measurements. A criterion, which uses the local in-situ measurements of Âi = Âi(£) and attributes the inner vortex edge to a rapid change (±-step) in the meridional pro¯le of the mixing ratio..., is developed in Chapter 6 to determine the (Antarctic) inner vortex edge. In turn, the outer vortex edge of a well-developed Antarctic vortex is proposed to attribute to the position of a local maximum of ...H2O in the polar vortex area. For a well-developed Antarctic vortex, the ...-parametrization of tracer-tracer correlations allows to distinguish the tracer-tracer inter-relationships in the vortex core, vortex boundary region and surf zone (Chapter 7). This is clearly illustrated by analyzing the tracer-tracer relationships Âi ¡ ÂN2O obtained from the in-situ data of the APE-GAIA campaign for i = CFCl3 (CFC-11), CF2Cl2 (CFC-12), CBrClF2 (H-1211) and SF6. The solitary anomalous points in the ...CFC11 ¡ ÂN2O correlation, observed in the Antarctic vortex core during the APE-GAIA and ASHOE campaigns, are interpreted in terms of small-scale localized differential descent. As detailed in Chapter 8, the quasi-Lagrangian coordinate system fÂN2O; ¢ÂN2Og is an effective tool for evaluation of the differential ozone losses in the polar vortex area. With this purpose, a two-parametric reference function ...O3 = F(...), which characterizes the unperturbed O3 distribution in the early winter polar vortex area, is introduced to separate and quantify in terms of the meridional coordinate ...2O the differential ozone losses in the vortex core and vortex boundary region. The method is applied to analyze the ozone depletion in the Antarctic stratosphere during the austral spring 1999 (APE-GAIA campaign). In Chapter 9, the main results of the thesis are summarized.
Bayesian Networks are computer-based environmental models that are frequently used to support decision-making under uncertainty. Under data scarce conditions, Bayesian Networks can be developed, parameterized, and run based on expert knowledge only. However, the efficiency of expert-based Bayesian Network modeling is limited by the difficulty in deriving model inputs in the time available during expert workshops. This thesis therefore aimed at developing a simple and robust method for deriving conditional probability tables from expert estimates in a time-efficient way. The design and application of this new elicitation and conversion method is demonstrated using a case study in Xinjiang, Northwest China. The key characteristics of this method are its time-efficiency and the approach to use different conversion tables based on varying levels of confidence. Although the method has its limitations, e.g. it can only be applied for variables with one conditioning variable; it provides the opportunity to support the parameterization of Bayesian Networks which would otherwise remain half-finished due to time constraints. In addition, a case study in the Murray-Darling Basin, Australia, is used to compare Bayesian Network types and software to improve the presentation clarity of large Bayesian Networks. Both case studies aimed at gaining insights on how to improve the applicability of Bayesian Networks to support environmental management.
One possible approach to study systematically the influence of the deformation regime on the geometry of geological structures like folds and boudins is analogue modelling. For a complete understanding of the resulting structures, consideration of the third dimension is required. This PhD study deals with scaled analogue modelling under constriction and plane-strain conditions to improve our knowledge of folding and boudinage of lower crustal rocks in space and time. Plasticine is an appropriate analogue material for rocks in the lower crust. Therefore, this material was used for the experiments. The macroscopic behaviour of most types of plasticine is quite similar to rocks undergoing strain-rate softening and strain hardening regardless of the different microscopic aspects of deformation. Therefore, if one is aware that the stress exponent and viscosity increase with increasing strain, the original plasticine types used with stress exponents ranging from 5.8 to 8.0 are adequate for modelling geologic structures. The same holds for plasticine/oil mixtures. Thus, plasticine and plasticine/oil mixtures can be used to model the viscous flow of different rock types in the lower crust. If climb-accommodated dislocation creep and associated steady-state flow is assumed for the natural rocks, the plasticine/oil mixtures should be used, which flow under steady-state conditions. Three different experimental studies of plane-strain coaxial deformation of stiff layers, with viscosity η2 and stress exponent n2, embedded in a weak matrix, with viscosity η1 and stress exponent n1, have been carried out. The undeformed samples (matrix plus layer) were cubes with an edge length of 12 cm. All experimental runs have been carried out at T = 25 ± 1°C and varying strain rates ė, ranging from 7.9 x 10 high -6 s high -1 to 1.7 x 10 high -2 s high -1, until a finite longitudinal strain of 30% – 40% was achieved. The first experimental study improved the understanding about the evolution of folds and boudins when the layer is oriented perpendicular to the Y-axis of the finite strain ellipsoid. The rock analogues used were Beck’s green plasticine (matrix) and Beck’s black plasticine (competent layer), both of which are strain-rate softening modelling materials with stress exponent n = ca. 8. The effective viscosity η of the matrix plasticine was changed by adding different amounts of oil to the original plasticine. At a strain rate ė of 10 high -3 s high -1 and a finite strain e of 10%, the effective viscosity of the matrix ranges from 1.2 x 10 high 6 to 7.2 x 10 high 6 Pa s. The effective viscosity of the competent layer has been determined as 4.2 x 10 high 7 Pa s. If the viscosity ratio is large (> ca. 20) and the initial thickness of the competent layer is small, both folds and boudins develop simultaneously. Although the growth rate of the folds seems to be higher than the growth rate of the boudins, the wavelength of both structures is approximately the same as is suggested by analytical solutions. A further unexpected, but characteristic, aspect of the deformed competent layer is a significant increase in thickness, which can be used to distinguish plane-strain folds and boudins from constrictional folds and boudins. In the second experimental study, the impact of varying strain rates on growing folds and boudins under plane strain have been investigated. The strain rates used range from 7.9 x 10 high -6 s high -1 to 1.7 x 10 high -2 s high -1. The stiff layer and matrix consist of non-linear viscous Kolb grey and Beck’s green plasticine, respectively, both of which are strain-rate softening modelling materials with power law exponents (n) and apparent viscosities (η) ranging from 6.5 to 7.9 and 8.5 x 10 high 6 to 7.2 x 10 high 6 Pa s, respectively. The effective viscosity (η) of the matrix plasticine was partly modified by adding oil to the original plasticine. At the strain rates used in the experiments the viscosity ratio between layer and matrix ranges between 3 and 10. Different runs have been carried out where the layer was oriented perpendicular to the principal strain axes (X>Y>Z). The results suggest a considerable influence of the strain rate on the geometry of the deformed stiff layer including its thickness. This holds for every type of layer orientation (S ┴ X, S ┴ Y, S ┴ Z). If the stiff layer is oriented perpendicular to the short axis Z of the finite strain ellipsoid, the number of the resulting boudins and the thickness of the stiff layer increase, whereas the length of boudins decreases with increasing strain rate. If the stiff layer is oriented perpendicular to the long axis, X, of the finite strain ellipsoid, enlargement of the strain rate results in increasing wavelength of folds, whereas the number of folds and the degree of thickening of the stiff layer decreased. If the stiff layer is oriented perpendicular to the intermediate Y-axis of the finite strain ellipsoid enlargement of the strain rate results in a decreasing number of boudins and folds associated with increasing wavelengths of both structures. The wavelength of folds is approximately half of the boudins wavelength. This is true for the case where folds and boudins develop simultaneously (S ┴ Y) and for cases where both structures develop independently (folds at S ┴ X and boudins at S ┴ Z). In the third experimental study, scaled analogue experiments have been carried out to demonstrate the growth of plane-strain folds and boudins through space and time. Previous 3D-studies are based only on finite deformation structures. Their results can therefore not be used to prove if both structures grew simultaneously or in sequence. Plane strain acted on a single stiff layer that was embedded in a weak matrix, with the layer oriented perpendicular to the intermediate Y-axis of the finite strain ellipsoid. Two different experimental runs have been carried out using computer tomography (CT) to analyse the results. The first run was carried out without interruption. During the second run, the deformation was stopped in each case at longitudinal strain increments of 10%. Every experiment was carried out at a temperature T of 25°C and a strain rate, ė, of ca. 4 x 10 high -3 s high -1 until a finite longitudinal strain of 40% was achieved with a viscosity contrast m of 18.6 between the non-linear viscous layer (Kolb brown plasticine) and the matrix (Beck’s green plasticine with 150 ml oil kg high -1). The apparent viscosity, η, and the stress exponent, n, for the layer at a strain rate ė = ca. 10 high -3 s high -1 and a finite strain e = 10% are 2.23 x 10 high 7 Pa s and n = 5.8 and for the matrix 1.2 x 10 high 6 Pa s and 10.5. These new data that result from incremental analogue modelling corroborate previous suggestions that folds and boudins are coeval structures in cases of plane-strain coaxial deformation with the stiff layer oriented perpendicular to the intermediate Y-axis of the finite strain ellipsoid. They will be of interest for all workers who are dealing with plane-strain boudins and folds, where the fold axes are parallel to the major axis (X) of the finite strain ellipsoid. As has been demonstrated by the first experimental study, coeval folding and boudinage under plane strain, with S ┴ Y, are associated with a significant increase in the thickness of the competent layer. The latter phenomenon does not occur in other cases of simultaneous folding and boudinage, such as bulk pure constriction. To study the impact of layer thickness on the geometry of folds and boudins under pure constriction, we carried out additional experiments using different types of plasticine for a stiff layer and a weaker matrix to model folding and boudinaging under pure constriction, with the initially planar layer oriented parallel to the Xaxis of the finite strain ellipsoid. The stiff layer and matrix consist of non-linear viscous Kolb brown and Beck’s green plasticine, respectively, both of which are strain-rate softening modelling materials. Six runs have been carried out using thicknesses of the stiff layer of 1, 2, 4, 6, 8 and 10 ± 0.2 mm. All experimental runs were carried out at a temperature T of 30 ± 2°C and a strain rate, ė, of ca. 1.1 x 10 high -4 s high -1 until a finite longitudinal strain of 40% was achieved with a viscosity contrast m of 3.1 between the stiff layer (Kolb brown plasticine) and the matrix (Beck’s green plasticine). The apparent viscosity, η, and the stress exponent, n, for the layer at a strain rate ė = ca. 10 high -3 s high -1 and a finite strain e = 10% are 2.23 x 10 high 7 Pa s and n = 5.8 and for the matrix 7.2 x 10 high 6 Pa s and 7.9. Our results suggest a considerable influence of the initial thickness of the stiff layer on the geometry of the deformed stiff layer. There is no evidence for folding in XY=XZ-sections if the initial thickness of the competent layer is larger than ca. 8 mm. If the initial thickness of the competent layer is set at ca. 10 ± 0.2 mm, both folds and boudins develop simultaneously. However, the growth rate of the boudins seems to be higher than the growth rate of the folds. A further expected, but characteristic, aspect of the deformed competent layer is no change in thickness of the competent layer, which can be used to distinguish plane-strain folds and boudins from constrictional folds and boudins. The model results are important for the analysis and interpretation of deformation structures in rheologically stratified rocks undergoing dislocation creep under bulk constriction. Tectonic settings where constrictional folds and boudins may develop simultaneously are stems of salt diapirs, subduction zones or thermal plumes. To make (paleo) viscosimetric statements possible, the rheological data of the different plasticine types were related to the geometrical data. When comparing the normalized dominant wavelength Wd obtained from the deformed layer of the models with the theoretical dominant wavelength (Ld) calculated using the Smith equation (1977, 1979), the latter probably also holds when folding and boudinage develop simultaneously (S ┴ Y) and when boudins develop independently (S ┴ Z), but can obviously not be applied at very low viscosity ratios as is indicated by the low-strain-rate experiments.
To reconstruct ocean circulation changes during specific periods of Earth history, benthic and planktic foraminifera were used as proxies in the different parts of this thesis. Both studied time periods, the Late Cretaceous and the early Pleistocene, are characterized by long-term climate cooling and major changes in ocean circulation. The first part of this thesis concentrated in the Late Cretaceous. During the Late Cretaceous long-term cooling phase, benthic foraminiferal δ18O values show a positive shift lasting about 1.5 Myr (71.5–70 Ma). This shift can be observed on a global scale and has become known as the Campanian-Maastrichtian Boundary Event (CMBE). It is proposed that this δ18O excursion is influenced either by changing intermediate- to deep-water circulation or by temporal build-up of Antarctic ice sheets. Benthic foraminiferal assemblage counts from a southern high-latitudinal site near Antarctica (ODP Site 690) are analyzed to test if the influence of the CMBE on the benthic species composition. One of the two discussed hypotheses for the causation of the δ18O transition is a change in intermediate- to deep-water circulation from low-latitude to high-latitude water masses. This change would result in cooler temperatures, higher oxygen concentration, and possibly lower organic-matter flux at the seafloor, causing a major benthic foraminiferal assemblage change. Another possible explanation of the δ18O transition of the CMBE is significant ice formation on Antarctica. However no major benthic foraminiferal assemblage change would be expected in this case. The benthic foraminiferal assemblage of Site 690 shows a separation of the studied succession into two parts with significantly different species composition. The older part (73.0–70.5 Ma) is dominated by species, which are typical for lower bottom water oxygen concentration and more common in low-latitude assemblages. Species dominating the younger part (70.0–68.0 Ma) are indicators for well-oxygenated bottom waters and more common in high-latitude assemblages. This change in the benthic foraminiferal assemblages is interpreted to represent a shift of low-latitude toward high-latitude dominated intermediateto deep-water sources. A change in oceanic circulation was therefore at least a major component of the CMBE. The Pacific Ocean contributed significantly to the climatic development during the Late Cretaceous cooling period. The contribution of ocean circulation changes in the Pacific Ocean to the Late Cretaceous climatic development in general and the CMBE and Mid-Maastrichtian Event (MME) in particular, however, is poorly understood. Previously measured high resolution planktic and benthic stable isotope data and a neodymium (Nd) isotope record from the Pacific ODP Site 1210 (Shatsky Rise, tropical Pacific Ocean) for the Campanian to Maastrichtian (69.5 to 72.5 Ma) are used to reconstruct changes in surface- and bottom water temperatures as well as changes in the source region of deep- to intermediate waters [see Appendix 4; Jung et al. 2013]. The results of the benthic foraminiferal δ18O and Nd isotope records in combination with Nd isotope records from other studies indicate changes in the intensity of intermediate- to deep ocean circulation in the tropical Pacific across the Campanian-Maastrichtian interval [see Appendix 4; Jung et al. 2013]. During the early Maastrichtian (72.5 to 69.5 Ma), a three-million-year-long period of cooler conditions and a simultaneous change towards less radiogenic Nd isotope signatures is interpreted to represent a period of increased admixture and northward flow of deep waters from the Southern Ocean (Southern Component Water, SCW). This change was probably caused by an intensified formation of deep waters in the Southern Ocean. This was reduced again during the MME (69.5 to 68.5 Ma). This early Maastrichtian cold interval is similar to the CMBEδ13C fall and succeeding δ13C rise towards the MME and is therefore also interpreted to represent tectonically forced, long-term changes in the global carbon cycle and thus a tectonic forcing of the early Maastrichtian climate cooling. Overall, the Campanian-Maastrichtian Nd and stable isotope records of Shatsky Rise indicate changes in ocean circulation that are paralleled by global warming and cooling periods. The fluctuating strength of SCW contribution in the tropical Pacific points towards an increased respectively weakened ocean circulation, which is probably related to the strength of deep-water formation in the Southern Ocean [see Appendix 4; Jung et al. 2013]. For this study, the analysis of benthic foraminiferal assemblages of Site 1210 is carried out for the same time interval (69.5 to 72.5 Ma) as Nd and stable isotopes to evaluate the influence of intermediate- to deep ocean circulation changes on the benthic foraminiferal community. The possible reaction of benthic foraminiferal assemblages is compared to the results of stable isotope and neodymium isotopes. The observed changes in species abundances only partly reflect the circulation changes reconstructed with Nd and stable oxygen istopes. For example, Stensioina spp., Aragonia spp. and Lenticulina spp., cold-water preferring species, start to be increasingly abundant at the beginning of enhanced influence of SCW. However, their abundance pattern does not follow the varying strength of the cold SCW influence at Shatsky Rise. Other species prefer lesser oxygen concentrations and warmer bottom water, e.g. Paralabamina spp. and Globorotalites spp. Paralabamina spp. has its highest relativ abundance at the beginning of the studied succession, where the influence of SCW is small. However, this taxa occurs throughout the record, even though the influence of SCW increases. Globorotalites spp. is even most abundance after the CMBE, where bottom waters are till cold and influenced by SCW. This leads to the conclusion that the varying strength of SCW in the tropical Pacific at Shatsky Rise through the studied interval is not facilitating a significant faunal turnover as has been observed at the South Atlantic Site 690 (Chapter 3). These results of the benthic foraminiferal assemblage analysis suggest a rather minor influence of the SCW on the major environmental factors that are generally influencing benthic foraminiferal communities (e.g., oxygen concentration, organic matter flux to the sea floor, bottom-water temperature). The second major part of this thesis focused on the late Pliocene-earliest Pleistocene. The late Pliocene is characterized by a long-term global cooling trend resulting in a major increase of Arctic ice sheets from around 3 Ma onwards, culminating in the Plio-Pleistocene intensification of the Northern Hemisphere glaciation. At around 2.7 Ma, large amplitude glacial-interglacial excursions (~1‰ δ18O in benthic foraminiferal calcite) in benthic oxygen isotopes can be observed. Marine isotope stage (MIS) 100 at around 2.55 Ma is the first glacial, when widespread ice rafted debris has been found in sediments in the North Atlantic Ocean. To gain a deeper understanding of the climatic evolution of the latest Pliocene-early Pleistocene, it is necessary to improve the reconstructions of North Atlantic paleohydrography, as the North Atlantic provides a key region for global climate. The consequences of the intensification of Northern Hemisphere on the early Pleistocene North Atlantic thermocline stratification and intermediate waters are still poorly understood. However, surface hydrography, the history of the thermocline and development of North Atlantic intermediate waters are well-studied for the Last Glacial Maximum (LGM). These well-known mechanisms responsible for the LGM in comparison with the present-day interglacial North Atlantic are used as an analogue for te early Pleistocene glacialinterglacials cycles. In this study, suborbitally resolved stable oxygen and carbon isotope and Mg/Ca records are measured from a deep-dwelling planktic foraminifera (Globorotaliacrassaformis) from Integrated Ocean Drilling Program Site U1313 (North Atlantic, 41°N) covering marine oxygen isotope stages MIS 103 to 95 (early Pleistocene, 2.6 to 2.4 Ma). The results are interpreted to represent a change in intermediate-water masses on glacialinterglacial timescales. During glacials geochemical records in G. crassaformis (~500–1000 m) bear the imprint of Glacial North Atlantic Intermediate Water (GNAIW), while during interglacials this species reflects the signature of the influence of Mediterranean Outflow Water (MOW) in combination with the subtropical gyre. The comparison of this data with the published records from G. ruber from the same samples facilitates the reconstruction of glacial-interglacial stratification changes of the upper water column at Site U1313. The results show that larger gradients of temperature, salinity and δ13C prevailed during glacials, suggesting a stronger stratification of the upper water column. This can be seen to indicate glacial-interglacial changes in ntermediate water masses in the North Atlantic similar to those reconstructed for the latest Pleistocene. As an additional proxy, the clumped isotope paleothermometer is applied for the Late Cretaceous study as well as for the early Pleistocene. This proxy is commonly assumed to be independent of other factors than temperature. Clumped isotopes are measured for the Late Cretaceous Site 690 on the planktic foraminiferal species Archaeoglobigerina australis and compared to already existing stable oxygen isotopes of this species. This is assumed to enable the reconstruction of paleotemperature independent of ice volume and therefore contribute to the long-lasting discussion whether there was a temporal ice build-up on Antarctic during the Campanian-Maastrichtian cooling period. For the early Pleistocene, the planktic foraminiferal species G. crassaformis is used from Site U1313 from MIS 99 (interglacial) and MIS 98 (glacial). This provides the opportunity to separate ice volume, salinity and temperature effects on the measured δ18O record of G. crassaformis. The results of the clumped isotope measurements reveal comparatively large standard errors. For the Late Cretaceous the standard error of the clumped isotope measurements proved too large to allow any conclusions on the temperature component on the δ18O record of A. australis. For the early Pleistocene, the temperature difference is also too small to be reconstructed with the standard error of the clumped isotope measurements in this study. Measuring many replicates of one sample would minimize the standard error considerably. However, the amount necessary to measure replicates cannot be gained for either time period, as almost all foraminifera were picked from the respective samples. It is concluded that the respective questions may be solved with a different method of clumped isotope analysis requiring less sample material. This method is, for example, available at the ETH Zurich.
Trace elemental concentrations of bivalve shells content a wealthy of environmental and climatic information of the past, and therefore the studies of trace elemental distributions in bivalve shells gained increasing interest lately. However, after more than half century of research, most of the trace elemental variations are still not well understood and trace elemental proxies are far from being routinely applicable. This dissertation focuses on a better understanding of the trace elemental chemistry of Arctica islandica shells from Iceland, and paving the way for the application of the trace elemental proxies to reconstruct the environmental and climatic changes. Traits of trace elemental concentrations on A. islandica shells were explored and evaluated. Then based the geochemical traits of the shells, four non-environmental/climatic controlling is indentified. (1) Trace elemental concentrations of bivalve shells are effected by early diagenesis by the leach or exchange of elemental ions, especially in shell tip part, even with the protection of periostrucum; (2) The analytical methods also affect the results of trace elemental concentrations, especially for the element, such as Mg, which is highly enriched in organic matrices; (3) Shell organic matrices are found play a dominating role on the concentration of trace elements on A. islandica shells. Most trace elements only occurred in insoluble organic matrices (IOM), although others are only found in the carbonate fraction. IOM of A. islandica shells is significantly enriched in Mg, while Li and Na are more deplete in IOM, but enriched in shell carbonate. Ba is more or less even contented in IOM and shell carbonate. The concentrations of certain elements vary between primary layer and secondary layer; (4) The vital /physiological controlling on trace elemental distributions of bivalve shells is also confirmed. Six elemental (B, Na, Mg, Mn, Sr, and Ba) concentrations show significant correlation (exponential functions) with ontogenetic age and shell grow rates (logarithmic equations). It is worthy to remark that B, Mg, Sr and Ba concentrations are negatively correlated with shell growth rate, positive with ontogenetic age, while the concentrations of Na and Mn show the opposite trends. At last, all the controlling described above can be taken into account and corrected to extract the environmental and climatic signal by a kind of standardization. The derived six exponential functions of the high correlations between six trace elemental concentrations and ontogenetic year are applied to make the standardization of these element-Ca ratios. The gotten standardized indices are compared with the variations of environmental and climatic parameters in this region, and many correlations are found. Standardized indices of Sr/Ca ratios are strongly related to the sun spot number, autumn NAO, autumn Europe surface air temperature (SAT) and Arctic sea surface temperature anomaly (TA), and those of Mg/Ca ratios are strongly associated with Arctic TA, Europe SAT and Solar variation (irradiance). The variations of autumn Europe SAT demonstrated more similarity with standardized indices of B/Ca than other parameters. Except for the SAT index of Arctic, the standardized indices of Na/Ca showed no distinct relation to temperature. European precipitation and the Arctic sea level pressure index compared well the Na/Ca ratios of the shells, and so did the autumn NAO. Standardized indices of Mn/Ca were correlated with the number of hurricanes in the North Atlantic, Northern Europe SAT and sun spot number.
This work analyses several granitic bodies of the Variscan Orogen of Central and Western Europe in order to improve our knowledge about different aspects of their evolution, regarding their ascent and emplacement mechanisms, as well as their deformation history. In the Iberian Massif two granitoid bodies, namely the La Bazana pluton and the Nisa-Alburquerque batholith, were studied in order to decipher their ascent and emplacement history. The La Bazana pluton is a small, sub-circular body in map view that intruded into rocks of the Ossa-Morena Zone in the core of a late upright antiform. Its three-dimensional drop-pipe shape, its internal dome foliation pattern and the structure of the host rock suggest that the magma ascended and emplaced diapirically. The Nisa-Alburquerque batholith is a large body that intruded into rocks of the Central Iberian Zone, the Central Unit, and the Ossa-Morena Zone. Its cartographic shape is elongate and parallel to the NW—SE to WNW—ESE Variscan structures. In the light of the available structural data and the gravimetric models, the intrusion is viewed as a continuous lateral magma flow from the eastern root guided towards the west through the southern limb of a kilometre-scale antiform. As mass-transfer mechanisms, a combination of rigid translation of the country rocks, stoping, and possibly ballooning is proposed. In the Bohemian Massif several small granitoid bodies showing a strong solid-state deformation were studied in order to integrate their tectonometamorphic history in the geotectonic framework of the south-western Bohemian Massif, focusing principally on the deformation phase referred to as D3. Four ductile deformation phases are proposed for the study area. D1 produced high-temperature fabrics under upper amphibolite to granulite facies conditions. Its kinematics is unknown. D2 occurred under amphibolite to upper greenschist facies conditions under N—S to NNW—SSE compression. It is responsible for a subvertical NW—SE striking foliation in migmatites developed under dextral simple shear and for the deformation at the Bayerischer Pfahl shear-zone system at its earlier stages. Many granitoid dykes and stocks were found to be affected by sinistral shear along subvertical planes trending ENE to ESE. Since this deformation, which is called D3 in the present work, is not compatible with a N—S to NNW—SSE compression, it is proposed that these sinistral shear zones in granites do not belong to the Bayerischer Pfahl shear-zone system and constitute themselves a separated one, which is called “D3 shear-zone system”. D3 took place under upper greenschist to lower amphibolite facies conditions (~480-550°C). Both the intrusion and the deformation of the granites affected by D3 occurred at deep to intermediate levels of the crust, whereas the deformation took place under NE—SW compression. Datings on two of the deformed granites yielded 324.4 ± 0.8 Ma and 315.0 ± 1.0 Ma: Thus, the age of D3 is most probably ~315 Ma. The intrusion of most of the sheared granitoids was pre-kinematic with respect to D3. After D3 the N—S to NNW—SSE compression which governed D2 was restored, giving way to the next deformation phase D4, which was linked to further deformation at and next to the principal shears of the Bayerischer Pfahl shear-zone system under greenschist facies conditions. The causes for the change of the stress field leading to a NE—SW compression during D3 might be related to (1) global changes in the dynamics of the tectonic plates in late Variscan times, (2) orogenic collapse leading to the sinking of the Teplá-Barrandian and lateral extrusion of the surrounding Moldanubian rocks, (3) distortion of the regional stress field by local intrusion of large stocks, such as the Saldenburg granite of the Fürstenstein Massif, or (4) distortion of the regional stress field due to the existence of ephemeral releasing bends in the Bayerischer Pfahl shear zone during its early evolution.
Gridded maps of meteorological variables are needed for the evaluation of weather and climate models and for climate change monitoring. In order to produce them, values at locations where no observing stations are available need to be estimated from point-wise observations. For the interpolation of meteorological observations deterministic and stochastic methods are often combined. Deterministic methods can account for ancillary information such as elevation, continentality or satellite observations. Stochastic methods such as kriging reproduce observed values at the station locations and also account for spatial variability. In the first two studies of this thesis, a flexible interpolation method for the gridding of locally observed daily extreme temperatures is developed that also provides an optimal estimate of the interpolation ncertainty. In the third study, an observational dataset is created using this interpolation method and then applied to evaluate a climate simulation for Africa.
In the first study, the Regression-Kriging-Kriging (RKK) method is tested for the interpolation of daily minimum and maximum temperatures (Tmin and Tmax) in different regions in Europe. RKK accounts for elevation, continentality index and zonal mean temperature and is applicable in regions of differing station density and climate. The accuracy of RKK is compared to Inverse Distance Weighting, a common deterministic interpolation method, and to Ordinary Kriging, a common stochastic interpolation method. The first step in RKK is to use regression kriging, in which multiple linear regression accounts for topographical effects on the temperature field and kriging minimizes the regression error, to interpolate climatological means. In the second step daily deviations from the monthly climatology are interpolated using simple kriging. Owing to the large climatological differences across the investigation area the interpolation is performed in homogeneous subregions defined according to the Köppen-Geiger climate classification. Cross validation demonstrates the superiority of RKK over the simpler algorithms in terms of accuracy and preservation of spatial variability. The interpolation performance however strongly varies across Europe, being considerably higher over Central Europe (highest station density) than over Greenland (few stations along the coast line). This illustrates the strong impact of the station density on the accuracy of the interpolation result. Satellites provide comprehensive observations of climate variables such as land surface temperature (LST) and cloud cover (CC). However, LST is associated with high uncertainty (standard error ~ 1-2°C), preventing its direct application in meteorology and climatology. The second study investigates the usefulness of LST and CC as predictors for the gridding of daily Tmin and Tmax. The RKK algorithm is compared with similar interpolation methods that apply LST and CC in addition to the predictors used with the RKK algorithm. The investigation is conducted in two regions, Central Europe and the Iberian Peninsula, which differ strongly in average cloud cover (Central Europe is approximately 30% cloud free and the Iberian Peninsula approximately 60 % cloud free). RKKLST (in which monthly mean LST is used as an additional predictor) yields for Central Europe no clear improvement over RKK, yet it reduces the interpolation error over the Iberian Peninsula. This finding can be explained by the higher percentage of cloud free pixels over that region in summer which enables a more robust determination of monthly mean LST. Adding a regression step for daily anomalies (using the predictor CC) yields the RKRK method and improves the preservation of spatial variability over the Iberian Peninsula. Moreover, a successive reduction of the station number (from 140 to 10 stations) reveals an increasing superiority of RKKLST and RKRK over RKK in both regions.
The application of a gridded observational dataset for climate monitoring or climate model validation requires knowledge of the uncertainties associated with the dataset. The estimation of the interpolation uncertainty, here the inter quartile range is the used uncertainty measure, is therefore an important issue within the frame of this thesis. By means of cross validation it is shown that the largest uncertainties occur in regions of low station density (e.g. Greenland), in mountainous regions and along coastlines (in these regions model evaluation results should be interpreted carefully). The magnitude of the interpolation error mainly depends on the station density, while the complexity of terrain has substantially less influence. On average over all regions and investigation days the target precision of the uncertainty estimate is reached. However, on local scales and for single days it can be clearly over- or underestimated. The application of satellite-derived predictors (LST and CC) yields no noteworthy improvement of the uncertainty estimate.
In the last study two regional climate simulations for Africa using the ERA-Interim driven COSMO-CLM (CCLM) model at two different horizontal resolutions (0.22° and 0.44°) are validated. It is assessed whether observed patterns and statistical properties of daily Tmin and Tmax are correctly represented in the model. The ERA-Interim reanalysis and a specially created observational dataset are used as reference. The observational dataset is generated by applying the RKRK algorithm (developed within the second study). The investigations show an occasionally large bias in Tmin and Tmax. The hemispheric summers are generally too warm and the temporal variability in temperature is too high, particularly over extra tropical Africa. The diurnal temperature range is overestimated by about 2°C in the northern subtropics but underestimated by about 2°C over large parts of the African tropics. CCLM reproduces the observed frequency distribution of daily Tmin and Tmax in all African climate regions, and the extreme values in the lower percentiles (5, 10, 20%) for Tmin are well simulated. The higher percentiles (80, 90, 95%) for Tmax are however overestimated by 2-5°C. For both Tmin and Tmax the 0.22° simulation is on average 0.5°C warmer than the 0.44° simulation. Additionally, the higher percentiles are about 1°C warmer for both Tmin and Tmax in the higher resolution run, while the lower percentiles in both runs match very well. Although the temperature pattern is represented in more detail along the coastlines and in topographically complex regions, the higher resolution simulation yields no qualitative improvement.
To summarize, the choice of the appropriate algorithm mainly depends on the interpolation conditions. In cases where the station density is high across the target region and the predictor space is adequately covered by observing stations, the computationally less demanding RK algorithm should be preferred. In regions where the station density is low the more robust RKRK algorithm should be the first choice. Due to the strong physical relation of both CC and LST to Tmin and Tmax the missing information is at least partially compensated for. The estimation of the interpolation uncertainty could be improved by applying a normal score transformation to the data prior to a kriging step. This is because the kriging assumption that the increments of the variable of interest are second order stationary can be approximately met by a normal score transformation.
High field strength element systematics and Lu-Hf & Sm-Nd garnet geochronology of orogenic eclogites
(2008)
Concerning the Bulk Silicate Earth (BSE), the depleted mantle and the continental crust are thought to balance the budget of refractory and lithophile elements, resulting in complementary trace element patterns. However, the two high field strength elements (HFSE) Niob and Tantal appear to contradict this mass balance. All reservoirs of the silicate Earth exhibit subchondritic Nb/Ta ratios, possibly as a result of Nb depletion. The two HFSE Zr and Hf on the other hand seem not to be fractionated between the silicate reservoirs. They show more or less chondritic Zr/Hf ratios. In this study a series of orogenic eclogites from different localities was analyzed to determine their HFSE concentrations and to contribute to the question if eclogites could form a hidden reservoir to account for the mass imbalance of the BSE. The results show that the orogenic eclogites have subchondritic Nb/Ta ratios and near chondritic Zr/Hf ratios. The investigated eclogites show no fractionation of Nb/Ta ratios and no enrichment of Nb compared to e.g. MOR-basalts, the likely precursor of these rocks. With an average Nb/Ta ratio of 14.9 these eclogites could not balance the differences between BSE and chondrite. Additionally, with an average Nb/Ta ≈ MORB they also cannot balance the small differences in the Nb/Ta of the crust and the mantle. LA-ICPMS analyses of rutiles in these eclogites reveal a zonation of Nb/Ta ratios in this mineral, with rutile cores having higher Nb/Ta than rutile rims. As a consequence, Laser Ablation data of rutiles have to be evaluated carefully and cannot necessarily reflect a bulk rock Nb and Ta composition, although over 90% of these elements reside in rutile.
In Belize, which is well known for the Belize Barrier Reef and its offshore atolls, coastal lagoons are frequent morphological features along the coast. They represent transitional environments between siliciclastic and carbonate settings. In order to shed light into the Holocene evolution of coastal lagoon environments, five localities along the central coast of Belize were selected as coring sites. These include four coastal lagoons and one marsh area, namely Mantatee Lagoon, Mullins River Beach, Colson Point Lagoon, Commerce Bight Lagoon and Sapodilla Lagoon. A total of 26 sediment cores with core lengths ranging from 109 cm to 500 cm, were drilled using a Lanesky-vibracorer. Overall, 73 m of Holocene sediments and Pleistocene soil were recovered. Together with 58 radiocarbon dates the sediments reveal details on the sediment architecture and depositional features of the localities.
Marine inundation of the mainland and coastal lagoon formation started around 6 kyrs cal BP.
As a response to sea-level rise during the Holocene transgression, facies retrograded towards the coast, as seen in marginal marine overlying brackish mollusk faunas. Evidence for late Holocene progradation of facies due to sea-level stagnation is largely lacking. The occurrence of landward thinning sand beds, hiatuses and marine fauna in lagoonal successions are indications of event (overwash) sedimentation. Sediments recovered are largely of Holocene age (<7980 cal BP), overlying Pleistocene sections. Analyses of sediment composition and texture, radiocarbon dating and mollusk shell identification were used to describe and correlate sedimentary facies.
XRD analyses have identified quartz as the dominant mineral, with the Maya Mountains as main source of coastal lagoon sediments. The most common sedimentary facies include peat and peaty sediment, mud, sand, and poorly sorted sediments. Pleistocene soil forms the basement of Holocene sediments. Holocene mud represents lagoon background permanent sedimentation.
Peats and peat-rich sequences were deposited in mangrove swamp environments, whereas sandy facies mainly occur in the shoreface, beach, barriers, bars, barrier spits and overwash deposits. Facies successions could be identified for each locality, but it has proven difficult to correlate the stratigraphic sequences, especially among localities. These differences among the five locations studied suggest that apart from regional influence such as sea-level rise, local environmental factors such as small-scale variation in geomorphology and resulting facies heterogeneity, connectivity of the lagoon with the sea, antecedent topography and river discharge, were responsible for coastal sedimentation and lagoon development in the Holocene of Belize.
Faunal composition and distribution patterns of mollusk assemblages from 20 shell concentrations in cores collected in coastal lagoons, a mangrove-fringed tidal inlet and the marginal marine area (shallow subtidal) show considerable variation due to environmental heterogeneity and the interplay of several environmental factors in the course of the mid-late Holocene (ca. 6000 cal BP to modern). The investigated fauna ≥2 mm comprises 2246 bivalve, gastropod and 11 scaphopod specimens. Fifty-three mollusk species, belonging to 42 families, were identified. The bivalve Anomalocardia cuneimeris and cerithid gastropods are the dominant species and account for 78% of the total fauna. Diversity indices are low in concentrations from lagoons and relatively high in the marginal marine and tidal inlet areas.
Based on cluster analysis and nonmetric multidimensional scaling (NMDS), seven lagoonal assemblages and three marginal marine/tidal inlet assemblages were defined. A separation between lagoonal and marginal marine/tidal inlet assemblages seen in ordination indicates a lagoon-onshore gradient. The statistical separation among lagoonal assemblages demonstrates environmental changes during the Holocene evolution of the coastal lagoons, which is probably related to the formation of barriers and spits. The controlling factors of species distribution patterns are difficult to figure out, probably due to the heterogeneity of the barrier-lagoon systems and the interaction of paleoecological and paleoenvironmental factors. In addition to the taxonomic analysis, a taphonomic analysis of 1827 valves of A. cuneimeris from coastal lagoons was carried out. There is no relationship between depth and age of shells and their taphonomic condition. Size-frequency distributions and right-left valve ratios of A. cuneimeris suggest that valves were not transported over long distances but were deposited parautochthonously in their original habitat. Shells from tidal inlet and marginal marine environments were also predominantly deposited in their original habitats.
Since the Belize coast was repeatedly affected by hurricanes and the paleohurricane record for this region is poor, the sediment cores have been examined in order to identify storm deposits.
The paleohurricane record presented in this study spans the past 8000 years and exhibits three periods with increased evidences of hurricane strikes occurring at 6000-4900 cal yr BP, 4200-3600 cal yr BP and 2200-1500 cal yr BP. Two earlier events around 7100 and 7900 cal yr BP and more recent events around 180 cal yr BP and during modern times have been detected. Sand layers, redeposited corals and lagoon shell concentrations have been used as proxies for storm deposition. Additionally, hiatuses and reversed ages may indicate storm influence. While sand layers and corals represent overwash deposits, the lagoon shell concentrations, which mainly comprise the bivalve Anomalocardia cuneimeris and cerithid gastropods, have been deposited due to changes in lagoon salinity during and after storm landfalls. Comparison with other studies reveals similarities with one record from Belize, but hardly any matches with other published records. The potential for paleotempestology reconstructions of the barrier-lagoon complexes along the central Belize coast differs depending on geomorphology, and deposition of washovers in the lagoon basins is limited, probably due to the interplay of biological, geological and geomorphological processes.
This study describes the Holocene sedimentary lagoonal deposition history, including event sedimentation and benthic foraminiferal analyzes, from about 10 kyrs BP until today. This is the first study describing the sedimentation of a Maldivian atoll lagoon in such detail. Thirty-nine sediment cores have been recovered from the deep Rasdhoo Atoll lagoon of the Maldives (4°N/73°W). Seventeen sediment cores were opened, described, and 296 sediment samples have been collected and analyzed. Different methods have been used to evaluate the coarse- and fine-grained carbonate components and a total of fifty-eight samples have been dated radiometrically by Beta Analytic Inc., Miami, Florida. In general, the Rasdhoo Atoll lagoon sediments can be divided into (1) a Late Pleistocene soil, (2) an early Holocene peat layer composed of mangrove deposits which mark the beginning inundation of the atoll lagoon by the rising Holocene sea-level at 10,320 ± 100 yrs BP, and (3) carbonate sediments starting to fill up the lagoon 7850 ± 140 yrs BP until today. The transition from peat to carbonate is characterized by a considerable hiatus. Six different carbonate sediment facies are classified by statistical analyses, listed in decreasing abundance:
(1) mollusk-coral-algal floatstone to rudstone (30%)
(2) mollusk-coral-red algae rudstone (23%)
(3) mollusk-coral-algal wackestone to floatstone (23%)
(4) mollusk-coral wackestone (13%)
(5) mollusk-coral mudstone to wackestone (9%)
(6) mollusk mudstone (2%)
Based on grain-sizes in combination with coral identification, the facies represent both lagoonal background sedimentation (mostly fine-grained sediments (matrix >50%)) and event sedimentation (coarse-grained sediment layers composing reefal components).
Six coarser grained layers in muddy background sediments of the Rasdhoo Atoll lagoon were interpreted as Holocene tsunami events, based on the increase of allochthonous skeletal material with shallow-water reef affinity such as fragments of shallow-water coral species, coralline red algae, and reef-dwelling foraminifera in these layers, as well as AMS dating:
• Event 1: 420 - 890 yrs BP (655 yrs BP)
• Event 2: 890 - 1560 yrs BP (1225 yrs BP)
• Event 3: 2040 - 2340 yrs BP (2190 yrs BP)
• Event 4: 2420 - 3380 yrs BP (2900 yrs BP)
• Event 5: 3890 - 4330 yrs BP (4110 yrs BP)
• Event 6: 5480 - 5760 yrs BP (5620 yrs BP)
Five of the six layers may be correlated to previously published tsunami events at adjacent coastal research sites. The mid-late Holocene atoll lagoon archive is incomplete though based on the assumption that major earthquakes at the Indonesian subduction zone generated more than six major tsunamis during the past 6.5 kyrs.
According to Gischler (2006), the sediments of the Rasdhoo Atoll lagoon can be divided into two areas: (1) a central to marginal deep lagoon with a lateral west-to-east gradient of sediment facies distribution, visible in sections <4 kyrs BP with sedimentary facies of mudstone to wackestone in the western part (e.g., cores 16, 18, and 34) and coarse-grained coral and algal-rich sediments in the eastern part of the lagoon (e.g., cores 30 and 31). (2) A northern enclosed and shallow area between the sand apron and the sand spit accumulating “sandy” sediments of wackestone facies (cores 2, 19, 25, and 26).
Comparing the sediment accumulation data of the lagoon with two reconstructed local sea-level curves, three different sequence-stratigraphical systems tracts are visible: (1) a lowstand systems tract (LST) >10 kyrs BP. Pleistocene brownish soil superposing subaerially exposed Pleistocene reef limestone. (2) A transgressive systems tract (TST) 10-6.5 kyrs BP. A peat layer marks the beginning of the inundation, and the carbonate sedimentation starts with very low sedimentation rates of 0.02 m/kyr. (3) A highstand systems tract (HST) 6.5-0 kyrs BP, further divided into three stages (6.5-3, 3-1, 1-0 kyrs BP). The sea-level rise slowed down, sedimentation rates are increasing continuously up to a maximum of 1.4 m/kyr, the sand spit developed some 4 kyrs BP, the lagoonal circulation got restricted, and the lateral west-to-east gradient of grain-size accumulation started. From 1-0 kyrs BP the sedimentation rates slowed down to modern mean sedimentation rates of 0.6 m/kyr.
Two cores, one core from the center of the lagoon (core 16) and one core from the northern margin of the lagoon (core 19), have been analyzed on diversity and assemblages of benthic foraminifera in high-resolution. The transitions of Ammonia spp. to a more even and diverse fauna marks a significant environmental change at 7.0 kyrs BP in core 16 (onset of a stable environment in the deep lagoon after the sea-level rise slowed down at HST stage 1) and at 4.0 kyrs BP in core 19. A continuing environmental change after 1.4 kyrs BP in core 16 caused the fauna to become more even, a recovery of diversity and a permanent decline of foraminiferal accumulation rate. The changes in the faunas at 4.0 kyrs BP and at 1.4 kyrs BP could be explained with the sand spit formation in the northwestern and western lagoon. The sand spit has apparently acted as an obstacle in lagoonal circulation and might have caused unstable environmental conditions due to a more rapid circulation at the shallow marine site of core 19 and a slowdown of bottom water circulation in the main lagoon (core 16) leading to higher residence times and to lower oxygen and higher nutrient concentrations.
This research was conducted in the Rwenzori Region of the Western Branch, East African Rift System (EARS). The EARS is a tectonic structure extending over a length of more than 3000 km from the Afar Triple Junction, in Ethiopia, to Lake Malawi in the south. The Western Rift System is a roughly NE to ENE trending sector of the EARS, which runs along the western boundary of Uganda and the neighboring Democratic Republic of Congo (D.R.C). It stretches 2100 km from Nimule, NW on Uganda-Sudan border, extending to Lake Malawi in the SE of Africa. The unusual uplift of the Rwenzori Mountains within an extensional regime and the mechanisms associated with the high frequency of seismic activity in the region was hardly understood and therefore, had remained a subject of contention that needed to be critically addressed in detail. To my knowledge, this was probably the first study to be performed and documented in great depth within the domains of seismic noise variation, seismic anisotropy and b value analyses beneath the Rwenzori Region. After about six years of operation (2006-2012), the seismology group of the RIFTLINK Research Project (www.riftlink.org) acquired a vast amount of high-quality, digital data that were collected using a seismic network of well calibrated seismic equipment. The project was divided into two phases. Phase I, that operated between February 2006 - September 2007, consisted of thirty-two temporary seismic stations, which were selectively spread out in the Rwenzori Region on the Ugandan side, to detect and record extremely weak as well as strong naturally occurring earthquakes. The seismic equipment used included EDL and REFTEK digitizers, which were coupled with Güralp and MARK sensors respectively (REFTEKS: only short-period MARK sensors, EDLs: short-period MARK plus few broadband Güralp Sensors). Exactly 22375 earthquakes were recorded. The data were processed using the SEISAN software package. About 14413 earthquakes were carefully localized using the velocity model of Bram (1975) that implements a Vp=Vs ratio fixed at 1.74. Phase II, that extended between 2009-2012 consisted of thirty-two seismic stations, which were spread out around the Rwenzori Mountains, both on the Ugandan side and the neighboring D.R.C. Only Taurus digitizers that were coupled with Trillium sensors were used in the D.R.C. On the Ugandan side however, both EDL and Taurus digitizers, which were coupled with Trillium and Güralp sensors were used. ...
The present PhD-thesis was prepared within subproject B8 of the DFG-Sonderforschungsbereich (SFB) 641 “The Tropospheric Ice Phase”. The subproject B8 was entitled “Interactions of volatile organic compounds with airborne ice crystals”. Results of previous studies have shown that various volatile organic compounds (VOC) and semivolatile organic compounds (SVOC) are incorporated into the atmospheric ice phase and several uptake mechanisms are discussed in the literature. The aim of this study was to identify the dominating VOC and SVOC in airborne snow collected at Jungfraujoch in the Swiss Alps (3580 m asl) and to study in laboratory experiments the uptake mechanism of organic compounds into snow and ice. For this purpose an analytical method to analyse freshly fallen snow samples was developed and evaluated in a first step. The method consists of headspace (HS) solid phase dynamic extraction (SPDE) followed by gas chromatography combined with mass spectrometry (GC/MS). During the extraction process a new cooling device was successfully integrated into the HS-SPDE-GC/MS method to enhance the extraction yield. Extraction and desorption parameters such as the number of extraction cycles, extraction temperature, desorption volume and desorption flow rate have been optimized. Detection limits for benzene, toluene, ethylbenzene, m-, p-, o- xylene (BTEX) ranged from 19 ng L-1 (benzene) to 30 ng L-1 (m/p-xylene), while those for C6-C10 n-aldehydes ranged from 21 ng L-1 (n-heptanal) to 63 ng L-1 (n-hexanal). Furthermore, freshly fallen snow samples were collected at the High Altitude Research Station Jungfraujoch (3580 m asl, Switzerland) during the field campaigns “Cloud and aerosol characterization experiment” (CLACE) 4 and 5 in February and March 2005 and 2006, respectively. Freshly fallen snow samples collected directly in-cloud on a high altitude remote location were used as approximation of airborne ice crystals since sampling of airborne ice crystals in quantities sufficient for analysis of individual organic compounds is not yet possible. In the collected snow samples a wide range of organic compounds were identified, namely BTEX, n-aldehydes (C6-C10), terpenes, chlorinated hydrocarbons and alkylated monoaromatics. The most abundant organic compounds in snow samples from Jungfaujoch during CLACE 4 and 5 were n-hexanal with a median concentration of 1.324 μg L-1 (CLACE 5) followed by n-nonanal (CLACE 5) with a median concentration of 1.239 μg L-1. High concentration variations of the analytes in snow samples collected at the same time at the same place argue for a heterogeneous composition of snow and ice. Several indicators were found that the origin of the n-aldehydes in the snow can be attributed to direct biogenic emissions from vegetation and indirect biogenic emissions through photochemical oxidation of fatty acids and alkenes. In a second step laboratory experiments were carried out to clarify the uptake mechanism of volatile and semivolatile organic compounds into snow/ice. Organic compounds can be incorporated into the atmospheric ice phase either by the process of gas scavenging, liquid scavenging (riming) or particle scavenging. Gas scavenging (incorporation of the organic compounds from the gas phase during growing of ice crystals) revealed to be ineffective based on previous laboratory experiments in which ice crystals were growing in the presence of aromatic hydrocarbons (BTEX) in the gas phase. In the present study the process of liquid scavenging (riming) was investigated in the laboratory using aqueous standard solutions containing BTEX, naldehydes (C6-C10), methyl tert-butyl ether (MTBE) and ethyl tert-butyl ether (ETBE). The headspace above the standard solution was sampled after adjusting the aqueous solutions to definite temperatures by use of a thermostat. Measurement were carried out at 25°C, 15°C and 5°C (water), -5°C and -15°C (supercooled water) and -25°C (ice). Results have shown that the known trend of lower gas phase concentrations over water concomitant with lower temperatures (Henry’s Law) is only valid for temperatures above 0°C. At temperature below 0°C, increasing concentrations of the analytes (BTEX, MTBE, ETBE and n-aldehydes) were determined in the gas phase together with decreasing temperatures. Dimensionless Henry’s law coefficients (KAW) were calculated from the concentrations of the organic compounds in the headspace above the standard solutions at temperatures between 25°C and -25°C. The observed inversion of Henry’s law coefficients of volatile and semivolatile organic compounds at a water temperature of approximately 0°C is explained by the formation of ordered zones of H2O molecules in supercooled water called “ice-like-clusters”. Together with decreasing temperatures the degree of formation of ordered zones increases which results in the removal of the organic molecules from the liquid phase and transfer into the gas phase. At a temperature of -25°C the supercooled water is converted into ice and a further significant increase of the gas phase concentrations of hydrophobic compounds such as BTEX is observed. In comparison, less hydrophobic compounds such as MTBE, ETBE and n-aldehydes are detected in lower amounts in the gas phase above the water/ice phase due to the higher water solubility and lower Henry coefficients compared to BTEX. The results show that in the absence of particles the uptake of BTEX MTBE, ETBE and C6-C10-naldehydes into ice not enhanced during freezing of a supercooled liquid, since at -25°C for these analytes the concentrations in the gas phase are higher at -25°C (ice) compared with -15°C (supercooled liquid). The heterogeneous distribution of BTEX and n-aldehydes concentrations in snow samples collected during the CLACE field campaigns suggests that adsorption of the organic compounds to particles followed by incorporation of the particles into the snow and ice might play a major role in the uptake process of organic compounds into snow and ice. To increase the knowledge about uptake processes of organic compounds into snow and ice further experiments are required with should include aerosol particles in the experimental setup to evaluate the influence of particle scavenging in the uptake processes.
In der hier vorliegenden Arbeit wurde der troposphärische Kreislauf von Carbonylsulfid (COS) untersucht. COS ist ein Quellgas des stratosphärischen SulfatAerosols, das die Strahlungsbilanz beeinflussen und den chemischen Abbau des stratosphärischen Ozons beschleunigen kann. Trotz zahlreicher Studien sind die Quellen und Senken des atmosphärischen COS bisher nur unzulänglich quantifiziert. Insbesondere bestehen große Unsicherheiten in den Abschätzungen der Beiträge des Ozeans und der anthropogenen Quellen, sowie der Senkenstärke der Landvegetation. Schiffs und flugzeuggetragene Messungen des atmosphärischen COS ergaben kein einheitliches interhemisphärisches Verhältnis (IHR=MNH /M SH ). Während die Messungen von Bingemer et al. (1990), StaubesDiederich (1992) und Johnson et al. (1993) ein IHR zwischen 1.10 und 1.25 zeigten, fanden die Messungen von Torres et al. (1980), StaubesDiederich (1992), Weiss et al. (1995) und Thornton et al. (1996) keinen oder nur einen geringfügigen N/SGradienten. Die Untersuchung von Chin und Davis (1993) zeigt ein N/SVerhältnis der COS Quellstärke von 2.3, das hauptsächlich auf die stärkeren anthropogenen Quellen auf der Nordhalbkugel zurückzuführen ist. Es ist unklar, ob der zeitweilige Konzentrationsüberschuß der Nordhemisphäre Zeichen anthropogener Quellen dort oder Teil eines durch die Senkenfunktion der Landpflanzen verursachten saisonalen Signals ist. Die Konsistenz der Breitenverteilung des COSMischungsverhältnisses mit den geographischen bzw. saisonalen Variationen der COSQuellen und Senken muß überprüft werden. Dazu werden genaue Kenntnissen der Quell und Senkenstärken des atmosphärischen COS und ihrer raumzeitlichen Variabilität benötigt. Vor dem obigen Hintergrund ergeben sich als Schwerpunkte dieser Arbeit: (1) der Austausch von COS zwischen Atmosphäre und Ozean sowie (2) zwischen Atmosphäre und terrestrischer Vegetation und (3) die raumzeitliche Variabilität des atmosphärischen COS. Zur Untersuchung des Austausches von COS zwischen Atmosphäre und Ozean wurde das KonzentrationsUngleichgewicht von COS zwischen Ozean und Atmosphäre durch Messungen des COS im Seewasser und in der Meeresluft ermittelt und die resultierenden Austauschflüsse mit einem Modell berechnet. Die Messungen fanden an Bord des Forschungsschiffs Polarstern während der Fahrten ANT/XV1 (15.10.6.11.1997, BremerhavenKapstadt) und ANT /XV5 (26.5.6.20.1998, KapstadtBremerhaven) statt. Die Konzentration des gelösten COS und das Sättigungsverhältnis von COS zwischen Ozean und Atmosphäre zeigen ausgeprägte Tagesgänge und saisonale und geographische Variationen. Die mittlere Konzentration von COS im Seewasser beträgt 14.7 pmol L -1 für die HerbstFahrt bzw. 18.1 pmol L -1 für die SommerFahrt. Höchste COSKonzentrationen werden in der jeweiligen SommerHemisphäre und in Gebieten mit hoher biologischer Produktivität beobachtet, d.h. im BenguelaStrom im November, im NordostAtlantik im Juni und in den Auftriebgebieten vor Westafrika im Oktober bzw. Juni. In den übrigen Gebieten sind die Konzentrationen um eine Größenordnung niedriger. Die Konzentration von COS im Seewasser steigt frühmorgens von ihrem tiefsten Stand an. Um ca. 15 Uhr Ortszeit erreicht sie ihr Maximum, danach nimmt sie ab. Der Tagesgang unterstützt die Theorie, daß COS im Seewasser photochemisch produziert wird. Während der Tagesstunden wird eine Übersättigung des offenen Ozean für COS gefunden. Dagegen ist eine Untersättigung des Ozeans in den späten Nachtstunden zu beobachten. Der Ozean wirkt in den Tagesstunden als COSQuelle, in der späten Nacht als COSSenke. Die Untersättigung tritt sogar im Sommer in produktiven Meeresgebieten regelmäßig auf. Eine Konsequenz dieser Beobachtung ist die weitere Reduzierung der ozeanischen Quelle von COS gegenüber bisher publizierten Abschätzungen. Methylmercaptan (CH 3 SH) ist in allen Seewasserproben zu beobachten. Der Tagesmittelwert der CH 3 SHKonzentration variiert zwischen 29 und 303 pm L -1 und ist 316 fach größer als der der COSKonzentration. Der Tagesgang der CH 3 SHKonzentration zeigt ein Minimum um die Mittagszeit. Die Tagesmittel der CH 3 SH und COSKonzentrationen sind signifikant miteinander korreliert. Diese Daten liefern den Beweis dafür, daß CH 3 SH eine der wichtigen Vorgängersubstanzen von COS ist. Die Regressionslinie der Korrelation zwischen den mittleren COS und CH 3 SHKonzentrationen weist nur einen geringfügigen Achsenabschnitt auf. Somit kann die CH 3 SHKonzentration als ein Indikator der Konzentration von COSVorgängern benutzt werden. Es besteht außerdem eine Korrelation zwischen der CH 3 SHKonzentration und dem Logarithmus der Konzentration des gelösten Chlorophyll a. Diese Korrelation deutet darauf hin, daß der Gehalt von CH 3 SH im Seewasser eine enge Beziehung zur marinen Primärproduktion hat. COS wird im Seewasser durch Hydrolyse abgebaut. Die Abbaurate hängt von der Temperatur des Seewassers ab. Je wärmer das Seewasser ist, desto schneller wird COS abgebaut, und um so kürzer ist die Lebenszeit von COS im Seewasser. Die Lebenszeit kann einerseits durch das ReaktionsgeschwindigkeitsGesetz von Arrhenius berechnet werden, andererseits läßt sie sich durch exponentielle Anpassung an den nächtlichen Konzentrationsverlauf (d.h. bei Abwesenheit von Photoproduktion) abschätzen. Eine solche Anpassung des exponentiellen Abklingens wurde anhand von dicht gestaffelten Messungen während einiger Nächte vorgenommen. Die gefitteten Lebenszeiten stimmen mit den theoretischen Werten gut überein, obwohl die gefittete Lebenszeit neben Hydrolyse noch von anderen Prozessen (z.B. Transport nach unten, AirSeaAustausch, usw.) beeinflußt wird. Diese gute Übereinstimmung unterstützt die Aussage, daß die Hydrolyse eine bedeutende Rolle beim Abbau von COS im Seewasser spielt. Die berechnete HydrolyseLebenszeit ist mit dem Tagesmittel der COSKonzentration korreliert. Da die Tagesmittelwerte sowohl zeitliche wie auch räumliche Mittelwerte der COSKonzentrationen darstellen, zeigt diese Korrelation, daß Hydrolyse eine bedeutende Rolle in der raumzeitlichen Variabilität der COSKonzentration einnimmt. Da die Konzentration des gelösten COS von mehreren Faktoren abhängig ist, scheint eine multivariable Betrachtung sinnvoll. Hierfür wurde eine "Multiple Linear Regression Analysis'' (MLRA) ausgeführt. Diese Analyse ergibt ein empirisches Modell der folgenden Form für die Berechnung des Tagesmittels der COSKonzentration: [COS] = 1.8# 13log[Chl] - 1.5W s 0.057G - 0.73, mit [COS] = mittlere Konzentration von COS in pmol L -1 # = HydrolyseLebenszeit in Stunde [Chl] = mittlere Konzentration von Chlorophyll a in mg m -3 W s = Windgeschwindigkeit in m s -1 G = Intensität der Globalstrahlung in W m -2 . Die Parameter auf der rechten Seite der Gleichung können direkt oder indirekt von Satelliten aus gemessen werden, deshalb kann dieses Modell für die Abschätzung der Konzentration von COS im Seewasser anhand von Satelliten Daten verwendet werden. Das empirische Modell soll noch durch weitere Messungen bestätigt bzw. verbessert werden. Der Austauschfluß von COS zwischen der Atmosphäre und dem offenen Ozean wurde mit dem AirSeaFlußModell von Liss and Slater (1974) zusammen mit dem Modell von Erickson (1993) f
The timing and duration of leaf deployment strongly regulate earth-atmosphere interactions and biotic processes. Leaf dynamics therefore have major implications for life on earth, including the global energy balance, carbon and water cycles, feedbacks to climate, species extinction risk and agriculture. Evidence of shifts in the timing of leaf deployment and senescence (leaf phenology) as a result of climate change has been accumulating over the past decades, particularly in relation to spring phenology in the northern hemisphere. However, leaf phenological change in other parts of the world has received less attention. This thesis quantifies global phenological change over the past three decades using remotely sensed data. Phenological change was found to be widespread and severe, also in the southern hemisphere. While the detected change testifies of the phenological plasticity of many plant species, it is not clear if the duration of leaf deployment (leaf habit) is equally sensitive to environmental change. Since evergreen and deciduous leaf habits are often distinctly sorted along environmental gradients, ecologists have hypothesised that these patterns result from natural selection for an optimal leaf habit, under a given environmental regime. Such evolutionary convergence can be examined by testing if the physiological niche that is occupied by a particular leaf habit (evergreen or deciduous) is similar among regions with distinct evolutionary histories. Using a process-based model of plant growth and a constructed map of evergreen and deciduous vegetation, the physiological niche of leaf habits was quantified in four global biogeographic realms. Substantial niche overlap was found between the same leaf habit in different realms, suggesting evolutionary convergence of the physiological niche. This implies a sensitivity of leaf habit to environmental change, as environmental variables determine the geographic space where the physiological niche allows a positive carbon balance, and therefore occurrence of the leaf habit. Since the physiological niche consists of the integrated effects of physiological traits and trade-offs, environmental dependencies and leaf habit and phenology, an understanding of the carbon economy of individual plants requires decomposing the physiological niche into its components. Using empirical data on leaf phenology, leaf habit and physiological processes from woody species in a seasonally dry African savanna, a simple carbon balance model was parametrised. Carbon gain varied considerably between species as a result of substantial variation in leaf habit, leaf phenology and physiological traits. The multiple lines of evidence in this thesis therefore suggest that, while convergent selective forces may determine the dominant leaf habit in a particular environment, inter-specific variation is substantial, potentially as a consequence of historical contingencies or competitive interactions.
An eclogite barometer has profound importance in the study of upper mantle processes and potential application to diamond prospecting. Studies on the partitioning of Li between clinopyroxene (cpx) and garnet (grt) in natural samples have shown that this particular element is very sensitive to changes in pressure and could be calibrated as the barometer demanded for bimineralic eclogites. Experiments were performed from 4 to 13 GPa and 1100 to 1400°C in the CMAS (CaO, MgO, Al2O3, SiO2) system with Li added as Li3PO4 to quantify this pressure dependence into a barometer expressed in the following equation: P= (0.00255*T-lnKd)/0.2351 where P is in GPa, T is in °C and Kd is defined as the partition coefficient of Li (in ppm) between clinopyroxene and garnet. The experimental pressures are reproduced to ± 0.38 GPa (1σ) by this equation. This barometer is strictly applicable only to CMAS. Experiments at 1300°C, 8-12 GPa showed that Henry’s Law is fulfilled for Li partitioning between cpx and grt in the concentration range of approximately 0.01 – 1 wt% Li. Direct application of the equation to experiments in natural systems performed at 1300°C from 4 GPa to 13 GPa consistently overestimates pressures by approximately 2 GPa. Our previous experiments in the system CaO-MgO-Al2O3-SiO2 + Li3PO4 showed that the partitioning of Li between garnet and clinopyroxene is pressure dependent in eclogitic bulk compositions. This supports experimentally the hypothesis of Seitz et al. (2003), based on the analysis of Li in eclogitic xenoliths and inclusions in diamond, that the partitioning of this particular element between clinopyroxene and garnet is very sensitive to changes in pressure and could be calibrated as a barometer for bimineralic eclogites. In order to calibrate this pressure dependence into a barometer, experiments were performed in natural systems using starting materials sourced from a well preserved eclogitic xenolith from the Roberts Victor kimberlite pipe (South Africa) to extrapolate our findings in CMAS to natural systems. Sixteen multianvil experiments were performed from 4-13 GPa and 1100-1500°C. Our findings reinforced the general trend we observed in the CMAS system, that KdLi cpx-grt for Li decreases with increasing P, and that at P ≥ 12 GPa, garnet is able to incorporate more Li than clinopyroxene. Multiple linear regression was applied to our experimental results to create the barometer: P = (0.000963*T – ln KdLi cpx-grt + 1.581) / 0.252 Where P is pressure in GPa, T is temperature in °C and KdLi cpx-grt is defined as the partitioning coefficient of Li obtained by dividing the concentration of Li in cpx by the concentration of Li in garnet. This barometer reproduces the experimental conditions to ± 0.2 GPa. It is applicable to eclogitic xenoliths, to garnet pyroxenites and to peridotitic and eclogitic inclusions in diamond. Application of the barometer to diamond bearing xenoliths results in pressures in the diamond stability field. Clinopyroxene is easily corrupted in xenoliths and also preferentially takes in Li during short lived metasomatic processes. Care must be taken therefore to analyse primary, unaltered clinopyroxene. Our preliminary application to natural samples shows that the barometer can be applied beyond the experimental range to pressures down to 3 GPa. Seventeen eclogitic xenoliths were chosen from a sample set of greater than 200 for their fresh microscopic and macroscopic appearance and were analyzed for Li content in coexisting garnet (grt) and clinopyroxene (cpx). These samples can be subdivided into two groups on the basis of Mg in cpx (cpfu: cations per formula unit, based on 6 oxygens): Group 1 with Mg > 0.75, and Group 2 with Mg < 0.75. Group 1 xenoliths show lower Li contents in both grt and cpx compared to Group 2. The Li barom ter calibrated in Hanrahan et al. (2009b)/Chapter 3 was applied to these samples as well as available literature data to obtain pressures of provenance - Group 2 xenoliths often provide pressures that appear unrealistic for eclogitic xenoliths. In light of observed crystal chemical relations in the natural samples, a new fitting procedure was applied to the experimental data presented in Chapter 3. This new fit appears to be more realistic than the previous fit, although a strong relationship with Mg# remains present, suggesting that Li-barometry is, at present, only applicable to Mg-rich eclogites. Inclusions in diamond, with the exception of eclogitic inclusions of coexisting majorite and cpx, often yield pressures that are inconsistent with the pressures required for diamond formation. Although an interesting observation when comparing all of the data is that inclusions in diamond have significantly higher average Li concentrations compared to xenoliths, which suggests that Li is highly present in the fluids from which diamonds form in the mantle, an observation which was previously made for the deep mantle as a result of high Li in ferropericlase inclusions in diamond (Seitz et al. 2003).
Many natural minerals exist in the form of a solid solution. The systematic changes in structural and physical properties of oxide solid solutions are of geological importance and allow for wide applications. In order to understand the composition-structureproperty relations, substitutional solid solutions of CuxZn2−xTiO4, ZnxMg1−xTi2O5 and CuxMg1−xTi2O5 have been synthesised by mechanochemical activation assisted solid state synthesis. Self-propagating high-temperature synthesis has been employed to achieve the interstitial solid solutions of Ti5Si3Zx (Z refers to the element boron or oxygen).
The changes in the crystal structure and physical properties due to the formation of solid solutions are investigated by employing X-ray diffraction, neutron diffraction, Raman spectroscopy, low-temperature heat capacity, thermal expansion, scanning electron microscopy, UV-vis spectroscopy, plane-wave ultrasound spectroscopy and density functional theory calculations.
This thesis presents microstructural investigations of rock salt from the central part of the Gorleben salt dome (Northern Germany). The main emphasis was to characterize the rock salt microfabrics, to identify operating deformation mechanisms in halite and anhydrite and to decipher the macro- and microstructural distribution of hydrocarbons, which have been encountered during the underground exploration of the salt dome. The microfabrics of the Knäuel- and the Streifensalz formation indicate that strain-induced grain boundary migration has been active during deformation of halite. Crystal plastic deformation of halite is further documented by lattice bending, subgrain formation and minor subgrain rotation. Evidence for pressure solution of halite has not been found, but cannot be excluded because of the small grain size, the lack of LPO and the low differential stress (1.1 - 1.3 MPa) as deduced from subgrain-size piezometry. Solution precipitation creep was proven for intercalated anhydrite layers and clusters, which have been deformed in the brittle-ductile regime. Brittle deformation of anhydrite in terms of boudinage and fracturing was counteracted by viscous creep of halite which caused a re-sealing of fractures and a reestablishing of the characteristic sealing capacity of rock salt. Hydrocarbons are mainly located along cross cut 1 West of the Gorleben exploration mine and are heterogeneously distributed in the rock salt. They are incorporated in the rock salt foliation in the form of streaks, dispersed clouds, clusters and isolated patches. On the micro-scale, hydrocarbons are trapped along grain boundaries of halite and/or anhydrite, in micro-capillary tubes of anhydrite and in pore space of the rare rock salt with elevated porosity (< 1.26 vol.-%). Such elevated porosities correlate with elevated hydrocarbon concentrations of several hundred ppm. The overall concentrations of hydrocarbons, however, are very low (< 0.05 wt.-%). Elevated porosity is depicted to be a remnant originating from an early stage of salt uplift when fluid and hydrocarbons have migrated and spread from the Staßfurt Karbonat (z2SK) into the superjacent Gorleben Hauptsalz. During halokinesis and the strong reworking of the salt body hydrocarbons have been redistributed and dismembered resulting in the isolated present-day occurrences. The distribution of hydrocarbons shows no relation to local variations in the rock salt fabric. The microstructures of hydrocarbon-bearing and hydrocarbon-free Gorleben rock salt are not distinguishable from each other. Likewise, the presence of hydrocarbons should not have influenced the mechanical behavior or the rock salt as indicated by the microfabrics studied and by geomechanical data. The pure amounts of hydrocarbons are too low for any detectable impact on the barrier properties of this part of rock salt. Although hydrocarbons have migrated into the Gorleben Hauptsalz during an early stage of salt uplift when the sealing capacity of rock salt was diminished, the major implication of their isolated distribution patterns is that the Gorleben rock salt was able to regain its sealing capacity during subsequent deformation and re-equilibration. Former migration pathways for fluid and hydrocarbons have been healed and do not exist anymore. The application of X-ray computed tomography (CT) allows the 3D visualization and quantification of anhydrite, pore space and fluid phases located along grain-boundaries or trapped as intracrystalline inclusions. The 3D reconstruction of anhydrite clusters and pore space for the same sample reveals different spatial distribution patterns. This fact implies that anhydrite is not responsible for such elevated pore space in the rock salt studied, which has been largely closed during the polyphase deformation history of the Gorleben salt dome. High-resolution nanoCT scans (≤ 1 μm voxel size) of single intra- and intercrystalline fluid inclusions in rock salt enable a characterization of gaseous, solid and liquid phases inside single fluid inclusions and give exact information on morphology and shape. The 3D reconstruction of grain boundary fluid inclusions allows the amount, volumes, surface areas or diameters of various types to be determined. Non-destructive X-ray CT imaging is presented as very useful tool to characterize the structural inventory of rock salt. This non-destructive technique offers new perspectives for microstructural studies and for a wide range of research in structural geology, in general.
Semi-arid African ecosystems influence trends and variability in global terrestrial carbon dynamics. However, there are uncertainties in potential effects of future climates for semi-arid ecosystems, especially for niche ecosystems. At the same time, African ecosystems provide the livelihoods and ecosystem services for around 1.4 billion people. Future population growth and associated changes in land use pose a challenge for the protection of African biodiversity. Therefore, this work focussed on future impacts of climate change on African ecosystems and carbon dynamics and also for African protected areas (PAs), where they may cooccur with other global change factors. Another focus was on uncertainties associated with future projections and with modelling the Nama Karoo, as an example of a semi-arid niche ecosystem. Dynamic vegetation models (DVMs) were the main research tool.
In Chapter 2, we analysed climate change impacts on African ecosystems and carbon pools until the end of the 21st century and associated uncertainties based on an ensemble of vegetation simulations with the DVM adaptive dynamic vegetation model (aDGVM). We investigated the impact of increased atmospheric CO2 concentrations and two climate change scenarios (medium (RCP4.5) and high emissions (RCP8.5); RCP - representative concentration pathway) on vegetation changes. Differences in the simulated vegetation were primarily driven by assumptions about the influence of CO2 on plants. Elevated CO2 concentrations led to increased total aboveground vegetation biomass and shrub encroachment into grasslands and savannas for both climate scenarios. In simulations without the direct influence of CO2 on plants, there was hardly any shrub encroachment and vegetation biomass decreased or varied between a slight decrease in some cases and a slight increase in others. Based on these results, biome changes due to climate change are likely in Africa in the future. Due to the large uncertainties in future projections, strategies to adapt to climate change must be flexible.
The simulated vegetation in Chapter 2 represented potential, natural vegetation and is particularly suitable to investigate PAs. However, PAs do not exist isolated from their environment and social developments. In Chapter 3, the vegetation projections with CO2 effect from Chapter 2 were combined with projections for population density and land use. Except for many PAs in North Africa, most PAs were adversely affected by at least one of the three drivers by the end of the 21st century in both investigated scenarios ("middle-of-the-road" and "fossil-fuelled development"). Cooccurrence of the drivers varied by region and scenario for PAs. Both scenarios implied increasing challenges for the conservation of African biodiversity in PAs. The impact of climate change on vegetation is likely to be exacerbated by socio-economic change for most African PAs. Strong mitigation of future climate change together with equitable societal development may facilitate successful ecosystem conservation.
The simulations in Chapters 2 and 3 showed large-scale patterns of vegetation change, but their low resolution makes them unsuitable for local analyses. In Chapter 4, the challenges of simulating smaller scale, semi-arid ecosystems and their carbon cycle were analysed for the Nama Karoo with the aDGVM2 and its shrub module. The aDGVM2 is based on the aDGVM, but represents plants more flexibly. In all tested aDGVM2 configurations, the carbon fluxes improved compared to initial simulations but still overestimated them. The measured morphology of the dwarf shrubs and soil water dynamics were not reproduced in aDGVM2. Semi-arid soil water dynamics and coping strategies of semi-arid dwarf shrubs under drought stress are not adequately implemented in the aDGVM2. Further field research on semi-arid water and carbon dynamics of vegetation is necessary to parameterise the aDGVM2 for dwarf shrubs. If these challenges are overcome, DVMs can be a powerful tool for much-needed research on the impacts of climate change on the Nama Karoo.
The analyses have shown that climate change under medium to high emission scenarios is likely to lead to large-scale changes in ecosystems and the carbon balance in Africa. Because lower emissions scenarios come with less uncertainty, climate change adaptation strategies likely need to be less complex or extensive if climate change is minimised. For African PAs, the challenges of climate change may be exacerbated by socio-economic factors to a regionally varying extent. This research suggests that successful ecosystem conservation depends on climate change mitigation measures and ensuring equitable, sustainable development. The shown uncertainties, e.g., in the implementation of the CO2 effect on plants or vegetation dynamics in more niche ecosystems, help to focus future research efforts and increase our understanding of the range of plausible futures we may need to adapt to.
This thesis deals with the analysis of “presolar” silicates and oxides by high resolution mass spectrometry and electron microscopy techniques. This “stardust” was identified by its extreme oxygen isotopic anomalies, which point to nucleosynthetic reactions in stellar interiors, in the carbonaceous chondrite Acfer 094. Isotopic, chemical and mineralogical studies on these stardust grains therefore allow the testing of astrophysical questions on Earth, which are otherwise only accessible by spectroscopy and theoretical models. The class of presolar silicates has been identified only six years ago in 2002, although it was known already from spectroscopic observations that silicates represent the most abundant type of dust in the galaxy. The development of the “NanoSIMS” was a crucial step in this respect, because this ion probe with its superior spatial resolution of only 50 nm allowed the detection of the typically 300 nm sized presolar silicates. A total of 142 presolar silicates and 20 presolar oxides were identified within Acfer 094, whose matrix therefore contains 163 ± 14 ppm presolar silicates and 26 ± 6 ppm presolar oxides. This is among the highest amounts reported so far for any primitive solar system material. The majority of detected stardust grains derive from asymptotic giant branch stars of 1 – 2.5 Msun and close-to-solar or slightly lower-than-solar metallicity. However, by measuring the Si isotopic compositions of some enigmatic grains, it could be shown that there is a sub-class of presolar silicates characterized by an extreme enrichment of 17O and a moderate enhancement of 30Si relative to solar, whose origins might be explained by formation in binary stellar systems. About 10% of all grains exhibit an enrichment in 18O and some of them also of 28Si relative to solar, which most likely point to an origin in type II supernova explosions. The Si isotopic measurements also allowed to quantify the effect of the s-process on the Si isotopes in low-mass asymptotic giant branch stars. The results agree well with theoretical predictions. The grains were furthermore characterized by SEM and the chemistries of about half of the grains were determined by Auger electron spectroscopy. The majority of grain morphologies are consistent with what is expected from condensation experiments. However, a lot of grains are altered by Fe-rich minerals, which are either of primary condensation or of secondary ISM or solar nebula origin. Furthermore, complex presolar grains consisting of refractory Al-rich grains attached to silicate material could be identified, which have been predicted by condensation theory and observational evidence. Nine presolar silicates were analyzed by combined NanoSIMS/TEM studies. The majority of grains are Mg-rich and amorphous, which is in contrast to astrophysical evidence, which mainly postulate crystalline Mg-rich and amorphous Fe-rich circumstellar condensates. However, the grains might have been rendered amorphous by secondary processes in the ISM or could have condensed under non-equilibrium, low-temperature conditions in the circumstellar outflow. The grains are more likely characterized by a variable, pyroxene-like chemistry, which could be a result of sputtering in the ISM, which preferentially removes Mg. The detected crystalline presolar silicates in this study and in other work are all olivines, whereas grains with a pyroxene stoichiometry are all amorphous except one. This supports astrophysical models which point to different formation pathways for these two types of grains and therefore different crystallinity. However, the relatively high Fe content of three detected presolar olivines in this study and in other work is in contrast to astrophysical evidence and theoretical considerations, which predict essentially Fe-free crystalline grains. It is therefore possible that the infrared spectra might also be compatible with less Mg-rich olivines. The only crystalline presolar silicate with a pyroxene-like stoichiometry is the unusual grain 1_07: although it is chemically enstatite, the electron diffraction pattern could only be indexed to silicate perovskite, which is stable above ~23 GPa. The discovery of a high-pressure phase of presolar origin shows that dust grains encountering interstellar shocks might not necessarily be completely destroyed. In astrophysical models it is in principle also possible that a fraction of larger grains might survive such a shock wave encounter as a high-pressure modification, which is supported by this discovery.
Das Thema dieser Arbeit war die Untersuchung der natürlichen Variationen von den zwei primordialen Uranisotopen (238U und 235U) mit einem Schwerpunkt auf Proben, die (1) die kontinentale Kruste und ihre Verwitterungsprodukte (d.h. Granite, Shales und Flusswasser) repräsentieren, (2) Produkte der hydrothermalen Alteration vom mittelozeanischen Rücken widerspiegeln (d.h. alterierte Basalte, Karbonatgänge und hydrothermales Wasser) und (3) aus abgegrenzten euxinischen Becken (d.h. Proben aus der Wassersäule und den dazugehörigen Sedimenten) stammen. Das allgemeine Ziel war das Verständnis, unter welchen Bedingungen und Mechanismen eine Fraktionierung der zwei häufigsten Uranisotope (238U und 235U) in der Natur erfolgt, zu verbessern.
Die untersuchten Haupt- und Nebenflüsse unterscheiden sich sowohl in Ihrer Urankonzentration (c(U)) als auch in Ihrer Uranisotopenzusammensetzung (δ238U), wobei die Nebenflüsse eine geringere Urankonzentration (0.87 nmol/kg bis 3.08 nmol/kg) und eine schwerere Uranisotopenzusammensetzung aufweisen (-0.29 ‰ bis +0.01 ‰ im δ238U) im Vergleich zu den Hauptflüssen (c(U) = 5.19 nmol/kg bis 11.69 nmol/kg und d238U = -0.31 ‰ bis +0.13 ‰) aufweisen. Die untersuchten Gesteinsproben fallen alle in einen recht schmalen Bereich von δ238U, zwischen -0.45 ‰ und -0.21 ‰, mit einem Durchschnittswert von -0.30 ‰ ± 0.04 ‰ (doppelte Standardabweichung). Deren Uranisotopenvariationen sind unabhängig von der Urankonzentration (11.8 µg/g bis 1.3 µg/g), dem Alter (3.80 Ga bis 328 Ma), der Probenlokalität und Grad der Differenzierung. Basierend auf den Ergebnissen der Hauptflüsse, die die Uranhauptquelle für den Ozean darstellen, schlagen wir für zukünftige Berechnungen in der Massenbilanz des Urans einen neuen Wert als beste Abschätzung für die Quelle des Urans im Ozean vor, δ238U = -0.23 ‰.
Die Produkte der hydrothermalen Alteration, alterierte Basalte und Kalziumkarbonatgänge, zeigten etwas stärkere Isotopenvariationen (δ238U zwischen -0.63 ‰ und +0.27 ‰) als erwartet und die hydrothermalen Fluide wiesen eine etwas leichtere Uranisotopenzusammensetzung als Meerwasser ((-0.43 ± 0.25) ‰ vs. (-0.37 ± 0.03) ‰) auf. Diese Ergebnisse sind in Übereinstimmung mit einem Modell, dass annimmt, dass die beobachtete Isotopenfraktionierung hauptsächlich ein Ergebnis von Redoxprozessen ist, z.B. die partielle Reduktion von löslichem UVI aus dem Meerwasser während der hydrothermalen Alteration, was zu einer Anreicherung der schweren Uranisotope in der reduzierten Uranspezies (UIV) führt und 2) das bevorzugte Entfernen von UIV aus den hydrothermalen Fluid und der Einbau in die alterierte ozeanische Kruste. Durch diesen Prozess wird das hydrothermale Fluid an schweren Uranisotopen verarmt und somit würden auch die alterierten Basalte und Karbonate ein niedriges δ238U aufweisen, wenn sie mit dem isotopisch leichten hydrothermalen Fluid in Kontakt gekommen sind.
Die Untersuchung von Wasser- und Sedimentproben aus der Ostsee und dem anoxischen Kyllaren Fjord (Norwegen) auf deren Uran- und Mo-Isotopenzusammensetzung zeigte, dass die Uranisotopenzusammensetzung der Sedimente abhängt von (1) dem Ausmaß des Uranaustrags aus der Wassersäule (in einer ähnlichen Art und Weise wie bei den Molybdänisotopen) und (2) der Sedimentationsrate, d.h. der Fraktion von authigenem- relativ zum dedritischen Uran in den Sedimenten. Aufgrund der hohen Sedimentationsrate zeigen die Sedimente aus dem Kyllaren Fjord nur eine moderate authigene Urananreicherung und eine leichtere Uranisotopenzusammensetzung als Sedimente aus dem Schwarzen Meer. In den anoxischen Becken der Ostsee erfolgt dagegen eine starke Mo- und schwache U-Isotopenfraktionierung zwischen Wasser und Sediment. Durch die regelmäßigen auftretenden Spülereignisse mit sauerstoffreichem Wasser wurden vermutlich die ursprünglichen anoxischen Mo- und U-Isotopensignaturen der Sedimente verändert. Demzufolge müssen die Sedimente durchgehend anoxischen Bedingungen ausgesetzt sein, um eine Mo- und U-Isotopensignatur von den Redoxbedingungen während der Ablagerungen zu speichern.
Der Vergleich zwischen Molybdän- und Uranisotopen in der Ostsee und dem anoxischen Kyllaren Fjord zeigte, dass sich Uran- und Molybdänisotope in stark euxinischen Wassersäulen (c(H2S) > 11 µmol/L) entgegengesetzt verhalten. Dementsprechend ergänzen sich die beiden Isotopensysteme und können genutzt werden, um die Ablagerungsbedingungen in abgeschlossenen Becken und die Redoxentwicklung des Paläoozeans zu untersuchen.
Occurrence and sources of 2,4,7,9-tetramethyl-5-decyne-4,7-diol (TMDD) in the aquatic environment
(2011)
The aim of the present study was to identify the sources of 2,4,7,9-tetramethyl-5-decyne-4,7-diol (TMDD) into the aquatic environment and to investigate its occurrence in rivers and wastewater treatment plants (WWTPs). Therefore, TMDD was analyzed in 441 wastewater samples from influents and effluents of 27 municipal WWTPs, in 6 sludge samples, in 52 wastewater samples from 3 sewage systems of municipal WWTPs, in 489 surface samples from 24 rivers, in 9 wastewater samples of 3 paper-recycling industries and in 65 groundwater samples. TMDD was also analyzed in household paper products, in 23 samples of toilet
papers, in 5 types of paper towels and in 12 types of paper tissues. The samples were collected between 2007 and 2011. The water samples were extracted with solid phase extraction (SPE) and the household paper samples with Soxhlet extraction. Gas chromatography-mass spectrometry (GC-MS) was used for quantification purposes. Between November 2007 and January 2008, TMDD was detected in the river Rhine at Worms with permanent high concentrations (up to 1330 ng/L). The results showed that TMDD is uniformly distributed across the river at Worms. An increase of the mean TMDD concentration from approximately 500 ng/L to 1000 ng/L was registered in January 2008. Due to the minor fluctuations of the TMDD concentration during the sampling period it is expected that the input of TMDD into the river is continuous. Therefore, TMDD might rather originate from effluents of municipal WWTPs than from temporal sources. The mean TMDD load based on the analysis of 147 water samples collected in the River Rhine was 62.8 kg/d which is equivalent to 23 t/a suggesting that TMDD must be used and/or produced in high quantities in order to be found in those high concentrations. To determine if TMDD is discharged by effluents of municipal WWTPs into the rivers, 24 hours influent and effluent samples of four municipal WWTPs in the Frankfurt/Rhine-Main metropolitan region were collected during November 2008 and February 2010 and analyzed for TMDD. The TMDD influent concentrations varied between 134 ng/L and 5846 ng/L and the effluent concentrations between <LOQ (limit of quantitation) and 3539 ng/L. The TMDD elimination rates in the four WWTPs varied between 33% and 68%. The results showed that effluents of municipal WWTPs are an important source of TMDD in the aquatic environment because TMDD is not completely removed from the sewage during the wastewater treatment. Weekly and daily variations of the TMDD concentration in the influents of two municipal WWTPs indicated that both private households and indirect industrial dischargers contribute to the introduction of TMDD into the municipal sewage systems. A more detailed study of the TMDD elimination rate in the different wastewater treatment stages was carried out in the WWTP Niederrad/Griesheim in Frankfurt am Main. The results showed that the removal of TMDD is mainly carried out during the aerobic biological treatments, where the elimination rate was 46%. In contrast, during the anoxic treatment the removal efficiency was only 1.4% and during the mechanical treatment the elimination rate was 19%. To determine the sources of TMDD in the sewage, household paper products (paper tissues, toilet papers and paper towels) were analyzed for TMDD using Soxhlet extraction. TMDD was detected in 83% of the samples (n=40). The highest mean TMDD concentrations were found in recycled toilet paper (0.20 μg/g) and in paper towels (0.11 μg/g). In paper tissues and non-recycled toilet paper the mean TMDD concentrations were lower 0.080 μg/g and 0.025 μg/g respectively. According to these results the high TMDD influent concentrations found previously in municipal WWTPs (mean 1.20 μg/L) cannot be explained due to migration of TMDD from the household paper products into the sewage. Thus indirect industrial dischargers are the cause of the high influent TMDD concentrations. Effluents of municipal WWTPs with different indirect industrial dischargers (textile-, metal processing-, food processing-, electroplating-, paper-recycling- and printing ink factories) were analyzed. The highest mean TMDD concentrations were found in the effluents of municipal WWTPs that have paper-recycling (71.3 μg/L) and printing ink factories (138 μg/L) as indirect industrial dischargers. These results were confirmed by analyzing process wastewater of three paper-recycling factories located in Germany. High TMDD concentrations were detected and fluctuated between 1.83 μg/L and 113 μg/L. TMDD was also analyzed in the wastewater of a non-recycling-paper factory but its concentration was much lower (0.066 μg/L) indicating that TMDD is introduced into the processing water during the papermaking process due to the use of waste paper. Analyses of wastewater samples from different parts of the sewage pipes of a municipal WWTP in Hesse, which receives the wastewater from a printing ink factory, were carried out. The TMDD concentration in the wastewater sample from the sewage pipe of the printing ink factory was much higher (3,300 μg/L) than the TMDD concentration detected in the other wastewater samples from the sewage system (0.030 μg/L – 0.89 g/L). These results confirm the printing ink production as one of the principal sources of TMDD in the sewage. Analysis of surface water samples of the River Modau downstream from the effluent of the WWTP Nieder-Ramstadt showed TMDD concentrations of up to 28.0 μg/L. These high TMDD concentrations might be caused by the indirect wastewater discharges of a paint factory connected to the municipal sewage system. These results indicate that TMDD is introduced into the municipal WWTPs principally by indirect industrial dischargers and they are mainly paint and printing ink factories. The paper-recycling factories also represent an important source of TMDD in municipal WWTPs but indirectly. According to statements given by the representatives of two paper recycling factories neither TMDD or any other TMDD containing product is used or added during the papermaking process. Therefore, TMDD is washed out from the printing inks of the coloured waste paper and concentrated in the process wastewater in the closed water circuits of paper-recycling factories reaching rivers and municipal WWTPs. The occurrence and distribution of TMDD in surface waters in Germany was also studied. The results showed that TMDD is widely distributed across different rivers systems in the federal states of Hesse, North-Rhine-Westphalia, Bavaria, Baden-Wuerttemberg and Rhineland-Palatinate. In Hesse, TMDD was detected in the some of main rivers with mean concentrations of 812 ng/L (Schwarzbach, Hessian Ried), 374 ng/L (Kinzig), 393 ng/L (Main, at Frankfurt), 539 ng/L (Werra), 326 ng/L (Fulda), 151 ng/L (Emsbach) and 161 ng/L (Nidda). In small rivers (creeks) the mean TMDD concentrations varied between <LOQ (Diemel, Urselbach) and 1890 ng/L (Darmbach). The results showed that the TMDD concentrations in creeks are highly influenced by both effluents of WWTPs and by the distance between the sampling point and the nearest WWTP. Surface samples from sampling locations downstream from WWTPs dischargers showed higher TMDD concentrations (mean 518 ng/L) than sampling locations upstream from WWTPs dischargers (mean 35.1 ng/L). The behavior of TMDD during bank filtration was investigated at two locations, at a water utility company at the Lower River Rhine (urban area) and at the Oderbruch polder (rural area). The results indicated that TMDD is removed from the surface water by bank filtration at both sampling locations. The removal process is probably carried out in the first meters of the aquifer (hyporheic zone) by biodegradation processes, since TMDD does not tend to be absorbed by sediments and it was not found in the groundwater of monitoring wells. In groundwater samples from the Hessian Ried (n=23) TMDD was found only in five samples and the highest TMDD concentration was 135 ng/L. According to these results, TMDD does not represent a concern for drinking water in Germany, since it does not reach the groundwater with high concentrations and it has a low toxicity potential. The input of TMDD into the North Sea was estimated to be 60.7 t/a by considering the mean transported loads of TMDD by the River Rhine at Wesel (58.3 t/a) and Meuse in the Netherlands (2.40 t/a). The estimated discharge of TMDD by German municipal WWTPs (8.19 t/a) and paper-recycling factories (9.24 t/a) into rivers seems to be too low considering that the mean TMDD load in the River Rhine downstream from Wesel is 58.3 t/a. However, due to the high density of population and industries at the Lower Rhine it is expected that more relevant sources of TMDD are located along the Rhine River increasing the transported load. According to the results of this PhD project TMDD is a non-ionic surfactant contained in products, which are applied on surfaces (printing inks and paints) and has the potential to reach the aquatic environment. Therefore, TMDD should fulfill the requirement of a biodegradability of 80% established by the “Law on the Environmental Impact of Detergents and Cleaning Products” in Germany. However, due to the partial elimination rates of TMDD obtained in municipal WWTPs (between 33% and 68%) and to the absence of information about the execution of the biodegradation test on TMDD, it is unknown if TMDD is in accordance with this law. Otherwise, its use as surfactant in such products is questionable.
The objective of the present doctoral thesis was to investigate the occurrence, distribution, and behaviour of six hydrophilic ethers: ethyl tert-butyl ether (ETBE), 1,4-dioxane, ethylene glycol dimethyl ether (monoglyme), diethylene glycol dimethyl ether (diglyme), triethylene glycol dimethyl ether (triglyme), and tetraethylene glycol dimethyl ether (tetraglyme) in surface-, waste-, ground- and drinking water samples. Solid phase extraction and gas chromatography/mass spectrometry were used to analyze the six hydrophilic ethers. Altogether more than 150 surface water samples, almost 100 of each groundwater and wastewater samples, and 10 raw and drinking water samples were analyzed during the research project.
Initially, the method was validated in order to simultaneously determine the analytes of interest in various aquatic environments. A solid phase extraction method that uses coconut charcoal (Resprep® activated coconut charcoal, Restek) or carbon molecular sieve material (SupelcleanTM Envi-CarbTM Plus, Supelco) for analyte absorption were found suitable for determination of ETBE, 1,4-dioxane, and glymes in surface-, drinking-, ground- and wastewater samples. Precision and accuracy of both methods was demonstrated for all analytes of interest. The recovery of target compounds from the ultrapure water spiked at 1.0 µg L−1 was between 86.8 % and 98.2 %, with relative standard deviation below 6 %. The samples spiked at 10.0 µg L−1 gave slightly higher recovery of 90.6 % to 112.2 % with a relative standard deviation below 3.4 % for each analyte. Detection and quantification limits in ultrapure water and surface waters were furthermore established. The limit of quantitation (LOQ) in ultrapure water ranged between 0.024 µg L−1 to 0.057 µg L−1 using Restek cartridges, and 0.030 µg L−1 to 0.069 µg L−1 using Supelco cartridges. In the surface water samples the calculated LOQ was 0.032 µg L−1 to 0.067µg L−1 using coconut charcoal material and 0.032 µg L−1 to 0.052 µg L−1 using the carbon molecular sieve material. Moreover, stability of the unpreserved and preserved water samples as well as the extracts was determined. Preservation of samples with sodium bisulfate (at 1 gram per Liter) resulted in much better stability of the ethers in water samples. Subsequently, 27 samples obtained from seven surface water bodies in Germany (Rivers Rhine, Lippe, Main, Oder, Rur, Schwarzbach and Wesel-Datteln Canal) were analyzed for the six hydrophilic ethers. ETBE was present in only two surface waters (Rhine River and Wesel-Datteln Canal) with concentrations close to the LOQ (up to 0.065 µg L−1). 1,4-Dioxane was detected in all of the water samples at concentrations reaching 1.93 µg L–1. Monoglyme was identified only in the Main and Rhine Rivers at the maximum concentration of 0.114 µg L–1 and 0.427 µg L–1, respectively. Very high concentrations (up to 1.73 µg L−1) of diglyme, triglyme, and tetraglyme were detected in the samples from the Oder River. These glymes were also detected in the Rhine River; however the concentrations did not exceed 0.200 µg L–1. Furthermore, tetraglyme was detected in the Main River at an average concentration of 0.409 µg L–1 (n = 6) and in one sample from the Rur River at 0.192 µg L–1.
Four sampling campaigns were conducted at the Oderbruch polder between October 2009 and May 2012, in order to study the behavior of the hydrophilic ethers and organophosphates during riverbank filtration and in the anoxic aquifer. Moreover the suitability of these target compounds was assessed for their use as groundwater organic tracers. At the time of each sampling campaign, concentrations of triglyme and tetraglyme in the Oder River were between 20–185 ng L–1 (n = 4) and 273¬–1576 ng L–1 (n = 4). Monoglyme, diglyme, and 1,4-dioxane were analyzed only during the two last sampling campaigns. At that time, the concentration of diglyme in Oder River was 65¬–94 ng L-1 (n = 2) and 1,4-dioxane 1610¬–3290 ng L–1 (n = 2). In the drainage ditch, following bank filtration, concentrations of ethers ranged between 1090 ng L–1 and 1467 ng L–1 for 1,4-dioxane, 23¬ng L–1 and 41 ng L–1 for diglyme, 37 ng L–1 and 149 ng L–1 for triglyme, and 496 ng L–1 and 1403 ng L–1 for tetraglyme. In the anoxic aquifer, 1,4-dioxane showed the greatest persistence during the groundwater passage. At the distance of 1150 m from the river and an estimated groundwater age of 41.9 years, a concentration above 200 ng L−1 was detected. A positive correlation was found for the inorganic tracer chloride (Cl−) with 1,4-dioxane and tetraglyme. Similarities in the behavior of Cl− and the organic compound suggested that 1,4-dioxane and tetraglyme are controlled by the same hydraulic process and therefore can be used as additional tracers to study the dynamics of the groundwater system. These results show that high concentrations of ethers are present in the surface water and are not removed during bank filtration processes. Moreover, the hydrophilic ethers persist in the anoxic aquifer and little or no degradation is expected, supporting, their possible application as organic tracers.
A separate sampling project was conducted for 1,4-dioxane that focused primarily on its fate in the aquatic environment. This study provided missing information on the extent of water pollution with 1,4-dioxane is Germany. Numerous waste-, surface-, ground- and drinking water samples were collected in order to determine the persistence of 1,4-dioxane in the aquatic environment. The occurrence of 1,4-dioxane was determined in wastewater samples from four municipal sewage treatment plants (STP). The influent and effluent samples were collected during weekly campaigns. The average influent concentrations in all four plants ranged from 262 ± 32 ng L−1 to 834 ± 480 ng L−1, whereas the average effluents concentrations were between 267 ± 35 ng L−1 and 62,260 ± 36,000 ng L−1. The source of increased 1,4-dioxane concentrations in one of the effluents was identified to originate from impurities in the methanol used in the postanoxic denitrification process. Spatial and temporal distribution of 1,4-dioxane in the river Main, Rhine, and Oder was also examined. Concentrations reaching 2,200 ng L−1 in the Oder River, and 860 ng L−1 in both Main and Rhine River were detected. The average load during the sampling was estimated to be 6.5 kg d−1 in the Main, 34.1 kg d−1 in the Oder, and 134.5 kg d−1 in the Rhine River. In all of the sampled rivers, concentrations of 1,4-dioxane increased with distance from the mouth of the river and were found to negatively correlate with the discharge of the river. In order to determine if 1,4-dioxane can reach drinking water supplies, samples from a Rhine River bank filtration site and potable water from two drinking water production facilities were analyzed for the presence of 1,4-dioxane in the raw water and finished potable water. The raw water (following bank filtration) contained 650 ng L−1 to 670 ng L−1 of 1,4-dioxane, whereas the concentration in the finished drinking water fell only to 600 ng L−1 and 490 ng L−1, respectively.
During the final project, investigations of the source identification of high glyme concentrations in the Oder River were carried out. During four sampling campaigns between January, 2012 and April, 2013, 50 samples from the Oder River in the Oderbruch region and Poland were collected. During the first two samplings in the Oderbruch polder, glymes were detected at concentration reaching 0.07 µg L-1 (diglyme), 0.54 µg L−1 (triglyme) and 1.73 µg L−1 (tetraglyme) in the Oder River. The extensive sampling campaign of the Oder River (about 500 km) in Poland helped to identify the area of possible glyme entry into the river. During that sampling the maximum concentrations of triglyme and tetraglyme were 0.46 µg L−1 and 2.21 µg L−1, respectively. A closer investigation of the identified area of pollution, helped to determine the possible sources of glymes in the Oder River. Hence, the final sampling focused on the Kaczawa River, a left tributary of the Oder River and Czarna Woda, a left tributary of Kaczawa River. Moreover, samples from an industrial wastewater treatment plant were collected. Samples from Czarna Woda stream and Kaczawa River contained even higher concentrations of diglyme, triglyme, and tetraglyme, reaching 5.18 µg L−1, 12.87 µg L−1 and 80.81 µg L−1, respectively. Finally, three water samples from a wastewater treatment plant receiving influents from a copper smelter were analyzed. Diglyme, triglyme, and tetraglyme were present at an average concentration of 569 µg L−1, 4300 µg L−1, and 65900 µg L−1, respectively in the wastewater. Further research helped to identify the source of the glymes in the wastewater. The gas desulfurization process – Solinox implemented in the nearby copper smelter uses glymes as physical absorption medium for sulfur dioxide.
Results of this doctoral research provide important information about the occurrence, distribution, and behavior of hydrophilic ethers: 1,4-dioxane, monoglyme, diglyme, triglyme, and tetraglyme in the aquatic environment. A method capable of analyzing a wide range of ether compounds: from a volatile ETBE to a high molecular weight tetraglyme was validated. 1,4-Dioxane and tetraglyme were found to be applicable as organic tracers, since they are not easily attenuated during bank filtration and the anoxic groundwater passage. The extent of water pollution with 1,4-dioxane was shown in waste-, surface-, ground-, and drinking waters. One source of extremely high concentrations of 1,4-dioxane in a municipal sewage treatment plant applying postanoxic denitrification was identified, however more information is needed on the entry of 1,4-dioxane into surface waters. Moreover, 1,4-dioxane was present in drinking water samples from river bank filtration, which demonstrates its persistence in the aquatic environment and its low degradation potential during bank filtration and subsequent water treatment. Furthermore, this was the first study that focused primarily on identifying sources of glymes in surface waters. Glymes find a widespread use in industrial sectors, hence establishing their origin in the surface water is difficult (as with 1,4-dioxane). In this work, a gas desulphurization process was identified to be a dominating source of glyme pollution in the Oder River.
Despite mounting evidence of the anthropogenic influence on the Earth's climate, underlying mechanisms of climate change often remain elusive. The investigation of periods of rapid climate change from geological archives may provide crucial information about magnitude, duration, teleconnections of and regional responses to global and hemispheric scale climate perturbations. Thus, paleoclimate reconstructions may help in mitigating and adapting to the challenges of the coming decades. The '8.2 kyr B.P. climatic event' has previously been proposed as a possible analogue for the future climatic scenario of a reduced Atlantic Meridional Overturning Circulation (AMOC). The catastrophic drainage of the Laurentide meltwater lakes through the Hudson Bay and into the Labrador Sea, that occurred ca. 8.47 kyr B.P., caused the slowdown of the AMOC around 8.2 kyr B.P.. Subsequently, reduced heat transfer towards Europe triggered a substantial decline in (winter) temperature and pronounced changes in atmospheric circulation patterns in many regions of the northern hemisphere, especially the North Atlantic realm and Europe. Among the regions affected by the 8.2 kyr B.P. climatic event, the Eastern Mediterranean region is of particular interest for both past and future climate developments. Traditionally characterized as a region highly sensitive to variations in the climate systems of the high and low latitudes, abrupt climate changes have the potential to strongly alter atmospheric circulation patterns and thus precipitation distribution in the region that may have severe socioeconomical consequences. The analysis of stable hydrogen (δD) and oxygen isotopes (δ18O) in precipitation is an excellent tool to trace changes in atmospheric circulation. Here, we present a comparative study of δD and δ18O in precipitation from the Eastern Mediterranean region both in a present day scenario and during the 8.2 kyr B.P. climatic event. We analyze the influences of topography, air mass trajectory, climate and seasonality among others the stable isotopic compositions of meteoric waters from the Central Anatolian Plateau (CAP), Turkey, in order to create a first-order template which may serve as a reference against which paleoenvironmental proxy data may be more accurately interpreted and tested. Further, we employ a multiproxy approach on the early Holocene peat deposits of the classical site of Tenaghi Philippon (TP), NE Greece, to investigate paleoenvironmental responses to northern hemisphere cooling during the 8.2 kyr B.P. climatic event and aim to determine changes atmospheric circulation from δD of leaf wax n-alkanes (δDwax).
Based on δD and δ18O data from more than 480 surface water samples from the CAP, we characterize moisture sources affecting the net isotopic budget of precipitation, manifesting in a systematic north-south difference in near-sea level moisture compositions. Rainout, induced by the major orographic barriers of the plateau, the Pontic Mountains to the north and the Taurus Mountains to the south, strongly shape the modern patterns of δD and δ18O. Stable isotope data from the semi-arid plateau interior provide clear evidence for an evaporitic regime that drastically affects surface water compositions. Strong evaporative enrichment contrasts rainfall patterns along the plateau margins, in part obfuscating the effects of topography and air mass trajectory.
Consequently, in order to address possible influences of evaporation on δD and δ18O in paleoprecipitation from TP, we analyze n-alkane abundances and distributions along with stable carbon isotope compositions of total organic carbon (δ13CTOC) and palynological data to estimate surface moisture conditions during the early Holocene (ca. 8.7 - 7.5 kyr B.P.) and especially during the 8.2 kyr B.P. climatic event. A period of relatively dry surface conditions from ca. 8.7 to 8.2 kyr B.P., indicated by low values of the 'aquatic index' (Paq ) and by elevated Average Chain Length (ACL) values, in concert with elevated δ13CTOC values, precedes the 8.2 kyr B.P. climatic event. The event itself is characterized by slightly wetter, more humid conditions, as suggested by an increase in Paq values as well as reduced ACL and δ13CTOC values between ca. 8.2 and 7.9 kyr B.P.. In the upper section of the core, a distinct change in paleohydrology becomes. A steep increase in Paq and a decrease in ACL values as well as variations in δ13CTOC from 7.9 kyr B.P. onwards imply considerably elevated surface moisture levels, likely caused by the increased activity of the karstic system of the surrounding mountains. Collectively, the biomarker proxies presented here, reveal a concise picture of changing moisture conditions at TP that is consistent with palynological data and provide detailed paleoenvironmental information for the analysis of δDwax as a paleoprecipitation proxy. The long-term decline in δDwax values characterizes the lower section of the core until ca. 8.2 kyr B.P.. The 8.2 kyr B.P. climatic event itself is connected to two distinct positive hydrogen isotope excursions: a minor shift in δDwax around 8.2 kyr B.P. and a major shift in δDwax between ca. 8.1 and 8.0 kyr B.P.. The upper part of the section shows a progressive trend towards higher δDwax values. With no indication of increased evaporitic conditions at TP during the 8.2 kyr B.P. climatic event, as evident in biomarker proxies and pollen data, we link shifts in δDwax to changes in Mediterranean air mass trajectories supplying precipitation to northeastern Greece, with variations in the relative contributions of northerly derived, D-depleted moisture and southerly-derived, D-enriched moisture. Possible control mechanisms include changes in the influence of the Siberian High and differences in the influence of the African and Asian monsoon circulations on anticyclonic conditions in the Mediterranean region as well as regional inflow of moist air from the Aegean Sea.
Die Wechselwirkung zwischen zwei verschiedenartigen Wellenphänomenen in einer Höhe von ca. 10 bis 100 km, der mittleren Atmosphäre, ist das zentrale Thema der vorliegenden Arbeit. Schwerewellen entstehen durch Oszillationen der Luft in einer stabil geschichteten Atmosphäre. Durch die Vielzahl von Schwerewellen-Paketen, die in der Troposphäre durch Gebirge, Gewitter, Fronten und andere dynamische Prozesse angeregt werden, wird Energie und Impuls in die mittleren Atmosphäre transportiert. Durch den turbulenten Zerfall von brechenden Schwerewellen wird auf die mittlere Strömung eine Kraft ausgeübt, welche im Bereich der Mesopause bei ca. 90 km maximal wird. Daraus resultiert die sogenannte interhemispherische residuelle Zirkulation, die in der Mesosphäre den Sommer- mit dem Winterpol verbindet und die beeindruckend kalte Sommer-Mesopause mit Temperaturen von unter −140°C verursacht. Thermische Gezeiten sind ein weiterer wichtiger Teil in der Dynamik der mittleren Atmosphäre. Sie werden durch die Erwärmung der Tagseite der Erde angeregt und sind globale Schwingungen mit Perioden von 24 Stunden und harmonischen Vielfachen. Mit Wind- und Temperatur-Amplituden von bis zu 50 m/s und 30 K dominieren sie die Tagesvariabilität im Mesopausen-Bereich.
In der Mesosphäre wird die Wechselwirkung zwischen Schwerewellen und thermischen Gezeiten wichtig. Dort wird durch die Gezeitenwinde das Brechen von Schwerewellen zeitlich moduliert und eine periodische Kraft erzeugt, welche auf die Gezeiten rückwirkt. Doch selbst unter Zuhilfenahme modernster Hochleistungsrechner kann in komplexen Zirkulationsmodellen nur ein Bruchteil des turbulenten sowie des Wellen-Spektrums aufgelöst werden. Der Effekt der nichtaufgelösten Skalen, wie Turbulenz und Schwerewellen, muss somit in effizienter Weise parametrisiert werden. Üblicherweise wird in Schwerewellen-Parametrisierungen die horizontale und zeitliche Variation des Hintergrundmediums vernachlässigt. Es entsteht eine vertikale Säule, in der sich stationäre Schwerewellen-Züge instantan nach oben ausbreiten. Es ist jedoch äußerst fraglich, inwieweit eine solche Beschreibung, auf der ein Großteil früherer Untersuchungen basiert, für das Ergründen der Schwerewellen-Gezeiten-Wechselwirkung hinreicht. Für diese Arbeit wurde deswegen das Ziel gesetzt, die Defizite der konventionellen Beschreibung der Schwerewellen-Ausbreitung in realistischen Gezeiten zu quantifizieren.
Die "Ray Tracing"-Methode wird auf die Problemstellung der Schwerewellen-Gezeiten-Wechselwirkung angewendet. In der "Ray Tracing"-Methode werden Schwerewellen-Pakete entlang ihrer Ausbreitungspfade explizit verfolgt und Veränderungen der Schwerewellen-Eigenschaften durch den Einfluss der Hintergrundströmung berücksichtigt. Vom Autor wurde das globale "Ray Tracing"-Modell RAPAGI (RAy PArameterization of Gravity-wave Impacts) entwickelt und mit realistischen Gezeitenfeldern aus dem Zirkulationsmodell HAMMONIA (HAmburg MOdel of the Neutral and Ionized Atmosphere) betrieben. In verschiedenen "Ray Tracing"-Experimenten wird für ein einfaches Schwerewellen-Ensemble gezeigt, wie horizontale Gradienten des Hintergrundmediums sowie dessen Zeitabhängigkeit wesentlichen Einfluss auf die Ausbreitung und Dissipation von Schwerewellen nehmen. Zum einen führt die durch Gezeitenwellen hervorgerufene Transienz zu einer tageszeitlichen Modulation der absoluten Schwerewellen-Frequenz.
Die dadurch induzierten Variationen der horizontalen Phasengeschwindigkeit der Schwerewellen können die anfängliche Phasengeschwindigkeit um bis zu eine Größenordnung übertreffen und folgen dem Verlauf des Hintergrundwindes. Die kritische Filterung von Schwerewellen wird durch diese Modulation abgeschwächt, was im Vergleich zu konventionellen Schwerewellen-Parametrisierungen zu einer im Mittel um 30 % geringeren Kraftwirkung auf die Gezeiten führt. Zum anderen werden durch horizontale Gradienten in der gesamten Hintergrundströmung Schwerewellen-Pakete horizontal abgelenkt. Wellen, die gegen die Hintergrundströmung laufen, werden in der Stratosphäre in die Maxima der Wind-Jets hineingeführt. Durch dieses Verhalten wird analog zum Fermatschen Prinzip der geometrischen Optik die Laufzeit der Schwerewellen in der mittleren Atmosphäre minimiert. Es entsteht eine Fokussierung von Schwerewellen-Feldern, bei gleichzeitiger Zunahme der horizontalen Wellenzahl in den Experimenten im Mittel um ca. 10 %. Dadurch reduziert sich der Schwerewellen-Impulsfluss und die mittlere und ebenfalls die periodische Kraft auf die Hintergrundströmung im Mittel um weitere 20 % bis 30 %. Konventionelle Schwerewellen-Parametrisierungen scheinen somit die Kraftwirkung von brechenden Schwerewellen zu uberschätzen. Aus den Ergebnissen der Arbeit wird klar, dass Schwerewellen-Parametrisierungen nicht "blind" für jede Untersuchung genutzt werden können. Alle Annahmen und Näherungen in Parametrisierungen müssen je nach Zielstellung neu getestet werden.
Physical soil properties feature high spatial variabilities which are known to affect geophysical measurements. However, these variations are not considered in most cases. The challenging task is to quantify the influence of soil heterogeneities on geophysical data. This question is analysed for DC resistivity and GPR measurements which are frequently used for near-surface explorations. To determine the pattern of electric soil properties in situ with the required high spatial resolution, geophysical measuring techniques are methodically enhanced. High-resolution dipole-dipole resistivity measurements are used to determine the electric conductivity distribution of the topsoil. Due to the small electrode separations, the actual electrode geometry has to be considered and an analytic expression for geometric factors is derived instead of assuming point electrodes. Two methods are used to determine soil permittivity with GPR:(i) the coefficient of reflection at the interface air-soil is measured with an air-launched horn antenna, (ii) the velocity of the groundwave is measured with a new setup using two receiver antennas enhancing the lateral resolution from in the best case 0.5 m for standard techniques to approximately 0.1 m with the new technique. With the optimised measuring techniques, the electric properties of sandy soils are determined in the field. Conductivity and permittivity show high spatial variability with correlation lengths of a few decimetres. Geostatistical simulation techniques are used to generate synthetic random media featuring the same statistical properties as in the field. FD calculations are carried out with this media to provide realistic synthetic data of resistivity and GPR measurements. Conductivity variations as determined in the field generate significant variations of simulated Schlumberger sounding curves resulting in uncertainties of the inverted models. Even in pedologically homogeneous sandy soil, moisture pattern and resulting permittivity variations cause strong GPR diffractions as demonstated by FD calculations. This influences the detectability of small objects such as e.g. landmines or of large reflectors as e.g. the groundwater table. Conductivity variations as typical for soils showed to have a minor effect on GPR measurements than variations of permittivity. In summary, geostatistical analysis and simulation provide a powerful tool to simulate geophysical measurements under field conditions including soil heterogeneity which can be used to quantify the uncertainty of field measurements by geologic noise.