Description: The Van (Eastern Anatolia, Turkey) earthquake occurred on Sunday, October 23, 2011 with a moment magnitude of 7.2. The tectonics of this region is characterized by strike–slip faulting on the Bitlis Suture Zone, and thrusting in the Zagros fold and thrust belt. Using high-rate (1 second) GPS data from permanent GNSS stations from the CORS-TR network, co-seismic displacements of eleven stations were determined using precise point positioning during this earthquake. We used the time series of coordinate changes for fourteen CORS-TR stations, and calculated the crust movements before and after the earthquake. According to the PPP solutions computed using high frequency GPS data to determine the co-seismic motions of stations, we conclude for the Van earthquake an occurrence time of 10:41:22 (UTC). No pre-seismic horizontal movement of stations at the level more than 5 mm before the earthquake could be observed. That means that no kinematic warning or prediction before the earthquake exists. Along an east–west horizontal line north of the Van Sea with a length of about 100 km, the northern part of this line experienced extension of 0.2–1 ppm in a NW–SE direction. The southern part experienced N–S shortening of 0.5–1.5 ppm. The N–S shortening we estimated geodetically matches well with the N–S shortening and thrust focal mechanism derived independently using seismic data by the USGS. Co-seismic surface displacements derived from the GPS data are consistent with the teleseismic source model given by the USGS. The geodetic source model derived from the GPS data reproduces the same moment magnitude and centroid as the teleseismic model, but shows a higher spatial resolution of the slip distribution. We also analyzed the post-seismic surface displacements derived from the GPS data within the first two weeks after the mainshock. No reasonable slip distribution on the co-seismic fault plane could be found, indicating that the sources for the early post-seismic deformation might come from the widely scattered aftershocks.

Description: We determined a high-resolution 3-D S-wave velocity model for a 26 km × 12 km area in the northern part of the basin of Santiago de Chile. To reach this goal, we used microtremor recordings at 125 sites for deriving the horizontal-to-vertical (H/V) spectral ratios that we inverted to retrieve local S-wave velocity profiles. In the inversion procedure, we used additional geological and geophysical constraints and values of the thickness of the sedimentary cover already determined by gravimetric measurements, which were found to vary substantially over short distances in the investigated area. The resulting model was derived by interpolation with a kriging technique between the single S-wave velocity profiles and shows locally good agreement with the few existing velocity profile data, but allows the entire area, as well as deeper parts of the basin, to be represented in greater detail. The wealth of available data allowed us to check if any correlation between the S-wave velocity in the uppermost 30 m (v30S) and the slope of topography, a new technique recently proposed by Wald and Allen, exists on a local scale. We observed that while one lithology might provide a greater scatter in the velocity values for the investigated area, almost no correlation between topographic gradient and calculated v30S exists, whereas a better link is found between v30S and the local geology. Finally, we compared the v30S distribution with the MSK intensities for the 1985 Valparaiso event, pointing out that high intensities are found where the expected v30S values are low and over a thick sedimentary cover. Although this evidence cannot be generalized for all possible earthquakes, it indicates the influence of site effects modifying the ground motion when earthquakes occur well outside of the Santiago basin.

Description: Amplitude ratio of 30 short-period conspicuous P5KP and PKPab phases from five intermediate depth or deep events in Fiji-Tonga recorded at European stations around 150° distance shows a mean value two to three times the ratio of the synthetic amplitudes obtained by the normal-mode theory (and ak135 model) or by full-wave theory (and PREM). There is a large variance in the results, also observed in five amplitude ratios from one event in Argentina observed at temporary stations in China around 156°. Global recordings of three major deep earthquakes in Fiji, Bonin, and Western Brazil observed at ASAR, WRA, and ZRNK arrays, at 59 North America stations and at six South Pole stations displayed conspicuous P4KP and PcP (or ScP) phases. The amplitude ratio values of P4KP vs P(S)cP are sometimes almost one order of magnitude larger than the corresponding values of the synthetics. In both cases, arrival times and slowness values (corrected for ellipticity and station elevation) at the distances up to 23° beyond the A cutoff point predicted by ray theory match both the synthetics, suggesting the observations are the AB branch of PmKP (m = 4, 5) around 1 Hz. In disagreement to ray theory, no reliable BC branch is observed neither on the recordings nor on the normal-mode synthetics. The high amplitude ratio values cannot be explained by realistic perturbations of the velocity or attenuation values of the global models in the proximity of the core-to-mantle boundary (CMB). We speculate that the focusing effects and/or strong scattering most likely associated to some anomalous velocity areas of the lowermost mantle are responsible for that. The results suggest limitations of the previous evaluations of the short-period attenuation in the outer core from PmKP amplitudes (m ≥ 3), irrespective of the fact that they are obtained by using ray theory, normal-mode or full-wave synthetics. Attempts to use PmKP arrival times in order to refine velocity structure in the proximity of CMB should be also regarded with care if the propagation times have been computed with ray theory.

Description: We propose a new rapid procedure for determining the energy magnitude Me for shallow events from broadband teleseismic P-wave signals within the distance range 20°–98°. To accomplish this task, we compute spectral amplitude decay functions for different periods using numerical simulations based on the reference Earth model AK135Q. By means of these functions, we correct the spectra of the teleseismic recordings for the propagation path effects, and calculate the radiated seismic energy ES, and hence Me. We use cumulative P-wave windows for simulating a real- or near real-time procedure and test it for 61 shallow earthquakes. The results show that our approach is able to provide a rapid and reliable Me determination within 7–15 minutes after the earthquake origin time, and is therefore suitable for implementation in rapid response systems.

