Transfrontier macroseismic observations of the Ml = 5.4 earthquake of February 22, 2003 at Rambervillers,France

The 2003 Ml = 5.4 Rambervillers earthquake, north-east of France, is the largest seismic event recorded north of the Alps since the 1992 Ms = 5.3, I₀ = VII, Roermond earthquake, Netherlands. With a maximum macroseismic intensity of VI-VII EMS-98, the 2003 event was broadly felt to a distance of 300 km from the epicentre. It provides a unique opportunity to test and compare the different procedures used in France, Germany and Switzerland when evaluating macroseismic intensities. The main purpose of this paper is to present a common transfrontier macroseismic map based on the EMS-98 intensity scale. Maximum horizontal accelerations recorded in the area are compared to the intensity values, and we propose to use a differential technique to re-estimate the magnitude of the 1682 Remiremont, I₀ = VIII, earthquake, which occurred 40 km south of Rambervillers.     

Linking ground motion measurements and macroseismic observations in France: a case study based on accelerometric and macroseismic databases

Comparison between accelerometric and macroseismic observations is made for three M _w?=?4.5 earthquakes, which occurred in north-eastern France and south-western Germany in 2003 and 2004. Scalar and spectral instrumental parameters are processed from the accelerometric data recorded by nine accelerometric stations located between 29 and 180 km from the epicentres. Macroseismic data are based on French Internet reports. In addition to the single questionnaire intensity, analysis of the internal correlation between the encoded answers highlights four predominant fields of questions bearing different physical meanings: (1) “vibratory motions of small objects”, (2) “displacement and fall of objects”, (3) “acoustic noise” and (4) “personal feelings”. Best correlations between macroseismic and instrumental observations are obtained when the macroseismic parameters are averaged over 10-km-radius circles around each station. Macroseismic intensities predicted by published peak ground velocity (PGV)–intensity relationships agree with our observed intensities, contrary to those based on peak ground acceleration (PGA). Correlation between the macroseismic and instrumental data for intensities between II and V (EMS-98) is better for PGV than for PGA. Correlation with the response spectra exhibits clear frequency dependence for all macroseismic parameters. Horizontal and vertical components are significantly correlated with the macroseismic parameters between 1 and 10 Hz, a range corresponding to both natural frequencies of most buildings and high energy content in the seismic ground motion. Between 10 and 25 Hz, a clear lack of correlation between macroseismic and instrumental observations exists. It could be due to a combination of the decrease in the energy signal above 10 Hz, a high level of anthropogenic noise and an increase in variability in soil conditions. Above 25 Hz, the correlation coefficients between the acceleration response spectra and the macroseismic parameters are close to the PGA correlation level.     

Macroseismic field of the October 27, 2004 Vrancea (Romania) moderate subcrustal earthquake

On October 27, 2004, a moderate size earthquake occurred in the Vrancea seismogenic region (Romania). The Vrancea seismic zone is an area of concentrated seismicity at intermediate depths beneath the bending area of the southeastern Carpathians. The 2004 M _w?=?6 Vrancea subcrustal earthquake is the largest seismic event recorded in Romania since the 1990 earthquakes. With a maximum macroseismic intensity of VII Medvedev–Sponheuer–Kárník (MSK-64) scale, the seismic event was felt to a distance of 600 km from the epicentre. This earthquake caused no serious damage and human injuries. The main purpose of this paper is to present the macroseismic map of the earthquake based on the MSK-64 intensity scale. After the evaluation of the macroseismic effects of this earthquake, an intensity dataset has been obtained for 475 sites in the Romanian territory. Also, the maximum horizontal accelerations recorded in the area by the K2 network are compared to the intensity values.     

Evidence of nonlinear site response in horizontal-to-vertical spectral ratio from near-field earthquakes

