A Multilevel Procedure at Urban Scale to Assess the Vulnerability and the Exposure of Residential Masonry Buildings: The Case Study of Pordenone, Northeast Italy

: More than the 60% of the Italian residential building stock had already been built by 1974, when seismic codes were enforced on a minimal part of the country. Unreinforced masonry buildings represent most of that share, but they are typical for each region, in terms of both materials and structural conﬁgurations. The deﬁnition of ‘regional’, i.e., more speciﬁc, vulnerability and exposure models are required to improve existing forecast models. The research presents a new geographic information system (GIS)-based multilevel procedure for earthquake disaster prevention planning at urban scale; it includes multicriteria analysis, such as architectural types, structural vulnerability analysis, microzonation studies, and socio-economic aspects. The procedure has been applied to the municipality of Pordenone (PN), a district town of the Friuli–Venezia–Giulia region, in Northeast Italy. To assess the urban seismic risk, more than 5000 masonry residential buildings were investigated and common types within sub-municipal areas and exposure data were collected. Simpliﬁed mechanical analysis provided a ‘regional’ vulnerability model through typological fragility curves. The integration of results into GIS tool permitted the deﬁnition of cross-mapping among vulnerability, damage scenarios (conditional and unconditional) and exposure (seismic losses, casualties, impact), with respect to various earthquake intensities expected in the town. These results are presented at di ﬀ erent scales: from the single building, to submunicipal area and to the entire town.

Lastly, the implementation of this framework into a GIS tool permits the visualization of the developed seismic scenarios at different levels of the urban scale.
The procedure here proposed was applied on the masonry-built heritage of Pordenone district town, Friuli-Venezia-Giulia (Northeast Italy). This region has one of the highest seismic hazards in Italy (according to the Italian hazard map [71]). Moreover, as a pilot case study, the municipality of Pordenone represents a significant medium-sized city of about 50,000 people, with a small medieval old town around which the city developed.
The research was carried out taking into consideration the following steps ( Figure 1): (a) identification of civil building types of Pordenone; (b) update and calibration of regional mechanical fragility and exposure models (masonry-built heritage); (c) validation of a new GIS-based multilevel and multicriteria analysis for urban earthquake disaster prevention planning; (d) damage scenarios analysis (conditional and unconditional) and impacts (economic losses, accessible buildings, victims, etc.); (e) definition of local seismic risk indicators for the development of regional mitigation strategies.

Multicriteria and Multilevel Framework Design and Implementation in QGIS Environment
Urban scale seismic risk assessment requires multicriteria framework that well fits with the GIS environment, where different information can be layered and integrated, such as general urban data, census statistics, building types, structural vulnerability analysis, microzonation studies and socioeconomic aspects, fragility and exposure models and scenarios, etc. (Table 1).

Multicriteria and Multilevel Framework Design and Implementation in QGIS Environment
Urban scale seismic risk assessment requires multicriteria framework that well fits with the GIS environment, where different information can be layered and integrated, such as general urban data, census statistics, building types, structural vulnerability analysis, microzonation studies and socio-economic aspects, fragility and exposure models and scenarios, etc. (Table 1).
The GIS database permits an effective management of all information and the implementation of damage and loss estimation models (probability functions). To automate and optimize the procedure, all computations were performed inside the GIS environment through query, field calculator, array, buffer, join, etc. The analysis was carried out by mapping all outputs through the open-source software Quantum GIS 3.0.3 (QGIS), released by the Open Source Geospatial Foundation [72].
Then, the implementation of information identified by the CARTIS approach, allowed for the developing of accurate fragility and exposure models for seismic risk predictions at urban scale. The first (2014) and second level (2016) CARTIS forms ( Table 2) both permit an expeditious survey of building types within the town districts (first level) and a detailed investigation of vulnerability of residential buildings (second level).
The new GIS-based procedure presents a flexible structure based on four different levels of detail (building, census unit, district, city), to provide local and global assessment of the urban seismic risk. This multilevel approach aims at identifying buildings, streets, and district areas with high level of risk based on the evaluation of various criteria. By mapping working at various administrative levels, it represents a tool for a better management and interpretation of urban data, which can be useful to allocate the resources for territorial seismic risk mitigation plans [73] (Table 3) (Figure 2). Table 1. Multicriteria framework for QGIS-based urban scale seismic risk assessment.

General urban data
Google Satellite, CTR (Regional technical map), urban plans, territorial information, census statistics (ISTAT), geotechnical studies, etc.

The Municipality of Pordenone: Urban Expansion and Conformation
The Friulian territory is characterized by high seismic activity [71]. Located in the Northeast of the country, Friuli-Venezia-Giulia region borders with Austria and the Carnic Dolomites (North), Slovenia (East), the Adriatic Sea (Southeast) and Veneto (Southwest). This area was struck by the dramatic earthquake of 1976 that caused about 17,000 collapses, 200,000 displaced, 965 victims and 3,000 injured [74]. In these areas, high seismic hazard combines with the vulnerability and exposure of the existing building stock, thus defining an ideal case for the validation and calibration of the procedure herein presented.
Pordenone is a medium-sized town of about 50,000 inhabitants and 10,000 buildings (9171 residential [41]), located in the southwest of the region (low Friulian plain, south of the Carnic Prealps), and crossed by the Noncello river. Its consolidated old town developed through a main axis in north-south direction and is composed of arcade buildings with 3-4 floors of low medieval origin. The building fabric of this area is characterized by masonry clusters that underwent frequent transformations over the last century, which were generally aimed at adapting the buildings to commercial uses.
Since the 16th century, the town has begun to expand in all directions and small urban villages arose around the old district, such as: S. Giorgio, S. Antonio, S. Giuliano and Trinità, Borgo Colonna, Torre, Borgo Meduna, S. Gregorio, Villanova, Noncello, and Rorai Grande.

