Source characteristics of the deep Philippine earthquake cluster, 23 and 24 July 2010

We determine the rupture velocity, rupture area, stress drop and duration of four strong deep-focus earthquakes in the Philippines by back-projecting the teleseismic P waves. Four deep-focus earthquakes occurred in a totally consumed Molucca microplate; their focal depths were greater than 550 km and their moment magnitudes were between M _w 6.6 and M _w 7.6. By studying this deep-focus cluster, we are able to estimate the rupture velocity, rupture area and stress drop which would assist in constraining the physical mechanism for earthquakes deeper than 500 km. Since the Molucca microplate is totally consumed, little evidence is left on the surface for us to do research. This deep-focus cluster provides us the opportunity to reveal the properties of this totally consumed microplate by using seismic method for the first time. Four earthquakes in this deep-focus cluster all have multiple rupture subevents. The M _w 7.3 event ruptures in two subevents, the M _w 7.6 and M _w 7.4 events both have three subevents. The M _w 6.6 event has single peak on the amplitude as a function of time; however, its energy releases at two spatially separated areas. Our results show that this deep-focus cluster has a slow rupture velocity which is about 0.27 to 0.43 of the shear wave velocity, long-scaled duration, concentrated energy release area, and high stress drop. These source properties are similar to those of other deep earthquakes occurring in warm slabs and indicate that the totally consumed Molucca microplate possibly is a warm plate.     

Validation of ground‐motion simulations for historical events using MDoF systems

The study presented in this paper addresses the issue of engineering validation of Graves and Pitarka's (2010) hybrid broadband ground motion simulation methodology with respect to some well‐recorded historical events and considering the response of multiple degrees of freedom (MDoF) systems. Herein, validation encompasses detailed assessment of how similar is, for a given event, the seismic response due to comparable hybrid broadband simulated records and real records. In the first part of this study, in order to investigate the dynamic response of a wide range of buildings, MDoF structures are modeled as elastic continuum systems consisting of a combination of a flexural cantilever beam coupled with a shear cantilever beam. A number of such continuum systems are selected including the following: (1) 16 oscillation periods between 0.1 and 6 s; (2) three shear to flexural deformation ratios to represent respectively shear‐wall structures, dual systems, and moment‐resisting frames; and (3) two stiffness distributions along the height of the systems, that is, uniform and linear. Demand spectra in terms of generalized maximum interstory drift ratio (IDR) and peak floor acceleration (PFA) are derived using simulations and actual recordings for four historical earthquakes, namely, the 1979 M_w 6.5 Imperial Valley earthquake, 1989 M_w 6.8 Loma Prieta earthquake, 1992 M_w 7.2 Landers earthquake, and 1994 M_w 6.7 Northridge earthquake. In the second part, for two nonlinear case study structures, the IDR and PFA distributions over the height and their statistics, are obtained and compared for both recorded and simulated time histories. These structures are steel moment frames designed for high seismic hazard, 20‐story high‐rise and 6‐story low‐rise buildings. The results from this study highlight the similarities and differences between simulated and real records in terms of median and intra‐event standard deviation of logs of seismic demands for MDoF building systems. This general agreement, in a broad range of moderate and long periods, may provide confidence in the use of the simulation methodology for engineering applications, whereas the discrepancies, statistically significant only at short periods, may help in addressing improvements in generation of synthetic records. Copyright © 2013 John Wiley & Sons, Ltd.     

Estimation of strong ground motions from hypothetical earthquakes on the Cascadia subduction zone,Pacific Northwest

Strong ground motions are estimated for the Pacific Northwest assuming that large shallow earthquakes, similar to those experienced in southern Chile, southwestern Japan, and Colombia, may also occur on the Cascadia subduction zone. Fifty-six strong motion recordings for twenty-five subduction earthquakes ofM _s7.0 are used to estimate the response spectra that may result from earthquakesM _w<81/4. Large variations in observed ground motion levels are noted for a given site distance and earthquake magnitude. When compared with motions that have been observed in the western United States, large subduction zone earthquakes produce relatively large ground motions at surprisingly large distances. An earthquake similar to the 22 May 1960 Chilean earthquake (M _w 9.5) is the largest event that is considered to be plausible for the Cascadia subduction zone. This event has a moment which is two orders of magnitude larger than the largest earthquake for which we have strong motion records. The empirical Green's function technique is used to synthesize strong ground motions for such giant earthquakes. Observed teleseismicP-waveforms from giant earthquakes are also modeled using the empirical Green's function technique in order to constrain model parameters. The teleseismic modeling in the period range of 1.0 to 50 sec strongly suggests that fewer Green's functions should be randomly summed than is required to match the long-period moments of giant earthquakes. It appears that a large portion of the moment associated with giant earthquakes occurs at very long periods that are outside the frequency band of interest for strong ground motions. Nevertheless, the occurrence of a giant earthquake in the Pacific Northwest may produce quite strong shaking over a very large region.     

