No description available.

Der Bericht dokumentiert Ort, Datum, Magnitude und Tiefe der Erdbeben in Sachsen über einen Zeitraum von drei Jahren. Bei der Auswertung der Registrierungen werden zusätzlich zu den Online-Stationen des seismologischen Sachsennetzes auch die Stationen der Nachbarländer einbezogen, deren Standorte in der Veröffentlichung dargestellt sind. Ein besonderes Ereignis stellt dabei der Erdbebenschwarm von Oktober bis Dezember 2008 dar. Das Kapitel »Induzierte Seismizität« befasst sich mit den Untersuchungen zur bergbauinduzierten Seismizität der Uranerzgrube Schlema-Alberoda und zur Böschungsbewegung des ehemaligen Braunkohletagebaues Nachterstedt.

Im Jahr 1806 verfasste Heinrich von Kleist die Kurzgeschichte Das Erdbeben in Chili, in der nach einem Erdbeben in Santiago de Chile ein junges Paar aus religiösen Gründen gelyncht wird. Das Problem des Leidens wird im Text deutlich thematisiert. In diesem Aufsatz wird versucht, durch eine chronologische Textanalyse auszulegen, ob Kleist in der Novelle das logische oder induktive Problem des Leidens oder beides thematisiere. Es wird argumentiert, dass das große Gewicht des unnötigen Leidens in der Geschichte und die fehlende Rolle Gottes als moralischer Agent von der Perspektive des Erzählers eine Thematisierung des induktiven Problems des Leidens andeuten.

