Magnitude–intensity and intensity–attenuation relationships for atlas region and Algerian earthquakes

This paper presents the results of an investigation of the magnitude–intensity and intensity–attenuation relationships for earthquakes in the Atlas block and Algeria using macroseismic data. This work is based on a selected sample of isoseismal maps from 32 events which were recently revised. Surface-wave magnitudes, M_s, are recalculated using the Prague formula and range from 4·2 to 7·45. Because the Atlas mountains block is in a collision zone, earthquakes occur in general within a layer 15 km deep. Expressions of general form for the magnitude–intensity and intensity–attenuation correlations are adopted and are, respectively, and where R² = d² + h², d the source distance in km, h the focal depth in km, M_s the revised surface-wave magnitude, M_sc the predicted surface-wave magnitude, I_i the intensity at isoseismal i, I the predicted intensity, σ the standard deviation and P is zero for 50-percentile values and one for 84-percentile, and the coefficients A's and B's are determined by regression analysis. The results of this study show that the intensity–attenuation models are adequate to predict quite well the die-out of intensity with distance in the Atlas zone and coastal Algeria; it is also found that magnitude can be predicted accurately by calibrating isoseismal radii against revised instrumental surface-wave magnitude. Such magnitude–intensity relationships may be used to evaluate the magnitude of historical earthquakes in the region under survey, with no instrumental data, for which isoseismal radii and intensities are available.     

The French seismic CATalogue (FCAT-17)

In regions that undergo low deformation rates, as is the case for metropolitan France (i.e. the part of France in Europe), the use of historical seismicity, in addition to instrumental data, is necessary when dealing with seismic hazard assessment. This paper presents the strategy adopted to develop a parametric earthquake catalogue using moment magnitude M_w, as the reference magnitude scale to cover both instrumental and historical periods for metropolitan France. Work performed within the framework of the SiHex (SIsmicité de l’HEXagone) (Cara et al. Bull Soc Géol Fr 186:3–19, 2015. doi: 10.2113/qssqfbull.186.1.3) and SIGMA (SeIsmic Ground Motion Assessment; EDF-CEA-AREVA-ENEL) projects, respectively on instrumental and historical earthquakes, have been combined to produce the French seismic CATalogue, version 2017 (FCAT-17). The SiHex catalogue is composed of ~40,000 natural earthquakes, for which the hypocentral location and M_w magnitude are given. In the frame of the SIGMA research program, an integrated study has been realized on historical seismicity from intensity prediction equations (IPE) calibration in M_w detailed in Baumont et al. (submitted) companion paper to their application to earthquakes of the SISFRANCE macroseismic database (BRGM, EDF, IRSN), through a dedicated strategy developed by Traversa et al. (Bull Earthq Eng, 2017. doi: 10.1007/s10518-017-0178-7) companion paper, to compute their M_w magnitude and depth. Macroseismic data and epicentral location and intensity used both in IPE calibration and inversion process, are those of SISFRANCE without any revision. The inversion process allows the main macroseismic field specificities reported by SISFRANCE to be taken into account with an exploration tree approach. It also allows capturing the epistemic uncertainties associated with macroseismic data and to IPEs selection. For events that exhibit a poorly constrained macroseismic field (mainly old, cross border or off-shore earthquakes), joint inversion of M_w and depth is not possible, and depth needs to be fixed to calculate M_w. Regional a priori depths have been defined for this purpose based on analysis of earthquakes with a well constrained macroseismic field where joint inversion of M_w and depth is possible. As a result, 27% of SISFRANCE earthquake seismological parameters have been jointly inverted and for the other 73% M_w has been calculated assuming a priori depths. The FCAT-17 catalogue is composed of the SIGMA historical parametric catalogue (magnitude range between 3.5 up to 7.0), covering from AD463 to 1965, and of the SiHex instrumental one, extending from 1965 to 2009. Historical part of the catalogue results from an automatic inversion of SISFRANCE data. A quality index is estimated for each historical earthquake according to the way the events are processed. All magnitudes are given in M_w which makes this catalogue directly usable as an input for probabilistic or deterministic seismic hazard studies. Uncertainties on magnitudes and depths are provided for historical earthquakes following calculation scheme presented in Traversa et al. (2017). Uncertainties on magnitudes for instrumental events are from Cara et al. (J Seismol 21:551–565, 2017. doi: 10.1007/s10950-016-9617-1).     

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.     

