Automated inversion of long period signals from shallow volcanic and induced seismic events

Vulkanische seismische Ereignisse, die langperiodische Signale erzeugen (LP-Events), sind wichtige Indikatoren fur im Erdinneren stattfindende Fluidbewegungen. Ihre Untersuchung ist essentiell, um Gefahrdungen fur die Umgebung von Vulkangebieten abschatzen zu konnen. LP-Events treten nicht nur an Vulkanen auf, sie sind auch in der Umgebung von Fluid-Lagerstatten und Bergwerken zu beobachten. Dort werden sie als Anzeichen fur das Auftreten langsam ablaufender Deformationsprozesse gedeutet, welche in Zusammenhang mit strukturellen Inhomogenitaten stehen konnen und durch ausere Einwirkungen ausgelost werden (induzierte Ereignisse). Die Interpretation solcher Ereignisse kann einen Beitrag zu seismischen Gefahrdungsanalysen liefern. Bei der Untersuchung von LP-Events hat man auf der einen Seite die lang erforschten lokalen vulkanischen Ereignisse, welche vulkantypische Charakteristika aufweisen konnen und in von starker Topographie und heterogenen Untergrundstrukturen gepragten Regionen vorkommen. Auf der anderen Seite findet man induzierte LP-Events, welche noch nicht der Gegenstand systematischer Untersuchungen waren. Es existiert einerseits Software welche zur kontinuierlichen Uberwachung und Auswertung von seismischen Signalen bei vulkanischen Krisen gedacht ist. Andererseits gibt es Programme, welche gros- und intermediarskalige seismische Ereignisse fur bestimmte Regionen analysieren und auftretende Ereignisse katalogisieren. Bis zum aktuellen Zeitpunkt gibt es keinen systematischen Ansatz, die Untersuchungen von vulkanischen und induzierten seismischen LP-Events zu kombinieren. Im Rahmen der vorliegenden Arbeit wurde ein Programm entwickelt, welches es erlaubt, langperiodische seismische Signale schnell und unabhangig von der vorliegenden Geometrie der betreffenden Region zu analysieren. Der zugrundeliegende Algorithmus basiert auf einer algebraischen Inversion seismischer Daten in einem Einschrittverfahren in Kombination mit einer dreidimensionalen raumlichen Gittersuche bezuglich moglicher Quellpositionen. Um die Inversion ausfuhren zu konnen, werden 18 Green’s-Funktionen (GFn) benotigt. Die GFn beinhalten Informationen uber das lokale oder regional Geschwindigkeitsmodell. Sie werden fur die betreffende Region im Voraus berechnet und in einer Dateistruktur bereitgestellt, was eine sehr schnelle Inversion der Daten ermoglicht. Im Gegensatz zu anderen Methoden fur die Inversion seismischer Daten ist im vorliegenden Ansatz die vollstandige dreidimensionale Untergrundstruktur in den GFn enthalten, dies beinhaltet auch Topographie oder Anisotropien. Die seismischen Daten konnen zum einen ereignisbasiert analysiert werden. Hierbei invertiert der Algorithmus Daten fur eine durch ein externes (Detektions-)Programm bestimmte Quellzeit, Die Bestimmung dieser Quellzeit kann durch eine weitere Gittersuche innerhalb eines gegebenen Zeitintervalls verfeinert werden. Unter Verzicht auf externe Detektion von Ereignissen kann das entwickelte Programm zum anderen auch zur schnellen Untersuchung eines kontinuierlichen Datenstroms verwendet werden. Dies wird durch die schrittweise Anwendung in einem gleitenden Zeitfenster ermoglicht. Die fur jedes Zeitfenster bestimmte Losung erlaubt eine unabhangige Detektion von Ereignissen durch den Vergleich mit einem anwendungsfallabhangigen Schwellenwert. Der neuentwickelte Algorithmus wurde mit synthetischen und echten Daten getestet, Daten fur ein- und dreidimensionale Geschwindigkeitsmodelle wurden verwendet. Der Code ist lauffahig und auf Daten anwendbar. Die Ergebnisse sind konsistent, sowohl mit Resultaten von bekannten methodischen Untersuchungen als auch mit Parametern, die von anderen Programmen im Fall der Echtdatenstudien ermittelt wurden. Es zeigte sich, das die Topographie einen starken Einflus auf den Signalgehalt der GFn hat. Darum mussen vor nachfolgenden Anwendungen auf Falle mit komplizierter Geometrie weitere Tests durchgefuhrt werden, welche die genauen Auswirkungen verschiedener Topographien quantifizieren konnen. In dieser Arbeit wurden Fallstudien fur vulkanische und induzierte LP-Events bearbei- tet. Die Ergebnisse dieser Studien legen nahe, das numerische Verbesserungen sowie die Implementierung weiterer Regularisierungsparameter angeraten sind, um robuste Losun- gen zu gewahrleisten. Die allgemeine Funktionsfahigkeit des Algorithmus’ wurde demonstriert, sowohl fur die statische Inversion eines zeitlich begrenzten Datensatzes als auch fur die Anwendung auf kontinuierliche synthetische und reale Datensatze als selbstandige Inversions- und Detektionsroutine. Unter Einbeziehung der aufgefuhrten Verbesserungsmoglichkeiten kann der vorgestellte Algorithmus im Rahmen eines automatisierten, echtzeitnahen Monitoring-Systems verwendet werden. Volcanic seismic events that emit long period signals (LP-events) are considered import- ant indicators for current subsurface fluid or magma movement. The analysis of such events is crucial for the monitoring of volcanic regions in order to evaluate the evolution of volcanic unrest and the relative hazard assessment. In addition, induced LP-events oc- cur near fluid reservoirs and mining sites, where they are interpreted as evidence for slow rupturing or deformation processes. The continuous long-term evaluation of LP-signals in such regions can allow to infer local or regional structural inhomogeneities, which again can contribute to seismic hazard estimations. So far, the analyses of volcanic and induced seismic LP-events are strictly separated. Software tools are available, which are either specialised on the monitoring of volcanic unrest or restricted to the continuous analysis of large and medium scaled seismic events on regions with a flat topography. A new software tool has been developed in the scope of this thesis. It allows the fast and simple analysis of seismic LP-events, independent of the geometry and size of the source region. The algorithm is based on an algebraic one-step time-domain inversion of full waveform data, combined with a gridsearch for possible source locations. For conducting the inversion, a set of 18 Green’s functions (GFs) is used that take into account the local or regional seismic velocity model. The GFs are generated in advance to allow a fast data processing. In contrast to other approaches, the complete infomation about the three-dimensional structure of the underground is contained within the GFs – this includes topography as well as anisotropies. Seismic data can be analysed event-based, in which the source time can be determined by an additional gridsearch in time. Without the need for an external event detection routine it is furthermore possible to analyse continuous seismic data by applying a moving time-window approach. Data and estimated parameters are compared for successive overlapping time intervals. The level of agreement between them is expressed in one value and a pre-defined threshold is utilised for the detection of LP-events. The new algorithm has been tested under several considerations using synthetic data sets for one- and three-dimensional geometries. The code is operational and yields results that are consistent with former findings existing event analyses. Topography has a large influence on the GFs, therefore special care has to be taken of the generation of GFs in case of complex geometries. Case studies have been carried out for both volcanic and induced seismic LP-events. The results suggest that numerical improvements and the implementation of regularisa- tion parameters are advisable. The general capability of the inversion algorithm to act as an automated inversion detection routine on continuous data is demonstrated for both synthetic and real data. These tests show that further optimisation with respect to computational speed is re- quired. Overcoming these limitations will allow to use the developed algorithm within the framework of an automated monitoring system.