Die relativ hohe seismische Aktivität der Tote-Meer-Störungszone (Dead Sea Fault Zone - DSFZ) ist mit einem hohen Gefahrenpotential verbunden, welches zu einem erheblichen Erdbebenrisiko für die Ballungszentren in den Ländern Syrien, Libanon, Palästina, Jordanien und Israel führt. Eine Vielzahl massiver, zerstörerischer Erdbeben hat sich in diesem Raum in den letzten zwei Jahrtausenden ereignet. Ihre Wiederholungsrate zeigt Anzeichen für eine zeitliche Abhängigkeit, insbesondere wenn lange Zeiträume in Betracht gezogen werden. Die Berücksichtigung der zeitlichen Abhängigkeit des Auftretens von Erdbeben ist für eine realistische seismische Gefährdungseinschätzung von großer Bedeutung. Ziel der vorliegenden Arbeit ist es, anhand des zeitabhängigen Auftretens von Erdbeben eine robuste wahrscheinlichkeitstheoretische seismische Gefährdungseinschätzung am Beispiel der DSFZ zu entwickeln. Mittels dieser Methode soll die zeitliche Abhängigkeit des Auftretens von großen Erdbeben (Mw ≥ 6) untersucht und somit eine Gefährdungseinschätzung für das Untersuchungsgebiet getroffen werden. Primär gilt es zu prüfen, ob das Auftreten von großen Erdbeben tatsächlich einer zeitlichen Abhängigkeit unterliegt und wenn ja, inwiefern diese bestimmt werden kann. Zu diesem Zweck werden insgesamt vier zeitabhängige statistische Verteilungen (Weibull, Gamma, Lognormal und Brownian Passage Time (BPT)) sowie die zeitunabhängige Exponentialverteilung (Poisson-Prozess) getestet. Zur Abschätzung der jeweiligen Modellparameter wird eine modifizierte Methode der gewichteten Maximum-Likelihood-Schätzung (MLE) verwendet. Um einzuschätzen, ob die Wiederholungsrate von Erdbeben einer unimodalen oder multimodalen Form folgt, wird ein nichtparametrischer Bootstrap-Test für Multimodalität durchgeführt. Im Falle einer multimodalen Form wird neben der MLE zusätzlich eine Erwartungsmaximierungsmethode (EM) herangezogen. Zur Auswahl des am besten geeigneten Modells wird zum einem das Bayesschen Informationskriterium (BIC) und zum anderen der modifizierte Kolmogorow-Smirnow-Goodness-of-Fit-Test angewendet. Abschließend werden mittels der Bootstrap-Methode die Konfidenzintervalle der geschätzten Parameter berechnet. Als Datengrundlage werden Erdbeben mit Mw ≥ 6 seit dem Jahre 300 n. Chr. herangezogen. Das Untersuchungsgebiet erstreckt sich von 29.5° N bis 37° N und umfasst ein ca. 40 km breites Gebiet entlang der DSFZ. Aufgrund der seismotektonischen Situation im Untersuchungsgebiet wird zwischen einer südlichen, zentralen und nördlichen Subzone unterschieden. Dabei kann die südliche Subzone aus Mangel an Daten nicht für die Analysen herangezogen werden. Die Ergebnisse für die zentrale Subzone zeigen keinen signifikanten multimodalen Verlauf der Wiederholungsrate von Erdbeben. Des Weiteren ist kein signifikanter Unterschied zwischen den zeitabhängigen und dem zeitunabhängigem Modell zu verzeichnen. Da das zeitunabhängige Modell vergleichsweise einfach interpretierbar ist, wird die Wiederholungsrate von Erdbeben in dieser Subzone unter Annahme der Exponentialverteilungs-Hypothese abgeschätzt. Sie wird demnach als zeitunabhängig betrachtet und beträgt 9.72 * 10-3 Erdbeben (mit Mw ≥ 6) pro Jahr. Einen besonderen Fall stellt die nördliche Subzone dar. In diesem Gebiet tritt im Durchschnitt alle 51 Jahre ein massives Erdbeben (Mw ≥ 6) auf. Das letzte Erdbeben dieser Größe ereignete sich 1872 und liegt somit bereits 137 Jahre zurück. Somit ist in diesem Gebiet ein Erdbeben dieser Stärke überfällig. Im statistischen Mittel liegt die Zeit zwischen zwei Erdbeben zu 96% unter 137 Jahren. Zudem wird eine deutliche zeitliche Abhängigkeit der Erdbeben-Wiederauftretensrate durch die Ergebnisse der in der Arbeit neu entwickelten statistischen Verfahren bestätigt. Dabei ist festzustellen, dass die Wiederholungsrate insbesondere kurz nach einem Erdbeben eine sehr hohe zeitliche Abhängigkeit aufweist. Am besten repräsentiert werden die seismischen Bedingungen in der genannten Subzone durch ein bi-modales Weibull-Weibull-Modell. Die Wiederholungsrate ist eine glatte Zeitfunktion, welche zwei Häufungen von Datenpunkten in der Zeit nach dem Erdbeben zeigt. Dabei umfasst die erste Häufung einen Zeitraum von 80 Jahren, ausgehend vom Zeitpunkt des jeweiligen Bebens. Innerhalb dieser Zeitspanne ist die Wiederholungsrate extrem zeitabhängig. Die Wiederholungsrate direkt nach einem Beben ist sehr niedrig und steigert sich in den folgenden 10 Jahren erheblich bis zu einem Maximum von 0.024 Erdbeben/Jahr. Anschließend sinkt die Rate und erreicht ihr Minimum nach weiteren 70 Jahren mit 0.0145 Erdbeben/Jahr. An dieses Minimum schließt sich die zweite Häufung von Daten an, dessen Dauer abhängig von der Erdbebenwiederholungszeit ist. Innerhalb dieses Zeitfensters nimmt die Erdbeben-Wiederauftretensrate annähernd konstant um 0.015 Erdbeben/Jahr zu. Diese Ergebnisse bilden die Grundlage für eine zeitabhängige probabilistische seismische Gefährdungseinschätzung (PSHA) für die seismische Quellregion, die den nördlichen Raum der DSFZ umfasst. / The seismicity of the Dead Sea fault zone (DSFZ) during the last two millennia is characterized by a number of damaging and partly devastating earthquakes. These events pose a considerable seismic hazard and seismic risk to Syria, Lebanon, Palestine, Jordan, and Israel. The occurrence rates for large earthquakes along the DSFZ show indications to temporal changes in the long-term view. The aim of this thesis is to find out, if the occurrence rates of large earthquakes (Mw ≥ 6) in different parts of the DSFZ are time-dependent and how. The results are applied to probabilistic seismic hazard assessments (PSHA) in the DSFZ and neighboring areas. Therefore, four time-dependent statistical models (distributions), including Weibull, Gamma, Lognormal and Brownian Passage Time (BPT), are applied beside the exponential distribution (Poisson process) as the classical time-independent model. In order to make sure, if the earthquake occurrence rate follows a unimodal or a multimodal form, a nonparametric bootstrap test of multimodality has been done. A modified method of weighted Maximum Likelihood Estimation (MLE) is applied to estimate the parameters of the models. For the multimodal cases, an Expectation Maximization (EM) method is used in addition to the MLE method. The selection of the best model is done by two methods; the Bayesian Information Criterion (BIC) as well as a modified Kolmogorov-Smirnov goodness-of-fit test. Finally, the confidence intervals of the estimated parameters corresponding to the candidate models are calculated, using the bootstrap confidence sets. In this thesis, earthquakes with Mw ≥ 6 along the DSFZ, with a width of about 20 km and inside 29.5° ≤ latitude ≤ 37° are considered as the dataset. The completeness of this dataset is calculated since 300 A.D. The DSFZ has been divided into three sub zones; the southern, the central and the northern sub zone respectively. The central and the northern sub zones have been investigated but not the southern sub zone, because of the lack of sufficient data. The results of the thesis for the central part of the DSFZ show that the earthquake occurrence rate does not significantly pursue a multimodal form. There is also no considerable difference between the time-dependent and time-independent models. Since the time-independent model is easier to interpret, the earthquake occurrence rate in this sub zone has been estimated under the exponential distribution assumption (Poisson process) and will be considered as time-independent with the amount of 9.72 * 10-3 events/year. The northern part of the DSFZ is a special case, where the last earthquake has occurred in 1872 (about 137 years ago). However, the mean recurrence time of Mw ≥ 6 events in this area is about 51 years. Moreover, about 96 percent of the observed earthquake inter-event times (the time between two successive earthquakes) in the dataset regarding to this sub zone are smaller than 137 years. Therefore, it is a zone with an overdue earthquake. The results for this sub zone verify that the earthquake occurrence rate is strongly time-dependent, especially shortly after an earthquake occurrence. A bimodal Weibull-Weibull model has been selected as the best fit for this sub zone. The earthquake occurrence rate, corresponding to the selected model, is a smooth function of time and reveals two clusters within the time after an earthquake occurrence. The first cluster begins right after an earthquake occurrence, lasts about 80 years, and is explicitly time-dependent. The occurrence rate, regarding to this cluster, is considerably lower right after an earthquake occurrence, increases strongly during the following ten years and reaches its maximum about 0.024 events/year, then decreases over the next 70 years to its minimum about 0.0145 events/year. The second cluster begins 80 years after an earthquake occurrence and lasts until the next earthquake occurs. The earthquake occurrence rate, corresponding to this cluster, increases extremely slowly, such as it can be considered as an almost constant rate about 0.015 events/year. The results are applied to calculate the time-dependent PSHA in the northern part of the DSFZ and neighbouring areas.

