Analysis of macroseismic fields using statistical data depth functions: considerations leading to attenuation probabilistic modelling

Modelling seismic attenuation is one of the most critical points in the hazard assessment process. In this article we consider the spatial distribution of the effects caused by an earthquake as expressed by the values of the macroseismic intensity recorded at various locations surrounding the epicentre. Considering the ordinal nature of the intensity, a way to show its decay with distance is to draw curves—isoseismal lines—on maps, which bound points of intensity not smaller than a fixed value. These lines usually take the form of closed and nested curves around the epicentre, with highly different shapes because of the effects of ground conditions and of complexities in rupture propagation. Forecasting seismic attenuation of future earthquakes requires stochastic modelling of the decay on the basis of a common spatial pattern. The aim of this study is to consider a statistical methodology that identifies a general shape, if it exists, for isoseismal lines of a set of macroseismic fields. Data depth is a general nonparametric method for analysis of probability distributions and datasets. It has arisen as a statistical method to order points of a multivariate space, e.g., Euclidean space \({\mathbb {R}}^{p}\), \(p \ge 1\), according to the centrality with respect to a distribution or a given data cloud. Recently, this method has been extended to the ordering of functions and trajectories. In our case, for a fixed intensity decay \(\varDelta I\), we build a set of convex hulls that enclose the sites of felt intensity \(I_s \ge I_0 -\varDelta I\), one for each macroseismic field of a set of earthquakes that are considered as similar from the attenuation point of view. By applying data depth functions to this functional dataset, it is possible to identify the most central curve, i.e., the attenuation pattern, and to consider other properties like variability, outlyingness, and possible clustering of such curves. Results are shown for earthquakes that occurred on the Central Po Plain in May 2012, and on the eastern flank of Mt. Etna since 1865.