The Municipality of Pordenone: Urban Expansion and Conformation
The Friulian territory is characterized by high seismic activity [71]. Located in the Northeast of the country, Friuli-Venezia-Giulia region borders with Austria and the Carnic Dolomites (North), Slovenia (East), the Adriatic Sea (Southeast) and Veneto (Southwest). This area was struck by the dramatic earthquake of 1976 that caused about 17,000 collapses, 200,000 displaced, 965 victims and 3000 injured [74]. In these areas, high seismic hazard combines with the vulnerability and exposure of the existing building stock, thus defining an ideal case for the validation and calibration of the procedure herein presented.
Pordenone is a medium-sized town of about 50,000 inhabitants and 10,000 buildings (9171 residential [41]), located in the southwest of the region (low Friulian plain, south of the Carnic Prealps), and crossed by the Noncello river. Its consolidated old town developed through a main axis in north-south direction and is composed of arcade buildings with 3-4 floors of low medieval origin. The building fabric of this area is characterized by masonry clusters that underwent frequent transformations over the last century, which were generally aimed at adapting the buildings to commercial uses.
Since the 16th century, the town has begun to expand in all directions and small urban villages arose around the old district, such as: S. Giorgio, S. Antonio, S. Giuliano and Trinità, Borgo Colonna, Torre, Borgo Meduna, S. Gregorio, Villanova, Noncello, and Rorai Grande.
In the 19th century, the industrial development of the town led to the construction of many cotton and paper mills in its suburbs. These factories gradually became as many small urban centers, characterized by an irregular growth of low-density residential masonry buildings for the working class. After the Second World War, the process of saturating the empty spaces between the town and the suburbs led the various urban built up areas to merge into a vast organically unitary city. In this regard, Pordenone transformed from a 'big village' into a real town: the buildings involved in this expansion were mostly brickwork houses from 2 to 5 floors and high-density reinforced concrete buildings. Lastly, in recent decades, new constructions have been limited, but extension works on existing buildings are still frequent.

Civil Building Types Characterization-CARTIS Form Survey
A georeferenced (QGIS) multilevel database was created to assess the vulnerability of building types at urban scale. The creation of specific attribute table (shapefile) with the fields of the CARTIS forms permitted an efficient and rapid insertion of the data into the folder and the integration of different layers (vector, raster), which are quite useful when dealing with this type of analysis: Google Satellite, CTR, historical cartographies, seismic microzonation maps, ShakeMaps, etc.
Secondly, according to the CARTIS approach, the municipality of Pordenone was divided into nine sections (districts), i.e., areas characterized by homogeneity of the built fabric through aspects, such as age of original installation, construction, and structural techniques ( Figure 3 and Table 4).
Heritage 2020, 3 FOR PEER REVIEW 8 In the 19th century, the industrial development of the town led to the construction of many cotton and paper mills in its suburbs. These factories gradually became as many small urban centers, characterized by an irregular growth of low-density residential masonry buildings for the working class. After the Second World War, the process of saturating the empty spaces between the town and the suburbs led the various urban built up areas to merge into a vast organically unitary city. In this regard, Pordenone transformed from a 'big village' into a real town: the buildings involved in this expansion were mostly brickwork houses from 2 to 5 floors and high-density reinforced concrete buildings. Lastly, in recent decades, new constructions have been limited, but extension works on existing buildings are still frequent.

Civil Building Types Characterization-CARTIS Form Survey
A georeferenced (QGIS) multilevel database was created to assess the vulnerability of building types at urban scale. The creation of specific attribute table (shapefile) with the fields of the CARTIS forms permitted an efficient and rapid insertion of the data into the folder and the integration of different layers (vector, raster), which are quite useful when dealing with this type of analysis: Google Satellite, CTR, historical cartographies, seismic microzonation maps, ShakeMaps, etc.
Secondly, according to the CARTIS approach, the municipality of Pordenone was divided into nine sections (districts), i.e., areas characterized by homogeneity of the built fabric through aspects, such as age of original installation, construction, and structural techniques ( Figure 3 and Table 4).   In total, five masonry-built residential types were recognized (MUR1-district 1, MUR1-districts 3-9, MUR2, MUR3, MUR4), depending on significant parameters, such as the construction period, the number of floors and the characteristics of vertical and horizontal structural components.
According to census data [41], Pordenone counts 3532 load-bearing masonry buildings (39%), 3760 r.c. structures (41%) and 1879 other types (20%). Table 5 shows the masonry types according to the construction periods (pre-1919; 1919-1945; 1946-1970; post-1970). It is noted that buildings built before 1919 or at most in the early 1900s prevail in the central zones (districts 1 and 2) [75]. They are typically arranged in clusters, with more or less regular texture, from rough stones and poor mortar (MUR2) to ashlars and clay bricks with lime mortar (MUR1), and with traditional timber floors and roofs (i.e., flexible) ( Table 6). The most common types of buildings in the suburbs (from district 3 to district 9) are the ones built in 1946-1970, made of clay brick walls and precast r.c. horizontal structures, such as precast joists ('Varese' type) or lightweight hollow clay bricks with steel rebars ('SAP' type) (both are semirigid types). This stock was divided into two classes: low-rise buildings with 3 floors or less (MUR3), which mostly represent the single-family residential buildings in the municipal area; and high-rise buildings with more than three floors (MUR4), that is the typical for popular housing of those years [76] (Figure 4).

