ALBERT

Description: Seismological agencies play an important role in seismological research and seismic hazard mitigation by providing source parameters of seismic events (location, magnitude, mechanism), and keeping these data accessible in the long term. The availability of catalogues of seismic source parameters is an important requirement for the evaluation and mitigation of seismic hazards, and the catalogues are particularly valuable to the research community as they provide fundamental long-term data in the geophysical sciences. This work is well motivated by the rising demand for developing more robust and efficient methods for routine source parameter estimation, and ultimately generation of reliable catalogues of seismic source parameters. Specifically, the aim of this work is to develop some methods to determine hypocentre location and temporal evolution of seismic sources based on regional and teleseismic observations more accurately, and investigate the potential of these methods to be integrated in near real-time processing. To achieve this, a location method that considers several events simultaneously and improves the relative location accuracy among nearby events has been provided. This method tries to reduce the biasing effects of the lateral velocity heterogeneities (or equivalently to compensate for limitations and inaccuracies in the assumed velocity model) by calculating a set of timing corrections for each seismic station. The systematic errors introduced into the locations by the inaccuracies in the assumed velocity structure can be corrected without explicitly solving for a velocity model. Application to sets of multiple earthquakes in complex tectonic environments with strongly heterogeneous structure such as subduction zones and plate boundary region demonstrate that this relocation process significantly improves the hypocentre locations compared to standard locations. To meet the computational demands of this location process, a new open-source software package has been developed that allows for efficient relocation of large-scale multiple seismic events using arrival time data. Upon that, a flexible location framework is provided which can be tailored to various application cases on local, regional, and global scales. The latest version of the software distribution including source codes, a user guide, an example data set, and a change history, is freely available to the community. The developed relocation algorithm has been modified slightly and then its performance in a simulated near real-time processing has been evaluated. It has been demonstrated that applying the proposed technique significantly reduces the bias in routine locations and enhance the ability to locate the lower magnitude events using only regional arrival data. Finally, to return to emphasis on global seismic monitoring, an inversion framework has been developed to determine the seismic source time function through direct waveform fitting of teleseismic recordings. The inversion technique can be systematically applied to moderate- size seismic events and has the potential to be performed in near real-time applications. It is exemplified by application to an abnormal seismic event; the 2017 North Korean nuclear explosion. The presented work and application case studies in this dissertation represent the first step in an effort to establish a framework for automatic, routine generation of reliable catalogues of seismic event locations and source time functions.

Description: Seismologische Dienste spielen eine wichtige Rolle in der seismologischen Forschung und Gefährdungsanalyse, indem sie Quellparameter für seismische Ereignisse (Lokalisierung, Stärke, Mechanismus) bereitstellen und diese Daten langfristig verfügbar machen. Die Verfügbarkeit von Katalogen seismischer Quellparameter ist eine wichtige Voraussetzung für die Bewertung der seismischen Gefährdung und Minderung der Auswirkungen von Erdbeben. Die seismischen Kataloge sind für die Forschungsgemeinschaft besonders wertvoll, da sie grundlegende Hintergrundsdaten über lange Zeiträume liefern. Die vorliegende Arbeit ist motiviert durch die steigende Nachfrage nach der Entwicklung robusterer und effizienterer Methoden für die routinemäßige Schätzung von Quellparametern und schließlich die Verbesserung der Qualität und Quantität der seismischen Kataloge. Ziel dieser Arbeit ist es insbesondere Methoden zu entwickeln, um die Lokalisierung von Hypozentren und die zeitliche Entwicklung der Momentfreisetzung bei einem Erdbeben auf der Grundlage regionaler und teleseismischer Beobachtungen genauer zu bestimmen, und das Potenzial dieser Methoden für eine zeitnahe Verarbeitung zu untersuchen. Um dies zu erreichen, wurde eine Lokalisierungsmethode bereitgestellt, die mehrere Ereignisse gleichzeitig berücksichtigt und die relative Genauigkeit der Hypozentren zwischen benachbarten Ereignissen verbessert. Dieses Verfahren versucht, die Verzerrungen durch Geschwindigkeits-heterogenitäten zu reduzieren, oder äquivalent Einschränkungen und Ungenauigkeiten in dem angenommenen Geschwindigkeitsmodell auszugleichen, indem richtungsabhängige Zeitkorrekturen für jede seismische Station berechnet werden. Die systematischen Fehler, die durch die Ungenauigkeiten in der angenommenen Geschwindigkeitsstruktur in die Lokalisierung projiziert werden, können korrigiert werden, ohne explizit nach einem Geschwindigkeitsmodell zu lösen. Die Anwendung auf Ensemble von Erdbeben in tektonischen Situationen mit stark heterogener, komplexer Struktur wie Subduktionszonen und Plattengrenzregionen zeigt, dass diese Relokalisierung die Hypozentren im Vergleich zu Standardlokalisierungen signifikant verbessert. Um den rechnerischen Anforderungen der Lokalisierung gerecht zu werden, wurde ein neues Open-Source-Softwarepaket entwickelt, das eine effiziente Relokalisierung einer großen Anzahl von Erdbeben unter Verwendung von Ankunftszeitdaten ermöglicht. Daraufhin wird ein flexibler Lokalisierungsrahmen bereitgestellt, der an verschiedene Fälle auf lokaler, regionaler und globaler Ebene angepasst werden kann. Die neueste Version der Software einschließlich Quellcodes, Benutzerhandbuch, Beispieldatensatz und Änderungsverlauf ist für die Öffentlichkeit frei verfügbar. Der entwickelte Relokalisierungs-Algorithmus wurde geringfügig verändert und anschließend seine Leistung in einem simulierten Echtzeit Prozess bewertet. Es zeigte sich, dass die Anwendung der vorgeschlagenen Technik die Verzerrung der Hypozentren der Routineauswertung erheblich verringert und die Fähigkeit zum Lokalisieren der Ereignisse mit niedriger Magnitude unter Verwendung nur regionaler Ankunftsdaten verbessert. Um sich wieder auf die globale seismische Überwachung zu konzentrieren, wurde schließlich ein Inversionssystem entwickelt, um die Zeitfunktion der seismischen Quelle durch direkte Wellenformanpassung von teleseismischen Aufzeichnungen zu bestimmen. Die Inversionstechnik kann systematisch auf seismische Ereignisse mittlerer Größe angewendet werden und hat das Potenzial, in nahezu Echtzeitanwendungen ausgeführt zu werden. Dies wird beispielhaft durch die Anwendung auf ein abnormales seismisches Ereignis veranschaulicht: die nordkoreanische Atomexplosion 2017. Die in dieser Dissertation vorgestellten Arbeits- und Anwendungsfallstudien stellen den ersten Schritt dar, um einen Rahmen für die automatische, routinemäßige Erzeugung zuverlässiger Kataloge von seismischen Ereignisorten und Quellzeitfunktion zu schaffen.

Description: This cumulative thesis is concerned with the evolution of geomagnetic activity since the beginning of the 20th century, that is, the time-dependent response of the geomagnetic field to solar forcing. The focus lies on the description of the magnetospheric response field at ground level, which is particularly sensitive to the ring current system, and an interpretation of its variability in terms of the solar wind driving. Thereby, this work contributes to a comprehensive understanding of long-term solar-terrestrial interactions. The common basis of the presented publications is formed by a reanalysis of vector magnetic field measurements from geomagnetic observatories located at low and middle geomagnetic latitudes. In the first two studies, new ring current targeting geomagnetic activity indices are derived, the Annual and Hourly Magnetospheric Currents indices (A/HMC). Compared to existing indices (e.g., the Dst index), they do not only extend the covered period by at least three solar cycles but also constitute a qualitative improvement concerning the absolute index level and the ~11-year solar cycle variability. The analysis of A/HMC shows that (a) the annual geomagnetic activity experiences an interval-dependent trend with an overall linear decline during 1900–2010 of ~5 % (b) the average trend-free activity level amounts to ~28 nT (c) the solar cycle related variability shows amplitudes of ~15–45 nT (d) the activity level for geomagnetically quiet conditions (Kp〈2) lies slightly below 20 nT. The plausibility of the last three points is ensured by comparison to independent estimations either based on magnetic field measurements from LEO satellite missions (since the 1990s) or the modeling of geomagnetic activity from solar wind input (since the 1960s). An independent validation of the longterm trend is problematic mainly because the sensitivity of the locally measured geomagnetic activity depends on geomagnetic latitude. Consequently, A/HMC is neither directly comparable to global geomagnetic activity indices (e.g., the aa index) nor to the partly reconstructed open solar magnetic flux, which requires a homogeneous response of the ground-based measurements to the interplanetary magnetic field and the solar wind speed. The last study combines a consistent, HMC-based identification of geomagnetic storms from 1930–2015 with an analysis of the corresponding spatial (magnetic local time-dependent) disturbance patterns. Amongst others, the disturbances at dawn and dusk, particularly their evolution during the storm recovery phases, are shown to be indicative of the solar wind driving structure (Interplanetary Coronal Mass Ejections vs. Stream or Co-rotating Interaction Regions), which enables a backward-prediction of the storm driver classes. The results indicate that ICME-driven geomagnetic storms have decreased since 1930 which is consistent with the concurrent decrease of HMC. Out of the collection of compiled follow-up studies the inclusion of measurements from high-latitude geomagnetic observatories into the third study’s framework seems most promising at this point.

Description: The Cheb Basin (CZ) is a shallow Neogene intracontinental basin located in the western Eger Rift. The Cheb Basin is characterized by active seismicity and diffuse degassing of mantle-derived CO2 in mofette fields. Within the Cheb Basin, the Hartoušov mofette field shows a daily CO2 flux of 23–97 tons. More than 99% of CO2 released over an area of 0.35 km2. Seismic active periods have been observed in 2000 and 2014 in the Hartoušov mofette field. Due to the active geodynamic processes, the Cheb Basin is considered to be an ideal region for the continental deep biosphere research focussing on the interaction of biological processes with geological processes. To study the influence of CO2 degassing on microbial community in the surface and subsurface environments, two 3-m shallow drillings and a 108.5-m deep scientific drilling were conducted in 2015 and 2016 respectively. Additionally, the fluid retrieved from the deep drilling borehole was also recovered. The different ecosystems were compared regarding their geochemical properties, microbial abundances, and microbial community structures. The geochemistry of the mofette is characterized by low pH, high TOC, and sulfate contents while the subsurface environment shows a neutral pH, and various TOC and sulfate contents in different lithological settings. Striking differences in the microbial community highlight the substantial impact of elevated CO2 concentrations and high saline groundwater on microbial processes. In general, the microorganisms had low abundance in the deep subsurface sediment compared with the shallow mofette. However, within the mofette and the deep subsurface sediment, the abundance of microbes does not show a typical decrease with depth, indicating that the uprising CO2-rich groundwater has a strong influence on the microbial communities via providing sufficient substrate for anaerobic chemolithoautotrophic microorganisms. Illumina MiSeq sequencing of the 16S rRNA genes and multivariate statistics reveals that the pH strongly influences the microbial community composition in the mofette, while the subsurface microbial community is significantly influenced by the groundwater which motivated by the degassing CO2. Acidophilic microorganisms show a much higher relative abundance in the mofette. Meanwhile, the OTUs assigned to family Comamonadaceae are the dominant taxa which characterize the subsurface communities. Additionally, taxa involved in sulfur cycling characterizing the microbial communities in both mofette and CO2 dominated subsurface environments. Another investigated important geo–bio interaction is the influence of the seismic activity. During seismic events, released H2 may serve as the electron donor for microbial hydrogenotrophic processes, such as methanogenesis. To determine whether the seismic events can potentially trigger methanogenesis by the elevated geogenic H2 concentration, we performed laboratory simulation experiments with sediments retrieved from the drillings. The simulation results indicate that after the addition of hydrogen, substantial amounts of methane were produced in incubated mofette sediments and deep subsurface sediments. The methanogenic hydrogenotrophic genera Methanobacterium was highly enriched during the incubation. The modeling of the in-situ observation of the earthquake swarm period in 2000 at the Novy Kostel focal area/Czech Republic and our laboratory simulation experiments reveals a close relation between seismic activities and microbial methane production via earthquake-induced H2 release. We thus conclude that H2 – which is released during seismic activity – can potentially trigger methanogenic activity in the deep subsurface. Based on this conclusion, we further hypothesize that the hydrogenotrophic early life on Earth was boosted by the Late Heavy Bombardment induced seismic activity in approximately 4.2 to 3.8 Ga.

