Macroseismic intensity hazard maps for Italy based on a recent grid source model

Seismic hazard maps from probabilistic seismic hazard analysis or PSHA collect, at different sites, the values of the (site-specific) ground motion intensity measures of interest that, taken individually, have the same exceedance return period. For large-scale analyses, a widely used intensity measure is the macroseismic (MS) intensity, that provides an assessment of the earthquake effect based on the observed consequences in the hit area. Hazard maps can be developed in terms of MS intensity, and some examples exist in this respect. In the case of Italy, the last MS hazard map is based on the same seismic source model (known as MPS04) adopted to derive the design seismic actions of the current building code, a study dating more than ten years ago. It provides results in terms of countrywide Mercalli–Cancani–Sieberg (MCS) intensity level with 475 years return period. This short paper presents and discusses MCS probabilistic seismic hazard maps for Italy based on a recent grid-seismicity source model, herein named MPS19, synthetizing the large effort of a wide scientific community. The results, which are obtained by means of classical PSHA, are given in the form of maps referring to the 475 years return period, and also others of earthquake engineering interest. Moreover, it is discussed that the return period does not univocally identifies the MS intensity because, although MS is, by definition, a discrete random variable, it is modelled, in a given earthquake, by means of a normal distribution, that is, treated as continuous. Thus, the maps of the minimum return period causing the occurrence or exceedance of different MCS intensities are also provided. Finally, the comparison between the 475 years return period hazard map presented and the one which is currently the point of reference in Italy, that is, computed using MPS04, is briefly discussed. All the computed maps are made available to the reader as supplemental material.