The phenomenon of induced seismicity is caused by anthropogenic activity such as: underground and opencast mining, extraction of conventional and unconventional hydrocarbons, construction of water reservoirs and production of geothermal energy. In recent years, interest in induced seismicity increased due to the fact that it causes increasingly stronger earthquakes, even above 4 on the Richter scale. Thus, it poses a threat to people, technical and urban infrastructure. This study analyzed the seismic event of M = 4.6, which occurred on the 15 September 2018 in the Rudna copper mine area in SW Poland. For this purpose, Sentinel 1 satellite data and DInSAR processing method were used to determine the ground movement values in the satellite line of sight. Based on the results for four image pairs, the area disturbed by the seismic event was determined. The maximum values of subsidence ranged from -65 mm to -75 mm depending on the analysed dataset and the area of deformation was determined at approx. 4 km sq. The results indicate the usefulness of the adopted method to determine ground deformation caused by induced seismicity in an underground mining area.

The main goal of CIC-BREL project (Cracked and Inelastic Calculation of BRacing Elements) is to develop an analytical method to distribute horizontal forces on bracing elements, in this case reinforced masonry shear walls, of a building considering the cracked and inelastic state of material. The moment curvature curve of the wall section is created first depending on the section geometry and material properties of both the masonry units and steel reinforcement. This curve will start with an elastic material behavior, then continue in inelastic material behavior where the masonry crushes and the steel start to yield, until the maximum bending moment M_p is reached. Due to reinforced masonry wall ductility, post maximum capacity is also considered assuming a maximum curvature of 0.1%. From the moment curvature curve, the force displacement curve could be extracted depending on the wall height and wall boundary conditions. Matrix formulation has been developed for both elastic and damaged stiffness matrix, considering different boundary conditions. Fixed-fixed boundary condition which usually exists at the middle stories or last story with strong top diaphragm, fixed-pinned which is the case of the last story that has a relatively soft top diaphragm, and pinned-fixed in the first story case. Other boundary conditions could be considered depending on the degree of fixation on the wall both ends at the top and the bottom. The matrix formulation combined with the force-displacement curve which considers different material stages (elastic, inelastic, ductile post peak force) is used to define forces in each bracing element even after elastic behavior. After elastic phase of each wall the stiffness of the element will degrade leading to a less portion of the total lateral force; other elastic walls, i.e., stronger walls, will receive more portion of the total force leading to a redistribution of the total force. This process will be iterated until the total force is distributed on each bracing element depending on the wall section state: elastic, inelastic and ductile post-peak capacity. Flowcharts clearly will show this process. Finally, a Fortran code is developed to show examples using this method. The developed analytical method will be verified by the results of shake table tests held at the University of California in San Diego, USA. Last test performed in the year 2018 uses T-section reinforced masonry walls, subjected to shakings with increased intensity. The total applied force for each shaking could be defined depending on the structural weight and shaking intensity (acceleration). The damage and displacement at each intensity has been recorded and evaluated. Depending on these test results, the results of the analytically developed method will be compared and evaluated. Total system displacement at different lateral load values has been compared for analytical calculations and shake table tests; furthermore, each wall state at increased load has been compared, good agreement could be noticed.:Acknowledgement 5 1. Introduction 7 1.1. State of the Art 9 1.2. Elastic Formulae 9 1.3. Example, Elastic Calculation 12 1.3.1. Stiffnesses of the System 13 1.3.2. Torsion due to Eccentric Lateral Loading 14 1.3.3. Distribution of the Lateral Load on Wall “j” and Floor “i” 15 2. Force Displacement Curve of RM Shear Wall 19 2.1. Introduction 19 2.2. Cantilever Wall 19 2.2.1. Cantilever Elastic Wall 19 2.2.2. Cantilever Inelastic Wall 21 2.2.3. Cantilever Post-Peak Wall 22 2.3. Fixed-Fixed Wall 23 2.3.1. Fixed-Fixed Elastic Wall 23 2.3.2. Fixed-Fixed Inelastic Wall 24 2.3.3. Fixed-Fixed Post-Peak Wall 26 2.4. Moment – Curvature Analysis 26 2.5. Example, Rectangle Cross Section, Cantilever 29 a) Moment Curvature Curve 29 b) Force Displacement Curve 32 2.6. Example, Rectangle Cross Section, Fixed-Fixed 33 a) Moment Curvature Curve 33 b) Force Displacement Curve 33 2.7. Example, T Cross Section, Cantilever 35 a) Moment Curvature Curve 35 b) Force Displacement Curve 41 2.8. Example, T Cross Section, Fixed-Fixed 43 a) Moment Curvature Curve 43 b) Force Displacement Curve 43 3. Matrix Formulation 47 3.1. Procedure 47 3.2. Structure Discretization 47 3.3. Element, i.e.; Wall, Local Stiffness Matrix 48 3.4. Stiffness Matrix of Fixed-Pinned Beam 52 3.4.1. Elastic 52 3.4.2. Pre-Peak Inelastic 54 3.4.3. Post-Peak Inelastic 55 3.4.4. Normal Force Part in the Stiffness Matrix 56 3.5. Stiffness Matrix of Pinned-Fixed Beam 57 3.5.1. Elastic 57 3.5.2. Post-Peak Inelastic 57 3.6. Stiffness Matrix of Fixed-Fixed Beam 58 3.6.1. Elastic 58 3.6.2. Post-Peak Inelastic 60 3.7. Summary of Stiffness Matrices 61 3.7.1. Fixed-Fixed 61 3.7.2. Fixed-Pinned 62 3.7.3. Pinned-Fixed 63 3.8. Transformation Matrix 63 3.9. Assemble the Structure Stiffness Matrix 65 3.10. Assemble the Structure Nodal Vector 66 3.11. Solve, Get Nodal Displacements and Forces 66 4. Matrix Formulation and Deformation Based Method 69 4.1. Elastic Method in Distributing Lateral Force 69 4.2. Elastic and Inelastic Method in Distributing Lateral Force 70 5. Shake Table Tests 73 5.1. Introduction 73 5.2. Design of Test Structure 73 5.3. Material Properties 75 5.4. Tests and Observations 75 5.4.1. Tests up to Mul-90% 76 5.4.2. Tests with Mul-120% 76 5.4.3. Tests with Mul-133% 76 5.5. Deformations 77 6. Verification 81 6.1. T Cross Section, Dimensions, Reinforcement and Materials 81 6.2. Moment Curvature Curve 82 6.3. Force Displacement Curve 85 6.4. Force Displacement Curve of the Structure 88 7. Conclusions and Suggestions 91 8. References 93 Appendix 1, Timoshenko Beam 95 • Fixed-Fixed 95 • Fixed-Pinned 95 • Pinned-Fixed 96 Appendix 2, Bernoulli Beam 97 • Fixed-Fixed 97 • Fixed-Pinned 97 • Pinned-Fixed 98

