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 - 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 - Ramos, C. A1 - Mechie, James A1 - Feng, M. T1 - Shear wave velocity and Poisson's ratio models across the southern Chile convergent margin at 38°15′S JF - Geophysical journal international N2 - Using active and passive seismology data we derive a shear (S) wave velocity model and a Poisson's ratio (σ) model across the Chilean convergent margin along a profile at 38°15′S, where the Mw 9.5 Valdivia earthquake occurred in 1960. The derived S-wave velocity model was constructed using three independently obtained velocity models that were merged together. In the upper part of the profile (0–2 km depth), controlled source data from explosions were used to obtain an S-wave traveltime tomogram. For the middle part (2–20 km depth), data from a temporary seismology array were used to carry out a dispersion analysis. The resulting dispersion curves were used to obtain a 3-D S-wave velocity model. In the lower part (20–75 km depth, depending on the longitude), an already existent local earthquake tomographic image was merged with the other two sections. This final S-wave velocity model and already existent compressional (P) wave velocity models along the same transect allowed us to obtain a Poisson's ratio model. The results of this study show that the velocities and Poisson's ratios in the continental crust of this part of the Chilean convergent margin are in agreement with geological features inferred from other studies and can be explained in terms of normal rock types. There is no requirement to call on the existence of measurable amounts of present-day fluids, in terms of seismic velocities, above the plate interface in the continental crust of the Coastal Cordillera and the Central Valley in this part of the Chilean convergent margin. This is in agreement with a recent model of water being transported down and released from the subduction zone. KW - Controlled source seismology KW - Seismic tomography KW - Computational seismology KW - South America Y1 - 2016 U6 - https://doi.org/10.1093/gji/ggv541 SN - 0956-540X SN - 1365-246X VL - 204 SP - 1620 EP - 1635 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Hobiger, M. A1 - Cornou, C. A1 - Wathelet, M. A1 - Di Giulio, G. A1 - Knapmeyer-Endrun, B. A1 - Renalier, F. A1 - Bard, Pierre-Yves A1 - Savvaidis, Alexandros A1 - Hailemikael, S. A1 - Le Bihan, N. A1 - Ohrnberger, Matthias A1 - Theodoulidis, N. T1 - Ground structure imaging by inversions of Rayleigh wave ellipticity sensitivity analysis and application to European strong-motion sites JF - Geophysical journal international N2 - The knowledge of the local soil structure is important for the assessment of seismic hazards. A widespread, but time-consuming technique to retrieve the parameters of the local underground is the drilling of boreholes. Another way to obtain the shear wave velocity profile at a given location is the inversion of surface wave dispersion curves. To ensure a good resolution for both superficial and deeper layers, the used dispersion curves need to cover a wide frequency range. This wide frequency range can be obtained using several arrays of seismic sensors or a single array comprising a large number of sensors. Consequently, these measurements are time-consuming. A simpler alternative is provided by the use of the ellipticity of Rayleigh waves. The frequency dependence of the ellipticity is tightly linked to the shear wave velocity profile. Furthermore, it can be measured using a single seismic sensor. As soil structures obtained by scaling of a given model exhibit the same ellipticity curve, any inversion of the ellipticity curve alone will be ambiguous. Therefore, additional measurements which fix the absolute value of the shear wave velocity profile at some points have to be included in the inversion process. Small-scale spatial autocorrelation measurements or MASW measurements can provide the needed data. Using a theoretical soil structure, we show which parts of the ellipticity curve have to be included in the inversion process to get a reliable result and which parts can be omitted. Furthermore, the use of autocorrelation or high-frequency dispersion curves will be highlighted. The resulting guidelines for inversions including ellipticity data are then applied to real data measurements collected at 14 different sites during the European NERIES project. It is found that the results are in good agreement with dispersion curve measurements. Furthermore, the method can help in identifying the mode of Rayleigh waves in dispersion curve measurements. KW - Inverse theory KW - Surface waves and free oscillations KW - Site effects KW - Computational seismology KW - Wave propagation Y1 - 2013 U6 - https://doi.org/10.1093/gji/ggs005 SN - 0956-540X VL - 192 IS - 1 SP - 207 EP - 229 PB - Oxford Univ. Press CY - Oxford ER -