The issue is addressed as to whether the horizontal-to-vertical spectral ratio (HVSR) method is sensitive to the amplitude of ground motion from near-field earthquakes. Twenty-one three-component accelerograms from two closely located similar soil sites in the town of Lefkas are used. The recordings represent 17 earthquakes covering a wide range of magnitudes, epicentral distances and azimuths. Peak horizontal accelerations (PGA) and velocities (PGV) lie in the ranges 20–540 cm/s² and 1.4–55.2 cm/s. For each HVS ratio, the site's fundamental-resonance frequency, f_res, is determined visually. Linear correlation analysis shows that f_res is strongly (negatively) correlated to PGA and PGV (r between −0.7 and −0.8); no correlation is found with resonance amplitude or epicentral distance. We show that the observed correlation is attributable to soil nonlinearity and indicate how weak-motion estimates of f_res can be corrected for use in assessing site response during strong shaking.     

Macroseismic estimates of magnitude in Italy

A multi-parametric study of empirical relationships between macroseismic data and magnitude is presented for the Italian region by the analysis of a new extended data set concerning 146 earthquakes. The available magnitude determinations include all of the most intense earthquakes which occurred in Italy in the last century and have been obtained by an accurate revision of original instrumental data. Intensity data have been revised and upgraded on the basis of the most recent studies: only local intensities directly documented have been used. Macroseismic determinations ofM _s ,m _B andM _wa magnitudes have been performed. The empirical relationships between maximum felt intensity (I _max ) and magnitude have been determined by the use of a distribution-free approach and a linear regression analysis. This last parameterization allows for the explanation of more than 60% of the variation in magnitude. In order to improve these results, the linear dependence between magnitude,I _max and average distances (in logarithm) corresponding to fixed attenuation values has been explored. The comparison between instrumental magnitudes and corresponding macroseismic estimates obtained from empirical relationships shows that the respective uncertainties are comparable.     

Intensity attenuation in the Campania region,Southern Italy

Ground motion prediction equations (GMPE) in terms of macroseismic intensity are a prerequisite for intensity-based shake maps and seismic hazard assessment and have the advantage of direct relation to earthquake damage and good data availability also for historical events. In this study, we derive GMPE for macroseismic intensity for the Campania region in southern Italy. This region is highly exposed to the seismic hazard related to the high seismicity with moderate- to large-magnitude earthquakes in the Appenninic belt. The relations are based on physical considerations and are easy to implement for the user. The uncertainties in earthquake source parameters are accounted for through a Monte Carlo approach and results are compared to those obtained through a standard regression scheme. One relation takes into account the finite dimensions of the fault plane and describes the site intensity as a function of Joyner–Boore distance. Additionally, a relation describing the intensity as a function of epicentral distance is derived for implementation in cases where the dimensions of the fault plane are unknown. The relations are based on an extensive dataset of macroseismic intensities for large earthquakes in the Campania region and are valid in the magnitude range M _w = 6.3–7.0 for shallow crustal earthquakes. Results indicate that the uncertainties in earthquake source parameters are negligible in comparison to the spread in the intensity data. The GMPE provide a good overall fit to historical earthquakes in the region and can provide the intensities for a future earthquake within 1 intensity unit.     

Estimation of Strong Ground Motions for 2001 Bhuj (M w 7.6), India Earthquake

