Murru, Maura

In this study, we compute the effect of stress change due to previous historical earthquakes on the probability of occurrence of future earthquakes on neighboring faults. Following a methodology developed in the last decade, we start from the estimate of the probability of occurrence in the next 50 years for a characteristic earthquake on known seismogenic structures, based on a time-dependent renewal model. Then a physical model for the Coulomb stress change caused by previous earthquakes on these structures is applied. The influence of this stress change on the occurrence rate of characteristic earthquakes is computed, taking into account both permanent (clock advance) and temporary (rate-and-state) perturbations. We apply this method to the computation of earthquake hazard of the main seismogenic structures recognized in the Central and Southern Apennines region, for which both historical and paleoseismological data are available. This study provides the opportunity of reviewing the problems connected with the estimate of the parameters of a renewal model in case of characteristic earthquakes characterized by return times longer than the time spanned by the available catalogues and the applicability of the concept of characteristic earthquake itself. The results show that the estimated effect of earthquake interaction in this region is small compared with the uncertainties affecting the statistical model used for the basic time-dependent hazard assessment.

The phenomenon of earthquake clustering, i.e., the increase of occurrence probability for seismic events close in space and time to other previous earthquakes, has been modeled both by statistical and physical processes. From a statistical viewpoint the so-called epidemic model (ETAS) introduced by Ogata in 1988 and its variations have become fairly well known in the seismological community. Tests on real seismicity and comparison with a plain time-independent Poissonian model through likelihood-based methods have reliably proved their validity. On the other hand, in the last decade many papers have been published on the so-called Coulomb stress change principle, based on the theory of elasticity, showing qualitatively that an increase of the Coulomb stress in a given area is usually associated with an increase of seismic activity. More specifically, the rate-and-state theory developed by Dieterich in the ′90s has been able to give a physical justification to the phenomenon known as Omori law. According to this law, a mainshock is followed by a series of aftershocks whose frequency decreases in time as an inverse power law. In this study we give an outline of the above-mentioned stochastic and physical models, and build up an approach by which these models can be merged in a single algorithm and statistically tested. The application to the seismicity of Japan from 1970 to 2003 shows that the new model incorporating the physical concept of the rate-and-state theory performs not worse than the purely stochastic model with two free parameters only. The numerical results obtained in these applications are related to physical characters of the model as the stress change produced by an earthquake close to its edges and to the A and σ parameters of the rateand- state constitutive law.

The recent Amatrice strong event (Mw6.0) occurred on August 24, 2016 in Central Apennines (Italy) in a seismic gap zone, motivated us to study and provide better understanding of the seismic hazard assessment in the macro area defined as “Central Italy”. The area affected by the sequence is placed between the Mw6.0 1997 Colfiorito sequence to the north (Umbria-Marche region) the Campotosto area hit by the 2009 L’Aquila sequence Mw6.3 (Abruzzo region) to the south. The Amatrice earthquake occurred while there was an ongoing effort to update the 2004 seismic hazard map (MPS04) for the Italian territory, requested in 2015 by the Italian Civil Protection Agency to the Center for Seismic Hazard (CPS) of the Istituto Nazionale di Geofisica e Vulcanologia INGV. Therefore, in this study we brought to our attention new earthquake source data and recently developed ground-motion prediction equations (GMPEs). Our aim was to validate whether the seismic hazard assessment in this area has changed with respect to 2004, year in which the MPS04 map was released. In order to understand the impact of the recent earthquakes on the seismic hazard assessment in central Italy we compared the annual seismic rates calculated using a smoothed seismicity approach over two different periods; the Parametric Catalog of the Historical Italian earthquakes (CPTI15) from 1871 to 2003 and the historical and instrumental catalogs from 1871 up to 31 August 2016. Results are presented also in terms of peak ground acceleration (PGA), using the recent ground-motion prediction equations (GMPEs) at Amatrice, interested by the 2016 sequence.

We develop an ensemble earthquake rate model that provides spatially variable time-independent (Poisson) long-term annual occurrence rates of seismic events throughout Italy, for magnitude bin of 0.1 units from Mw ≥4.5 in spatial cells of 0.1° x 0.1°. We weighed seismic activity rates of smoothed seismicity and fault-based inputs to build our earthquake rupture forecast model, merging it into a single ensemble model. Both inputs adopt a tapered Gutenberg-Richter relation with a single b-value and a single corner magnitude estimated by earthquakes catalog. The spatial smoothed seismicity was obtained using the classical kernel smoothing method with the inclusion of magnitude dependent completeness periods applied to the Historical (CPTI15) and Instrumental seismic catalogs. For each seismogenic source provided by the Database of the Individual Seismogenic Sources (DISS), we computed the annual rate of the events above Mw4.5, assuming that the seismic moments of the earthquakes generated by each fault are distributed according to the tapered Gutenberg-Richter relation with the same parameters of the smoothed seismicity models. Comparing seismic annual rates of the catalogs with those of the seismogenic sources, we realized that there is a good agreement between these rates in Central Apennines zones, whereas the seismogenic rates are higher than those of the catalogs in the north east and south of Italy. We also tested our model against the strong Italian earthquakes (Mw5.5+), in order to check if the total number (N-test) and the spatial distribution (S-test) of these events was compatible with our model, obtaining good results, i.e. high p-values in the test. The final model will be a branch of the new Italian seismic hazard map.