Globale Urbanisierungsprozesse und eine Zunahme der räumlichen Konzentration von exponierten Elementen wie Menschen, Gebäude, Infrastruktur und ökonomische Werte induzieren ein ungekanntes Risiko in erdbebengefährdeten Regionen. Wenn keine Abschwächung des Risikos erfolgt werden dramatische Folgen in der Zukunft erwartet. Diese umfassen eine beispiellose Anzahl an Todesopfer, enorme ökonomische und ökologische Verluste und Ausfälle bezüglich kritischer Infrastruktur und Versorgung etc. Um derartige Gefährdungen abzuschwächen sind detaillierte Informationen über seismisches Risiko notwendig. Die seismische Verwundbarkeit von Siedlungsarealen ist dabei als zentrale, konstituierende Komponente von seismischem Risiko zu berücksichtigen. In diesem Zusammenhang ist es von besonderem Interesse das Verhalten von Gebäudeinventaren unter einem bestimmten Erdbebeneinfluss abschätzen zu können. Das Hauptziel der Arbeit war es maßgeschneiderte Methoden zu entwickeln, die eine Bewertung der seismischen Vulnerabilität von Siedlungsräumen, basierend auf Fernerkundungsdaten, durchführbar machen. Es wurden Methoden aus dem Bereich des maschinellen Lernens adaptiert, um Verwundbarkeitsstufen von Gebäuden und homogenen Siedlungsstrukturen zu bestimmen. Hierfür wurden Merkmale aus Fernerkundungsdaten abgeleitet und mit in situ Informationen verknüpft. Wir verwenden verschiedene Ensembles von Fernerkundungssensoren, um die urbane Morphologie umfassend zu charakterisieren. Empirische Ergebnisse, die für die erdbebengefährdeten Städte Padang (Indonesien) und Istanbul (Türkei) generiert werden konnten, bestätigen die Durchführbarkeit der entwickelten Verfahren. Zukünftige Arbeiten können daran anknüpfen und beispielsweise empirische Erkenntnisse in weiteren Fallstudien anzweifeln, eine Verbesserung der Methodik vornehmen, Konzepte und Ansätze auf andere Sensorsysteme oder Datenquellen übertragen oder Daten und Methoden im Rahmen von holistischen Risikobewertungsstrategien anwenden. / Global urbanization processes and increasing spatial concentration of exposed elements such as people, buildings, infrastructure, and economic values in earthquake prone regions induce seismic risk at a uniquely high level. This situation, when left unmitigated, is expected to cause unprecedented death tolls, enormous economic and ecological losses, and critical infrastructure and service failures, etc., in the future. To mitigate those perils requires detailed knowledge about seismic risks. As an important constituent element of seismic risk, the seismic vulnerability of the built environment has to be assessed. In particular, it is crucial to know about the behavior of the building inventory under a certain level of ground shaking. The main goal of the thesis was to develop and evaluate tailored methods and procedures that allow for a viable seismic vulnerability assessment of the built environment with remote sensing data. In particular, methods from the machine learning domain were adapted to estimate vulnerability levels of buildings and homogeneous urban structures based on features derived from remote sensing and by incorporation of in situ knowledge. To this purpose we deploy ensembles of earth observation sensors to exhaustively characterize the urban morphology. Empirical results, obtained for the earthquake prone cities Padang (Indonesia) and Istanbul (Turkey), confirm the viability of the approaches. Overall, this thesis provides some promising results, which show that remote sensing has a high capability to contribute to a rapid screening assessment of the seismic vulnerability of buildings and urban structures. Further work can build upon these results and may challenge empirical findings in further case studies, enhance developed and applied methods, transfer concepts and approaches to other sensor systems and data sources, or apply data and methodologies within integrative and holistic risk assessment strategies.