Large earthquake doublets and fault plane heterogeneity in the northern Solomon Islands subduction zone

In the Solomon Islands and New Britain subduction zones, the largest earthquakes commonly occur as pairs with small separation in time, space and magnitude. This doublet behavior has been attributed to a pattern of fault plane heterogeneity consisting of closely spaced asperities such that the failure of one asperity triggers slip in adjacent asperities. We analyzed body waves of the January 31, 1974,M _w =7.3, February 1, 1974,M _w =7.4, July 20, 1975 (1437)M _w =7.6 and July 20, 1975 (1945),M _w =7.3 doublet events using an iterative, multiple station inversion technique to determine the spatio-temporal distribution of seismic moment release associated with these events. Although the 1974 doublet has smaller body wave moments than the 1975 events, their source histories are more complicated, lasting over 40 seconds and consisting of several subevents located near the epicentral regions. The second 1975 event is well modeled by a simple point source initiating at a depth of 15 km and rupturing an approximate 20 km region about the epicenter. The source history of the first 1975 event reveals a westerly propagating rupture, extending about 50 km from its hypocenter at a depth of 25 km. The asperities of the 1975 events are of comparable size and do not overlap one another, consistent with the asperity triggering hypothesis. The relatively large source areas and small seismic moments of the 1974 doublet events indicate failure of weaker portions of the fault plane in their epicentral regions. Variations in the roughness of the bathymetry of the subducting plate, accompanying subduction of the Woodlark Rise, may be responsible for changes in the mechanical properties of the plate interface.To understand how variations in fault plane coupling and strength affect the interplate seismicity pattern, we relocated 85 underthrusting earthquakes in the northern Solomon Islands Are since 1964. Relatively few smaller magnitude underthrusting events overlap the Solomon Islands doublet asperity regions, where fault coupling and strength are inferred to be the greatest. However, these asperity regions have been the sites of several previous earthquakes withM _s 7.0. The source regions of the 1974 doublet events, which we infer to be mechanically weak, contain many smaller magnitude events but have not generated any otherM _s 7.0 earthquakes in the historic past. The central portion of the northern Solomon Islands Arc between the two largest doublet events in 1971 (studied in detail bySchwartz et al., 1989a) and 1975 contains the greatest number of smaller magnitude underthrusting earthquakes. The location of this small region sandwiched between two strongly coupled portions of the plate interface suggest that it may be the site of the next large northern Solomon Islands earthquake. However, this region has experienced no known earthquakes withM _s 7.0 and may represent a relatively aseismic portion of the subduction zone.     

Wavelet-Hilbert transform-based simulation of pulse-like ground motion

In nonlinear dynamic structural analysis, a suite of pulse-like ground motions is required for the performance-based design of structures near active faults. The dissimilarity in the amplitude and frequency content of the earthquake time series referred to nonstationary properties in temporal and spectral, respectively. An approach is proposed based on the nonstationary properties of the far-field records and the seismological information in an event for simulating pulse-like records. The pulse-like earthquake time history is estimated via the superposition of the residual part of the earthquake with the estimated pulse. The wavelet-based Hilbert transform is utilized to characterize the nonstationary properties, the instantaneous amplitude, and frequencies of far-field records to model residual part. The effects of near-fault and pulse are estimated based on the seismological properties of the region. The validation of the procedure is indicated by comparing simulated time-series, response spectra, and Arias intensity with recorded pulse-like records in two different earthquakes in California; the M_w 6.7 1994 Northridge and the M_w 6.5 1979 Imperial valley.

     

Magnitude reassessment of New Zealand earthquakes