The Problem of Regionalization in the Study of Magnitude-intensity Correspondences in the Mediterranean Area: A Distribution-free Approach

--A distribution-free statistical approach is proposed for tackling the problem of regionalization devoted to the study of magnitude-intensity relationships in the Mediterranean area. The training-set constitutes a compilation of more than 2000 earthquakes which occurred in the Mediterranean region since the end of the nineteenth century, when instrumental data became available, with epicentral or maximum intensity S VI and M_S magnitude values. Firstly an empirical magnitude-intensity correspondence has been computed for each intensity class by using the entire data set in the assumption of homogeneity at the regional scale. Residuals of such relation have been analyzed by a distribution-free statistical approach in order to evaluate the opportunity of a regionalization able to locally improve the performances of magnitude-intensity relationships. The analysis indicates that data concerning larger earthquakes (intensity S VII) do not suggest the opportunity of zonation, and that unbiased estimates of macroseismic magnitude can also be obtained in the assumption that magnitude-intensity correspondences are uniform over the entire Mediterranean area. Therefore, better constrained relations determined for the entire Mediterranean region should be preferred to ill-defined local ones. As concerns smaller events (intensity VI), the procedure suggests that medium/small-scale lateral variations (on a wavelength lower than 10² km) should be taken into account if an efficient estimate of magnitudes from maximum observed intensity is pursued, but that data presently available are not sufficient to suggest any reliable zonation of the area under study.     

A ground motion prediction equation for JMA instrumental seismic intensity for shallow crustal earthquakes in active tectonic regimes

The JMA (Japan Meteorological Agency) seismic intensity scale has been used in Japan as a measure of earthquake ground shaking effects since 1949. It has traditionally been assessed after an earthquake based on the judgment of JMA officials. In 1996 the scale was revised as an instrumental seismic intensity measure (I_JMA) that could be used to rapidly assess the expected damage after an earthquake without having to conduct a survey. Since its revision, Japanese researchers have developed several ground motion prediction equations (GMPEs) for I_JMA using Japanese ground motion data. In this paper, we develop a new empirical GMPE for I_JMA based on the strong motion database and functional forms used to develop similar GMPEs for peak response parameters as part of the PEER (Pacific Earthquake Engineering Research Center) Next Generation Attenuation (NGA) project. We consider this relationship to be valid for shallow crustal earthquakes in active tectonic regimes for moment magnitudes ( M ) ranging from 5.0 up to 7.5–8.5 (depending on fault mechanism) and rupture distances ranging from 0 to 200 km. A comparison of this GMPE with relationships developed by Japanese researchers for crustal and shallow subduction earthquakes shows relatively good agreement among all of the relationships at M 7.0 but relatively poor agreement at small magnitudes. Our GMPE predicts the highest intensities at small magnitudes, which together with research on other ground motion parameters, indicates that it provides conservative or upwardly biased estimates of I_JMA for M <5.5. Copyright © 2010 John Wiley & Sons, Ltd.     

A Tectonophysical Analysis of Earthquake Frequency—Size Relationship Types for Catastrophic Earthquakes in Central Asia