Introduction
Macroseismic (MS) intensity provides a proxy for the impact of an earthquake based on its implications (macroscopically) observed on the affected communities, built environment and structures (e.g., Grünthal 1998). Different MS intensity scales exist (see Musson et al. 2010, for an overview); in general, they consist of grades, identified by ordinal numbers, that differ each other for the earthquake effects that they describe. Although site-specific ground motion intensity measures (IMs) are recognized far more appropriate to characterize earthquakes in the case of analysis and design of a specific structure (e.g., CEN 2004; C.S.LL.PP. 2018) at the construction site, MS intensity is still widely used. For example, right after an earthquake, ShakeMap (Wald et al. 1999a) elaborates the map of MS intensity for the hit area (e.g., Michelini et al. 2019), which may be used for rapid loss assessment. Also, MS intensity may serve to express the vulnerability of existing structures in an area of interest, based on empirical observations (e.g., Giovinazzi and Lagomarsino 2004). Furthermore, in those areas where the seismic monitoring network is sparse or even absent, the observed MS intensity can be used to infer ground motion IMs (e.g., peak ground acceleration; Faenza and Michelini 2010;Gomez-Capera et al. 2020) or to estimate the occurred event magnitude (e.g., Sibol et al. 1987;Azzaro et al. 2011). The Italian seismic catalog, so-called Catalogo Parametrico dei Terremoti Italiani or CPTI (Rovida et al. 2020), contains information, such as magnitude and location, about earthquakes occurred from 1000 AD to the end of 2017, most of which are retrieved based on macroseismic intensity estimates (i.e., effects assessed according to the historical evidences as provided by dedicated studies; e.g., Boschi et al. 2000) being no instrumental data available at the time of occurrence. The first earthquake of the catalog for which early instrumental data are available dates back to 1918 (Sandron et al. 2014) but, until the second half of sixties of the last century, the majority of the earthquakes' magnitude estimations are based on historical information.
For all the cited reasons, it may be worth to compute countrywide probabilistic maps, for fixed exceedance return periods T r , in terms of MS intensity, especially for countries, such as Italy, where most of the built environment is very slowly renovated, to say the least. Note that such maps probabilistically describe earthquake consequences, then they could be also seen as risk maps; however, they are usually identified in literature (and hereafter) as hazard maps. On the other hand, if the required risk metric is the damage state of a specific structural typology, the macroseismic intensity hazard needs to be combined with specific vulnerability models, that is, damage probability matrices (see, for example, Iervolino et al. 2015, among others).
In fact, macroseismic hazard assessment studies for Italy have been proposed in the past by Slejko et al. (1998) and Albarello et al. (2000). They used the source model of Meletti et al. (2000) and derived seismicity rates from the earthquake catalog of Camassi and Stucchi (1997). A few years later, another study, dealing with the hazard assessment in terms of MS intensity at the national scale, was published by Gomez Capera et al. (2010). It used the models that, still at the time of writing this paper, are at the basis of the hazard assessment, commonly named as MPS04, which is adopted by the current Italian building code and is considered as the reference MS intensity hazard assessment study for Italy. All the cited works refer to the Mercalli-Cancani- Sieberg (MCS;Sieberg 1930) scale, and present hazard maps collecting, at the national scale, the MCS intensities with ten percent probability in fifty years of occurring or being exceeded, that is 475 years (yr) return period.
All the hazard maps provided in the cited studies were computed by implementing the classical probabilistic seismic hazard analysis or PSHA framework (e.g., Reiter 1990;Mc Guire 2004). Site-specific PSHA first requires to identify the seismic sources in the region of interest and to characterize them in terms of mean annual number (i.e., the rate) of earthquakes (mainshocks, in fact) above a minimum magnitude (M) of interest, and magnitude distribution given the occurrence of one earthquake event (e.g., Gutenberg and Richter 1944). Finally, (at least) one ground motion prediction equation (sometimes referred to as attenuation relationship) is necessary to model, for any possible earthquake magnitude and location on the source, the probabilistic distribution of the ground motion intensity at the site of interest. The characterization of the sources is based on the available information, historical or instrumental, about past earthquakes, which are typically collected in the earthquake catalogs.
In fact, it is to mention that there are methodologies, alternative to classical PSHA, that may be employed to assess seismic hazard in terms of MS intensity. For example, with reference to Italy, Mucciarelli et al. (2000) and Albarello et al. (2002) provide a countrywide macroseismic hazard assessment profiting of the so-called site approach (Magri et al. 1994). The input of the analysis is defined for each site individually, and consists of a catalog collecting local macroseismic data; e.g., observations according to documentary sources or estimations by means of attenuation relationships. Gomez Capera et al. (2010) compare the results obtained using the site approach to those obtained by means of PSHA.
When PSHA is performed in terms of MS intensity, the implicit assumption is that the resolution of the assessment is such that the effects at a given location, identified by a point (i.e., a longitude-latitude pair) in the analysis, represent the large-scale effect on the built environment at that location. Moreover, the earthquakes magnitude distribution characterizing the seismic sources can be replaced by the distribution of the epicentral intensity I 0 , that is directly related to the historical information provided by the seismic catalogs. Indeed, for a given earthquake, I 0 is generally assumed as the largest MS intensity level observed in the epicentral area of the seismic event (e.g., Gruppo di Lavoro 2004). Alternatively, if a source model is already available in term of magnitude distribution, it is also possible to convert each magnitude value into the expected (i.e., not the observed) value of the MS intensity at the epicenter, I e (e.g., Albarello et al. 2007). Regardless of whether the input is in terms of I 0 or I e , once the propagation model for MS intensity is chosen (e.g., Ambraseys 1985), the PSHA allows to compute, for a given point, a MS intensity value with a given return period and, extending the analysis to a grid of points, the sought map of MS intensity values characterized by the same return period can be computed; this is despite the fact that MS intensity is defined by means of discrete grades, and some issues may be raised in this respect.
Recently, a set of new seismic hazard models for Italy, named MPS19, was developed as a possible update of the previous source model for the country Visini et al. 2021). The objective of the simple study presented herein is to discuss the seismic hazard assessment for Italy in terms of MS intensity, considering the MCS scale and adopting a source model, not coincident with, yet based on, MPS19 (and retaining its name herein). Four hazard maps, referred to as IMPS19, derived via classical PSHA, are discussed, that is, those showing the MS intensity with return periods equal to 50, 475, 975 and 2475 yr, corresponding to 63%, 10%, 5% and 2% exceedance probability in 50 yr (a typical exposure time in earthquake engineering), respectively.
The remainder of the paper is structured such that the basics of PSHA are recalled, first. After describing the source and the propagation models used in the analysis, the hazard maps for the four selected return periods are presented. Subsequently, for each MCS grade, the lowest return period causing its occurrence or exceedance is mapped to discuss some issues deriving from the fact that attenuation relationship models the MS intensity as a continuous random variable. Finally, the IMPS19 hazard map for T r = 475 yr is briefly compared to the counterpart provided of Gomez Capera et al. (2010), referred to as IMPS04. Some final remarks close the study. Hazard maps developed in this study are provided as supplemental material.