Statistics and Results of Typological Study (GIS-Based Inventory)
The application of the second level CARTIS form (2016) to about 1000 buildings in Pordenone permitted to define the masonry types described above (Section 3.2) and the following extension of type codes to each building within the municipal area (GIS layer-CARTIS buildings).
The GIS layers 'Census Sections' and 'districts' contain results for the 554 sections and the nine districts of Pordenone, expressed in terms of amount of buildings, percentage of building types, and average of floors.
The results of the typological analysis show that in the central districts (i.e., districts 1 and 2) pre-1945 types (MUR1-MUR2) prevail (72%), while from district 3 to district 9 there is a strong increase of the MUR3 type (on average 70%). Moreover, the MUR1 type is around 30% of the masonry buildings in districts 5, 7, 9, i.e., suburbs having several historical clusters developed along the main old roads ( Figures 5 and 6).
The multilevel inventory provided the automatic calculation (fields calculator tool) and the visualization of the results through maps and graphics. As an example, Figures 7-9 present the maps of the old town (district 1) and suburbs (district 3) in terms of number of building for each CARTIS type, for the three levels of the GIS database are reported.

Statistics and Results of Typological Study (GIS-based Inventory)
The application of the second level CARTIS form (2016) to about 1000 buildings in Pordenone permitted to define the masonry types described above (Section 3.2) and the following extension of type codes to each building within the municipal area (GIS layer-CARTIS buildings).
The GIS layers 'Census Sections' and 'districts' contain results for the 554 sections and the nine districts of Pordenone, expressed in terms of amount of buildings, percentage of building types, and average of floors.
The results of the typological analysis show that in the central districts (i.e., districts 1 and 2) pre-1945 types (MUR1-MUR2) prevail (72%), while from district 3 to district 9 there is a strong increase of the MUR3 type (on average 70%). Moreover, the MUR1 type is around 30% of the masonry buildings in districts 5, 7, 9, i.e., suburbs having several historical clusters developed along the main old roads (Figures 5 and 6).
The multilevel inventory provided the automatic calculation (fields calculator tool) and the visualization of the results through maps and graphics. As an example, Figures 7-9 present the maps of the old town (district 1) and suburbs (district 3) in terms of number of building for each CARTIS type, for the three levels of the GIS database are reported. The application of the second level CARTIS form (2016) to about 1000 buildings in Pordenone permitted to define the masonry types described above (Section 3.2) and the following extension of type codes to each building within the municipal area (GIS layer-CARTIS buildings).
The GIS layers 'Census Sections' and 'districts' contain results for the 554 sections and the nine districts of Pordenone, expressed in terms of amount of buildings, percentage of building types, and average of floors.
The results of the typological analysis show that in the central districts (i.e., districts 1 and 2) pre-1945 types (MUR1-MUR2) prevail (72%), while from district 3 to district 9 there is a strong increase of the MUR3 type (on average 70%). Moreover, the MUR1 type is around 30% of the masonry buildings in districts 5, 7, 9, i.e., suburbs having several historical clusters developed along the main old roads (Figures 5 and 6).
The multilevel inventory provided the automatic calculation (fields calculator tool) and the visualization of the results through maps and graphics. As an example, Figures 7-9 present the maps of the old town (district 1) and suburbs (district 3) in terms of number of building for each CARTIS type, for the three levels of the GIS database are reported.

The Old Town subtypes
The old town of Pordenone (district 1) develops from the ancient town hall through a main axis in the north-south direction along which more than 800 structural units shape either simple or complex clusters. The complexity of this area required a specific typological study to be carried out, to conveniently detail the database better.
The information of 850 US (Structural Units) (33 nonresidential; 73 R.C. 624 MUR1; 7 MUR2; 102 MUR3; 11 MUR4) identified within the district were inserted in the GIS database. Moreover, new subtypes of the MUR1 and MUR2 (district 1) were defined for these clusters, to improve representativeness of the vulnerability classification and the following seismic risk analysis ( Figure  10).
Units and subtypes were defined according to the ground plan of the old town buildings, by taking into consideration various factors, such as surface, plan and elevation regularity, aggregation, presence of arcades, and warping of floors. Four subtypes were identified inside the old town (Table  7), as follows: • MUR1-T1: <150 m² surface; pre-1900 age; close rectangular shape; unit in connection; loadbearing masonry in x-y direction; presence of arcades; timber floors and roof.
• MUR1-T3: 80-300 m² surface; pre-1900 age; large sized quadrangular units with several internal partitions or with the presence of isolated masonry pillars; load-bearing masonry; timber or semirigid floors and roof.