Description: Das Eger-Becken (CZ) ist ein flaches, intrakontinentales neogenes Becken im westlichen Eger-Graben. Das Eger-Becken zeichnet sich durch aktive Seismizität und die diffuse Entgasung von aus dem Mantel stammenden CO2 in Mofettenfeldern aus. Das Mofettenfeld von Hartoušov weist einen täglichen CO2-Fluss von 23-97 Tonnen auf. Mehr als 99% des CO2 werden auf einer Fläche von 0,35 km2 freigesetzt. Im Untersuchungsgebiet wurden in den Jahren 2000 und 2014 seismisch aktive Perioden beobachtet. Aufgrund der aktiven geodynamischen Prozesse gilt das Egerer Becken als ideale Region für die kontinentale Tiefenbiosphärenforschung, die sich auf die Wechselwirkung von biologischen Prozessen mit geologischen Prozessen konzentriert. Zur Untersuchung des Einflusses der CO2-Entgasung auf die mikrobielle Gemeinschaft in der ober- und unterirdischen Umwelt wurden 2015 und 2016 zwei 3 m tiefe Flachbohrungen und eine 108,5 m tiefe wissenschaftliche Bohrung durchgeführt. Zusätzlich wurde auch aus dem Tiefbohrloch Flüssigkeit gewonnen. Die verschiedenen Ökosysteme wurden hinsichtlich ihrer geochemischen Eigenschaften, der mikrobiellen Abundanzen und der mikrobiellen Gemeinschaftsstrukturen verglichen. Die Geochemie der Mofetten zeichnet sich durch einen niedrigen pH-Wert und hohe TOC- und Sulfatgehalte aus, während das unterirdische Milieu einen neutralen pH-Wert und verschiedene TOC- und Sulfatgehalte in unterschiedlichen lithologischen Umgebungen aufweist. Auffällige Unterschiede in der mikrobiellen Gemeinschaft unterstreichen den erheblichen Einfluss erhöhter CO2-Konzentrationen und stark salzhaltigen Grundwassers auf mikrobielle Prozesse. Generell waren die mikrobiellen Abundanzen in dem tiefen Untergrundsediment im Vergleich zur flachen Mofette gering. Innerhalb der Mofette und des tiefen unterirdischen Sediments zeigt die Häufigkeit der Mikroorganismen jedoch keine typische Abnahme mit der Tiefe, was darauf hinweist, dass das aufsteigende CO2-reiche Grundwasser einen starken Einfluss auf die mikrobiellen Gemeinschaften hat, indem es genügend Substrat für anaerobe chemolithoautotrophe Mikroorganismen bietet. Die Illumina-MiSeq-Sequenzierung der 16S rRNA-Gene und die multivariate Statistik zeigen, dass der pH-Wert die Zusammensetzung der mikrobiellen Gemeinschaft in der Mofette signifikant bestimmt, während die unterirdische mikrobielle Gemeinschaft signifikant vom Grundwasser beeinflusst wird, das durch das ausgasende CO2 geprägt ist. Azidophile Mikroorganismen zeigen eine viel höhere relative Abundanz in der Mofette, wohingegen die der Familie Comamonadaceae zugeordneten OTUs die dominierenden Taxa der unterirdischen Gemeinschaften darstellen. Zusätzlich charakterisieren Taxa, die am Schwefelzyklus beteiligt sind, die mikrobiellen Gemeinschaften sowohl in der Mofette als auch in der CO2-dominierten unterirdischen Umwelt. Eine weitere wichtige Untersuchung der Geo-Bio-Interaktion ist der Einfluss der seismischen Aktivität. Während seismischer Ereignisse kann freigesetztes H2 als Elektronendonator für mikrobielle hydrogenotrophe Prozesse, wie z.B. die Methanogenese, dienen. Um zu bestimmen, ob die seismischen Ereignisse durch die erhöhten geogenen H2-Konzentrationen möglicherweise methanogene Prozesse auslösen können, führten wir Laborsimulationsexperimente mit Sedimenten durch, die aus den Bohrungen gewonnen wurden. Die Simulationsexperimente weisen darauf hin, dass nach der Zugabe von Wasserstoff beträchtliche Mengen an Methan in inkubierten Mofettensedimenten und tiefen unterirdischen Sedimenten produziert wurden. Die methanogene hydrogenotrophe Gattung Methanobacterium wurde während der Inkubation stark angereichert. Die Modellierung der in-situ-Beobachtung der Erdbeben-Schwarmzeit im Jahr 2000 im Schwerpunktgebiet Novy Kostel/Tschechische Republik und unsere Laborsimulationsexperimente zeigen einen engen Zusammenhang zwischen seismischen Aktivitäten und der biotischen Methanproduktion durch erdbebeninduzierte H2-Freisetzung. Wir kommen daher zu dem Schluss, dass H2 - dass bei seismischer Aktivität freigesetzt wird - möglicherweise methanogene Aktivität im tiefen Untergrund auslösen kann. Basierend auf dieser Schlussfolgerung gehen wir weiter davon aus, dass das frühe hydrogenotrophe Leben, durch die durch Late Heavy Bombardment induzierte seismische Aktivität in etwa 4,2 bis 3,8 Ga verstärkt wurde.

Description: Rivers have always flooded their floodplains. Over 2.5 billion people worldwide have been affected by flooding in recent decades. The economic damage is also considerable, averaging 100 billion US dollars per year. There is no doubt that damage and other negative effects of floods can be avoided. However, this has a price: financially and politically. Costs and benefits can be estimated through risk assessments. Questions about the location and frequency of floods, about the objects that could be affected and their vulnerability are of importance for flood risk managers, insurance companies and politicians. Thus, both variables and factors from the fields of hydrology and sociol-economics play a role with multi-layered connections. One example are dikes along a river, which on the one hand contain floods, but on the other hand, by narrowing the natural floodplains, accelerate the flood discharge and increase the danger of flooding for the residents downstream. Such larger connections must be included in the assessment of flood risk.

Description: Flüsse haben seit jeher ihre Auen überflutet. In den vergangenen Jahrzehnten waren weltweit über 2,5 Milliarden Menschen durch Hochwasser betroffen. Auch der ökonomische Schaden ist mit durchschnittlich 100 Milliarden US Dollar pro Jahr erheblich. Zweifelsohne können Schäden und andere negative Auswirkungen von Hochwasser vermieden werden. Allerdings hat dies einen Preis: finanziell und politisch. Kosten und Nutzen lassen sich durch Risikobewertungen abschätzen. Dabei werden in der Wasserwirtschaft, von Versicherungen und der Politik Fragen nach dem Ort und der Häufigkeit von Überflutungen, nach den Dingen, die betroffen sein könnten und deren Anfälligkeit untersucht. Somit spielen sowohl Größen und Faktoren aus den Bereichen der Hydrologie und Sozioökonmie mit vielschichtigen Zusammenhängen eine Rolle. Ein anschauliches Beispiel sind Deiche entlang eines Flusses, die einerseits in ihrem Abschnitt Überflutungen eindämmen, andererseits aber durch die Einengung der natürlichen Vorländer den Hochwasserabfluss beschleunigen und die Gefährdung für die Anlieger flussab verschärfen. Solche größeren Zusammenhänge müssen in der Bewertung des Hochwasserrisikos einbezogen werden. In derzeit gängigen Verfahren geht dies mit vereinfachenden Annahmen einher. Risikoabschätzungen sind daher unscharf und mit Unsicherheiten verbunden. Diese Arbeit untersucht den Nutzen und die Möglichkeiten neuer Datensätze für eine Verbesserung der Hochwasserrisikoabschätzung. Es werden neue Methoden und Modelle entwickelt, die die angesprochenen Zusammenhänge stärker berücksichtigen und auch die bestehenden Unsicherheiten der Modellergebnisse beziffern und somit die Verlässlichkeit der getroffenen Aussagen einordnen lassen. Dafür werden Daten zu Hochwasserereignissen aus verschiedenen Quellen erfasst und ausgewertet. Dazu zählen neben Niederschlags-und Durchflussaufzeichnungen an Messstationen beispielsweise auch Bilder aus sozialen Medien, die mit Ortsangaben und Bildinhalten helfen können, die Überflutungsflächen abzugrenzen und Hochwasserschäden zu schätzen. Verfahren des Maschinellen Lernens wurden erfolgreich eingesetzt, um aus vielfältigen Daten, Zusammenhänge zwischen Hochwasser und Auswirkungen zu erkennen, besser zu verstehen und verbesserte Modelle zu entwickeln. Solche Risikomodelle helfen bei der Entwicklung und Bewertung von Strategien zur Minderung des Hochwasserrisikos. Diese Werkzeuge ermöglichen darüber hinaus Einblicke in das Zusammenspiel verschiedener Faktoren sowie Aussagen zu den zu erwartenden Folgen auch von Hochwassern, die das bisher bekannte Ausmaß übersteigen. Diese Arbeit verzeichnet Fortschritte in Bezug auf eine verbesserte Bewertung von Hochwasserrisiken durch die Nutzung vielfältiger Daten aus unterschiedlichen Quellen mit innovativen Verfahren sowie der Weiterentwicklung von Modellen. Das Hochwasserrisiko unterliegt durch wirtschaftliche Entwicklungen und klimatische Veränderungen einem steten Wandel. Um das Wissen über Risiken aktuell zu halten sind robuste, leistungs- und anpassungsfähige Verfahren wie sie in dieser Arbeit vorgestellt werden von zunehmender Bedeutung.

Description: Flooding is a vast problem in many parts of the world, including Europe. It occurs mainly due to extreme weather conditions (e.g. heavy rainfall and snowmelt) and the consequences of flood events can be devastating. Flood risk is mainly defined as a combination of the probability of an event and its potential adverse impacts. Therefore, it covers three major dynamic components: hazard (physical characteristics of a flood event), exposure (people and their physical environment that being exposed to flood), and vulnerability (the elements at risk). Floods are natural phenomena and cannot be fully prevented. However, their risk can be managed and mitigated. For a sound flood risk management and mitigation, a proper risk assessment is needed. First of all, this is attained by a clear understanding of the flood risk dynamics. For instance, human activity may contribute to an increase in flood risk. Anthropogenic climate change causes higher intensity of rainfall and sea level rise and therefore an increase in scale and frequency of the flood events. On the other hand, inappropriate management of risk and structural protection measures may not be very effective for risk reduction. Additionally, due to the growth of number of assets and people within the flood-prone areas, risk increases. To address these issues, the first objective of this thesis is to perform a sensitivity analysis to understand the impacts of changes in each flood risk component on overall risk and further their mutual interactions. A multitude of changes along the risk chain are simulated by regional flood model (RFM) where all processes from atmosphere through catchment and river system to damage mechanisms are taken into consideration. The impacts of changes in risk components are explored by plausible change scenarios for the mesoscale Mulde catchment (sub-basin of the Elbe) in Germany. A proper risk assessment is ensured by the reasonable representation of the real-world flood event. Traditionally, flood risk is assessed by assuming homogeneous return periods of flood peaks throughout the considered catchment. However, in reality, flood events are spatially heterogeneous and therefore traditional assumption misestimates flood risk especially for large regions. In this thesis, two different studies investigate the importance of spatial dependence in large scale flood risk assessment for different spatial scales. In the first one, the “real” spatial dependence of return period of flood damages is represented by continuous risk modelling approach where spatially coherent patterns of hydrological and meteorological controls (i.e. soil moisture and weather patterns) are included. Further the risk estimations under this modelled dependence assumption are compared with two other assumptions on the spatial dependence of return periods of flood damages: complete dependence (homogeneous return periods) and independence (randomly generated heterogeneous return periods) for the Elbe catchment in Germany. The second study represents the “real” spatial dependence by multivariate dependence models. Similar to the first study, the three different assumptions on the spatial dependence of return periods of flood damages are compared, but at national (United Kingdom and Germany) and continental (Europe) scales. Furthermore, the impacts of the different models, tail dependence, and the structural flood protection level on the flood risk under different spatial dependence assumptions are investigated. The outcomes of the sensitivity analysis framework suggest that flood risk can vary dramatically as a result of possible change scenarios. The risk components that have not received much attention (e.g. changes in dike systems and in vulnerability) may mask the influence of climate change that is often investigated component. The results of the spatial dependence research in this thesis further show that the damage under the false assumption of complete dependence is 100 % larger than the damage under the modelled dependence assumption, for the events with return periods greater than approximately 200 years in the Elbe catchment. The complete dependence assumption overestimates the 200-year flood damage, a benchmark indicator for the insurance industry, by 139 %, 188 % and 246 % for the UK, Germany and Europe, respectively. The misestimation of risk under different assumptions can vary from upstream to downstream of the catchment. Besides, tail dependence in the model and flood protection level in the catchments can affect the risk estimation and the differences between different spatial dependence assumptions. In conclusion, the broader consideration of the risk components, which possibly affect the flood risk in a comprehensive way, and the consideration of the spatial dependence of flood return periods are strongly recommended for a better understanding of flood risk and consequently for a sound flood risk management and mitigation.