The surface-wave magnitudes of a selection of New Zealand earthquakes have been determined on a consistent basis using the ‘Prague formula’ and station corrections. The earthquakes range in magnitude from about 5 to 7.8, covering the instrumental period 1901–1988. Magnitudes for many of the earlier events had not been properly determined previously; and some significant discrepancies from the traditional magnitudes were found. The use of European station data (160° < D < 175°) is important to New Zealand because of its geographical isolation. These distant data were found to give consistently slightly higher M_s than closer stations, but could be used without bias through the station correction procedure. The relationship between M_s and M_L was found for 31 ‘shallow’ New Zealand events and much of the scatter was explained as a function of depth. Significant differences in M_s/M_L expressions from Europe and California were also found. The limited New Zealand data for M_w and M₀ related well to Californian and global relationships with M_s.     

A note on the M w 6.3 earthquake in Iceland on 29 May 2008 at 15:45 UTC

This preliminary study aims to investigate a M _w 6.3 earthquake that occurred in South Iceland on Thursday 29 May 2008 at 15:45 UTC. The epicentre was in the Olfus District between the towns of Selfoss and Hveragerdi. This study examines the data recorded and the damage observed immediately after the event. Horizontal accelerations of up to 80%g were recorded in the epicentral region and there is visual evidence that the vertical acceleration exceeded 1 g. The PGA data is compared to a ground motion estimation model developed for the South Iceland earthquakes in June 2000. In general the basic properties of this event are found to be similar to the characteristics of the South Iceland earthquakes in June 2000. The duration of strong-motion is short and the intensity attenuates rapidly with increasing distance. The earthquake action resisted by buildings in the near fault area is inspected through evaluation of elastic as well as inelastic response spectra. The vast majority of structures seemed to withstand the strong-motion fairly competently and without significant visual damage due firstly to the low-rise, predominantly reinforced concrete or timber, style of buildings. Secondly, the short duration of strong-motion contributed to the endurance of structures.     

Possible precursory accelerated Benioff strain in the region of Sistan Suture Zone of Eastern Iran

We investigated whether accelerated seismic strain release precedes large earthquakes occurring in and around the Sistan Suture Zone, Eastern Iran. Online catalogs of teleseismic events occurring post-1960 within the region 27.0°–37.0°N, 55.0°–65.0°E, report five M _w > 7.0 earthquakes, namely, 1968 Dasht-e-Bayaz, 1978 Tabas, 1979 Khuli-Buniabad, 1981 Sirch and 1997 Zirkuh-e-Q’aenat events. We defined four earthquake test episodes, 1968–1978, 1978–1981, 1979–1981, and 1981–1997, with all catalogued intermediate events having magnitudes within 2.0 units that of the final large event. Using the 1968 event as the starting point, we investigated possible increased moderate earthquake activity patterns prior to the large events of 1978, 1981 and 1997 by examining if the cumulative Benioff strain released from such preceding events followed a power law time-to-failure. Our investigation seem to suggest that the 1978, 1981 and 1997 events (i) followed a period of accelerated moderate earthquake activity and (ii) the radius of their optimal critical region, R, scaled with their magnitude, M, according to the scaling law log R ∝ 0.36 M. Our suggestions conform to those proposed by similar investigations in varied seismotectonic regimes.     

Earth's rotation-triggered earthquakes before the 2018 MS5.7 Xingwen earthquake

For earthquakes (M_L≥2.0) that occurred from January 2006 to October 2018 around the M_S5.7 Xingwen earthquake occurred on 16 December 2018 in Xingwen, Sichuan province, China, we statistically investigated the correlation between the phase of Earth's rotation and the occurrence of earthquakes via Schuster's test to determine the signals that triggered earthquakes before the M_S5.7 Xingwen event. The results were evaluated based on the P-value where a smaller P-value corresponded to a higher correlation between the occurrence of an earthquake and Earth's rotation. We investigated the spatial distribution of P-values in the region around the epicenter of the M_S5.7 Xingwen event, and obtained a result exhibiting a extremely low-P-value region. The M_S5.7 event occurred inside near the northern boundary of this region. Furthermore, we analyzed the temporal evolution of P-values for earthquakes that occurred within the extremely low-P-value region and found that some extremely low P-values (less that 0.1%), i.e., significant correlation, were calculated for earthquakes that occurred before the M_S5.7 Xingwen earthquake. Among sixty-one earthquakes with the lowest P-value, occurred from May 2014 to April 2018, a vast majority of them occurred during the acceleration of Earth's rotation. The lower P-value obtained in this study reveals that the Xingwen source body probably was extremely unstable prior to the occurrence of the M_S5.7 Xingwen earthquake.     

