Refine
Year of publication
Document Type
- Doctoral Thesis (11)
- Article (3)
- Other (3)
- Master's Thesis (2)
- Bachelor Thesis (1)
- Habilitation Thesis (1)
Keywords
- Seismologie (21) (remove)
Institute
- Institut für Geowissenschaften (20)
- Extern (2)
- Historisches Institut (1)
Seismological and seismotectonic analysis of the northwestern Argentine Central Andean foreland
(2020)
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.
A comprehensive study on seismic hazard and earthquake triggering is crucial for effective mitigation of earthquake risks. The destructive nature of earthquakes motivates researchers to work on forecasting despite the apparent randomness of the earthquake occurrences. Understanding their underlying mechanisms and patterns is vital, given their potential for widespread devastation and loss of life. This thesis combines methodologies, including Coulomb stress calculations and aftershock analysis, to shed light on earthquake complexities, ultimately enhancing seismic hazard assessment.
The Coulomb failure stress (CFS) criterion is widely used to predict the spatial distributions of aftershocks following large earthquakes. However, uncertainties associated with CFS calculations arise from non-unique slip inversions and unknown fault networks, particularly due to the choice of the assumed aftershocks (receiver) mechanisms. Recent studies have proposed alternative stress quantities and deep neural network approaches as superior to CFS with predefined receiver mechanisms. To challenge these propositions, I utilized 289 slip inversions from the SRCMOD database to calculate more realistic CFS values for a layered-half space and variable receiver mechanisms. The analysis also investigates the impact of magnitude cutoff, grid size variation, and aftershock duration on the ranking of stress metrics using receiver operating characteristic (ROC) analysis. Results reveal the performance of stress metrics significantly improves after accounting for receiver variability and for larger aftershocks and shorter time periods, without altering the relative ranking of the different stress metrics.
To corroborate Coulomb stress calculations with the findings of earthquake source studies in more detail, I studied the source properties of the 2005 Kashmir earthquake and its aftershocks, aiming to unravel the seismotectonics of the NW Himalayan syntaxis. I simultaneously relocated the mainshock and its largest aftershocks using phase data, followed by a comprehensive analysis of Coulomb stress changes on the aftershock planes. By computing the Coulomb failure stress changes on the aftershock faults, I found that all large aftershocks lie in regions of positive stress change, indicating triggering by either co-seismic or post-seismic slip on the mainshock fault.
Finally, I investigated the relationship between mainshock-induced stress changes and associated seismicity parameters, in particular those of the frequency-magnitude (Gutenberg-Richter) distribution and the temporal aftershock decay (Omori-Utsu law). For that purpose, I used my global data set of 127 mainshock-aftershock sequences with the calculated Coulomb Stress (ΔCFS) and the alternative receiver-independent stress metrics in the vicinity of the mainshocks and analyzed the aftershocks properties depend on the stress values. Surprisingly, the results show a clear positive correlation between the Gutenberg-Richter b-value and induced stress, contrary to expectations from laboratory experiments. This observation highlights the significance of structural heterogeneity and strength variations in seismicity patterns. Furthermore, the study demonstrates that aftershock productivity increases nonlinearly with stress, while the Omori-Utsu parameters c and p systematically decrease with increasing stress changes. These partly unexpected findings have significant implications for future estimations of aftershock hazard.
The findings in this thesis provides valuable insights into earthquake triggering mechanisms by examining the relationship between stress changes and aftershock occurrence. The results contribute to improved understanding of earthquake behavior and can aid in the development of more accurate probabilistic-seismic hazard forecasts and risk reduction strategies.
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.
"Terrae motus factus est magnus". In diesen und ähnlichen Worten erinnern mittelalterliche Geschichtsschreiber stets an das verspürte Eintreten von Erdbeben. Für die ereignisgeschichtliche Rekonstruktion der historischen Seismizität besitzt das Verständnis, dieser seit dem Frühmittelalter zunehmend standardisiert gebrauchten Narrativen, einen hohen Wert. Daher ist es wichtig, mit den bislang nahezu unerkannt geblieben Intentionen, Vorstellungsstrukturen und Argumentationsstrategien früh- und hochmittelalterlicher Geschichtsschreiber bekannt zu werden. Ausgehend von den antiken Ursprüngen ermittelt diese Arbeit die Bandbreite einer auf "terrae motus" aufbauenden, spezifisch mittelalterlichen Traditionsbildung und setzt sie in den Kontext zum Wissens- und Erfahrungshorizont früh- und hochmittelalterlicher Gelehrter. Erdbeben besaßen ein außerordentliches hermeneutisches Potential für das mittelalterliche Weltverständnis. Somit sind mittelalterliche Erdbebenbeschreibungen hinsichtlich ihrer deskriptiven Qualität und argumentativen Wertigkeit verschieden. Die Historiographie- und Ideengeschichte sowie die seismologische Parametrisierung von mittelalterlichen Erdbeben wird von diesem Wissen gleichermaßen profitieren.
Tensile source components of swarm events in West Bohemia in 2000 by considering seismic anisotropy
(2006)
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.
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.
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.
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.
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.
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.