@article{DahmHeimannMetzetal.2021, author = {Dahm, Torsten and Heimann, Sebastian and Metz, Malte and Isken, Marius Paul}, title = {A self-similar dynamic rupture model based on the simplified wave-rupture analogy}, series = {Geophysical journal international / the Royal Astronomical Society, the Deutsche Geophysikalische Gesellschaft and the European Geophysical Society}, volume = {225}, journal = {Geophysical journal international / the Royal Astronomical Society, the Deutsche Geophysikalische Gesellschaft and the European Geophysical Society}, number = {3}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggab045}, pages = {1586 -- 1604}, year = {2021}, abstract = {The investigation of stresses, faults, structure and seismic hazards requires a good understanding and mapping of earthquake rupture and slip. Constraining the finite source of earthquakes from seismic and geodetic waveforms is challenging because the directional effects of the rupture itself are small and dynamic numerical solutions often include a large number of free parameters. The computational effort is large and therefore difficult to use in an exploratory forward modelling or inversion approach. Here, we use a simplified self-similar fracture model with only a few parameters, where the propagation of the fracture front is decoupled from the calculation of the slip. The approximative method is flexible and computationally efficient. We discuss the strengths and limitations of the model with real-case examples of well-studied earthquakes. These include the M-w 8.3 2015 Illapel, Chile, megathrust earthquake at the plate interface of a subduction zone and examples of continental intraplate strike-slip earthquakes like the M-w 7.1 2016 Kumamoto, Japan, multisegment variable slip event or the M-w 7.5 2018 Palu, Indonesia, supershear earthquake. Despite the simplicity of the model, a large number of observational features ranging from different rupture-front isochrones and slip distributions to directional waveform effects or high slip patches are easy to model. The temporal evolution of slip rate and rise time are derived from the incremental growth of the rupture and the stress drop without imposing other constraints. The new model is fast and implemented in the open-source Python seismology toolbox Pyrocko, ready to study the physics of rupture and to be used in finite source inversions.}, language = {en} } @article{JamalreyhaniRezapourCescaetal.2022, author = {Jamalreyhani, Mohammadreza and Rezapour, Mehdi and Cesca, Simone and Dahm, Torsten and Heimann, Sebastian and Sudhaus, Henriette and Isken, Marius Paul}, title = {Insight into the 2017-2019 Lurestan arc seismic sequence (Zagros, Iran); complex earthquake interaction in the basement and sediments}, series = {Geophysical journal international}, volume = {230}, journal = {Geophysical journal international}, number = {1}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggac057}, pages = {114 -- 130}, year = {2022}, abstract = {Despite its high-seismogenic potential, the details of the seismogenic processes of Zagros Simply Folded Belt (SFB) remains debated. Three large earthquakes (M-w 7.3, 5.9 and 6.3) struck in the Lurestan arc of the Zagros SFB in 2017 and 2018. The sequence was recorded by seismic stations at regional, and teleseismic distances. Coseismic surface displacements, measured by Sentinel-1A/B satellites, provide additional data and a unique opportunity to study these earthquakes in detail. Here, we complement previous studies of the coseismic slip distribution of the 12 November 2017 M-w 7.3 Ezgeleh earthquake by a detailed analysis of its aftershocks, and we analysed the rupture process of the two interrelated earthquakes (25 August 2018 M-w 5.9 Tazehabad and the 25 November 2018 M-w 6.3 Sarpol-e Zahab earthquakes). We model the surface displacements obtained from Interferometric Synthetic Aperture Radar (InSAR) measurements and seismic records. We conduct non-linear probabilistic optimizations based on joint InSAR and seismic data to obtain finite-fault rupture of these earthquakes. The Lurestan arc earthquakes were followed by a sustained aftershock activity, with 133 aftershocks exceeding M-n 4.0 until 30 December 2019. We rely on the permanent seismic networks of Iran and Iraq to relocate similar to 700 M-n 3 + events and estimate moment tensor solutions for 85 aftershocks down to M-w 4.0. The 2017 Ezgeleh earthquake has been considered to activate a low-angle (similar to 17 degrees) dextral-thrust fault at the depth of 10-20 km. However, most of its aftershocks have shallow centroid depths (8-12 km). The joint interpretation of finite source models, moment tensor and hypocentral location indicate that the 2018 Tazehabad and Sarpol-e Zahab earthquakes ruptured different strike-slip structures, providing evidence for the activation of the sinistral and dextral strike-slip faults, respectively. The deformation in the Lurestan arc is seismically accommodated by a complex fault system involving both thrust and strike-slip faults. Knowledge about the deformation characteristics is important for the understanding of crustal shortening, faulting and hazard and risk assessment in this region.}, language = {en} } @article{CescaGrigoliHeimannetal.2016, author = {Cesca, Simone and Grigoli, Francesco and Heimann, Sebastian and Dahm, Torsten and Kriegerowski, Marius and Sobiesiak, M. and Tassara, C. and Olcay, M.}, title = {The M-w 8.1 2014 Iquique, Chile, seismic sequence: a tale of foreshocks and aftershocks}, series = {Geophysical journal international}, volume = {204}, journal = {Geophysical journal international}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggv544}, pages = {1766 -- 1780}, year = {2016}, abstract = {The 2014 April 1, M-w 8.1 Iquique (Chile) earthquake struck in the Northern Chile seismic gap. With a rupture length of less than 200 km, it left unbroken large segments of the former gap. Early studies were able to model the main rupture features but results are ambiguous with respect to the role of aseismic slip and left open questions on the remaining hazard at the Northern Chile gap. A striking observation of the 2014 earthquake has been its extensive preparation phase, with more than 1300 events with magnitude above M-L 3, occurring during the 15 months preceding the main shock. Increasing seismicity rates and observed peak magnitudes accompanied the last three weeks before the main shock. Thanks to the large data sets of regional recordings, we assess the precursor activity, compare foreshocks and aftershocks and model rupture preparation and rupture effects. To tackle inversion challenges for moderate events with an asymmetric network geometry, we use full waveforms techniques to locate events, map the seismicity rate and derive source parameters, obtaining moment tensors for more than 300 events (magnitudes M-w 4.0-8.1) in the period 2013 January 1-2014 April 30. This unique data set of fore- and aftershocks is investigated to distinguish rupture process models and models of strain and stress rotation during an earthquake. Results indicate that the spatial distributions of foreshocks delineated the shallower part of the rupture areas of the main shock and its largest aftershock, well matching the spatial extension of the aftershocks cloud. Most moment tensors correspond to almost pure double couple thrust mechanisms, consistent with the slab orientation. Whereas no significant differences are observed among thrust mechanisms in different areas, nor among thrust foreshocks and aftershocks, the early aftershock sequence is characterized by the presence of normal fault mechanisms, striking parallel to the trench but dipping westward. These events likely occurred in the shallow wedge structure close to the slab interface and are consequence of the increased extensional stress in this region after the largest events. The overall stress inversion result suggests a minor stress rotation after the main shock, but a significant release of the deviatoric stress. The temporal change in the distribution of focal mechanisms can also be explained in terms of the spatial heterogeneity of the stress field: under such interpretation, the potential of a large megathrust earthquake breaking a larger segment offshore Northern Chile remains high.}, language = {en} } @article{KulikovaSchurrKruegeretal.2016, author = {Kulikova, Galina and Schurr, Bernd and Kr{\"u}ger, Frank and Brzoska, Elisabeth and Heimann, Sebastian}, title = {Source parameters of the Sarez-Pamir earthquake of 1911 February 18}, series = {Geophysical journal international}, volume = {205}, journal = {Geophysical journal international}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggw069}, pages = {1086 -- 1098}, year = {2016}, abstract = {The Ms \&\#8764; 7.7 Sarez-Pamir earthquake of 1911 February 18 is the largest instrumentally recorded earthquake in the Pamir region. It triggered one of the largest landslides of the past century, building a giant natural dam and forming Lake Sarez. As for many strong earthquakes from that time, information about source parameters of the Sarez-Pamir earthquake is limited due to the sparse observations. Here, we present the analysis of analogue seismic records of the Sarez-Pamir earthquake. We have collected, scanned and digitized 26 seismic records from 13 stations worldwide to relocate the epicentre and determine the event's depth (\&\#8764;26 km) and magnitude (mB7.3 and Ms7.7). The unusually good quality of the digitized waveforms allowed their modelling, revealing an NE-striking sinistral strike-slip focal mechanism in accordance with regional tectonics. The shallow depth and magnitude (Mw7.3) of the earthquake were confirmed. Additionally, we investigated the possible contribution of the landslide to the waveforms and present an alternative source model assuming the landslide and earthquake occurred in close sequence.}, language = {en} } @article{CescaGrigoliHeimannetal.2014, author = {Cesca, Simone and Grigoli, Francesco and Heimann, Sebastian and Gonzalez, Alvaro and Buforn, Elisa and Maghsoudi, Samira and Blanch, Estefania and Dahm, Torsten}, title = {The 2013 September-October seismic sequence offshore Spain: a case of seismicity triggered by gas injection?}, series = {Geophysical journal international}, volume = {198}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggu172}, pages = {941 -- 953}, year = {2014}, abstract = {A spatially localized seismic sequence originated few tens of kilometres offshore the Mediterranean coast of Spain, close to the Ebro river delta, starting on 2013 September 5, and lasting at least until 2013 October. The sequence culminated in a maximal moment magnitude M-w 4.3 earthquake, on 2013 October 1. The most relevant seismogenic feature in the area is the Fosa de Amposta fault system, which includes different strands mapped at different distances to the coast, with a general NE-SW orientation, roughly parallel to the coastline. However, no significant known historical seismicity has involved this fault system in the past. The epicentral region is also located near the offshore platform of the Castor project, where gas is conducted through a pipeline from mainland and where it was recently injected in a depleted oil reservoir, at about 2 km depth. We analyse the temporal evolution of the seismic sequence and use full waveform techniques to derive absolute and relative locations, estimate depths and focal mechanisms for the largest events in the sequence (with magnitude mbLg larger than 3), and compare them to a previous event (2012 April 8, mbLg 3.3) taking place in the same region prior to the gas injection. Moment tensor inversion results show that the overall seismicity in this sequence is characterized by oblique mechanisms with a normal fault component, with a 30A degrees low-dip angle plane oriented NNE-SSW and a subvertical plane oriented NW-SE. The combined analysis of hypocentral location and focal mechanisms could indicate that the seismic sequence corresponds to rupture processes along shallow low-dip surfaces, which could have been triggered by the gas injection in the reservoir, and excludes the activation of the Amposta fault, as its known orientation is inconsistent with focal mechanism results. An alternative scenario includes the iterated triggering of a system of steep faults oriented NW-SE, which were identified by prior marine seismics investigations.}, language = {en} }