We performed a tectonophysical analysis of earthquake frequency–size relationship types for large Central Asian earthquakes in the regions of dynamical influence due to major earthquake-generating faults based on data for the last 100 years. We identified four types of frequency–size curves, depending on the presence/absence of characteristic earthquakes and the presence or absence of a downward bend in the tail of the curve. This classification by the shape of the tail in frequency–size relationships correlates well with the values of the maximum observed magnitude. Thus, faults of the first type (there are characteristic earthquakes, but no downward bend) with M_max ≥ 8.0 are classified as posing the highest seismic hazard; faults with characteristic earthquakes and a bend, and with M_max = 7.5–7.9, are treated as rather hazardous; faults of the third type with M_max = 7.1–7.5 are treated as posing potential hazard; and lastly, faults with a bend, without characteristic earthquakes, and with a typical magnitude M_max ≤ 7.0, are classified as involving little hazard. The tail types in frequency–size curves are interpreted using the model of a nonlinear multiplicative cascade. The model can be used to treat different tail types as corresponding to the occurrence/nonoccurrence of nonlinear positive and negative feedback in earthquake rupture zones, with this feedback being responsible for the occurrence of earthquakes with different magnitudes. This interpretation and clustering of earthquake-generating faults by the behavior the tail of the relevant frequency–size plot shows raises the question about the physical mechanisms that underlie this behavior. We think that the occurrence of great earthquakes is related to a decrease in effective strength (viscosity) in the interblock space of faults at a scale appropriate to the rupture zone size.     

Calibration of magnitude scales for earthquakes of the Mediterranean

In order to provide the tools for uniform size determination for Mediterranean earthquakes over the last 50-year period of instrumental seismology, we have regressed the magnitude determinations for 220 earthquakes of the European-Mediterranean region over the 1977–1991 period, reported by three international centres, 11 national and regional networks and 101 individual stations and observatories, using seismic moments from the Harvard CMTs. We calibrate M(M₀) regression curves for the magnitude scales commonly used for Mediterranean earthquakes (M_L, M_WA, m_b, M_S, M_LH, M_LV, M_D, M); we also calibrate static corrections or specific regressions for individual observatories and we verify the reliability of the reports of different organizations and observatories. Our analysis shows that the teleseismic magnitudes (m_b, M_S) computed by international centers (ISC, NEIC) provide good measures of earthquake size, with low standard deviations (0.17–0.23), allowing one to regress stable regional calibrations with respect to the seismic moment and to correct systematic biases such as the hypocentral depth for M_S and the radiation pattern for m_b; while m_b is commonly reputed to be an inadequate measure of earthquake size, we find that the ISC m_b is still today the most precise measure to use to regress M_W and M₀ for earthquakes of the European-Mediterranean region; few individual observatories report teleseismic magnitudes requiring specific dynamic calibrations (BJI, MOS). Regional surface-wave magnitudes (M_LV, M_LH) reported in Eastern Europe generally provide reliable measures of earthquake size, with standard deviations often in the 0.25–0.35 range; the introduction of a small (±0.1–0.2) static station correction is sometimes required. While the Richter magnitude M_L is the measure of earthquake size most commonly reported in the press whenever an earthquake strikes, we find that M_L has not been computed in the European-Mediterranean in the last 15 years; the reported local magnitudes M_WA and M_L do not conform to the Richter formula and are of poor quality and little use, with few exceptions requiring ad hoc calibrations similar to the M_S regression (EMSC, ATH). The duration magnitude M_D used by most seismic networks confirms that its use requires accurate station calibrations and should be restricted only to events with low seismic moments.     

Empirical Global Relations Converting M S and m b to Moment Magnitude

