TY - JOUR A1 - Nooshiri, Nima A1 - Saul, Joachim A1 - Heimann, Sebastian A1 - Tilmann, Frederik A1 - Dahm, Torsten T1 - Revision of earthquake hypocentre locations in global bulletin data sets using source-specific station terms JF - Geophysical journal international N2 - Global earthquake locations are often associated with very large systematic travel-time residuals even for clear arrivals, especially for regional and near-regional stations in subduction zones because of their strongly heterogeneous velocity structure. Travel-time corrections can drastically reduce travel-time residuals at regional stations and, in consequence, improve the relative location accuracy. We have extended the shrinking-box source-specific station terms technique to regional and teleseismic distances and adopted the algorithm for probabilistic, nonlinear, global-search location. We evaluated the potential of the method to compute precise relative hypocentre locations on a global scale. The method has been applied to two specific test regions using existing P- and pP-phase picks. The first data set consists of 3103 events along the Chilean margin and the second one comprises 1680 earthquakes in the Tonga-Fiji subduction zone. Pick data were obtained from the GEOFON earthquake bulletin, produced using data from all available, global station networks. A set of timing corrections varying as a function of source position was calculated for each seismic station. In this way, we could correct the systematic errors introduced into the locations by the inaccuracies in the assumed velocity structure without explicitly solving for a velocity model. Residual statistics show that the median absolute deviation of the travel-time residuals is reduced by 40-60 per cent at regional distances, where the velocity anomalies are strong. Moreover, the spread of the travel-time residuals decreased by similar to 20 per cent at teleseismic distances (>28 degrees). Furthermore, strong variations in initial residuals as a function of recording distance are smoothed out in the final residuals. The relocated catalogues exhibit less scattered locations in depth and sharper images of the seismicity associated with the subducting slabs. Comparison with a high-resolution local catalogue reveals that our relocation process significantly improves the hypocentre locations compared to standard locations. KW - Seismicity and tectonics KW - Computational seismology KW - Subduction zone processes KW - Pacific Ocean KW - South America Y1 - 2016 U6 - https://doi.org/10.1093/gji/ggw405 SN - 0956-540X SN - 1365-246X VL - 208 IS - 2 SP - 589 EP - 602 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Letort, Jean A1 - Retailleau, Lise A1 - Boue, Pierre A1 - Radiguet, Mathilde A1 - Gardonio, Blandine A1 - Cotton, Fabrice Pierre A1 - Campillo, Michel T1 - Lateral variations of the Guerrero-Oaxaca subduction zone (Mexico) derived from weak seismicity (M(b)3.5+) detected on a single array at teleseismic distance JF - Geophysical journal international N2 - Detections of pP and sP phase arrivals (the so-called depth phases) at teleseismic distance provide one of the best ways to estimate earthquake focal depth, as the P-pP and the P-sP delays are strongly dependent on the depth. Based on a new processing workflow and using a single seismic array at teleseismic distance, we can estimate the depth of clusters of small events down to magnitude M-b 3.5. Our method provides a direct view of the relative variations of the seismicity depth from an active area. This study focuses on the application of this new methodology to study the lateral variations of the Guerrero subduction zone (Mexico) using the Eielson seismic array in Alaska (USA). After denoising the signals, 1232 M-b 3.5 + events were detected, with clear P, pP, sP and PcP arrivals. A high-resolution view of the lateral variations of the depth of the seismicity of the Guerero-Oaxaca area is thus obtained. The seismicity is shown to be mainly clustered along the interface, coherently following the geometry of the plate as constrained by the receiver-function analysis along the Meso America Subduction Experiment profile. From this study, the hypothesis of tears on the western part of Guerrero and the eastern part of Oaxaca are strongly confirmed by dramatic lateral changes in the depth of the earthquake clusters. The presence of these two tears might explain the observed lateral variations in seismicity, which is correlated with the boundaries of the slow slip events. KW - North America KW - Time-series analysis KW - Body waves KW - Earthquake source observations KW - Seismicity and tectonics Y1 - 2018 U6 - https://doi.org/10.1093/gji/ggy035 SN - 0956-540X SN - 1365-246X VL - 213 IS - 2 SP - 1002 EP - 1012 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Hendriyana, Andri A1 - Bauer, Klaus A1 - Muksin, Umar A1 - Weber, Michael T1 - AIC-based diffraction stacking for local earthquake locations at the Sumatran Fault (Indonesia) JF - Geophysical journal international N2 - We present a new workflow for the localization of seismic events which is based on a diffraction stacking approach. In order to address the effects from complex source radiation patterns, we suggest to compute diffraction stacking from a characteristic function (CF) instead of stacking the original waveform data. A new CF, which is called in the following mAIC (modified from Akaike Information Criterion) is proposed. We demonstrate that both P- and S-wave onsets can be detected accurately. To avoid cross-talk between P and S waves due to inaccurate velocity models, we separate the P and S waves from the mAIC function by making use of polarization attributes. Then, the final image function is represented by the largest eigenvalue as a result of the covariance analysis between P-and S-image functions. Results from synthetic experiments show that the proposed diffraction stacking provides reliable results. The workflow of the diffraction stacking method was finally applied to local earthquake data from Sumatra, Indonesia. Recordings from a temporary network of 42 stations deployed for nine months around the Tarutung pull-apart basin were analysed. The seismic event locations resulting from the diffraction stacking method align along a segment of the Sumatran Fault. A more complex distribution of seismicity is imaged within and around the Tarutung basin. Two lineaments striking N-S were found in the centre of the Tarutung basin which support independent results from structural geology. KW - Time-series analysis KW - Body waves KW - Computational seismology KW - Earthquake source observations KW - Seismicity and tectonics Y1 - 2018 U6 - https://doi.org/10.1093/gji/ggy045 SN - 0956-540X SN - 1365-246X VL - 213 IS - 2 SP - 952 EP - 962 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Krüger, Frank A1 - Kulikova, Galina A1 - Landgraf, Angela T1 - Magnitudes for the historical 1885 (Belovodskoe), the 1887 (Verny) and the 1889 (Chilik) earthquakes in Central Asia determined from magnetogram recordings JF - Geophysical journal international N2 - Six large magnitude earthquakes in Central Asia which occurred at the end of the 19th century were recorded on early magnetographs in Great Britain. Scalar seismic moment estimates of the 1911 Chon-Kemin, the 1902 Atushi and the 1907 Karatag earthquakes in Central Asia were recently determined by historical seismogram modelling. For those events, we find agreement between moment magnitudes estimated from seismograms and from magnetograms. This supports the assumption of linear scaling of magnetogram amplitudes as function of M-0, which we then use to estimate the moment magnitudes for earlier large-magnitude events, that is, the 1885 Belovodskoe, 1887 Verny and 1889 Chilik earthquakes. The magnetometer data imply that the Chilik earthquake had M(W)7.9, slightly smaller than the Chon-Kemin event with M(W)8.0. The Verny earthquake, however, for which we estimate M(W)7.7, is likely larger than listed in catalogues (M7.3). Similarly, we find a larger magnitude M(W)7.6 (instead of the previous M6.9) for the Belovodskoe earthquake, but this remains uncertain due to measurement imprecision. KW - Earthquake source observations KW - Seismicity and tectonics KW - Intraplate processes Y1 - 2018 U6 - https://doi.org/10.1093/gji/ggy377 SN - 0956-540X SN - 1365-246X VL - 215 IS - 3 SP - 1824 EP - 1840 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Kulikova, Galina A1 - Schurr, Bernd A1 - Krüger, Frank A1 - Brzoska, Elisabeth A1 - Heimann, Sebastian T1 - Source parameters of the Sarez-Pamir earthquake of 1911 February 18 JF - Geophysical journal international N2 - The Ms ∼ 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 (∼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. KW - Earthquake source observations KW - Seismicity and tectonics KW - Body waves KW - Theoretical seismology Y1 - 2016 U6 - https://doi.org/10.1093/gji/ggw069 SN - 0956-540X SN - 1365-246X VL - 205 SP - 1086 EP - 1098 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Kulikova, Galina A1 - Krüger, Frank T1 - Source process of the 1911 M8.0 Chon-Kemin earthquake: investigation results by analogue seismic records JF - Geophysical journal international N2 - Several destructive earthquakes have occurred in Tien-Shan region at the beginning of 20th century. However, the detailed seismological characteristics, especially source parameters of those earthquakes are still poorly investigated. The Chon-Kemin earthquake is the strongest instrumentally recorded earthquake in the Tien-Shan region. This earthquake has produced an approximately 200 km long system of surface ruptures along Kemin-Chilik fault zone and killed about similar to 400 people. Several studies presented the different information on the earthquake epicentre location and magnitude, and two different focal mechanisms were also published. The reason for the limited knowledge of the source parameters for the Chon-Kemin earthquake is the complexity of old analogue records processing, digitization and analysis. In this study the data from 23 seismic stations worldwide were collected and digitized. The earthquake epicentre was relocated to 42.996NA degrees and 77.367EA degrees, the hypocentre depth is estimated between 10 and 20 km. The magnitude was recalculated to m(B) 8.05, M-s 7.94 and M-w 8.02. The focal mechanism, determined from amplitude ratios comparison of the observed and synthetic seismograms, was: str = 264A degrees, dip = 52A degrees, rake = 98A degrees. The apparent source time duration was between similar to 45 and similar to 70 s, the maximum slip occurred 25 s after the beginning of the rupture. Two subevents were clearly detected from the waveforms with the scalar moment ratio between them of about 1/3, the third subevent was also detected with less certainty. Taking into account surface rupture information, the fault geometry model with three patches was proposed. Based on scaling relations we conclude that the total rupture length was between similar to 260 and 300 km and a maximum rupture width could reach similar to 70 km. KW - Earthquake source observations KW - Seismicity and tectonics KW - Body waves KW - Theoretical seismology Y1 - 2015 U6 - https://doi.org/10.1093/gji/ggv091 SN - 0956-540X SN - 1365-246X VL - 201 IS - 3 SP - 1891 EP - 1911 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Passarelli, Luigi A1 - Hainzl, Sebastian A1 - Cesca, Simone A1 - Maccaferri, Francesco A1 - Mucciarelli, Marco A1 - Rößler, Dirk A1 - Corbi, Fabio A1 - Dahm, Torsten A1 - Rivalta, Eleonora T1 - Aseismic transient driving the swarm-like seismic sequence in the Pollino range, Southern Italy JF - Geophysical journal international N2 - Tectonic earthquake swarms challenge our understanding of earthquake processes since it is difficult to link observations to the underlying physical mechanisms and to assess the hazard they pose. Transient forcing is thought to initiate and drive the spatio-temporal release of energy during swarms. The nature of the transient forcing may vary across sequences and range from aseismic creeping or transient slip to diffusion of pore pressure pulses to fluid redistribution and migration within the seismogenic crust. Distinguishing between such forcing mechanisms may be critical to reduce epistemic uncertainties in the assessment of hazard due to seismic swarms, because it can provide information on the frequency-magnitude distribution of the earthquakes (often deviating from the assumed Gutenberg-Richter relation) and on the expected source parameters influencing the ground motion (for example the stress drop). Here we study the ongoing Pollino range (Southern Italy) seismic swarm, a long-lasting seismic sequence with more than five thousand events recorded and located since October 2010. The two largest shocks (magnitude M-w = 4.2 and M-w = 5.1) are among the largest earthquakes ever recorded in an area which represents a seismic gap in the Italian historical earthquake catalogue. We investigate the geometrical, mechanical and statistical characteristics of the largest earthquakes and of the entire swarm. We calculate the focal mechanisms of the M-l > 3 events in the sequence and the transfer of Coulomb stress on nearby known faults and analyse the statistics of the earthquake catalogue. We find that only 25 per cent of the earthquakes in the sequence can be explained as aftershocks, and the remaining 75 per cent may be attributed to a transient forcing. The b-values change in time throughout the sequence, with low b-values correlated with the period of highest rate of activity and with the occurrence of the largest shock. In the light of recent studies on the palaeoseismic and historical activity in the Pollino area, we identify two scenarios consistent with the observations and our analysis: This and past seismic swarms may have been 'passive' features, with small fault patches failing on largely locked faults, or may have been accompanied by an 'active', largely aseismic, release of a large portion of the accumulated tectonic strain. Those scenarios have very different implications for the seismic hazard of the area. KW - Seismicity and tectonics KW - Statistical seismology KW - Dynamics: seismotectonics Y1 - 2015 U6 - https://doi.org/10.1093/gji/ggv111 SN - 0956-540X SN - 1365-246X VL - 201 IS - 3 SP - 1553 EP - 1567 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Ghods, Abdolreza A1 - Shabanian, Esmaeil A1 - Bergman, Eric A1 - Faridi, Mohammad A1 - Donner, Stefanie A1 - Mortezanejad, Gholamreza A1 - Aziz-Zanjani, Asiyeh T1 - The Varzaghan-Ahar, Iran, Earthquake Doublet (M-w 6.4, 6.2): implications for the geodynamics of northwest Iran JF - Geophysical journal international N2 - On 2012 August 11, a pair of large, damaging earthquakes struck the Varzaghan-Ahar region in northwest Iran, in a region where there was no major mapped fault or any well-documented historical seismicity. To investigate the active tectonics of the source region we applied a combination of seismological methods (local aftershock network, calibrated multiple event relocation and focal mechanism studies), field observations (structural geology and geomorphological) and inversions for the regional stress field. The epicentral region is north of the North Tabriz Fault. The first main shock is characterized by right-lateral strike-slip motion on an almost E-W fault plane of about 23 km length extending from the surface to a depth of about 14 km. The second main shock occurred on an ENE-striking fault that dips at 60-70A degrees to the NW. Independent inversions of focal mechanisms and geologically determined fault kinematic data for the active stress state yield a transpressional tectonic regime with sigma(1) oriented N132E. For the region northeast of the North Tabriz Fault, the presence of rigid lithosphere of the South Caspian Basin implies the kinematic adjustment by northward transferring of the contracted masses through both distributed deformation and structural deflections. Our results suggest that the kinematic adjustment inside a contracting wedge may occur along interacting crosswise or conjugate faults to accommodate low rates of internal deformation. At a global scale, our results indicate that despite the basic assumption of 'rigid blocks' in geodetic plate modelling, internal deformation of block-like regions could control the kinematics of deformation and the level of seismic hazard within and around such regions of low deformation rate. KW - Earthquake source observations KW - Seismicity and tectonics KW - Continental neotectonics KW - Continental tectonics: compressional KW - Dynamics: seismotectonics KW - Asia Y1 - 2015 U6 - https://doi.org/10.1093/gji/ggv306 SN - 0956-540X SN - 1365-246X VL - 203 IS - 1 SP - 522 EP - 540 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Grigoli, Francesco A1 - Cesca, Simone A1 - Amoroso, Ortensia A1 - Emolo, Antonio A1 - Zollo, Aldo A1 - Dahm, Torsten T1 - Automated seismic event location by waveform coherence analysis JF - Geophysical journal international N2 - Automated location of seismic events is a very important task in microseismic monitoring operations as well for local and regional seismic monitoring. Since microseismic records are generally characterized by low signal-to-noise ratio, automated location methods are requested to be noise robust and sufficiently accurate. Most of the standard automated location routines are based on the automated picking, identification and association of the first arrivals of P and S waves and on the minimization of the residuals between theoretical and observed arrival times of the considered seismic phases. Although current methods can accurately pick P onsets, the automatic picking of the S onset is still problematic, especially when the P coda overlaps the S wave onset. In this paper, we propose a picking free earthquake location method based on the use of the short-term-average/long-term-average (STA/LTA) traces at different stations as observed data. For the P phases, we use the STA/LTA traces of the vertical energy function, whereas for the S phases, we use the STA/LTA traces of a second characteristic function, which is obtained using the principal component analysis technique. In order to locate the seismic event, we scan the space of possible hypocentral locations and origin times, and stack the STA/LTA traces along the theoretical arrival time surface for both P and S phases. Iterating this procedure on a 3-D grid, we retrieve a multidimensional matrix whose absolute maximum corresponds to the spatial coordinates of the seismic event. A pilot application was performed in the Campania-Lucania region (southern Italy) using a seismic network (Irpinia Seismic Network) with an aperture of about 150 km. We located 196 crustal earthquakes (depth < 20 km) with magnitude range 1.1 < M-L < 2.7. A subset of these locations were compared with accurate manual locations refined by using a double-difference technique. Our results indicate a good agreement with manual locations. Moreover, our method is noise robust and performs better than classical location methods based on the automatic picking of the P and S waves first arrivals. KW - Time-series analysis KW - Inverse theory KW - Earthquake source observations KW - Seismicity and tectonics KW - Early warning Y1 - 2014 U6 - https://doi.org/10.1093/gji/ggt477 SN - 0956-540X SN - 1365-246X VL - 196 IS - 3 SP - 1742 EP - 1753 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Alinaghi, Alireza A1 - Kruger, Frank T1 - Seismic array analysis and redetermination of depths of earthquakes in Tien-Shan: implications for strength of the crust and lithosphere JF - Geophysical journal international N2 - We have redetermined focal depths of moderate and major earthquakes with reported lower-crust and upper-mantle depths that have occurred in Tien-Shan, since the availability of broad-band array data. Records of earthquakes at global arrays have been used for identification and modelling of depth phases in order to make accurate estimation of focal depths. Our results show that half of the purportedly deep earthquakes are indeed originating from depths attributable to middle-crust and lower-crust regions. Also one exceptional event in the northern foreland of Tien-Shan in Junggar Basin is located in the upper mantle at the depth of 64 km. Such unusually deep earthquakes for intraplate continental tectonic domain are all located at the margin of Tien-Shan with its adjacent stable blocks and at least some of them have occurred where the brittle behaviour of continental rocks is not highly expected. The reverse mechanisms of all these earthquakes and their proximity to formerly subducting and later colliding and underplating stable blocks and their interactions with overlying Tien-Shan are clues to explain this extremity. KW - Earthquake source observations KW - Seismicity and tectonics KW - Body waves Y1 - 2014 U6 - https://doi.org/10.1093/gji/ggu141 SN - 0956-540X SN - 1365-246X VL - 198 IS - 2 SP - 1111 EP - 1129 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Lange, Dietrich A1 - Bedford, J. R. A1 - Moreno, M. A1 - Tilmann, F. A1 - Báez, Juan Carlos A1 - Bevis, M. A1 - Krüger, Frank T1 - Comparison of postseismic afterslip models with aftershock seismicity for three subduction-zone earthquakes: Nias 2005, Maule 2010 and Tohoku 2011 JF - Geophysical journal international N2 - We focus on the relation between seismic and total postseismic afterslip following the Maule M-w 8.8 earthquake on 2010 February 27 in central Chile. First, we calculate the cumulative slip released by aftershock seismicity. We do this by summing up the aftershock regions and slip estimated from scaling relations. Comparing the cumulative seismic slip with afterslip modelswe showthat seismic slip of individual aftershocks exceeds locally the inverted afterslip model from geodetic constraints. As the afterslip model implicitly contains the displacements from the aftershocks, this reflects the tendency of afterslip models to smear out the actual slip pattern. However, it also suggests that locally slip for a number of the larger aftershocks exceeds the aseismic slip in spite of the fact that the total equivalent moment of the afterslip exceeds the cumulative moment of aftershocks by a large factor. This effect, seen weakly for the Maule 2010 and also for the Tohoku 2011 earthquake, can be explained by taking into account the uncertainties of the seismicity and afterslip models. In spite of uncertainties, the hypocentral region of the Nias 2005 earthquake is suggested to release a large fraction of moment almost purely seismically. Therefore, these aftershocks are not driven solely by the afterslip but instead their slip areas have probably been stressed by interseismic loading and the mainshock rupture. In a second step, we divide the megathrust of the Maule 2010 rupture into discrete cells and count the number of aftershocks that occur within 50 km of the centre of each cell as a function of time. We then compare this number to a time-dependent afterslip model by defining the 'afterslip to aftershock ratio' (ASAR) for each cell as the slope of the best fitting line when the afterslip at time t is plotted against aftershock count. Although we find a linear relation between afterslip and aftershocks for most cells, there is significant variability in ASAR in both the downdip and along-strike directions of the megathrust. We compare the spatial distribution of ASAR with the spatial distribution of seismic coupling, coseismic slip and Bouguer gravity anomaly, and in each case we find no significant correlation. KW - Creep and deformation KW - Earthquake dynamics KW - Seismicity and tectonics KW - Continental margins: convergent Y1 - 2014 U6 - https://doi.org/10.1093/gji/ggu292 SN - 0956-540X SN - 1365-246X VL - 199 IS - 2 SP - 784 EP - 799 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Zöller, Gert A1 - Ben-Zion, Yehuda T1 - Large earthquake hazard of the San Jacinto fault zone, CA, from long record of simulated seismicity assimilating the available instrumental and paleoseismic data JF - Pure and applied geophysics N2 - We investigate spatio-temporal properties of earthquake patterns in the San Jacinto fault zone (SJFZ), California, between Cajon Pass and the Superstition Hill Fault, using a long record of simulated seismicity constrained by available seismological and geological data. The model provides an effective realization of a large segmented strike-slip fault zone in a 3D elastic half-space, with heterogeneous distribution of static friction chosen to represent several clear step-overs at the surface. The simulated synthetic catalog reproduces well the basic statistical features of the instrumental seismicity recorded at the SJFZ area since 1981. The model also produces events larger than those included in the short instrumental record, consistent with paleo-earthquakes documented at sites along the SJFZ for the last 1,400 years. The general agreement between the synthetic and observed data allows us to address with the long-simulated seismicity questions related to large earthquakes and expected seismic hazard. The interaction between m a parts per thousand yen 7 events on different sections of the SJFZ is found to be close to random. The hazard associated with m a parts per thousand yen 7 events on the SJFZ increases significantly if the long record of simulated seismicity is taken into account. The model simulations indicate that the recent increased number of observed intermediate SJFZ earthquakes is a robust statistical feature heralding the occurrence of m a parts per thousand yen 7 earthquakes. The hypocenters of the m a parts per thousand yen 5 events in the simulation results move progressively towards the hypocenter of the upcoming m a parts per thousand yen 7 earthquake. KW - Earthquake dynamics KW - Earthquake interaction KW - forecasting KW - prediction KW - Statistical seismology KW - Seismicity and tectonics Y1 - 2014 U6 - https://doi.org/10.1007/s00024-014-0783-1 SN - 0033-4553 SN - 1420-9136 VL - 171 IS - 11 SP - 2955 EP - 2965 PB - Springer CY - Basel ER - TY - JOUR A1 - Braeuer, B. A1 - Asch, Günter A1 - Hofstetter, Rami A1 - Haberland, Christian A1 - Jaser, Darweesh A1 - El-Kelani, Radwan J.. A1 - Weber, Michael H. T1 - Microseismicity distribution in the southern Dead Sea basin and its implications on the structure of the basin JF - Geophysical journal international N2 - While the Dead Sea basin has been studied for a long time, the available knowledge about the detailed seismicity distribution in the area, as well as the deeper structure of the basin, is limited. Therefore, within the framework of the international project DESIRE (DEad Sea Integrated REsearch project), a dense temporary local seismological network was operated in the southern Dead Sea area. We use 530 local earthquakes, having all together 26 730 P- and S-arrival times for a simultaneous inversion of 1-D velocity models, station corrections and precise earthquake locations. Jackknife tests suggest an accuracy of the derived hypocentre locations of about 1 km. Thus, the result is the first clear image of the absolute distribution of the microseismicity of the area, especially in depth. The seismicity is concentrated in the upper crust down to 20 km depth while the lower limit of the seismicity is reached at 31 km depth. The seismic events at the eastern boundary fault (EBF) in the southern part of the study area represent the northward transform motion of the Arabian Plate along the Dead Sea Transform. North of the Boqeq fault the seismic activity represents the transfer of the motion in the pull-apart basin from the eastern to the western boundary. We find that from the surface downward the seismic events are tracing the boundary faults of the basin. The western boundary is mapped down to 12 km depth while the EBF reaches about 17 km depth, forming an asymmetric basin. One fifth of the data set is related to a specific cluster in time and space, which occurred in 2007 February at the western border fault. This cluster is aligned vertically, that is, it is perpendicular to the direction of the dominating left-lateral strike-slip movement at the main transform fault. KW - Seismicity and tectonics KW - Continental tectonics: strike-slip and transform KW - Asia Y1 - 2012 U6 - https://doi.org/10.1111/j.1365-246X.2011.05318.x SN - 0956-540X VL - 188 IS - 3 SP - 873 EP - 878 PB - Wiley-Blackwell CY - Malden ER - TY - JOUR A1 - Pesicek, J. D. A1 - Engdahl, E. R. A1 - Thurber, C. H. A1 - DeShon, H. R. A1 - Lange, Dietrich T1 - Mantle subducting slab structure in the region of the 2010 M8.8 Maule earthquake (30-40 degrees S), Chile JF - Geophysical journal international N2 - We present a new tomographic model of the mantle in the area of the 2010 M8.8 Maule earthquake and surrounding regions. Increased ray coverage provided by the aftershock data allows us to image the detailed subducting slab structure in the mantle, from the region of flat slab subduction north of the Maule rupture to the area of overlapping rupture between the 1960 M9.5 and the 2010 M8.8 events to the south. We have combined teleseismic primary and depth phase arrivals with available local arrivals to better constrain the teleseismic earthquake locations in the region, which we use to conduct nested regionalglobal tomography. The new model reveals the detailed structure of the flat slab and its transition to a more moderately dipping slab in the Maule region. South of the Maule region, a steeply dipping relic slab is imaged from similar to 200 to 1000 km depth that is distinct from the moderately dipping slab above it and from the more northerly slab at similar depths. We interpret the images as revealing both horizontal and vertical tearing of the slab at similar to 38 degrees S to explain the imaged pattern of slab anomalies in the southern portion of the model. In contrast, the transition from a horizontal to moderately subducting slab in the northern portion of the model is imaged as a continuous slab bend. We speculate that the tearing was most likely facilitated by a fracture zone in the downgoing plate or alternatively by a continental scale terrane boundary in the overriding plate. KW - Seismicity and tectonics KW - Seismic tomography KW - Subduction zone processes Y1 - 2012 U6 - https://doi.org/10.1111/j.1365-246X.2012.05624.x SN - 0956-540X VL - 191 IS - 1 SP - 317 EP - 324 PB - Wiley-Blackwell CY - Hoboken ER - TY - JOUR A1 - Hainzl, Sebastian A1 - Zöller, Gert A1 - Brietzke, Gilbert B. A1 - Hinzen, Klaus-G. T1 - Comparison of deterministic and stochastic earthquake simulators for fault interactions in the Lower Rhine Embayment, Germany JF - Geophysical journal international N2 - Time-dependent probabilistic seismic hazard assessment requires a stochastic description of earthquake occurrences. While short-term seismicity models are well-constrained by observations, the recurrences of characteristic on-fault earthquakes are only derived from theoretical considerations, uncertain palaeo-events or proxy data. Despite the involved uncertainties and complexity, simple statistical models for a quasi-period recurrence of on-fault events are implemented in seismic hazard assessments. To test the applicability of statistical models, such as the Brownian relaxation oscillator or the stress release model, we perform a systematic comparison with deterministic simulations based on rate- and state-dependent friction, high-resolution representations of fault systems and quasi-dynamic rupture propagation. For the specific fault network of the Lower Rhine Embayment, Germany, we run both stochastic and deterministic model simulations based on the same fault geometries and stress interactions. Our results indicate that the stochastic simulators are able to reproduce the first-order characteristics of the major earthquakes on isolated faults as well as for coupled faults with moderate stress interactions. However, we find that all tested statistical models fail to reproduce the characteristics of strongly coupled faults, because multisegment rupturing resulting from a spatiotemporally correlated stress field is underestimated in the stochastic simulators. Our results suggest that stochastic models have to be extended by multirupture probability distributions to provide more reliable results. KW - Earthquake interaction KW - forecasting KW - and prediction KW - Seismicity and tectonics KW - Statistical seismology Y1 - 2013 U6 - https://doi.org/10.1093/gji/ggt271 SN - 0956-540X SN - 1365-246X VL - 195 IS - 1 SP - 684 EP - 694 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Sen, Ali Tolga A1 - Cesca, Simone A1 - Bischoff, Monika A1 - Meier, Thomas A1 - Dahm, Torsten T1 - Automated full moment tensor inversion of coal mining-induced seismicity JF - Geophysical journal international N2 - Seismicity induced by coal mining in the Ruhr region, Germany, has been monitored continuously over the last 25 yr. In 2006, a dense temporary network (HAMNET) was deployed to locally monitor seismicity induced by longwall mining close to the town of Hamm. Between 2006 July and 2007 July, more than 7000 events with magnitudes M-L from -1.7 to 2.0 were detected. The spatiotemporal distribution of seismicity shows high correlation with the mining activity. In order to monitor rupture processes, we set up an automated source inversion routine and successfully perform double couple and full moment tensor (MT) inversions for more than 1000 events with magnitudes above M-L -0.5. The source inversion is based on a full waveform approach, both in the frequency and in the time domain, providing information about the centroid location, focal mechanism, scalar moment and full MT. Inversion results indicate a strong dominance of normal faulting focal mechanisms, with a steeper plane and a subhorizontal one. Fault planes are oriented parallel to the mining stopes. We classify the focal mechanisms based on their orientation and observe different frequency-magnitude distributions for families of events with different focal mechanisms; the overall frequency-magnitude distribution is not fitting the Gutenberg-Richter relation. Full MTs indicate that non-negligible opening tensile components accompanied normal faulting source mechanisms. Finally, extended source models are investigated for largest events. Results suggest that the rupture processes mostly occurred along the subvertical planes. KW - Geomechanics KW - Fracture and flow KW - Earthquake source observations KW - Seismicity and tectonics Y1 - 2013 U6 - https://doi.org/10.1093/gji/ggt300 SN - 0956-540X SN - 1365-246X VL - 195 IS - 2 SP - 1267 EP - 1281 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Muksin, Umar A1 - Haberland, Christian A1 - Bauer, Klaus A1 - Weber, Michael H. T1 - Three-dimensional upper crustal structure of the geothermal system in Tarutung (North Sumatra, Indonesia) revealed by seismic attenuation tomography JF - Geophysical journal international N2 - The geothermal potential in Tarutung is controlled by both the Sumatra Fault system and young arc volcanism. In this study we use the spatial distribution of seismic attenuation, calculated from local earthquake recordings, to image the 3-D seismic attenuation of the area and relate it with the temperature anomalies and the fluid distribution of the subsurface. A temporary seismic network of 42 stations was deployed around Tarutung and Sarulla (south of Tarutung) for a period of 10 months starting in 2011 May. Within this period, the network recorded 2586 local events. A high-quality subset of 229 events recorded by at least 10 stations was used for the attenuation inversion (tomography). Path-average attenuation (t(p)*) was calculated by using a spectral inversion method. The spread function, the contour lines of the model resolution matrix and the recovery test results show that our 3-D attenuation model (Q(p)) has good resolution around the Tarutung Basin and along the Sarulla graben. High attenuation (low Q(p)) related to the geothermal system is found in the northeast of the Tarutung Basin suggesting fluid pathways from below the Sumatra Fault. The upper part of the studied geothermal system in the Tarutung district seems to be mainly controlled by the fault structure rather than by magmatic activities. In the southwest of the Tarutung Basin, the high attenuation zone is associated with the Martimbang volcano. In the Sarulla region, a low-Q(p) anomaly is found along the graben within the vicinity of the Hopong caldera. KW - Seismicity and tectonics KW - Body waves KW - Seismic attenuation KW - Seismic tomography Y1 - 2013 U6 - https://doi.org/10.1093/gji/ggt383 SN - 0956-540X SN - 1365-246X VL - 195 IS - 3 SP - 2037 EP - 2049 PB - Oxford Univ. Press CY - Oxford ER -