The Attenuation of Seismic Intensity in Italy, Part II: Modeling and Validation

Abstract Several different attenuation models have recently been proposed for the Italian region to characterize the decay of macroseismic intensity with the distance from the source. The significant scatter between these relationships and some significant drawbacks that seem to characterize previous approaches (described in a companion article by Pasolini et al. , 2008) suggest that the problem needs to be reconsidered. As a first step toward more detailed analyses in the future, this study aimed at developing an isotropic attenuation relationship for the Italian area. Because this attenuation relationship has to be used primarily in probabilistic seismic hazard assessment, major attention was given to evaluating the attenuation relationship in its complete probabilistic form. Another important aspect of this study was the preliminary evaluation of the intrinsic (i.e., independent of the specific attenuation relationship to be used) scattering of data, which represents the lowest threshold for the residual variance that cannot be explained by the attenuation relationship. Furthermore, the peculiar formal features of intensity data and relevant uncertainties were considered carefully. To reduce possible biases, the completeness of the available database was checked and a suitable data selection procedure was applied. Since epicentral intensity cannot be defined unambiguously from the experimental point of view, the attenuation relationship was scaled with a new variable that is more representative of the earthquake dimension. Several criteria were considered when evaluating competing attenuation formulas (explained variance, Bayesian information criteria, Akaike information criteria, etc.). Statistical uncertainty about empirical parameters was evaluated by using standard approaches and bootstrap simulations. The performance of the selected relationship with respect to a control sample was analyzed by using a distribution-free approach. The resulting equation for the expected intensity I at a site located at epicentral distance R is where , h =(3.91±0.27) km, and I E is the average expected intensity at the epicenter for a given earthquake that can be computed from the intensity data (when available) or by using empirical relationships with the moment magnitude M w or the epicentral intensity I 0 reported by the Italian seismic catalog Comparison of the model standard deviation (S.D.) (0.69 intensity degrees) with the intrinsic one (0.62) indicates that this attenuation equation is not far from being optimal.

[1]  G. Schwarz.  Estimating the Dimension of a Model , 1978 .

[2]  David M. Perkins,et al.  SEISRISK III; a computer program for seismic hazard estimation , 1987 .

[3]  N. Metropolis,et al.  The Monte Carlo method. , 1949 .

[4]  Carl Allin Cornell,et al.  Probabilistic Analysis of Damage to Structures under Seismic Loads , 1971 .

[5]  Gianluca Valensise,et al.  Defining seismogenic sources from historical earthquake felt reports , 1999, Bulletin of the Seismological Society of America.

[6]  Yosihiko Ogata,et al.  Statistical relations between the parameters of aftershocks in time, space, and magnitude , 1997 .

[7]  M. S. Barbano,et al.  The attenuation of seismic intensity in the Etna region and comparison with other Italian volcanic districts , 2006 .

[8]  D. Albarello,et al.  Attenuation relationship of macroseismic intensity in Italy for probabilistic seismic hazard assessment , 2004 .

[9]  U. Chandra,et al.  Attenuation of intensities in Iran , 1979, Bulletin of the Seismological Society of America.

[10]  Jung-Ki Kim,et al.  Intensity Attenuation in the Sino-Korean Craton , 2002 .

[11]  O. Scotti,et al.  Regional intensity attenuation models for France and the estimation of magnitude and location of historical earthquakes , 2006 .

[12]  Roger M.W. Musson,et al.  Intensity attenuation in the U.K. , 2005 .

[13]  Paolo Gasperini,et al.  The Attenuation of Seismic Intensity in Italy: A Bilinear Shape Indicates the Dominance of Deep Phases at Epicentral Distances Longer than 45 km , 2001 .

[14]  Indra N. Gupta,et al.  Spatial attenuation of intensities for central U.S. earthquakes , 1976, Bulletin of the Seismological Society of America.

[15]  D. E. Goldberg,et al.  Genetic Algorithms in Search , 1989 .