The strong ground motions for the 2001 Bhuj (M _w 7.6) India earthquake have been estimated on hard rock and B/C boundary (NEHRP) levels using a recently modified version of stochastic finite fault modeling based on dynamic corner frequency (Motazedian and Atkinson in Bull Seismol Soc Am 95, 995–1010 2005). Incorporation of dynamic corner frequency removes the limitations of earlier stochastic methods. Simulations were carried out at 13 sites in Gujarat where structural response recorder (SRR) recordings are available. In addition, accelerograms were simulated at the B/C boundary at a large number of points distributed on a grid. The corresponding response spectra have also been estimated. The values of peak ground accelerations and spectral accelerations at three periods (0.4, 0.75 and 1.25 s) are presented in the form of contour maps. The maximum value of peak ground acceleration (PGA) in the center of meizoseismal zone is 550 cm/s². The response spectral acceleration in same zone is 900 cm/s² (T = 0.4 s), 600 cm/s² (T = 0.75 s) and 300 cm/s² (T = 1.25 s). The innermost PGA contour is on the fault plane. A comparison of the PGA values obtained at 13 sites in this study with those obtained in earlier studies on the same sites, but employing different methods, show that the present PGA values are comparable at most of the sites. The rate of decay of PGA values is fast at short distances as compared to that at longer distances. The PGA values obtained here put some constraints on the expected values from a similar earthquake in the region. A synthetic intensity map has been prepared from the estimated values of PGA using an empirical relation. A comparison with the reported intensity map of the earthquake shows the synthetic MMI values, as expected, are lower by 1 unit compared to reported intensity map. The contour map of PGA along with the contour maps of spectral acceleration at various periods permit the assessment of damage potential to various categories of houses and other structures. Such information will be quite important in planning of mitigation and disaster management programs in the region.     

Macroseismic evaluation of the impact of the Olyutorskii earthquake and its larger aftershocks

This paper presents the results from a macroseismic survey of the impact and consequences of the M _W = 7.6 April 20(21), 2006 Olyutorskii earthquake in the area of the Koryak Autonomous Okrug and the adjacent areas in Kamchatka and Magadan regions. The earthquake was felt over an area of about 400000 sq. km with intensities of II to IX–X on the MSK-64 scale. Information was gathered from 37 population centers situated in this area and was used to present a summary of felt effects, to construct an isoseismals map, and to determine the macroseismic magnitude.     

Probabilistic seismic hazard assessment in SE-Spain based on macroseismic site histories

In many countries such as Spain earthquake databases still mainly comprise macroseismic data from felt effects. The full exploit of this information is of basic importance for seismic risk assessment and emergency planning, given the strict link between macroseismic intensity and damage. A probabilistic procedure specifically developed to handle macroseismic data, mostly relying on site information and seismogenic-source free, has been applied to evaluate seismic hazard in SE-Spain (Alicante-Murcia region). Present seismicity is moderate-low with largest magnitudes slightly over Mw5.0. The historical record includes very destructive earthquakes, maximum EMS98 intensities reaching IX–X and X in the nineteenth century (e.g., Torrevieja 1829 earthquake). Very recently, two events in the area on 11 May 2011 (Mw4.5, Mw5.2) killed nine people, injured 300, and produced important damage in the city of Lorca. Regional hazard maps for the area together with specific hazard curves at selected localities are obtained. Results are compared with the maximum observed intensities in the period 1300–2012, and with the values in the seismic hazard map from the Spanish Building Code in force. In general, the maximum felt intensity values are closer to the hazard values calculated for 2 % probability of exceedance in 50 years, using felt and expected intensity. The intensity-based probabilistic hazard maps obtained through the applied approach reduce the inherent smoothing of those based on standard probabilistic seismic hazard assessment approaches for the region, allowing identifying possible over- or sub-estimates of site hazard values, providing very valuable information for risk reduction strategies or for future updates of the building code hazard maps.     

Field survey around strong motion stations and its implications on the seismic intensity in the Lushan earthquake on April 20, 2013

The M _s7.0 Lushan earthquake on April 20, 2013 is another destructive event in China since the M _s8.0 Wenchuan earthquake in 2008 and M _s7.1 Yushu earthquake in 2010. A large number of strong motion recordings were accumulated by the National Strong Motion Observation Network System of China. The maximum peak ground acceleration (PGA) at Station 51BXD in Baoxing Country is recorded as ?1,005.3 cm/s², which is even larger than the maximum one in the Wenchuan earthquake. A field survey around three typical strong motion stations confirms that the earthquake damage is consistent with the issued map of macroseismic intensity. For the oscillation period 0.3–1.0 s which is the common natural period range of the Chinese civil building, a comparison shows that the observed response spectrums are considerably smaller than the designed values in the Chinese code and this could be one of the reasons that the macroseismic intensity is lower than what we expected despite the high amplitude of PGAs. The Housner spectral intensities from 16 stations are also basically correlated with their macroseismic intensities, and the empirical distribution of spectral intensities from Lushan and Wenchuan Earthquakes under the Chinese scale is almost identical with those under the European scale.     