Paleoseismic study of the Kamishiro Fault on the northern segment of the Itoigawa–Shizuoka Tectonic Line,Japan

The M_w 6.2 (Mj 6.8) Nagano (Japan) earthquake of 22 November 2014 produced a 9.3-km long surface rupture zone with a thrust-dominated displacement of up to 1.5 m, which duplicated the pre-existing Kamishiro Fault along the Itoigawa–Shizuoka Tectonic Line (ISTL), the plate-boundary between the Eurasian and North American plates, northern Nagano Prefecture, central Japan. To characterize the activity of the seismogenic fault zone, we conducted a paleoseismic study of the Kamishiro Fault. Field investigations and trench excavations revealed that seven morphogenic paleohistorical earthquakes (E2–E8) prior to the 2014 M_w 6.2 Nagano earthquake (E1) have occurred on the Kamishiro Fault during the last ca. 6000 years. Three of these events (E2–E4) are well constrained and correspond to historical earthquakes occurring in the last ca. 1200 years. This suggests an average recurrence interval of ca. 300–400 years on the seismogenic fault of the 2014 Kamishiro earthquake in the past 1200 years. The most recent event prior to the 2014 earthquakes (E1) is E2 and the penultimate and antepenultimate faulting events are E3 and E4, respectively. The penultimate faulting event (E3) occurred during the period of AD 1800–1400 and is associated with the 1791 M_w 6.8 earthquake. The antepenultimate faulting event (E4) is inferred to have occurred during the period of ca. AD 1000–700, likely corresponding to the AD 841 M_w 6.5 earthquake. The oldest faulting event (E8) in the study area is thought to have occurred during the period of ca. 5600–6000 years. The throw rate during the early Holocene is estimated to be 1.2–3.3 mm/a (average, 2.2 mm/a) with an average amount of characteristic offset of 0.7–1.1 m produced by individual event. When compared with active intraplate faults on Honshu Island, Japan, these slip rates and recurrence interval estimated for morphogenic earthquakes on the Kamishiro Fault along the ISTL appear high and short, respectively. This indicates that present activity on this fault is closely related to seismic faulting along the plate boundary between the Eurasian and North American plates.     

Ground motion prediction in Beirut: a multi-step procedure coupling empirical Green’s functions,ground motion prediction equations and instrumental transfer functions

The strike slip Yammouneh fault is the longest fault in Lebanon, crossing the territory from South to North. It was responsible for major historical earthquakes like the 1202 A.D. earthquake, estimated to M_s7.6. This paper presents a site-specific estimation of the ground motion caused by a potential M_w7.5 earthquake on the Yammouneh fault, similar to the 1202 event, for various sites within the Beirut area. The empirical Green’s function technique EGF is used to estimate the median and the standard deviations of the seismic ground motion at the reference station BHL, taking into account epistemic and aleatory uncertainties related to source parameters. These uncertainties were quantified through a sensitivity analysis of the position of the rupture nucleation X_nuc, the slip roughness parameter K, the corner frequency f_c and the magnitude M_c of the EGF. The rock ground motion is then transferred to various other sites within the Beirut area, using instrumental Fourier transfer functions. Site amplification factors are next deduced by computing the ratio between response spectra at sediment sites and at a reference rock station. Considering the limits of the EGF method in the near field of extended sources, the EGF approach is considered only up to a magnitude M_w of 6.5. Selected Ground Motion Predictive Equations are then used to simulate a M_w7.5 event at a reference station. By applying the amplification factors, the response spectra at the different sites of Beirut are also calculated and compared with the actual response spectra used in the Lebanese regulations.     

Quantification of major earthquake tsunamis of the Japan Sea

The size of major tsunamigenic earthquakes which occurred in the Japan Sea is quantified on the basis of seismic and tsunamigenic source parameters. The tsunami magnitude M_t is determined from the instrumental tsunami-wave amplitudes. The M_t values thus obtained are on average 0.2 units larger than the values of moment magnitude M_w, though the M_t scale has originally been adjusted to agree with M_w. Moreover, the volume of displaced water at the source is on average 2.3 times as large as that for the Pacific events with a comparable M_w. Nevertheless, the observed height of the sea-level disturbance at the source is found consistent with the amount of crustal deformation computed for the seismic fault models. These results indicate that the tsunami source potential itself is large for M_w in comparison with the Pacific events. The large source potential is explained in terms of the effective difference both in the rigidity of the source medium and in the geometry of the fault motion. For the Japan Sea events, the M_t scale still provides the physical measure of the tsunami potential, and M_t minus 0.2 corresponds to M_w. This predicts that the maximum amplitude of tsunami waves from Japan Sea earthquakes is at least two times as large as that from Pacific earthquakes with a comparable M_w.     