A systematic analysis is made of static Coulomb stress changes and earthquake occurrence in the area of the North Aegean Sea, Greece, in order to assess the prospect of using static stress changes to construct a regional earthquake likelihood model. The earthquake data set comprises all events of magnitude M ≥ 5.2 which have occurred since 1964. This is compared to the evolving stress field due to constant tectonic loading and perturbations due to coseismic slip associated with major earthquakes (M ≥ 6.4) over the same period. The stress was resolved for sixteen fault orientation classes, covering the observed focal mechanisms of all earthquakes in the region. Analysis using error diagrams shows that earthquake occurrence is better correlated with the constant tectonic loading component of the stress field than with the total stress field changes since 1964, and that little, if any, information on earthquake occurrence is lost if only the maximum of the tectonic loading over the fault orientation classes is considered. Moreover, the information on earthquake occurrence is actually increased by taking the maximum of the evolving stress field since 1964, and of its coseismic–slip component, over the fault orientation classes. The maximum, over fault orientation classes, of linear combinations of the tectonic loading and the evolving stress field is insignificantly better correlated with earthquake occurrence than the maximum of the tectonic loading by itself. A composite stress–change variable is constructed from ordering of the maximum tectonic loading component and the maximum coseismic–slip component, in order to optimize the correlation with earthquake occurrence. The results indicate that it would be difficult to construct a time–varying earthquake likelihood model from the evolving stress field that is more informative than a time–invariant model based on the constant tectonic loading.

One important aspect of the seismicity is the spatiotemporal clustering; hence, the distinction between independent and triggered events is a critical part of the analysis of seismic catalogs. Stochastic declustering of seismicity allows a probabilistic distinction between these two kinds of events. Such an approach, usually performed with the epidemic‐type aftershock sequence (ETAS) model, avoids the bias in the estimation of the frequency–magnitude distribution parameters if we consider a subset of the catalog, that is, only the independent or the triggered events. In this article, we present a framework to properly include the probabilities of any event to be independent (or triggered) both in the temporal variation of the seismic rates and in the estimation of the b‐value of the Gutenberg–Richter law. This framework is then applied to a high‐definition seismic catalog in the central part of Italy covering the period from April 2010 to December 2015. The results of our analysis show that the seismic activity from the beginning of the catalog to March 2013 is characterized by a low degree of clustering and a relatively high b‐value, whereas the following period exhibits a higher degree of clustering and a smaller b‐value.

In this paper, we gather and take stock of the results produced by the Operational Earthquake Forecasting (OEF) system in Italy, during its first 10 yr of operativity. The system is run in real-time: every midnight and after each ML 3.5 + event, it produces the weekly forecast of earthquakes expected by an ensemble model in each cell of a spatial grid covering the entire Italian territory. To e v aluate the performance skill of the OEF-Italy forecasts, we consider here standard tests of the Collaboratory for the Study of Earthquake Predictability, which have been opportunely adapted to the case of the overlapped weekly OEF forecasts; then we also adopt new performance measures borrowed from other research fields, like meteorology, specific to validate alarm-based systems by a binary criterion (forecast: yes/no; occurrence: yes/no). Our final aim is to: (i) investigate possible weaknesses and room for improvements in the OEF-Italy stochastic modelling, (ii) provide performance measures that could be helpful for stakeholders who act through a boolean logic (making an action or not) and (iii) highlight possible features in the Italian tectonic seismic activity.