Subduction zones are regions of intense earthquake activity up to great depth. Sources are located inside the subducting lithosphere and, as a consequence, seismic radiation from subduction zone earthquakes is strongly affected by the interior slab structure. The wave field of these intraslab events observed in the forearc region is profoundly influenced by a seismically slow layer atop the slab surface. This several kilometer thick low-velocity channel (wave guide) causes the entrapment of seismic energy producing strong guided wave phases that appear in P onsets in certain regions of the forearc. Observations at the Chile-Peru subduction zone presented here, as well as observations at several other circum-pacific subduction zones show such signals. Guided wave analysis contributes details of immense value regarding the processes near the slab surface, such as layering of subducted lithosphere, source locations of intraslab seismicity and most of all, range and manner of mineralogical phase transitions. <br><br> Seismological data stem from intermediate depth events (depth range 70 km - 300 km) recorded in northern Chile near 21 Grad S during the collaborative research initiative " Deformation Processes in the Andes" (SFB 267). A subset of stations - all located within a slab-parallel transect close to 69 Grad W - show low-frequency first arrivals (2 Hz), sometimes followed by a second high-frequency phase. <br><br> We employ 2-dimensional finite-difference simulations of complete P-SV wave propagation to explore the parameter space of subduction zone wave guides and explain the observations. Key processes underlying the guided wave propagation are studied: Two distinct mechanisms of decoupling of trapped energy from the wave guide are analyzed - a prerequisite to observe the phases at stations located at large distances from the wave guide (up to 100 km). Variations of guided wave effects perpendicular to the strike of the subduction zone are investigated, such as the influence of phases traveling in the fast slab. Further, the merits and limits of guided wave analysis are assessed. Frequency spectra of the guided wave onsets prove to be a robust quantity that captures guided wave characteristics at subduction zones including higher mode excitation. They facilitate the inference of wave guide structure and source positioning: The peak frequency of the guided wave fundamental mode is associated with a certain combination of layer width and velocity contrast. The excitation strength of the guided wave fundamental mode and higher modes is associated with source position and orientation relative to the low-velocity layer. <br><br> The guided wave signals at the Chile-Peru subduction zone are caused by energy that leaks from the subduction zone wave guide. On the one hand, the bend shape of the slab allows for leakage at a depth of 100 km. On the other, equalization of velocities between the wave guide and the host rocks causes further energy leakage at the contact zone between continental and oceanic crust (70 km depth). Guided waves bearing information on deep slab structure can therefore be recorded at specific regions in the forearc. These regions are determined based on slab geometry, and their locations coincide with the observations. <br><br> A number of strong constraints on the structure of the Chile-Peru slab are inferred: The deep wave guide for intraslab events is formed by a layer of 2 km average width that remains seismically slow (7 percent velocity reduction compared to surrounding mantle). This low-velocity layer at the top of the Chile-Peru slab is imaged from a depth of 100 km down to at least 160 km. Intermediate depth events causing the observed phases are located inside the layer or directly beneath it in the slab mantle. The layer is interpreted as partially eclogized lower oceanic crust persisting to depth beyond the volcanic arc. / Subduktionszonen sind bis in große Tiefen von intensiver Erdbebentätigkeit geprägt. Die Erdbebenquellen befinden sich in der subduzierten Lithosphäre (Slab), ihr Wellenfeld wird deshalb stark von der internen Slab-Struktur beeinflusst. Eine Schicht mit reduzierter seismischer Geschwindigkeit im oberen Bereich der Platte kann als Wellenleiter für diese Signale fungieren. In der nur wenige Kilometer dicken Schicht entstehen sogenannte geführte Wellen, die in Teilen des Forearc beobachtet werden. Diese Phasen bergen wertvolle Informationen über die Struktur nahe der Slab-Oberfläche, wie zum Beispiel Dicke der Schichtung, Herdlokationen und vor allem Tiefe und Art mineralogischer Umsetzungen. <br><br> Die Beobachtungen stammen von mitteltiefen Beben (70 km - 300 km) im Untersuchungsgebiet in Nord-Chile und wurden im Rahmen des Sonderforschungsbereich 267 " Deformationsprozesse in den Anden" aufgezeichnet. Stationen in einem Streifen um 69 Grad W, der sich parallel zum Streichen der Subduktionszone erstreckt, zeigen niederfrequente Ersteinsätze, denen teilweise höherfrequente Phasen folgen. <br><br> Mit Hilfe eines 2-dimensionalen Finite-Differenzen-Algorithmus werden die P-SV Wellenausbreitung simuliert, und die Beobachtungen erklärt. Zentrale Fragestellungen zu Wellenleitern in Subduktionszonen werden untersucht: Es werden zwei Mechanismen, die das Auskoppeln seismischer Energie aus dem Wellenleiter ermöglichen beschrieben - eine Grundvoraussetzung für das Auftreten von geführten Wellen in großen Entfernungen vom Wellenleiter (bis zu 100 km). Des weiteren werden Stärken und Grenzen der Analyse von geführten Wellen erörtert. <br><br> Die Spektren der geführten Wellenzüge erweisen sich als robuste Messgröße, um die Charakteristika des Wellenleiters zu bestimmt. Struktur des Wellenleiters und Quellpositionen können so für festgelegte Quell-Empfänger-Geometrien abgeleitet werden. Die Peak-Frequenz der Grundmode wird durch eine Kombination aus Dicke der Schicht und Geschwindigkeitskontrast bestimmt. Die Stärke der Anregung der Grundmode und höherer Moden lässt auf die Lage und Orientierung der Erdbebenquelle relativ zur Schicht schließen. Geschwindigkeitskontrast, Schichtdicke und Quellposition sind von herausragender Bedeutung, um mineralogische Interpretationen des Wellenleiters zu überprüfen. <br><br> Aufbauend auf die Simulationen werden die Beobachtungen interpretiert und Auskunft über die Struktur der Chile-Peru Subduktionszone erhalten: Eine dünne Schicht an der Slab-Oberfläche (durchschnittlich 2 km dick) trägt geringere seismische Geschwindigkeiten als der umgebende Mantel und fungiert als Wellenleiter für intra-platten Ereignisse in Tiefen von 100 bis mindestens 160 km. Ereignisse, die geführte Wellen hervorrufen, liegen in dieser Schicht oder direkt darunter im subduzierten Mantel. Um zu den Stationen in der Forearc-Region zu gelangen, entkoppelt ein Teil der geführten Wellen in einer Tiefe von circa 100 km aus der Niedergeschwindigkeitsschicht. Die Krümmung des Slab erlaubt das Austreten der Wellen und nimmt auch Einfluss auf die Pulsformen. <br><br> Der Wellenleiter in der Chile-Peru Subduktionszone ergibt sich als unregelmäßige Schicht mit reduzierter seismischer Geschwindigkeit, in der geführte Wellen entstehen, in unterschiedlichen Tiefen wieder austreten, und an die freie Oberfläche gelangen. Die Beobachtungsgebiete befinden sich im Forearc und werden durch die Geometrie und Struktur der subduzierten Platte festgelegt. <br><br> Die nur wenige Kilometer dicke, seismisch langsame Schicht an der Oberfläche des Chile-Peru Slab legt nahe, dass die Unterkruste der subduzierten Platte bis in große Tiefen besteht und nicht vollständig eklogitisiert ist. Abgeleitete Schichtdicke, Geschwindigkeitskontrast

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. / 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.