The Old Town Subtypes
The old town of Pordenone (district 1) develops from the ancient town hall through a main axis in the north-south direction along which more than 800 structural units shape either simple or complex clusters. The complexity of this area required a specific typological study to be carried out, to conveniently detail the database better.
The information of 850 US (Structural Units) (33 nonresidential; 73 R.C. 624 MUR1; 7 MUR2; 102 MUR3; 11 MUR4) identified within the district were inserted in the GIS database. Moreover, new subtypes of the MUR1 and MUR2 (district 1) were defined for these clusters, to improve representativeness of the vulnerability classification and the following seismic risk analysis ( Figure 10).
Units and subtypes were defined according to the ground plan of the old town buildings, by taking into consideration various factors, such as surface, plan and elevation regularity, aggregation, presence of arcades, and warping of floors. Four subtypes were identified inside the old town (Table 7), as follows: • MUR1-T1: <150 m 2 surface; pre-1900 age; close rectangular shape; unit in connection; load-bearing masonry in x-y direction; presence of arcades; timber floors and roof. • MUR1-T2: <80 m 2 surface; pre-1900 age; small regular square shape with few internal partitions; load-bearing masonry; timber floors and roof.    The 48% of units belong to the MUR1-T1 subtype, characterized by terraced arcaded buildings; these are generally the most precious buildings of the old town, often showing frescoed facades. The 24% of the sample presents small regular inclusive units (MUR1-T2); the 15% is made up of large size quadrangular buildings (MUR1-T3) with an irregular internal distribution, as result of many aggregation processes; lastly, the 13% (MUR1-T4) corresponds to corner units of building clusters ( Figure 11). The 48% of units belong to the MUR1-T1 subtype, characterized by terraced arcaded buildings; these are generally the most precious buildings of the old town, often showing frescoed facades. The 24% of the sample presents small regular inclusive units (MUR1-T2); the 15% is made up of large size quadrangular buildings (MUR1-T3) with an irregular internal distribution, as result of many aggregation processes; lastly, the 13% (MUR1-T4) corresponds to corner units of building clusters ( Figure 11). The typological surveys and the information gathered in the context of the CARTIS project were used to produce a local fragility model for the different types of masonry buildings of the town.

The Vulnus VB 4.0 Procedure
The calculation of fragility curves for the Pordenone masonry building stock was carried out by a mechanical procedure, which is based on simplified modeling applied to unreinforced masonry (URM) buildings through the software Vulnus VB 4.0. This permits vulnerability analyses based on few information about geometry, material properties, construction features and some qualitative data [34,50]. The simplified analysis performed by Vulnus VB 4.0 is mainly based on resistance checks and on the calculation of the possible collapse mechanisms that the various parts of the structure may suffer. This follows the kinematic approach for local verification, which is also suggested by the Italian legislation [77]. Specifically, in-plane (IP) and out-of-plane (OOP) analyses can be carried out through resistance and linear kinematic calculations, and the horizontal accelerations (a) that activate The typological surveys and the information gathered in the context of the CARTIS project were used to produce a local fragility model for the different types of masonry buildings of the town.

The Vulnus VB 4.0 Procedure
The calculation of fragility curves for the Pordenone masonry building stock was carried out by a mechanical procedure, which is based on simplified modeling applied to unreinforced masonry (URM) buildings through the software Vulnus VB 4.0. This permits vulnerability analyses based on few information about geometry, material properties, construction features and some qualitative data [34,50]. The simplified analysis performed by Vulnus VB 4.0 is mainly based on resistance checks and on the calculation of the possible collapse mechanisms that the various parts of the structure may suffer. This follows the kinematic approach for local verification, which is also suggested by the Italian legislation [77]. Specifically, in-plane (IP) and out-of-plane (OOP) analyses can be carried out through resistance and linear kinematic calculations, and the horizontal accelerations (a) that activate the main mechanisms can be registered. Then, based on the acceleration (a) normalized to gravity g (a/g), the indexes I1 and I2 are computed, for the IP and OOP resistances, respectively.
In more detail, I1 represents the shear resistance of the building normalized to its total weight, evaluated in its weak direction as the sum of the shear strengths of parallel wall systems analyzed in their average plane as rigidly coupled. In the case of irregularities in plan and elevation, this index is corrected, to consider the non-uniform distributions of normal and tangential stresses. The I2 index, instead, is a more complex parameter. Vulnus VB 4.0 outputs the value of the triggering accelerations of the possible OOP mechanisms associated with the vertical (i.e., tilting of the overall walls and tilting and flexural collapse of the walls at the top story) and horizontal (i.e., bending and arching effect with failure at the top story, tilting and flexural collapse of the arch effect with failure of shoulders and detachment of the transverse walls, both at the top story) masonry portions of each wall, and combines them to obtain only one index (I2) that can be representative of the OOP behavior of the whole building.
All other relevant vulnerability aspects that are not directly computable (i.e., qualitative information regarding the types of floors, roof and foundations, the configuration and regularity of the building, the state of preservation, the quality of connections and construction details, etc.) are taken in consideration through a further index, I3. It is calculated as a weighted average (from 0 to 1, where 0 represents a building designed according to anti-seismic regulations) of the scores assigned to the parameters identified by the second level GNDT form.
Lastly, using these three indices and applying the fuzzy set theory [51], three fragility curves are calculated by the software. They represent the central probability (White curve) and two extreme probabilities (Upper-and Lower-Bounds curves) of exceeding a moderate-severe damage state, associated to a DS2-3 (damage scale according to EMS-98), as a function of the expected peak ground acceleration (PGA).
The association of the fragility curve calculated by Vulnus VB 4.0 with a damage DS2-3 is a reasonable assumption considering that the triggering acceleration of a certain mechanism, assessed by linear analysis, is a necessary condition for the mechanism activation (DS2, corresponding to a Limited Damage-LD, according to [78]), can be very close to the maximum system capacity (DS3, corresponding to a Significant Damage-SD), but it is not yet a sufficient condition to turn the mechanism into a partial collapse (DS4, corresponding to a Near Collapse-NC).