Methodology
This section describes the way in which PSHA is implemented in the case the location of interest is affected by s point-like seismic sources and the seismicity of each source is characterized in terms of I e . The analysis ultimately provides the (usually annual) rate of earthquakes that cause MS intensity larger than a certain threshold (ms) at the location, indicated as MS>ms , as per Eq. (1): The exceedance rate is the reciprocal of the exceedance return period. The equation is written assuming that, for each source, seismicity is defined in terms of rate of earthquakes associated to bins, say n in number, of I e . Index k = {1, 2, … , s} identifies each of the sources affecting the seismic hazard at the location of interest and k,j is the rate of earthquakes, with a specific epicentral intensity I e = i e,j , j = {1, 2, … , n} , occurring at the k-th source. The term P MS > ms | | | I e = i e,j , R = r k , provided by the propagation model, denotes the conditional probability that MS is larger than ms given that the considered location is subjected to an earthquake with I e = i e,j and that the epicentral distance from the k-th source is R = r k . It should be noted that propagation models usually consider MS as a continuous random variable (RV), although macroseismic scales are usually discrete. This is an issue that typically arises with the analysis of macroseismic data for developing propagation models or also relationships between MS intensity and other ground motion parameters (e.g., Wald et al. 1999b;Faenza and Michelini 2010). In PSHA calculations, threating MS as continuous RV allows to obtain a unique MS value for any return period of interest. A side effect is that, when hazard maps are presented discretizing the MS intensity in accordance with the MCS scale, at a given grid point, the same MS intensity could be associated to different return periods. Due to this issue, the maps of the minimum return period associated to each MS intensity grade are also discussed in the following.