The existence of several magnitude scales used by seismological centers all over the world and the compilation of earthquake catalogs by many authors have rendered globally valid relations connecting magnitude scales a necessity. This would allow the creation of a homogeneous global earthquake catalog, a useful tool for earthquake research. Of special interest is the definition of global relations converting different magnitude scales to the most reliable and useful scale of magnitude, the moment magnitude, M _W. In order to accomplish this, a very large sample of data from international seismological sources (ISC, NEIC, HRVD, etc.) has been collected and processed. The magnitude scales tested against M _W are the surface wave magnitude, M _S, the body wave magnitude, m _b, and the local magnitude, M _L. The moment magnitudes adopted have been taken from the CMT solutions of HRVD and USGS. The data set used in this study contains 20,407 earthquakes, which occurred all over the world during the time period 1.1.1976–31.5.2003, for which moment magnitudes are available. It is shown that well-defined relations hold between M _W and m _b and M _S and that these relations can be reliably used for compiling homogeneous, with respect to magnitude, earthquake catalogs.     

The macroseismic field of the Balkan area

A catalogue of 356 macroseimic maps which are available for the Balkan area was compiled, including information on the source parameters of the corresponding earthquakes, the macroseismic parameters of their strength and their macroseismic field. The data analysis of this catalogue yields new empirical relations for attenuation, which can be applied for the calibration of historical events, modelling of isoseismals and seismic hazard assessment. An appropriate analysis allowed the separation and estimation of the average values of the geometrical spreading, n, and anelastic attenuation factor, c, for the examined area which were found equal to –3.227 ± 0.112 and –0.0033 ± 0.0010. Scaling relations for the focal macroseismic intensity, I_f, and the epicentral intensity I₀, versus the earthquake moment magnitude were also determined for each Balkan country. A gradual decrease of the order of 0.5 to 1 intensity unit is demonstrated for recent (after 1970) earthquakes in Greece. Finally the depths of the examined earthquakes as they robustly determined (error <5 km) on the basis of macroseismic data were found to have small values ( 10 km). However large magnitude earthquakes show higher focal depths ( 25 km), in accordance with an increase of the seismic fault dimensions for such events.     

Remarks about magnitudes and maximum intensities of the earthquakes in Friuli, 1976

Summary The magnitudesM and maximum intensitiesI ₀ of 66 earthquakes which occurred in the period from 1901 to 1973 in the eastern alpine area are compared with those of 62 earthquakes occurred in Friuli, northern Italy between 6 May and 15 September 1976. The average focal depth of the Friuli events is about 10 km. The existing empirical relation betweenM andI ₀ for this area is improved by many new data.     

Estimation of the probability of strongest seismic disasters based on the extreme value theory

This paper presents the review of the experience in applying the approach based on the limiting distributions of the extreme value theory (the generalized Pareto distribution, GPS, and generalized extreme value distribution, GEV) for deriving the distributions of maximal magnitudes and related ground accelerations from the earthquakes on the future time intervals of a given duration. The results of analyzing the global and regional earthquake catalogs and the ground peak accelerations during the earthquakes are described. It is shown that the magnitude of the strongest possible earthquake M _max (and analogous characteristics for other types of data), which is often used in seismic risk assessment, is potentially unstable. We suggest a stable alternative for M _max in the form of quantiles Q _q (τ) of the maximal possible earthquake, which could occur during the future time interval of length τ. The quantity of the characteristic maximal event M _c, which has been introduced in our previous publications, is another helpful robust scalar parameter. All the cases of approximation of the tails of empirical distributions, which were studied in our works, turned out to be finite (bounded); however, the rightmost point of these distributions, M _max, is often poorly detectable and unstable. Therefore, the M _max parameter has a low practical value.     

The Kaliningrad, West Russia earthquakes on the 21st of September 2004—Surprise events in a very low-seismicity area

The earthquakes in Kaliningrad, West Russia on the 21st of September 2004 were unexpected in a seismically quiet area. The main shock of magnitude m_b = 5 was widely felt around the Baltic Sea. A comparison with some historic earthquakes in Northern Europe shows that its perceptibility area is smaller than that of the 1904 Oslo Graben earthquake of an estimated magnitude (M_L) 5.4 but larger than those of the 1759 Kattegat and 1819 Lurøy earthquakes. The latter are claimed to have had magnitudes (M_S) in the range of 5.7-6.0. An analysis of the Lg phase of the Kaliningrad earthquakes as recorded at a number of European stations accords only weakly with the macroseismic intensity pattern that shows fast attenuation towards west-northwest and southeast. The strike-slip focal mechanism of the main shock is discussed in the context of remnant glacial rebound stresses in generating present-day seismicity in N. Europe.     