Application of Vulnus VB 4.0 to the Sample of Buildings
The mechanical analysis of seismic vulnerability was carried out on a sample of 60 masonry buildings, previously detected through the CARTIS second level form (2016) and selected for their representativeness and completeness of available data. Table 8 shows the corresponding types ascribed to the buildings. Plans and main elevations were available for all the buildings. In addition, for buildings built from 1946, sections and some information about the materials and construction techniques were also provided. Design documents were generally available for more recent buildings. These aspects affect the quality of the information available in the GNDT form, which roughly follows the quality of the second level CARTIS one. Therefore, the reliability of information is higher with regard to the geometric aspects, whereas it is less accurate for the materials properties (especially in older buildings) and the quality of the connections between structural elements, since no results of onsite investigation tests were available. Table 9 shows the average values of the three indices provided by Vulnus VB 4.0 for the types of the sample of 60 buildings analyzed here. The MUR2 type is the most vulnerable, having the lowest I1 and I2 indexes and the higher I3 one, mainly due to the poor quality of the masonry. With regard to the subtypes of the old town, the most vulnerable is the subtype T4, which is characterized by irregular plan and large surface area, while the subtype T2 is the least vulnerable, due to the small size and regularity of the plan (simple quadrangular geometry). Lastly, for MUR3 and MUR4 types, there is a reduction in vulnerability, due to the progress of construction techniques over time, but not as significant as expected. However, the behavior of MUR3 and MUR4 types is strongly influenced by the presence of reduced thickness clay brick walls (25-30 cm), highly vulnerable horizontal structures (e.g., Varese and SAP types), and either poor materials or structural details. All these aspects strongly influenced the results in terms of both vulnerability indices and fragility curves. The Vulnus VB 4.0 procedure also provides the fragility curves that define the probability of exceeding a certain level of damage with respect to the PGA. Fragility curves are processed by points: therefore, the first step was to fit them to obtain cumulative log-normal distribution curves, by finding the mean and the standard deviation ( Figure 12).
For each of the identified masonry types, i.e., the four subtypes of the old town MUR1-T1/T2/T3/T4, and the remaining MUR2, MUR3, MUR4, the mean curves for White (Figure 12), Upper-Bounds and Lower-Bounds were obtained. Figure 13 shows the White fragility curve of all the building types: the most vulnerable type is MUR2, characterized by a poor masonry quality; MUR1-T2 is less vulnerable, whereas among the others (MUR-T1, MUR-T3 and MUR-T4) there are no substantial differences, since the structural details used in these historical buildings are the same. MUR3 and MUR4 types, which represent more recent buildings, have overall lower vulnerability but still respectively comparable to the MUR1-T2 and MUR1-T3 subtypes, due to the onset of the above-mentioned aspects that affect the seismic behavior. Moreover, MUR4 is slightly more vulnerable than MUR3; this is probably due to the presence of buildings that are taller, less regular, and with a larger surface area than the others are.     Figure 14 shows the complete model, with the White, Upper-and Lower-Bound fragility curves, for each masonry type. As it can be seen, the dispersion range defined by the Upper-and Lower-Bound curves increases for the type MUR3 according to its decrease of vulnerability. This feature is due to the reduction in the Lower-Bound fragility related to a better seismic behavior.   Figure 14 shows the complete model, with the White, Upper-and Lower-Bound fragility curves, for each masonry type. As it can be seen, the dispersion range defined by the Upper-and Lower-Bound curves increases for the type MUR3 according to its decrease of vulnerability. This feature is due to the reduction in the Lower-Bound fragility related to a better seismic behavior.   Figure 14 shows the complete model, with the White, Upper-and Lower-Bound fragility curves, for each masonry type. As it can be seen, the dispersion range defined by the Upper-and Lower-Bound curves increases for the type MUR3 according to its decrease of vulnerability. This feature is due to the reduction in the Lower-Bound fragility related to a better seismic behavior.

Extension and Calibration of Fragility Curves
The fragility curves obtained by Vulnus VB 4.0 analysis for the described sample were defined for a damage level DS2-3, i.e., 'moderate-severe', according to the classification provided by the EMS-98 [55]. To have a complete picture and, therefore, to calculate seismic damage scenarios for the ordinary built at urban scale, it was necessary to extend the curves to all five damage levels DS1-DS5.
A reference model, derived from the conversion in PGA of the model of Lagomarsino and Cattari [36] through the law of correlation according to Margottini et al. [79], was calibrated on the DS2-3 fragility curves. The following steps [54] were applied: 1.
for each vulnerability class (from A-high to F-low) of the macroseismic model a fragility curve corresponding to a DS2-3 damage grade was calculated, to be compared with the mechanical ones obtained through the Vulnus VB 4.0 approach; 2.
for each mechanical fragility curve, an optimal linear combination between two curves (DS2-3) of the macroseismic model was made, thanks to the genetic algorithm NSGA-II (i.e., Non-dominated Sorting Genetic Algorithm [54]), aimed at minimizing the absolute and relative errors; 3. the combination coefficients associated to the different classes of vulnerability were used to generate a further set of fragility curves, associated to the five levels of damage (DS1-DS5) for the building types here analyzed. Table 10 reports the results of that calibration, i.e., the percentages of linear combination of the classes of the macroseismic model, able to obtain the set of fragility curves for the five levels of damage. According to the EMS-98 classification, unreinforced masonry buildings (URM) belong to the first four classes (A-D), based on type of structure from rubble stone (A) to massive stone (C) or URM with RC floors (C-D). Nevertheless, it should be noted that solid or hollow clay brick masonry (with or without mortar repointing) which characterize the MUR1, MUR3, and MUR4 types are not included among the specific cases considered by the EMS-98 [6].
With regard to the 'White' combinations, almost all building types (except for MUR2), belong to the C-D range. Consistently, the Upper-Bounds curves are shifted to the left to the B-C range, whereas the Lower Bounds curves are more displaced to the right to D class. The comparison between the Vulnus VB 4.0 mechanical curves (local mechanical-heuristic model) and the calibrated macroseismic fragility curves DS2-3 shows some differences in terms of dispersion, due to their different nature (the first derived in a mechanical way and the second obtained by expert judgment). The standard deviation of the calibrated curves is, in fact, higher, thus reflecting the greater dispersion typical of the macroseismic approach and consistent with empirical observations on the seismic behavior of pseudo-plasticity of masonry structures, but different from that of mechanical curves. This effect is more evident with Lower-Bounds ( Figure 15).

