TY - THES A1 - Gassner, Alexandra Carina T1 - The character of the core-mantle boundary : a systematic study using PcP N2 - Assuming that liquid iron alloy from the outer core interacts with the solid silicate-rich lower mantle the influence on the core-mantle reflected phase PcP is studied. If the core-mantle boundary is not a sharp discontinuity, this becomes apparent in the waveform and amplitude of PcP. Iron-silicate mixing would lead to regions of partial melting with higher density which in turn reduces the velocity of seismic waves. On the basis of the calculation and interpretation of short-period synthetic seismograms, using the reflectivity and Gauss Beam method, a model space is evaluated for these ultra-low velocity zones (ULVZs). The aim of this thesis is to analyse the behaviour of PcP between 10° and 40° source distance for such models using different velocity and density configurations. Furthermore, the resolution limits of seismic data are discussed. The influence of the assumed layer thickness, dominant source frequency and ULVZ topography are analysed. The Gräfenberg and NORSAR arrays are then used to investigate PcP from deep earthquakes and nuclear explosions. The seismic resolution of an ULVZ is limited both for velocity and density contrasts and layer thicknesses. Even a very thin global core-mantle transition zone (CMTZ), rather than a discrete boundary and also with strong impedance contrasts, seems possible: If no precursor is observable but the PcP_model /PcP_smooth amplitude reduction amounts to more than 10%, a very thin ULVZ of 5 km with a first-order discontinuity may exist. Otherwise, if amplitude reductions of less than 10% are obtained, this could indicate either a moderate, thin ULVZ or a gradient mantle-side CMTZ. Synthetic computations reveal notable amplitude variations as function of the distance and the impedance contrasts. Thereby a primary density effect in the very steep-angle range and a pronounced velocity dependency in the wide-angle region can be predicted. In view of the modelled findings, there is evidence for a 10 to 13.5 km thick ULVZ 600 km south-eastern of Moscow with a NW-SE extension of about 450 km. Here a single specific assumption about the velocity and density anomaly is not possible. This is in agreement with the synthetic results in which several models create similar amplitude-waveform characteristics. For example, a ULVZ model with contrasts of -5% VP , -15% VS and +5% density explain the measured PcP amplitudes. Moreover, below SW Finland and NNW of the Caspian Sea a CMB topography can be assumed. The amplitude measurements indicate a wavelength of 200 km and a height of 1 km topography, previously also shown in the study by Kampfmann and Müller (1989). Better constraints might be provided by a joined analysis of seismological data, mineralogical experiments and geodynamic modelling. N2 - Unter der Annahme, dass flüssiges Eisen aus dem äußeren Erdkern mit dem festen, silikat-reichen Unteren Mantel reagiert, wird eine Einflussnahme auf die Kern-Mantel Reflexionsphase PcP erwartet. Ist die Kern-Mantel Grenze aufgeweicht, und nicht wie bislang angenommen ein diskreter Übergang, so zeichnet sich dies in der Wellenform und Amplitude von PcP ab. Die Interaktion mit Eisen führt zu teilweise aufgeschmolzenen Bereichen höherer Dichte, welche die seismischen Wellengeschwindigkeiten herabsetzen. Basierend auf den Berechnungen von kurzperiodischen synthetischen Seismogrammen, mittels der Reflektivitäts- und Gauss Beam Methode, soll ein möglicher Modellraum dieser Niedriggeschwindigkeitszonen ermittelt werden. Das Ziel dieser Arbeit ist es das Verhalten von PcP im Distanzbereich von 10° bis 40° unter dem Einfluss dieser Modelle mit diversen Geschwindigkeits- und Dichtekontrasten zu untersuchen. Ferner wird das Auflösungsvermögen hinsichtlich seismischer Daten diskutiert. Entscheidende Parameter wie Anomaliedicke, Quellfrequenz und Topographie werden hierbei analysiert. Tiefe Erdbeben und Kernexplosionen, die sich im entsprechenden Entfernungsbereich zum Gräfenberg und NORSAR Array befinden, werden anschließend im Hinblick auf PcP ausgewertet. Das seismische Auflösungsvermögen von Niedriggeschwindigkeitszonen ist stark begrenzt sowohl in Bezug auf Geschwindigkeits- und Dichtekontraste als auch hinsichtlich der Mächtigkeit. Es besteht sogar die Möglichkeit einer dünnen, globalen Kern-Mantel Übergangszone, selbst mit großen Impedanzkontrasten, ohne dass dies mit seismologischen Methoden detektiert werden könnte: Wird kein precursor zu PcP beobachtet aber das PcPmodel /PcPsmooth Amplitudenverhältnis zeigt gleichzeitig eine Reduktion von mehr als 10%, dann könnte eine sehr dünne Niedriggeschwindigkeitszone von ca. 5 km Mächtigkeit und einer Diskontinuität erster Ordnung vorliegen. Andererseits, ist PcP um weniger als 10% reduziert, könnte dies entweder auf eine dünne, moderate Niedriggeschwindigkeitszone oder einen graduellen Kern-Mantel Übergang hindeuten. Die synthetischen Berechnungen ergeben starke Amplitudenvariationen als Funktion der Distanz, welche auf den Impedanzkontrast zurückzuführen sind. Dabei ergibt sich ein primärer Dichteeffekt im extremen Steilwinkelbereich und ein maßgeblicher Geschwindigkeitseinfluss im Weitwinkelbereich. Im Hinblick auf die modellierten Resultate lässt sich eine 10 - 13.5 km mächtige Niedriggeschwindigkeitszone 600 km südöstlich von Moskau mit einer NW-SE Ausdehnung von mindestens 450 km folgern, wobei eine exakte Aussage über Geschwindigkeiten und Dichte nicht möglich ist. Dies ist im Konsens mit den synthetischen Berechnungen, wonach viele unterschiedliche Modelle ähnliche Amplituden- und Wellenformcharakteristiken erzeugen. Zum Beispiel erklärt ein Modell mit Kontrasten von -5% VP , -15% VS and +5% Dichte die gemessenen PcP Amplituden. Darüber hinaus können unterhalb des südwestlichen Finnlands und nord-nordwestlich des Kaspischen Meeres Undulationen an der Kern-Mantel Grenze selbst vermutet werden. Unter Berücksichtigung früherer Studien, z. B. von Kampfmann and Müller (1989), deuten die Messergebnisse auf eine laterale Topographie von 200 km und eine Höhe von 1 km hin. Eine Eingrenzung der potentiellen Anomaliemodelle kann nur durch eine gemeinsame Auswertung mit mineralogischen Experimenten und geodynamischen Modellierungen erfolgen. KW - Kern-Mantel Grenze KW - Seismologie KW - Ultra-Niedriggeschwindigkeitszonen KW - Steilwinkel-Analyse von PcP KW - Tiefbeben und Kernexplosionen KW - core-mantle boundary KW - seismology KW - ultra-low velocity zones KW - steep-angle analysis of PcP KW - deep earthquakes and nuclear explosions Y1 - 2012 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-63590 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 - GEN A1 - Krüger, Frank A1 - Ohrnberger, Matthias A1 - Rößler, Dirk T1 - Rupture imaging of large earthquakes with a poststack isochrone migration method N2 - Rapid and robust characterization of large earthquakes in terms of their spatial extent and temporal duration is of high importance for disaster mitigation and early warning applications. Backtracking of seismic P-waves was successfully used by several authors to image the rupture process of the great Sumatra earthquake (26.12.2004) using short period and broadband arrays. We follow here an approach of Walker et al. to backtrack and stack broadband waveforms from global network stations using traveltimes for a global Earth model to obtain the overall spatio-temporal development of the energy radiation of large earthquakes in a quick and robust way. We present results for selected events with well studied source processes (Kokoxili 14.11.2001, Tokachi-Oki 25.09.2003, Nias 28.03.2005). Further, we apply the technique in a semi-real time fashion to broadband data of earthquakes with a broadband magnitude >= 7 (roughly corresponding to Mw 6.5). Processing is based on first automatic detection messages from the GEOFON extended virtual network (GEVN). KW - Seismologie KW - Erdbeben KW - Array Seismologie KW - Migration KW - Seismology KW - Earthquake KW - Array Seismology KW - Migration Y1 - 2008 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-18395 ER - TY - THES A1 - Kulikova, Galina T1 - Source parameters of the major historical earthquakes in the Tien-Shan region from the late 19th to the early 20th century T1 - Herdparameter von starken, historischen Erdbeben in der Tien-Shan Region zur Jahrhundertwende (19./20. Jahrhundert) N2 - The Tien-Shan and the neighboring Pamir region are two of the largest mountain belts in the world. Their deformation is dominated by intermontane basins bounded by active thrust and reverse faulting. The Tien-Shan mountain belt is characterized by a very high rate of seismicity along its margins as well as within the Tien-Shan interior. The study area of the here presented thesis, the western part of the Tien-Shan region, is currently seismically active with small and moderate sized earthquakes. However, at the end of the 19th beginning of the 20th century, this region was struck by a remarkable series of large magnitude (M>7) earthquakes, two of them reached magnitude 8. Those large earthquakes occurred prior to the installation of the global digital seismic network and therefore were recorded only by analog seismic instruments. The processing of the analog data brings several difficulties, for example, not always the true parameters of the recording system are known. Another complicated task is the digitization of those records - a very time-consuming and delicate part. Therefore a special set of techniques is developed and modern methods are adapted for the digitized instrumental data analysis. The main goal of the presented thesis is to evaluate the impact of large magnitude M≥7.0 earthquakes, which occurred at the turn of 19th to 20th century in the Tien-Shan region, on the overall regional tectonics. A further objective is to investigate the accuracy of previously estimated source parameters for those earthquakes, which were mainly based on macroseismic observations, and re-estimate them based on the instrumental data. An additional aim of this study is to develop the tools and methods for faster and more productive usage of analog seismic data in modern seismology. In this thesis, the ten strongest and most interesting historical earthquakes in Tien-Shan region are analyzed. The methods and tool for digitizing and processing the analog seismic data are presented. The source parameters of the two major M≥8.0 earthquakes in the Northern Tien-Shan are re-estimated in individual case studies. Those studies are published as peer-reviewed scientific articles in reputed journals. Additionally, the Sarez-Pamir earthquake and its connection with one of the largest landslides in the world, Usoy landslide, is investigated by seismic modeling. These results are also published as a research paper. With the developed techniques, the source parameters of seven more major earthquakes in the region are determined and their impact on the regional tectonics was investigated. The large magnitudes of those earthquakes are confirmed by instrumental data. The focal mechanism of these earthquakes were determined providing evidence for responsible faults or fault systems. N2 - Der Tien-Shan und die angrenzende Pamir Region sind zwei der größten Gebirgszüge der Welt. Deformation findet hier hauptsächlich an aktiven Auf- und Abschiebungszonen statt, welche intermontane Becken umschließen. Der Tien-Shan Gebirgszug ist sowohl an den Störungszonen als auch innerhalb der Becken durch eine hohe Seismizitätsrate charakterisiert. Das Untersuchungsgebiet der hier präsentierten Dissertation, der westliche Bereich der Tien-Shan Region, ist in den letzten Jahrzehnten durch das Auftreten kleiner und mittlerer Erdbeben gekennzeichnet. Jedoch wurde diese Region an der Wende zum 20. Jahrhundert von einer Reihe außergewöhnlich starker Erdbeben (M>7) heimgesucht. Zwei von ihnen erreichten sogar die Magnitude 8. Diese starken Erdbeben ereigneten sich vor der Installation eines globalen, digitalen seismischen Netzwerks und wurden daher nur von analogen seismischen Instrumenten aufgezeichnet. Die Bearbeitung von analogen Daten birgt mehrere Schwierigkeiten, z.B. sind die wahren Werte der Instrumentencharakterisitik nicht immer bekannt. Ein weiterer komplizierter Teil ist die Digitalisierung dieser Aufzeichnungen, die sehr zeitaufwändig und diffizil ist. Um diesen und weiteren Schwierigkeiten zu begegnen, wurden in der vorliegenden Arbeit spezielle Techniken zur Digitalisierung analoger, seismischer Daten und moderne Methoden der Datenanalyse speziell für digitalisierte, analoge Daten adaptiert. Das Hauptziel der hier präsentierten Dissertation ist die Auswertung der Auswirkungen von starken Erdbeben (Magnitude 7.0) auf die regionale Tektonik, welche zur Jahrhundertwende zum 20. Jahrhundert in der Tian Shan Region stattgefunden haben. Eine weitere Zielsetzung ist die Überprüfung der Genauigkeit von früher bestimmten Herdparametern dieser Erdbeben, welche hauptsächlich auf makroseismischen Untersuchungen beruhen, und deren erneute Bestimmung mithilfe instrumenteller Daten. Außerdem sollen in dieser Arbeit die notwendigen Werkzeuge und Methoden für die schnellere und produktivere Nutzung von analogen, seismischen Daten für zukünftige Studien in der modernen Seismologie entwickelt werden. In dieser Arbeit werden die zehn stärksten und interessantesten historischen Erdbeben der Tien-Shan Region analysiert. Es werden Methoden und Werkzeuge für die Digitalisierung und Bearbeitung von analogen seismischen Daten vorgestellt. Die Herdparameter der zwei bedeutendsten Erdbeben mit M~8 im nördlichen Tien-Shan werden erneut bestimmt und in separaten, detaillierten Fallstudien behandelt. Diese Studien sind als wissenschaftlich begutachtete Artikel in renommierten Fachzeitschriften publiziert. Zusätzlich wurde das Sarez-Pamir Erdbeben und seine Verbindung mit einem der größten Erdrutsche der Welt, dem Usoy Erdrutsch, basierend auf seismischer Modellierung untersucht. Diese Ergebnisse sind ebenfalls in einem wissenschaftlichen Aufsatz veröffentlicht. Mit den hier entwickelten Methoden wurden die Herdparameter der sieben stärksten Erdbeben in der Region bestimmt und ihre Auswirkung auf die regionale Tektonik untersucht. Die Magnituden dieser Erdbeben wurden mit Hilfe instrumental aufgezeichneter Daten bestätigt. Die Herdflächenlösungen dieser Erdbeben wurden bestimmt und geben Hinweise auf die möglichen verantwortlichen Störungen. KW - Seismology KW - historical earthquakes KW - Seismologie KW - historische Erdbeben Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-88370 ER - TY - THES A1 - Köhler, Andreas T1 - Recognition and investigation of temporal patterns in seismic wavefields using unsupervised learning techniques T1 - Unüberwachte Erkennung und Untersuchung von zeitlichen Mustern in seismischen Wellenfeldern N2 - Modern acquisition of seismic data on receiver networks worldwide produces an increasing amount of continuous wavefield recordings. Hence, in addition to manual data inspection, seismogram interpretation requires new processing utilities for event detection, signal classification and data visualization. Various machine learning algorithms, which can be adapted to seismological problems, have been suggested in the field of pattern recognition. This can be done either by means of supervised learning using manually defined training data or by unsupervised clustering and visualization. The latter allows the recognition of wavefield patterns, such as short-term transients and long-term variations, with a minimum of domain knowledge. Besides classical earthquake seismology, investigations of temporal patterns in seismic data also concern novel approaches such as noise cross-correlation or ambient seismic vibration analysis in general, which have moved into focus within the last decade. In order to find records suitable for the respective approach or simply for quality control, unsupervised preprocessing becomes important and valuable for large data sets. Machine learning techniques require the parametrization of the data using feature vectors. Applied to seismic recordings, wavefield properties have to be computed from the raw seismograms. For an unsupervised approach, all potential wavefield features have to be considered to reduce subjectivity to a minimum. Furthermore, automatic dimensionality reduction, i.e. feature selection, is required in order to decrease computational cost, enhance interpretability and improve discriminative power. This study presents an unsupervised feature selection and learning approach for the discovery, imaging and interpretation of significant temporal patterns in seismic single-station or network recordings. In particular, techniques permitting an intuitive, quickly interpretable and concise overview of available records are suggested. For this purpose, the data is parametrized by real-valued feature vectors for short time windows using standard seismic analysis tools as feature generation methods, such as frequency-wavenumber, polarization, and spectral analysis. The choice of the time window length is dependent on the expected durations of patterns to be recognized or discriminated. We use Self-Organizing Maps (SOMs) for a data-driven feature selection, visualization and clustering procedure, which is particularly suitable for high-dimensional data sets. Using synthetics composed of Rayleigh and Love waves and three different types of real-world data sets, we show the robustness and reliability of our unsupervised learning approach with respect to the effect of algorithm parameters and data set properties. Furthermore, we approve the capability of the clustering and imaging techniques. For all data, we find improved discriminative power of our feature selection procedure compared to feature subsets manually selected from individual wavefield parametrization methods. In particular, enhanced performance is observed compared to the most favorable individual feature generation method, which is found to be the frequency spectrum. The method is applied to regional earthquake records at the European Broadband Network with the aim to define suitable features for earthquake detection and seismic phase classification. For the latter, we find that a combination of spectral and polarization features favor S wave detection at a single receiver. However, SOM-based visualization of phase discrimination shows that clustering applied to the records of two stations only allows onset or P wave detection, respectively. In order to improve the discrimination of S waves on receiver networks, we recommend to consider additionally the temporal context of feature vectors. The application to continuous recordings of seismicity close to an active volcano (Mount Merapi, Java, Indonesia) shows that two typical volcano-seismic events (VTB and Guguran) can be detected and distinguished by clustering. In contrast, so-called MP events cannot be discriminated. Comparable results are obtained for selected features and recognition rates regarding a previously implemented supervised classification system. Finally, we test the reliability of wavefield clustering to improve common ambient vibration analysis methods such as estimation of dispersion curves and horizontal to vertical spectral ratios. It is found, that in general, the identified short- and long-term patterns have no significant impact on those estimates. However, for individual sites, effects of local sources can be identified. Leaving out the corresponding clusters, yields reduced uncertainties or allows for improving estimation of dispersion curves. N2 - Die Anzahl der weltweit kontinuierlich aufzeichnenden seismischen Messstationen ist in den vergangenen Jahren immer weiter angestiegen. Aus diesem Grund steht eine große Menge von seismischen Datensätzen zu Forschungszwecken zur Verfügung. Insbesondere betrifft dies passive Verfahren zur geologischen Strukturerkundung entweder mittels transienter Ereignisse wie Erdbeben oder unter der Verwendung der permanent vorhandenen natürlichen seismischen Bodenunruhe. Die Bearbeitung dieser Daten erfordert neben der klassischen manuellen Seismogrammanalyse verstärkt auch den Einsatz automatischer Detektionssysteme. Mit Hilfe von überwachten Lernverfahren, d.h. unter Verwendung von seismischen Signalen deren Auftreten bekannt ist, ist es möglich, unbekannte Muster zu klassifizieren. Im Gegensatz dazu hatte die vorliegende Arbeit zum Ziel, ein allgemeines, unüberwachtes Verfahren zur quantitativen Zerlegung seismischer Wellenfelder zu entwickeln. Dies wird mittels einer automatischen Clusterung von Seismogrammzeitfenstern bzw. über die Visualisierung von zeitlichen Mustern auf unterschiedlichen Zeitskalen erreicht. Als unüberwachtes Lernverfahren, das neben der Clusterung auch eine einfach interpretierbare Visualisierung hoch-dimensionaler Datensätze über eine zweidimensionale Darstellung ermöglicht, wurde der Self-organizing-map Algorithmus (SOM) gewählt. Für automatische Lernverfahren ist die Parametrisierung der Seismogramme mittels Merkmalsvektoren erforderlich. Im vorliegenden Fall wurden möglichst viele potentielle Wellenfeldmerkmale unter Verwendung von verschiedenen seismischen Einzel- und Mehrstationsanalyseverfahren für aufeinanderfolgende kurze Zeitfenster berechnet. Um eine datenadaptive und effiziente Parametrisierung zu erreichen, wurde darüberhinaus ein quantitatives Auswahlverfahren für geeignete Merkmale entwickelt, das über einen mehrstufigen Filter bestehend aus einem Signifikanztest und einer SOM-basierenden Korrelationsanalyse redundante und irrelevante Eigenschaften aussortiert. Mit den neu implementierten Techniken wurden verschiedene Arten von seismischen Datensätzen unter Berücksichtigung verschiedener seismologischer Fragestellungen bearbeitet. Die Algorithmen und deren Parameter wurden zunächst intensiv und quantitativ mit Hilfe synthetischer Daten getestet und optimiert. Anschließend wurden reale Aufzeichnungen regionaler Erdbeben und vulkanischer Seismizität verwendet. Im ersten Fall konnten geeignete Merkmale zur Detektion und Klassifizierung von Erdbebenwellenphasen gefunden und die Diskriminierung dieser Signale mit Hilfe der SOM-Darstellung untersucht werden. Unter Verwendung des zweiten Datensatzes wurden Cluster typischer vulkano-seismischer Signale am Vulkan Mount Merapi (Java, Indonesien) detektiert, die sich zur Vorhersage von Eruptionen eignen. Beide Anwendungen haben gezeigt, dass, verglichen mit einzelnen Methoden, automatisch gefundene Kombinationen von Merkmalen verschiedener Parametrisierungsverfahren deutlich bessere Klassifizierungsraten zur Folge haben. Zudem können die Erkenntnisse über die Clusterung von seismischen Signalen dazu verwendet werden, verbesserte automatische Klassifizierungssysteme zu entwickeln. Abschließend wurden Aufzeichnungen der natürlichen seismischen Bodenunruhe bearbeitet. Insbesondere konnte der Einfluss kurzzeitiger und längerfristiger Variationen im Wellenfeld auf Methoden zur passiven Strukturerkundung untersucht werden. Es hat sich gezeigt, dass in einzelnen Fällen tageszeitabhängige Muster und lokale seismische Quellen die Ergebnisse negativ beeinflussen können. Die Wellenfeldzerlegung mittels Clusterung hat es erlaubt, diese Signale zu identifizieren und somit von der Analyse auszuschließen. KW - Seismologie KW - Mustererkennung KW - Unüberwachtes Lernen KW - Seismology KW - Pattern Recognition KW - Unsupervised Learning Y1 - 2009 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-29702 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 - GEN A1 - Lipke, Katrin A1 - Krüger, Frank A1 - Rößler, Dirk T1 - Subduction zone structure along Sumatra from receiver functions N2 - Receiver functions are a good tool to investigate the seismotectonic structure beneath the a seismic station. In this study we apply the method to stations situated on or near Sumatra to find constraints on a more detailed velocity model which should improve earthquake localisation. We estimate shallow Moho-depths (~ 21 km) close to the trench and depths of ~30 km at greater distances. First evidences for the dip direction of the slab of ~60° are provided. Receiver functions were calculated for 20 stations for altogether 110 earthquakes in the distance range between 30° and 95° from the receiver. However the number of receiver functions per station is strongly variable as it depends on the installation date, the signal-to-noise-ratio of the station and the reliability of the acquisition. N2 - Receiver Funkttion stellen eine gut Methode zur Untersuchung von Seismotektonischen Strukturen unterhalb einer seismischen Station dar. In dieser Arbeit wenden wir die Methode auf Station auf oder nahe Sumatra an um Hinweise für ein detaillierteres Geschwindigkeitsmodell zu erhalten, welches die Lokalisierung von Erdbeben verbessern sollte. Wir ermitteln flache Moho-Tiefen (~21 km) in der Nähe des Trenchs und Tiefen um die 30 km in größeren Distanzen. Erste Hinweise für eine Einfallsrichtung des Slabs von ~60° konnten gefunden werden. Receiver Funktionen wurden für 20 Stationen für insgesamt 110 Erdbeben im Distanzbereich zwischen 30° und 95° berechnet. allerdings ist die Anzahl von Receiver Funktionen pro Station sehr variabel, da sie vom Installationszeitpunkt, dem Signal-Rausch-Verhältnis und der Zuverlässigkeit der Datenaufnahme an der Station abhängt. KW - Receiver Funktionen KW - Sumatra KW - Seismologie KW - Subduktion KW - Receiver Functions KW - Sumatra KW - Seismology KW - Subduction Y1 - 2008 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-18260 ER - TY - JOUR A1 - Lipke, Katrin A1 - Zitzmann, Max A1 - Amberger, Manuel A1 - Ehlert, Carsten A1 - Rößler, Dirk A1 - Krüger, Frank A1 - Ohrnberger, Matthias T1 - Traveltime residuals at regional and teleseismic distances for SE-Asia N2 - Traveltime residuals for worldwide seismic stations are calculated. We use P and S waves from earthquakes in SE-Asia at teleseismic and regional distances. The obtained station residuals help to enhance earthquake localisation. Furthermore we calculated regional source dependent station residuals. They show a systematic dependence of the locality of the source. These source dependent residuals reflect heterogenities along the path and can be used for a refinement of earthquake localisation. N2 - Laufzeitresiduen für weltweite seismische Stationen werden berechnet. Wir nutzen P - und S-Wellen von Erdbeben in Südostasien in teleseismischen und regionalen Distanzen. Die so erhaltenen Stationsresiduen helfen, die Lokaliesierung von Erdbeben zu verbessern. Außerdem berechnen wir regional quellabhängige Stationsresiduen. Diese zeigen eine systematische Abhänbgigkeit vom Ort der Quelle. Sie spiegeln Heterogenitäten entlang des Strahlweges wieder und können für eine Verfeinerung der Ersbebenlokaliesierung genutzt werden. KW - Seismologie KW - Südostasien KW - Laufzeitresiduen KW - GITEWS KW - Wellengeschwindigkeiten KW - seismology KW - Southeast Asia KW - traveltime KW - residuals KW - GITEWS Y1 - 2007 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-14117 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 - Metzger, Sabrina T1 - Neotectonic deformation over space and time as observed by space-based geodesy T1 - Über die Vermessung neotektonischer Deformation in Raum und mit Hilfe von satellitengestützter Geodäsie N2 - Alfred Wegeners ideas on continental drift were doubted for several decades until the discovery of polarization changes at the Atlantic seafloor and the seismic catalogs imaging oceanic subduction underneath the continental crust (Wadati-Benioff Zone). It took another 20 years until plate motion could be directly observed and quantified by using space geodesy. Since then, it is unthinkable to do neotectonic research without the use of satellite-based methods. Thanks to a tremendeous increase of instrumental observations in space and time over the last decades we significantly increased our knowledge on the complexity of the seismic cycle, that is, the interplay of tectonic stress build up and release. Our classical assumption, earthquakes were the only significant phenomena of strain release previously accumulated in a linear fashion, is outdated. We now know that this concept is actually decorated with a wide range of slow and fast processes such as triggered slip, afterslip, post-seismic and visco-elastic relaxation of the lower crust, dynamic pore-pressure changes in the elastic crust, aseismic creep, slow slip events and seismic swarms. On the basis of eleven peer-reviewed papers studies I here present the diversity of crustal deformation processes. Based on time-series analyses of radar imagery and satellited-based positioning data I quantify tectonic surface deformation and use numerical and analytical models and independent geologic and seismologic data to better understand the underlying crustal processes. The main part of my work focuses on the deformation observed in the Pamir, the Hindu Kush and the Tian Shan that together build the highly active continental collision zone between Northwest-India and Eurasia. Centered around the Sarez earthquake that ruptured the center of the Pamir in 2015 I present diverse examples of crustal deformation phenomena. Driver of the deformation is the Indian indenter, bulldozing into the Pamir, compressing the orogen that then collapses westward into the Tajik depression. A second natural observatory of mine to study tectonic deformation is the oceanic subduction zone in Chile that repeatedly hosts large earthquakes of magnitude 8 and more. These are best to study post-seismic relaxation processes and coupling of large earthquake. My findings nicely illustrate how complex fashion and how much the different deformation phenomena are coupled in space and time. My publications contribute to the awareness that the classical concept of the seismic cycle needs to be revised, which, in turn, has a large influence in the classical, probabilistic seismic hazard assessment that primarily relies on statistically solid recurrence times. N2 - Alfred Wegeners Thesen des Kontinentaldrifts fanden erst in den 1960er und 1970er Jahren Akzeptanz, als die krustalen Polarisationswechsel auf dem atlantischen Meeresboden entdeckt wurden und Erdbebenkataloge das Abtauchen von ozeanischer Kruste unter kontinentale Kruste abbildeten (Wadati-Benioff-Zone). Es dauerte jedoch weitere 20 Jahre, bis die Geodäsie erstmals Plattenbewegung sicht- und quantifizierbar machte. Seit dann sind satellitengestützte Messmethoden aus der neotektonischen Forschung nicht mehr wegzudenken. Dank einer stetig (zeitlich und räumlich) wachsenden Anzahl instrumenteller Beobachtungsdaten wird unser Verständnis des Erdbebenzyklus—des Wechselspiels zwischen tektonischem Spannungsauf- und -abbau—immer komplexer. Das klassische Konzept, nur Erdbeben setzten die zuvor linear aufgebaute Spannungsenergie instantan frei, wird heutzutage durch eine Vielzahl von zusätzlichen schnelleren und langsameren Prozessen ergänzt. Beispiele dafür sind getriggerte Versätze (triggered slip), Nachbeben (afterslip), postseismische und visko-elastische Relaxation der tieferen Kruste, dynamische, elastische Veränderungen des Gesteins-Porendrucks, aseismisches Kriechen sowie Spannungsabbau durch kleine Erdbebenschwärme. Anhand von elf begutachteten und bereits veröffentlichten Arbeiten präsentiere ich in meiner Habilitationsschrift die Diversität krustaler Deformationsprozesse. Ich analysiere Zeitreihen von Radar-Satellitenaufnahmen und satellitengestützten Positionierungssystemen um die tektonische Oberflächenbewegung zu quantifizieren. Der Vergleich von kinematischen Beobachtungen mit geologischen und seismischen Indizien sowie die Simulation ebenjener durch rechnergestützte Modelle ermöglichen mir, die verursachenden krustalen Prozesse besser verstehen. Der Hauptteil meiner Arbeiten beschreibt rezente, krustale Bewegungen im Pamir, Hindu Kush und Tien Shan, welche zusammen das westliche Ende der kontinentalen Kollisionszone zwischen dem indischen und eurasischen Kontinent bilden. Rund um ein starkes Erdbeben, welches 2015 den Zentralpamir erschüttert hat, zeige ich vielseitige Beispiele von hochaktiver krustaler Deformation. Verursacht werden diese Bewegungen durch den nordwestindischen Kontinentalsporn, welcher (fast) ungebremst in den Pamir hineinrammt, ihn auftürmt, zusammenquetscht, und ihn gravitationsbedingt gegen Westen ins tadschikische Becken kollabieren lässt. Der zweite thematische Schwerpunkt liegt auf Prozessen, welche durch Megathrust-Erdbeben, also Beben mit einer Magnitude>8, hervorgerufen werden. Diese Anwendungen fokussieren sich auf die ozeanischen Subduktionszone von Chile und zeigen die Wichtigkeit vertikaler Hebungsdaten um, beispielsweise, den Einfluss tektonischer Prozesse auf den Gesteins-Porendruck zu verstehen. Zusammenfassend veranschaulichen und bestätigen meine Arbeiten, wie stark und komplex die oben beschriebenen Prozesse räumlich und zeitlich korrelieren, und dass das klassische Konzept des Erdbebenzyklus überholt ist. Letztere Einsicht hat grossen Einfluss auf probabilistische seismische Gefährdungsanalysen, welche grundsätzlich statistische Vorhersagbarkeit annehmen. KW - radar satellite interferometry KW - tectonics KW - geodesy KW - seismology KW - earthquakes KW - InSAR KW - InSAR KW - Erdbeben KW - Geodäsie KW - Radar-Satelliteninterferometrie KW - Seismologie KW - Tektonik Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-599225 ER - TY - THES A1 - Petersen, Gesa Maria T1 - Source array and receiver array analysis of Vogtland/ West Bohemia earthquake clusters T1 - Untersuchung von Erdbebenclustern im Vogtland/ Westböhmen durch Quell- und Empfängerarrays N2 - Die Region Vogtland/ West Böhmen im Grenzgebiet zwischen Deutschland und Tschechien ist bekannt für ihre geologische Aktivität. Holozäner Vulkanismus, Gasaustritte an Mofetten und Quellen und wiederkehrende Erdbebenschwärme sind Ausdruck geodynamischer Prozesse im Untergrund. Während des Erdbebenschwarms 2008/2009 in Nový Kostel installierte die Universität Potsdam ein temporäres Array in Rohrbach, in einer Epizentraldistanz von etwa 10 km und mit einer Aperatur von etwa 0.75 km. 22 Erdbeben wurden für Quellarray- beam forming ausgewählt. Quellarrays sind örtliche Cluster von Erdbeben, die von einer Empfängerstation aufgezeichnet werden. Wegen der Reziprozität der Green’schen Funktionen können diese in ähnlicher Weise genutzt werden wie Empfängerarrays, bei denen mehrere Stationen ein einzelnes Beben aufzeichnen. Die Kreuzkorrelationskoeffizienten aller Beben des Quellarrays, aufgezeichnet an einer einzelnen Station, sind in der Regel höher als für einzelne Ereignisse, die an allen Stationen des Empfängerarrays aufgezeichnet wurden. Dies deutet hinsichtlich der aufgelösten Frequenzen auf einen heterogenen Untergrund unter den Array-Stationen und ein vergleichsweise homogenes Quellarray-Volumen hin. Beam forming wurde mit den horizontalen und vertikalen Spuren aller Quellarray-Beben, aufgezeichnet auf allen 11 Stationen des Empfängerarrays, durchgeführt. Die Ergebnisse wurden im Hinblick auf Konversionen und reflektierte Phasen analysiert. Während die theoretische Richtung der direkten P-Welle im Falle der Quellarray-Aufzeichnungen gut übereinstimmt, wird eine Empfängerarray-Missweisung von 15° bis 25° beobachtet. Eine PS Phase, die der direkten P Phase folgt und eine mögliche SP Phase, die kurz vor der direkten S-Phase ankommt, wurden auf den summierten Spuren mehrerer Stationen interpretiert. Aus der Betrachtung der Laufzeiten resultiert eine Konversionstiefe von 0.6-0.9 km Tiefe. Ein zweites Quellarray, bestehend aus 12 tieferen Beben wurde zusätzlich analysiert, um eine nach ca. 0.85 s ausschließlich auf den Aufzeichnungen tieferer Beben auftretende Strukturphase zu deuten. Zusätzlich zum beam forming wurden zwei Lokalisierungsmethoden von Reflexionen und Konversionen für einfach reflektierte/konvertierte Phasen entwickelt und zur Auswertung verwendet. Während die erste, analytische Methode eine homogene Geschwindigkeit entlang des Laufwegs annimmt, wird in der zweiten Methode eine 3-D-Rastersuche ausgeführt, in der ein 1-D-Geschwindigkeitsmodell verwendet wird. Auf Grund der eindeutigen beam forming Ergebnisse und der hohen Ähnlichkeit der Wellenformen der Erdbeben, die für das Quellarray genutzt wurden, bieten Quellarrays bestehend aus Mikrobeben aus dem untersuchten Gebiet gute Möglichkeiten zur Untersuchung von Krustenstrukturen. N2 - The Vogtland, located at the border region between the Czech Republic and Germany, is known for Holocene volcanism, gas and fluid emissions as well as for reoccurring earthquake swarms, pointing towards a high geodynamic activity. During the earthquake swarm in 2008/2009, a temporary array was installed close to Rohrbach (Germany), at an epicentral distance of about 10 km from the Nový Kostel focal zone (aperture ~0.75 km). 22 events of the recorded swarm were selected to set up a source array. Source arrays are spatially clustered earthquakes, which can be used in a similar manner as receiver array recordings of single events (Green’s functions reciprocity). The application of array seismology techniques like beam forming requires similar waveforms and precisely known origin times and locations. The resemblance of waveforms was assured by visual selection of events and quantified with the calculation of cross-correlation coefficients. We observed that the different events recorded at a single station generally show greater resemblances than the recordings of one event at all stations of the receiver array. This indicates a heterogeneous subsurface beneath the receiver array and a comparably homogeneous source array volume with respect to the frequency-dependent resolution of both arrays. Beam forming was applied on the Z, N and E component recordings of the source array events at 11 stations, and the results were analysed with respect to converted or reflected crustal phases. While the theoretical back azimuth of the direct phases match the beam forming results in case of the source array analysis, in case of receiver array beam forming derivations of 15°-25° are observed. PS phases, closely following the direct P phase and presumably SP phases, arriving shortly before the direct S phase can be observed on several stations. Based on the time differences to the direct P and S phases we inferred a conversion depth of about 0.6-0.9 km. A second deeper source array was set up in order to interpret a structural phase arriving 0.85 s after the direct P phase on records of deeper events only. Additionally to the source array beam forming method an analytical method with a fixed medium velocity and a grid search method, both for determining conversion/ reflection locations of phases traveling off the direct line between source and receiver array, were developed and applied to other observed phases. In conclusion, we think that the distinct beam forming results along with the striking waveform resemblance reveal the opportunities of using source arrays consisting of small swarm events for the analysis of crustal structures. KW - seismology KW - Vogtland KW - West Bohemia KW - source array KW - beam forming KW - earthquake swarms KW - Seismologie KW - Vogtland KW - Westböhmen KW - Quellarray KW - Beamforming KW - Erdbebenschwärme KW - Schwarmbeben Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-406671 ER - TY - THES A1 - Rößler, Dirk T1 - Retrieval of earthquake source parameters in inhomogeneous anisotropic mediawith application to swarm events in West Bohemia in 2000 T1 - Bestimmung von Erdbebenparametern in inhomogenen anisotropen Medien mit Anwendung auf Schwarmbeben im Vogtland im Jahr 2000 N2 - Earthquakes form by sudden brittle failure of rock mostly as shear ruptures along a rupture plane. Beside this, mechanisms other than pure shearing have been observed for some earthquakes mainly in volcanic areas. Possible explanations include complex rupture geometries and tensile earthquakes. Tensile earthquakes occur by opening or closure of cracks during rupturing. They are likely to be often connected with fluids that cause pressure changes in the pore space of rocks leading to earthquake triggering. Tensile components have been reported for swarm earthquakes in West Bohemia in 2000. The aim and subject of this work is an assessment and the accurate determination of such tensile components for earthquakes in anisotropic media. Currently used standard techniques for the retrieval of earthquake source mechanisms assume isotropic rock properties. By means of moment tensors, equivalent forces acting at the source are used to explain the radiated wavefield. Conversely, seismic anisotropy, i.e. directional dependence of elastic properties, has been observed in the earth's crust and mantle such as in West Bohemia. In comparison to isotropy, anisotropy causes modifications in wave amplitudes and shear-wave splitting. In this work, effects of seismic anisotropy on true or apparent tensile source components of earthquakes are investigated. In addition, earthquake source parameters are determined considering anisotropy. It is shown that moment tensors and radiation patterns due to shear sources in anisotropic media may be similar to those of tensile sources in isotropic media. In contrast, similarities between tensile earthquakes in anisotropic rocks and shear sources in isotropic media may exist. As a consequence, the interpretation of tensile source components is ambiguous. The effects that are due to anisotropy depend on the orientation of the earthquake source and the degree of anisotropy. The moment of an earthquake is also influenced by anisotropy. The orientation of fault planes can be reliably determined even if isotropy instead of anisotropy is assumed and if the spectra of the compressional waves are used. Greater difficulties may arise when the spectra of split shear waves are additionally included. Retrieved moment tensors show systematic artefacts. Observed tensile source components determined for events in West Bohemia in 1997 can only partly be attributed to the effects of moderate anisotropy. Furthermore, moment tensors determined earlier for earthquakes induced at the German Continental Deep Drilling Program (KTB), Bavaria, were reinterpreted under assumptions of anisotropic rock properties near the borehole. The events can be consistently identified as shear sources, although their moment tensors comprise tensile components that are considered to be apparent. These results emphasise the necessity to consider anisotropy to uniquely determine tensile source parameters. Therefore, a new inversion algorithm has been developed, tested, and successfully applied to 112 earthquakes that occurred during the most recent intense swarm episode in West Bohemia in 2000 at the German-Czech border. Their source mechanisms have been retrieved using isotropic and anisotropic velocity models. Determined local magnitudes are in the range between 1.6 and 3.2. Fault-plane solutions are similar to each other and characterised by left-lateral faulting on steeply dipping, roughly North-South oriented rupture planes. Their dip angles decrease above a depth of about 8.4km. Tensile source components indicating positive volume changes are found for more than 60% of the considered earthquakes. Their size depends on source time and location. They are significant at the beginning of the swarm and at depths below 8.4km but they decrease in importance later in the course of the swarm. Determined principle stress axes include P axes striking Northeast and Taxes striking Southeast. They resemble those found earlier in Central Europe. However, depth-dependence in plunge is observed. Plunge angles of the P axes decrease gradually from 50° towards shallow angles with increasing depth. In contrast, the plunge angles of the T axes change rapidly from about 8° above a depth of 8.4km to 21° below this depth. By this thesis, spatial and temporal variations in tensile source components and stress conditions have been reported for the first time for swarm earthquakes in West Bohemia in 2000. They also persist, when anisotropy is assumed and can be explained by intrusion of fluids into the opened cracks during tensile faulting. N2 - Erdbeben entstehen durch plötzlichen Sprödbruch des Gesteins, meist als Scherbruch entlang einer Bruchfläche. Daneben werden für einige Beben v.a. in vulkanischen Gebieten auch Mechanismen beobachtet, die scheinbar vom Modell des Scherbruches abweichen. Ursachen dafür beinhalten komplexe Bruchgeometrien und tensile Erdbeben. Bei tensilen Erdbeben kommt es während des Bruchs zum Öffnen oder Schließen der Bruchfläche und damit zu Volumenänderungen. Erdbeben mit tensilen Anteilen stehen wahrscheinlich oft im Zusammenhang mit Fluiden, welche zur Durckänderung im Porenraum von Gesteinen und damit zum Auslösen des Bebens führen. Sie wurden auch im Vogtland während eines Erdbebenschwarms im Jahr 1997 beobachtet. Die Beurteilung und sichere Bestimmung tensiler Anteile von Erdbeben sind Ziel und Gegenstand dieser Arbeit. Bei Standardverfahren zur Bestimmung von Erdbebenmechanismen werden isotrope Gesteinseigenschaften angenommen. Momententensoren beschreiben dabei Kräfte, die das abgestrahlte Wellenfeld erklären. Allerdings wird seismische Anisotropie, d.h. Richtungsabhängigkeit elastischer Eigenschaften, in der Erdkruste und im Mantel wie z.B. im Vogtland beobachtet. Anisotropie bewirkt im Vergleich zu isotropen Medien Veränderungen der Wellenamplituden und -polariserungen sowie das Aufspalten von Scherwellen. In der vorliegenden Arbeit werden daher der Einfluss seismischer Anisotropie auf wahre oder scheinbar auftretende tensile Quellanteile untersucht und Erdbebenmechanismen unter Berücksichtigung seismischer Anisotropie bestimmt. Es wird gezeigt, dass Momententensoren und Abstrahlmuster von Scherbrüchen in anisotropen Medien denen von tensilen Brüchen in isotropen Medien ähneln können. Umgekehrt treten Ähnlichkeiten tensiler Beben in anisotropen Gesteinen mit Scherbrüchen in isotropen Medien auf. Damit existieren Mehrdeutigkeiten beobachteter tensiler Quellanteile. Die Effekte von Anisotropie hängen von der Orientierung des Bruches und vom Grad der Anisotropie ab. Außerdem beeinflusst Anisotropie das Moment eines Bebens. Herdflächenorientierungen können auch dann verlässlich bestimmt werden, wenn man Isotropie statt Anisotropie annimmt und die Spektren von Kompressionswellen verwendet. Bei Hinzunahme der Spektren von Scherwellen können Uneindeutigkeiten auftreten. Abgeleitete Momententensoren zeigen systematische Artefakte. Beobachtungen tensiler Quellanteile von Beben im Vogtland im Jahr 1997 können nicht allein durch moderate Anisotropie erklärt werden. Weiterhin wurden früher bestimmte Momententensoren induzierter Beben nahe der Kontinentalen Tiefbohrung, Bayern, unter Annahme anisotroper Parameter reinterpretiert. Die Beben werden einheitlich als Scherbrüche charakterisiert, obwohl deren Momententensoren tensile Bestandteile enthalten, die als scheinbar angesehen werden. Die Resultate unterstreichen die Notwendigkeit, seismische Anisotropie zu berücksichtigen, um tensile Komponenten von Erdbeben eindeutig zu bestimmen. Ein daher neu entwickelter Inversionsalgorithmus wurde getestet und erfolgreich auf 112 Erdbeben der letzten intensiven Schwarmepisode im Jahr 2000 im Vogtland an der deutsch-tschechischen Grenze angewandt. Die Herdparameter wurden unter Verwendung isotroper und anisotroper Geschwindigkeitsmodelle ermittelt. Die Beben zeigen Lokalmagnituden zwischen 1,6 und 3,2. Sie weisen zueinander ähnliche Herdflächenlösungen mit linkslateralem Versatz auf steil einfallenden, etwa Nord-Süd orientierten Bruchflächen auf. Die Fallwinkel nehmen oberhalb 8,4km Tiefe ab. Für über 60% der betrachteten Erdbeben werden tensile Quellanteile mit Volumenvergrößerung beobachtet. Die tensilen Komponenten zeigen Abhängigkeiten von Herdzeit und -ort. Sie sind zu Beginn des Schwarms sowie in Tiefen unterhalb 8,4km besonders signifikant und nehmen später an Bedeutung ab. Abgeleitete Hauptspannungsachsen enthalten P Achsen mit nordwestlicher und T Achsen mit südwestlicher Streichrichtung. Sie ähneln denen in Mitteleuropa. Es werden tiefenabhängige Fallwinkel beobachtet. Die Änderungen erfolgen für die P Achsen graduell von 50° hin zu flacheren Fallwinkeln bei tieferen Beben. Sie erfolgen jedoch abrupt für die T Achsen von etwa 8° oberhalb einer Tiefe von etwa 8,4km zu 21° einfallend unterhalb dessen. Mit dieser Arbeit werden erstmals zeitliche und räumliche Veränderungen tensiler Quellanteile und Spannungszustände im Vogtland für Erdbeben im Jahr 2000 beobachtet. Diese haben auch dann Bestand, wenn seismische Anisotropie berücksichtigt wird. Sie können durch Fluide erklärt werden, die in die Bruchflächen eindringen. KW - Seismologie KW - Momententensor KW - Anisotropie KW - Erdbeben KW - Wellenausbreitung KW - Vogtland KW - Schwarmbeben KW - tensile Anteile KW - Hauptspannungsachse KW - earthquake swarm KW - anisotropy KW - tensile earthquakes KW - moment tensor KW - West Bohemia Y1 - 2006 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-7758 ER - TY - GEN A1 - Rößler, Dirk A1 - Krüger, Frank A1 - Ohrnberger, Matthias T1 - Rupture propagation of recent large TsE off-coast Sumatra and Java N2 - The spatio-temporal evolution of the three recent tsunamogenic earthquakes (TsE) off-coast N-Sumatra (Mw9.3), 28/03/2005 (Mw8.5) off-coast Nias, on 17/07/2006 (Mw7.7) off-coast Java. Start time, duration, and propagation of the rupture are retrieved. All parameters can be obtained rapidly after recording of the first-arrival phases in near-real time processing. We exploit semblance analysis, backpropagation and broad-band seismograms within 30°-95° distance. Image enhancement is reached by stacking the semblance of arrays within different directions. For the three events, the rupture extends over about 1150, 150, and 200km, respectively. The events in 2004, 2005, and 2006 had source durations of at least 480s, 120s, and 180s, respectively. We observe unilateral rupture propagation for all events except for the rupture onset and the Nias event, where there is evidence for a bilateral start of the rupture. Whereas average rupture speed of the events in 2004 and 2005 is in the order of the S-wave speed (≈2.5-3km/s), unusually slow rupturing (≈1.5 km/s) is indicated for the July 2006 event. For the July 2006 event we find rupturing of a 200 x 100 km wide area in at least 2 phases with propagation from NW to SE. The event has some characteristics of a circular rupture followed by unilateral faulting with change in slip rate. Fault area and aftershock distribution coincide. Spatial and temporal resolution are frequency dependent. Studies of a Mw6.0 earthquake on 2006/09/21 and one synthetic source show a ≈1° limit in resolution. Retrieved source area, source duration as well as peak values for semblance and beam power generally increase with the size of the earthquake making possible an automatic detection and classification of large and small earthquakes. KW - Tsunami KW - Erdbeben KW - Indischer Ozean KW - Bruchausbreitung KW - Seismologie KW - Tsunami KW - Earthquake KW - Indonesia KW - Indian Ocean KW - Rupture Propagation KW - Seismology Y1 - 2007 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-13039 ER - TY - JOUR A1 - Rößler, Dirk A1 - Krüger, Frank A1 - Ohrnberger, Matthias T1 - Rupture propagation of the TsE (Mw7.7) on 17 July 2006 off-coast Java N2 - The Mw=7.7 tsunamogenic earthquake (TsE) on 17 July 2006, 08:19:28 shock the Indian Ocean at about 15 km depth off-coast Java, Indonesia. It caused a local tsunami with wave heights exceeding 2 m. The death toll reached several hundred. Thousands of people were displaced. By means of standard array methods, we have investigated the propagation and the extent of the rupture front of the causative earthquake. Waveform similarity is expressed by means of the semblance. We back-propagate the semblance for first-arrival phases recorded at broad-band stations within teleseismic distances (30°-95°). Image enhancement is realised by stacking the semblance of 8 arrays within different epicentral and azimuthal directions. From teleseismic observations we find rupturing of a 200 x 100 km wide area in at least 2 phases with propagation from NW to SE and source duration >125 s. The event has some characteristics of a circular rupture followed by unilateral faulting with change in slip rate. Unusually slow rupturing (≈1.5 km/s) is indicated. Fault area and aftershock distribution coincide. Spatial and temporal resolution are frequency dependent. Studies of a Mw6.0 earthquake on 2006/09/21 and one synthetic source show a ≈1° limit in resolution. Retrieved source area, source duration as well as peak values for semblance and beam power increase with the size of the earthquake making possible an automatic detection and classification of large and small earthquakes. KW - Seismologie KW - Erdbeben KW - Tsunami KW - Indischer Ozean Y1 - 2007 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-12964 ER - TY - JOUR A1 - Rößler, Dirk A1 - Krüger, Frank A1 - Rümpker, Georg A1 - Psencik, Ivan T1 - Tensile source components of swarm events in West Bohemia in 2000 by considering seismic anisotropy N2 - Earthquake swarms occur frequently in West Bohemia, Central Europe. Their occurrence is correlated with and propably triggered by fluids that escape on the earth's surface near the epicentres. These fluids raise up periodically from a seemingbly deep-seated source in the upper mantle. Moment tensors for swarm events in 1997 indicate tensile faulting. However, they were determined under assumption of seismic isotropy although anisotropy can be observed. Anisotropy may obscure moment tensors and their interpretation. In 2000, more than 10,000 swarm earthquakes occurred near Novy Kostel, West Bohemia. Event triggering by fluid injection is likely. Activity lasted from 28/08 until 31/12/00 (9 phases) with maximum ML=3.2. High quality P-wave seismograms were used to retrieve the source mechanisms for 112 events between 28/08/00 and 30/10/00 using > 20 stations. We determine the source geometry using a new algorithm and different velocity models including anisotropy. From inversions of P waves we observe ML<3.2, strike-slip events on steep N-S oriented faults with additional normal or reverse components. Tensile components seem to be evident for more than 60% of the processed swarm events in West Bohemia during the phases 1-7. Being most significant at great depths and at phases 1-4 during the swarm they are time and location dependent. Although tensile components are reduced when anisotropy is assumed they persist and seem to be important. They can be explained by pore-pressure changes due to the injection of fluids that raise up. Our findings agree with other observations e.g. correlation of fluid transport and seismicity, variations in b-value, forcing rate, and in pore pressure diffusion. Tests of our results show their significance. KW - Seismologie KW - Erbeben KW - Momententensor KW - Anisotropie KW - Vogtland KW - Seismology KW - tensile earthquake KW - moment tensor KW - anisotropy KW - West Bohemia Y1 - 2006 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-12975 ER - TY - THES A1 - Sen, Ali Tolga T1 - Inversion of seismic source parameters for weak mining-induced and natural earthquakes T1 - Inversion von seismischen Quellparametern für schwache Bergbau-induzierte und natürliche Erdbeben N2 - The purpose of this thesis is to develop an automated inversion scheme to derive point and finite source parameters for weak earthquakes, here intended with the unusual meaning of earthquakes with magnitudes at the limit or below the bottom magnitude threshold of standard source inversion routines. The adopted inversion approaches entirely rely on existing inversion software, the methodological work mostly targeting the development and tuning of optimized inversion flows. The resulting inversion scheme is tested for very different datasets, and thus allows the discussion on the source inversion problem at different scales. In the first application, dealing with mining induced seismicity, the source parameters determination is addressed at a local scale, with source-sensor distance of less than 3 km. In this context, weak seismicity corresponds to event below magnitude MW 2.0, which are rarely target of automated source inversion routines. The second application considers a regional dataset, namely the aftershock sequence of the 2010 Maule earthquake (Chile), using broadband stations at regional distances, below 300 km. In this case, the magnitude range of the target aftershocks range down to MW 4.0. This dataset is here considered as a weak seismicity case, since the analysis of such moderate seismicity is generally investigated only by moment tensor inversion routines, with no attempt to resolve source duration or finite source parameters. In this work, automated multi-step inversion schemes are applied to both datasets with the aim of resolving point source parameters, both using double couple (DC) and full moment tensor (MT) models, source duration and finite source parameters. A major result of the analysis of weaker events is the increased size of resulting moment tensor catalogues, which interpretation may become not trivial. For this reason, a novel focal mechanism clustering approach is used to automatically classify focal mechanisms, allowing the investigation of the most relevant and repetitive rupture features. The inversion of the mining induced seismicity dataset reveals the repetitive occurrence of similar rupture processes, where the source geometry is controlled by the shape of the mined panel. Moreover, moment tensor solutions indicate a significant contribution of tensile processes. Also the second application highlights some characteristic geometrical features of the fault planes, which show a general consistency with the orientation of the slab. The additional inversion for source duration allowed to verify the empirical correlation for moment normalized earthquakes in subduction zones among a decreasing rupture duration with increasing source depth, which was so far only observed for larger events. N2 - Ziel dieser Dissertation ist es ein automatisches Inversionsschema zur Ableitung von Parametern für Punkt- und finite Quellen von "schwachen" Erdbeben im Sinne von aussergewöhnlich kleinen Erdbeben an oder unterhalb der Magnitudenschwelle für Standard-Inversionsroutinen zu entwickeln. Der verwendete Inversionsansatz baut auf bestehender Inversionssoftware auf, die methodische Arbeit zielt hauptsächlich auf die Entwicklung und das Tuning des optimierten Inversionsablaufs. Das resultierende Inversionsverfahren ist für verschiedene Datensätze getestet und erlaubt somit eine Diskussion über das Herdinversionsproblem auf unterschiedlichen Skalierungen. In der ersten Anwendung, die sich mit induzierter Seismizität in Minen beschäftigt, werden die Herdparameter im lokalen Bereich mit einer Herd-Empfänger Entfernung von weniger als 3 km bestimmt. In diesem Zusammenhang bezieht sich die "schwache" Seismizität auf Ereignisse mit Magnituden unter Mw 2.0, welche sehr selten das Ziel von automatisierten Herdinversionsroutinen sind. Die zweite Anwendung berücksichtigt den regionalen Datensatz der Nachbebesequenz des Maule (Chile) Erdbebens von 2010, aufgezeichnet durch Breitband-Stationen in einer maximalen Entfernung von 300 km. In diesem Fall umfasst der Magnitudenbereich Ereignisse mit Magnituden kleiner als Mw 4.0. Auch dieser Datensatz wird hier als "schwache" Seismizität interpretiert, da die Analyse solch moderater Seismizität normalerweise nur durch Routinen für Momententensorinversion erfolgt, welche keinerlei Auflösung für die Herddauer oder finite Herdparameter beinhaltet. In dieser Arbeit wird ein automatisiertes, mehrstufiges Inversionsverfahren auf beide Datensätze angewendet mit dem Ziel sowohl Punktquellen-Herdparameter für double-couple (DC) und volle Momententensor (MT) Modelle als auch für Herddauer und finite Herdparameter zu bestimmen. Ein Hauptergebnis dieser Analyse von schwachen Ereignissen ist die Vergrößerung von Momententensor-Katalogen, deren Interpretation wichtig werden könnte. Aus diesem Grund wurde ein neuartiger Cluster-Ansatz für Herdmechanismen angewendet um diese automatisch zu klassifizieren, der somit die Untersuchung der wichtigsten und sich wiederholenden Bruchprozessen erlaubt. Die Inversion des Datensatzes von induzierter Seismizität in Minen enthüllt das sich wiederholende Auftreten von ähnlichen Bruchprozessen, wobei die Herdgeometrie von der Form der Minentäfelung kontrolliert wird. Darüber hinaus deuten die Momententensorlösungen einen signifikanten Beitrag von tensilen Prozessen an. Auch die zweite Anwendung zeigt einige charackteristische, geometrische Besonderheiten der Bruchflächen auf, welche prinzipiell konsistent zur Orientierung der subduzierten Platte sind. Die zusätzliche Inversion nach der Herdzeit erlaubt die Verifizierung von empirischen Korrelationen von momentennormalisierten Erdbeben in Subduktionszonen zwischen einer abnehmenden Bruchdauer mit zunehmender Herdtiefe, welche bisher nur für größere Ereignisse beobachtet wurde. KW - Seismologie KW - induzierten Seismizität KW - Quellenumkehr KW - Momententensor KW - Aufenthaltsdauer KW - seismology KW - induced seismicity KW - source inversion KW - moment tensor KW - source duration Y1 - 2014 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-71914 ER - TY - THES A1 - Zeckra, Martin T1 - Seismological and seismotectonic analysis of the northwestern Argentine Central Andean foreland N2 - After a severe M W 5.7 earthquake on October 17, 2015 in El Galpón in the province of Salta NW Argentina, I installed a local seismological network around the estimated epicenter. The network covered an area characterized by inherited Cretaceous normal faults and neotectonic faults with unknown recurrence intervals, some of which may have been reactivated normal faults. The 13 three-component seismic stations recorded data continuously for 15 months. The 2015 earthquake took place in the Santa Bárbara System of the Andean foreland, at about 17km depth. This region is the easternmost morphostructural region of the central Andes. As a part of the broken foreland, it is bounded to the north by the Subandes fold-and-thrust belt and the Sierras Pampeanas to the south; to the east lies the Chaco-Paraná basin. A multi-stage morphotectonic evolution with thick-skinned basement uplift and coeval thin-skinned deformation in the intermontane basins is suggested for the study area. The release of stresses associated with the foreland deformation can result in strong earthquakes, as the study area is known for recurrent and historical, destructive earthquakes. The available continuous record reaches back in time, when the strongest event in 1692 (magnitude 7 or intensity IX) destroyed the city of Esteco. Destructive earthquakes and surface deformation are thus a hallmark of this part of the Andean foreland. With state-of-the-art Python packages (e.g. pyrocko, ObsPy), a semi-automatic approach is followed to analyze the collected continuous data of the seismological network. The resulting 1435 hypocenter locations consist of three different groups: 1.) local crustal earthquakes (nearly half of the events belong to this group), 2.) interplate activity, of regional distance in the slab of the Nazca-plate, and 3.) very deep earthquakes at about 600km depth. My major interest focused on the first event class. Those crustal events are partly aftershock events of the El Galpón earthquake and a second earthquake, in the south of the same fault. Further events can be considered as background seismicity of other faults within the study area. Strikingly, the seismogenic zone encompass the whole crust and propagates brittle deformation down, close to the Moho. From the collected seismological data, a local seismic velocity model is estimated, using VELEST. After the execution of various stability tests, the robust minimum 1D-velocity model implies guiding values for the composition of the local, subsurface structure of the crust. Afterwards, performing a hypocenter relocation enables the assignment of individual earthquakes to aftershock clusters or extended seismotectonic structures. This allows the mapping of previously unknown seismogenic faults. Finally, focal mechanisms are modeled for events with acurately located hypocenters, using the newly derived local velocity model. A compressive regime is attested by the majority of focal mechanisms, while the strike direction of the individual seismogenic structures is in agreement with the overall north – south orientation of the Central Andes, its mountain front, and individual mountain ranges in the southern Santa-Bárbara-System. N2 - Nach einem schweren Erdbeben der Magnitude M W 5.7 am 17. Oktober 2015 in El Galpón, in der Provinz Salta im Nordwesten Argentiniens, habe ich ein lokales seismologisches Netzwerk, um das vermutete Epizentrum herum, aufgebaut. Dabei haben 13 Stationen kontinuierlich für 15 Monate gemessen. Das Netzwerk wurde in einem Gebiet installiert, welches durch tektonische Störungen charakterisiert ist, die entweder in der Kreidezeit zunächst als Abschiebungen initiiert und später als Aufschiebungen reaktiviert wurden oder in der geologischen jüngeren Vergangenheit erst entstanden sind. Die Intervallzeiten zwischen zwei Erdbeben sind dabei häufig unbekannt. Das Erdbeben von 2015 trat im Santa-Barbara-System im Argentinischen Vorland, 17 km unter der Erdoberfläche auf. Diese Region ist die östlichste strukturgeologische Provinz der Zentralanden und dem broken-foreland-Typus zuzuordnen. Im Norden schließt sich der Bolivianische Faltengürtel (Sierras Subandinas) und im Süden die Sierras Pampeanas an; im Osten liegt das Chaco-Paraná Becken. Eine mehrstufige morphotektonische Entwicklung wird hier vermutet, bei der das Grundgebirge durch als thick-skinned bezeichnete Deformation herausgehoben wurde und die dazwischen liegenden Intermontanbecken gleichzeitig Deformation des Typs thin-skinned erfahren haben. Die plötzliche Spannungsfreisetzung, die mit dieser Vorlanddeformation einhergeht, kann zu starken Erdbeben führen. Das Untersuchungsgebiet ist für wiederkehrende und historische, zerstörerische Erdbeben bekannt. Der zur Verfügung stehenden Aufzeichnungen reichen bis in das Jahr 1692 zurück, als ein Erdbeben der Magnitude M 7 (oder Intensität IX) die Stadt Esteco zerstörte. Daher sind zerstörerische Erdbeben ein besonderes Kennzeichen in diesem Teil des Andenvorlands. Für die Analyse der im seismologischen Netzwerk aufgezeichneten kontinuierlichen Daten wurde ein semiautomatischer Ansatz verfolgt, der mittels hochmoderner Python-Bibliotheken informationstechnisch umgesetzt wurde. Die resultierenden 1435 Erdbeben bestehen aus drei verschiedenen Gruppen: 1.) lokale Erdbeben in der Erdkruste (die etwa die Hälfte aller Events ausmachen), 2.) weiter entfernte Interplattenaktivität, die durch die Subduktion der Nazca-Platte unter den Südamerikanischen Kontinent hervorgerufen wird und 3.) sehr tiefen Erdbeben in etwa 600 km Tiefe. Mein Hauptaugenmerk lag dabei auf der ersten Gruppe. Diese krustalen Ereignisse sind teilweise Nachbeben des El Galpón Erdbebens und eines weiteren Bebens, welches sich weiter im Süden an der gleichen Störung ereignete. Die restlichen Beben können der allgemeinen Hintergrundaktivität entlang weiterer Störungen im Untersuchungsgebiet zugeschrieben werden. Beachtenswert ist dabei der Umstand, dass die Erdbebenaktivität in der gesamten Kruste beobachtet werden kann und sich dadurch die Deformation bis a fast an den Erdmantel ausbreitet. Mit den gesammelten Daten kann, unter der Verwendung der VELEST Software, ein lokales seismisches Geschwindigkeitsmodell bestimmt werden. Nach der Durchführung verschiedener Stabilitätstests, können aus dem robusten eindimensionalen Modell Richtwerte für die Zusammensetzung und den Aufbau der Erdkruste gewonnen werden. Dieanschließende Relokalisierung von Erdbebenherden erlaubt die Zuordnung einzelner Erdbeben zu Erdbebenclustern oder ausgedehnten seismotektonischen Strukturen. Dadurch können sogar zuvor unbekannte seismogene Störungen kartiert werden. Schlussendlich, werden Herdflächenlösungen für Beben mit präzise lokalisierten Erdbebenherden und unter der Einbeziehung des neu bestimmten lokalen Geschwindigkeitsmodells modelliert. Der Großteil der resultierenden Lösungen bestätigt das vorherrschende kompressive Regime. Das Streichen der einzelnen seismogenen Strukturen stimmt größtenteils mit der allgemeinen Nord – Süd Ausrichtugn der Zentralanden, ihrer Gebirgsfront und den einzelnen Höhenzügen im Santa-Barbará-System überein. T2 - Seismologische und Seismotektonische Analyse des Vorlandsystems der nordwestargentinischen Zentralanden KW - Seismology KW - Seismotektonik KW - Geophysics KW - Andes KW - Geosciences KW - Argentina KW - Anden KW - Seismologie KW - Geophysik KW - Geowissenschaften KW - Argentinien Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-473240 ER -