Description: Hochwasser sind ein großes Problem und treten hauptsächlich aufgrund extremer Wetterbedingungen (z. B. starker Regen und Schneeschmelze) auf. Die Folgen von Hochwasserereignissen können verheerend sein. Das Konzept des Hochwasserrisikos beinhaltet die drei Komponenten: Gefahr, Exposition und Vulnerabilität. Hochwasser sind natürliche Phänomene und können nicht sicher verhindert werden. Das Risiko kann jedoch gesteuert und gemindert werden. Für ein solides Hochwasserrisikomanagement und die Minderung des Risikos ist eine ordnungsgemäße Risikobewertung und ein klares Verständnis der Hochwasserrisikodynamik erforderlich. Beispielsweise verursacht der anthropogene Klimawandel eine höhere Intensität der Niederschläge und einen Anstieg des Meeresspiegels und damit eine Zunahme des Ausmaßes und der Häufigkeit von Hochwasserereignissen. Andererseits können unangemessene strukturelle Schutzmaßnahmen, das Anwachsen von Vermögenswerten und eine steigende Anzahl betroffener Personen in den hochwassergefährdeten Gebieten das Risiko erhöhen. Um diese Probleme zu adressieren, besteht ein Ziel dieser Arbeit aus der Durchführung einer Sensitivitätsanalyse, um die Auswirkungen von Änderungen in jeder Hochwasserrisikokomponente auf das Gesamtrisiko und deren Wechselwirkungen untereinander zu verstehen. Eine angemessene Risikobewertung wird auch durch die korrekte k Darstellung des realen Hochwasserereignisses erreicht. Traditionell wird das Hochwasserrisiko bewertet, indem homogene Wiederkehrintervalle von Hochwasserspitzen im gesamten Einzugsgebiet angenommen werden. In der Realität sind Hochwasserereignisse jedoch räumlich heterogen, weshalb die traditionelle Annahme von Homogenität das Hochwasserrisiko insbesondere für große Einzugsgebiete falsch einschätzt. In dieser Arbeit wird die Bedeutung der räumlichen Abhängigkeit bei der Bewertung des Hochwasserrisikos in großem Maßstab in zwei Studien für verschiedene räumliche Skalen untersucht. In der ersten Untersuchung wird die „reale“ räumliche Abhängigkeit durch einen kontinuierlichen Risikomodellierungsansatz dargestellt. Zusätzlich werden die Risikoabschätzungen unter dieser modellierten Abhängigkeitsannahme mit zwei weiteren Annahmen zur räumlichen Abhängigkeit der Wiederkehrintervalle von Hochwasser verglichen: vollständige Abhängigkeit und Unabhängigkeit für das Elbeeinzugsgebiet in Deutschland. Die zweite Studie repräsentiert die „reale“ räumliche Abhängigkeit durch ein copula-basiertes Abhängigkeitsmodell. In ähnlicher Weise werden die drei verschiedenen Annahmen zur räumlichen Abhängigkeit der Wiederkehrintervalle von Hochwasser auf nationaler und kontinentaler Ebene verglichen. Außerdem wird der Einfluss von „Tail-dependences“ im Modell sowie von Hochwasserschutzmaßnahmen auf die räumliche Abhängigkeit untersucht. Die Ergebnisse dieser Arbeit unter Anwendung des Sensitivitätsanalyse-Frameworks zeigen, dass das Hochwasserrisiko aufgrund möglicher Änderungsszenarien dramatisch variieren kann. Der Einfluss des Klimawandels kann durch Änderungen anderer Risikokomponenten (z. B. Änderungen der Deichsysteme und der Vulnerabilität) überdeckt werden. Die Untersuchung zur räumlichen Abhängigkeit zeigen, dass der Schaden unter der Annahme vollständiger Abhängigkeit für Ereignisse mit Wiederkehrintervalle von mehr als ungefähr 200 Jahren im Elbeeinzugsgebiet 100 % größer als der Schaden unter modellierter Abhängigkeit. Die Annahme vollständiger Abhängigkeit überschätzt den 200-jährigen Hochwasserschaden, einen Referenzindikator für die Versicherungsbranche, um 139 %, 188 % und 246 % für Vereinigte Königreich, Deutschland und Europa. Die Fehleinschätzung des Hochwasserrisikos kann unter verschiedenen Annahmen von Abhängigkeit zwischen Oberlauf und Unterlauf eines Einzugsgebietes stark variieren. Zudem können „Tail-dependences“ im Modell sowie der Hochwasserschutz im Einzugsgebiet die Ergebnisse der Risikoabschätzung, unter verschiedenen Annahmen der räumlichen Abhängigkeit, beeinflussen. Abschließend wird eine umfangreiche Berücksichtigung der Risikokomponenten und insbesondere der räumlichen Abhängigkeit von Wiederkehrintervallen stark empfohlen, um das Hochwasserrisiko und damit dessen Management und Minderung besser verstehen zu können.

Description: This dissertation was carried out as part of the international and interdisciplinary graduate school StRATEGy. This group has set itself the goal of investigating geological processes that take place on different temporal and spatial scales and have shaped the southern central Andes. This study focuses on claystones and carbonates of the Yacoraite Fm. that were deposited between Maastricht and Dan in the Cretaceous Salta Rift Basin. The former rift basin is located in northwest Argentina and is divided into the sub-basins Tres Cruces, Metán-Alemanía and Lomas de Olmedo. The overall motivation for this study was to gain new knowledge about the evolution of marine and lacustrine conditions during the Yacoraite Fm. Deposit in the Tres Cruces and Metán-Alemanía sub-basins. Other important aspects that were examined within the scope of this dissertation are the conversion of organic matter from Yacoraite Fm. into oil and its genetic relationship to selected oils produced and natural oil spills. The results of my study show that the Yacoraite Fm. began to be deposited under marine conditions and that a lacustrine environment developed by the end of the deposition in the Tres Cruces and Metán-Alemanía Basins. In general, the kerogen of Yacoraite Fm. consists mainly of the kerogen types II, III and II / III mixtures. Kerogen type III is mainly found in samples from the Yacoraite Fm., whose TOC values are low. Due to the adsorption of hydrocarbons on the mineral surfaces (mineral matrix effect), the content of type III kerogen with Rock-Eval pyrolysis in these samples could be overestimated. Investigations using organic petrography show that the organic particles of Yacoraite Fm. mainly consist of alginites and some vitrinite-like particles. The pyrolysis GC of the rock samples showed that the Yacoraite Fm. generates low-sulfur oils with a predominantly low-wax, paraffinic-naphthenic-aromatic composition and paraffinic wax-rich oils. Small proportions of paraffinic, low-wax oils and a gas condensate-generating facies are also predicted. Here, too, mineral matrix effects were taken into account, which can lead to a quantitative overestimation of the gas-forming character. The results of an additional 1D tank modeling carried out show that the beginning (10% TR) of the oil genesis took place between ≈10 Ma and ≈4 Ma. Most of the oil (from ≈50% to 65%) was generated prior to the development of structural traps formed during the Plio-Pleistocene Diaguita deformation phase. Only ≈10% of the total oil generated was formed and potentially trapped after the formation of structural traps. Important factors in the risk assessment of this petroleum system, which can determine the small amounts of generated and migrated oil, are the generally low TOC contents and the variable thickness of the Yacoraite Fm. Additional risks are associated with a low density of information about potentially existing reservoir structures and the quality of the overburden.

Description: Diese Dissertation wurde im Rahmen des internationalen und interdisziplinären Graduiertenkollegs StRATEGy durchgeführt. Diese Gruppe hat sich zum Ziel gesetzt, geologische Prozesse zu untersuchen, die auf unterschiedlichen zeitlichen und räumlichen Skalen ablaufen und die südlichen Zentralanden geprägt haben. Diese Studie konzentriert sich auf die Tonsteine und Karbonate der Yacoraite Fm., die zwischen Maastricht und Dan im kreidezeitlichen Salta-Grabenbecken abgelagert wurden. Das ehemalige Riftbecken liegt im Nordwesten Argentiniens und gliedert sich in die Teilbecken Tres Cruces, Metán-Alemanía und Lomas de Olmedo. Die übergreifende Motivation für diese Studie war es, neue Erkenntnisse über die Entwicklung der marinen und lakustrinen Bedingungen während der Ablagerung der Yacoraite Fm. in den Tres Cruces und Metán-Alemanía Sub-Becken zu gewinnen. Weitere wichtige Aspekte, die im Rahmen dieser Dissertation untersucht wurden, sind die Umwandlung von organischer Materie der Yacoraite Fm. in Öl sowie deren genetische Beziehung zu ausgewählten produzierten Ölen und natürlichen Ölaustritten. Die Ergebnisse meiner Studie zeigen, dass die Ablagerung der Yacoraite Fm. unter marinen Bedingungen begann und sich bis zum Ende der Ablagerung in den Tres Cruces- und Metán-Alemanía-Becken ein lakustrines Milieu entwickelte. Im Allgemeinen besteht das Kerogen der Yacoraite Fm. überwiegend aus den Kerogentypen II, III und II/III-Mischungen. Der Kerogentyp III findet sich vor allem in Proben aus der Yacoraite Fm., deren TOC-Werte niedrig sind. Aufgrund der Adsorption von Kohlenwasserstoffen an den Mineraloberflächen (Mineralmatrixeffekt) könnte der Gehalt an Typ-III-Kerogen mit der Rock-Eval-Pyrolyse in diesen Proben überschätzt werden. Untersuchungen mittels organischer Petrographie zeigen, dass die organischen Partikel der Yacoraite Fm. hauptsächlich aus Alginiten und einigen Vitrinit-artigen Partikeln bestehen. Die Pyrolyse-GC der Gesteinsproben zeigte, dass die Yacoraite Fm. schwefelarme Öle mit einer überwiegend Wachs-armen paraffinisch-naphthenisch-aromatischen Zusammensetzung und paraffinische Wachs-reiche Öle generiert. Geringe Anteile paraffinischer, Wachs-armer Öle und einer Gaskondensat-generierenden Fazies werden ebenfalls vorausgesagt. Auch hier wurden Mineralmatrixeffekte berücksichtigt, die zu einer quantitativen Überschätzung des gasbildenden Charakters führen können. Die Ergebnisse einer zusätzlich durchgeführten 1D-Beckenmodellierung zeigen, dass der Beginn (10 %TR) der Ölgenese zwischen ≈10 Ma und ≈4 Ma stattfand. Der größte Teil des Öls (von ≈50 % bis 65 %) wurde vor der Entwicklung struktureller Fallen gebildet, die während der plio-pleistozänen Diaguita Deformatiosphase gebildet wurden. Nur ≈10 % des insgesamt generierten Öls wurde nach der Entstehung struktureller Fallen gebildet und potentiell darin gefangen. Wichtige Faktoren in der Risikobewertung dieses Erdölsystems, welche die geringen Mengen an generiertem und migriertem Öl bestimmen können, stellen die allgemein niedrigen TOC-Gehalte und die variable Mächtigkeit der Yacoraite Fm. dar. Weitere Risiken sind mit einer niedrigen Informationsdichte über potentiell vorhandene Reservoirstrukturen und die Qualität der Deckgesteine verbunden.

