TY - THES A1 - Mester, Benedikt T1 - Modeling flood-induced human displacement risk under global change T1 - Modellierung des hochwasserbedingten Risikos der Vertreibung von Menschen unter globalen Veränderungen N2 - Extreme flooding displaces an average of 12 million people every year. Marginalized populations in low-income countries are in particular at high risk, but also industrialized countries are susceptible to displacement and its inherent societal impacts. The risk of being displaced results from a complex interaction of flood hazard, population exposed in the floodplains, and socio-economic vulnerability. Ongoing global warming changes the intensity, frequency, and duration of flood hazards, undermining existing protection measures. Meanwhile, settlements in attractive yet hazardous flood-prone areas have led to a higher degree of population exposure. Finally, the vulnerability to displacement is altered by demographic and social change, shifting economic power, urbanization, and technological development. These risk components have been investigated intensively in the context of loss of life and economic damage, however, only little is known about the risk of displacement under global change. This thesis aims to improve our understanding of flood-induced displacement risk under global climate change and socio-economic change. This objective is tackled by addressing the following three research questions. First, by focusing on the choice of input data, how well can a global flood modeling chain reproduce flood hazards of historic events that lead to displacement? Second, what are the socio-economic characteristics that shape the vulnerability to displacement? Finally, to what degree has climate change potentially contributed to recent flood-induced displacement events? To answer the first question, a global flood modeling chain is evaluated by comparing simulated flood extent with satellite-derived inundation information for eight major flood events. A focus is set on the sensitivity to different combinations of the underlying climate reanalysis datasets and global hydrological models which serve as an input for the global hydraulic model. An evaluation scheme of performance scores shows that simulated flood extent is mostly overestimated without the consideration of flood protection and only for a few events dependent on the choice of global hydrological models. Results are more sensitive to the underlying climate forcing, with two datasets differing substantially from a third one. In contrast, the incorporation of flood protection standards results in an underestimation of flood extent, pointing to potential deficiencies in the protection level estimates or the flood frequency distribution within the modeling chain. Following the analysis of a physical flood hazard model, the socio-economic drivers of vulnerability to displacement are investigated in the next step. For this purpose, a satellite- based, global collection of flood footprints is linked with two disaster inventories to match societal impacts with the corresponding flood hazard. For each event the number of affected population, assets, and critical infrastructure, as well as socio-economic indicators are computed. The resulting datasets are made publicly available and contain 335 displacement events and 695 mortality/damage events. Based on this new data product, event-specific displacement vulnerabilities are determined and multiple (national) dependencies with the socio-economic predictors are derived. The results suggest that economic prosperity only partially shapes vulnerability to displacement; urbanization, infant mortality rate, the share of elderly, population density and critical infrastructure exhibit a stronger functional relationship, suggesting that higher levels of development are generally associated with lower vulnerability. Besides examining the contextual drivers of vulnerability, the role of climate change in the context of human displacement is also being explored. An impact attribution approach is applied on the example of Cyclone Idai and associated extreme coastal flooding in Mozambique. A combination of coastal flood modeling and satellite imagery is used to construct factual and counterfactual flood events. This storyline-type attribution method allows investigating the isolated or combined effects of sea level rise and the intensification of cyclone wind speeds on coastal flooding. The results suggest that displacement risk has increased by 3.1 to 3.5% due to the total effects of climate change on coastal flooding, with the effects of increasing wind speed being the dominant factor. In conclusion, this thesis highlights the potentials and challenges of modeling flood- induced displacement risk. While this work explores the sensitivity of global flood modeling to the choice of input data, new questions arise on how to effectively improve the reproduction of flood return periods and the representation of protection levels. It is also demonstrated that disentangling displacement vulnerabilities is feasible, with the results providing useful information for risk assessments, effective humanitarian aid, and disaster relief. The impact attribution study is a first step in assessing the effects of global warming on displacement risk, leading to new research challenges, e.g., coupling fluvial and coastal flood models or the attribution of other hazard types and displacement events. This thesis is one of the first to address flood-induced displacement risk from a global perspective. The findings motivate for further development of the global flood modeling chain to improve our understanding of displacement vulnerability and the effects of global warming. N2 - Durch extreme Überschwemmungen werden jedes Jahr durchschnittlich 12 Millionen Menschen vertrieben. Vor allem marginalisierte Bevölkerungsgruppen in Ländern mit niedrigem Einkommen sind stark gefährdet, aber auch Industrieländer sind anfällig für Vertreibungen und die damit verbundenen gesellschaftlichen Auswirkungen. Das Risiko der Vertreibung ergibt sich aus einer komplexen Wechselwirkung zwischen der Hochwassergefahr, der Exposition der in den Überschwemmungsgebieten lebenden Bevölkerung und der sozioökonomischen Vulnerabilität. Die fortschreitende globale Erderwärmung verändert die Intensität, Häufigkeit und Dauer von Hochwassergefahren und untergräbt die bestehenden Schutzmaßnahmen. Gleichzeitig hat die Besiedlung attraktiver, aber gefährdeter Überschwemmungsgebiete zu einem höheren Maß an Exposition der Bevölkerung geführt. Schließlich wird die Vulnerabilität für Vertreibungen durch den demografischen und sozialen Wandel, die Verlagerung der Wirtschaftskräfte, die Urbanisierung und die technologische Entwicklung verändert. Diese Risikokomponenten wurden im Zusammenhang mit dem Verlust von Menschenleben und wirtschaftlichen Schäden intensiv untersucht, über das Risiko der Vertreibung im Rahmen des globalen Wandels ist jedoch nur wenig bekannt. Diese Arbeit zielt darauf ab, unser Verständnis des durch Überschwemmungen verursachten Vertreibungsrisikos unter dem Einfluss des globalen Klimawandels und des sozioökonomischen Wandels zu verbessern. Dieses Ziel wird durch die Beantwortung der folgenden drei Forschungsfragen erreicht. Erstens: Wie gut kann eine globale Hochwassermodellierungskette die Hochwassergefahren historischer Ereignisse, die zu Vertreibung geführt haben, reproduzieren, wobei ein Fokus auf die Wahl der Eingangsdaten gelegt wird? Zweitens: Welches sind die sozioökonomischen Merkmale, die die Vulnerabilität für Vertreibung beeinflussen? Und schließlich, inwieweit hat der Klimawandel möglicherweise zu den jüngsten hochwasserbedingten Vertreibungsereignissen beigetragen? Zur Beantwortung der ersten Frage wird eine globale Hochwassermodellierungskette durch den Vergleich der simulierten Überschwemmungsfläche mit satellitengestützten Überschwemmungsdaten für acht große Hochwasserereignisse überprüft. Der Schwerpunkt liegt dabei auf der Sensitivität gegenüber verschiedenen Kombinationen der zugrunde liegenden Klimareanalysedatensätzen und globalen hydrologischen Modellen, die als Input für das globale Hydraulikmodell dienen. Ein Bewertungsschema von Leistungsindikatoren zeigt, dass die simulierte Überschwemmungsfläche ohne Berücksichtigung des Hochwasserschutzes meist überschätzt wird und nur bei wenigen Ereignissen von der Wahl der globalen hydrologischen Modelle abhängt. Die Ergebnisse sind empfindlicher gegenüber dem zugrunde liegenden Climate Forcing, wobei sich zwei Datensätze erheblich von einem dritten unterscheiden. Im Gegensatz dazu führt die Einbeziehung von Hochwasserschutznormen zu einer Unterschätzung der Überschwemmungsfläche, was auf mögliche Mängel bei der Schätzung des Schutzniveaus oder der Hochwasserhäufigkeitsverteilung innerhalb der Modellierungskette hinweist. Nach der Analyse des physikalischen Hochwassergefahrenmodells werden in einem nächsten Schritt die sozioökonomischen Triebkräfte für die Vulnerabilität für Vertreibungen untersucht. Zu diesem Zweck wird eine satellitengestützte, globale Sammlung von Hochwasseüberschwemmungsflächen mit zwei Katastrophendatenbänken verknüpft, um die gesellschaftlichen Auswirkungen mit der entsprechenden Hochwassergefahr zusammenzuführen. Für jedes Ereignis werden die Anzahl der betroffenen Menschen, Vermögenswerte und kritischen Infrastrukturen sowie sozioökonomische Indikatoren berechnet. Die daraus resultierenden Datensätze werden öffentlich zugänglich gemacht und enthalten 335 Vertreibungsereignisse und 695 Todesopfer-/Schadensereignisse. Auf der Grundlage dieses neuen Datenprodukts werden ereignisspezifische Vertreibungsvulnerabilitäten bestimmt und vielfältige (nationale) Abhängigkeiten mit den sozioökonomischen Prädiktoren abgeleitet. Die Ergebnisse deuten darauf hin, dass wirtschaftlicher Wohlstand nur teilweise die Anfälligkeit für Vertreibungen beeinflusst; Urbanisierung, Kindersterblichkeitsrate, der Anteil älterer Menschen, Bevölkerungsdichte und kritische Infrastrukturen weisen eine stärkere funktionale Beziehung auf, was den Schluss zulässt, dass ein höheres Entwicklungsniveau im Allgemeinen mit einer geringeren Vulnerabilität verbunden ist. Neben der Untersuchung der kontextabhängigen Faktoren der Vulnerabilität wird auch die Rolle des Klimawandels im Zusammenhang mit der Vertreibung von Menschen untersucht. Am Beispiel des Zyklons Idai und den damit verbundenen extremen Küstenüberschwemmungen in Mosambik wird ein Ansatz zur Attribution der Auswirkungen angewandt. Eine Kombination aus Küstenüberflutungsmodellierung und Satellitenbildern wird verwendet, um faktische und kontrafaktische Überschwemmungsereignisse zu konstruieren. Diese Storyline-artige Attributionsmethode ermöglicht die Untersuchung der isolierten oder kombinierten Auswirkungen des Meeresspiegelanstiegs und der Intensivierung der Windgeschwindigkeiten von Zyklonen auf die Küstenüberflutung. Die Ergebnisse deuten darauf hin, dass das Vertreibungsrisiko durch die Gesamtwirkung des Klimawandels auf Küstenüberschwemmungen um 3,1 bis 3,5 % gestiegen ist, wobei die Auswirkungen der zunehmenden Windgeschwindigkeit der dominierende Faktor sind. Zusammenfassend zeigt diese Arbeit die Potentiale und Herausforderungen der Modellierung von hochwasserbedingten Vertreibungsrisiken auf. Während diese Arbeit die Sensitivität der globalen Hochwassermodellierung in Bezug auf die Wahl der Eingabedaten untersucht, ergeben sich neue Fragen, wie die Reproduktion von Wiederkehrintervallen und die Darstellung von Schutzniveaus effektiv verbessert werden kann. Die Ergebnisse liefern nützliche Informationen für Risikobewertungen, effektive humanitäre Hilfe und Katastrophenhilfe. Die Studie zur Auswirkungs-Attribution ist ein erster Schritt zur Bewertung der Effekte der globalen Erwärmung auf das Vertreibungsrisiko und führt zu neuen Forschungsherausforderungen, z. B. zur Kopplung von Fluss- und Küstenhochwassermodellen oder zur Untersuchung anderer Gefahrenarten. Diese Arbeit ist eine der ersten, die das durch Überschwemmungen verursachte Vertreibungsrisiko aus einer globalen Perspektive heraus betrachtet. Die Ergebnisse motivieren dazu, die globale Hochwassermodellierungskette weiterzuentwickeln, um unser Verständnis der Vertreibungsvulnerabilität und der Auswirkungen der globalen Erderwärmung zu vertiefen. KW - displacement KW - flooding KW - remote sensing KW - vulnerability KW - global flood model KW - Vertreibung KW - Überschwemmungen KW - Fernerkundung KW - Vulnerabilität KW - globales Überschwemmungsmodell Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-609293 ER - TY - THES A1 - Manu, Evans T1 - Hydrogeochemical characterization of water resources in the Pra Basin (Ghana) for quality assessment and water management T1 - Hydrogeochemische Charakterisierung der Wasserressourcen im Pra-Becken (Ghana) zur Qualitätsbewertung und Wasserbewirtschaftung BT - field observations and geochemical modelling BT - Feldbeobachtungen und geochemische Modellierung N2 - Watershed management requires an understanding of key hydrochemical processes. The Pra Basin is one of the five major river basins in Ghana with a population of over 4.2 million people. Currently, water resources management faces challenges due to surface water pollution caused by the unregulated release of untreated household and industrial waste into aquatic ecosystems and illegal mining activities. This has increased the need for groundwater as the most reliable water supply. Our understanding of groundwater recharge mechanisms and chemical evolution in the basin has been inadequate, making effective management difficult. Therefore, the main objective of this work is to gain insight into the processes that determine the hydrogeochemical evolution of groundwater quality in the Pra Basin. The combined use of stable isotope, hydrochemistry, and water level data provides the basis for conceptualizing the chemical evolution of groundwater in the Pra Basin. For this purpose, the origin and evaporation rates of water infiltrating into the unsaturated zone were evaluated. In addition, Chloride Mass Balance (CMB) and Water Table Fluctuations (WTF) were considered to quantify groundwater recharge for the basin. Indices such as water quality index (WQI), sodium adsorption ratio (SAR), Wilcox diagram, and salinity (USSL) were used in this study to determine the quality of the resource for use as drinking water and for irrigation purposes. Due to the heterogeneity of the hydrochemical data, the statistical techniques of hierarchical cluster and factor analysis were applied to subdivide the data according to their spatial correlation. A conceptual hydrogeochemical model was developed and subsequently validated by applying combinatorial inverse and reaction pathway-based geochemical models to determine plausible mineral assemblages that control the chemical composition of the groundwater. The interactions between water and rock determine the groundwater quality in the Pra Basin. The results underline that the groundwater is of good quality and can be used for drinking water and irrigation purposes. It was demonstrated that there is a large groundwater potential to meet the entire Pra Basin’s current and future water demands. The main recharge area was identified as the northern zone, while the southern zone is the discharge area. The predominant influence of weathering of silicate minerals plays a key role in the chemical evolution of the groundwater. The work presented here provides fundamental insights into the hydrochemistry of the Pra Basin and provides data important to water managers for informed decision-making in planning and allocating water resources for various purposes. A novel inverse modelling approach was used in this study to identify different mineral compositions that determine the chemical evolution of groundwater in the Pra Basin. This modelling technique has the potential to simulate the composition of groundwater at the basin scale with large hydrochemical heterogeneity, using average water composition to represent established spatial groupings of water chemistry. N2 - Die Bewirtschaftung von Wassereinzugsgebieten erfordert ein Verständnis der wichtigsten hydrochemischen Prozesse. Das Pra-Becken ist eines der fünf großen Flusseinzugsgebiete Ghanas mit einer Bevölkerung von über 4,2 Millionen Menschen. Die Bewirtschaftung der Wasserressourcen wird derzeit durch die Verschmutzung der Oberflächengewässer erschwert, die durch die unkontrollierte Einleitung von unbehandelten Haushalts- und Industrieabfällen in die aquatischen Ökosysteme und durch illegale Bergbauaktivitäten entsteht. Dies hat den Bedarf an Grundwasser als zuverlässigste Wasserversorgung erhöht. Unser Verständnis der Mechanismen der Grundwasserneubildung und der chemischen Entwicklung im Einzugsgebiet ist bislang unzureichend, was eine wirksame Bewirtschaftung erschwert. Daher ist das Hauptziel dieser Arbeit Einblicke in die Prozesse zu bekommen, welche die hydrogeochemische Entwicklung der Grundwasserqualität im Pra-Becken bestimmen. Die kombinierte Verwendung von Daten stabiler Isotope, der Hydrochemie und von Wasserständen bildet die Grundlage für die Konzeption der chemischen Entwicklung des Grundwassers im Pra-Becken. Dafür wurden die Herkunft und die Verdunstungsraten des in die ungesättigte Zone infiltrierenden Wassers bewertet. Darüber hinaus wurden die Chlorid-Massenbilanz und die Wasserspiegelschwankungen betrachtet, um die Grundwasserneubildung für das Einzugsgebiet zu quantifizieren. Indizes wie der Wasserqualitätsindex (WQI), das Natriumadsorptionsverhältnis (SAR), das Wilcox-Diagramm und der Salzgehalt (USSL) wurden in dieser Studie verwendet, um die Qualität der Ressource für die Verwendung als Trinkwasser und zu Bewässerungszwecken zu bestimmen. Aufgrund der Heterogenität der hydrochemischen Daten wurden die statistischen Verfahren der hierarchischen Cluster- und Faktorenanalyse angewandt, um die Daten entsprechend ihrer räumlichen Korrelation zu unterteilen. Ein konzeptionelles hydrogeochemisches Modell wurde entwickelt und anschließend durch Anwendung kombinatorischer inverser und reaktionspfadbasierter geochemischer Modelle validiert, um plausible mineralische Assemblagen zu bestimmen, welche die chemische Zusammensetzung des Grundwassers kontrollieren. Die Wechselwirkungen zwischen Wasser und Gestein bestimmen die Grundwasserqualität im Pra-Becken. Die Ergebnisse unterstreichen, dass das Grundwasser eine gute Qualität aufweist und als Trinkwasser und für Bewässerungszwecke genutzt werden kann. Es wurde nachgewiesen, dass ein großes Grundwasserpotenzial vorhanden ist, um den derzeitigen und künftigen Wasserbedarf des gesamten Pra-Beckens zu decken. Als Hauptneubildungsgebiet wurde die nördliche Zone im Gebiet identifiziert, während die südliche Zone das Abflussgebiet ist. Der vorherrschende Einfluss der Verwitterung von Silikatmineralen spielt bei der chemischen Entwicklung des Grundwassers eine zentrale Rolle. Die hier vorgestellte Arbeit gibt grundlegende Einblicke in die Hydrochemie des Pra-Beckens und liefert für das Wassermanagement wichtige Daten für eine fundierte Entscheidungsfindung bei der Planung und Zuweisung von Wasserressourcen für verschiedene Zwecke. Ein neuartiger Ansatz zur inversen Modellierungwurde in dieser Studie eingesetzt, um unterschiedliche Mineralzusammensetzungen zu ermitteln, welche die chemische Entwicklung des Grundwassers im Pra-Becken bestimmen. Diese Modellierungstechnik hat das Potenzial, die Zusammensetzung eines Grundwassers auf der Skala eines Beckens mit großer hydrochemischer Heterogenität zu simulieren, wobei die durchschnittliche Wasserzusammensetzung zur Darstellung der etablierten räumlichen Gruppierungen der Wasserchemie verwendet wird. KW - PHREEQC KW - PHREEQC KW - combinatorial inverse modelling KW - kombinatorische inverse Modellierung KW - reaction path modelling KW - Reaktionspfadmodellierung KW - groundwater evolution KW - Grundwasserentwicklung KW - water quality KW - Wasserqualität KW - rock-water interaction KW - Gestein-Wasser-Wechselwirkung KW - silicate weathering KW - Silikatverwitterung Y1 - 2024 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-628062 ER - TY - THES A1 - Erbello Doelesso, Asfaw T1 - Cenozoic magma-assisted continental rifting and crustal block rotations in an extensional overlap zone between two rift segments, Southwest Ethiopia, East Africa N2 - Continental rifts are key geodynamic regions where the complex interplay of magmatism and faulting activity can be studied to understand the driving forces of extension and the formation of new divergent plate boundaries. Well-preserved rift morphology can provide a wealth of information on the growth, interaction, and linkage of normal-fault systems through time. If rift basins are preserved over longer geologic time periods, sedimentary archives generated during extensional processes may mirror tectonic and climatic influences on erosional and sedimentary processes that have varied over time. Rift basins are furthermore strategic areas for hydrocarbon and geothermal energy exploration, and they play a central role in species dispersal and evolution as well as providing or inhibiting hydrologic connectivity along basins at emerging plate boundaries. The Cenozoic East African rift system (EARS) is one of the most important continental extension zones, reflecting a range of evolutionary stages from an early rift stage with isolated basins in Malawi to an advanced stage of continental extension in southern Afar. Consequently, the EARS is an ideal natural laboratory that lends itself to the study of different stages in the breakup of a continent. The volcanically and seismically active eastern branch of the EARS is characterized by multiple, laterally offset tectonic and magmatic segments where adjacent extensional basins facilitate crustal extension either across a broad deformation zone or via major transfer faulting. The Broadly Rifted Zone (BRZ) in southern Ethiopia is an integral part of the eastern branch of the EARS; in this region, rift segments of the southern Ethiopian Rift (sMER) and northern Kenyan Rift (nKR) propagate in opposite directions in a region with one of the earliest manifestations of volcanism and extensional tectonism in East Africa. The basin margins of the Chew-Bahir Basin and the Gofa Province, characterized by a semi-arid climate and largely uniform lithology, provide ideal conditions for studying the tectonic and geomorphologic features of this complex kinematic transfer zone, but more importantly, this area is suitable for characterizing and quantifying the overlap between the propagating structures of the sMER and nKR and the resulting deformation patterns of the BRZ transfer zones. In this study, I have combined data from thermochronology, thermal modeling, morphometry, paleomagnetic analysis, geochronology, and geomorphological field observations with information from published studies to reconstruct the spatiotemporal relationship between volcanism and fault activity in the BRZ and quantify the deformation patterns of the overlapping rift segments. I present the following results: (1) new thermochronological data from the en-échelon basin margins and footwall blocks of the rift flanks and morphometric results verified in the field to link different phases of magmatism and faulting during extension and infer geomorphological landscape features related to the current tectonic interaction between the nKR and the sMER; (2) temporally constrained paleomagnetic data from the BRZ overlap zone between the Ethiopian and Kenyan rifts to quantitatively determine block rotation between the two segments. Combining the collected data, time-temperature histories of thermal modeling results from representative samples show well-defined deformation phases between 25–20 Ma, 15–9Ma, and ~5 Ma to the present. Each deformation phase is characterized by the onset of rapid cooling (>2°C/Ma) of the crust associated with uplift or exhumation of the rift shoulder. After an initial, spatially very diffuse phase of extension, the rift has gradually evolved into a system of connected structures formed in an increasingly focused rift zone during the last 5 Ma. Regarding the morphometric analysis of the rift structures, it can be shown that normalized slope indices of the river courses, spatial arrangement of knickpoints in the river longitudinal profiles of the footwall blocks, local relief values, and the average maximum values of the slope of the river profiles indicate a gradual increase in the extension rate from north (Sawula basin: mature) to south (Chew Bahir: young). The complexity of the structural evolution of the BRZ overlap zone between nKR and sMER is further emphasized by the documentation of crustal blocks around a vertical axis. A comparison of the mean directions obtained for the Eo-Oligocene (Ds=352.6°, Is=-17.0°, N=18, α95=5.5°) and Miocene (Ds=2.9°, Is=0.9°, N=9, α95=12.4°) volcanics relative to the pole for stable South Africa and with respect to the corresponding ages of the analyzed units record a significant counterclockwise rotation of ~11.1°± 6.4° and insignificant CCW rotation of ~3.2° ± 11.5°, respectively. N2 - Wie Kontinentalverschiebungen entstehen und sich im Laufe der Zeit entwickeln, ist nach wie vor eine offene Frage in den Geowissenschaften. Hier untersuche ich das Rifting entlang der EARS, insbesondere in der BRZ zwischen dem nKR und dem sMER, einem aktiven Rift auf der Erde. Einige Studien legen nahe, dass ein Mantelquellenauftrieb von geschmolzenem Material und die damit verbundenen Mantelkonvektionsströme eine Krustenschwächung und -dehnung verursachen, die zur Bildung des vulkanisch und seismisch aktiven Rifts in Ostafrika führt. Im Gegensatz dazu argumentieren andere, dass eine Fernfeldbelastung im Zusammenhang mit einer benachbarten Plattenkonfiguration zwischen der eurasischen und der arabischen Platte das Rifting in EARS ausgelöst haben könnte. In diesem Projekt habe ich mit Hilfe von Tieftemperatur-Thermochronologie, Paläomagnetismus und Flusseinzugsgebietsmorphometrie die Beziehungen zwischen Magmatismus, Verwerfungen, einer Blockdeformation zwischen den Rifts und der Entwicklung der Oberflächentopographie während des Riftings untersucht. Basierend auf den Ergebnissen der Modellierung der Zeit-Temperatur-Geschichte stelle ich fest, dass die frühe Deformation in Verbindung mit Verwerfungen bei ~22 Ma begann, nach dem Ende des massiven Vulkanismus, der zwischen 45-27 Ma stattfand. Die räumlich verteilte anfängliche Deformation wanderte auf eine schmale Zone zu und breitete sich entlang der Gofa-Provinz nach Norden aus, bevor sich die aktuelle Verformung in der südlichen Gofa-Provinz und im Chew Bahir-Becken während einer späten Phase des Riftings bei ~15-9 Ma bzw. in der Gegenwart lokalisierte. Die räumliche Verschiebung der Deformation entspricht zeitlich einem erneuten Magmatismus. Kombinierte paläomagnetische und geochronologische Ergebnisse zeigen eine gegen den Uhrzeigersinn gerichtete vertikale Achsenblockrotation zwischen der sich ausbreitenden nKR und sMER, die im Laufe der Zeit abnahm. Die Abnahme der vertikalen Blockrotation gegen den Uhrzeigersinn in Verbindung mit der räumlichen und zeitlichen Veränderung der Deformation könnte auf eine charakteristische Rift-Evolutionsphase hinweisen, die von einer verteilten Anfangsdeformation und dem seitlichen Wachstum normaler Verwerfungen bis hin zur Lokalisierung der Dehnungen in der Spätphase des Rifting reicht, bevor sich ein durchgehendes größeres Riftbecken entwickelt. Die Ergebnisse meiner Studie, die zeigen, dass die anfängliche Deformation nach dem Ende des massiven Vulkanismus auftrat, unterstützen die Vorstellung einer Schwächung der Kruste durch massiven Vulkanismus vor dem anfänglichen Rifting als Vorläufer einer weit verbreiteten Dehnungsverwerfung und bieten somit weitere Einblicke in magmengestützte Deformationsprozesse im ostafrikanischen Riftsystem. KW - magma assisted continental rifting KW - rift segments interaction KW - counterclockwise block rotation between overlapping rift segments KW - magmagestütztes kontinentales Rifting KW - Interaktion zwischen sich ausbreitenden Riftsegmenten KW - gegen den Uhrzeigersinn gerichtete Rotation von Krustenblöcken zwischen zwei überlappenden Riftsegmenten Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-610968 ER - TY - THES A1 - Freisleben, Roland T1 - Deciphering the mechanisms of permanent forearc deformation based on marine terraces T1 - Untersuchung der Mechanismen dauerhafter Forearc-Deformation anhand von marinen Terrassen N2 - The Andes reflect Cenozoic deformation and uplift along the South American margin in the context of regional shortening associated with the interaction between the subducting Nazca plate and the overriding continental South American plate. Simultaneously, multiple levels of uplifted marine terraces constitute laterally continuous geomorphic features related to the accumulation of permanent forearc deformation in the coastal realm. However, the mechanisms responsible for permanent coastal uplift and the persistency of current/decadal deformation patterns over millennial timescales are still not fully understood. This dissertation presents a continental-scale database of last interglacial terrace elevations and uplift rates along the South American coast that provides the basis for an analysis of a variety of mechanisms that are possibly responsible for the accumulation of permanent coastal uplift. Regional-scale mapping and analysis of multiple, late Pleistocene terrace levels in central Chile furthermore provide valuable insights regarding the persistency of current seismic asperities, the role of upper-plate faulting, and the impact of bathymetric ridges on permanent forearc deformation. The database of last interglacial terrace elevations reveals an almost continuous signal of background-uplift rates along the South American coast at ~0.22 mm/yr that is modified by various short- to long-wavelength changes. Spatial correlations with crustal faults and subducted bathymetric ridges suggest long-term deformation to be affected by these features, while the latitudinal variability of climate forcing factors has a profound impact on the generation and preservation of marine terraces. Systematic wavelength analyses and comparisons of the terrace-uplift rate signal with different tectonic parameters reveal short-wavelength deformation to result from crustal faulting, while intermediate- to long-wavelength deformation might indicate various extents of long-term seismotectonic segments on the megathrust, which are at least partially controlled by the subduction of bathymetric anomalies. The observed signal of background-uplift rate is likely accumulated by moderate earthquakes near the Moho, suggesting multiple, spatiotemporally distinct phases of uplift that manifest as a continuous uplift signal over millennial timescales. Various levels of late Pleistocene marine terraces in the 2015 M8.3 Illapel-earthquake area reveal a range of uplift rates between 0.1 and 0.6 mm/yr and indicate decreasing uplift rates since ~400 ka. These glacial-cycle uplift rates do not correlate with current or decadal estimates of coastal deformation suggesting seismic asperities not to be persistent features on the megathrust that control the accumulation of permanent forearc deformation over long timescales of 105 years. Trench-parallel, crustal normal faults modulate the characteristics of permanent forearc-deformation; upper-plate extension likely represents a second-order phenomenon resulting from subduction erosion and subsequent underplating that lead to regional tectonic uplift and local gravitational collapse of the forearc. In addition, variable activity with respect to the subduction of the Juan Fernández Ridge can be detected in the upper plate over the course of multiple interglacial periods, emphasizing the role of bathymetric anomalies in causing local increases in terrace-uplift rate. This thesis therefore provides new insights into the current understanding of subduction-zone processes and the dynamics of coastal forearc deformation, whose different interacting forcing factors impact the topographic and geomorphic evolution of the western South American coast. N2 - Die känozoischen Anden resultieren aus der regionalen Verkürzung und Hebung der kontinentalen Kruste entlang des südamerikanischen Kontinentalrandes infolge der Interaktion zwischen der subduzierenden Nazca-Platte und der südamerikanischen Platte. Zahlreiche, durch tektonische Prozesse angehobene marine Abrasionsplattformen entlang der Küsten am westlichen Kontinentalrand Südamerikas bilden lateral kontinuierliche Stufen, welche die Akkumulation dauerhafter Deformation im küstennahen Forearc-Bereich widerspiegeln. Die Mechanismen, welche für die dauerhafte Hebung der Küste und die Beständigkeit derzeitiger/dekadischer Deformationsmuster auch über Jahrtausende hinweg verantwortlich sind, sind jedoch noch nicht vollständig geklärt. Vor diesem Hintergrund wird in dieser Dissertation zunächst eine kontinentale Datenbank der gehobenen marinen Terrassen des letzten Interglazials (125,000 Jahre), ihrer Höhenverteilung sowie der daraus abgeleiteten Hebungsraten für die Westküste von Südamerika vorgestellt. Diese Datenbank bildet die Grundlage für die Analyse einer Vielzahl von Mechanismen, welche möglicherweise für die Akkumulation dauerhafter Küstenhebung und Deformation verantwortlich sind. In einem weiteren Schritt wurden zusätzlich spätpleistozäne Terrassenniveaus in Zentralchile regionalmaßstäblich kartiert und hinsichtlich ihrer räumlich-zeitlichen Entwicklung analysiert. Diese Informationen liefern wertvolle Erkenntnisse über die Beständigkeit von seismischen Asperitäten, den Einfluss bathymetrischer Rücken auf die dauerhafte Deformation der Oberplatte sowie die Rolle von krustalen Störungen hinsichtlich der lokalen topographischen Differenzierung seismotektonischer Segmente. Die Datenbank der im letzten Interglazial entstandenen marinen Terrassen und der abgeleiteten Hebungsraten zeigt ein nahezu kontinuierliches Signal von Hintergrund-Hebungsraten entlang der südamerikanischen Küste von ~0,22 mm/Jahr, das durch verschiedene kurz- bis langwellige Modifikationen gekennzeichnet ist. Räumliche Übereinstimmungen mit krustalen Störungen und subduzierten bathymetrischen Rücken deuten darauf hin, dass die langfristige Deformation von diesen Parametern beeinflusst wird, während die breitengradabhängige Variabilität von Klimafaktoren tiefgreifende Auswirkungen auf die Entstehung und Erhaltung von marinen Terrassen hat. Systematische Wellenlängenanalysen und Vergleiche des Signals der Terrassenhebungsrate mit verschiedenen tektonischen Parametern zeigen, dass Deformation im kurzwelligen Bereich auf krustale Störungen zurückzuführen ist. Die Deformation im mittel- bis langwelligen Bereich könnte hingegen auf verschiedene Ausdehnungen langfristig aktiver seismotektonischer Segmente der Subduktionszone hindeuten, die zumindest teilweise durch die Subduktion bathymetrischer Anomalien kontrolliert werden. Das beobachtete Signal der Hintergrund-Hebungsrate wird wahrscheinlich durch Erdbeben mittlerer Magnituden in der Nähe der Moho generiert, was auf mehrere, räumlich und zeitlich getrennte Hebungsphasen hindeutet, die sich über Jahrtausende hinweg als kontinuierliches Hebungssignal manifestieren. Gehobene spätpleistozäne marine Terrassen im Illapel-Erdbebengebiet von 2015 weisen eine Hebungsrate zwischen 0,1 und 0,6 mm/Jahr auf und deuten insgesamt auf abnehmende tektonische Hebungsraten seit ~400,000 Jahren hin. Diese langfristigen Hebungsraten auf glazialen Zeitskalen korrelieren nicht mit den Hebungsraten auf kurzen Zeitskalen, was darauf schließen lässt, dass seismische Asperitäten im Bereich der Subduktionszone keine dauerhaften Merkmale sind, welche die Deformation und Hebung der Forearc-Regionen über lange Zeiträume von mehreren 105 Jahren steuern. Wie auch an anderen aktiven konvergenten Plattenrändern existieren in Zentralchile krustale Abschiebungen parallel zum Plattenrand, die das Erscheinungsbild dauerhafter Forearc-Deformation stark beeinflussen. Die Extensionsprozesse der Oberplatte sind dabei wahrscheinlich ein Phänomen zweiter Ordnung, welches durch Subduktionserosion und anschließende Akkretion hervorgerufen wird und sich in regionaler tektonischer Hebung sowie einem gravitativen Kollaps der Forearc-Region äußert. Darüber hinaus lassen sich über verschiedene Glazial- und Interglazialzyklen hinweg Änderungen in der Bedeutung des Juan-Fernández-Rückens hinsichtlich der Deformations-prozesse in der Oberplatte nachweisen; dies unterstreicht die transiente Rolle bathymetrischer Anomalien bei der Beeinflussung lokaler Änderungen der Terrassenhebungsrate. Die vorgelegte Arbeit liefert daher neue Einblicke für ein besseres Verständnis der Deformationsprozesse in Subduktionszonen und den küstennahen Forearc-Regionen sowie ihrer unterschiedlich interagierenden Antriebsfaktoren, welche die topographische und geomorphologische Entwicklung der westlichen südamerikanischen Küste beeinflussen. KW - South America KW - marine terrace KW - tectonic geomorphology KW - subduction KW - Südamerika KW - marine Terrassen KW - Subduktion KW - tektonische Geomorphologie Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-610359 ER - TY - THES A1 - Metz, Malte T1 - Finite fault earthquake source inversions T1 - Ausgedehnte Erdbebenquellinversion BT - implementation and testing of a novel physics-based rupture model BT - Implementierung und Validierung eines neuen selbst-ähnlichen Bruchmodels N2 - Earthquake modeling is the key to a profound understanding of a rupture. Its kinematics or dynamics are derived from advanced rupture models that allow, for example, to reconstruct the direction and velocity of the rupture front or the evolving slip distribution behind the rupture front. Such models are often parameterized by a lattice of interacting sub-faults with many degrees of freedom, where, for example, the time history of the slip and rake on each sub-fault are inverted. To avoid overfitting or other numerical instabilities during a finite-fault estimation, most models are stabilized by geometric rather than physical constraints such as smoothing. As a basis for the inversion approach of this study, we build on a new pseudo-dynamic rupture model (PDR) with only a few free parameters and a simple geometry as a physics-based solution of an earthquake rupture. The PDR derives the instantaneous slip from a given stress drop on the fault plane, with boundary conditions on the developing crack surface guaranteed at all times via a boundary element approach. As a side product, the source time function on each point on the rupture plane is not constraint and develops by itself without additional parametrization. The code was made publicly available as part of the Pyrocko and Grond Python packages. The approach was compared with conventional modeling for different earthquakes. For example, for the Mw 7.1 2016 Kumamoto, Japan, earthquake, the effects of geometric changes in the rupture surface on the slip and slip rate distributions could be reproduced by simply projecting stress vectors. For the Mw 7.5 2018 Palu, Indonesia, strike-slip earthquake, we also modelled rupture propagation using the 2D Eikonal equation and assuming a linear relationship between rupture and shear wave velocity. This allowed us to give a deeper and faster propagating rupture front and the resulting upward refraction as a new possible explanation for the apparent supershear observed at the Earth's surface. The thesis investigates three aspects of earthquake inversion using PDR: (1) to test whether implementing a simplified rupture model with few parameters into a probabilistic Bayesian scheme without constraining geometric parameters is feasible, and whether this leads to fast and robust results that can be used for subsequent fast information systems (e.g., ground motion predictions). (2) To investigate whether combining broadband and strong-motion seismic records together with near-field ground deformation data improves the reliability of estimated rupture models in a Bayesian inversion. (3) To investigate whether a complex rupture can be represented by the inversion of multiple PDR sources and for what type of earthquakes this is recommended. I developed the PDR inversion approach and applied the joint data inversions to two seismic sequences in different tectonic settings. Using multiple frequency bands and a multiple source inversion approach, I captured the multi-modal behaviour of the Mw 8.2 2021 South Sandwich subduction earthquake with a large, curved and slow rupturing shallow earthquake bounded by two faster and deeper smaller events. I could cross-validate the results with other methods, i.e., P-wave energy back-projection, a clustering analysis of aftershocks and a simple tsunami forward model. The joint analysis of ground deformation and seismic data within a multiple source inversion also shed light on an earthquake triplet, which occurred in July 2022 in SE Iran. From the inversion and aftershock relocalization, I found indications for a vertical separation between the shallower mainshocks within the sedimentary cover and deeper aftershocks at the sediment-basement interface. The vertical offset could be caused by the ductile response of the evident salt layer to stress perturbations from the mainshocks. The applications highlight the versatility of the simple PDR in probabilistic seismic source inversion capturing features of rather different, complex earthquakes. Limitations, as the evident focus on the major slip patches of the rupture are discussed as well as differences to other finite fault modeling methods. N2 - Erdbebenmodelle sind der Schlüssel zu einem detaillierten Verständnis der zugrunde liegenden Bruchprozesse. Die kinematischen oder dynamischen Brucheigenschaften werden mit Hilfe von ausgedehnten Bruchmodellen bestimmt. Dadurch können Details, wie z.B. die Bruchrichtung und -geschwindigkeit oder die Verschiebungsverteilung, aufgelöst werden. Häufig sind ausgedehnte Bruchmodelle durch sehr viele freie Parameter definiert, etwa individuelle Verschiebungen und Verschiebungsrichtungen auf den diskretisierten Bruchflächenelementen. Die große Anzahl an Parametern sorgt dafür, dass Inversionsprobleme hochgradig unterbestimmt sind. Um daraus resultierende numerische Instabilitäten zu verhinden, werden diese Modelle häufig mit zusätzlichen eher geometrischen als physikalischen Annahmen stabilisiert, z.B. im Bezug auf die Rauigkeit der Verschiebung auf der Bruchfläche. Die Basis für die Inversionsmethode in dieser Dissertaton bildet das von uns entwickelete pseudo-dynamische Bruchmodel (PDR). Die PDR basiert auf wenigen freien Parametern und einer simplen, planaren Geometrie und ergibt eine physik-gestützte Lösung für Erdbebenbrüche. Die PDR bestimmt die instantane Verschiebung basierend auf gegebenen Spannungsänderungen auf der Bruchfläche. Die Randbedingung der Spannungsänderung wird dabei zu jedem Zeitpunkt der Bruchentwicklung über eine Randelementmethode eingehalten. Als Nebenprodukt dessen kann die Herdzeitfunktion an jedem Punkt der Bruchfläche als Ergebnis des Models bestimmt werden, und muss daher nicht vorher definiert werden. Der PDR-Modellierungsansatz wurde mit anderen Modellen anhand verschiedener Erdbeben verglichen. Am Beispiel des Mw 7,1 2016 Kumamoto, Japan, Bebens konnte der Effekt einer gekrümmten Bruchfläche auf die daraus resultierenden Verschiebungsverteilung und Verschiebungsraten durch eine Projezierung der Spannungsvektoren reproduziert werden. Für das Mw 7,5 2018 Palu, Indonesien, Beben haben wir die Bruchausbreitung auf Grundlage der 2D-Eikonalgleichung und basierend auf einem angenommenen linearen Zusammenhang zwischen Bruch- und Scherwellengeschwindigkeit modelliert. Dadurch konnten wir die beobachtete Supershear-Bruchausbreitung als Ergebnis einer möglichen tiefen und daher schnelleren Bruchfront mit einer Abstrahlung an die Erdoberfläche erklären. Der PDR-Vorwärtsmodellierungs-Code wurde in den Open-Source Python Paketen Pyrocko und Grond veröffentlicht. Meine Dissertation beleuchtet drei Aspekte der Erdbebeninversion unter Zuhilfenahme der PDR: (1) Ist eine Implementation eines simplen Bruchmodels mit wenigen Parametern in ein probabalistisches Bayesisches Inversionsprogramm möglich? Kann dies schnelle und robuste Ergebnisse für weitere Folgeanwendungen, wie Bodenbeschleunigungsvorhersagen, liefern? (2) Wie hilft die Kombination aus seismischen Breitband- und Accelerometerdaten mit Nahfelddeformationsdaten, Inversionsergebnisse mit der PDR zu verbessern? (3) Können komplexe Brüche über einen multiplen PDR-Quellinversionsansatz aufgelöst werden und wenn ja, wann ist dies möglich? Ich habe den PDR-Inversionsansatz entwickelt und auf zwei Erdbeben-Sequenzen in verschiedenen tektonischen Umgebungen angewandt. Mit Hilfe von verschiedenen Datensätzen in mehreren Frequenzbändern innerhalb von einfachen und multiplen Bruchflächeninversionen konnte ich das multi-modale Mw 8,2 2021 South Sandwich Erdbeben characterisieren. Dieses bestand aus einem langen, flachen, langsam brechenden Beben entlang der gekrümmten Subduktionszone, welches durch zwei kleinere, tiefere Brüche mit schnelleren Bruchgeschwindigkeiten begrenzt wurde. Die Validierung mit Ergebnissen aus einer P-Wellen Back-Projection, der Clusteranalyse von Nachbeben und einer Tsunami-Modelierung zeigten eine hohe Konsistenz mit den PDR-Resultaten. Die Kombination von seismischen Daten und Oberflächendeformationen in einer multiplen PDR-Inversion habe ich auch zur Analyse eines Beben-Triplets vom Juni 2022 im Südosten des Irans genutzt. Die Inversionen konnten im Zusammenspiel mit relokalisierten Nachbeben einen neuen Fall von vertikaler Haupt-/Nachbebenseparation auflösen. Während die großen Hauptbeben im flachen Sediment stattfanden, sind die Nachbeben hauptsächlich entlang der tieferen Grenzfläche zwischen Sediment und kristallinem Grundgebirge aufgetreten. Eine Erklärung dafür ist das duktile Fließen einer vorhandenen Salzschicht auf der Grenzfläche, ausgelöst durch Spannungsänderungen im Zuge der Hauptbeben. Die Anwendungen konnten die Vielseitigkeit der PDR als simples Quellmodel innerhalb von seismischen Quellinversionen zeigen. Limitierungen der Inversion, wie der augenscheinliche Fokus auf den Hauptverschiebungsbereich eines Bebens, werden in dieser Arbeit genauso diskutiert wie die Einordnung der PDR im Vergleich zu anderen ausgedehnten Quellmodellen. KW - seismology KW - inversion KW - source model KW - Seismologie KW - Inversion KW - Bruchmodel Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-619745 ER - TY - THES A1 - Mantiloni, Lorenzo T1 - Modeling stress and dike pathways in calerdas T1 - Modellierung von Spannung und Magmagängen in Calderen BT - towards a physics-based forecast of eruptive vent locations BT - Entwicklung physikbasierter Vorhersagen von Eruptionsschlot-Lokationen N2 - Volcanic hazard assessment relies on physics-based models of hazards, such as lava flows and pyroclastic density currents, whose outcomes are very sensitive to the location where future eruptions will occur. On the contrary, forecast of vent opening locations in volcanic areas typically relies on purely data-driven approaches, where the spatial density of past eruptive vents informs the probability maps of future vent opening. Such techniques may be suboptimal in volcanic systems with missing or scarce data, and where the controls on magma pathways may change over time. An alternative approach was recently proposed, relying on a model of stress-driven pathways of magmatic dikes. In that approach, the crustal stress was optimized so that dike trajectories linked consistently the location of the magma chamber to that of past vents. The retrieved information on the stress state was then used to forecast future dike trajectories. The validation of such an approach requires extensive application to nature. Before doing so, however, several important limitations need to be removed, most importantly the two-dimensional (2D) character of the models and theoretical concepts. In this thesis, I develop methods and tools so that a physics-based strategy of stress inversion and eruptive vent forecast in volcanoes can be applied to three dimensional (3D) problems. In the first part, I test the stress inversion and vent forecast strategy on analog models, still within a 2D framework, but improving on the efficiency of the stress optimization. In the second part, I discuss how to correctly account for gravitational loading/unloading due to complex 3D topography with a Boundary-Element numerical model. Then, I develop a new, simplified but fast model of dike pathways in 3D, designed for running large numbers of simulations at minimal computational cost, and able to backtrack dike trajectories from vents on the surface. Finally, I combine the stress and dike models to simulate dike pathways in synthetic calderas. In the third part, I describe a framework of stress inversion and vent forecast strategy in 3D for calderas. The stress inversion relies on, first, describing the magma storage below a caldera in terms of a probability density function. Next, dike trajectories are backtracked from the known locations of past vents down through the crust, and the optimization algorithm seeks for the stress models which lead trajectories through the regions of highest probability. I apply the new strategy to the synthetic scenarios presented in the second part, and I exploit the results from the stress inversions to produce probability maps of future vent locations for some of those scenarios. In the fourth part, I present the inversion of different deformation source models applied to the ongoing ground deformation observed across the Rhenish Massif in Central Europe. The region includes the Eifel Volcanic Fields in Germany, a potential application case for the vent forecast strategy. The results show how the observed deformation may be due to melt accumulation in sub-horizontal structures in the lower crust or upper mantle. The thesis concludes with a discussion of the stress inversion and vent forecast strategy, its limitations and applicability to real volcanoes. Potential developments of the modeling tools and concepts presented here are also discussed, as well as possible applications to other geophysical problems. N2 - Die Analyse vulkanischer Gefahren basiert auf physikalischen Gefahrenmodellen wie Lavaströmen und pyroklastischen Dichteströmen, deren Ergebnisse von den Lokationen der Eruptionen abhängt. Die Vorhersage der Lokationen von Spaltenöffungen (eruptive vent opening) in Vulkangebieten basiert typischerweise auf rein datengesteuerten Ansätzen, bei denen die räumliche Verteilung vergangener Eruptionsspalten die Wahrscheinlichkeitskarten zukünftiger Spaltöffnungen beeinflusst. Solche Techniken sind in Vulkansystemen mit fehlenden oder knappen Daten suboptimal, ebenso wie dort, wo sich physikalische Bedingungen der Magmapfade im Laufe der Zeit verändern können. Kürzlich wurde ein alternativer Ansatz vorgeschlagen, der auf einem Modell spannungsbedingter Magmagänge basiert. Bei diesem Ansatz wurde die Krustenspannung so optimiert, dass die Verläufe der Magmagänge- sowohl mit der Lage der Magmakammer als auch mit der Lage früherer Spaltenöffnungen konsistent sind Die erlangten Informationen über den Spannungszustand werden dann zur Vorhersage zukünftiger Magmagänge verwendet. Die Validierung eines solchen Ansatzes erfordert eine umfassende Anwendung auf natürliche Vulkansysteme. Zuvor müssen jedoch einige wichtige Einschränkungen beseitigt werden, insbesondere der zweidimensionale (2D) Charakter der Modelle und theoretischen Konzepte. In dieser Arbeit entwickle ich Methoden und Software, um physikbasierte Strategien der Spannungsinversion und der Vorhersage von Lokationen von Magmagängen in Vulkanen auf dreidimensionale (3D) Probleme anwenden zu können. Im ersten Teil teste ich die Spannungsinversion und die Vorhersage von eruptiven Magmagängen an analogen Modellen, immer noch innerhalb eines 2D-Rahmens, aber mit verbesserter Effizienz der Spannungsoptimierung. Im zweiten Teil diskutiere ich, wie die gravitative Belastung/Entlastung aufgrund einer komplexen 3D-Topographie mit einem numerischen Grenzelementmodell berücksichtigt werden kann. Anschließend entwickle ich ein neues, vereinfachtes, aber effizientes Modell der Magmagang-Ausbreitung in 3D Dieses ist für die Durchführung einer großen Anzahl von Simulationen mit minimalem Rechenaufwand konzipiert und dabei in der Lage, Magmagang-Ausbreitungen an der Oberfläche zurückzuverfolgen. Schließlich kombiniere ich die Spannungs- und Magma-Ausbreitungsmodelle, um Magmagänge in synthetischen Calderas zu simulieren. Im dritten Teil erarbeite ich ein Rahmenmodell für Spannungsinversion und dreidimensionale Magmagang-Vorhersagen für Calderas. Die Spannungsinversion beruht auf der Annahme einer Magmaspeicherung unterhalb einer Caldera, deren Geometrie durch eine Wahrscheinlichkeitsdichtefunktion beschrieben wird. Vergangene Eruptionsstandorte werden durch die Erdkruste zurück propagiert, dabei werden die Spannungsmodelle so optimiert, dass die Endpunkte der Optimierung (Starpunkte der Magmagänge) mit der angenommenen Magmakammer übereinstimmen. Ich verwende die synthetischen Szenarien aus dem zweiten Teil, um die Spannungsinversion in kontrollieren Szenarien zu testen. Im vierten Teil analysieren wir die Bodenverformung im Rheinischen Massiv, Mitteleuropa. Die Region umfasst die Eifel-Vulkanfelder in Deutschland, die einen potenziellen Anwendungsfall für die Strategie zur Vorhersage von Spalteneruptionen darsetellen. Dabei nutze ich in den Inversionen unterschiedliche Deformations-Quellmodelle, um jährliche Magmavolumenänderungen aufzulösen, unter der Annahme, das Magma in die Lithosphäre injiziert wird. Die Ergebnisse zeigen, dass die beobachtete Verformung auf die Ansammlung von Schmelze in subhorizontalen Strukturen in der unteren Kruste oder im oberen Mantel zurückzuführen sein könnte. Abschließend diskutiere ich die Limitationen der erarbeiteten Konzepte und die Möglichkeiten für Anwendungen sowohl in verschiedenen Hochrisko-Vulkanen, als auch in anderen Kontexten, wie besispielsweise in tektonischen und geothermischen Modellierungen. KW - dike pathways KW - stress modeling KW - volcanic hazard assessment KW - Calderas KW - Magmagänge KW - Stressmodellierung KW - Risikobewertung von Vulkanausbrüchen KW - Calderas Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-612621 ER - TY - THES A1 - Illien, Luc T1 - Time-dependent properties of the shallow subsurface BT - groundwater and earthquake damage dynamics from seismic interferometry N2 - The shallow Earth’s layers are at the interplay of many physical processes: some being driven by atmospheric forcing (precipitation, temperature...) whereas others take their origins at depth, for instance ground shaking due to seismic activity. These forcings cause the subsurface to continuously change its mechanical properties, therefore modulating the strength of the surface geomaterials and hydrological fluxes. Because our societies settle and rely on the layers hosting these time-dependent properties, constraining the hydro-mechanical dynamics of the shallow subsurface is crucial for our future geographical development. One way to investigate the ever-changing physical changes occurring under our feet is through the inference of seismic velocity changes from ambient noise, a technique called seismic interferometry. In this dissertation, I use this method to monitor the evolution of groundwater storage and damage induced by earthquakes. Two research lines are investigated that comprise the key controls of groundwater recharge in steep landscapes and the predictability and duration of the transient physical properties due to earthquake ground shaking. These two types of dynamics modulate each other and influence the velocity changes in ways that are challenging to disentangle. A part of my doctoral research also addresses this interaction. Seismic data from a range of field settings spanning several climatic conditions (wet to arid climate) in various seismic-prone areas are considered. I constrain the obtained seismic velocity time-series using simple physical models, independent dataset, geophysical tools and nonlinear analysis. Additionally, a methodological development is proposed to improve the time-resolution of passive seismic monitoring. N2 - Die oberflächennahen Erdschichten sind ein Zusammenspiel vieler physikalischer Prozesse: Einige werden durch atmosphärische Einflüsse (Niederschlag, Temperatur...) angetrieben, während andere ihren Ursprung im Untergrund haben, beispielsweise Bodenerschütterungen aufgrund seismischer Aktivität. Diese Kräfte führen zu sich ständig ändernden mechanische Eigenschaften des Untergrunds, und damit auch Änderungen der Festigkeit der Geomaterialien der Oberfläche, sowie der Wasserflüsse. Als Gesellschaft leben wir auf der Erdoberfläche, und sind auf diese zeitabhängigen Eigenschaften angewiesen. Für unsere zukünftige, geographische Entwicklung ist es daher essentiell, die hydromechanischen Dynamiken des flachen Untergrunds zu verstehen. Eine Möglichkeit, die sich stetig ändernden physikalischen Eigenschaften unter unseren Füßen zu untersuchen, ist die Ableitung seismischer Geschwindigkeitsänderungen aus dem Umgebungslärm, mit Hilfe einer Technik namens seismische Interferometrie. In dieser Dissertation verwende ich diese Methode, um die Entwicklung der Grundwasserspeicherung, sowie durch Erdbeben verursachte Schäden zu untersuchen. Es werden zwei Forschungslinien untersucht, diese umfassen die Hauptkontrollen der Grundwasserneubildung in steilen Landschaften und die Vorhersagbarkeit und Dauer der transienten physikalischen Eigenschaften aufgrund von Erdbebenerschütterungen. Beide Dynamiken beeinflussen sich gegenseitig, sowie die Geschwindigkeitsänderungen im Untergrund in einer Weise, die schwer zu entschlüsseln ist. Ein Teil dieser Doktorarbeit behandelt daher auch diese Wechselwirkung. Wir verwenden seismische Daten, die unter verschiedenen klimatischen Bedingungen (feuchtes bis trockenes Klima) in mehreren, seismisch gefährdeten Gebieten erhoben wurden. Die gewonnenen seismischen Geschwindigkeitszeitreihen werden mit Hilfe einfacher, physikalischer Modelle, unabhängiger Datensätze, geophysikalischer Instrumente und nichtlinearer Analysen untersucht. Darüber hinaus wird eine methodische Entwicklung vorgeschlagen, um die zeitliche Auflösung des passiven, seismischen Monitorings zu verbessern. KW - seismology KW - seismic noise KW - mountain hydrology KW - earthquake damage KW - geomechanics KW - geophysics KW - Erdbebenschäden KW - Geomechanik KW - Geophysik KW - Gebirgshydrologie KW - seismisches Rauschen KW - Seismologie Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-599367 ER - TY - THES A1 - Gómez Zapata, Juan Camilo T1 - Towards unifying approaches in exposure modelling for scenario-based multi-hazard risk assessments N2 - This cumulative thesis presents a stepwise investigation of the exposure modelling process for risk assessment due to natural hazards while highlighting its, to date, not much-discussed importance and associated uncertainties. Although “exposure” refers to a very broad concept of everything (and everyone) that is susceptible to damage, in this thesis it is narrowed down to the modelling of large-area residential building stocks. Classical building exposure models for risk applications have been constructed fully relying on unverified expert elicitation over data sources (e.g., outdated census datasets), and hence have been implicitly assumed to be static in time and in space. Moreover, their spatial representation has also typically been simplified by geographically aggregating the inferred composition onto coarse administrative units whose boundaries do not always capture the spatial variability of the hazard intensities required for accurate risk assessments. These two shortcomings and the related epistemic uncertainties embedded within exposure models are tackled in the first three chapters of the thesis. The exposure composition of large-area residential building stocks is studied on the scope of scenario-based earthquake loss models. Then, the proposal of optimal spatial aggregation areas of exposure models for various hazard-related vulnerabilities is presented, focusing on ground-shaking and tsunami risks. Subsequently, once the experience is gained in the study of the composition and spatial aggregation of exposure for various hazards, this thesis moves towards a multi-hazard context while addressing cumulative damage and losses due to consecutive hazard scenarios. This is achieved by proposing a novel method to account for the pre-existing damage descriptions on building portfolios as a key input to account for scenario-based multi-risk assessment. Finally, this thesis shows how the integration of the aforementioned elements can be used in risk communication practices. This is done through a modular architecture based on the exploration of quantitative risk scenarios that are contrasted with social risk perceptions of the directly exposed communities to natural hazards. In Chapter 1, a Bayesian approach is proposed to update the prior assumptions on such composition (i.e., proportions per building typology). This is achieved by integrating high-quality real observations and then capturing the intrinsic probabilistic nature of the exposure model. Such observations are accounted as real evidence from both: field inspections (Chapter 2) and freely available data sources to update existing (but outdated) exposure models (Chapter 3). In these two chapters, earthquake scenarios with parametrised ground motion fields were transversally used to investigate the role of such epistemic uncertainties related to the exposure composition through sensitivity analyses. Parametrised scenarios of seismic ground shaking were the hazard input utilised to study the physical vulnerability of building portfolios. The second issue that was investigated, which refers to the spatial aggregation of building exposure models, was investigated within two decoupled vulnerability contexts: due to seismic ground shaking through the integration of remote sensing techniques (Chapter 3); and within a multi-hazard context by integrating the occurrence of associated tsunamis (Chapter 4). Therein, a careful selection of the spatial aggregation entities while pursuing computational efficiency and accuracy in the risk estimates due to such independent hazard scenarios (i.e., earthquake and tsunami) are discussed. Therefore, in this thesis, the physical vulnerability of large-area building portfolios due to tsunamis is considered through two main frames: considering and disregarding the interaction at the vulnerability level, through consecutive and decoupled hazard scenarios respectively, which were then contrasted. Contrary to Chapter 4, where no cumulative damages are addressed, in Chapter 5, data and approaches, which were already generated in former sections, are integrated with a novel modular method to ultimately study the likely interactions at the vulnerability level on building portfolios. This is tested by evaluating cumulative damages and losses after earthquakes with increasing magnitude followed by their respective tsunamis. Such a novel method is grounded on the possibility of re-using existing fragility models within a probabilistic framework. The same approach is followed in Chapter 6 to forecast the likely cumulative damages to be experienced by a building stock located in a volcanic multi-hazard setting (ash-fall and lahars). In that section, special focus was made on the manner the forecasted loss metrics are communicated to locally exposed communities. Co-existing quantitative scientific approaches (i.e., comprehensive exposure models; explorative risk scenarios involving single and multiple hazards) and semi-qualitative social risk perception (i.e., level of understanding that the exposed communities have about their own risk) were jointly considered. Such an integration ultimately allowed this thesis to also contribute to enhancing preparedness, science divulgation at the local level as well as technology transfer initiatives. Finally, a synthesis of this thesis along with some perspectives for improvement and future work are presented. N2 - Diese kumulative Diplomarbeit stellt eine schrittweise Untersuchung des Expositionsmodellierungsprozesses für die Risikobewertung durch Naturgefahren dar und weist auf seine bisher wenig diskutierte Bedeutung und die damit verbundenen Unsicherheiten hin. Obwohl sich „Exposition“ auf einen sehr weiten Begriff von allem (und jedem) bezieht, der für Schäden anfällig ist, wird er in dieser Arbeit auf die Modellierung von großräumigen Wohngebäudebeständen eingeengt. Klassische Gebäudeexpositionsmodelle für Risikoanwendungen wurden vollständig auf der Grundlage unbestätigter Expertenerhebungen über Datenquellen (z. B. veraltete Volkszählungsdatensätze) erstellt und wurden daher implizit als zeitlich und räumlich statisch angenommen. Darüber hinaus wurde ihre räumliche Darstellung typischerweise auch vereinfacht, indem die abgeleitete Zusammensetzung geografisch auf grobe Verwaltungseinheiten aggregiert wurde, deren Grenzen nicht immer die räumliche Variabilität der Gefahrenintensitäten erfassen, die für genaue Risikobewertungen erforderlich sind. Diese beiden Mängel und die damit verbundenen epistemischen Unsicherheiten, die in Expositionsmodellen eingebettet sind, werden in den ersten drei Kapiteln der Dissertation verfolgt. Die Exposure-Zusammensetzung von großflächigen Wohngebäudebeständen wird im Rahmen szenariobasierter Erdbebenschadenmodelle untersucht. Anschließend wird der Vorschlag optimaler räumlicher Aggregationsbereiche von Expositionsmodellen für verschiedene gefahrenbezogene Anfälligkeiten präsentiert, wobei der Schwerpunkt auf Bodenerschütterungs- und Tsunami-Risiken liegt. Anschließend, sobald die Erfahrung in der Untersuchung der Zusammensetzung und räumlichen Aggregation der Exposition für verschiedene Gefahren gesammelt wurde, bewegt sich diese Arbeit in Richtung eines Kontextes mit mehreren Gefahren, während sie sich mit kumulativen Schäden und Verlusten aufgrund aufeinanderfolgender Gefahrenszenarien befasst. Dies wird erreicht, indem eine neuartige Methode vorgeschlagen wird, um die bereits bestehenden Schadensbeschreibungen an Gebäudeportfolios als Schlüsseleingabe für die Berücksichtigung einer szenariobasierten Multi-Risiko-Bewertung zu berücksichtigen. Abschließend zeigt diese Arbeit, wie die Integration der oben genannten Elemente in der Risikokommunikation genutzt werden kann. Dies erfolgt durch eine modulare Architektur, die auf der Untersuchung quantitativer Risikoszenarien basiert, die mit der sozialen Risikowahrnehmung der direkt von Naturgefahren betroffenen Gemeinschaften kontrastiert werden. In Kapitel 1 wird ein bayesianischer Ansatz vorgeschlagen, um die früheren Annahmen zu einer solchen Zusammensetzung (d. h. Anteile pro Gebäudetypologie) zu aktualisieren. Dies wird erreicht, indem hochwertige reale Beobachtungen integriert und dann die intrinsische Wahrscheinlichkeitsnatur des Expositionsmodells erfasst wird. Solche Beobachtungen werden sowohl aus Feldbegehungen (Kapitel 2) als auch aus frei verfügbaren Datenquellen zur Aktualisierung bestehender (aber veralteter) Expositionsmodelle (Kapitel 3) als echte Beweise gewertet. In diesen beiden Kapiteln wurden Erdbebenszenarien mit parametrisierten Bodenbewegungsfeldern transversal verwendet, um die Rolle solcher epistemischen Unsicherheiten in Bezug auf die Expositionszusammensetzung durch Sensitivitätsanalysen zu untersuchen. Parametrisierte Szenarien seismischer Bodenerschütterungen waren der Gefahreneingang, der verwendet wurde, um die physische Anfälligkeit von Gebäudeportfolios zu untersuchen. Das zweite untersuchte Problem, das sich auf die räumliche Aggregation von Gebäudeexpositionsmodellen bezieht, wurde in zwei entkoppelten Vulnerabilitätskontexten untersucht: durch seismische Bodenerschütterungen durch die Integration von Fernerkundungstechniken (Kapitel 3); und innerhalb eines Multi-Hazard-Kontextes durch Einbeziehung des Auftretens assoziierter Tsunamis (Kapitel 4). Darin wird eine sorgfältige Auswahl der räumlichen Aggregationseinheiten bei gleichzeitigem Streben nach Recheneffizienz und Genauigkeit bei den Risikoschätzungen aufgrund solcher unabhängiger Gefahrenszenarien (d. h. Erdbeben und Tsunami) diskutiert. Daher wird in dieser Arbeit die physische Vulnerabilität von großen Gebäudeportfolios durch Tsunamis durch zwei Hauptrahmen betrachtet: Berücksichtigung und Nichtberücksichtigung der Wechselwirkung auf der Vulnerabilitätsebene, durch aufeinanderfolgende bzw. entkoppelte Gefahrenszenarien, die dann gegenübergestellt wurden. Im Gegensatz zu Kapitel 4, wo keine kumulativen Schäden angesprochen werden, werden in Kapitel 5 Daten und Ansätze, die bereits in früheren Abschnitten generiert wurden, mit einer neuartigen modularen Methode integriert, um letztendlich die wahrscheinlichen Wechselwirkungen auf der Schwachstellenebene beim Aufbau von Portfolios zu untersuchen. Dies wird getestet, indem kumulative Schäden und Verluste nach Erdbeben mit zunehmender Magnitude gefolgt von den jeweiligen Tsunamis bewertet werden. Eine solche neuartige Methode basiert auf der Möglichkeit, bestehende Fragilitätsmodelle innerhalb eines probabilistischen Rahmens wiederzuverwenden. Derselbe Ansatz wird in Kapitel 6 verfolgt, um die wahrscheinlichen kumulativen Schäden zu prognostizieren, denen ein Gebäudebestand ausgesetzt sein wird, der sich in einer vulkanischen Umgebung mit mehreren Gefahren (Aschefall und Lahare) befindet. In diesem Abschnitt wurde besonderes Augenmerk auf die Art und Weise gelegt, wie die prognostizierten Verlustmetriken an lokal exponierte Gemeinden kommuniziert werden. Koexistierende quantitative wissenschaftliche Ansätze (d. h. umfassende Expositionsmodelle; explorative Risikoszenarien mit Einzel- und Mehrfachgefahren) und semiqualitative soziale Risikowahrnehmung (d. h. Grad des Verständnisses, das die exponierten Gemeinschaften über ihr eigenes Risiko haben) wurden gemeinsam berücksichtigt. Eine solche Integration ermöglichte es dieser Arbeit schließlich auch, zur Verbesserung der Bereitschaft, der wissenschaftlichen Verbreitung auf lokaler Ebene sowie zu Technologietransferinitiativen beizutragen. Abschließend wird eine Zusammenfassung dieser These zusammen mit einigen Perspektiven für Verbesserungen und zukünftige Arbeiten präsentiert. KW - exposure KW - multi-hazard KW - risk analysis KW - risk communication KW - uncertainty analysis KW - tsunami risk KW - spatial aggregation KW - seismic risk KW - Expositionsmodellen KW - Multi-Hazard KW - Risikoanalyse KW - Risikokommunikation KW - seismisches Risiko KW - räumliche Aggregation KW - Tsunami-Risiko KW - Unsicherheitsanalyse Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-586140 ER - TY - THES A1 - Zali, Zahra T1 - Volcanic tremor analysis based on advanced signal processing concepts including music information retrieval (MIR) strategies N2 - Volcanoes are one of the Earth’s most dynamic zones and responsible for many changes in our planet. Volcano seismology aims to provide an understanding of the physical processes in volcanic systems and anticipate the style and timing of eruptions by analyzing the seismic records. Volcanic tremor signals are usually observed in the seismic records before or during volcanic eruptions. Their analysis contributes to evaluate the evolving volcanic activity and potentially predict eruptions. Years of continuous seismic monitoring now provide useful information for operational eruption forecasting. The continuously growing amount of seismic recordings, however, poses a challenge for analysis, information extraction, and interpretation, to support timely decision making during volcanic crises. Furthermore, the complexity of eruption processes and precursory activities makes the analysis challenging. A challenge in studying seismic signals of volcanic origin is the coexistence of transient signal swarms and long-lasting volcanic tremor signals. Separating transient events from volcanic tremors can, therefore, contribute to improving our understanding of the underlying physical processes. Some similar issues (data reduction, source separation, extraction, and classification) are addressed in the context of music information retrieval (MIR). The signal characteristics of acoustic and seismic recordings comprise a number of similarities. This thesis is going beyond classical signal analysis techniques usually employed in seismology by exploiting similarities of seismic and acoustic signals and building the information retrieval strategy on the expertise developed in the field of MIR. First, inspired by the idea of harmonic–percussive separation (HPS) in musical signal processing, I have developed a method to extract harmonic volcanic tremor signals and to detect transient events from seismic recordings. This provides a clean tremor signal suitable for tremor investigation along with a characteristic function suitable for earthquake detection. Second, using HPS algorithms, I have developed a noise reduction technique for seismic signals. This method is especially useful for denoising ocean bottom seismometers, which are highly contaminated by noise. The advantage of this method compared to other denoising techniques is that it doesn’t introduce distortion to the broadband earthquake waveforms, which makes it reliable for different applications in passive seismological analysis. Third, to address the challenge of extracting information from high-dimensional data and investigating the complex eruptive phases, I have developed an advanced machine learning model that results in a comprehensive signal processing scheme for volcanic tremors. Using this method seismic signatures of major eruptive phases can be automatically detected. This helps to provide a chronology of the volcanic system. Also, this model is capable to detect weak precursory volcanic tremors prior to the eruption, which could be used as an indicator of imminent eruptive activity. The extracted patterns of seismicity and their temporal variations finally provide an explanation for the transition mechanism between eruptive phases. N2 - Vulkane gehören zu den dynamischsten Zonen der Erde und sind für viele Veränderungen auf unserem Planeten verantwortlich. Die Vulkanseismologie zielt darauf ab, physikalischen Prozesse in Vulkansystemen besser zu verstehen und die Art und den Zeitpunkt von Eruptionen durch die Analyse der seismischen Aufzeichnungen vorherzusagen. Die Signale vulkanischer Tremore werden normalerweise vor oder während Vulkanausbrüchen beobachtet und müssen überwacht werden, um die vulkanische Aktivität zu bewerten. Die Untersuchung vulkanischer Tremore ist ein wichtiger Teil der Vulkanüberwachung, die darauf abzielt, Anzeichen für das Erwachen oder Wiedererwachen von Vulkanen zu erkennen und möglicherweise Ausbrüche vorherzusagen. Mehrere Dekaden kontinuierlicher seismischer Überwachung liefern nützliche Informationen für die operative Eruptionsvorhersage. Die ständig wachsende Menge an seismischen Aufzeichnungen stellt jedoch eine Herausforderung für die Analyse, Informationsextraktion und Interpretation für die zeitnahe Entscheidungsfindung während Vulkankrisen dar. Darüber hinaus erschweren die Komplexität der Eruptionsprozesse und Vorläuferaktivitäten die Analyse. Eine Herausforderung bei der Untersuchung seismischer Signale vulkanischen Ursprungs ist die Koexistenz von transienten Signalschwärmen und lang anhaltenden vulkanischen Tremoren. Die Trennung dieser beiden Signaltypen kann daher dazu beitragen, unser Verständnis der zugrunde liegenden physikalischen Prozesse zu verbessern. Einige ähnliche Probleme (Datenreduktion, Quellentrennung, Extraktion und Klassifizierung) werden im Zusammenhang mit Music Information Retrieval (MIR, dt. Etwa Musik-Informationsabruf) behandelt. Die Signaleigenschaften von akustischen und seismischen Aufzeichnungen weisen eine Reihe von Gemeinsamkeiten auf. Ich gehe über die klassischen Signalanalysetechniken hinaus, die normalerweise in der Seismologie verwendet werden, indem ich die Ähnlichkeiten von seismischen und akustischen Signalen und das Fachwissen aus dem Gebiet der MIR zur Informationsgewinnung nutze. Inspiriert von der Idee der harmonisch-perkussiven Trennung (HPS) in der musikalischen Signalverarbeitung habe ich eine Methode entwickelt, mit der harmonische vulkanische Erschütterungssignale extrahiert und transiente Ereignisse aus seismischen Aufzeichnungen erkannt werden können. Dies liefert ein sauberes Tremorsignal für die Tremoruntersuchung, sowie eine charakteristischen Funktion, die für die Erdbebenerkennung geeignet ist. Weiterhin habe ich unter Verwendung von HPS-Algorithmen eine Rauschunterdrückungstechnik für seismische Signale entwickelt. Diese kann zum Beispiel verwendet werden, um klarere Signale an Meeresbodenseismometern zu erhalten, die sonst durch zu starkes Rauschen überdeckt sind. Der Vorteil dieser Methode im Vergleich zu anderen Denoising-Techniken besteht darin, dass sie keine Verzerrung in der Breitbandantwort der Erdbebenwellen einführt, was sie für verschiedene Anwendungen in der passiven seismologischen Analyse zuverlässiger macht. Um Informationen aus hochdimensionalen Daten zu extrahieren und komplexe Eruptionsphasen zu untersuchen, habe ich ein fortschrittliches maschinelles Lernmodell entwickelt, aus dem ein umfassendes Signalverarbeitungsschema für vulkanische Erschütterungen abgeleitet werden kann. Mit dieser Methode können automatisch seismische Signaturen größerer Eruptionsphasen identifizieren werden. Dies ist nützlich, um die Chronologie eines Vulkansystems zu verstehen. Außerdem ist dieses Modell in der Lage, schwache vulkanische Vorläuferbeben zu erkennen, die als Indikator für bevorstehende Eruptionsaktivität verwendet werden könnten. Basierend auf den extrahierten Seismizitätsmustern und ihren zeitlichen Variationen liefere ich eine Erklärung für den Übergangsmechanismus zwischen verschiedenen Eruptionsphasen. KW - seismic signal processing KW - machine learning KW - volcano seismology KW - music information retrieval KW - noise reduction Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-610866 ER - TY - THES A1 - Arboleda Zapata, Mauricio T1 - Adapted inversion strategies for electrical resistivity data to explore layered near-surface environments N2 - The electrical resistivity tomography (ERT) method is widely used to investigate geological, geotechnical, and hydrogeological problems in inland and aquatic environments (i.e., lakes, rivers, and seas). The objective of the ERT method is to obtain reliable resistivity models of the subsurface that can be interpreted in terms of the subsurface structure and petrophysical properties. The reliability of the resulting resistivity models depends not only on the quality of the acquired data, but also on the employed inversion strategy. Inversion of ERT data results in multiple solutions that explain the measured data equally well. Typical inversion approaches rely on different deterministic (local) strategies that consider different smoothing and damping strategies to stabilize the inversion. However, such strategies suffer from the trade-off of smearing possible sharp subsurface interfaces separating layers with resistivity contrasts of up to several orders of magnitude. When prior information (e.g., from outcrops, boreholes, or other geophysical surveys) suggests sharp resistivity variations, it might be advantageous to adapt the parameterization and inversion strategies to obtain more stable and geologically reliable model solutions. Adaptations to traditional local inversions, for example, by using different structural and/or geostatistical constraints, may help to retrieve sharper model solutions. In addition, layer-based model parameterization in combination with local or global inversion approaches can be used to obtain models with sharp boundaries. In this thesis, I study three typical layered near-surface environments in which prior information is used to adapt 2D inversion strategies to favor layered model solutions. In cooperation with the coauthors of Chapters 2-4, I consider two general strategies. Our first approach uses a layer-based model parameterization and a well-established global inversion strategy to generate ensembles of model solutions and assess uncertainties related to the non-uniqueness of the inverse problem. We apply this method to invert ERT data sets collected in an inland coastal area of northern France (Chapter~2) and offshore of two Arctic regions (Chapter~3). Our second approach consists of using geostatistical regularizations with different correlation lengths. We apply this strategy to a more complex subsurface scenario on a local intermountain alluvial fan in southwestern Germany (Chapter~4). Overall, our inversion approaches allow us to obtain resistivity models that agree with the general geological understanding of the studied field sites. These strategies are rather general and can be applied to various geological environments where a layered subsurface structure is expected. The flexibility of our strategies allows adaptations to invert other kinds of geophysical data sets such as seismic refraction or electromagnetic induction methods, and could be considered for joint inversion approaches. N2 - Die ERT-Methode (Electrical Resistivity Tomography) wird häufig zur Untersuchung geologischer, geotechnischer und hydrogeologischer Probleme im Binnenland und in Gewässern wie beispielsweise Seen, Flüssen oder dem Meer eingesetzt. Das Ziel der ERT-Methode ist es, zuverlässige Widerstandsmodelle des Untergrunds zu erhalten, die in Bezug auf die Struktur des Untergrundes und dessen petrophysikalischer Eigenschaften interpretiert werden können. Die Zuverlässigkeit der resultierenden Widerstandsmodelle hängt nicht nur von der Qualität der erfassten Daten ab, sondern auch von der angewendeten Inversionsstrategie. Die Inversion von ERT-Daten führt zu mehreren Lösungen, die die gemessenen Daten gleich gut erklären. Typische Inversionsansätze basieren auf verschiedenen deterministischen (lokalen) Strategien, die verschiedene Glättungs- und Dämpfungsstrategien berücksichtigen, um die Inversion zu stabilisieren. Diese Strategien haben jedoch den Nachteil, möglicherweise auftretende scharfe Grenzflächen zu verwischen. Es gibt jedoch Szenarien, in denen der Untergrund durch Schichten mit scharfen Grenzflächen gekennzeichnet ist, die Schichten mit hohem Widerstandskontrast (z. B. bis zu mehreren Größenordnungen) voneinander trennen. Wenn Vorwissen (z. B. aus Aufschlüssen, Bohrungen oder anderen geophysikalischen Untersuchungen) auf scharfe Widerstandsvariationen hindeutet, kann es von Vorteil sein, die Parametrisierungs- und Inversionsstrategien anzupassen, um stabilere und geologisch zuverlässige Modelllösungen zu erhalten. Anpassungen traditioneller lokaler Inversionen, beispielweise durch die Verwendung verschiedener struktureller und/oder geostatistischer Bedingungen, können helfen, schärfere Modelllösungen zu erhalten. Zusätzlich kann eine schichtbasierte Modellparametrisierung in Kombination mit lokalen oder globalen Inversionsansätzen verwendet werden, um Modelle mit scharfen Grenzen zu erhalten. In dieser Arbeit habe ich drei typische geschichtete oberflächennahe Umgebungen untersucht, in denen Vorabinformationen verwendet werden, um 2D-Inversionsstrategien so anzupassen, dass geschichtete Untergrundlösungen bevorzugt werden. In Zusammenarbeit mit den Co-Autoren der Kapitel 2-4 habe ich zwei allgemeine Strategien in Betracht gezogen. Unser erster Ansatz verwendet eine schichtbasierte Modellparametrisierung und eine gut etablierte globale Inversionsstrategie. Diese Strategie erzeugt Ensembles von Modelllösungen mithilfe derer die Unsicherheiten im Zusammenhang der Nicht-Eindeutigkeit des inversen Problems bewertet werden können. Wir wenden diese Methode an, um ERT-Datensätze zu invertieren, die in einem Binnenküstengebiet in Nordfrankreich (Kapitel 2) und vor der Küste zweier arktischer Regionen (Kapitel 3) gesammelt wurden. Unser zweiter Ansatz besteht darin, geostatistische Regularisierungen mit unterschiedlichen Korrelationslängen zu verwenden. Wir wenden diese Strategie auf ein komplexeres Untergrundszenario an, das sich auf einen lokalen Schwemmfächer in einem Mittelgebirge im Südwesten Deutshclands umfasst (Kapitel 4). Insgesamt ermöglichen uns unsere Inversionsansätze, Widerstandsmodelle zu erhalten, die mit dem allgemeinen geologischen Verständnis der untersuchten Feldstandorte übereinstimmen. Diese Strategien sind allgemeingültig und können in verschiedenen geologischen Umgebungen angewandt werden, in denen eine geschichtete Struktur des Untergrunds zu erwarten ist. Zudem erlaubt es die Flexibilität unserer Strategien, dass diese an die Inversion anderer geophysikalischer Datensätze wie seismischer Refraktionsmessungen oder elektromagentischer Induktionsverfahren angepasst werden können. Außerdem könnten solche Strategien für gemeinsame Inversionsansätze in Betracht gezogen werden. KW - Near-surface geophysics KW - Electrical resistivity tomography KW - Non-uniqueness KW - Global inversion KW - Particle swarm optimization KW - Ensemble analysis KW - Oberflächennahe Geophysik KW - Tomographie des elektrischen Widerstands KW - Nicht-Einmaligkeit KW - Globale Inversion KW - Partikelschwarm-Optimierung KW - Ensemble-Analyse Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-581357 ER - TY - THES A1 - Schifferle, Lukas T1 - Optical properties of (Mg,Fe)O at high pressure T1 - Optische Eigenschaften von (Mg,Fe)O unter Hochdruck N2 - Large parts of the Earth’s interior are inaccessible to direct observation, yet global geodynamic processes are governed by the physical material properties under extreme pressure and temperature conditions. It is therefore essential to investigate the deep Earth’s physical properties through in-situ laboratory experiments. With this goal in mind, the optical properties of mantle minerals at high pressure offer a unique way to determine a variety of physical properties, in a straight-forward, reproducible, and time-effective manner, thus providing valuable insights into the physical processes of the deep Earth. This thesis focusses on the system Mg-Fe-O, specifically on the optical properties of periclase (MgO) and its iron-bearing variant ferropericlase ((Mg,Fe)O), forming a major planetary building block. The primary objective is to establish links between physical material properties and optical properties. In particular the spin transition in ferropericlase, the second-most abundant phase of the lower mantle, is known to change the physical material properties. Although the spin transition region likely extends down to the core-mantle boundary, the ef-fects of the mixed-spin state, where both high- and low-spin state are present, remains poorly constrained. In the studies presented herein, we show how optical properties are linked to physical properties such as electrical conductivity, radiative thermal conductivity and viscosity. We also show how the optical properties reveal changes in the chemical bonding. Furthermore, we unveil how the chemical bonding, the optical and other physical properties are affected by the iron spin transition. We find opposing trends in the pres-sure dependence of the refractive index of MgO and (Mg,Fe)O. From 1 atm to ~140 GPa, the refractive index of MgO decreases by ~2.4% from 1.737 to 1.696 (±0.017). In contrast, the refractive index of (Mg0.87Fe0.13)O (Fp13) and (Mg0.76Fe0.24)O (Fp24) ferropericlase increases with pressure, likely because Fe Fe interactions between adjacent iron sites hinder a strong decrease of polarizability, as it is observed with increasing density in the case of pure MgO. An analysis of the index dispersion in MgO (decreasing by ~23% from 1 atm to ~103 GPa) reflects a widening of the band gap from ~7.4 eV at 1 atm to ~8.5 (±0.6) eV at ~103 GPa. The index dispersion (between 550 and 870 nm) of Fp13 reveals a decrease by a factor of ~3 over the spin transition range (~44–100 GPa). We show that the electrical band gap of ferropericlase significantly widens up to ~4.7 eV in the mixed spin region, equivalent to an increase by a factor of ~1.7. We propose that this is due to a lower electron mobility between adjacent Fe2+ sites of opposite spin, explaining the previously observed low electrical conductivity in the mixed spin region. From the study of absorbance spectra in Fp13, we show an increasing covalency of the Fe-O bond with pressure for high-spin ferropericlase, whereas in the low-spin state a trend to a more ionic nature of the Fe-O bond is observed, indicating a bond weakening effect of the spin transition. We found that the spin transition is ultimately caused by both an increase of the ligand field-splitting energy and a decreasing spin-pairing energy of high-spin Fe2+. N2 - Geodynamische Prozesse werden von den physikalischen Materialeigenschaften unter den extremen Druck- und Temperaturbedingungen des Erdinneren gesteuert, gerade diese Areale sind aber faktisch nicht für direkte Beobachtungen zugänglich. Umso wichtiger ist es, die physikalischen Eigenschaften unter Bedingungen des Erdinneren zu untersuchen. Mit diesem Ziel vor Augen erlaubt das Studium der optischen Eigenschaften von Mineralen des Erdmantels, eine große Bandbreite an physikalischen Materialeigenschaften, in einer einfachen, reproduzierbaren und effizienten Art und Weise zu bestimmen. Dadurch bieten sich wichtige Einblicke in die physikalischen Prozessen des Erdinneren. Die vorliegende Arbeit konzentriert sich auf das System Mg-Fe-O, im Speziellen auf Periklas (MgO) und seine Eisen-haltige Variante Ferroperiklas ((Mg,Fe)O), ein wichtiger Baustein planetarer Körper. Das Hauptziel der Arbeit besteht darin Verbindungen zwischen optischen Eigenschaften und physikalischen Materialeigenschaften zu finden. Gerade der Spin-Übergang in Ferroperiklas, der zweithäufigsten Phase des unteren Erdmantels, ist dabei von Bedeutung, da damit Veränderungen in den physikalischen Materialeigenschaften einhergehen. Obwohl sich der Spinübergangsbereich vermutlich bis zur Kern-Mantel-Grenze erstreckt, sind die Auswirkungen des gemischten Spin-Zustandes, bei dem sowohl Hoch- als auch Tief-Spin präsent sind, nur unzureichend untersucht. Die hier vorgestellten Studien zeigen, wie optische Eigenschaften mit anderen wichtigen physikalischen Eigenschaften wie elektrischer und thermischer Leitfähigkeit, Viskosität oder auch mit der chemischen Bindung verbunden sind. Daraus lässt sich auch ableiten wie der Spin-Übergang in Ferroperiklas diese Eigenschaften beeinflusst. Von Raumbedingungen bis zu ~140 GPa sinkt der Brechungsindex von MgO um ~2.4 % von 1.737 auf 1.696 (±0.017). Im Gegensatz dazu steigt der Brechungsindex von (Mg0.87Fe0.13)O (Fp13) und (Mg0.76Fe0.24)O (Fp24) Ferroperiklas mit dem Druck an. Dies ist auf Fe-Fe Wechselwirkungen zwischen benachbarten Eisenpositionen zurückzuführen, die eine starke Verringerung der Polarisierbarkeit, wie im Falle von reinem MgO mit zunehmender Dichte, behindern. Eine Analyse der Dispersion des Brechungsindexes von MgO (Abnahme um ~23 % von 1 Atm zu ~103 GPa) offenbart eine Verbreiterung der Bandlücke von ~7.4 eV bei 1 Atm zu ~8.5 (±0.6) eV bei ~103 GPa. Die Messung der Dispersion (zwischen 550 und 870 nm) in Fp13 zeigt eine starke Abnahme über den Bereich des Spin-Überganges (~44–100 GPa) bis zu einem Faktor von ~3. Die Bandlücke nimmt in der Region des gemischten Spin-Zustandes signifikant auf bis zu ~4.7 eV zu (entspricht einer Zunahme um den Faktor ~1.7). Dies deutet auf eine Verringerung der Elektronen-Mobilität zwischen benachbarten Fe2+-Positionen mit unterschiedlichem Spin-Zustand hin, was die bereits in früheren Arbeiten beobachtete Abnahme der elektrischen Leitfähigkeit im Bereich des gemischten Spin-Zustandes erklärt. Absorptionsspektren an Fp13 zeigen eine Druck-bedingte Zunahme der Kovalenz der Fe-O Bindung für Ferroperiklas im Hoch-Spin Zustand, wohingegen Tief-Spin Ferroperiklas einen Trend zu einer mehr ionischen Fe-O Bindung auf-weist, was auf einen Bindungs-schwächenden Effekt des Spin-Wechsels hinweist. Der Übergang von Hoch- zu Tiefspin ist letztlich auf eine Zunahme der Ligandenfeldaufspaltungsenergie sowie eine abnehmende Spinpaarungsenergie von Hoch-Spin Fe2+ zurückzuführen. KW - optical properties KW - optische Eigenschaften KW - high pressure KW - Hochdruck KW - earth mantle KW - Erdmantel KW - diamond anvil cell KW - Diamantstempelzelle KW - ferropericlase KW - Ferroperiklas KW - spectroscopy KW - Spektroskopie Y1 - 2024 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-622166 ER - TY - THES A1 - Antonoglou, Nikolaos T1 - GNSS-based remote sensing: Innovative observation of key hydrological parameters in the Central Andes T1 - GNSS-basierte Fernerkundung: Innovative Beobachtung der wichtigsten hydrologischen Parameter in den zentralen Anden N2 - The Central Andean region is characterized by diverse climate zones with sharp transitions between them. In this work, the area of interest is the South-Central Andes in northwestern Argentina that borders with Bolivia and Chile. The focus is the observation of soil moisture and water vapour with Global Navigation Satellite System (GNSS) remote-sensing methodologies. Because of the rapid temporal and spatial variations of water vapour and moisture circulations, monitoring this part of the hydrological cycle is crucial for understanding the mechanisms that control the local climate. Moreover, GNSS-based techniques have previously shown high potential and are appropriate for further investigation. This study includes both logistic-organization effort and data analysis. As for the prior, three GNSS ground stations were installed in remote locations in northwestern Argentina to acquire observations, where there was no availability of third-party data. The methodological development for the observation of the climate variables of soil moisture and water vapour is independent and relies on different approaches. The soil-moisture estimation with GNSS reflectometry is an approximation that has demonstrated promising results, but it has yet to be operationally employed. Thus, a more advanced algorithm that exploits more observations from multiple satellite constellations was developed using data from two pilot stations in Germany. Additionally, this algorithm was slightly modified and used in a sea-level measurement campaign. Although the objective of this application is not related to monitoring hydrological parameters, its methodology is based on the same principles and helps to evaluate the core algorithm. On the other hand, water-vapour monitoring with GNSS observations is a well-established technique that is utilized operationally. Hence, the scope of this study is conducting a meteorological analysis by examining the along-the-zenith air-moisture levels and introducing indices related to the azimuthal gradient. The results of the experiments indicate higher-quality soil moisture observations with the new algorithm. Furthermore, the analysis using the stations in northwestern Argentina illustrates the limits of this technology because of varying soil conditions and shows future research directions. The water-vapour analysis points out the strong influence of the topography on atmospheric moisture circulation and rainfall generation. Moreover, the GNSS time series allows for the identification of seasonal signatures, and the azimuthal-gradient indices permit the detection of main circulation pathways. N2 - Die Zentralanden sind eine Region, in der verschiedene Klimazonen nur durch kurze Übergänge gekennzeichnet sind. Der geographische Schwerpunkt dieser Arbeit liegt in den südlichen Zentralanden im Grenzgebiet zwischen Argentinien, Bolivien und Chile, und der wissenschaftliche Schwerpunkt ist in der Überwachung der Bodenfeuchtigkeit und des Wasserdampfs mit Fernerkundungsmethoden des Globales Navigationssatellitensystem (Global Navigation Satellite System - GNSS) angesiedelt. Wegen der raschen zeitlichen und räumlichen Schwankungen des Wasserdampfs und den damit häufig verbundenen Niederschlägen und der Feuchtigkeitszirkulation ist die Beobachtung dieses Teils des hydrologischen Zyklus von entscheidender Bedeutung für das Verständnis des lokalen Klimas. Darüber hinaus haben GNSS-gestützte Techniken in anderen Studien bereits ein hohes Potenzial gezeigt, erfordern aber in einigen Bereichen weitere Untersuchungen. Diese Studie umfasst sowohl logistischen Aufwand als auch Datenanalyse. Dazu wurden drei GNSS-Bodenstationen in abgelegenen Orten im Nordwesten Argentiniens installiert, um Beobachtungen zu sammeln, da dort keine externen Daten verfügbar waren. Die methodische Entwicklung für die Beobachtung der Klimavariablen Bodenfeuchtigkeit und Wasserdampfs ist unabhängig voneinander. Die Messung der Bodenfeuchte mit Hilfe der GNSS-Reflektometrie ist eine Annäherung, die vielversprechende Ergebnisse erbracht hat, aber bisher noch nicht operationell eingesetzt wurde. Daher wurde ein fortschrittlicherer Algorithmus entwickelt, der Beobachtungen von mehreren Satellitenkonstellationen nutzt und unter anderem Daten von zwei Pilotstationen in Deutschland verwendet. Außerdem wurde dieser Algorithmus leicht modifiziert und in einer Meeresspiegelmesskampagne eingesetzt. Obwohl diese Andwendung nicht direkt mit der Überwachung hydrologischer Parameter zusammenhängt, basiert die Methodik auf denselben Prinzipien und hilft bei der Bewertung des entwickelten Algorithmus. Auf der anderen Seite ist die Überwachung des Wasserdampfs mit GNSS-Beobachtungen eine anerkannte Technik, die in der Praxis bereits seit mehreren Jahren eingesetzt wird. Diese Studie befasst sich daher mit der Durchführung einer meteorologischen Analyse der Luftfeuchtigkeitswerte entlang des Zenits und der Entwicklung von klimatischen Indizes, die sich auf den azimutalen Gradienten beziehen. Die Ergebnisse der Experimente zeigen, dass die Qualität der Bodenfeuchtebeobachtungen mit dem neuen Algorithmus vielversprechend und besser sind. Darüber hinaus zeigt die Analyse anhand der Stationen im nordwesten Argentiniens die Grenzen dieser Technologie aufgrund der sehr unterschiedlichen Bodenbedingungen auf und gibt mögliche zukünftige Forschungsrichtung an. Die Wasserdampfanalyse verdeutlicht den Einfluss der Topographie auf die Luftfeuchtigkeit und der Regenmenge. Außerdem ermöglichen die GNSS-Zeitreihen die Identifizierung der jahreszeitlichen Signaturen, und Messungen der azimutal Gradienten erlauben die Erkennung der wichtigsten Zirkulationswege. KW - remote sensing KW - GNSS KW - GPS KW - water vapour KW - soil moisture KW - Central Andes KW - zentrale Anden KW - globales Navigationssatellitensystem KW - globales Positionsbestimmungssystem KW - Fernerkundung KW - Bodenfeuchtigkeit KW - Wasserdampf Y1 - 2024 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-628256 ER - TY - THES A1 - Windirsch-Woiwode, Torben T1 - Permafrost carbon stabilisation by recreating a herbivore-driven ecosystem T1 - Stabilisierung von Permafrostkohlenstoff durch die Wiedereinführung eines Herbivor-geprägten Ökosystems N2 - With Arctic ground as a huge and temperature-sensitive carbon reservoir, maintaining low ground temperatures and frozen conditions to prevent further carbon emissions that contrib-ute to global climate warming is a key element in humankind’s fight to maintain habitable con-ditions on earth. Former studies showed that during the late Pleistocene, Arctic ground condi-tions were generally colder and more stable as the result of an ecosystem dominated by large herbivorous mammals and vast extents of graminoid vegetation – the mammoth steppe. Characterised by high plant productivity (grassland) and low ground insulation due to animal-caused compression and removal of snow, this ecosystem enabled deep permafrost aggrad-ation. Now, with tundra and shrub vegetation common in the terrestrial Arctic, these effects are not in place anymore. However, it appears to be possible to recreate this ecosystem local-ly by artificially increasing animal numbers, and hence keep Arctic ground cold to reduce or-ganic matter decomposition and carbon release into the atmosphere. By measuring thaw depth, total organic carbon and total nitrogen content, stable carbon iso-tope ratio, radiocarbon age, n-alkane and alcohol characteristics and assessing dominant vegetation types along grazing intensity transects in two contrasting Arctic areas, it was found that recreating conditions locally, similar to the mammoth steppe, seems to be possible. For permafrost-affected soil, it was shown that intensive grazing in direct comparison to non-grazed areas reduces active layer depth and leads to higher TOC contents in the active layer soil. For soil only frozen on top in winter, an increase of TOC with grazing intensity could not be found, most likely because of confounding factors such as vertical water and carbon movement, which is not possible with an impermeable layer in permafrost. In both areas, high animal activity led to a vegetation transformation towards species-poor graminoid-dominated landscapes with less shrubs. Lipid biomarker analysis revealed that, even though the available organic material is different between the study areas, in both permafrost-affected and sea-sonally frozen soils the organic material in sites affected by high animal activity was less de-composed than under less intensive grazing pressure. In conclusion, high animal activity af-fects decomposition processes in Arctic soils and the ground thermal regime, visible from reduced active layer depth in permafrost areas. Therefore, grazing management might be utilised to locally stabilise permafrost and reduce Arctic carbon emissions in the future, but is likely not scalable to the entire permafrost region. N2 - Mit dem arktischen Boden als riesigem und temperatursensiblen Kohlenstoffspeicher ist die Aufrechterhaltung niedriger Bodentemperaturen und gefrorener Bedingungen zur Verhinde-rung weiterer Kohlenstoffemissionen, die zum globalen Klimawandel beitragen, ein Schlüs-selelement im Kampf der Menschheit, die Erde weiterhin bewohnbar zu halten. Vorangehen-de Studien ergaben, dass die Bodenbedingungen in der Arktis während des späten Pleisto-zäns im Allgemeinen kälter und dadurch stabiler waren, als Ergebnis eines Ökosystems, das von großen pflanzenfressenden Säugetieren und weiten Flächen grasartiger Vegetation do-miniert wurde - der Mammutsteppe. Gekennzeichnet durch hohe Pflanzenproduktivität (Gras-land) und geringe Bodenisolierung aufgrund von Kompression und Schneeräumung durch Tiere, ermöglichte dieses Ökosystem eine tiefreichende Entwicklung des Permafrosts. Heut-zutage, mit der vorherrschenden Tundra- und Strauchvegetation in der Arktis, sind diese Ef-fekte nicht mehr präsent. Es scheint aber möglich, dieses Ökosystem lokal durch künstliche Erhöhung der Tierbestände nachzubilden und somit den arktischen Boden kühl zu halten, um den Abbau von organischem Material und die Freisetzung von Kohlenstoff in die Atmosphäre zu verringern. Durch Messungen der Auftautiefe, des Gesamtgehalts des organischen Kohlenstoffs und Stickstoffs, des stabilen Kohlenstoff-Isotopenverhältnisses, des Radiocarbonalters, der n-Alkan- und Alkoholcharakteristika sowie durch Bestimmung der vorherrschenden Vegetati-onstypen entlang von Beweidungsgradienten in zwei unterschiedlichen arktischen Gebieten habe ich festgestellt, dass die Schaffung ähnlicher Bedingungen wie in der Mammutsteppe möglich sein könnte. Für durch Permafrost beeinflusste Böden konnte ich zeigen, dass eine intensive Beweidung im direkten Vergleich mit unbeweideten Gebieten die Tiefe der Auftau-schicht verringert und zu höheren Gehalten an organischem Kohlenstoff im oberen Bodenbe-reich führt. Für im Winter nur oberflächlich gefrorene Böden konnte kein Anstieg des organi-schen Kohlenstoffgehalts mit zunehmender Beweidungsintensität festgestellt werden, höchstwahrscheinlich aufgrund von Störfaktoren wie vertikalen Wasser- und Kohlenstoffbe-wegungen, die nicht durch eine undurchlässige Schicht wie beim Permafrost begrenzt sind. In beiden Gebieten führte eine hohe Tieraktivität zu einer Umwandlung der Vegetation hin zu artenarmen, von Gräsern dominierten Landschaften mit weniger Sträuchern. Die Analyse von Lipid-Biomarkern ergab, dass das verfügbare organische Material zwar zwischen den Unter-suchungsgebieten unterschiedlich war, aber sowohl in Permafrostgebieten als auch in saiso-nal gefrorenen Böden in Bereichen mit hoher Tieraktivität weniger stark zersetzt war als unter geringerer Beweidungsintensität. Zusammenfassend beeinflusst eine hohe Tieraktivität die Zersetzungsvorgänge in arktischen Böden und das thermische Regime des Bodens, was sich in einer reduzierten Tiefe der Auftauschicht in Permafrostgebieten widerspiegelt. Daher könn-te das Beweidungsmanagement in Zukunft aktiv eingesetzt werden, um den Permafrost lokal zu stabilisieren und gefroren zu halten sowie die Kohlenstoffemissionen in der Arktis zu ver-ringern. Aufgrund der Größe der Fläche, die in der terrestrischen Arktis von Permafrost be-einflusst ist, wird ein solches Beweidungsmanagement aber nicht als Maßnahme auf die ge-samte Permafrostregion ausgedehnt werden können. KW - permafrost KW - carbon KW - climate change KW - grazing KW - Arctic KW - Arktis KW - Kohlenstoff KW - Klimawandel KW - Beweidung KW - Permafrost Y1 - 2024 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-624240 ER - TY - THES A1 - Müller, Daniel T1 - Identification and monitoring of structures, controls, and evolution dynamics of hydrothermal systems and associated alteration through high-resolution remote sensing and in situ analysis T1 - Identifizierung und Überwachung von Strukturen, Kontrollmechanismen und Entwicklungsdynamiken hydrothermaler Systeme und damit verbundener Alteration durch hochauflösende Fernerkundung und In-situ-Analyse N2 - Volcanic hydrothermal systems are an integral part of most volcanoes and typically involve a heat source, adequate fluid supply, and fracture or pore systems through which the fluids can circulate within the volcanic edifice. Associated with this are subtle but powerful processes that can significantly influence the evolution of volcanic activity or the stability of the near-surface volcanic system through mechanical weakening, permeability reduction, and sealing of the affected volcanic rock. These processes are well constrained for rock samples by laboratory analyses but are still difficult to extrapolate and evaluate at the scale of an entire volcano. Advances in unmanned aircraft systems (UAS), sensor technology, and photogrammetric processing routines now allow us to image volcanic surfaces at the centimeter scale and thus study volcanic hydrothermal systems in great detail. This thesis aims to explore the potential of UAS approaches for studying the structures, processes, and dynamics of volcanic hydrothermal systems but also to develop methodological approaches to uncover secondary information hidden in the data, capable of indicating spatiotemporal dynamics or potentially critical developments associated with hydrothermal alteration. To accomplish this, the thesis describes the investigation of two near-surface volcanic hydrothermal systems, the El Tatio geyser field in Chile and the fumarole field of La Fossa di Vulcano (Italy), both of which are among the best-studied sites of their kind. Through image analysis, statistical, and spatial analyses we have been able to provide the most detailed structural images of both study sites to date, with new insights into the driving forces of such systems but also revealing new potential controls, which are summarized in conceptual site-specific models. Furthermore, the thesis explores methodological remote sensing approaches to detect, classify and constrain hydrothermal alteration and surface degassing from UAS-derived data, evaluated them by mineralogical and chemical ground-truthing, and compares the alteration pattern with the present-day degassing activity. A significant contribution of the often neglected diffuse degassing activity to the total amount of degassing is revealed and constrains secondary processes and dynamics associated with hydrothermal alteration that lead to potentially critical developments like surface sealing. The results and methods used provide new approaches for alteration research, for the monitoring of degassing and alteration effects, and for thermal monitoring of fumarole fields, with the potential to be incorporated into volcano monitoring routines. N2 - Vulkanische Hydrothermalsysteme sind ein integraler Bestandteil vieler Vulkane und erfordern allgemein eine ausreichende Wärmequelle, eine ausreichende Flüssigkeitszufuhr und Kluft- oder Porensysteme, durch die Flüssigkeiten zirkulieren können. Damit verbunden sind subtile, aber wirksame Prozesse, welche die Entwicklung der vulkanischen Aktivität oder die Stabilität des oberflächennahen Vulkansystems durch mechanische Schwächung, Verringerung der Durchlässigkeit und Versiegelung des betroffenen vulkanischen Gesteins erheblich beeinflussen können. Solche Prozesse sind für Gesteinsproben durch Laboranalysen gut definiert, aber es ist immer noch schwierig sie auf der Skala eines ganzen Vulkans zu bewerten. Fortschritte bei unbemannten Flugsystemen (UAS), Sensortechnologie und photogrammetrischen Prozessierungs-routinen ermöglichen es uns heute, Vulkanoberflächen im Zentimeterbereich abzubilden und damit vulkanische Hydrothermalsysteme sehr detailliert zu untersuchen. Ziel dieser Arbeit ist es, deren Potenzial für die Untersuchung der Strukturen, Prozesse und Dynamik solcher Systeme zu erforschen, aber auch methodische Ansätze zu finden, um in den Daten verborgene Sekundärinformationen zu analysieren, die auf raum-zeitliche Dynamiken oder potenziell kritische Entwicklungen im Zusammenhang mit hydrothermaler Alteration hinweisen können. Wir haben zwei oberflächennahe vulkanische hydrothermale Systeme analysiert, das Geysirfeld von El Tatio in Chile und das Fumarolenfeld von La Fossa di Vulcano in Italien, die beide zu den am besten untersuchten Systemen ihrer Art gehören. Durch Bildanalyse, statistische und räumliche Analysen konnten wir das bisher detaillierteste Abbild des strukturellen Aufbaus beider Standorte erstellen, neue Einblicke in die oft faszinierende Systematik solcher Systeme geben, aber auch neue potenzielle Kontrollfaktoren aufzeigen. Die Ergebnisse werden in konzeptionellen Modellen zusammengefasst. Darüber hinaus haben wir methodische Ansätze der Fernerkundung zur Erkennung, Klassifizierung und räumlichen Eingrenzung hydrothermaler Alteration aus UAS-Daten untersucht, durch mineralogisches und chemisches „Ground-Truthing“ bewertet und die Alterationsmuster mit der aktuellen Entgasungsaktivität verglichen. Wir zeigen dass die oft nicht berücksichtigte diffuse Aktivität einen signifikanten Beitrag zur Gesamtaktivität liefert, aber auch Bereiche in denen sekundäre Prozesse hydrothermaler Alteration scheinbar zu potenziell kritischen Entwicklungen wie Oberflächenversiegelung führen. Die Ergebnisse und die verwendeten Methoden bieten neue Ansätze für die Alterationsforschung, für die Überwachung von Entgasungs- und Alterationseffekten und für die thermische Überwachung von Fumarolenfeldern, und haben Potential in Vulkanüberwachungsroutinen integriert zu werden. KW - volcano remote sensing KW - Fernerkundung an Vulkanen KW - volcanic hydrothermal systems KW - vulkanische Entgasungs-und Hydrothermalsysteme KW - Vulkanüberwachung KW - hydrothermale Alteration KW - vulkanische Entgasungssysteme KW - Drohnen-Fernerkundung Y1 - 2024 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-626683 ER - TY - THES A1 - Lilienkamp, Henning T1 - Enhanced computational approaches for data-driven characterization of earthquake ground motion and rapid earthquake impact assessment T1 - Fortgeschrittene Berechnungsansätze für die datengestützte Charakterisierung von Erdbeben-Bodenbewegungen und die schnelle Einschätzung von Erdbebenauswirkungen N2 - Rapidly growing seismic and macroseismic databases and simplified access to advanced machine learning methods have in recent years opened up vast opportunities to address challenges in engineering and strong motion seismology from novel, datacentric perspectives. In this thesis, I explore the opportunities of such perspectives for the tasks of ground motion modeling and rapid earthquake impact assessment, tasks with major implications for long-term earthquake disaster mitigation. In my first study, I utilize the rich strong motion database from the Kanto basin, Japan, and apply the U-Net artificial neural network architecture to develop a deep learning based ground motion model. The operational prototype provides statistical estimates of expected ground shaking, given descriptions of a specific earthquake source, wave propagation paths, and geophysical site conditions. The U-Net interprets ground motion data in its spatial context, potentially taking into account, for example, the geological properties in the vicinity of observation sites. Predictions of ground motion intensity are thereby calibrated to individual observation sites and earthquake locations. The second study addresses the explicit incorporation of rupture forward directivity into ground motion modeling. Incorporation of this phenomenon, causing strong, pulse like ground shaking in the vicinity of earthquake sources, is usually associated with an intolerable increase in computational demand during probabilistic seismic hazard analysis (PSHA) calculations. I suggest an approach in which I utilize an artificial neural network to efficiently approximate the average, directivity-related adjustment to ground motion predictions for earthquake ruptures from the 2022 New Zealand National Seismic Hazard Model. The practical implementation in an actual PSHA calculation demonstrates the efficiency and operational readiness of my model. In a follow-up study, I present a proof of concept for an alternative strategy in which I target the generalizing applicability to ruptures other than those from the New Zealand National Seismic Hazard Model. In the third study, I address the usability of pseudo-intensity reports obtained from macroseismic observations by non-expert citizens for rapid impact assessment. I demonstrate that the statistical properties of pseudo-intensity collections describing the intensity of shaking are correlated with the societal impact of earthquakes. In a second step, I develop a probabilistic model that, within minutes of an event, quantifies the probability of an earthquake to cause considerable societal impact. Under certain conditions, such a quick and preliminary method might be useful to support decision makers in their efforts to organize auxiliary measures for earthquake disaster response while results from more elaborate impact assessment frameworks are not yet available. The application of machine learning methods to datasets that only partially reveal characteristics of Big Data, qualify the majority of results obtained in this thesis as explorative insights rather than ready-to-use solutions to real world problems. The practical usefulness of this work will be better assessed in the future by applying the approaches developed to growing and increasingly complex data sets. N2 - Das rapide Wachstum seismischer und makroseismischer Datenbanken und der vereinfachte Zugang zu fortschrittlichen Methoden aus dem Bereich des maschinellen Lernens haben in den letzen Jahren die datenfokussierte Betrachtung von Fragestellungen in der Seismologie ermöglicht. In dieser Arbeit erforsche ich das Potenzial solcher Betrachtungsweisen im Hinblick auf die Modellierung erdbebenbedingter Bodenerschütterungen und der raschen Einschätzung von gesellschaftlichen Erdbebenauswirkungen, Disziplinen von erheblicher Bedeutung für den langfristigen Erdbebenkatastrophenschutz in seismisch aktiven Regionen. In meiner ersten Studie nutze ich die Vielzahl an Bodenbewegungsdaten aus der Kanto Region in Japan, sowie eine spezielle neuronale Netzwerkarchitektur (U-Net) um ein Bodenbewegungsmodell zu entwickeln. Der einsatzbereite Prototyp liefert auf Basis der Charakterisierung von Erdbebenherden, Wellenausbreitungspfaden und Bodenbeschaffenheiten statistische Schätzungen der zu erwartenden Bodenerschütterungen. Das U-Net interpretiert Bodenbewegungsdaten im räumlichen Kontext, sodass etwa die geologischen Beschaffenheiten in der Umgebung von Messstationen mit einbezogen werden können. Auch die absoluten Koordinaten von Erdbebenherden und Messstationen werden berücksichtigt. Die zweite Studie behandelt die explizite Berücksichtigung richtungsabhängiger Verstärkungseffekte in der Bodenbewegungsmodellierung. Obwohl solche Effekte starke, impulsartige Erschütterungen in der Nähe von Erdbebenherden erzeugen, die eine erhebliche seismische Beanspruchung von Gebäuden darstellen, wird deren explizite Modellierung in der seismischen Gefährdungsabschätzung aufgrund des nicht vertretbaren Rechenaufwandes ausgelassen. Mit meinem, auf einem neuronalen Netzwerk basierenden, Ansatz schlage ich eine Methode vor, umdieses Vorhaben effizient für Erdbebenszenarien aus dem neuseeländischen seismischen Gefährdungsmodell für 2022 (NSHM) umzusetzen. Die Implementierung in einer seismischen Gefährdungsrechnung unterstreicht die Praktikabilität meines Modells. In einer anschließenden Machbarkeitsstudie untersuche ich einen alternativen Ansatz der auf die Anwendbarkeit auf beliebige Erdbebeszenarien abzielt. Die abschließende dritte Studie befasst sich mit dem potenziellen Nutzen der von makroseismischen Beobachtungen abgeleiteten pseudo-Erschütterungsintensitäten für die rasche Abschätzung von gesellschaftlichen Erdbebenauswirkungen. Ich zeige, dass sich aus den Merkmalen solcher Daten Schlussfolgerungen über die gesellschaftlichen Folgen eines Erdbebens ableiten lassen. Basierend darauf formuliere ich ein statistisches Modell, welches innerhalb weniger Minuten nach einem Erdbeben die Wahrscheinlichkeit für das Auftreten beachtlicher gesellschaftlicher Auswirkungen liefert. Ich komme zu dem Schluss, dass ein solches Modell, unter bestimmten Bedingungen, hilfreich sein könnte, um EntscheidungsträgerInnen in ihren Bestrebungen Hilfsmaßnahmen zu organisieren zu unterstützen. Die Anwendung von Methoden des maschinellen Lernens auf Datensätze die sich nur begrenzt als Big Data charakterisieren lassen, qualifizieren die Mehrheit der Ergebnisse dieser Arbeit als explorative Einblicke und weniger als einsatzbereite Lösungen für praktische Fragestellungen. Der praktische Nutzen dieser Arbeit wird sich in erst in Zukunft an der Anwendung der erarbeiteten Ansätze auf wachsende und zunehmend komplexe Datensätze final abschätzen lassen. KW - seismology KW - machine learning KW - deep learning KW - ground motion modeling KW - seismic hazard KW - rapid earthquake impact assessment KW - geophysics KW - Deep Learning KW - Geophysik KW - Bodenbewegungsmodellierung KW - maschinelles Lernen KW - schnelle Einschätzung von Erdbebenauswirkungen KW - seismische Gefährdung KW - Seismologie Y1 - 2024 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-631954 ER - TY - THES A1 - Sharma, Shubham T1 - Integrated approaches to earthquake forecasting T1 - Integrierte Ansätze zur Vorhersage von Erdbeben BT - insights from Coulomb stress, seismotectonics, and aftershock sequences BT - Erkenntnisse aus Coulomb-Stress, Seismotektonik und Nachbebenfolgen N2 - A comprehensive study on seismic hazard and earthquake triggering is crucial for effective mitigation of earthquake risks. The destructive nature of earthquakes motivates researchers to work on forecasting despite the apparent randomness of the earthquake occurrences. Understanding their underlying mechanisms and patterns is vital, given their potential for widespread devastation and loss of life. This thesis combines methodologies, including Coulomb stress calculations and aftershock analysis, to shed light on earthquake complexities, ultimately enhancing seismic hazard assessment. The Coulomb failure stress (CFS) criterion is widely used to predict the spatial distributions of aftershocks following large earthquakes. However, uncertainties associated with CFS calculations arise from non-unique slip inversions and unknown fault networks, particularly due to the choice of the assumed aftershocks (receiver) mechanisms. Recent studies have proposed alternative stress quantities and deep neural network approaches as superior to CFS with predefined receiver mechanisms. To challenge these propositions, I utilized 289 slip inversions from the SRCMOD database to calculate more realistic CFS values for a layered-half space and variable receiver mechanisms. The analysis also investigates the impact of magnitude cutoff, grid size variation, and aftershock duration on the ranking of stress metrics using receiver operating characteristic (ROC) analysis. Results reveal the performance of stress metrics significantly improves after accounting for receiver variability and for larger aftershocks and shorter time periods, without altering the relative ranking of the different stress metrics. To corroborate Coulomb stress calculations with the findings of earthquake source studies in more detail, I studied the source properties of the 2005 Kashmir earthquake and its aftershocks, aiming to unravel the seismotectonics of the NW Himalayan syntaxis. I simultaneously relocated the mainshock and its largest aftershocks using phase data, followed by a comprehensive analysis of Coulomb stress changes on the aftershock planes. By computing the Coulomb failure stress changes on the aftershock faults, I found that all large aftershocks lie in regions of positive stress change, indicating triggering by either co-seismic or post-seismic slip on the mainshock fault. Finally, I investigated the relationship between mainshock-induced stress changes and associated seismicity parameters, in particular those of the frequency-magnitude (Gutenberg-Richter) distribution and the temporal aftershock decay (Omori-Utsu law). For that purpose, I used my global data set of 127 mainshock-aftershock sequences with the calculated Coulomb Stress (ΔCFS) and the alternative receiver-independent stress metrics in the vicinity of the mainshocks and analyzed the aftershocks properties depend on the stress values. Surprisingly, the results show a clear positive correlation between the Gutenberg-Richter b-value and induced stress, contrary to expectations from laboratory experiments. This observation highlights the significance of structural heterogeneity and strength variations in seismicity patterns. Furthermore, the study demonstrates that aftershock productivity increases nonlinearly with stress, while the Omori-Utsu parameters c and p systematically decrease with increasing stress changes. These partly unexpected findings have significant implications for future estimations of aftershock hazard. The findings in this thesis provides valuable insights into earthquake triggering mechanisms by examining the relationship between stress changes and aftershock occurrence. The results contribute to improved understanding of earthquake behavior and can aid in the development of more accurate probabilistic-seismic hazard forecasts and risk reduction strategies. N2 - Ein umfassendes Verständnis der seismischen Gefahr und Erdbebenauslösung ist wichtig für eine Minderung von Erdbebenrisiken. Die zerstörerische Natur von Erdbeben motiviert Forscher dazu, trotz der scheinbaren Zufälligkeit der Erdbebenereignisse an Vorhersagen zu arbeiten. Das Verständnis der den Beben zugrunde liegenden Mechanismen und Muster ist angesichts ihres Potenzials für weitreichende Verwüstung und den Verlust von Menschenleben von entscheidender Bedeutung. Diese Arbeit kombiniert Methoden, einschließlich der Berechnung der Coulombschen Spannung und der Analyse von Nachbeben, um die Komplexitäten von Erdbeben besser zu verstehen und letztendlich die Bewertung der seismischen Gefahr zu verbessern. Das Coulomb Spannungskriterium (CFS) wird oft verwendet, um die räumliche Verteilung von Nachbeben nach großen Erdbeben vorherzusagen. Jedoch ergeben sich Unsicherheiten bei der Berechnung von CFS aus nicht eindeutigen slip-inversion und der unbekannten Störungsnetzwerken, insbesondere aufgrund der Unsicherheit bezüglich der Nachbebenmechanismen (Empfänger). Neueste Studien deuten darauf hin dass alternative Spannungsgrößen und Deep-Learning-Ansätze gegenüber CFS mit vordefinierten Empfängermechanismen. Um diese Ergebnisse zu hinterfragen, habe ich 289 Slip-inversion uberlegensind aus der SRCMOD-Datenbank verwendet, um realistischere CFS-Werte für einen geschichteten Halbraum und variable Empfängermechanismen zu berechnen. Dabei habe ich auch den Einfluss von Magnitudenschwellenwerten, Gittergrößenvariationen und der Nachbeben-Dauer auf die vorhersagemöglichkeiten der Spannungsmetriken unter Verwendung der ROC-Analyse (Receiver Operating Characteristic) untersucht. Die Ergebnisse zeigen, dass die berudzsidtizung von variablen Empfangermechanism und größere Nachbeben und kürzere Zeiträume die vorhersagekraft steigern, wobei die relative Rangfolge der verschiedenen Spannungsmetriken nicht geändert wird. Um die Coulomb Spannungsberechnungen genauer mit den Ergebnissen von Erdbebenstudien abzugleichen, habe ich die Quelleneigenschaften des Erdbebens von Kaschmir aus dem Jahr 2005 und seiner Nachbeben mit dem ziel, die Seismotektonik des NW-Himalaya Syntaxis zu entschlüsseln, detailliert untersucht. Ich habe gleichzeitig das Hauptbeben und seine größten Nachbeben unter Verwendung von seismischen Phaseneinsetzen relokalisiert und anschließend eine umfassende Analyse der Coulomb Spannungsänderungen auf den Bruchflächen der Nachbeben durchgeführt. Durch die Berechnung der Coulomb Spannungsänderungen an den während der Nachbeben aktivierten Störungen konnte ich herausfinden, dass alle großen Nachbeben in Regionen mit positiven Spannungsänderungen liegen, was auf eine Auslösung durch entweder ko-seismische oder post-seismische Verschiebungen des Hauptbebens hinweist. Schließlich habe ich die Beziehung zwischen den durch Hauptbeben verursachten Spannungsänderungen und den damit verbundenen seismischen Parametern untersucht, insbesondere denen der Häufigkeits-Magnituden (Gutenberg-Richter) Verteilung und des zeitlichen Nachbebenabklingens (Omori-Utsu-Gesetz). Zu diesem Zweck habe ich meinen globalen Datensatz von 127 Hauptbeben-Nachbeben-Sequenzen mit den in der Umgebung der Hauptbeben berechneten Coulomb Spannungen ($\Delta$CFS) zusammen mit den alternativen, empfänger-unabhängigen Spannungsmetriken, verwendet und die Eigenschaften in Abhängigkeit der Spannungswerte analysiert. Überraschenderweise zeigen die Ergebnisse eine klar positive Korrelation zwischen dem $b$-Wert der Gutenberg-Richter-Verteilung und der induzierten Spannung, was im Kontrast zu den Erwartungen aus Laborexperimenten steht. Diese Beobachtung unterstreicht die Bedeutung struktureller Heterogenitäten und Festigkeitsvariationen in seismischen Mustern. Darüber hinaus zeigt die Studie, dass die Anzahl von Nachbeben nichtlinear mit der Spannung zunimmt, während die Omori-Utsu-Parameter $c$ und $p$ systematisch mit zunehmenden Spannungsänderungen abnehmen. Diese teilweise unerwarteten Ergebnisse haben bedeutende Auswirkungen auf zukünftige Abschätzungen der Nachbebengefahr. Die Ergebnisse dieser Arbeit liefern wertvolle Einblicke in die Mechanismen der Erdbebenauslösung, indem sie die Beziehung zwischen Spannungsänderungen und dem Auftreten von Nachbeben untersuchen. Die Ergebnisse tragen zu einem besseren Verständnis des Verhaltens von Erdbeben bei und können bei der Entwicklung genauerer probabilistischer, seismischer Gefahreneinschätzungen und Risikominderungsstrategien helfen. KW - earthquake KW - forecasting KW - hazards KW - seismology KW - Erdbeben KW - Vorhersage KW - Gefahren KW - Seismologie Y1 - 2024 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-636125 ER - TY - THES A1 - Koyan, Philipp T1 - 3D attribute analysis and classification to interpret ground-penetrating radar (GPR) data collected across sedimentary environments: Synthetic studies and field examples T1 - 3D Attributanalyse und -klassifizierung zur Interpretation von Georadar-Daten in sedimentären Umgebungen: Synthetische Studien und Feldbeispiele N2 - Die Untersuchung des oberflächennahen Untergrundes erfolgt heutzutage bei Frage- stellungen aus den Bereichen des Bauwesens, der Archäologie oder der Geologie und Hydrologie oft mittels zerstörungsfreier beziehungsweise zerstörungsarmer Methoden der angewandten Geophysik. Ein Bereich, der eine immer zentralere Rolle in Forschung und Ingenieurwesen einnimmt, ist die Untersuchung von sedimentären Umgebungen, zum Beispiel zur Charakterisierung oberflächennaher Grundwassersysteme. Ein in diesem Kontext häufig eingesetztes Verfahren ist das des Georadars (oftmals GPR - aus dem Englischen ground-penetrating radar). Dabei werden kurze elektromagnetische Impulse von einer Antenne in den Untergrund ausgesendet, welche dort wiederum an Kontrasten der elektromagnetischen Eigenschaften (wie zum Beispiel an der Grundwasseroberfläche) reflektiert, gebrochen oder gestreut werden. Eine Empfangsantenne zeichnet diese Signale in Form derer Amplituden und Laufzeiten auf. Eine Analyse dieser aufgezeichneten Signale ermöglicht Aussagen über den Untergrund, beispielsweise über die Tiefenlage der Grundwasseroberfläche oder die Lagerung und Charakteristika oberflächennaher Sedimentschichten. Dank des hohen Auflösungsvermögens der GPR-Methode sowie stetiger technologischer Entwicklungen erfolgt heutzutage die Aufzeichnung von GPR- Daten immer häufiger in 3D. Trotz des hohen zeitlichen und technischen Aufwandes für die Datenakquisition und -bearbeitung werden die resultierenden 3D-Datensätze, welche den Untergrund hochauflösend abbilden, typischerweise von Hand interpretiert. Dies ist in der Regel ein äußerst zeitaufwendiger Analyseschritt. Daher werden oft repräsentative 2D-Schnitte aus dem 3D-Datensatz gewählt, in denen markante Reflektionsstrukuren markiert werden. Aus diesen Strukturen werden dann sich ähnelnde Bereiche im Untergrund als so genannte Radar-Fazies zusammengefasst. Die anhand von 2D-Schnitten erlangten Resultate werden dann als repräsentativ für die gesamte untersuchte Fläche angesehen. In dieser Form durchgeführte Interpretationen sind folglich oft unvollständig sowie zudem in hohem Maße von der Expertise der Interpretierenden abhängig und daher in der Regel nicht reproduzierbar. Eine vielversprechende Alternative beziehungsweise Ergänzung zur manuellen In- terpretation ist die Verwendung von so genannten GPR-Attributen. Dabei werden nicht die aufgezeichneten Daten selbst, sondern daraus abgeleitete Größen, welche die markanten Reflexionsstrukturen in 3D charakterisieren, zur Interpretation herangezogen. In dieser Arbeit wird anhand verschiedener Feld- und Modelldatensätze untersucht, welche Attribute sich dafür insbesondere eignen. Zudem zeigt diese Arbeit, wie ausgewählte Attribute mittels spezieller Bearbeitungs- und Klassifizierungsmethoden zur Erstellung von 3D-Faziesmodellen genutzt werden können. Dank der Möglichkeit der Erstellung so genannter attributbasierter 3D-GPR-Faziesmodelle können zukünftige Interpretationen zu gewissen Teilen automatisiert und somit effizienter durchgeführt werden. Weiterhin beschreiben die so erhaltenen Resultate den untersuchten Untergrund in reproduzierbarer Art und Weise sowie umfänglicher als es bisher mittels manueller Interpretationsmethoden typischerweise möglich war. N2 - Today, near-surface investigations are frequently conducted using non-destructive or minimally invasive methods of applied geophysics, particularly in the fields of civil engineering, archaeology, geology, and hydrology. One field that plays an increasingly central role in research and engineering is the examination of sedimentary environments, for example, for characterizing near-surface groundwater systems. A commonly employed method in this context is ground-penetrating radar (GPR). In this technique, short electromagnetic pulses are emitted into the subsurface by an antenna, which are then reflected, refracted, or scattered at contrasts in electromagnetic properties (such as the water table). A receiving antenna records these signals in terms of their amplitudes and travel times. Analysis of the recorded signals allows for inferences about the subsurface, such as the depth of the groundwater table or the composition and characteristics of near-surface sediment layers. Due to the high resolution of the GPR method and continuous technological advancements, GPR data acquisition is increasingly performed in three-dimensional (3D) fashion today. Despite the considerable temporal and technical efforts involved in data acquisition and processing, the resulting 3D data sets (providing high-resolution images of the subsurface) are typically interpreted manually. This is generally an extremely time-consuming analysis step. Therefore, representative 2D sections highlighting distinctive reflection structures are often selected from the 3D data set. Regions showing similar structures are then grouped into so-called radar facies. The results obtained from 2D sections are considered representative of the entire investigated area. Interpretations conducted in this manner are often incomplete and highly dependent on the expertise of the interpreters, making them generally non-reproducible. A promising alternative or complement to manual interpretation is the use of GPR attributes. Instead of using the recorded data directly, derived quantities characterizing distinctive reflection structures in 3D are applied for interpretation. Using various field and synthetic data sets, this thesis investigates which attributes are particularly suitable for this purpose. Additionally, the study demonstrates how selected attributes can be utilized through specific processing and classification methods to create 3D facies models. The ability to generate attribute-based 3D GPR facies models allows for partially automated and more efficient interpretations in the future. Furthermore, the results obtained in this manner describe the subsurface in a reproducible and more comprehensive manner than what has typically been achievable through manual interpretation methods. KW - ground-penetrating radar KW - sedimentary environments KW - 3D KW - applied geophysics KW - near-surface geophysics KW - Georadar KW - sedimentäre Systeme KW - angewandte Geophysik KW - oberflächennahe Geophysik KW - Attribute KW - attributes KW - geophysics KW - Geophysik Y1 - 2024 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-639488 ER -