TY - THES A1 - Cattania, Camilla T1 - Improvement of aftershock models based on Coulomb stress changes and rate-and-state dependent friction T1 - Verbesserte Nachbebenmodelle durch Berücksichtigung von Coulombspannungsänderungen und Rate-State abhängiger Reibung N2 - Earthquake clustering has proven the most useful tool to forecast changes in seismicity rates in the short and medium term (hours to months), and efforts are currently being made to extend the scope of such models to operational earthquake forecasting. The overarching goal of the research presented in this thesis is to improve physics-based earthquake forecasts, with a focus on aftershock sequences. Physical models of triggered seismicity are based on the redistribution of stresses in the crust, coupled with the rate-and-state constitutive law proposed by Dieterich to calculate changes in seismicity rate. This type of models are known as Coulomb- rate and-state (CRS) models. In spite of the success of the Coulomb hypothesis, CRS models typically performed poorly in comparison to statistical ones, and they have been underepresented in the operational forecasting context. In this thesis, I address some of these issues, and in particular these questions: (1) How can we realistically model the uncertainties and heterogeneity of the mainshock stress field? (2) What is the effect of time dependent stresses in the postseismic phase on seismicity? I focus on two case studies from different tectonic settings: the Mw 9.0 Tohoku megathrust and the Mw 6.0 Parkfield strike slip earthquake. I study aleatoric uncertainties using a Monte Carlo method. I find that the existence of multiple receiver faults is the most important source of intrinsic stress heterogeneity, and CRS models perform better when this variability is taken into account. Epistemic uncertainties inherited from the slip models also have a significant impact on the forecast, and I find that an ensemble model based on several slip distributions outperforms most individual models. I address the role of postseismic stresses due to aseismic slip on the mainshock fault (afterslip) and to the redistribution of stresses by previous aftershocks (secondary triggering). I find that modeling secondary triggering improves model performance. The effect of afterslip is less clear, and difficult to assess for near-fault aftershocks due to the large uncertainties of the afterslip models. Off-fault events, on the other hand, are less sensitive to the details of the slip distribution: I find that following the Tohoku earthquake, afterslip promotes seismicity in the Fukushima region. To evaluate the performance of the improved CRS models in a pseudo-operational context, I submitted them for independent testing to a collaborative experiment carried out by CSEP for the 2010-2012 Canterbury sequence. Preliminary results indicate that physical models generally perform well compared to statistical ones, suggesting that CRS models may have a role to play in the future of operational forecasting. To facilitate efforts in this direction, and to enable future studies of earthquake triggering by time dependent processes, I have made the code open source. In the final part of this thesis I summarize the capabilities of the program and outline technical aspects regarding performance and parallelization strategies. N2 - Die örtliche und zeitlich Häufung von Erdbeben ist geeignet, um Änderungen in Seismizitätsraten auf kurzen bis mittleren Zeitskalen (Stunden bis Monate) zu prognostizieren. Kürzlich wurden vermehrt Anstrengungen unternommen, den Umfang solcher Modelle auf Operationelle Erdbebenvorhersage auszudehnen, welche die Veröffentlichung von Erdbebenwahrscheinlichkeiten beinhaltet mit dem Ziel, die Bevölkerung besser auf mögliche Erdbeben vorzubereiten. Das vorrangige Ziel dieser Dissertation ist die Verbesserung von kurz- und mittelfristiger Erdbebenprognose basierend auf physikalischen Modellen. Ich konzentriere mich hier auf Nachbebensequenzen. Physikalische Modelle, die getriggerte Seimizität erklären, basieren auf der Umverteilung von Spannungen in der Erdkruste. Berechnung der Coulomb Spannung können kombiniert werden mit dem konstituivem Gesetz von Dieterich, welches die Berechnung von Änderungen in der Seismizitätsrate ermöglicht. Diese Modelle sind als Coulomb-Rate-and-State (CRS) Modelle bekannt. Trotz der erfolgreichen Überprüfung der Coulomb-Hypothese, schneiden CRS-Modelle im Vergleich mit statistischen Modellen schlecht ab, und wurden deshalb bisher kaum im Kontext operationeller Erdbenbenvorhersage genutzt. In dieser Arbeit, gehe ich auf einige der auftretenden Probleme ein. Im Besonderen wende ich mich den folgenden Fragen zu: (1) Wie können wir die Unsicherheiten und die Heterogenität des Spannungsfeldes infolge des Hauptbebens realistisch modellieren? (2)Welche Auswirkungen haben zeitlich variable Spannungsänderungen in der postseismischen Phase? Ich konzentriere mich hierbei auf zwei Beispiele in unterschiedlichen tektonischen Regionen: die Aufschiebung des Mw9.0 Tohoku Erdbeben und die Blattverschiebung des Mw6.0 Parkfield Erdbeben. Ich untersuche aleotorische Unsicherheiten der Coulomb-Spannung durch Variabilität in der Orientierung der betroffenen Bruchflächen und durch Spannungsgradienten innerhalb von Modellzellen. Ich zeige, dass die Existenz der unterschiedlichen Bruchflächen die bedeutenste Quelle für intrinsiche Spannungheterogenität ist und das CRS-Modelle deutlich besser abschneiden, wenn diese Variabilität berücksichtigt wird. Die epistemischen Unsicherheiten aufgrund von unterschiedlichen Ergebnissen von Inversionen von Daten für die Verschiebung entlang der Bruchfläche haben ebenso erhebliche Auswirkungen auf die Vorhersage. Ich gehe dann auf die Rolle von postseismischen Spannung ein, insbesondere auf zwei Prozesse: aseismische Verschiebung entlang der Störungsfläche des Hauptbebens (Afterslip) und die Veränderung von Spannungen durch vorhergehende Nachbeben (sekundäres Triggern). Ich demonstriere, dass das Modellieren von sekundärem Triggern die Modellvorhersage in beiden Fallbeispielen verbessert. Die Einbeziehung von Afterslip verbessert die Qualität der Vorhersage nur für die Nachbebensequenz des Parkfield Erdbebens. Dagegen kann ich nachweisen, dass Afterslip infolge des Tohoku Bebens eine höhere Seismizität auf Abschiebungsflächen im Hangenden begünstigt. Die dargestellten Verbesserungen des CRS-Modells sind sehr vielversprechend im Kontext operationeller Erdbebenvorhersage, verlangen aber nach weiterer Überprüfung. Ich stelle die vorläufigen Ergebnisse eines gemeinschaftlichen Tests für die Erdbebenfolge von Canterbury 2010-2012 vor, welcher von CSEP durchgeführt wurde. Die physikalischen Modelle schneiden hier im Vergleich mit statistischen Modellen gut ab. Daher scheint eine Anwendung von CSR-Modellen, die Unsicherheiten und sekundäres Triggering berücksichtigen, in zukünftigen operationellen Erdbebenvorhersagen empfehlenswert. Um die Bemühungen in dieser Richtung zu unterstützen und weitere Studien zum Triggern von Erdbeben durch zeitabhängige Prozesse zu ermöglichen, habe ich meinen Open Source Code öffentlich zugänglich gemacht. Im letzen Teil dieser Arbeit fasse ich die Leistungsfähigkeit des Programms zusammen und skizziere die technischen Aspekte bezüglich der Effiziens und der Parallelisierung des Programmes. KW - earthquake forecasting KW - earthquake interaction KW - Coulomb stress KW - rate-state friction KW - Erdbebenvorhersage KW - Coulombspannung KW - Erdbebenwechselwirkung Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-87097 ER - TY - THES A1 - Hainzl, Sebastian T1 - Earthquake triggering and interaction T1 - Erdbebenentstehung und Wechselwirkungen N2 - Earthquake faults interact with each other in many different ways and hence earthquakes cannot be treated as individual independent events. Although earthquake interactions generally lead to a complex evolution of the crustal stress field, it does not necessarily mean that the earthquake occurrence becomes random and completely unpredictable. In particular, the interplay between earthquakes can rather explain the occurrence of pronounced characteristics such as periods of accelerated and depressed seismicity (seismic quiescence) as well as spatiotemporal earthquake clustering (swarms and aftershock sequences). Ignoring the time-dependence of the process by looking at time-averaged values – as largely done in standard procedures of seismic hazard assessment – can thus lead to erroneous estimations not only of the activity level of future earthquakes but also of their spatial distribution. Therefore, it exists an urgent need for applicable time-dependent models. In my work, I aimed at better understanding and characterization of the earthquake interactions in order to improve seismic hazard estimations. For this purpose, I studied seismicity patterns on spatial scales ranging from hydraulic fracture experiments (meter to kilometer) to fault system size (hundreds of kilometers), while the temporal scale of interest varied from the immediate aftershock activity (minutes to months) to seismic cycles (tens to thousands of years). My studies revealed a number of new characteristics of fluid-induced and stress-triggered earthquake clustering as well as precursory phenomena in earthquake cycles. Data analysis of earthquake and deformation data were accompanied by statistical and physics-based model simulations which allow a better understanding of the role of structural heterogeneities, stress changes, afterslip and fluid flow. Finally, new strategies and methods have been developed and tested which help to improve seismic hazard estimations by taking the time-dependence of the earthquake process appropriately into account. N2 - Erdbeben interagieren in vielfältiger Weise miteinander, weshalb sie nicht als einzelne, unabhängige Ereignisse behandelt werden können. Obwohl diese Erdbebenwechselwirkungen in der Regel zu einer komplexen Entwicklung des Spannungsfelds führen, bedeutet dies nicht zwangsläufig, dass Erdbeben rein zufällig und völlig unberechenbar auftreten. Insbesondere kann das Zusammenspiel zwischen Erdbeben zu ausgeprägten Charakteristiken wie Phasen beschleunigter Aktivität, seismischer Ruhe sowie raumzeitlichen Erdbebenanhäufungen (Schwärme und Nachbebensequenzen) führen. Die Vernachlässigung der Zeitabhängigkeit des Erdbebenprozesses kann somit zu fehlerhaften Einschätzungen nicht nur des zukünftigen Aktivitätsniveaus, sondern auch der räumlichen Verteilung führen. Daher besteht ein dringender Bedarf an geeigneten zeitabhängigen Seismizitätsmodellen. Meine Arbeit zielt auf ein verbessertes Verständnis und Charakterisierung der Interaktionen von Erdbeben ab, um Abschätzungen der Erdbebengefährdung zu verbessern. Zu diesem Zweck untersuche ich Seismizitätsmuster auf den räumlichen Skalen von hydraulisch induzierten Öffnungsbrüchen (Meter bis Kilometer) bis zu Verwerfungssystemen (Hunderte von Kilometern), während die zeitlichen Skalen von Nachbebenaktivität (Minuten bis Monate) bis zu seismischen Zyklen (bis zu mehrere tausendend Jahre) reichen. Meine Studien ergeben eine Reihe neuer Merkmale von Fluid- und Spannungs-induzierten Erdbeben. Ergänzend zur reinen Datenanalyse der Erdbeben- und Deformationsdaten liefern statistische und Physik-basierte Modellsimulationen ein besseres Verständnis der Rolle von strukturellen Heterogenitäten, Spannungsänderungen und postseismischen Prozessen. Schließlich konnten neue Strategien und Methoden entwickelt und getestet werden, mit denen die Erdbebengefährdung besser eingeschätzt werden kann, indem die Zeitabhängigkeit des Erdbebens Prozess angemessen berücksichtigt wird. KW - Erdbebeninteraktion KW - Erdbebengefährdung KW - earthquake interaction KW - seismic hazard Y1 - 2011 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-50095 ER -