Anisotropic radiation modelling of macroseismic intensities for estimation of the attenuation structure of the upper crust in Greece

A method is suggested for the analysis of macroseismic intensity data in order to accurately determine an average attenuation structure of the upper part of the crust in an area. The method is based on a model which assumes that the observed intensities depend on source properties (radiation pattern, size, focal depth), geometrical spreading and anelastic attenuation. The method is applied to 13,008 intensity values, observed in corresponding sites of Greece and grouped (in 4228 groups), according to their spatial clustering in order to diminish observational errors and site effects. An average intensity attenuation coefficient,c=–0.0039±0.0016, corresponding to a quality factor, Q=350±140, is determined for the upper 20 km of the crust in this area. This value is relatively low, in good agreement with the relatively high heat flow and high seismic activity of this area. A byproduct of the present study is the determination, for each earthquake, of a macroseismic focal depth and of a macroseismic size, which is strongly correlatted with both the earthquake's magnitude and its seismic moment determined by independent methods.