In our paper we analyze the data base obtained from the observations of the Italian Seismological Network from 1975 to 1994 by using a simple algotithm to determine the rate of occurrence of seismic events condi- tioned by the occurrence of previous events after a period of quiescence. The number of observed pairs of earthquakes depends on several parameters: the magnitude threshold of the two events, the spatial and tempo- ral ranges of the quiescence period preceding the first (non aftershock) event, the time elapsed between the first and the second events and the spatial dimension of the alarm area. The Akaike information criterion was adopted to assess the optimal set of space-time parameters used in the definition of non-aftershock (events not related to a stronger previous one). In Central Italy, the rate of M ³3.8 earthquakes preceded by at least one M ~ 3.3 foreshock within 14.1 km and 2 days is 30%, while the rate of M ~ 3.3 earthquakes followed by a M ~ 3.8 mainshock in the same space time range is 7%. We observed that the probability that an earthquake of magnitude MI will be followed by an earthquake of magnitude M2 (success rate) fits the law log À = a+b (Mi -M2) with b approximately equal to l. By computing the success rate for given values of magnitude threshold of the first and the second events over a dense grid of spatial coordinates, we obtained maps of this feature over the investigated area. The results of this process document variations larger than a factor of five in the success rate over the Italian territory.

The aim of this work is to quantitatively set up a simple hypothesis for occurrence of earthquakes conditioned by prior events, on the basis of a previously existing model and the use of recent instrumental observations. A simple procedure is presented in order to determine the conditional probability of pairs of events (foreshock-mainshock, mainshock-aftershock) with short time and space separation. The first event of a pair should not be an aftershock, i.e., itmust not be related to a stronger previous event. The Italian earthquake catalog of the Istituto Nazionale di Geofisica (ING) (1975–1995, M > 3:4), the earthquake catalog of the Japan Meteorological Agency (JMA) (1983– 1994, M > 3:0) and that of the National Observatory of Athens (NOA) (1982–1994, M > 3:8) were analyzed. The number of observed pairs depends on several parameters: the size of the space-time quiescence volume defining nonaftershocks, the inter event time, the minimum magnitude of the two events, and the spatial dimension of the alarm volume after the first event. The Akaike information criterion has been adopted to assess the optimum set of space-time parameters used in the definition of the pairs, assuming that the occurrence rate of subsequent events may be modeled by two Poisson processes with different rates: the higher rate refers to the space-time volume defined by the alarms and the lower one simulates earthquakes that occur in the nonalarm space-time volume. On the basis of the tests carried out on the seismic catalog of Italy, the occurrence rate of M > 3:8 earthquakes followed by a M > 3:8 mainshock within 10 km and 10 days (validity) is 0.459. We have observed, for all three catalogs, that the occurrence rate density for the second event of a couple (mainshock or aftershock) of magnitude M2 subsequent to a non aftershock of magnitude M1 in the time range T can be modeled by the following relationship: .T , M2/ D 10a 0Cb.M1􀀀M2/ with b varying from 0.74 (Japan) to 1.09 (Greece). The decrease of the occurrence rate in time for a mainshock after a foreshock or for large aftershocks after a mainshock, for all three databases, obeys the Omori’s law with p changing from 0.94 (Italy) to 2.0 (Greece).

The stability of seismicity rate in central and southern Italy was ex- amined in two data sets: from 1975.0 to 1995.0 with M = 3.4 and from 1987.5 to 1996.0 with M _--> 2.5. These are the approximate minimum magnitudes of complete reporting for the respective periods. The first set was used to evaluate the possibility that the 23 November 1980 Irpinia (M 6.9) earthquake was preceded by precursory seismic quiescence; the second was used to evaluate the conditions under which a current seismic quiescence could be detected in central or southern Italy. During the years before the Irpinia earthquake, the seismicity rate in the northern half of the rupture area and north of it was low. Whether this was a case of precursory quiescence or not is subject to interpretation because the background rate cannot be established for the years before 1975. If we accept the relatively constant seismicity rate in the Irpinia volume during the decade after the mainshock as representative for the back- ground rate, we have a clear case of precursory quiescence lasting at least 1.3 yr up to the mainshock. Alternatively, it can be postulated that the seismicity rate during the decade following this shock was elevated regionally because of the stress redis- tribution and that the low rates seen before it represent the normal background rate. Even if this reasoning is accepted, the fact remains that a volume including the northern part of the rupture produced no M _-> 3.4 earthquakes during 1.3 yr before the Irpinia earthquake but produced 10 earthquakes during the 4.7 previous yr. Given these facts, we favor the interpretation that the Irpinia 1980 earthquake was preceded by precursory seismic quiescence. In the entire data set, there are three other cases of quiescence of higher significance without a mainshock following. Since no other mainshock with M > 6 exists in the data set, no missed event exists. We propose that in Italy precursory seismic quiescence may precede mainshocks and that it may be detected in the future by the improved catalog of the modern data set beginning in 1987.5. Major magnitude scale changes give the mistaken appearance that fewer large-magnitude earthquakes occur in Italy now than in years before 1987. We pos- tulate that the rate of earthquakes has not changed and that the magnitude scale should be adjusted.