Inter-event times between the moderate-large earthquakes in Iran

Persian territory, which is dividable into major seismotectonic provinces, always suffers from damages of moderate and large earthquakes from ancient era to modern time. Therefore, temporal prediction of earthquake occurrence in this kind of area is an important topic. For this purpose, 628 moderate-large (5.5 ≤M_S≤ 8.2) earthquakes occurred in Persia during the period from 400 B.C. to 2015 C.E. were used. Considering the magnitudes of events preceding main shocks and the annual seismic moment release in seismic source areas in the provinces, we calibrated equations predicting inter-event time of occurrence of moderate and large earthquakes (M_W>5.5) in Iran. In each source area, inter-event times between moderate and large shocks with magnitudes equal to or larger than a certain cut-off magnitude (M_W5.5) were calculated. The inter-event times between the earthquakes were used to compute the relationships using multiple regression technique. Calculated relationships express the basic idea of the time predictable model predicting the occurrence time of the future main shock in a certain seismogen area. However, despite of unavoidable scatter in observations and uncertainties in the results, occurrence times of main shocks during the next years and decades in some source areas in Iran were determined.     

Relationship of Tsunami Intensity to Source Earthquake Magnitude as Retrieved from Historical Data

Operational prediction of near-field tsunamis in all existing Tsunami Warning Systems (TWSs) is based on fast determination of the position and size of submarine earthquakes. Exceedance of earthquake magnitude above some established threshold value, which can vary over different tsunamigenic zones, results in issuing a warning signal. Usually, a warning message has several (from 2 to 5) grades reflecting the degree of tsunami danger and sometimes contains expected wave heights at the coast. Current operational methodology is based on two main assumptions: (1) submarine earthquakes above some threshold magnitude can generate dangerous tsunamis and (2) the height of a resultant tsunami is, in general, proportional to the earthquake magnitude. While both assumptions are physically reasonable and generally correct, statistics of issued warnings are far from being satisfactory. For the last 55 years, up to 75% of warnings for regional tsunamis have turned out to be false, while each TWS has had at least a few cases of missing dangerous tsunamis. This paper presents the results of investigating the actual dependence of tsunami intensity on earthquake magnitude as it can be retrieved from historical observations and discusses the degree of correspondence of the above assumptions to real observations. Tsunami intensity, based on the Soloviev-Imamura scale is used as a measure of tsunami “size”. Its correlation with the M _s and M _w magnitudes is investigated based on historical data available for the instrumental period of observations (from 1900 to present).     

The University of Athens Hellenic Macroseismic Database (HMDB.UoA): historical earthquakes

A systematic study of historical earthquakes leading to the quantification of earthquake effects in terms of macroseismic data points (MDPs) and, consequently, earthquake parameters has been carried out in the last decade at the Laboratory of Seismology of the University of Athens. For each earthquake, the available background information was evaluated and the corresponding macroseismic intensities assessed in terms of the European Macroseismic Scale 1998. A considerable amount of these MDPs contributed to the Archive of Historical Earthquake Data inventory through European initiatives (NERIES and SHARE). Based on the structure of the European Database, the local version of the Hellenic Macroseismic Database (HMDB.UoA) was designed incorporating historical earthquakes of the period 1000–1899 from the eastern Aegean area, central Greece and Ionian Islands. In its present form, the HMDB.UoA includes 90 events with I _max?≥?7 (868 MDPs) and 1,088 events with I _max?<?7 (1,273 MDPs). The database is hosted on the website http://macroseismology.geol.uoa.gr/.     