Description: Rheology describes the flow of matter under the influence of stress, and - related to solids- it investigates how solids subjected to stresses deform. As the deformation of the Earth’s outer layers, the lithosphere and the crust, is a major focus of rheological studies, rheology in the geosciences describes how strain evolves in rocks of variable composition and temperature under tectonic stresses. It is here where deformation processes shape the form of ocean basins and mountain belts that ultimately result from the complex interplay between lithospheric plate motion and the susceptibility of rocks to the influence of plate-tectonic forces. A rigorous study of the strength of the lithosphere and deformation phenomena thus requires in-depth studies of the rheological characteristics of the involved materials and the temporal framework of deformation processes. This dissertation aims at analyzing the influence of the physical configuration of the lithosphere on the present-day thermal field and the overall rheological characteristics of the lithosphere to better understand variable expressions in the formation of passive continental margins and the behavior of strike-slip fault zones. The main methodological approach chosen is to estimate the present-day thermal field and the strength of the lithosphere by 3-D numerical modeling. The distribution of rock properties is provided by 3-D structural models, which are used as the basis for the thermal and rheological modeling. The structural models are based on geophysical and geological data integration, additionally constrained by 3-D density modeling. More specifically, to decipher the thermal and rheological characteristics of the lithosphere in both oceanic and continental domains, sedimentary basins in the Sea of Marmara (continental transform setting), the SW African passive margin (old oceanic crust), and the Norwegian passive margin (young oceanic crust) were selected for this study. The Sea of Marmara, in northwestern Turkey, is located where the dextral North Anatolian Fault zone (NAFZ) accommodates the westward escape of the Anatolian Plate toward the Aegean. Geophysical observations indicate that the crust is heterogeneous beneath the Marmara basin, but a detailed characterization of the lateral crustal heterogeneities is presented for the first time in this study. Here, I use different gravity datasets and the general non-uniqueness in potential field modeling, to propose three possible end-member scenarios of crustal configuration. The models suggest that pronounced gravitational anomalies in the basin originate from significant density heterogeneities within the crust. The rheological modeling reveals that associated variations in lithospheric strength control the mechanical segmentation of the NAFZ. Importantly, a strong crust that is mechanically coupled to the upper mantle spatially correlates with aseismic patches where the fault bends and changes its strike in response to the presence of high-density lower crustal bodies. Between the bends, mechanically weaker crustal domains that are decoupled from the mantle are characterized by creep. For the passive margins of SW Africa and Norway, two previously published 3-D conductive and lithospheric-scale thermal models were analyzed. These 3-D models differentiate various sedimentary, crustal, and mantle units and integrate different geophysical data, such as seismic observations and the gravity field. Here, the rheological modeling suggests that the present-day lithospheric strength across the oceanic domain is ultimately affected by the age and past thermal and tectonic processes as well as the depth of the thermal lithosphere-asthenosphere boundary, while the configuration of the crystalline crust dominantly controls the rheological behavior of the lithosphere beneath the continental domains of both passive margins. The thermal and rheological models show that the variations of lithospheric strength are fundamentally influenced by the temperature distribution within the lithosphere. Moreover, as the composition of the lithosphere significantly influences the present-day thermal field, it therefore also affects the rheological characteristics of the lithosphere. Overall my studies add to our understanding of regional tectonic deformation processes and the long-term behavior of sedimentary basins; they confirm other analyses that have pointed out that crustal heterogeneities in the continents result in diverse lithospheric thermal characteristics, which in turn results in higher complexity and variations of rheological behavior compared to oceanic domains with a thinner, more homogeneous crust.

Description: Die Rheologie ist die Wissenschaft, die sich mit dem Fließ- und Verformungsverhalten von Materie beschäftigt. Hierzu gehören neben Gasen und Flüssigkeiten vor allem auch Feststoffe, die einer Spannung ausgesetzt sind und einem daraus resultierenden Verformungsprozess unterliegen - entweder unter bruchhaften oder plastischen Bedingungen. In den Geowissenschaften umfasst die Rheologie die kombinierte Analyse tektonischer Spannungen und resultierender Deformationsphänomene in Gesteinen unter unterschiedlichen Temperatur- und Druckbedingungen sowie im Zusammenhang mit physikalischen Eigenschaften der Krusten- und Mantelgesteine. Die Verformung des lithosphärischen Mantels und der Kruste ist ein Schwerpunkt rheologischer Untersuchungen, denn in diesem Zusammenhang bilden sich Ozeanbecken und Gebirgsgürtel, die letztendlich aus dem komplexen Zusammenspiel der Bewegungen lithosphärischer Platten und der unterschiedlichen Deformierbarkeit von Krusten- und Mantelgesteinen unter dem Einfluss plattentektonischer Kräfte resultieren. Eine genaue Untersuchung der Festigkeit der Lithosphäre und der Deformationssphänomene erfordert daher eingehende Studien der rheologischen Eigenschaften der beteiligten Materialien. Vor diesem Hintergrund ist es das Ziel dieser Dissertation, die allgemeinen rheologischen Charakteristika der Lithosphäre in drei verschiedenen geodynamischen Bereichen zu analysieren, um unterschiedlich geprägte passive Kontinentalränder sowie das Verhalten von Transformstörungen innerhalb der Kontinente besser zu verstehen. Der wichtigste methodische Ansatz, der hierfür gewählt wurde, ist die numerische 3D-Modellierung, um eine Abschätzung des gegenwärtigen thermischen Feldes und der Festigkeit der Lithosphäre zu ermöglichen. Die räumliche Verteilung der Gesteinseigenschaften in Kruste und Mantel wird dabei durch 3-D-Strukturmodelle bereitgestellt, die als Grundlage für die thermische und rheologische Modellierung verwendet werden. Die Strukturmodelle basieren auf der Integration geophysikalischer und geologischer Daten, die zusätzlich durch eine 3D-Dichtemodellierung validiert werden. Um die thermischen und rheologischen Eigenschaften der Lithosphäre sowohl im ozeanischen als auch im kontinentalen Bereich zu entschlüsseln, wurden für diese Studie Sedimentbecken im Marmarameer im Bereich der kontinentalen Nordanatolischen Transformstörung sowie im Bereich der passiven Plattenränder von SW-Afrika (alte ozeanische Kruste) und vor Norwegen (junge ozeanische Kruste) ausgewählt. Das Marmarameer im Nordwesten der Türkei befindet sich in einer Region, wo die dextrale Nordanatolische Störung (NAFZ) die westwärts gerichtete Ausweichbewegung der Anatolischen Platte in Richtung Ägäis ermöglicht. Geophysikalische Beobachtungen deuten darauf hin, dass die kontinentale Kruste unter dem Marmara-Becken heterogen ist, allerdings stellt diese Arbeit zum ersten Mal eine detaillierte Charakterisierung dieser lateralen Krustenheterogenitäten vor. Hierzu verwende ich verschiedene Schweredaten und die Potenzialfeldmodellierung, um drei mögliche Szenarien zur Erklärung der Unterschiede im Charakter der Kruste vorzuschlagen. Die Modelle legen nahe, dass ausgeprägte Schwereanomalien im Becken von signifikanten Dichte-Heterogenitäten innerhalb der Kruste hervorgerufen werden. Die rheologische Modellierung zeigt, dass damit verbundene Unterschiede in der Festigkeit der Lithosphäre die mechanische Segmentierung der NAFZ steuern und sich auf seismogene Prozesse auswirken. Demnach korrelieren Krustenbereiche hoher Festigkeit, die mechanisch an den oberen Mantel gekoppelt sind, räumlich mit aseismischen Sektoren in der Region um die Störungszone, in denen sich das Streichen der NAFZ ändert. Zwischen den Bereichen mit den veränderten Streichrichtungen der Störung existieren dagegen mechanisch schwächere Krustenbereiche, die vom Mantel entkoppelt und durch Kriechbewegungen gekennzeichnet sind. Für die passiven Kontinentalränder von SW-Afrika und Norwegen wurden zwei veröffentlichte thermische 3-D-Modelle hinsichtlich des Einflusses der Temperaturverteilung auf die Festigkeit der Lithosphäre analysiert. Diese 3-D-Modelle differenzieren verschiedene Sediment-, Krusten- und Mantelbereiche und integrieren unterschiedliche geophysikalische Daten, wie zum Beispiel seismische Beobachtungen und Schwerefeldmessungen. Hier legt die rheologische Modellierung nahe, dass die derzeitige Lithosphärenfestigkeit im ozeanischen Bereich letztlich durch das Alter und vergangene thermische und tektonische Prozesse sowie die Tiefe der thermischen Grenze zwischen Lithosphäre und Asthenosphäre beeinflusst wird, während die Konfiguration der kristallinen Kruste das rheologische Verhalten der Lithosphäre in den kontinentalen Bereichen der beiden passiven Ränder dominiert. Die thermischen und rheologischen Modelle zeigen, dass die Variationen in der Festigkeit der Lithosphäre grundlegend von der Temperaturverteilung innerhalb der Lithosphäre selbst beeinflusst werden. Dabei steuert die Zusammensetzung der Lithosphäre das heutige thermische Feld entscheidend mit und darüber auch die rheologischen Eigenschaften der Lithosphäre. Diese Ergebnisse tragen somit zu einem besseren Verständnis regionaler tektonischer Deformationsprozesse und der dynamischen Langzeitentwicklung von Sedimentbecken bei; sie bestätigen außerdem frühere Analysen, die bereits darauf hingewiesen haben, dass die Heterogenität der Kruste in den Kontinenten mit unterschiedlichen thermischen Eigenschaften der Lithosphäre einhergeht, welches wiederum zu einer höheren Komplexität und Variabilität des rheologischen Verhaltens im Vergleich zu ozeanischen Gebieten mit einer geringer mächtigen, homogeneren Kruste führt.

Description: In this dissertation, I describe the mechanisms involved in magmatic plumbing system establishment and evolution. Magmatic plumbing systems play a key role in determining volcanic activity style and recognizing its complexities can help in forecasting eruptions, especially within hazardous volcanic systems such as calderas. I explore the mechanisms of dike emplacement and intrusion geometry that shape magmatic plumbing systems beneath caldera-like topographies and how their characteristics relate to precursory activity of a volcanic eruption. For this purpose, I use scaled laboratory models to study the effect of stress field reorientation on a propagating dike induced by caldera topography. I construct these models by using solid gelatin to mimic the elastic properties of the earth's crust with a caldera on the surface. I inject water as the magma analog and track the evolution of the experiments through qualitative (geometry and stress evolution) and quantitative (displacement and strain computation) descriptions. The results show that a vertical dike deviates towards and outside of the caldera-like margin due to stress field reorientation beneath the caldera-like topography. The propagating intrusion forms a circumferential-eruptive dike when the caldera-like size is small, whereas a cone sheet develops beneath the large caldera-like topography. To corroborate the results obtained from the experimental models, this thesis also describes the results of a case study utilizing seismic monitoring data associated with the unrest period of the 2015 phreatic eruption of Lascar volcano. Lascar has a crater with a small-scale caldera-like topography and exhibited long-lasting anomalous evolution of the number of long-period (LP) events preceding the 2015 eruption. I apply seismic techniques to constrain the hypocentral locations of LP events and characterize their spatial distribution, obtaining an image of Lascar's plumbing system. I observe an agreement in shallow hypocentral locations obtained through four different seismic techniques; nevertheless, the cross-correlation technique provides the best results. These results depict a plumbing system with a narrow sub-vertical deep conduit and a shallow hydrothermal system, where most LP events are located. These two regions are connected through an intermediate region of path divergence, whose geometry and orientation likely is influenced by stress reorientation due to topographic effects of the caldera-like crater. Finally, in order to further enhance the interpretations of the previous case study, the seismic data was analyzed in tandem with a complementary multiparametric monitoring dataset. This complementary study confirms that the anomalous LP activity occurred as a sign of unrest in the preparatory phase of the phreatic eruption. In addition, I show how changes observed in other monitored parameters enabled to detect further signs of unrest in the shallow hydrothermal system. Overall, this study demonstrates that detecting complex geometric regions within plumbing systems beneath volcanoes is fundamental to produce an effective forecast of eruptions that from a first view seem to occur without any precursory activity. Furthermore, through the development of this research I show that combining methods that include both observations and models allows one to obtain a more precise interpretation of the volcanic processes.