Historical intermediate-depth earthquakes in the southern Aegean Sea Benioff zone: modeling their anomalous macroseismic patterns with stochastic ground-motion simulations

We model the macroseismic damage distribution of four important intermediate-depth earthquakes of the southern Aegean Sea subduction zone, namely the destructive 1926 M?=?7.7 Rhodes and 1935 M?=?6.9 Crete earthquakes, the unique 1956 M?=?6.9 Amorgos aftershock (recently proposed to be triggered by a shallow event), and the more recent 2002 M?=?5.9 Milos earthquake, which all exhibit spatially anomalous macroseismic patterns. Macroseismic data for these events are collected from published macroseismic databases and compared with the spatial distribution of seismic motions obtained from stochastic simulation, converted to macroseismic intensity (Modified Mercalli scale, I_MM). For this conversion, we present an updated correlation between macroseismic intensities and peak measures of seismic motions (PGA and PGV) for the intermediate-depth earthquakes of the southern Aegean Sea. Input model parameters for the simulations, such as fault dimensions, stress parameters, and attenuation parameters (e.g. back-arc/along anelastic attenuation) are adopted from previous work performed in the area. Site-effects on the observed seismic motions are approximated using generic transfer functions proposed for the broader Aegean Sea area on the basis of V_S30 values from topographic slope proxies. The results are in very good agreement with the observed anomalous damage patterns, for which the largest intensities are often observed at distances >?100 km from the earthquake epicenters. We also consider two additional “prediction” but realistic intermediate-depth earthquake scenarios, and model their macroseismic distributions, to assess their expected damage impact in the broader southern Aegean area. The results suggest that intermediate-depth events, especially north of central Crete, have a prominent effect on a wide area of the outer Hellenic arc, with a very important impact on modern urban centers along northern Crete coasts (e.g. city of Heraklion), in excellent agreement with the available historical information.     

Probability of Peak Ground Horizontal and Peak Ground Vertical Accelerations at Tehran and Surrounding Areas

The active faults near Tehran are capable of producing Mw magnitudes of 6.62?C7.23; at epicentral distances of 25?km from the active faults, and Mw?=?7.23, the peak ground horizontal acceleration, PGH, is between 386 and 730?cm/s² and peak ground vertical acceleration, PGV, is between 192 and 261?cm/s²; the historic earthquake of the fourth century BC, Mw?=?7.16 produced the highest estimated PGH acceleration in Tehran, between 254 and 479?cm/s² and PGV acceleration between 127 and 173?cm/s². Earthquakes from 1909 to 2008, within 300?km from Tehran, are used for calculation of magnitude frequency relation, and results are applied to estimate PGH and PGV accelerations as a function of magnitudes for a set of fixed epicentral distance and site conditions; also as a function of epicentral distances for a set of fixed magnitudes and site conditions. Poisson??s distribution is used for probability calculation of PGH and PGV accelerations for several exposure times, site conditions and epicentral distances; accelerations with very high probability, near 1, are presented. At an epicentral distance of 10?km and exposure time of 450?years, in the northern part of Tehran, close to Mosha and the Northern Tehran faults, PGH acceleration is 800?C420?cm/s² and PGV is 400?C220?cm/s² with high probability. At an epicentral distance of 25?km and 1,000?years exposure time, PGH is 610?C320?cm/s² and PGV is 310?C160?cm/s² with high probability, where larger values are for soft soil and smaller values are for hard rock.     

Predicted ground motion after the L��Aquila 2009 earthquake (Italy, M w 6.3): input spectra for seismic microzoning