Input models
Recently, a new seismic hazard assessment for Italy, considering a large set of pseudospectral accelerations as the ground motion IM, has been released by Meletti et al. (2021). PSHA is developed via ninety-four seismic source models weighted according to statistical performances in describing the earthquake occurrences and taking into account an experts' elicitation process . The ensemble of the models are combined, via a logic tree approach, with a set of selected ground motion prediction equations . Thus, such a seismic model results in about six-hundred branches of the logic tree. (With respect to the specific purposes of the analyses herein discussed, it must be noted that the several seismic source models in the logic tree are characterized by methodologies, such as geodetic data, which do not allow a direct source parametrization in terms of MS intensity.) The direct application of this complex model is complicated, so that reproducibility of results is impaired, as it often happens for recent PSHA studies. In fact, to address these issues, a relatively easy-to-implement weighted average grid-seismicity model was derived from the ensemble of the logic tree, by the same working group (details are provided in Chioccarelli et al. 2021). It is made of a grid of about eleven thousand point-like seismic sources covering the whole national territory and the surrounding areas. For each point, seismicity is defined in terms of annual rate of earthquakes per moment magnitude bins (i.e., the activity rates), forty-six in number, the width of which is set to 0.1. The central value of the lowest magnitude interval is M4.5 across all Italy. The highest magnitude bin, with rate larger than zero, is centered at M9 for about the 85% of the sources, and at M8.3 for the others. Finally, the probabilistic distribution of the style-of-faulting is defined for each point-like source. Hereafter, this source model is referred to as MPS19.
To develop PSHA in terms of MS intensity using the MPS19 model, which is in terms of magnitude, the propagation model of Pasolini et al. (2008) was chosen. It allows to convert the seismicity rates for magnitude classes into seismicity rates in terms of epicentral intensity I e . This means converting the magnitude bin central value into I e according to the semi-empirical relationship of Albarello et al. (2007), and the rate of the magnitude bin is attributed to the obtained (converted) intensity. The conversion model considers I e as a continuous variable, thus, for each bin, the converted value was limited to twelve; i.e., the largest possible MCS grade. The propagation model also assumes that the MS intensity at the grid point of interest, in terms of MCS, is normally distributed, conditional on the epicentral distance and I e . It follows that the MS intensity is treated as a continuous RV, although it is, by definition, discrete. More specifically, the relationship between MS, I e and R is given by Eq. (2): where ⋅ is a normally distributed RV, with zero mean and standard deviation that is equal to 0.87 and represents the residual of MS at the considered point, given R (in km) and I e . In the study presented herein, the model was applied within its definition range of epicentral distances, that is, 0-300 km.
As shown, the model's covariates are I e and R , which means that the model does not allow to account in the analysis for the style of faulting, as other propagation models for Italy do (e.g., Gomez Capera 2007). However, to be employed within the PSHA, these models require the input in terms of I 0 (in lieu of I e ), that, as mentioned, is not available for the source model used in this study.
In Fig. 1, the MPS19 model is represented in terms of annual rates of I e larger than four values. (It should be noted that associating a value of I e to an earthquake with the epicenter in the sea could appear meaningless according to the definition of I e . However, in the context of PSHA, I e is a magnitude-equivalent measure required to apply Eq. (2)).
The rates tend to be larger in north-eastern Italy, along the Apennine mountain chain and in eastern Sicily. It can be observed that the largest annual rate corresponding to I e ≥ 6 (2) 1a) is about 2E−03 events per year, whereas it is about 8E−04 events per year for I e ≥ 7 (Fig. 1b). Considering I e ≥ 9 (Fig. 1c) and I e ≥ 11 (Fig. 1d), the largest rates are about 1.4E−04 and 2.9E−05 events per year, respectively. PSHA in terms of macroseismic intensity using this source characterization was performed via the REASSESS software (Chioccarelli et al. 2019), which implements the above discussed source and propagation model.