Das Bodenverhalten unter dynamischer Beanspruchung ist sehr komplex, wird jedoch in der Praxis häufig mit Hilfe von vereinfachten Modellen abgebildet. Die Gültigkeit solcher Modelle ist jedoch aufgrund des spannungs- und dehnungsabhängigen Bodenverhaltens sehr begrenzt. Alternativ dazu bieten sich dynamische numerische Berechnungen mit fortgeschrittenen Stofmodellen, die das Bodenverhalten in einem großen Dehnungs- und Spannungsbereich realitätsnah repräsentieren können. In dieser Arbeit wurde untersucht, inwieweit sich das komplexe Bodenverhalten unter dynamischen Einwirkungen mit Hilfe der Hypoplastizität abbilden lässt. Dabei wurde die entscheidende Rolle der Parameterermittlung veranschaulicht und zusätzlich ein angemessener Vorgang zur Bodenparameterbestimmung beschrieben. Zunächst wurde das Verhalten einer trockenen Sandschicht infolge von Erdbebenbeanspruchungen numerisch untersucht. Die Ergebnisse der Berechnungen zeigen, dass die Beschleunigungsamplifikation in der Nähe zur Bodenoberfläche von der Frequenz und der Amplitude der Grundbeschleunigung abhängt. Weiterhin nimmt die berechnete Eigenfrequenz und die entsprechende Amplifikation mit zunehmender Beschleunigungsamplitude ab. Des Weiteren wurde ein Zentrifugenversuch an einem im Sand eingebeteten Tunnel unter Erdbebeneinwirkungen nachgerechnet. Die berechneten Ergebnisse zeigen eine ausreichende Übereinstimmung mit dem Experiment. Mit der numerischen Nachrechnung wurde auch eine Abhängigkeit zwischen den Änderungen der Biegemomente in der Tunnelschale und der Oberflächensetzung im umliegenden Boden festgestellt. Die Standsicherheit von Böschungen unter Erdbebenbeanspruchungen stellt wegen des komplexen Bodenverhaltens eine weitere Herausforderung für die Berechnungen dar. Zunächst wurde überprüft, inwieweit sich das Böschungsverhalten mit der in der Praxis häufig eingesetzten pseudo-statischen Methode abbilden lässt. Hierfür wurde für eine in der Zentrifuge untersuchte Modellböschung die pseudostatische Analyse durchgeführt. Die im Zentrifugenversuch aufgetretenen oberflächennahen Gleitfläche lässt sich durch die pseudo-statische Methode nicht prognostizieren. Für eine oberflächennahe Gleitfläche wurde hingegen ein sehr hoher Standsicherheitsfaktor ermittelt. Mit einer numerischen Nachrechnung mit einem hypoplastischen Stoffmodell mit Betrachtung der intergranularen Dehnungen konnte das Verhalten der Modellböschung qualitativ und quantitativ sehr gut abgebildet werden. Somit wurden sowohl die oberflächennahe Gleitfläche als auch die Vertikal- und Horizontalverschiebungen realitätsnah wiedergegeben. In dieser Arbeit wurde des Weiteren ein Vorgang als Kombination zwischen den dynamischen numerischen Berechnungen und der pseudo-statischen Methode zur Bewertung der Standsicherheit von Böschungen unter dynamischer Einwirkung vorgeschlagen. Damit ließ sich ebenso ein realitätsnäher Stansicherheitsfaktor ermitteln. Da die Anwendung der pseudo-statischen Methode bei den Böschungen aus wassergesättigten kohäsionslosen Böden problematisch ist, lassen sich solche Böschungen entweder mit Zentrifugenmodellen oder numerisch mit fortgeschrittenen Stoffmodellen untersuchen. In dieser Arbeit wurden Nachrechnungen von Zentrifugenversuchen durchgeführt. Es handelt sich um einen Erddamm aus einem wassergesättigten, dicht gelagerten Nevada Sand unter Erdbebeneinwirkung. Mit der numerischen Berechnung wurde das Dammverhalten qualitativ und quantitativ sehr gut abgebildet. Sowohl die Dammverschiebungen als auch der Aufbau des Porenwasserdrucks zeigen eine sehr gute Übereinstimmung mit den Messungen. Weiterhin wurden mit den gleichen Bodenparametern zwei weitere Zentrifugenversuche unter Erdbebeneinwirkung nachgerechnet. Beide Modellversuche wurden mit einem locker gelagerten, wassergesättigten Nevada Sand durchgeführt. Bei einem Versuch wurde ein Erddamm und bei dem anderen eine Sandschicht untersucht. In den numerischen Nachrechnungen ließen sich sowohl die Verschiebungen als auch die Porenwasserdrücke in beiden Randwertproblemen realistisch abbilden. Weiterhin wurde die Wirkung von Schottersäulen zur Verhinderung der Bodenverflüssigung numerisch untersucht. Zunächst wurden die Dränage- und die Aussteifungswirkung der Schottersäulen unabhängig voneinander betrachtet. Die Dränagewirkung ist vernachlässigbar, da sich während eines Erdbebens der Porenwasserdruck sehr schnell aufbaut. Wegen der hohen Steifigkeit der Schottersäulen wird zwar weniger Porenwasserdruck in den Boden aufgebaut. Die effektive Spannung nimmt jedoch trotzdem unverhindert ab. Dies lässt sich damit begründen, dass die hohe Säulensteifigkeit zu einer Spannungsumlagerung in Richtung Säulen führt und ein Siloeffekt entsteht. Somit wird der Boden zum Teil von den Säulen getragen und die totale Spannung im Boden nimmt ab. In der 3D-Berechnungen ist dieser Siloeffekt deutlich geringer als in den 2D-Berechnungen. Nichtsdestotrotz zeigen sowohl die 2D- als auch die 3D-Berechnungen, dass die Säulensteifigkeit eine nur mäßige Wirkung zur Verhinderung der Bodenverflüssigung aufweist. In weiteren 3D-Berechnungen wurde der Einfluss der Säulenherstellung untersucht. Hierfür wurden Berechnungen mit erhöhter Bodendichte und Seitenspannung durchgeführt. Sowohl die Verdichtung als auch die Erhöhung der Seitenspannung verlangsamen den Porenwasserdruckaufbau bzw. die Abnahme der effektiven Spannung. Der Einfluss der Bodenverdichtung ist jedoch wesentlich höher. Weiterhin weist die Wirkung der Schottersäulen eine Abhängigkeit von der dynamischen Belastung auf. Die Bodenverflüssigung infolge eines kleinen Erdbebens wird verhindert, während sich die Verflüssigung infolge eines stärkeren Erdbebens nur um wenige Sekunden verzögert. / The soil behavior under dynamic loading is very complex. However, in daily use it is often illustrated by means of simplified models. The validity of these models is very limited due to the stress and strain-dependent soil behavior. Alternatively, dynamic numerical calculations can be performed with advanced constitutive models which can represent soil behavior in a wide range of strain and stress. In this work it was investigated, to which extent the complex soil behavior can be reproduced using hypoplasticity.Furthermore,the important role of parameter determination was illustrated. In addition, an appropriate procedure for determining soil parameters was described. First, the behavior of a dry sand layer under earthquake load was investigated numerically. The results of the calculations show that the acceleration amplification near the ground surface depends on the frequency and the amplitude of the basic acceleration. Furthermore, the calculated natural frequency and the corresponding amplification decrease with increasing acceleration amplitude. In addition, a centrifuge test on a tunnel embedded in sand under earthquake effects was numerically calculated. The calculated results show a satisfactory agreement with the experiment. The numerical calculation also revealed a dependency between the changes in the bending moments in the tunnel lining and the surface settlement of the surrounding soil. Due to the complex soil behavior, the stability of slopes under earthquake loads poses a further challenge for the calculations. Firstly, it was examined, to which extent slope behavior can be represented with the frequently used pseudo-static method. For this purpose the pseudo-static analysis was carried out for a model earth dam examined in the centrifuge. The pseudo-static method predicts a deep seated sliding surface in contrast to the shallow sliding surface in the centrifuge test. However, for a shallow sliding surface, a very high stability safety factor was determined. With a numerical calculation using a hypoplastic material model considering the intergranular strains, the behavior of the earth dam could be reproduced qualitatively and quantitatively very well. Thus, the shallow sliding surface as well as the vertical and horizontal displacements were reproduced realistically. In this thesis, a combination of the dynamic numerical calculation and the pseudo-static method for assessing the stability of slopes under dynamic influence was proposed. So, a realistic stability safety factor can be determined. The application of the pseudo-static method is problematic in case of slopes in saturated non-cohesive soil. These slopes can either be investigated with centrifuge models or numerically with advanced material models. In this work, numerical recalculations of centrifuge tests were carried out. It is an earth dam from a saturated Nevada sand under an earthquake effect. With the numerical calculation the dam behavior was reproduced qualitatively and quantitatively in a satisfactory manner. Both the dam displacements as well as the build-up of pore water pressure show a very good agreement with the measurements. Two further centrifuge tests were also carried out using the same soil parameters. Both model tests were conducted with a loose saturated Nevada sand. One test was carried out on an earth dam and the other on a sand layer. With the numerical calculations, both displacements and pore water pressures were reproduced realistically in both boundary value problems. In addition the effect of stone columns to prevent soil liquefaction was studied numerically. First, the drainage and stiffening effects of stone columns were examined separately. The drainage effect has no significant influence because of the very rapid build-up of pore water pressure during the earthquake. Due to the high stiffness of the stone columns, less pore water pressure builds up in the soil. However, the effective stress continues to decrease unhindered. The high stiffness of the columns leads to a stress redistribution in the direction of the columns and a silo effect arises. In 3D calculations, the silo effect is significantly lower than in 2D calculations. The 2D and 3D calculations show that the column stiffness has a moderate effect to prevent soil liquefaction. In further 3D calculations, the influence of column installation was investigated. Calculations with increased soil density and lateral stress were carried out for this purpose. Both the compaction and the increase of the lateral stress slow down the build-up of pore water pressure and the decrease in effective stress. However, the impact of soil compaction is much higher. Furthermore, the effect of stone columns depends on the dynamic load. The soil liquefaction due to a small earthquake is prevented, while liquefaction due to a stronger earthquake is delayed only by a few seconds.