Correlating Precursory Declines in Groundwater Radon with Earthquake Magnitude

Both studies at the Antung hot spring in eastern Taiwan and at the Paihe spring in southern Taiwan confirm that groundwater radon can be a consistent tracer for strain changes in the crust preceding an earthquake when observed in a low‐porosity fractured aquifer surrounded by a ductile formation. Recurrent anomalous declines in groundwater radon were observed at the Antung D1 monitoring well in eastern Taiwan prior to the five earthquakes of magnitude (M_w): 6.8, 6.1, 5.9, 5.4, and 5.0 that occurred on December 10, 2003; April 1, 2006; April 15, 2006; February 17, 2008; and July 12, 2011, respectively. For earthquakes occurring on the longitudinal valley fault in eastern Taiwan, the observed radon minima decrease as the earthquake magnitude increases. The above correlation has been proven to be useful for early warning local large earthquakes. In southern Taiwan, radon anomalous declines prior to the 2010 M_w 6.3 Jiasian, 2012 M_w 5.9 Wutai, and 2012 M_L 5.4 Kaohsiung earthquakes were also recorded at the Paihe spring. For earthquakes occurring on different faults in southern Taiwan, the correlation between the observed radon minima and the earthquake magnitude is not yet possible.     

The instability of the <Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript> and <Emphasis Type="Italic">M</Emphasis><Subscript>L</Subscript> comparison for earthquakes in Greece for the period 1969 to 2007

We use 576 earthquakes of magnitude, M _w, 3.3 to 6.8 that occurred within the region 33° N–42.5° N, 19° E–30° E in the time period 1969 to 2007 to investigate the stability of the relation between moment magnitude, M _w, and local magnitude, M _L, for earthquakes in Greece and the surrounding regions. We compare M _w to M _L as reported in the monthly bulletins of the National Observatory of Athens (NOA) and to M _L as reported in the bulletins of the Seismological Station of the Aristotle University of Thessaloniki. All earthquakes have been analyzed through regional or teleseismic waveform inversion, to obtain M _w, and have measured maximum trace amplitudes on the Wood–Anderson seismograph in Athens, which has been in operation since 1964. We show that the Athens Wood–Anderson seismograph performance has changed through time, affecting the computed by NOA M _L by at least 0.1 magnitude units. Specifically, since the beginning of 1996, its east–west component has been recording systematically much larger amplitudes compared to the north–south component. From the comparison between M _w and M _L reported by Thessaloniki, we also show that the performance of the sensors has changed several times through time, affecting the calculated M _L’s. We propose scaling relations to convert the M _L values reported from the two centers to M _w. The procedures followed here can be applied to other regions as well to examine the stability of magnitude calculations through time.     

Evidence for fault-related directionality and localized site effects from strong motion recordings of the 2003 Boumerdes (Algeria) earthquake: Consequences on damage distribution and the Algerian seismic code

The Algiers–Boumerdes region has been struck by a destructive magnitude 6.8 (M_w) earthquake on May 21, 2003. The study presented in this paper is based on main shock strong motions from 13 stations of the Algerian accelerograph network. A maximum 0.58g peak ground acceleration (PGA) has been recorded at 20 km from the epicenter, only about 150 m away from a PGA of 0.34g, with both a central frequency around 5 Hz, explained by a strong very localized site effect, confirmed by receiver function technique results showing peaks at 5 Hz with amplitudes changing by a factor of 2. Soil amplifications are also evidenced at stations located in the quaternary Mitidja basin, explaining the higher PGA values recorded at these stations than at stations located on firm soil at similar distances from the epicenter. A fault-related directionality effect observed on the strong motion records and confirmed by the study of the seismic movement anisotropy, in agreement with the N65 fault plan direction, explains the SW–NE orientation of the main damage zone. In the near field, strong motions present a high-frequency content starting at 3 Hz with a central frequency around 8 Hz, while in the far field their central frequency is around 3 Hz, explaining the high level of damage in the 3- to 4-story buildings in the epicentral zone. The design spectra overestimate the recorded mean response spectra, and its high corner frequency is less than the recorded one, leading to a re-examination of the seismic design code that should definitively integrate site-related coefficient, to account for the up to now neglected site amplification, as well as a re-modeling of the actual design spectra. Finally, both the proposed Algerian attenuation law and the worldwide laws usually used in Algeria underestimate the recorded accelerations of the 6.8 (M_w) Boumerdes earthquake, clearly showing that it is not possible to extrapolate the proposed Algerian law to major earthquakes.     