Results
To help the discussion of the results presented in this section, it is worth recalling that grades from I to V represent earthquakes that may be felt by human beings but not producing structural damage. The seismic damage, in macroseismic terms, observed on structures is taken into account by grades from VI onwards. Indeed, according to the definition of the MCS scale, degree VI (strong earthquake) corresponds to small cracks in buildings, whereas MCS VII (very strong) denotes significant cracks and falling chimneys. Grades VIII (severe) and IX (destroying) represent the partial and complete collapse of some buildings, respectively, whereas MCS X (completely destroying) refers to earthquakes that cause collapse of most buildings in the hit area. Grades XI and XII do not have differences in terms of the observed damage on structures, being both indicative of the destruction of entire urban settlements. It may be interesting to note that the maximum MS intensity reported so far in Italy is X-XI, according to the latest version of CPTI.
The IMPS19 hazard maps, computed considering a grid of about ten-thousands points and four return periods, that are 50, 475, 975 and 2475 yr, are reported in Fig. 2. The maps are represented using a continuous color scale, in accordance with the hypothesis of continuous RV adopted by the propagation model to describe the MS intensity in one seismic event. On the other hand, in order to comply with the discrete character of the MCS scale, Fig. 1 Maps of the annual rates of earthquakes for the grid-seismicity source model, considering epicentral intensities, in terms of I e , equal to or larger than 6 (a), 7 (b), 9 (c) and 11 (d) most of the previous works dealing with PSHA in MS intensity present hazard maps by rounding the intensity value to the closest integer (it is a common procedure in relevant literature; e.g., Gomez Capera et al. 2010). Herein, both the continuous and discrete scale are considered to present the results: more specifically, the interval for a certain grade, say y, contains the shades pertaining to the intensities between y − 0.5 and y + 0.5. It follows that the discrete values in each map of the figure are defined as the grades that, at the points of the considered grid, are observed or exceeded with the selected return period.
It can be observed that the pattern of each hazard map is somehow similar to that of the maps in Fig. 1, as expected. More specifically, for each return period, it is found that, apart from Sardinia, whose hazard is comparatively lower with respect to the inland Italy and Sicily, the grid points with the lowest MCS grade are in the northwest area of the country and part of Apulia (southeast Italy) and western Sicily; also, relatively moderate-to-high intensity grades are found in eastern Sicily, coastal areas and northeast Italy, whereas the locations with the largest MCS grade are along central and southern Apennine mountain chain.
For T r = 50 yr , the lowest grade across inland Italy (i.e., neglecting Sardinia) is equal to IV, even if the fraction of the country corresponding to such a grade is almost negligible. With reference to the same return period, MCS equal to or larger than V and VI is found in about the 30% and 60% of the territory, respectively. The largest grade for T r = 50 yr is MCS VII, and it pertains to a relatively small area, in central Italy, covering about 5% of the national territory. With reference to T r = 475 yr , the lowest MCS grade across the inland country is VI; however, it pertains to an area with negligible extension. The largest grade with 475 yr return period is MCS IX, and it is found in the most hazardous area of Italy; i.e., in central and southern Italy, covering about 15% of the considered grid points (taken individually). Looking at intermediate grades, MCS VII and VIII correspond to 30% and 44% of the points, respectively, including almost the whole northern Italy, coastal areas and Sicily. For T r = 975 yr , MCS VII corresponds to 13% of the grid, which refers to northern Italy. With reference to the same return period, grade equal to or larger than VIII and IX is associated to about 37% and 41% of the grid nodes, respectively, whereas the largest MCS is equal to X, yet it is found for a few points in southern Italy, covering about 5% of the grid. Finally, referring to T r = 2475 yr , grade equal to or larger than VIII, IX and X corresponds to about 27%, 40% and 26% of the grid points respectively, whereas the lowest across all the inland country is VII, yet it corresponds to a small area (i.e., lower than 1%).

Maps of the lowest return period associated to each MCS grade
As it was already anticipated, rounding the MS intensity level, for one return period, to the closest integer (i.e., the MCS grade) may arise an issue. In fact, at a given grid point, although the hazard result for one return period is different from that associated to another return period, being the MS intensity considered as a continuous RV according to Eq. (1), it may happen that they are both rounded to the same integer. Thus, it cannot be warranted that the MCS grade at one point of the grid monotonically increases with the increasing of T r . The consequence is that the hazard maps for two different return periods may show, at the same location, the same MCS grade. To give an example, Fig. 2b, c show that a fraction of northern Italy is classified as MCS VII for both T r = 475 yr and T r = 975 yr ; the same, yet in different regions of Italy, also happens for other MCS grades. In other words, the discrete character of the MCS scale implies that the association of the MS intensity for a given return period to one point is not unique, as happening for typical ground motion IMs instead. For this reason, it may be also useful to provide an alternative representation of hazard results with respect to the one in Fig. 2. Thus, Fig. 3 shows, for each grid point and MCS grade, the lowest associated return period: panels from (a) to (d) pertain to MCS from IV to VII, whereas panels from (e) to (h) refer to MCS from VIII to XII. The color scale in the figure is discrete because each map refers to one MCS grade and is intended to show, for any point in Italy, the lowest associated return period among the four considered so far. Moreover, if, for a given node of the grid, the T r associated to the MCS grade does The maps in the figure reveal that the earthquakes with MCS equal to or larger than IV and V (i.e., not damaging for structures) are relatively frequent across the whole inland Italy, being the lowest return period equal to or smaller than 50 yr anywhere. Conversely, those events classified as capable of destroying entire urban settlements, that is, grade XI and XII, are relatively rare, being their (lowest) return period even longer than 2475 yr anywhere in Italy. Results are less homogenous for MCS from VI to X. In fact, in each of the panels from (c) to (g), the lowest return period causing occurrence or exceedance of the grade declared in the map tends to be lower along the central and southern Apennine mountain chain: this is expected, being this area characterized by the largest seismic hazard, as discussed in the previous section. Accordingly, in the same area, destroying (IX) and completely destroying earthquakes (X) have return period equal to 475 and 2475 yr, respectively. The low-to-moderate seismic hazard regions are identified, in each panel, by larger return periods: to give an example, in northern Italy and coastal areas, MCS equal to or larger than X has a return period longer than 2475 yr.