LUW Fragility Model
The further extension to a territorial scale requires a model valid at a statistical level for the whole built stock examined. In such a connection, fragility curves with greater dispersion than the deterministic case (in accordance with empirical observations related to damage detected because of past seismic events) can be used.
To obtain a single model that best describes the behavior of the different types of buildings on a territorial scale, the three heuristic-mechanical models can be combined, so that the main information provided by the White, Lower and Upper Bounds fragility curves can be included.
The curves of the three models compose a single model called LUW (Lower-Upper-White), built as follows [54]: • between 0% and 2.5% of the White probability, the LUW is considered totally as the Upper-

LUW Fragility Model
The further extension to a territorial scale requires a model valid at a statistical level for the whole built stock examined. In such a connection, fragility curves with greater dispersion than the deterministic case (in accordance with empirical observations related to damage detected because of past seismic events) can be used.
To obtain a single model that best describes the behavior of the different types of buildings on a territorial scale, the three heuristic-mechanical models can be combined, so that the main information provided by the White, Lower and Upper Bounds fragility curves can be included.
The curves of the three models compose a single model called LUW (Lower-Upper-White), built as follows [54]: • between 0% and 2.5% of the White probability, the LUW is considered totally as the Upper-Bounds curve; • between 2.5% and 50% of the White probability, the LUW is assumed as a linear combination of the Upper-Bounds (from 100% to 0%) and the White (from 0% to 100%); • between 50% and 97.5% of the White probability, the LUW is obtained combining linearly the White (from 100% to 0%) and the Lower-Bounds (from 0% to 100%); • above 97.5% of the White probability, the LUW is given at 100% by the Lower-Bounds curve.
The reason underlying these choices was to increase the dispersion of the mechanical model, while still maintaining the mean value of the White one. Table 11 lists the parameters that define the curves thus obtained for all the building types and for the five levels of damage; Figure 16 shows the graphs of the LUW sets.
As already observed, the most vulnerable type for ordinary masonry buildings in Pordenone is MUR2. In general, for equal PGA and damage grade, the highest damage probability is obtained for the oldest buildings designed without seismic regulation. With regard to the four historical subtypes, the most vulnerable is MUR1-T1, with 20-40% that reach DS5, while the least vulnerable is MUR1-T2; MUR-T3 curves have a greater slope, while MUR-T4 reaches a greater percentage for DS5 level and has smaller intervals between one curve and another. Lastly, as for the most recent types MUR3 and MUR4 it can be seen that the first one is the least vulnerable, whereas MUR4 has higher slopes for lower damage levels, and shorter intervals, thus reaching higher probability for DS4 and DS5 damage levels.

Seismic Damage Maps
The fragility curves presented in Section 4 define a vulnerability model for the ordinary masonry buildings of Pordenone. The final step of this work would be the combination of this vulnerability with the exposure assessed through the CARTIS approach (Section 3) and the seismic hazard of the town, in order to evaluate seismic damage and risk, through GIS maps.
The damage maps represent the probability of reaching or overcoming a certain damage state, thus providing information about which buildings, areas or districts are more prone to seismic damage. Furthermore, risk maps can show the possible economic losses in terms of different indicators, such as repair or reconstruction costs, number of unusable buildings/dwellings, victims, injuries, and homeless people.

Seismic Damage Maps
The fragility curves presented in Section 4 define a vulnerability model for the ordinary masonry buildings of Pordenone. The final step of this work would be the combination of this vulnerability with the exposure assessed through the CARTIS approach (Section 3) and the seismic hazard of the town, in order to evaluate seismic damage and risk, through GIS maps.
The damage maps represent the probability of reaching or overcoming a certain damage state, thus providing information about which buildings, areas or districts are more prone to seismic damage. Furthermore, risk maps can show the possible economic losses in terms of different indicators, such as repair or reconstruction costs, number of unusable buildings/dwellings, victims, injuries, and homeless people.
This information can be very important for carrying out targeted investigations and defining priority criteria for seismic retrofit interventions. Specifically, these tools can be used for prevention and mitigation of seismic risk by authorities and institutions to manage resources in the aftermath of an earthquake and to select effective emergency measures and evacuation plans.

Conditional Damage
Conditional damage expresses the expected damage for a specific ground motion (e.g., PGA). The significant values of PGA for Pordenone were obtained from the parameters given by the Italian seismic code [80][81][82], i.e.,: values of a g (maximum horizontal acceleration), FO (maximum value of the amplification factor for the horizontal acceleration spectrum), and Tc* (reference value for determining the beginning of the plateau in the horizontal acceleration spectrum). Starting from a mesh of 10751 points that cover the whole Italian territory and computing nine return periods (Tr = 30, 50, 72, 101,  140, 201, 475, 975, 2475 years), the values expected for Pordenone were selected (Table 12). Table 12. Values of maximum horizontal acceleration (a g ), maximum value of the amplification factor for horizontal acceleration spectrum (FO), reference value for determining beginning of plateau in horizontal acceleration spectrum (Tc*) for Pordenone, according to Italian code [80]. These values were combined with the parameters of the morphological and stratigraphic characteristics that determine the local response (Table 13). Table 13. Description of soil types according to Italian code [80].