Description: In dieser Dissertation beschreibe ich die physikalischen Mechanismen, die maßgeblich an der Entstehung und Entwicklung von magmatischen Systemen beteiligt sind. Das magmatische System eines Vulkans spielt eine Schlüsselrolle bei der Bestimmung des vulkanischen Aktivitätsstils und die Erkennung seiner Komplexität kann bei der Vorhersage von Ausbrüchen helfen, insbesondere bei gefährlichen vulkanischen Systemen wie Calderas. Hier erforsche ich deshalb die Mechanismen der Platznahme von Eruptivgängen, sogenannten Dikes und die Geometrie der Intrusionen, welche die magmatischen Systeme unterhalb von Calderen formen, und beleuchte wie ihre Eigenschaften an der Vorläuferaktivität eines Vulkanausbruchs beteiligt sind. Zu diesem Zweck untersuche ich mit Hilfe von skalierten Labormodellen den Effekt, der durch die Caldera-Topographie induzierten Stressfeldneuausrichtung auf einen sich ausbreitenden Dike ausgeübt wird. Ich habe diese Modelle mit fester Gelatine durchgeführt, um die elastischen Eigenschaften der Erdkruste mit einer Caldera an der Oberfläche nachzuahmen. Ich injizierte Wasser als Magma-Analog und verfolgte die Entwicklung der Intrusion. Die Ergebnisse zeigen, dass ein vertikaler Dike aufgrund einer Umorientierung des oberflächennahen Spannungsfeldes unterhalb einer Caldera sowohl in Richtung, als auch nach außerhalb des Caldera-Randes abgelenkt wird. Die propagierende Intrusion bildet somit bei kleinen Calderen einen umlaufend eruptiven Dike, während sich unter großen Calderen bevorzugt sogenannte cone sheets bilden. Um die Ergebnisse aus den vorangegangenen experimentellen Modellen zu bestätigen, werden im Anschluß die Ergebnisse einer Fallstudie beschrieben, in welcher seismologische Monitoringdaten für den Zeitraum des Ausbruchs des Vulkans Lascar im Jahr 2015 herangezogen wurden. Lascar hat einen Krater mit einer kleinräumigen caldera-ähnlichen Topographie und zeigte vor dem phreatischen Ausbruch im Jahr 2015 eine lang anhaltende Periode anomaler Entwicklung bezüglich der Anzahl der täglich registrierten langperiodischen (LP) Erdbeben. Ich habe verschiedene seismische Techniken angewendet, um die Herdregionen (Hypozentren) dieser LP Erdbeben zu lokalisierenund ihre räumliche Verteilung zu charakterisieren, und somit ein Bild von Lascars magmatischen System zu erhalten. Ich beobachte eine Übereinstimmung in den flachen hypozentralen Positionen, die durch vier verschiedene seismische Techniken ermittelt wurden, allerdings lieferte die Kreuzkorrelationstechnik die besten Ergebnisse. Diese Ergebnisse zeigen ein magmatisches System mit einem schmalen tiefen subvertikalen Kanal und einem flachen hydrothermalen System, in dem sich die meisten LP-Ereignisse befinden. Diese beiden Bereiche sind durch einen Zwischenbereich der Pfaddivergenz miteinander verbunden, dessen Geometrie und Ausrichtung wahrscheinlich auf eine oberflächennahe Spannungsumorientierung aufgrund der topographischen Wirkung des caldera-ähnlichen Kraters zurückzuführen ist. Um die Interpretation der vorigen Fallstudie weiter zu vertiefen, habe ich die seismischen Daten in einer weiteren Studie zusammen mit einem ergänzenden multiparametrischen Überwachungsdatensatz analysiert, und wurde darin betstätigt, dass die anomale LP-Aktivität wahrscheinlich als Anzeichen für Unruhe in der Vorbereitungsphase der phreatischen Eruption gedeutet werden kann. Hier zeige ich u.a. auch, wie die Änderungen in den verschiedenen anderen überwachten Parametern weitere Anzeichen von Unruhe im flachen hydrothermalen System erkennen ließen. Insgesamt hat diese Studie gezeigt, wie grundlegend das Erkennen komplexer Geometrien des magmatischen Systems unter Vulkanen ist, um eine effektive Vorhersage von Eruptionen zu erstellen, die auf den ersten Blick ohne Vorläuferaktivität zu erfolgen scheinen. Darüber hinaus zeige ich durch die Entwicklung dieser Forschung, dass die Kombination von Methoden, die sowohl Beobachtungen als auch Modelle umfassen, eine genauere Interpretation der vulkanischen Prozesse ermöglicht.

Description: Near-Earth space represents a significant scientific and technological challenge. Particularly at magnetic low-latitudes, the horizontal magnetic field geometry at the dip equator and its closed field-lines support the existence of a distinct electric current system, abrupt electric field variations and the development of plasma irregularities. Of particular interest are small-scale irregularities associated with equatorial plasma depletions (EPDs). They are responsible for the disruption of trans-ionospheric radio waves used for navigation, communication, and Earth observation. The fast increase of satellite missions makes it imperative to study the near-Earth space, especially the phenomena known to harm space technology or disrupt their signals. EPDs correspond to the large-scale structure (i.e., tens to hundreds of kilometers) of topside F region irregularities commonly known as Spread F. They are observed as depleted-plasma density channels aligned with the ambient magnetic field in the post-sunset low-latitude ionosphere. Although the climatological variability of their occurrence in terms of season, longitude, local time and solar flux is well-known, their day to day variability is not. The sparse observations from ground-based instruments like radars and the few simultaneous measurements of ionospheric parameters by space-based instruments have left gaps in the knowledge of EPDs essential to comprehend their variability. In this dissertation, I profited from the unique observations of the ESA’s Swarm constellation mission launched in November 2013 to tackle three issues that revealed novel and significant results on the current knowledge of EPDs. I used Swarm’s measurements of the electron density, magnetic, and electric fields to answer, (1.) what is the direction of propagation of the electromagnetic energy associated with EPDs?, (2.) what are the spatial and temporal characteristics of the electric currents (field-aligned and diamagnetic currents) related to EPDs, i.e., seasonal/geographical, and local time dependencies?, and (3.) under what conditions does the balance between magnetic and plasma pressure across EPDs occur? The results indicate that: (1.) The electromagnetic energy associated with EPDs presents a preference for interhemispheric flows; that is, the related Poynting flux directs from one magnetic hemisphere to the other and varies with longitude and season. (2.) The field-aligned currents at the edges of EPDs are interhemispheric. They generally close in the hemisphere with the highest Pedersen conductance. Such hemispherical preference presents a seasonal/longitudinal dependence. The diamagnetic currents increase or decrease the magnetic pressure inside EPDs. These two effects rely on variations of the plasma temperature inside the EPDs that depend on longitude and local time. (3.) EPDs present lower or higher plasma pressure than the ambient. For low-pressure EPDs the plasma pressure gradients are mostly dominated by variations of the plasma density so that variations of the temperature are negligible. High-pressure EPDs suggest significant temperature variations with magnitudes of approximately twice the ambient. Since their occurrence is more frequent in the vicinity of the South Atlantic magnetic anomaly, such high temperatures are suggested to be due to particle precipitation. In a broader context, this dissertation shows how dedicated satellite missions with high-resolution capabilities improve the specification of the low-latitude ionospheric electrodynamics and expand knowledge on EPDs which is valuable for current and future communication, navigation, and Earth-observing missions. The contributions of this investigation represent several ’firsts’ in the study of EPDs: (1.) The first observational evidence of interhemispheric electromagnetic energy flux and field-aligned currents. (2.) The first spatial and temporal characterization of EPDs based on their associated field-aligned and diamagnetic currents. (3.) The first evidence of high plasma pressure in regions of depleted plasma density in the ionosphere. These findings provide new insights that promise to advance our current knowledge of not only EPDs but the low-latitude post-sunset ionosphere environment.

Description: To investigate the reliability and stability of spherical harmonic models based on archeo/-paleomagnetic data, 2000 Geomagnetic models were calculated. All models are based on the same data set but with randomized uncertainties. Comparison of these models to the geomagnetic field model gufm1 showed that large scale magnetic field structures up to spherical harmonic degree 4 are stable throughout all models. Through a ranking of all models by comparing the dipole coefficients to gufm1 more realistic uncertainty estimates were derived than the authors of the data provide. The derived uncertainty estimates were used in further modelling, which combines archeo/-paleomagnetic and historical data. The huge difference in data count, accuracy and coverage of these two very different data sources made it necessary to introduce a time dependent spatial damping, which was constructed to constrain the spatial complexity of the model. Finally 501 models were calculated by considering that each data point is a Gaussian random variable, whose mean is the original value and whose standard deviation is its uncertainty. The final model arhimag1k is calculated by taking the mean of the 501 sets of Gauss coefficients. arhimag1k fits different dependent and independent data sets well. It shows an early reverse flux patch at the core-mantle boundary between 1000 AD and 1200 AD at the location of the South Atlantic Anomaly today. Another interesting feature is a high latitude flux patch over Greenland between 1200 and 1400 AD. The dipole moment shows a constant behaviour between 1600 and 1840 AD. In the second part of the thesis 4 new paleointensities from 4 different flows of the island Fogo, which is part of Cape Verde, are presented. The data is fitted well by arhimag1k with the exception of the value at 1663 of 28.3 microtesla, which is approximately 10 microtesla lower than the model suggest.

Description: The current thesis contains the results from two experimental and one modelling study focused on the topic of ductile strain localization in the presence of material heterogeneities. Localization of strain in the high temperature regime is a well known feature of rock deformation occurring in nature at different scales and in a variety of lithologies. Large scale shear zones at the roots of major crustal fault zones are considered responsible for the activity of plate tectonics on our planet. A large number of mechanisms are suggested to be associated with strain softening and nucleation of localization. Among these, the presence of material heterogeneities within homogeneous host rocks is frequently observed in field examples to trigger shear zone development. Despite a number of studies conducted on the topic, the mechanisms controlling initiation and evolution of localization are not fully understood yet. We investigated, experimentally and by means of numerical modelling, phenomenological and microphysical aspects of high temperature strain localization in a homogeneous body containing single and paired inclusions of weaker material. A monomineralic carbonate system composed of Carrara marble (homogeneous, strong matrix) and Solnhofen limestone (weak planar inclusions) is selected for our studies based on its versatility as an experimental material and on the frequent occurrence of carbonate rocks at the core of natural shear zones. To explore the influence of different loading conditions on heterogeneity-induced high temperature shear zones we conducted torsion experiments under constant twist (deformation) rate and constant torque (stress) conditions in a Paterson-type deformation apparatus on hollow cylinders of marble containing single planar inclusions of limestone. At the imposed experimental conditions (900 ◦C temperature and 400 MPa confining pressure) both materials deform plastically and the marble is ≈ 9 times stronger than the limestone. The viscosity contrast between the two materials induces a perturbation of the stress field within the marble matrix at the tip of the planar inclusion. Early on along the deformation path (at bulk shear strains ≈ 0.3), heterogeneous distribution of strain can be observed under both loading conditions and a small area of incipient strain localization is formed at the tip of the weak limestone inclusion. Strongly deformed grains, incipient dynamic recrystallization and a weak crystallographic preferred orientation characterize the marble within an area a few mm in front of the inclusion. As the bulk strain is increased (up to γ ≈ 1), the area of microstructural modification is expanded along the inclusion plane, the texture strengthens and grain size refinement by dynamic recrystallization becomes pervasive. Locally, evidences for coexisting brittle deformation are also observed regardless of the imposed loading conditions. A shear zone is effectively formed within the deforming Carrara marble, its geometry controlled by the plane containing the thin plate of limestone. Thorough microstructural and textural analysis, however, do not reveal substantial differences in the mechanisms or magnitude of strain localization at the different loading conditions. We conclude that, in the presence of material heterogeneities capable of inducing strain softening, the imposed loading conditions do not affect ductile localization in its nucleating and transient stages. As the ultimate goal of experimental rock deformation is the extrapolation of results to geologically relevant time and space scales, we developed 2D numerical models reproducing (and benchmarked to) our experimental results. Our cm-scaled models have been implemented with a first-order strain-dependent softening law to reproduce the effect of rheological weakening in the deforming material. We successfully reproduced the local stress concentration at the inclusion tips and the strain localization initiated in the marble matrix. The heterogeneous distribution of strain and its evolution with imposed bulk deformation (i.e. the shape and extent of the nucleating shear zone) are observed to depend on the degree of softening imposed to the deforming matrix. When a second (artificial) softening step is introduced at elevated bulk strains in the model, the formation of a secondary high strain layer is observed at the core of the initial shear zone, analogous to the development of ultramylonite bands in high strain natural shear zones. Our results do not only reproduce the nucleation and transient evolution of a heterogeneity-induced high temperature shear zone with high accuracy, but also confirm the importance of introducing reliable softening laws capable of mimicking strain weakening to numerical models of crustal scale ductile processes. Material heterogeneities inducing strain localization in the field are often consisting of brittle precursors (joints and fractures). More generally, the interaction of brittle and ductile deformation mechanisms and its effect on the localization of strain have been a key topic in the structural geology community for a long time. The positive feedback between (micro)fracturing and ductile strain localization is a well recognized effect in a number of field examples. We experimentally investigated the influence of brittle deformation on the initiation and evolution of high temperature shear zones in a strong matrix containing pairs of weak material heterogeneities. Our Carrara marble-Solnhofen limestone inclusions system was tested in triaxial compression under constant strain rate and high temperature (900 ◦C) conditions in a Paterson deformation apparatus. The inclusion pairs were arranged in non-overlapping step-over geometries of either compressional or extensional nature. Experimental runs were conducted at different confining pressures (30, 50, 100 and 300 MPa) to induce various amounts of brittle deformation within the marble matrix. At low confinement (30 and 50 MPa) abundant brittle deformation is observed in all configurations, but the spatial distribution of cracks is dependent on the kinematics of the step-over region: concentrated along the shearing plane between the inclusions in the extensional samples, or broadly distributed around the inclusions but outside the step-over region in the compressional configuration. Accordingly, brittle-assisted ductile processes tend to localize deformation along the inclusions plane in the extensional geometry or to distribute widely across large areas of the matrix in the compressional step-over. At pressures of 100 and 300 MPa fracturing is mostly suppressed in both configurations and strain is accommodated almost entirely by viscous creep. In extensional samples this leads to progressive de-localization with increasing confinement. Our results show that, while ductile localization of strain is indeed more efficient where assisted by brittle processes, these latter are only effective if themselves heterogeneously distributed, ultimately a function of the local stress perturbations.