Comparison between different earthquake magnitudes determined by China Seismograph Network

By linear regression and orthogonal regression methods, comparisons are made between different magnitudes (lo-cal magnitude ML, surface wave magnitudes MS and MS7, long-period body wave magnitude mB and short-period body wave magnitude mb) determined by Institute of Geophysics, China Earthquake Administration, on the basis of observation data collected by China Seismograph Network between 1983 and 2004. Empirical relations between different magnitudes have been obtained. The result shows that: 1 As different magnitude scales reflect radiated energy by seismic waves within different periods, earthquake magnitudes can be described more objectively by using different scales for earthquakes of different magnitudes. When the epicentral distance is less than 1 000 km, local magnitude ML can be a preferable scale; In case M<4.5, there is little difference between the magnitude scales; In case 4.5MS, i.e., MS underestimates magnitudes of such events, therefore, mB can be a better choice; In case M>6.0, MS>mB>mb, both mB and mb underestimate the magnitudes, so MS is a preferable scale for deter-mining magnitudes of such events (6.08.5, a saturation phenomenon appears in MS, which cannot give an accurate reflection of the magnitudes of such large events; 2 In China, when the epicentral distance is less than 1 000 km, there is almost no difference between ML and MS, and thus there is no need to convert be-tween the two magnitudes in practice; 3 Although MS and MS7 are both surface wave magnitudes, MS is in general greater than MS7 by 0.2~0.3 magnitude, because different instruments and calculation formulae are used; 4 mB is almost equal to mb for earthquakes around mB4.0, but mB is larger than mb for those of mB≥4.5, because the periods of seismic waves used for measuring mB and mb are different though the calculation formulae are the same.     

Application of the time-predictable model in Peninsular India for future seismic hazard assessment

It has been the belief among Earth scientists that the Peninsular Shield is aseismic, as the region attained stability long ago. However, the earthquake at Koyna (10 December 1967), Bhadrachalam (13 April 1969), Broach (23 March 1970), Hyderabad (30 June 1983), Khillari (30 September 1993), Jabalpur (22 May 1997), Gujarat (26 January 2001), and additional ones of smaller magnitudes, altered this concept. This area has experienced many widely distributed shallow earthquakes, some of them having large magnitudes. It is now widely accepted that seismic activity still continues with moderate events. Therefore, a need has arisen to take into consideration recent seismological data to assess the future seismic status of Peninsular India. Earthquake generation model has been studied to develop the statistical relations with surface wave magnitude (M _S ≥ 4.5). Five seismogenic sources showing clustering of earthquakes and including at least three main shocks of magnitude 4.5 ≤ M _S ≤ 6.5 giving two repeat times, have been identified. It is mainly based on the so-called “regional time-predictable model”. For the considered region it is observed that the time interval between two consecutive main shocks depends on the preceding main shock magnitude (M _p ) not on the following main shocks magnitude M _f suggesting the validity of time predictable model in the region.     

Frequency analysis of annual maximum earthquakes within a geographical region

The risk formula, expressing the probability of at least one occurrence of earthquakes of greater-than-design-value magnitudes over the economic life of a structure, is modified taking into consideration the probability of no-earthquake years. The annual maximum earthquake magnitudes of three scales: Richter magnitude, also known as local magnitude (M_L), body-wave magnitude (M_b), and moment magnitude (M_M) in a geographical area encompassing the Bingöl Province in Turkey are taken from two sources: (1) report by Kalafat et al. (2007) [14] and (2) the web site reporting data by Kandilli Observatory which has been recording earthquakes occurring in and around Turkey since 1900. Statistical frequency analyses are applied on the three sample series using various probability distribution models, and magnitude versus average return period relationships are determined. The values of the M_L, M_b, and M_M series for 10% and 2% risk are computed to be around 7.2 and 8.3. The tectonic structure and seismic properties of the Bingöl region are also given briefly.     

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.     

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.