Comparison with MPS04 MCS hazard map
As discussed, Gomez Capera et al. (2010) provide a national hazard map, in terms of MCS, considering a return period of 475 yr. Such a map, herein indicated as IMPS04, is considered as a reference. In fact, as recalled above, the PSHA study was developed considering the same input and calculation procedure used for the seismic hazard assessment of Italy, denoted as MPS04, that is at basis of the current Italian building code (Stucchi et al. 2011). It relies on a source model featuring thirty-six areal source zones and no background seismicity. Intensity attenuation was modelled by means of regional propagation models which allow to account for the predominant rupture mechanism of the sources of Meletti et al. (2008). A logic tree accounting for alternative models of catalogue completeness, seismic parameters of the sources and propagation models was also considered; herein, the median map resulting from the logic tree is considered. The hazard map, whose data were provided by C. Meletti (personal communication, April 2021) is shown in Fig. 4a, together with the source zones that it is based on. Its pattern is quite similar to that of the map based on the grid-seismicity source model, for the same return period, shown in Fig. 2b. This is even more evident when comparing results considering the discrete MCS scale, as shown in Fig. 4b. More specifically, the map shows the absolute differences between IMPS04 results, in terms of MCS, and those denoted as IMPS19. It is shown that, in a significant fraction of Italy, which extends for about 65% of the territory (Sardinia Island is not taken into account, as IMPS04 does not give hazard results for the region) and includes almost all the MPS04 areal source zones, the MCS grade from Fig. 4a is the same as in Fig. 2b. The map of the differences also reveals that the differences are limited to one grade (except part of the small Giglio Island off-shore Tuscany, where the difference goes up to two grades, as shown in the figure). More specifically, the results of IMPS04 are lower than those of IMPS19 in the areas outside the source zones of Meletti et al. (2008) (e.g., northern Italy and Apulia region) and also in small areas enclosed by the MPS04 zones found in northern and central Italy; overall, these areas cover about 33% of the country. Finally, in the remaining 2% of Italy, which includes the grid points enclosed by the MPS04 zones found in a small fraction of northwestern Italy and in western Sicily, IMPS04 results are larger than those of IMPS19. This also happens in the Etna's volcanic region, even if a more recent study dealing with the hazard assessment for this area can be actually found in Azzaro et al. (2016). Figure 5 quantifies the size of the area, in the two hazard maps, associated to each MCS grade, in both percentage (i.e., relatively to the Italian inland surface, computed as the number of grid points corresponding to the MCS grade over the total number of points; left axis), and absolute terms (i.e., square kilometers; right axis). It is shown that, for each grade, the corresponding fraction of Italy according to the PSHA results based on the two source models is quite comparable. More specifically, earthquakes with grade equal to or larger than IX are found in about 15% of the national territory in the case of IMPS19 (as also discussed before) and 12% in the case of IMPS04. The same can be stated for the moderate-to-high hazardous regions: indeed, grade VII is found in about 30% of Italy in both the IMPS04 and IMPS19 hazard maps, whereas MCS VIII pertains to 35% and 44%, respectively. A slightly more evident difference is found in the case of MCS VI. Such a grade corresponds to about 13% of the country according to IMPS04, whereas it approaches to zero in the case of IMPS19. One of the reasons of these differences may be that MPS04 sources do not cover the whole Italy (i.e., there's no background seismicity), as the MPS19 grid-seismicity model does. For example, MCS VI is found in northern Italy and part of Apulia region, outside MPS04 sources (see Fig. 4a), whereas, in the same areas, it is MSC VII according to IMPS19 (Fig. 4b). What discussed so far shows that the IMPS19 hazard map is in a good agreement with the IMPS04 counterpart although source models, referred seismic catalogues and some propagation models 1 (i.e., input models) at the base of the two maps are different. However, it is worthwhile noting that such an agreement between the two models considered herein is not found when the seismic hazard in terms of spectral pseudo-accelerations is compared. 2 It is also to note that the comparison herein discussed should not be seen as any kind of validation of the results, yet it is intended only as informative; i.e., it by no means aims at providing a judgment. It has been extensively discussed elsewhere that the evaluation of hazard maps requires to compare them against ground motion observations in the years, an issue that, for large return periods (i.e., rare events) is hard, to say the least (e.g., Iervolino 2013;Iervolino et al. 2017).