Description: Die vorliegende Doktorarbeit umfasst Ergebnisse von zwei experimentellen und einer Modellierungsstudie. Diese befassen sich mit der Lokalisierung von duktilen Verformungen, hervorgerufen durch unterschiedliche Materialeigenschaften. Die Lokalisierung von Verformungen im Hochtemperaturbereich in unterschiedlichen Maßstäben und in einer Vielzahl von Lithologien ist ein bekanntes Merkmal der natürlichen Gesteinsdeformationen. So wird beispielsweise die Aktivität der Plattentektonik unseres Planeten durch weiträumige Scherzonen am Grund dieser Plattengrenzen verantwortlich gemacht. Dabei wird eine große Anzahl von Mechanismen mit der durch die Verformung hervorgerufenen Materialermüdung und der Ausbildung der Lokalisierung in Verbindung gebracht. Dabei wird unter diesen Mechanismen das Vorhandensein von Materialheterogenitäten innerhalb eines Gesteins häufig als Auslöser für die Ausbildung von Scherzonen beobachtet. Obwohl bereits Studien zu diesem Thema durchgeführt wurden, sind die kontrollierenden Mechanismen, die für die Initiierung und Entwicklung der Lokalisierung zuständig sind, bis heute nicht vollumfänglich verstanden. Aus diesem Grund wurden im Rahmen der vorgelegten Dissertation phänomenologische und mikrophysikalische Aspekte der Lokalisierung von Verformungen im Hochtemperaturbereich in einem homogenen Gesteinskörper, der mit einfachen und gepaarten Inklusionen aus weicheren Material versehen wurde, experimentell und unter Hilfenahme von numerischen Modellen untersucht. Da Karbonatgesteine häufig am Ursprung natürlicher Scherzonen auftreten und diese für ihre Vielseitigkeit als Experimentiermaterial bekannt sind, wurde ein monomineralisches Karbonatsystem, bestehend aus Carrara Marmor (homogene, starke Matrix) und Solnhofen Kalkstein (schwache Inklusionen) als zu untersuchendes Probenmaterial für diese Studie gewählt. Um den Einfluss unterschiedlicher Deformationsbedingungen auf die, durch die Materialheterogenität hervorgerufenen Scherzonen im Hochtemperaturbereich zu untersuchen, wurden Torsionsexperimente bei konstanter Torsionsrate und konstantem Drehmoment in einer Paterson-Deformationsapparatur an hohlen Carrara Marmorzylindern mit einer ebenen Inklusion bestehend aus Kalkstein durchgeführt. Unter den vorgegebenen Randbedingungen (Temperatur = 900 ◦C, Manteldruck = 400 MPa) verformten sich beide Materialien plastisch, wobei die Festigkeit des Marmors in etwa dem neunfachen der Kalksteinfestigkeit entspricht. An der Spitze der ebenen Kalksteininklusion wird durch den Viskositätskontrast der beiden Materialien dadurch eine Störung des Spannungsfeldes in der Marmormatrix hervorgerufen. In der frühen Phase der Deformation (Scherverformung γ ≈ 0.3) kann eine heterogene Verteilung der Verformungen in der gesamten Probe bei beiden Experimenttypen beobachtet werden. Zusätzlich beginnt sich an der Spitze der schwächeren Kalksteininklusion ein Bereich mit lokaler Verformung in der Marmormatrix auszubilden. Dieser ist durch stark deformierte Mineralkörner, beginnende dynamische Rekristallisation und einer schwach ausgeprägten kristallographisch bevorzugten Ausrichtung innerhalb einer Fläche weniger Millimeter charakterisiert. Mit ansteigender Gesamtverformung (γ ≈ 1) erweitert sich die Fläche der mikrostrukturellen Modifikationen entlang der Inklusionsebene. Zusätzlich konnten eine verfestigte Textur und eine Verfeinerung der Korngröße aufgrund dynamischer Rekristallisation beobachtet werden. Lokale Anzeichen für eine gleichzeitige spröde Verformung konnten, unabhängig von den Deformationsbedingungen, ebenfalls festgestellt werden. Innerhalb des deformierten Marmors bildete sich eine Scherzone aus, deren Geometrie maßgeblich durch die Ebene der Kalksteininklusion kontrolliert wird. Sorgfältig durchgeführte mikrostrukturelle Analysen zeigten jedoch keine wesentlichen Unterschiede der Mechanismen oder dem Ausmaß der Lokalisierung der Verformung bei unterschiedlichen Deformationsbedingungen. Daraus lässt sich schließen, dass bei dem Vorhandensein von Materialheterogenitäten, welche eine verformungsbedingte Materialermüdung hervorrufen können, die verwendeten Deformationsbedingungen keinen Einfluss auf die Lokalisierung duktiler Deformation in ihrer Entstehung und übergangsphasen haben. Da für gewöhnlich die Durchführung von Deformationsexperimenten auf die Extrapolation der gewonnenen Ergebnisse auf geologische Zeiträume abzielt, wurden zwei- dimensionale numerische Modelle entwickelt, welche in der Lage sind die aufgenommenen experimentellen Daten zu reproduzieren und zu bewerten. Diese, auf dem Zentimetermaßstab skalierten Modelle wurden mit einem verformungsbasierten Ermüdungsgesetz erster Ordnung umgesetzt, um den Effekt der rheologischen Materialermüdung nachzubilden. Die lokale Spannungskonzentration an der Spitze der Inklusionen und die in der Marmormatrix initiierten Lokalisierung der Verformung konnten mit diesen Modellen erfolgreich reproduziert werden. Dabei wurde festgestellt, dass die heterogene Verteilung der Verformung und deren Entwicklung mit zunehmender Gesamtverformung (z.B. Form und Umfang der sich ausbildenden Scherzone) abhängig vom Grad der Ermüdung der deformierten Matrix ist. Analog zu der Entwicklung von Ultramylonitbändern in natürlichen Scherzonen mit hoher Verformung, konnte bei der Einführung eines zweites (künstlichen) Ermüdungsschritts bei erhöhter Gesamtverformung, die Ausbildung einer zweiten Schicht mit großer Verformung am Kern der initialen Scherzone beobachtet werden. Mit den gewonnenen Ergebnissen der numerischen Simulationen wurde nicht nur die Ausbildung und transiente Entwicklung, der durch die Materialheterogenität hervorgerufenen Scherzone im Hochtemperaturbereich mit großer Genauigkeit reproduziert. Auch wurde die Wichtigkeit verlässliche Ermüdungsgesetze aufzustellen, die in der Lage sind die durch die Verformung hervorgerufenen Materialermüdung im geologischen Maßstab unter Hinzunahme von numerischen Modellen nachzuahmen, bestätigt. Die durch Materialheterogenitäten hervorgerufene Lokalisierung von Verformungen in der Natur bestehen häufig aus spröden Vorläufern, wie beispielsweise Klüften und Rissen. Die Interaktion von spröden und duktilen Deformationsmechanismen im Allgemeinen und ihr Effekt auf die Lokalisierung von Verformungen sind seit langem Schlüsselthema auf dem Gebiet der Strukturgeologie. Die Kopplung von spröden Bruchprozessen mit der Lokalisierung duktiler Verformungen ist oft Gegenstand der Untersuchung in einer Vielzahl von Feldstudien. Daher wurde experimentell der Einfluss von spröder Deformation auf die Initiierung und Entwicklung von Hochtemperatur Scherzonen in einer starken Matrix mit schwächeren Materialheterogenitäten untersucht. Dafür wurden Carrara Marmor-Solnhofen Kalksteininklusions-Systeme unter triaxialen Bedingungen bei konstanter axialer Verformungsrate und hoher Temperatur (T = 900 ◦C) in einer PatersonApparatur deformiert. Dabei wurden die Inklusionen in einer übereinander liegenden, aber nicht überlappenden Geometrie so angeordnet, dass sich bei axialer Probendeformation entweder Kompression oder Dehnung in der Fläche zwischen den Inklusionen einstellt. Die Versuche wurden bei verschiedenen Manteldrücken (P = 30, 50, 100 und 300 MPa) durchgeführt, um unterschiedliche Beträge an spröder Deformation in der Marmormatrix hervorzurufen. Bei geringen Manteldrücken (P = 30 und 50 MPa) kann ein hoher Anteil an spröder Deformation in allen Probenkonfigurationen beobachtet werden. Allerdings ist die räumliche Verteilung der Risse abhängig von der Kinematik der sich übereinanderliegenden Inklusionen. Bei der Dehnungskonfiguration sind die Risse entlang der Scherfläche zwischen den Inklusionen konzentriert, während sie bei der Kompressionskonfiguration außer halb der Inklusionen weit verteilt sind. Dementsprechend neigen duktile, durch spröde unterstützte Prozesse Deformationen entlang der Inklusionsebene bei Dehnung zu lokalisieren oder sich über weite Fläche der Matrix bei Kompression zu verteilen. Bei Manteldrücken von 100 und 300 MPa ist die Risserzeugung in beiden Konfigurationen weitestgehend unterdrückt und die Verformung wird fast ausschließlich durch viskoses Kriechen akkommodiert. Mit ansteigendem Manteldruck führt das bei Proben der Dehnungskonfiguration zu fortschreitender De-Lokalisierung. Die Ergebnisse zeigen, dass die Lokalisierung von duktilen Verformungen effizienter ist, wenn diese durch spröde Bruchprozesse assistiert werden. Allerdings sind diese spröden Prozesse nur dann effektiv, wenn sie heterogen verteilt sind, was letztendlich eine Funktion der lokalen Spannungsstörungen ist.

Description: Salt pans also termed playas are common landscape features of hydrologically closed basins in arid and semiarid zones, where evaporation significantly exceeds the local precipitation. The analysis and monitoring of salt pan environments is important for the evaluation of current and future impact of these landscape features. Locally, salt pans have importance for the ecosystem, wildlife and human health, and through dust emissions they influence the climate on regional and global scales. Increasing economic exploitation of these environments in the last years, e.g. by brine extraction for raw materials, as well as climate change severely affect the water, material and energy balance of these systems. Optical remote sensing has the potential to characterise salt pan environments and to increase the understanding of processes in playa basins, as well as to assess wider impacts and feedbacks that exist between climate forcing and human intervention in their regions. Remote sensing techniques can provide information for extensive regions on a high temporal basis compared to traditional field samples and ground observations. Specifically, for salt pans that are often challenging to study because of their large size, remote location, and limited accessibility due to missing infrastructure and ephemeral flooding. Furthermore, the availability of current and upcoming hyperspectral remote sensing data opened the opportunity for the analyses of the complex reflectance signatures that relate to the mineralogical mixtures found in the salt pan sediments. However, these new advances in sensor technology, as well as increased data availability currently have not been fully explored for the study of salt pan environments. The potential of new sensors needs to be assessed and state of the art methods need to be adapted and improved to provide reliable information for in depth analysis of processes and characterisation of the recent condition, as well as to support long-term monitoring and to evaluate environmental impacts of changing climate and anthropogenic activity. This thesis provides an assessment of the capabilities of optical remote sensing for the study of salt pan environments that combines the information of hyperspectral data with the increased temporal coverage of multispectral observations for a more complete understanding of spatial and temporal complexity of salt pan environments using the Omongwa salt pan located in the southwest Kalahari as a test site. In particular, hyperspectral data are used for unmixing of the mineralogical surface composition, spectral feature-based modelling for quantification of main crust components, as well as time-series based classification of multispectral data for the assessment of the long-term dynamic and the analysis of the seasonal process regime. The results show that the surface of the Omongwa pan can be categorized into three major crust types based on diagnostic absorption features and mineralogical ground truth data. The mineralogical crust iv types can be related to different zones of surface dynamic as well as pan morphology that influences brine flow during the pan inundation and desiccation cycles. Using current hyperspectral imagery, as well as simulated data of upcoming sensors, robust quantification of the gypsum component could be derived. For the test site the results further indicate that the crust dynamic is mainly driven by flooding events in the wet season, but it is also influenced by temperature and aeolian activity in the dry season. Overall, the scientific outcomes show that optical remote sensing can provide a wide range of information helpful for the study of salt pan environments. The thesis also highlights that remote sensing approaches are most relevant, when they are adapted to the specific site conditions and research scenario and that upcoming sensors will increase the potential for mineralogical, sedimentological and geomorphological analysis, and will improve the monitoring capabilities with increased data availability.