Soil Type Description
Soil A rocky or very rigid soils Soil B soft rocks and deposits of very dense coarse-grained soils or very consistent fine-grained soils Soil C deposits of medium dense coarse-grained soils or medium consistent fine-grained soils, with a depth of more than 30 m Soil D deposits of poorly dense coarse-grained soils or poorly consistent fine-grained soils, with a depth of more than 30 m Soil E soils with characteristics like those defined for categories C or D, with a depth not exceeding 30 m According to the Geological Report for Pordenone [83] and to the seismic microzonation study [84], Pordenone refers to three soil categories: the northern area can be associated with a type C soil (mostly gravel with depth of a few tens of meters); the historical center and the area south to the center to a type D (predominance of sands and clays); the remaining southernmost part of the town to a type E.
These variations can activate site effects of seismic action amplification and different PGAs for each area, or district, or even for each single ISTAT section. Now that the seismic hazard has been defined for all the different areas, i.e., every point in Pordenone can be associated with a particular value of PGA, it is possible to estimate the damage and, consequently, create damage maps. This can be done by associating to each residential masonry type the fragility results provided by the fragility set (Section 4) evaluated for that predicted PGA. Results are given in terms of probability of occurrence of the different damage states, and can be done for all the possible return periods Tr.
In this work, conditional damage maps for a Tr = 475 years were created, as reference for the Life Protection Limit State (SLV) (see Table 12). Figure 17 shows the results obtained for all the damage states, from DS1 (slight damage) to DS5 (collapse).
Heritage 2020, 3 FOR PEER REVIEW 23 In this work, conditional damage maps for a Tr = 475 years were created, as reference for the Life Protection Limit State (SLV) (see Table 12). Figure 17 shows the results obtained for all the damage states, from DS1 (slight damage) to DS5 (collapse).  Table 14 shows the conditional damage results (DS1-DS5) for the CARTIS district 1 (Tr=475 and soil C). The probability of exceedance of the five damage grades was calculated for each CARTIS building type and subtype of the old town. The DS1 (slight damage) values exceed 75% for all the sample, the DS3 (moderate-severe damage) reaches the highest percentages (more than 40%) for MUR1-T1 and MUR2 types, while the DS5 (collapse) is on average around 4-5%, except for the most vulnerable type (22%).  Table 14 shows the conditional damage results (DS1-DS5) for the CARTIS district 1 (Tr = 475 and soil C). The probability of exceedance of the five damage grades was calculated for each CARTIS building type and subtype of the old town. The DS1 (slight damage) values exceed 75% for all the sample, the DS3 (moderate-severe damage) reaches the highest percentages (more than 40%) for MUR1-T1 and MUR2 types, while the DS5 (collapse) is on average around 4-5%, except for the most vulnerable type (22%). This procedure can be made not only at a single-building level, but also considering wider areas (such as ISTAT sections or even CARTIS districts). A mean value of probability of reaching a specific damage state can be calculated with respect to the vulnerability, but also the number and the extension of buildings inside the area of interest. As an example, Figure 18 shows the results obtained for the damage state DS3, thus increasing the scale of the area from single buildings to ISTAT sections and, finally, to CARTIS districts.

Unconditional Damage
Unconditional damage represents the combination of multiple levels of ground motion (for various Tr values), taking into account the annual probability of reaching those levels. When computing unconditional damage, an observation time window is chosen, and all the possible

Unconditional Damage
Unconditional damage represents the combination of multiple levels of ground motion (for various Tr values), taking into account the annual probability of reaching those levels. When computing unconditional damage, an observation time window is chosen, and all the possible earthquake scenarios that can occur in the selected time are taken into consideration. Every scenario must be included with its own probability of occurrence in the observation time, as follows: where p is the probability that an earthquake with return period Tr occurs in the observation time To. In this work, observation times To of 10 and 50 years were applied. Tables 15 and 16 show the unconditional damage results for both years and soil type C. As for the conditional maps, different urban level can be represented, such as single buildings, ISTAT sections and CARTIS districts. Figure 19 shows the unconditional damage maps for a damage state DS3, for both To of 10 and 50 years, and for the three urban levels mentioned above.

Seismic Risk Maps
The damage maps combine seismic hazard and vulnerability. By adding information about exposure, such as number of buildings, reconstruction or other costs, and number of people, it is possible to produce risk maps that graphically represent losses and impact indicators, such as fatalities, injuries, economic losses, and impact on the usability of the building stock.
A framework for the calculation of the seismic risk is here proposed. It can be performed for each ISTAT section, and subsequently for each CARTIS district detected in Pordenone.

Seismic Risk Maps
The damage maps combine seismic hazard and vulnerability. By adding information about exposure, such as number of buildings, reconstruction or other costs, and number of people, it is possible to produce risk maps that graphically represent losses and impact indicators, such as fatalities, injuries, economic losses, and impact on the usability of the building stock.
A framework for the calculation of the seismic risk is here proposed. It can be performed for each ISTAT section, and subsequently for each CARTIS district detected in Pordenone.
In seismic risk maps, each damage state can turn into a risk indicator, thanks to different damage-to-risk matrices [60]. Table 17, in its first section (a), shows the matrix in which each damage state is associated with a different percentage of cost of repair or replacement (with a minimum and maximum value for each damage state, in order to have a possible range of economic losses). Section b shows the matrix that associates damage states with the percentage of casualties (fatalities or injuries). Lastly, section c shows the percentages of useable, not usable in the short or long-time span, and collapsed buildings associated with the damage states. Table 17. Matrix for: a. economic losses (% of cost of repair or replacement-minimum and maximum set), b. casualties (% of fatalities and injuries) and c. impact on buildings (% of useable, not usable in short or long-time span, and collapsed buildings). These partial results can be combined by multiplying the probability of occurrence of each damage state for the percentages shown in the matrices. Therefore, it is possible to obtain a single value that represents the most likely percentage of cost, people, or buildings that are associated with each indicator, as shown in Table 17.