Description: Near‐field ground‐motion data are available in semi‐real time either from modern strong‐motion or continuous Global Positioning System (GPS) networks, allowing robust solutions for earthquake source parameters, which are useful for rapid disaster assessment and early warning. These wide applications require the ground‐motion data to cover a very broad frequency band that, however, is usually not available. This paper presents a case study on the 2011 Mw 9.0 Tohoku earthquake, showing how the ground‐motion information from geodetic and seismic instrumentations is complementary, and suggesting the joint use of both types of data, particularly when the network coverage is sparse. First the strong‐motion records from the two Japanese networks, K‐NET and KiK‐Net, are analyzed using an automatic empirical baseline correction tool. The static coseismic displacement data are obtained by double integration and then used to derive the permanent slip distribution on the earthquake fault. Comparisons with the corresponding GPS‐based solutions yield a quantitative estimation of uncertainties of the empirical baseline correction. Furthermore, a dozen nearby GPS and strong‐motion station pairs are selected to demonstrate that the information in their time series agrees with each other. Finally, methods for combining both types of ground‐motion observation systems are discussed, and the wide applicability of this approach is highlighted.

Description: Mt. Merapi is one of the most dangerous volcanoes in Indonesia, located within the tectonically active region of south-central Java. This study investigates how Mt. Merapi affected - and was affected by - nearby tectonic earthquakes. In 2001, a Mw6.3 earthquake occurred in conjunction with an increase in fumarole temperature at Mt. Merapi. In 2006, another Mw6.3 earthquake took place, concomitant with an increase of magma extrusion and pyroclastic flows. Here, we develop theoretical models to study the amount of stress transfer between the earthquakes and the volcano, showing that dynamic, rather than static, stress changes are likely responsible for the temporal and spatial proximity of these events. Our examination of the 2001 and 2006 events implies that volcanic activity at Mt. Merapi is influenced by stress changes related to remote tectonic earthquakes, a finding that is important for volcano hazard assessment in this densely inhabited area.

Description: It is common practice in the seismological community to use, especially for large earthquakes, the moment magnitude Mw as a unique magnitude parameter to evaluate the earthquake’s damage potential. However, as a static measure of earthquake size, Mwdoes not provide direct information about the released seismicwave energy and its high frequency content,which is the more interesting information both for engineering purposes and for a rapid assessment of the earthquake’s shaking potential. Therefore, we recommend to provide to disaster management organizations besides Mw also sufficiently accurate energy magnitude determinations as soon as possible after large earthquakes. We developed and extensively tested a rapid method for calculating the energy magnitude Me within about 10–15 min after an earthquake’s occurrence. The method is based on pre-calculated spectral amplitude decay functions obtained from numerical simulations of Green’s functions. After empirical validation, the procedure has been applied offline to a large data set of 767 shallow earthquakes that have been grouped according to their type of mechanism (strike-slip, normal faulting, thrust faulting, etc.). The suitability of the proposed approach is discussed by comparing our rapid Me estimates with Mw published by GCMT as well as with Mw and Me reported by the USGS. Mw is on average slightly larger than our Me for all types of mechanisms. No clear dependence on source mechanism is observed for our Me estimates. In contrast, Me from the USGS is generally larger than Mw for strike-slip earthquakes and generally smaller for the other source types. For ∼67 per cent of the event data set our Me differs ≤ ±0.3 magnitude units (m.u.) from the respective Me values published by the USGS. However, larger discrepancies (up to 0.8 m.u.) may occur for strike-slip events. A reason of that may be the overcorrection of the energy flux applied by the USGS for this type of earthquakes. We follow the original definition of magnitude scales, which does not apply a priori mechanism corrections to measured amplitudes, also since reliable fault-plane solutions are hardly available within 10–15 min after the earthquake origin time. Notable is that our uncorrected Me data show a better linear correlation and less scatter with respect to Mw than Me of the USGS. Finally, by analysing the recordings of representative recent pairs of strong and great earthquakes, we emphasize the importance of combining Mw and Me in the rapid characterization of the seismic source. They are related to different aspects of the source and may differ occasionally even more than 1 m.u. This highlights the usefulness and importance of providing these two magnitude estimates together for a better assessment of an earthquake’s shaking potential and/or tsunamigenic potential.

Description: In 2007 a M7.7 earthquake occurred near the town of Tocopilla within the northern Chile seismic gap. Mainshock slip, derived from coseismic surface deformation, was confined to the depth range between 30-55 km. We relocated ~1100 events during six months before and one week after the mainshock. Aftershock seismicity is first congruent to the mainshock slip and then it spreads offshore west and northwest of Mejillones Peninsula (MP). Waveform modeling for 38 aftershocks reveals source mechanisms that are in the majority similar to the mainshock. However, a few events appear to occur in the upper plate, some with extensional mechanisms. Juxtaposing the Tocopilla aftershocks with those following the neighboring 1995 Antofagasta earthquake produces a striking symmetry across an EW axis in the center of MP. Events seem to skirt around MP, probably due to a shallower Moho there. We suggest that the seismogenic coupling zone in northern Chile changes its frictional behavior in the down-dip direction from unstable to mostly conditionally stable. For both earthquake sequences, aftershocks agglomerate in the conditionally stable region, whereas maximum inter-seismic slip deficit and co-seismic slip occurs in the unstable region. The boundary between the unstable and conditionally stable zones parallels the coastline. We identify a similar segmentation for other earthquakes in Chile and Peru, where the offshore segments break in great M〉8 earthquakes, and the onshore segments in smaller M〈8 earthquakes. Using critical taper analysis, we demonstrate a causal relationship between varying slip behavior on the interface and forearc wedge anatomy that can be attributed to spatial variations in the rate-dependency of friction.