Description: The Earth's inner magnetosphere is a very dynamic system, mostly driven by the external solar wind forcing exerted upon the magnetic field of our planet. Disturbances in the solar wind, such as coronal mass ejections and co-rotating interaction regions, cause geomagnetic storms, which lead to prominent changes in charged particle populations of the inner magnetosphere - the plasmasphere, ring current, and radiation belts. Satellites operating in the regions of elevated energetic and relativistic electron fluxes can be damaged by deep dielectric or surface charging during severe space weather events. Predicting the dynamics of the charged particles and mitigating their effects on the infrastructure is of particular importance, due to our increasing reliance on space technologies. The dynamics of particles in the plasmasphere, ring current, and radiation belts are strongly coupled by means of collisions and collisionless interactions with electromagnetic fields induced by the motion of charged particles. Multidimensional numerical models simplify the treatment of transport, acceleration, and loss processes of these particles, and allow us to predict how the near-Earth space environment responds to solar storms. The models inevitably rely on a number of simplifications and assumptions that affect model accuracy and complicate the interpretation of the results. In this dissertation, we quantify the processes that control electron dynamics in the inner magnetosphere, paying particular attention to the uncertainties of the employed numerical codes and tools. We use a set of convenient analytical solutions for advection and diffusion equations to test the accuracy and stability of the four-dimensional Versatile Electron Radiation Belt (VERB-4D) code. We show that numerical schemes implemented in the code converge to the analytical solutions and that the VERB-4D code demonstrates stable behavior independent of the assumed time step. The order of the numerical scheme for the convection equation is demonstrated to affect results of ring current and radiation belt simulations, and it is crucially important to use high-order numerical schemes to decrease numerical errors in the model. Using the thoroughly tested VERB-4D code, we model the dynamics of the ring current electrons during the 17 March 2013 storm. The discrepancies between the model and observations above 4.5 Earth's radii can be explained by uncertainties in the outer boundary conditions. Simulation results indicate that the electrons were transported from the geostationary orbit towards the Earth by the global-scale electric and magnetic fields. We investigate how simulation results depend on the input models and parameters. The model is shown to be particularly sensitive to the global electric field and electron lifetimes below 4.5 Earth's radii. The effects of radial diffusion and subauroral polarization streams are also quantified. We developed a data-assimilative code that blends together a convection model of energetic electron transport and loss and Van Allen Probes satellite data by means of the Kalman filter. We show that the Kalman filter can correct model uncertainties in the convection electric field, electron lifetimes, and boundary conditions. It is also demonstrated how the innovation vector - the difference between observations and model prediction - can be used to identify physical processes missing in the model of energetic electron dynamics. We computed radial profiles of phase space density of ultrarelativistic electrons, using Van Allen Probes measurements. We analyze the shape of the profiles during geomagnetically quiet and disturbed times and show that the formation of new local minimums in the radial profiles coincides with the ground observations of electromagnetic ion-cyclotron (EMIC) waves. This correlation indicates that EMIC waves are responsible for the loss of ultrarelativistic electrons from the heart of the outer radiation belt into the Earth's atmosphere.

Description: Die innere Magnetosphäre der Erde ist ein sehr dynamisches System, das hauptsächlich vom äußeren Sonnenwind beeinflusst wird, der auf das Magnetfeld unseres Planeten einwirkt. Störungen im Sonnenwind, wie z.B. koronale Massenauswürfe und sogenannte Korotierende Wechselwirkungsbereiche, verursachen geomagnetische Stürme, die zu deutlichen Veränderungen der Populationen geladener Teilchen in der inneren Magnetosphäre führen - Plasmasphäre, Ringstrom und Strahlungsgürtel. Satelliten, die in Regionen mit erhöhten energetischen und relativistischen Elektronenflüssen betrieben werden, können durch tiefe dielektrische Ladung oder Oberflächenladungen bei schweren Weltraumwetterereignissen beschädigt werden. Die Vorhersage der Dynamik der geladenen Teilchen und die Abschwächung ihrer Auswirkungen auf die Infrastruktur sind heutzutage von besonderer Bedeutung, insbesondere aufgrund unserer zunehmenden Abhängigkeit von Weltraumtechnologien. Die Dynamik von Teilchen in der Plasmasphäre, des Ringstrom und in den Strahlungsgürteln sind durch Kollisionen und kollisionsfreie Wechselwirkungen mit elektromagnetischen Feldern, die durch die Bewegung geladener Teilchen induziert werden, stark gekoppelt. Mehrdimensionale numerische Modelle vereinfachen die Betrachtung von Transport-, Beschleunigungs- und Verlustprozessen dieser Partikel und ermöglichen es uns, vorherzusagen, wie die erdnahe Weltraumumgebung auf Sonnenstürme reagiert. Die Modelle beruhen zwangsläufig auf einer Reihe von Vereinfachungen und Voraussetzungen, die sich auf die Modellgenauigkeit auswirken und die Interpretation der Ergebnisse erschweren. In dieser Dissertation quantifizieren wir die Prozesse, die die Dynamik der Elektronen in der inneren Magnetosphäre steuern. Dabei richten wir den Fokus insbesondere auch auf die Unsicherheiten der verwendeten numerischen Codes. Wir verwenden eine Reihe praktischer analytischer Lösungen für Advektions- und Diffusionsgleichungen, um die Genauigkeit und Stabilität des 4-dimensionalen ''Versatile Electron Radiation Belt'' Codes (VERB-4D Code) zu testen. Wir zeigen, dass die im Code implementierten numerischen Schemata zu den analytischen Lösungen konvergieren und der Code sich unabhängig vom angenommenen Zeitschritt stabil verhält. Wir demonstrieren, wie die Genauigkeit des numerischen Schemas für die Konvektionsgleichung die Ergebnisse von Ringstrom- und Strahlungsgürtelsimulationen beeinflussen kann, und dass es von entscheidender Beteutung ist, numerische Schemata höherer Ordnung zu verwenden, um numerische Fehler im Modell zu reduzieren. Mit dem ausführlich getesteten VERB-4D Code modellieren wir die Dynamik der Ringstromelektronen während des Sturms vom 17. März 2013. Wir zeigen, dass die Diskrepanzen zwischen dem Modell und Beobachtungen oberhalb von 4.5 Erdradien durch Unsicherheiten in den äußeren Randbedingungen erklärt werden können und dass die Elektronen durch die globalen elektrischen und magnetischen Felder von der geostationäre Umlaufbahn zur Erde transportiert wurden. Wir untersuchen weiterhin, wie die Simulationsergebnisse von den Eingabemodellen und Parametern abhängen. Wir zeigen, dass das Modell besonders empfindlich für das globale elektrische Feld und die Lebensdauer der Elektronen unterhalb von 4.5 Erdradien ist. Außerdem quantifizieren wir auch die Auswirkungen von radialer Diffusion und subauroralen Polarisationsströmen. Wir haben einen datenassimilativen Code entwickelt, der mithilfe des Kalman-Filters ein Konvektionsmodell für den Transport und den Verlust energetischer Elektronen mit den Satellitendaten der Van Allen Probes kombiniert. Wir zeigen, dass die Verwendung eines Kalman-Filters Modellunsicherheiten im elektrischen Konvektionsfeld, in der Lebensdauer der Elektronen und in den Randbedingungen korrigieren kann. Weiterhin zeigen wir, wie der Innovationsvektor - die Differenz zwischen Beobachtungen und Modellvorhersagen - verwendet werden kann, um physikalische Prozesse zu identifizieren, die im Modell der Dynamik der energetischen Elektronen fehlen. Außerdem berechnen wir radiale Profile der Phasenraumdichte ultrarelativistischer Elektronen mithilfe von Van Allen Probes-Messungen. Wir analysieren die Form der Profile und zeigen, dass die Entstehung neuer lokaler Minima in den radialen Profilen mit den Bodenbeobachtungen von EMIC-Wellen übereinstimmt. Diese Korrelation legt nahe, dass EMIC-Wellen für den Verlust ultrarelativistischer Elektronen vom Herzen des äußeren Strahlungsgürtels in die Erdatmosphäre verantwortlich sind.

Description: The plasmasphere is a dynamic region of cold, dense plasma surrounding the Earth. Its shape and size are highly susceptible to variations in solar and geomagnetic conditions. Having an accurate model of plasma density in the plasmasphere is important for GNSS navigation and for predicting hazardous effects of radiation in space on spacecraft. The distribution of cold plasma and its dynamic dependence on solar wind and geomagnetic conditions remain, however, poorly quantified. Existing empirical models of plasma density tend to be oversimplified as they are based on statistical averages over static parameters. Understanding the global dynamics of the plasmasphere using observations from space remains a challenge, as existing density measurements are sparse and limited to locations where satellites can provide in-situ observations. In this dissertation, we demonstrate how such sparse electron density measurements can be used to reconstruct the global electron density distribution in the plasmasphere and capture its dynamic dependence on solar wind and geomagnetic conditions. First, we develop an automated algorithm to determine the electron density from in-situ measurements of the electric field on the Van Allen Probes spacecraft. In particular, we design a neural network to infer the upper hybrid resonance frequency from the dynamic spectrograms obtained with the Electric and Magnetic Field Instrument Suite and Integrated Science (EMFISIS) instrumentation suite, which is then used to calculate the electron number density. The developed Neural-network-based Upper hybrid Resonance Determination (NURD) algorithm is applied to more than four years of EMFISIS measurements to produce the publicly available electron density data set. We utilize the obtained electron density data set to develop a new global model of plasma density by employing a neural network-based modeling approach. In addition to the location, the model takes the time history of geomagnetic indices and location as inputs, and produces electron density in the equatorial plane as an output. It is extensively validated using in-situ density measurements from the Van Allen Probes mission, and also by comparing the predicted global evolution of the plasmasphere with the global IMAGE EUV images of He+ distribution. The model successfully reproduces erosion of the plasmasphere on the night side as well as plume formation and evolution, and agrees well with data. The performance of neural networks strongly depends on the availability of training data, which is limited during intervals of high geomagnetic activity. In order to provide reliable density predictions during such intervals, we can employ physics-based modeling. We develop a new approach for optimally combining the neural network- and physics-based models of the plasmasphere by means of data assimilation. The developed approach utilizes advantages of both neural network- and physics-based modeling and produces reliable global plasma density reconstructions for quiet, disturbed, and extreme geomagnetic conditions. Finally, we extend the developed machine learning-based tools and apply them to another important problem in the field of space weather, the prediction of the geomagnetic index Kp. The Kp index is one of the most widely used indicators for space weather alerts and serves as input to various models, such as for the thermosphere, the radiation belts and the plasmasphere. It is therefore crucial to predict the Kp index accurately. Previous work in this area has mostly employed artificial neural networks to nowcast and make short-term predictions of Kp, basing their inferences on the recent history of Kp and solar wind measurements at L1. We analyze how the performance of neural networks compares to other machine learning algorithms for nowcasting and forecasting Kp for up to 12 hours ahead. Additionally, we investigate several machine learning and information theory methods for selecting the optimal inputs to a predictive model of Kp. The developed tools for feature selection can also be applied to other problems in space physics in order to reduce the input dimensionality and identify the most important drivers. Research outlined in this dissertation clearly demonstrates that machine learning tools can be used to develop empirical models from sparse data and also can be used to understand the underlying physical processes. Combining machine learning, physics-based modeling and data assimilation allows us to develop novel methods benefiting from these different approaches.