Large-scale aseismic creep in the areas of the strong earthquakes revealed from the GRACE data on the time variations of the Earth’s gravity field

The refinement of the accuracy and resolution of the monthly global gravity field models from the GRACE satellite mission, together with the accumulation of more than a decade-long series of these models, enabled us to reveal the processes that occur in the regions of large (M_w≥8) earthquakes that have not been studied previously. The previous research into the time variations of the gravity field in the regions of the giant earthquakes, such as the seismic catastrophes in Sumatra (2004) and Chile (2010), and the Tohoku mega earthquake in Japan (2011), covered the coseismic gravity jump followed by the long postseismic changes reaching almost the same amplitude. The coseismic gravity jumps resulting from the lower-magnitude events are almost unnoticeable. However, we have established a long steady growth of gravity anomalies after a number of such earthquakes. For instance, in the regions of the subduction earthquakes, the growth of the positive gravity anomaly above the oceanic trench was revealed after two events with magnitudes M_w=8.5 in the Sumatra region (the Nias earthquake of March 2005 and the Bengkulu event of September 2007 near the southern termination of Sumatra Island), after the earthquake with M_w=8.5 on Hokkaido in September 2007, a doublet Simushir earthquake with the magnitudes M_w = 8.3 and 8.1 in the Kuriles in November 2006 and January 2007, and after the earthquake off the Samoa Island in September 2009 (M_w=8.1). The steady changes in the gravity field have also been recorded after the earthquake in the Sichuan region (May 2008, M_w = 8.0) and after the doublet event with magnitudes 8.6 and 8.2, which occurred in the Wharton Basin of the Indian Ocean on April 11, 2012. The detailed analysis of the growth of the positive anomaly in gravity after the Simushir earthquake of November 2006 is presented. The growth started a few months after the event synchronously with the seismic activation on the downdip extension of the coseismically ruptured fault plane zone. The data demonstrating the increasing depth of the aftershocks since March 2007 and the approximately simultaneous change in the direction and average velocity of the horizontal surface displacements at the sites of the regional GPS network indicate that this earthquake induced postseismic displacements in a huge area extending to depths below 100 km. The total displacement since the beginning of the growth of the gravity anomaly up to July 2012 is estimated at 3.0 m in the upper part of the plate’s contact and 1.5 m in the lower part up to a depth of 100 km. With allowance for the size of the region captured by the deformations, the released total energy is equivalent to the earthquake with the magnitude M_w = 8.5. In our opinion, the growth of the gravity anomaly in these regions indicates a large-scale aseismic creep over the areas much more extensive than the source zone of the earthquake. These processes have not been previously revealed by the ground-based techniques. Hence, the time series of the GRACE gravity models are an important source of the new data about the locations and evolution of the locked segments of the subduction zones and their seismic potential.     

Harmonization check of M w within the central, northern, and northwestern European earthquake catalogue (CENEC)

Large data sets covering large areas and time spans and composed of many different independent sources raise the question of the obtained degree of harmonization. The present study is an analysis of the harmonization with respect to the moment magnitude M _w within the earthquake catalogue for central, northern, and northwestern Europe (CENEC). The CENEC earthquake catalogue (Grünthal et al., J Seismol, 2009) contains parameters for over 8,000 events in the time period 1000–2004 with magnitude M _w ≥ 3.5. Only about 2% of the data used for CENEC have original M _w magnitudes derived directly from digital data. Some of the local catalogues and data files providing data give M _w, but calculated by the respective agency from other magnitude measures or intensity. About 60% of the local data give strength measures other than M _w, and these have to be transformed by us using available formulae or new regressions based on original M _w data. Although all events are thus unified to M _w magnitude, inhomogeneity in the M _w obtained from over 40 local catalogues and data files and 50 special studies is inevitable. Two different approaches have been followed to investigate the compatibility of the different M _w sets throughout CENEC. The first harmonization check is performed using M _w from moment tensor solutions from SMTS and Pondrelli et al. (Phys Earth Planet Inter 130:71–101, 2002; Phys Earth Planet Inter 164:90–112, 2007). The method to derive the SMTS is described, e.g., by Braunmiller et al. (Tectonophysics 356:5–22, 2002) and Bernardi et al. (Geophys J Int 157:703–716, 2004), and the data are available in greater extent since 1997. One check is made against the M _w given in national catalogues and another against the M _w derived by applying different empirical relations developed for CENEC. The second harmonization check concerns the vast majority of data in CENEC related to earthquakes prior to 1997 or where no moment tensor based M _w exists. In this case, an empirical relation for the M _w dependence on epicentral intensity (I ₀) and focal depth (h) was derived for 41 master events, i.e., earthquakes, located all over central Europe, with high-quality data. To include also the data lacking h, the corresponding depth-independent relation for these 41 events was also derived. These equations are compared with the different sets of data from which CENEC has been composed, and the goodness of fit is demonstrated for each set. The vast majority of the events are very well or reasonably consistent with the respective relation so that the data can be said to be harmonized with respect to M _w, but there are exceptions, which are discussed in detail.     