After the April 6th 2009 L’Aquila earthquake (M _w 6.3), where 306 people died and a further 60,000 were displaced, seismic microzoning investigations have been carried out for towns affected by a macroseismic intensity equal to or greater than 7 MCS. Based upon seismotectonic data, historical seismicity and strong motion records, we defined input spectra to be used in the numerical simulations of seismic microzoning in four key municipalities, including the town of L’Aquila. We adopted two main approaches: uniform hazard response spectra are obtained by a probabilistic seismic hazard assessment introducing some time-dependency for individual faults on the study area; a deterministic design spectrum is computed from magnitude/distance pairs extracted by a stationary probabilistic analysis of historical intensities. The uniform hazard spectrum of the present Italian building code represents the third, less restrictive, response spectrum to be used for the numerical simulations in seismic microzoning. Strong motions recordings of the main shock of the L’Aquila sequence enlighten the critical role played by both the local response and distances metric for sites located above a seismogenic fault; however, these time-histories are compatible with the uncertainties of a deterministic utilization of ground motion predictive equations. As recordings at very near field are rare, they cannot be neglected while defining the seismic input. Disaggregation on the non-Possonian seismotectonic analysis and on the stationary site-intensity estimates reach very similar results in magnitude-distance pairs identification; we interpret this convergence as a validation of the geology-based model by historical observations.     

Exploration tree approach to estimate historical earthquakes Mw and depth,test cases from the French past seismicity

The estimation of the seismological parameters of historical earthquakes is a key step when performing seismic hazard assessment in moderate seismicity regions as France. We propose an original method to assess magnitude and depth of historical earthquakes using intensity data points. A flowchart based on an exploration tree (ET) approach allows to apply a consistent methodology to all the different configurations of the earthquake macroseismic field and to explore the inherent uncertainties. The method is applied to French test case historical earthquakes, using the SisFrance (BRGM, IRSN, EDF) macroseismic database and the intensity prediction equations (IPEs) calibrated in the companion paper (Baumont et al. Bull Earthq Eng, 2017). A weighted least square scheme allowing for the joint inversion of magnitude and depth is applied to earthquakes that exhibit a decay of intensity with distance. Two cases are distinguished: (1) a “Complete ET” is applied to earthquakes located within the metropolitan territory, while (2) a “Simplified ET” is applied to both, offshore and cross border events, lacking information at short distances but disposing of reliable data at large ones. Finally, a priori-depth-based magnitude computation is applied to ancient or poorly documented events, only described by single/sporadic intensity data or few macroseismic testimonies. Specific processing of “felt” testimonies allows exploiting this complementary information for poorly described earthquakes. Uncertainties associated to magnitude and depth estimates result from both, full propagation of uncertainties related to the original macroseismic information and the epistemic uncertainty related to the IPEs selection procedure.     

Seismic hazard assessment of the Province of Murcia (SE Spain): analysis of source contribution to hazard

A probabilistic seismic hazard assessment of the Province of Murcia in terms of peak ground acceleration (PGA) and spectral accelerations [SA(T)] is presented in this paper. In contrast to most of the previous studies in the region, which were performed for PGA making use of intensity-to-PGA relationships, hazard is here calculated in terms of magnitude and using European spectral ground-motion models. Moreover, we have considered the most important faults in the region as specific seismic sources, and also comprehensively reviewed the earthquake catalogue. Hazard calculations are performed following the Probabilistic Seismic Hazard Assessment (PSHA) methodology using a logic tree, which accounts for three different seismic source zonings and three different ground-motion models. Hazard maps in terms of PGA and SA(0.1, 0.2, 0.5, 1.0 and 2.0 s) and coefficient of variation (COV) for the 475-year return period are shown. Subsequent analysis is focused on three sites of the province, namely, the cities of Murcia, Lorca and Cartagena, which are important industrial and tourism centres. Results at these sites have been analysed to evaluate the influence of the different input options. The most important factor affecting the results is the choice of the attenuation relationship, whereas the influence of the selected seismic source zonings appears strongly site dependant. Finally, we have performed an analysis of source contribution to hazard at each of these cities to provide preliminary guidance in devising specific risk scenarios. We have found that local source zones control the hazard for PGA and SA(T ≤ 1.0 s), although contribution from specific fault sources and long-distance north Algerian sources becomes significant from SA(0.5 s) onwards.     