Description: Die Plasmasphäre ist eine die Erde umgebende dynamische Region aus kaltem, dichtem Plasma. Ihre Form und Größe sind sehr anfällig für Schwankungen der solaren und geomagnetischen Bedingungen. Ein präzises Modell der Plasmadichte in der Plasmasphäre ist wichtig für die GNSS-Navigation und für die Vorhersage gefährlicher Auswirkungen der kosmischen Strahlung auf Raumfahrzeuge. Die Verteilung des kalten Plasmas und seine dynamische Abhängigkeit vom Sonnenwind und den geomagnetischen Bedingungen sind jedoch nach wie vor nur unzureichend quantifiziert. Bestehende empirische Modelle der Plasmadichte sind in der Regel zu stark vereinfacht, da sie auf statistischen Durchschnittswerten statischer Parameter basieren. Das Verständnis der globalen Dynamik der Plasmasphäre anhand von Beobachtungen aus dem Weltraum bleibt eine Herausforderung, da vorhandene Dichtemessungen spärlich sind und sich auf Orte beschränken, an denen Satelliten In-situ-Beobachtungen liefern können. In dieser Dissertation zeigen wir, wie solche spärlichen Elektronendichtemessungen verwendet werden können, um die globale Elektronendichteverteilung in der Plasmasphäre zu rekonstruieren und ihre dynamische Abhängigkeit vom Sonnenwind und den geomagnetischen Bedingungen zu erfassen. Zunächst entwickeln wir einen automatisierten Algorithmus zur Bestimmung der Elektronendichte aus In-situ-Messungen des elektrischen Feldes der Van Allen Probes Raumsonden. Insbesondere entwerfen wir ein neuronales Netzwerk, um die obere Hybridresonanzfrequenz aus den dynamischen Spektrogrammen abzuleiten, die wir durch die Instrumentensuite „Electric and Magnetic Field Instrument Suite“ (EMFISIS) erhielten, welche dann zur Berechnung der Elektronenzahldichte verwendet wird. Der entwickelte „Neural-network-based Upper Hybrid Resonance Determination“ (NURD)-Algorithmus wird auf mehr als vier Jahre der EMFISIS-Messungen angewendet, um den öffentlich verfügbaren Elektronendichte-Datensatz zu erstellen. Wir verwenden den erhaltenen Elektronendichte-Datensatz, um ein neues globales Modell der Plasmadichte zu entwickeln, indem wir einen auf einem neuronalen Netzwerk basierenden Modellierungsansatz verwenden. Zusätzlich zum Ort nimmt das Modell den zeitlichen Verlauf der geomagnetischen Indizes und des Ortes als Eingabe und erzeugt als Ausgabe die Elektronendichte in der äquatorialebene. Dies wird ausführlich anhand von In-situ-Dichtemessungen der Van Allen Probes-Mission und durch den Vergleich der vom Modell vorhergesagten globalen Entwicklung der Plasmasphäre mit den globalen IMAGE EUV-Bildern der He+ -Verteilung validiert. Das Modell reproduziert erfolgreich die Erosion der Plasmasphäre auf der Nachtseite sowie die Bildung und Entwicklung von Fahnen und stimmt gut mit den Daten überein. Die Leistung neuronaler Netze hängt stark von der Verfügbarkeit von Trainingsdaten ab, die für Intervalle hoher geomagnetischer Aktivität nur spärlich vorhanden sind. Um zuverlässige Dichtevorhersagen während solcher Intervalle zu liefern, können wir eine physikalische Modellierung verwenden. Wir entwickeln einen neuen Ansatz zur optimalen Kombination der neuronalen Netzwerk- und physikbasierenden Modelle der Plasmasphäre mittels Datenassimilation. Der entwickelte Ansatz nutzt sowohl die Vorteile neuronaler Netze als auch die physikalischen Modellierung und liefert zuverlässige Rekonstruktionen der globalen Plasmadichte für ruhige, gestörte und extreme geomagnetische Bedingungen. Schließlich erweitern wir die entwickelten auf maschinellem Lernen basierten Werkzeuge und wenden sie auf ein weiteres wichtiges Problem im Bereich des Weltraumwetters an, die Vorhersage des geomagnetischen Index Kp. Der Kp-Index ist einer der am häufigsten verwendeten Indikatoren für Weltraumwetterwarnungen und dient als Eingabe für verschiedene Modelle, z.B. für die Thermosphäre, die Strahlungsgürtel und die Plasmasphäre. Es ist daher wichtig, den Kp-Index genau vorherzusagen. Frühere Arbeiten in diesem Bereich verwendeten hauptsächlich künstliche neuronale Netze, um Kurzzeit-Kp-Vorhersagen zu treffen, wobei deren Schlussfolgerungen auf der jüngsten Vergangenheit von Kp- und Sonnenwindmessungen am L1-Punkt beruhten. Wir analysieren, wie sich die Leistung neuronaler Netze im Vergleich zu anderen Algorithmen für maschinelles Lernen verhält, um kurz- und längerfristige Kp-Voraussagen von bis zu 12 Stunden treffen zu können. Zusätzlich untersuchen wir verschiedene Methoden des maschinellen Lernens und der Informationstheorie zur Auswahl der optimalen Eingaben für ein Vorhersagemodell von Kp. Die entwickelten Werkzeuge zur Merkmalsauswahl können auch auf andere Probleme in der Weltraumphysik angewendet werden, um die Eingabedimensionalität zu reduzieren und die wichtigsten Treiber zu identifizieren. Die in dieser Dissertation skizzierten Untersuchungen zeigen deutlich, dass Werkzeuge für maschinelles Lernen sowohl zur Entwicklung empirischer Modelle aus spärlichen Daten als auch zum Verstehen zugrunde liegender physikalischer Prozesse genutzt werden können. Die Kombination von maschinellem Lernen, physikbasierter Modellierung und Datenassimilation ermöglicht es uns, kombinierte Methoden zu entwickeln, die von unterschiedlichen Ansätzen profitieren.

Description: This thesis investigates how the permafrost microbiota responds to global warming. In detail, the constraints behind methane production in thawing permafrost were linked to methanogenic activity, abundance and composition. Furthermore, this thesis offers new insights into microbial adaptions to the changing environmental conditions during global warming. This was assesed by investigating the potential ecological relevant functions encoded by plasmid DNA within the permafrost microbiota. Permafrost of both interglacial and glacial origin spanning the Holocene to the late Pleistocene, including Eemian, were studied during long-term thaw incubations. Furthermore, several permafrost cores of different stratigraphy, soil type and vegetation cover were used to target the main constraints behind methane production during short-term thaw simulations. Short- and long-term incubations simulating thaw with and without the addition of substrate were combined with activity measurements, amplicon and metagenomic sequencing of permanently frozen and seasonally thawed active layer. Combined, it allowed to address the following questions. i) What constraints methane production when permafrost thaws and how is this linked to methanogenic activity, abundance and composition? ii) How does the methanogenic community composition change during long-term thawing conditions? iii) Which potential ecological relevant functions are encoded by plasmid DNA in active layer soils? The major outcomes of this thesis are as follows. i) Methane production from permafrost after long-term thaw simulation was found to be constrained mainly by the abundance of methanogens and the archaeal community composition. Deposits formed during periods of warmer temperatures and increased precipitation, (here represented by deposits from the Late Pleistocene of both interstadial and interglacial periods) were found to respond strongest to thawing conditions and to contain an archaeal community dominated by methanogenic archaea (40% and 100% of all detected archaea). Methanogenic population size and carbon density were identified as main predictors for potential methane production in thawing permafrost in short-term incubations when substrate was sufficiently available. ii) Besides determining the methanogenic activity after long-term thaw, the paleoenvironmental conditions were also found to influence the response of the methanogenic community composition. Substantial shifts within methanogenic community structure and a drop in diversity were observed in deposits formed during warmer periods, but not in deposits from stadials, when colder and drier conditions occurred. Overall, a shift towards a dominance of hydrogenotrophic methanogens was observed in all samples, except for the oldest interglacial deposits from the Eemian, which displayed a potential dominance of acetoclastic methanogens. The Eemian, which is discussed to serve as an analogue to current climate conditions, contained highly active methanogenic communities. However, all potential limitation of methane production after permafrost thaw, it means methanogenic community structure, methanogenic population size, and substrate pool might be overcome after permafrost had thawed on the long-term. iii) Enrichments with soil from the seasonally thawed active layer revealed that its plasmid DNA (‘metaplasmidome’) carries stress-response genes. In particular it encoded antibiotic resistance genes, heavy metal resistance genes, cold shock proteins and genes encoding UV-protection. Those are functions that are directly involved in the adaptation of microbial communities to stresses in polar environments. It was further found that metaplasmidomes from the Siberian active layer originate mainly from Gammaproteobacteria. By applying enrichment cultures followed by plasmid DNA extraction it was possible to obtain a higher average contigs length and significantly higher recovery of plasmid sequences than from extracting plasmid sequences from metagenomes. The approach of analyzing ‘metaplasmidomes’ established in this thesis is therefore suitable for studying the ecological role of plasmids in polar environments in general. This thesis emphasizes that including microbial community dynamics have the potential to improve permafrost-carbon projections. Microbially mediated methane release from permafrost environments may significantly impact future climate change. This thesis identified drivers of methanogenic composition, abundance and activity in thawing permafrost landscapes. Finally, this thesis underlines the importance to study how the current warming Arctic affects microbial communities in order to gain more insight into microbial response and adaptation strategies.

Description: Diese Dissertation untersucht die Reaktion der Permafrost-Mikrobiota auf die globale Erwärmung. Im Detail wurden mögliche Faktoren, die die Methanproduktion in tauendem Permafrost einschränken, im Zusammenhang methanogener Aktivität, Abundanz und Gemeinschaftszusammensetzung untersucht. Darüber hinaus bietet diese Dissertation neue Einblicke in mikrobielle Anpassungen an die sich ändernden Umweltbedingungen während der globalen Erwärmung. Dies wurde durch Untersuchung der potenziell ökologisch relevanten Funktionen bewertet, die von Plasmid-DNA innerhalb der Permafrost-Mikrobiota codiert werden. Permafrost, der seinen Ursprung in den Interglazialen und Glazialen aus dem Holozän bis zum späten Pleistozän, einschließlich des Eem, hat, wurde in Langzeit-Tau-Inkubationen untersucht. Darüber hinaus wurden mehrere Permafrostkerne mit unterschiedlicher Stratigraphie, Vegetationsbedeckung und unterschiedlichem Bodentyp verwendet, um die Faktoren, die die Methanproduktion während kurzfristiger Auftausimulationen bestimmen, zu ermitteln. Kurz- und Langzeitinkubationen, die das Auftauen mit und ohne Zugabe von Substrat in Kombination mit Aktivitätsmessungen, Amplikon- und Metagenom-Sequenzierung von permanent gefrorenem und saisonal aufgetautem Boden simulieren, ermöglichten die Beantwortung folgender Fragen: i) Welche Faktoren hemmen die Methanproduktion beim Auftauen des Permafrosts und wie hängt dies mit der Aktivität, Abundanz und Zusammensetzung methanogener Organismen zusammen? ii) Wie verändert sich die Gemeinschaftszusammensetzung methanogener Organismen unter langfristigen Auftaubedingungen? iii) Welche potenziell ökologisch relevanten Funktionen werden von Plasmid-DNA in saisonal getauten Böden kodiert? Die wichtigsten Ergebnisse dieser Arbeit können wie folgt zusammengefasst werden. i) Die Methanproduktion in langfristig getautem Permafrost wird hauptsächlich durch die Anzahl der methanogenen Archaeen und ihrem Anteil innerhalb der Archaeen bestimmt. Ablagerungen, die in wärmeren Perioden mit erhöhtem Niederschlag gebildet wurden, reagierten am stärksten auf das Tauen und enthielten eine von Methanogenen dominierte Archaeen-Gemeinschaft. In Kurzzeitinkubationen mit ausreichender Verfügbarkeit von Substrat wurden die Populationsgröße der methanogenen Organismen und die Kohlenstoffdichte als Hauptprädiktoren für die potenzielle Methanproduktion beim Auftauen von Permafrost identifiziert. ii) Auch die paläoökologischen Bedingungen beeinflussen die Reaktion der methanogenen Gemeinschaft und Aktivität, wenn Permafrost taut. Es wurden erhebliche Verschiebungen innerhalb der Gemeinschaftsstruktur und ein Rückgang der Diversität in Ablagerungen beobachtet, die in wärmeren Perioden gebildet wurden, jedoch nicht bei Ablagerungen aus kälteren und trockeneren Perioden. Insgesamt wurde in allen Proben eine Verschiebung hin zu einer Dominanz von hydrogenotrophen Methanogenen beobachtet, mit Ausnahme der ältesten interglazialen Ablagerungen aus dem Eem, die eine potenzielle Dominanz von acetoklastischen Methanogenen aufwiesen. Das Eem, das als Analogon zu den aktuellen Klimabedingungen diskutiert wird, enthielt hochaktive methanogene Gemeinschaften. iii) Anreicherungen aus Boden der saisonalen Auftauschicht zeigten, dass die enthaltene Plasmid-DNA („Metaplasmidom“) Stress-Reaktions-Gene trägt. Insbesondere codierte die Plasmid-DNA Antibiotikaresistenzgene, Schwermetallresistenzgene, Kälteschock-proteine und Gene, für den UV-Schutz, also Funktionen, die direkt an der Anpassung mikrobieller Gemeinschaften an Stress in polaren Umgebungen beteiligt sind. Weiterhin stammten die Metaplasmidome der saisonalen Auftauschicht Sibiriens hauptsächlich von Gammaproteobakterien. Durch die Anreicherung von Kulturen, gefolgt von einer Extraktion der Plasmid-DNA, war es möglich, eine höhere durchschnittliche Contig-Länge und eine signifikant höhere Wiederherstellung von Plasmidsequenzen zu erhalten als durch Extrahieren von Plasmidsequenzen aus Metagenomen. Der in dieser Arbeit etablierte Ansatz zur Analyse von „Metaplasmidomen“ ist ein geeigneter Ansatz zur Untersuchung der ökologischen Rolle von Plasmiden in polaren Regionen insgesamt. Diese Dissertation hebt hervor, wie wichtig es ist, die Abundanz, Zusammensetzung und Funktionen der mikrobiellen Gemeinschaft in Permafrost-Kohlenstoff-Projektionen einzubeziehen, und zwar nicht nur, da die mikrobiell vermittelte Methanfreisetzung aus Permafrostablagerungen das Potenzial hat, den zukünftigen Klimawandel erheblich zu beeinflussen. Vielmehr wurden in dieser Arbeit Abhängigkeiten methanogener Gemeinschaftsstrukturen, Abundanz und Aktivität identifiziert. Abschließend verdeutlicht diese Arbeit, wie wichtig es ist zu untersuchen, wie sich die derzeitige Erwärmung der Arktis auf mikrobielle Gemeinschaften auswirkt, um Einblicke in mikrobielle Reaktions- und Anpassungsstrategien zu erhalten.