A Preliminary Analysis of the Earthquake Disaster of Buildings during Two Destructive Earthquakes in Xinjiang

This paper introduces the characteristics of earthquake disasters in the 2013 Urumqi M_S5.1 and the 2015 Pishan M_S6.5 earthquakes, which are directly beneath the cities. Based on the discussion on regional tectonic background, site conditions, seismic fortification level of buildings, we preliminarily analyzed the causes of earthquake disaster for buildings caused by the two earthquakes.Finally,we give some advice on earthquake resistance for residential buildings and earthquake damage prevention in Xinjiang.     

Evidence for a Tang-Song Dynasty great earthquake along the Longmen Shan Thrust Belt prior to the 2008 <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> 7.9 Wenchuan earthquake,China

The magnitude (M _w) 7.9 Wenchuan earthquake occurred on 12 May 2008 in the Longmen Shan region of China, the transition zone between the Tibetan Plateau and the Sichuan Basin, resulting in widespread damage throughout central and western China. The steep, high-relief eastern margin of the Tibetan Plateau has undergone rapid Cenozoic uplift and denudation accompanied by folding and thrusting, yet no large thrust earthquakes are known prior to the 2008 M _w 7.9 Wenchuan earthquake. Field and excavation investigations reveal that a great historical earthquake occurred in the Sichuan region that ruptured a >200-km-long thrust fault within the Longmen Shan Thrust Belt, China, which also triggered the 2008 M _w 7.9 Wenchuan earthquake. The average co-seismic slip amount produced by this historical earthquake is estimated to be 2–3 m, comparable with that caused by the 2008 Wenchuan earthquake. Paleoseismic and archaeological evidence and radiocarbon dating results show that the penultimate great earthquake occurred in the Sichuan region during the late Tang-Song Dynasty, between AD 800 and 1000, suggesting a recurrence interval of ~1,000–1,200 years for Wenchuan-magnitude (M = ~8) earthquakes in the late Holocene within the Longmen Shan Thrust Belt. This finding is in contrast with previous estimates of 2,000–10,000 years for the recurrence interval of large earthquakes within the Longmen Shan Thrust Belt, as obtained from long-term slip rates based on the Global Positioning System and geological data, thereby necessitating substantial modifications to existing seismic-hazard models for the densely populated region at the eastern marginal zone of the Tibetan Plateau.     

Critical evaluation of strong motion in Kocaeli and Düzce (Turkey) earthquakes

Turkey was struck by two major events on August 17th and November 12th, 1999. Named Kocaeli (M_w=7.4) and Düzce (M_w=7.2) earthquakes, respectively, the two earthquakes provided the most extensive strong ground motion data set ever recorded in Turkey. The strong motion stations operated by the General Directorate of Disaster Affairs, the Kandilli Observatory and Earthquake Research Institute of Bogazici University and Istanbul Technical University have produced at least 27 strong motion records for the Kocaeli earthquake within 200 km of the fault. Kocaeli earthquake has generated six motions within 20 km of the fault adding significantly to the near-field database of ground motions for M_w>=7.0 strike–slip earthquakes. The paper discusses available strong motion data, studies their attenuation characteristics, analyses time domain, as well as spectral properties such as spectral accelerations with special emphasis on fault normal and fault parallel components and the elastic attenuation parameter, kappa. A simulation of the Kocaeli earthquake using code FINSIM is also presented.