Development of damage probability matrices based on Greek earthquake damage data

A comprehensive study is presented for empirical seismic vulnerability assessment of typical structural types, representative of the building stock of Southern Europe, based on a large set of damage statistics. The observational database was obtained from post-earthquake surveys carried out in the area struck by the September 7, 1999 Athens earthquake. After analysis of the collected observational data, a unified damage database has been created which comprises 180,945 damaged buildings from/after the near-field area of the earthquake. The damaged buildings are classified in specific structural types, according to the materials, seismic codes and construction techniques in Southern Europe. The seismic demand is described in terms of both the regional macroseismic intensity and the ratio αg/ao, where αg is the maximum peak ground acceleration (PGA) of the earthquake event and ao is the unique value PGA that characterizes each municipality shown on the Greek hazard map. The relative and cumulative frequencies of the different damage states for each structural type and each intensity level are computed in terms of damage ratio. Damage probability matrices (DPMs) and vulnerability curves are obtained for specific structural types. A comparison analysis is fulfilled between the produced and the existing vulnerability models.     

Selection and ranking of ground motion models for seismic hazard analysis in the Pyrenees

The issue addressed in this paper is the objective selection of appropriate ground motion models for seismic hazard assessment in the Pyrenees. The method of Scherbaum et al. (2004a) is applied in order to rank eight published ground motion models relevant to intraplate or to low deformation rate contexts. This method is based on a transparent and data-driven process which quantifies the model fit and also measures how well the underlying model assumptions are met. The method is applied to 15 accelerometric records obtained in the Pyrenees for events of local magnitude between 4.8 and 5.1, corresponding to moment magnitudes ranging from 3.7 to 3.9. Only stations at rock sites are considered. A total of 720 spectral amplitudes are used to rank the selected ground motion models. Some control parameters of these models, such as magnitude and distance definitions, may vary from one model to the other. It is thus important to correct the selected models for their difference with respect to the magnitude and distance definitions used for the Pyrenean data. Our analysis shows that, with these corrections, some of the ground motion models successfully fit the data. These are the Lussou et al. (2001) and the Berge-Thierry et al. (2003) models. According to the selected ground motion models, a possible scenario of a magnitude 6 event is proposed; it predicts response spectra accelerations of 0.08–0.1 g at 1 Hz at a hypocentral distance of 10 km.     

An ordered probit model for seismic intensity data

Seismic intensity, measured through the Mercalli–Cancani–Sieberg (MCS) scale, provides an assessment of ground shaking level deduced from building damages, any natural environment changes and from any observed effects or feelings. Generally, moving away from the earthquake epicentre, the effects are lower but intensities may vary in space, as there could be areas that amplify or reduce the shaking depending on the earthquake source geometry, geological features and local factors. Currently, the Istituto Nazionale di Geofisica e Vulcanologia analyzes, for each seismic event, intensity data collected through the online macroseismic questionnaire available at the web-page www.haisentitoilterremoto.it. Questionnaire responses are aggregated at the municipality level and analyzed to obtain an intensity defined on an ordinal categorical scale. The main aim of this work is to model macroseismic attenuation and obtain an intensity prediction equation which describes the decay of macroseismic intensity as a function of the magnitude and distance from the hypocentre. To do this we employ an ordered probit model, assuming that the intensity response variable is related through the link probit function to some predictors. Differently from what it is commonly done in the macroseismic literature, this approach takes properly into account the qualitative and ordinal nature of the macroseismic intensity as defined on the MCS scale. Using Markov chain Monte Carlo methods, we estimate the posterior probability of the intensity at each site. Moreover, by comparing observed and estimated intensities we are able to detect anomalous areas in terms of residuals. This kind of information can be useful for a better assessment of seismic risk and for promoting effective policies to reduce major damages.     

