TY - JOUR A1 - Flóvenz, Ólafur G. A1 - Wang, Rongjiang A1 - Hersir, Gylfi Páll A1 - Dahm, Torsten A1 - Hainzl, Sebastian A1 - Vassileva, Magdalena A1 - Drouin, Vincent A1 - Heimann, Sebastian A1 - Isken, Marius Paul A1 - Gudnason, Egill Á. A1 - Ágústsson, Kristján A1 - Ágústsdóttir, Thorbjörg A1 - Horálek, Josef A1 - Motagh, Mahdi A1 - Walter, Thomas R. A1 - Rivalta, Eleonora A1 - Jousset, Philippe A1 - Krawczyk, Charlotte M. A1 - Milkereit, Claus T1 - Cyclical geothermal unrest as a precursor to Iceland's 2021 Fagradalsfjall eruption JF - Nature geoscience N2 - Understanding and constraining the source of geodetic deformation in volcanic areas is an important component of hazard assessment. Here, we analyse deformation and seismicity for one year before the March 2021 Fagradalsfjall eruption in Iceland. We generate a high-resolution catalogue of 39,500 earthquakes using optical cable recordings and develop a poroelastic model to describe three pre-eruptional uplift and subsidence cycles at the Svartsengi geothermal field, 8 km west of the eruption site. We find the observed deformation is best explained by cyclic intrusions into a permeable aquifer by a fluid injected at 4 km depth below the geothermal field, with a total volume of 0.11 ± 0.05 km3 and a density of 850 ± 350 kg m–3. We therefore suggest that ingression of magmatic CO2 can explain the geodetic, gravity and seismic data, although some contribution of magma cannot be excluded. Y1 - 2022 U6 - https://doi.org/10.1038/s41561-022-00930-5 SN - 1752-0894 SN - 1752-0908 VL - 15 IS - 5 SP - 397 EP - 404 PB - Nature Research CY - Berlin ER - TY - JOUR A1 - Dahm, Torsten A1 - Heimann, Sebastian A1 - Metz, Malte A1 - Isken, Marius Paul T1 - A self-similar dynamic rupture model based on the simplified wave-rupture analogy JF - Geophysical journal international / the Royal Astronomical Society, the Deutsche Geophysikalische Gesellschaft and the European Geophysical Society N2 - 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. KW - Earthquake dynamics KW - Earthquake ground motions KW - Earthquake hazards KW - Earthquake source observations Y1 - 2021 U6 - https://doi.org/10.1093/gji/ggab045 SN - 0956-540X SN - 1365-246X VL - 225 IS - 3 SP - 1586 EP - 1604 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Jamalreyhani, Mohammadreza A1 - Rezapour, Mehdi A1 - Cesca, Simone A1 - Dahm, Torsten A1 - Heimann, Sebastian A1 - Sudhaus, Henriette A1 - Isken, Marius Paul T1 - Insight into the 2017-2019 Lurestan arc seismic sequence (Zagros, Iran); complex earthquake interaction in the basement and sediments JF - Geophysical journal international N2 - 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. KW - Joint Inversion KW - Waveform inversion KW - Earthquake source observations KW - Seismicity and tectonics Y1 - 2022 U6 - https://doi.org/10.1093/gji/ggac057 SN - 0956-540X SN - 1365-246X VL - 230 IS - 1 SP - 114 EP - 130 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Pilz, Marco A1 - Isken, Marius Paul A1 - Fleming, Kevin A1 - Orunbaev, Sagynbek A1 - Moldobekov, Bolot T1 - Long- and short-term monitoring of a dam in response to seasonal changes and ground motion loading BT - the test case of the Kurpsai Dam, Western Kyrgyz Republic JF - Pure and applied geophysics : PAGEOPH ; continuation of Geofisica pura e applicata N2 - An experimental multi-parameter structural monitoring system has been installed on the Kurpsai dam, western Kyrgyz Republic. This system consists of equipment for seismic and strain measurements for making longer- (days, weeks, months) and shorter- (minutes, hours) term observations, dealing with, for example seasonal (longer) effects or the response of the dam to ground motion from noise or seismic events. Fibre-optic strain sensors allow the seasonal and daily opening and closing of the spaces between the dam's segments to be tracked. For the seismic data, both amplitude (in terms of using differences in amplitudes in the Fourier spectra for mapping the modes of vibration of the dam) and their time-frequency distribution for a set of small to moderate seismic events are investigated and the corresponding phase variabilities (in terms of lagged coherency) are evaluated. Even for moderate levels of seismic-induced ground motion, some influence on the structural response can be detected, which then sees the dam quickly return to its original state. A seasonal component was identified in the strain measurements, while levels of noise arising from the operation of the dam's generators and associated water flow have been provisionally identified. KW - Structural health monitoring KW - Dam engineering KW - Operational and environmental effects KW - Strong-motion KW - Strain KW - Elastic response KW - Kurpsai dam Y1 - 2021 U6 - https://doi.org/10.1007/s00024-021-02861-5 SN - 0033-4553 SN - 1420-9136 VL - 178 IS - 10 SP - 4001 EP - 4020 PB - Birkhäuser CY - Basel ER - TY - JOUR A1 - Isken, Marius Paul A1 - Vasyura-Bathke, Hannes A1 - Dahm, Torsten A1 - Heimann, Sebastian T1 - De-noising distributed acoustic sensing data using an adaptive frequency-wavenumber filter JF - Geophysical journal international N2 - Data recorded by distributed acoustic sensing (DAS) along an optical fibre sample the spatial and temporal properties of seismic wavefields at high spatial density. Often leading to massive amount of data when collected for seismic monitoring along many kilometre long cables. The spatially coherent signals from weak seismic arrivals within the data are often obscured by incoherent noise. We present a flexible and computationally efficient filtering technique, which makes use of the dense spatial and temporal sampling of the data and that can handle the large amount of data. The presented adaptive frequency-wavenumber filter suppresses the incoherent seismic noise while amplifying the coherent wavefield. We analyse the response of the filter in time and spectral domain, and we demonstrate its performance on a noisy data set that was recorded in a vertical borehole observatory showing active and passive seismic phase arrivals. Lastly, we present a performant open-source software implementation enabling real-time filtering of large DAS data sets. KW - Fourier analysis KW - Image processing KW - Time-series analysis KW - Seismic noise KW - Distributed acoustic sensing Y1 - 2022 U6 - https://doi.org/10.1093/gji/ggac229 SN - 0956-540X SN - 1365-246X VL - 231 IS - 2 SP - 944 EP - 949 PB - Oxford University Press CY - Oxford ER - TY - JOUR A1 - Heimann, Sebastian A1 - Vasyura-Bathke, Hannes A1 - Sudhaus, Henriette A1 - Isken, Marius Paul A1 - Kriegerowski, Marius A1 - Steinberg, Andreas A1 - Dahm, Torsten T1 - A Python framework for efficient use of pre-computed Green's functions in seismological and other physical forward and inverse source problems JF - Solid earth N2 - The computation of such synthetic GFs is computationally and operationally demanding. As a consequence, the onthe-fly recalculation of synthetic GFs in each iteration of an optimisation is time-consuming and impractical. Therefore, the pre-calculation and efficient storage of synthetic GFs on a dense grid of source to receiver combinations enables the efficient lookup and utilisation of GFs in time-critical scenarios. We present a Python-based framework and toolkit - Pyrocko-GF - that enables the pre-calculation of synthetic GF stores, which are independent of their numerical calculation method and GF transfer function. The framework aids in the creation of such GF stores by interfacing a suite of established numerical forward modelling codes in seismology (computational back ends). So far, interfaces to back ends for layered Earth model cases have been provided; however, the architecture of Pyrocko-GF is designed to cover back ends for other geometries (e.g. full 3-D heterogeneous media) and other physical quantities (e.g. gravity, pressure, tilt). Therefore, Pyrocko-GF defines an extensible GF storage format suitable for a wide range of GF types, especially handling elasticity and wave propagation problems. The framework assists with visualisations, quality control, and the exchange of GF stores, which is supported through an online platform that provides many pre-calculated GF stores for local, regional, and global studies. The Pyrocko-GF toolkit comes with a well-documented application programming interface (API) for the Python programming language to efficiently facilitate forward modelling of geophysical processes, e.g. synthetic waveforms or static displacements for a wide range of source models. Y1 - 2019 U6 - https://doi.org/10.5194/se-10-1921-2019 SN - 1869-9510 SN - 1869-9529 VL - 10 IS - 6 SP - 1921 EP - 1935 PB - Copernicus CY - Göttingen ER -