[16]  G. Grünthal.  European Macroseismic Scale 1998 , 1998 .

[17]  David M. Perkins,et al.  SEISRISK II; a computer program for seismic hazard estimation , 1982 .

[18]  J. Bommer,et al.  PREDICTION OF HORIZONTAL RESPONSE SPECTRA IN EUROPE , 1996 .

[19]  Gianluca Valensise,et al.  Catalogue of Strong Italian Earthquakes from 461 B.C. to 1997 (Appendix to volume 43 N° 4, 2000) , 2000 .

[20]  John E. Dennis,et al.  Numerical methods for unconstrained optimization and nonlinear equations , 1983, Prentice Hall series in computational mathematics.

[21]  David E. Goldberg,et al.  Genetic Algorithms in Search Optimization and Machine Learning , 1988 .

[22]  William L. Goffe,et al.  SIMANN: FORTRAN module to perform Global Optimization of Statistical Functions with Simulated Annealing , 1992 .

[23]  P. Gasperini,et al.  Catalogo dei forti terremoti in Italia dal 461 a.C. , 1995 .

[24]  Dario Albarello,et al.  Estimates of site seismicity rates using ill-defined macroseismic data , 1994 .

[25]  P. Gasperini,et al.  Lateral variations of seismic intensity attenuation in Italy , 2003 .

[26]  P. Hall.  The Bootstrap and Edgeworth Expansion , 1992 .

[27]  D. Carroll.  Chemical laser modeling with genetic algorithms , 1996 .

[28]  A. Blake.  On the estimation of focal depth from macroseismic data , 1937 .

[29]  U. Chandra,et al.  Attenuation of intensities in the United States , 1979, Bulletin of the Seismological Society of America.

[30]  Goldberg,et al.  Genetic algorithms , 1993, Robust Control Systems with Genetic Algorithms.

[31]  D. Albarello,et al.  Validation of Intensity Attenuation Relationships , 2005 .

[32]  D. Albarello,et al.  The Attenuation of Seismic Intensity in Italy, Part I: Theoretical and Empirical Backgrounds , 2008 .

[34]  U. Chandra,et al.  Attenuation of intensities and effect of local site conditions on observed intensities during the Corinth, Greece, earthquakes of 24 and 25 February and 4 March 1981 , 1985 .

[35]  D. Wells,et al.  New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement , 1994, Bulletin of the Seismological Society of America.

[36]  Antonio Pugliese,et al.  Attenuation of peak horizontal acceleration and velocity from italian strong-motion records , 1987 .

[37]  B. F. Howell,et al.  Attenuation of Modified Mercalli intensity with distance from the epicenter , 1975, Bulletin of the Seismological Society of America.

[38]  D. Mayer-Rosa,et al.  The practice of earthquake hazard assessment , 1993 .

[39]  Paolo Gasperini,et al.  Deriving numerical estimates from descriptive information: the computation of earthquake parameters , 2000 .

[40]  Clifford M. Hurvich,et al.  Regression and time series model selection in small samples , 1989 .

[41]  F. Sabetta,et al.  Estimation of response spectra and simulation of nonstationary earthquake ground motions , 1996, Bulletin of the Seismological Society of America.

[42]  Claudio Margottini,et al.  Intensity versus ground motion: A new approach using Italian data , 1992 .

[43]  W. Bakun,et al.  MMI Attenuation and Historical Earthquakes in the Basin and Range Province of Western North America , 2006 .

[44]  Silvia Castellaro,et al.  Regression problems for magnitudes , 2006 .

[45]  Robert Tibshirani,et al.  Bootstrap Methods for Standard Errors, Confidence Intervals, and Other Measures of Statistical Accuracy , 1986 .

[46]  H. Akaike.  A new look at the statistical model identification , 1974 .

[47]  John H. Holland,et al.  Adaptation in Natural and Artificial Systems: An Introductory Analysis with Applications to Biology, Control, and Artificial Intelligence , 1992 .