Conclusions
Macroseismic intensity describes earthquakes' effects on communities, structures and environment. It may be used right after an earthquake for rapid loss assessment or, for example, in scenario-based risk analyses aimed at the evaluation of mitigations actions or for emergency planning purposes. Also, it may serve to characterize large-scale structural vulnerability on an empirical basis, or to have proxies for the event magnitude and the ground shaking in those areas where the seismic monitoring network is not dense or even absent. Therefore, it may be useful to derive, by means of PSHA, countrywide probabilistic seismic hazard maps providing the MS intensity (discrete) grade which is observed or exceeded with a return period of interest.
With reference to Italy, at least three studies dealing with PSHA in macroseismic intensity at the national scale can be found in literature, the most recent of which dates more than a decade ago. They developed the map of the intensity grades, according to the Mercalli-Cancani-Sieberg scale, that have ten percent probability of being exceeded in fifty years.
Herein a MS intensity hazard assessment study for Italy based on a recent grid-seismicity model was presented. The seismicity parameters were derived based on a conversion relationship between moment magnitude and the expected epicentral intensity. Results, that were obtained via PSHA and refer to the MCS grades corresponding to four return periods, were presented and, for each MCS grade, the map of the minimum return period causing its occurrence or exceedance was also discussed. Finally, the 475 yr return period hazard map based on the grid-seismicity source model was compared to the analogous map provided by the reference PSHA study in MS intensity for Italy. The main conclusions are recalled in the following.
• Given the return period, the lowest MCS grades pertain to Sardinia Island, part of Apulia, western Sicily and northern Italy. In the most hazardous area, that is, along the 1 The propagation model of Pasolini et al. (2008), used in this study, is also employed within the logic tree used by Gomez Capera et al. (2010). 2 Hazard maps for two spectral accelerations with the same return periods considered in this study based on MPS19 and MPS04 source model can be found in Chioccarelli et al. (2021) and Iervolino et al. (2018), respectively.
Apennine mountain chain in central and southern Italy, MCS IX has a return period of 475 yr. • Considering T r = 50 yr , a MCS equal to or larger than VI corresponds to about 60% of the considered grid points (taken individually). In the case of T r = 475 yr , MCS VII and VIII correspond to about 30% and 44% of the points, respectively. Similar percentages, yet for the larger grades, were found for T r = 975 yr , being 37% for MCS VIII and 41% for MCS IX. As pertaining to T r = 2475 yr , MCS IX was still found for about 40% of the grid points, whereas it is 26% for MCS X. • Anywhere in Italy, earthquakes of grade equal to or larger than IV have a return period lower than 50 yr, whereas those of MCS intensity equal to or larger than XI have a return period longer than 2475 yr. Along the Apennine mountain chain, earthquakes with grade equal to or larger than IX and X occur, on average, every 475 and 2475 yr, respectively. • The comparison between IMPS04 results and those presented in this paper shows that the MCS hazard maps from the two studies are generally comparable, and equivalent in more than half of Italy. IMPS04 results are lower than those found in this study mostly in the areas outside the zones of the source model used in IMPS04, whereas they are comparatively larger in part of Sicily and northwest area.
Hazard maps discussed in this study are provided with the paper as supplemental material.