Bayesian Analysis of the Magnitude of Earthquakes Located in a Seismic Region of Italy

: Bayesian statistical procedures are well known for estimating the probability distribution of the occurrence of an event. In this study, we applied the Bayesian statistical method to estimate the distribution of magnitudes of earthquakes that occurred in central Italy in two different time periods. Using the Monte Carlo sampling method, we recovered the real magnitude distribution using just a small amount of available data, randomly selected by the seismic catalogue.


Introduction
The seismic hazard assessment by using the Bayesian probability theory was introduced by Cornell [1] and was especially applied in earthquake engineering problems [2]. The Bayesian discrete approach was performed to forecast the interevent times for large earthquakes along the western and eastern Hellenic arc [3], and to predict probabilities of occurrence of moderate and strong earthquakes in the main seismogenic zones in Greece [4]. Stavrakakis and Drakopoulos [5] used the Bayesian extreme-value distribution to assess the seismic hazard in some seismogenic sources in Greece. Egozcue and Ruttener [6] applied Bayesian approaches for assessing seismic hazard with imprecise data; introducing the uncertainties inherent in the inaccuracy and heterogeneity of the measuring systems by which the data were recorded, they gave special attention to "priors" in the Bayesian estimation. Wang et al. [7] developed a Bayesian algorithm to predict the next major event induced by the Meishan fault in central Taiwan, based on one magnitude observation of 6.4 magnitude from the last event, along with the prior data including fault length of 14 km, rupture width of 15 km, rupture area of 216 km 2 , average displacement of 0.7 m, slip rate of 6 mm/yr, and five earthquake empirical models.
In this work, we apply the Bayesian probability theory to estimate the magnitude distribution in one of the most seismically active areas of Italy (Figure 1), struck recently by several strong earthquakes with local magnitude larger than 5.5 [8]. It is well known that the magnitude distribution of an earthquake sequence is given by the socalled Guttenberg-Richter law [9], which relates a threshold magnitude with the number of earthquakes with magnitude > in a form of semi-logarithmic power law as log ( ) = − , where is the number of earthquakes of magnitude > , a is the productivity of the earthquake, and b is a value that indicates the proportion of small events with respect to the large ones. The value b is a critical parameter that can provide information about the conditions of the stress state of the investigation area and is generally estimated by the method of maximum likelihood [10][11][12].
In this work, our aim is to estimate by using the Bayesian approach the magnitude distribution of earthquakes occurred in central Italy on the base of a small percentage of real data.

Bayesian Analysis
The inferences in the Bayesian analysis are based on the calculation of the conditional posterior probabilities calculated by the Bayes theorem [13,14]. The root of the calculation of Bayesian inference is the Bayes' theorem

Modeling the Magnitude of Earthquakes in Central Italy with Bayesian Analysis
We analyzed the seismic catalogue (magnitude larger or equal to 1.5) from 1995 and 2018 of the area shown in Figure 1. We selected two groups of events, between 1995 and 2008 and from 2009 to 2018. Figure 2 shows the two sequences of events.   Using the Bayesian inference, we obtained the probability distribution functions of the magnitudes of the two sets of data, using 100 Monte Carlo samples for the Bayesian Method.
The Bayesian method represented in Figure 4 reproduces the probability density functions of the two groups rather well. Figures 5 and 6 show the errors between the probability density function of the data and that estimated by the Bayesian inference method.   Error between the probability density function of the data (shown in Figure 3) and that estimated by the Bayesian method for the period from 1995 to 2008. Figure 6. Error between the probability density function of the data (shown in Figure 3) and that estimated by the Bayesian method for the period from 2009 to 2018.

Prediction of the Magnitude with Bayesian Method
In the following example, the data from 2009 to 2017 are considered to calculate the distribution of the magnitudes of the events and then to make a prediction of this distribution in 2018, using 100 samples for the Bayesian method. Figure 7 shows the distribution of magnitude for the period from 2009 to 2017 and the distribution obtained through the Bayesian prediction in 2018.

Conclusions
The Bayesian inference method used in this paper can correctly reproduce the probability distributions of magnitudes of earthquakes occurred in central Italy from 1995 to 2018. In addition, with just a few observations, the probability distributions of magnitudes of the events can be reproduced with a rather small error between the distribution of the data and that one obtained by Bayesian inference.