% of Damage Results
Then, exposure must be considered; for each ISTAT section, the following values were extracted from the GIS database: (i) number of buildings belonging to each MUR type; (ii) total floors area belonging to each MUR type; (iii) population associated with MUR buildings.
To calculate the economic losses, a value of 1350 EUR/m2 was chosen, according to Borzi et al. [60], as a reconstruction cost for ordinary masonry buildings. By multiplying this cost for the floor area in an ISTAT section, a total cost of the considered section can be found. By multiplying the percentage of losses by the total cost of the ISTAT section, an estimate of the economic losses can be obtained. The same calculation can be performed for casualties, multiplying the percentage of fatalities/injuries by the number of people, and also for impact, multiplying the percentage of usable, not usable or collapsed buildings by the total number of buildings in the ISTAT section. Figure 20 shows the general framework described above. The seismic risk expected losses/impact where calculated for observation times of 10 (Table 18) and 50 years (Table 19) for each district and for the whole municipality of Pordenone, and the associated maps were generated (Figures 21 and 22).
As for the damage maps, these analyses do not consider the entire built stock of the municipality, therefore they reflect a forecast of losses and impact for only the ordinary masonry-built stock. The seismic risk expected losses/impact where calculated for observation times of 10 (Table 18) and 50 years (Table 19) for each district and for the whole municipality of Pordenone, and the associated maps were generated (Figures 21 and 22).
As for the damage maps, these analyses do not consider the entire built stock of the municipality, therefore they reflect a forecast of losses and impact for only the ordinary masonry-built stock.

Conclusions
The seismic risk assessment of urban centers is an important issue due to the great vulnerability of masonry buildings and the frequent conditions of deterioration that characterize residential built

Conclusions
The seismic risk assessment of urban centers is an important issue due to the great vulnerability of masonry buildings and the frequent conditions of deterioration that characterize residential built

Conclusions
The seismic risk assessment of urban centers is an important issue due to the great vulnerability of masonry buildings and the frequent conditions of deterioration that characterize residential built heritage. When dealing with large scale assessments, simplified procedures based on expert judgement are required, since complete and detailed information for this type of analysis is often missing. Nevertheless, these procedures are generally based on census data that provide only few and poor information, to carry out accurate vulnerability assessments and seismic risk management of urban areas. Therefore, the updating of existing methods is suggested in order to: (a) create local inventory of building types (taxonomy); (b) develop mechanical fragility models for typical masonry-built types (with deep investigations to improve prediction models); (c) calculate damage scenarios (conditional and unconditional); (d) estimate the seismic risk of building types (losses and impact).
In this study, a new GIS-based multicriteria and multilevel procedure for the assessment of seismic risk of masonry-built types was presented. The typological and structural characterization of buildings (exposure model) and the following mechanical vulnerability analysis (fragility model) allowed for developing seismic risk maps. The results were produced with an automated process, thanks to the implementation of all the aspect mentioned before into a georeferenced and layered database of masonry-built types.
The procedure was applied to the case study of the municipality of Pordenone (Northeast Italy), focusing on ordinary load-bearing masonry buildings. The town was subdivided into nine districts, and various residential types were categorized by parameters, such as age of construction, number of floors and surface plan, aggregation, masonry quality, and type of horizontal structures.
All the collected data were implemented into a multilevel urban scale georeferenced (GIS) inventory of masonry building types, collecting specific information related to seismic vulnerability. Starting from these data it was possible to carry out mechanical vulnerability analyses according to building types configurations and derive the fragility sets for the eight types (MUR1-T1, MUR1-T2, MUR1-T3, MUR1-T4, and MUR2, MUR3, MUR4) of buildings.
Furthermore, the fragility models were produced using the Donà et al. [54] procedure, that consists of the calibration of a macroseismic fragility model on a mechanical one, derived by simplified analyses of masonry buildings by Vulnus VB 4.0. The specific fragility model obtained for Pordenone, highlights how the building types belonging to the most recent construction periods are less vulnerable, thanks to their better structural characteristics. Although buildings in the old town were the first ones to be built on the city, they are not the most vulnerable since they feature good quality structural details. Actually, the most vulnerable buildings are the ones belonging to MUR2 (pre-1919 and poor materials buildings).
Starting from this fragility model, a framework for damage and risk estimates was calculated. For what concerns damage assessment, both conditional damage (for a return period Tr = 475 years) and unconditional damage (for observation times To = 10 and 50 years) were calculated.
Furthermore, economic losses, casualties, and impact were evaluated. This procedure was implemented for single buildings, ISTAT sections, and CARTIS districts. In the next 50 years, an average economic loss of about EUR 250 million is expected for the town, as well as 22 fatalities, 28 injuries, and 777 unusable buildings (in a short time span, long time span, or in terms of collapses, respectively).
This pilot study provided the basis for the definition of a fragility and exposure 'regional' model for Northeast Italy, which can be taken as an example for further applications to other historical centres of the surrounding area. These applications prove that the methodology proposed can lead to practical outcomes. These results might be used by institutions not only to locate the most vulnerable areas, but also to know where to prioritize anti-seismic interventions, thus attaining a more targeted risk mitigation at different urban levels.