Seismic velocity and Poisson’s ratio tomography of the crust beneath southwest Anatolia: an insight into the occurrence of large earthquakes

The western part of Anatolia is one of the most seismically and tectonically active continental regions in the world, and much of it has been undergoing NS-directed extensional deformation since the Early Miocene. In this study, we determine 3-D tomographic images of the crust under the southwestern part of the North Anatolian Fault Zone by inverting a large number of arrival time data of P and S waves. From the obtained P- and S-wave velocity models, we estimated the Poisson’s ratio structures for a more reliable interpretation of the obtained anomalies. Our tomographic results confirmed the major tectonic features detected by previous studies and revealed new structural heterogeneities related to the active seismotectonics of the studied area. High P-wave velocity anomalies are recognized near the surface, while at deeper crustal layers, low P-wave velocities are widely distributed. The crustal S-wave velocity and Poisson’s ratio exhibit more structural heterogeneities compared to the P-wave velocity structure. Microearthquake activity is intense along highly heterogeneous zones in the southwestern part, which is characterized by low to high P-wave velocity, low S-wave velocity, and high Poisson’s ratio anomalies. Large earthquakes are also concentrated in zones dominated by low velocities and low to high Poisson’s ratios. Results of the checkerboard and synthetic tests indicate that the imaged anomalies are reliable features down to a depth of 25 km. Moreover, they are consistent with many geological and geophysical results obtained by other researchers along the southwestern part of the North Anatolian Fault Zone. An erratum to this article can be found at     

The ratio between corner frequencies of source spectra of <Emphasis Type="Italic">P</Emphasis>- and <Emphasis Type="Italic">S</Emphasis>-waves—a new discriminant between earthquakes and quarry blasts

This study analyzes and compares the P- and S-wave displacement spectra from local earthquakes and explosions of similar magnitudes. We propose a new approach to discrimination between low-magnitude shallow earthquakes and explosions by using ratios of P- to S-wave corner frequencies as a criterion. We have explored 2430 digital records of the Israeli Seismic Network (ISN) from 456 local events (226 earthquakes, 230 quarry blasts, and a few underwater explosions) of magnitudes Md?=?1.4–3.4, which occurred at distances up to 250 km during 2001–2013 years. P-wave and S-wave displacement spectra were computed for all events following Brune’s source model of earthquakes (1970, 1971) and applying the distance correction coefficients (Shapira and Hofstetter, Teconophysics 217:217–226, 1993; Ataeva G, Shapira A, Hofstetter A, J Seismol 19:389-401, 2015), The corner frequencies and moment magnitudes were determined using multiple stations for each event, and then the comparative analysis was performed.The analysis showed that both P-wave and especially S-wave displacement spectra of quarry blasts demonstrate the corner frequencies lower than those obtained from earthquakes of similar magnitudes. A clear separation between earthquake and explosion populations was obtained for ratios of P- to S-wave corner frequency f ₀(P)/f ₀(S). The ratios were computed for each event with corner frequencies f ₀ of P- and S-wave, which were obtained from the measured f ₀ ^I at individual stations, then corrected for distance and finally averaged. We obtained empirically the average estimation of f ₀(P)/f ₀(S)?=?1.23 for all used earthquakes, and 1.86 for all explosions. We found that the difference in the ratios can be an effective discrimination parameter which does not depend on estimated moment magnitude M _w .The new multi-station Corner Frequency Discriminant (CFD) for earthquakes and explosions in Israel was developed based on ratios P- to S-wave corner frequencies f ₀(P)/f ₀(S), with the empirical threshold value of the ratio for Israel as 1.48.