Project number: 265138 Project name: New ... - MATRIX - KIT

1 downloads 0 Views 3MB Size Report
volcano (e.g., Minakami, 1974; Shimizu et al., 1992; Wolfe, 1992; Legrand et al., ...... Spence, R.J.S., Kelman, I., Baxter, P.J., Zuccaro G., S. Petrazzuoli (2005), ...
Project number: Project name: Project acronym: Theme:

Start date:

265138 New methodologies for multi-hazard and multi-risk assessment methods for Europe MATRIX ENV.2010.6.1.3.4 Multi-risk evaluation and mitigation strategies

01.10.2010

End date: 31.12.2013 (39 months)

Deliverable: D7.3: Naples test case Version: Draft Responsible partner: AMRA Month due: M38 Month delivered: Primary authors(a): Alexander Garcia-Aristizabal Angela Di Ruocco Warner Marzocchi _______________________________ Signature (a)

12.2013 _________________ Date

see the Acknowledgements section for a full list of people contributing to this document

Reviewer: Kevin Fleming _______________________________ Signature

12.2013 _________________ Date

Authorised: Kevin Fleming _______________________________ Signature

12.2013 _________________ Date

Dissemination Level PU Public Restricted to other programme participants (including the Commission PP Services) Restricted to a group specified by the consortium (including the Commission RE Services) CO Confidential, only for members of the consortium (including the Commission Services)

X

1

2

Abstract

Considering only the natural sources of environmental hazards for the Naples city area, the main hazards that can be identified are (1) Volcanic hazards, which encompasses the active volcanoes of Somma-Vesuvius, Phlegrean Fields and the island of Ischia; Seismic hazards, due both to active faulting in the Campanian Apennines and to earthquakes arising from the three active volcanoes in the Naples area; Hydrogeological hazards, mainly due to flash floods, pyroclastic flows and rock falls due to slope instability; and Tsunami hazards, whose origin may be associated with earthquakes with epicentres in the sea or in coastal areas, coastal or submarine landslides, and volcanic eruptions (in this case, possible sources are volcanic activity on Ischia and Vesuvius, and submarine volcanoes in the low Tyrrhenian Sea). After collecting the data for the analyses, the Arenella was the target area in which it was possible to retrieve high resolution data, hence its selection as the pilot area for the analyses in the project. In this report, a summary of all the activities performed for the Naples test case within the framework of the MATRIX project is presented. It therefore constitutes a collection of the main results of different deliverables produced in the project. In particular, section 2 is a summary of the single hazard and risk assessments performed. Section 3 presents the results of the assessment of interactions and cascading effects. Finally, section 4 summarizes the results of interactions with local stakeholders in Naples.

Keywords: Naples test case, MATRIX project, single hazard and risk assessment, cascading effects and interactions, Social and institutional barriers for effective decision making using multi-risk data

3

4

Acknowledgments

The research leading to these results has received funding from the European Commission’s Seventh Framework Programme [FP7/2007-2013] under grant agreement n° 265138. Different partners from the MATRIX project have collaborated in the preparation of this report. Below we present a list of all people who participated in this document’s production (alphabetical order of institutions): - AMRA: Alexander Garcia-Aristizabal, Angela Di Ruocco, Warner Marzocchi, Jacopo Selva. - IIASA: Anna Scolobig.

5

6

Table of contents 1

Introduction ................................................................................................................. 12

2

Single hazards and risk assessment ........................................................................... 16 2.1

Volcanic ash-fall hazard and risk ......................................................................... 16

2.1.1

Nodes describing the volcanic hazard (Nodes 1 to 8) .................................... 16

2.1.2 Nodes describing Damage: Node D0 (inventory), node D1 (physical vulnerabilities) and Node D3 (Exposure) ..................................................................... 17 2.2 3

Seismic hazard and risk ...................................................................................... 22

Multi-hazard assessment for cascading effects Naples: considering interactions ....... 25 3.1

Identification of scenarios for the Naples test case.............................................. 25

3.2

Scenarios of cascading effects analysed in the Naples test case........................ 27

3.2.1 Volcanic earthquakes and seismic swarms triggered by volcanic unrest (scenario 1) ................................................................................................................. 27 Scenario description................................................................................................ 27 Kinds of interactions ................................................................................................ 28 Temporal aspects.................................................................................................... 28 Analysis ................................................................................................................... 29 Example considering the Campi Flegrei area ...................................................... 30 3.2.2 Scenario 2: Simultaneous occurrence of ash-fall and earthquake hazards (interaction at vulnerability level) ................................................................................. 34 Scenario description................................................................................................ 34 Kind of interactions.................................................................................................. 34 Temporal aspects.................................................................................................... 35 Analysis ................................................................................................................... 35 4 Interactions with stakeholders in Naples: Identification of social and institutional barriers to effective decision-making and governance in the case of multiple hazards ...... 44 5

Summary ..................................................................................................................... 46

The Naples test case provided the opportunity to test different methodologies developed in the MATRIX project in specific activities of single- and multi-risk assessments. In particular, we can summarize the following main results: .................................................. 46  In WP2, a procedure based on Bayesian analyses has been developed and applied to the Naples test case in order to calculate volcanic risk (considering ash fall). It is worth noting that the methodology proposed is general enough to be easily implemented for other kind of risks; .............................................................................................................. 46  In WP3, a conceptual framework to analyse interactions and cascade effects was developed. This framework was implemented in the Naples test case allowing us to

7

analyse the effects of interactions between two risk sources (seismic and volcanic) at different levels. ................................................................................................................... 46  The results of the interaction analysis and multi-risk assessment for this particular scenario were used in the WP6 for the identification of social and institutional barriers to effective decision-making and governance when considering multiple hazards. ............... 46 References ........................................................................................................................ 47

8

List of Figures Figure 1. Study area for the Naples test case. ................................................................... 12 Figure 2. Map of the target areas: Camaldoli, Arenella, and Furigrotta ............................. 14 Figure 3. Components of risk defining the event tree used in this study. ........................... 16 Figure 4. Hazard function for the Arenella district. The red area represents a 80% confidence interval (within the 10th and 90th percentiles). The red line represents the ‘best guess’ values (50th percentile) ................................................................................... 18 Figure 5. Fragility functions for all building typologies (TA k). The red area represents a 80% confidence interval. .................................................................................................... 21 Figure 6. (preliminary) risk curve, for an exposure time of 1 year. The red area represents a 80% confidence interval. ................................................................................................. 22 Figure 7. Seismic hazard curve for the Neapolitan area, obtained from the National seismic hazard map (Gruppo di Lavoro, 2004). The solid line represents the median value, whereas the dashed lines are the 16th and 84th percentiles and serve as an uncertainty measure for the results. ..................................................................................................... 23 Figure 8. Loss exceedance probability (a) and (mean) expected total loss (b) resulting from the seismic risk analysis for the Arenella area (Naples) .................................................... 24 Figure 9. Grid of points in which the PGA predictions have been performed following the methodology described before ........................................................................................... 31 Figure 10. Example of the PGA predictions for a single M2.2 event in the Campi Flegrei area using the methodology described in the previous paragraphs. .................................. 32 Figure 11. Hazard curve for a grid point located far from the seismic swarm source area (about 15 to the NE)........................................................................................................... 33 Figure 12. Hazard curve for a grid point located close to the Arenella area (one of the test areas defined in the MATRIX project). ............................................................................... 33 Figure 13. Hazard curve for a grid point located in an area with a high density of epicentral solutions............................................................................................................................. 34 Figure 14. Footprints of residential buildings in the Arenella area (Naples), with the most likely seismic vulnerability classes defined for each structure............................................ 37 Figure 15. Example curves used in the SAVE methodology for the calculation of fragility functions (in this example, for building class B and for damage levels D2 to D5). ............. 39 Figure 16. Histogram showing the cost distribution of rebuilding residential buildings in Arenella area (produced using a sample of 800 buildings out of about 1200 residential building in the Arenella area). ............................................................................................ 39 Figure 17. (a) Loss exceedance probability, and (b) Cumulative distribution function of the expected loss resulting from the seismic risk analysis for the Arenella area. Note that the plot in (b) represents the variability of the mean value as a result of the uncertainties considered. ........................................................................................................................ 40 Figure 18. Loss exceedance probability (a and c) and (mean) expected total loss (b and d) resulting from a risk analysis considering the combined effects of seismic and two different values of ash load: 3 kPa (a and b) and 5 kPa (c and d) for the Arenella area. ................. 41

9

Figure 19. Summary of results of the exceedance probabilities of loss for the three considered scenarios: (1) seismic, (2) seismic+3kPa ash load, and (3) seismic+5kPa ash load. ................................................................................................................................... 42 Figure 20. Summary of results of the expected losses (mean values) for the three considered scenarios: (1) seismic, (2) seismic+3kPa ash load, and (3) seismic+5kPa ash load. ................................................................................................................................... 42

10

List of Tables Table 1. Summary of some relevant parameters and data from the three selected areas of interest in Naples for multi-type hazard and risk analysis: Canaldoli, Arenella, and Fuorigrotta (LS – landslide and rockfalls, EQ – earthquakes, VE – volcanic eruptions, WF – wildfires).......................................................................................................................... 14 Table 2. Application’s taxonomy TA, from Spence et al. (2005) ........................................ 18 Table 3. Node D0, best estimate of p(TAi|TA(ISTAT)) ........................................................... 19 Table 4. Taxonomy relative to FM2, from Spence et al. (2005) ......................................... 19 Table 5. Node D1, association between input FMs and application’s TA .......................... 20 Table 6. Node D2, alpha values of (posterior) Dirichlet distribution. .................................. 21 Table 7. Possible cascade scenarios for the Naples test case (from MATRIX D3.3). ........ 26 Table 8. Time scales of interest for the volcanic earthquake hazard produced by triggered seismicity during volcanic unrests ...................................................................................... 29 Table 9. Time scales of interest for the cascading effect considering the simultaneous occurrence of volcanic ash-fall and earthquake hazards ................................................... 35 Table 10. Seismic-building structures classification (from Zuccaro et al., 2008) ................ 36 Table 11. Classification of damage to masonry buildings, according to the European Macroseismic Scale (EMS98, Grünthal, 1998) .................................................................. 38 Table 12. Seismic intensity increment for corresponding ash loads (from Zuccaro et al., 2008).................................................................................................................................. 39

11

1 Introduction Naples is located in south west Italy on a bay along the Tyrrhenian coast in the inner-most part of the Gulf of Naples. It is in the midst of two volcanic areas, the Phlegrean Fields and Mount Vesuvius. From the eastern hills of the Phlegrean Fields, the city slopes down northward to the Campanian Plain and south-eastward to the slopes of Mt. Vesuvius. The eastern area of Naples consists of a wide floodplain (the Sebeto river hollow) while the western urban area of Naples consists of the Bagnoli Plain and the remains of 4 ancient volcanic hollows generated by the Phlegrean Fields volcanism (Chiaia, Neapolis, Fuorigrotta and Soccavo, see Figure 1). The urban area of Naples is mainly composed of hills, most of them with heights ranging from 150 to 452 meters.

Figure 1. Study area for the Naples test case.

The geological history of the whole territory is complex: the latest substrate is mainly composed of volcanic debris. The Campania region is fragile and is exposed to several environmental risks, both natural and anthropogenic, that have in recent centuries caused disastrous events resulting in thousands of victims, massive loss of property and huge damage to the cultural heritage. Considering only the latter half of the 20th century, in Campania there were about 3300 victims of natural disasters (e.g., 2734 in the 1980 Irpinia earthquake, 40 in the 1944 eruption of Vesuvius and the rest due mainly to landslides and floods). In fact, Campania is the Italian region that has suffered the highest number of fatalities from natural catastrophes: more than 40% of the total for Italy. In particular, considering only the natural sources of environmental hazards for the Naples city area, the main hazards are:  Volcanic hazard, which encompasses the active volcanoes of Somma-Vesuvius, Phlegrean Fields and the island of Ischia; 12

 Seismic hazard, due both to active faulting in the Campanian Apennines and to earthquakes arising from the three active volcanoes in the Naples area;  Hydrogeological hazard, mainly due to flash floods, mudflows and rock falls due to slope instability;  Tsunami hazard, whose origin may be associated with earthquakes whose epicentres are in the sea or in coastal areas, coastal or submarine landslides, and volcanic eruptions (in this case, possible sources are volcanic activity on Ischia and Vesuvius, and submarine volcanoes in the low Tyrrhenian Sea).  Forest fires, which in the last few years have been reported in the city of Naples, mainly during the summer period. For example, more than 40 forest fires were registered in 2011 (Department of Agriculture, Campania Region 2012, and MATRIX deliverable D6.3 “Social and institutional barriers to effective multi-hazard decision making”). While the volcanic and seismic hazards can affect the whole city (with an intensity varying from place to place), the areas prone to the tsunami and hydrogeological hazards are respectively located on the coast and close to unstable slopes. In particular, in the Phlegrean Fields, the sectors mainly susceptible to landslides or rock falls are located along steep tuff slopes, in pyroclastic deposits and, in some places, in lavas. The recurrent morphologies are associated with the caldera rims, slopes of the coastal fault and the escarpments bordering river valleys. Considering the area of the city of Naples, the parts susceptible to landslides and rock falls are mainly located near the outcrops of Campanian Ignimbrite and Neapolitan Yellow Tuff on the Camaldoli hill (Soccavo and Pianura side). Similar events of lower intensity have been observed on the Posillipo hill, which is mainly composed of Neapolitan Yellow Tuff, along the slopes on the side of the Fuorigrotta district and along the tuff cliffs of Coroglio. In the MATRIX project, three target areas within the Naples city area have been selected for multi-type hazard and risk analysis. These areas have been selected according to the spatial distribution of natural hazards and the available information about the elements at risk. The selected target areas are: Camaldoli, Fuorigrotta and Arenella (Figure 2). Some characteristics of the three areas, namely their areas, numbers of inhabitants, numbers of residential and commercial buildings, etc., are summarized in Table 1.

13

Figure 2. Map of the target areas: Camaldoli, Arenella, and Furigrotta Table 1. Summary of some relevant parameters and data from the three selected areas of interest in Naples for multi-type hazard and risk analysis: Canaldoli, Arenella, and Fuorigrotta (LS – landslide and rockfalls, EQ – earthquakes, VE – volcanic eruptions, WF – wildfires).

Camaldoli 7.16 13,498

NAPLES Arenella Fuorigrotta 2.90 6.61 46,367 52,781

496,597,400 426

2,037,583,600 1,534,917,800 4,069,098,800 1,371 2,593 4,390

census section

census section

census section

census section

Hazard Resolution Hazards

LS, EQ, WF

EQ, VE

EQ, LS, VE

EQ, LS, VE, WF

# Total apartments # inhabitants/flat

5633 2.4

20237 2.3

17670 3.0

economic value/flat (k€)

88.2

100.7

86.9

Area (km^2) Residential population (number) Economic Value (€) Number of residential Buildings Exposure resolution

Total 16.67 112,646

2-3 inhabitants for flat circa 100 k€ for flat

Considering the scenario collection presented in the MATRIX MS3 report, the different hazards of interest for each area of interest are: - Camaldoli: Landslides, earthquakes and forest fires. - Arenella: Earthquakes and volcanic ash fall. - Fuorigrotta:

14

Earthquakes and volcanic hazards (such as vent opening, ash fall, and pyroclastic flows). After collecting the data for the analyses, the Arenella was selected as the target area, as it was possible to retrieve high resolution data. This report summarizes all the activities performed in the Naples test case within the framework of the MATRIX project. First, a summary of the results dealing with single hazard and risk assessment is presented. Then, section 3 presents the results of the assessment of interactions and cascading effects and finally, in section 4, the results of interactions with local stakeholders are discussed.

15

2 Single hazards and risk assessment In this section we will consider only the natural sources of environmental hazards for the Naples city area, as outlined in the previous section. In the framework of MATRIX project, for the test case in Naples two hazard sources have been considered for quantitative hazard and risk analyses: volcanic ash-fall and earthquakes. 2.1 Volcanic ash-fall hazard and risk The volcanic ash-fall hazard and risk assessment performed in the MATRIX project in the Naples test case was performed within the framework of WP2 “Single type risk assessment and comparability” and presented in the Deliverable D2.2 “Uncertainty Quantification”. This involved a quantitative application of a Bayesian event tree approach developed for the quantification and propagation of uncertainties in the case of volcanic risk assessment. For simplicity, the event tree is presented divided according to two general components of risk, namely hazard and damage, each of which can be regarded as a sub-event tree. The sub-event tree for the volcanic hazard assessment corresponds to the procedure described in the BET_VH model (Marzocchi et al., 2010, Selva et al., 2010), whereas the sub-event tree for the Damage is composed of three elements (nodes), relative to inventory definition (that here we call node D0), fragility functions (i.e., physical vulnerability, node D1), and losses (node D2). For simplicity, we describe the two components independently, but together they form a more general event tree (Figure 3):

Figure 3. Components of risk defining the event tree used in this study.

2.1.1 Nodes describing the volcanic hazard (Nodes 1 to 8) The Bayesian event tree approach for volcanic hazard assessment (BET_VH) is a methodology already described in the literature (e.g., Marzocchi et al., 2010, Selva et al., 2010) and is the procedure used as the reference for the volcanic hazard assessment in this document. A brief description is provided in the following paragraphs, but for a more detailed description the reader is invited to consult the original papers. Let’s assume that the function [θk] represents the probability density function of the random variable ,describing the probability of the jth event at the kth node. The distribution [θk] has to be unimodal and with domain [0,1] because the random variable is a probability. A suitable distribution with these requirements is the Beta distribution (e.g., Marzocchi et al., 2004). The key point of this strategy is that the available empirical and theoretical information can modify the mean and variance of this probability distribution. In general, we can consider all of the available information determines the mean of the distribution (which we can consider to be the aleatory uncertainty), while the variance of the distribution, can be considered as the epistemic uncertainty, is modified (usually 16

significantly reduced) by having more data available. The PDF associated with is calculated by following two steps (based on a common procedure of Bayesian analysis, e.g., Gelman et al., 1995): first, an “initial” (prior) distribution is defined (as a Beta/Dirichlet distribution, depending if the node has two or more possible outcomes) taking into account all the theoretical beliefs; then, this information is “updated” by considering past (observed) data (encoded as a Binomial/multinomial distribution) using Bayes’ theorem. Detailed descriptions of each node can be found in the Deliverable D2.2 and in the referenced works. 2.1.2 Nodes describing Damage: Node D0 (inventory), node D1 (physical vulnerabilities) and Node D3 (Exposure) The treatment of epistemic uncertainty in damage estimates is here done using the methodology whose development is still in progress within the on-going Italian project ByMuR (Bayesian Multi-Risk, ByMuR 2010-2013; Selva et al. 2012a). The method, for each typology of physical structure, can be summarized in the following logical steps:  Definition of the application’s Taxonomy (TA), Intensity Measure (IM), Damage states (DSs) and losses (L), referred to as the “application metrics”. In this description “application” is referred to the specific case analysed in the problem.  Collection of applicable Fragility Models (FMs), to be mapped in the application metrics (TA, IM and DSs). A weight to each model is assigned, based on its applicability in the target area, or due to its credibility as defined by expert opinion;  Collection of the input data for the target area (Inventory), to be mapped into the application metrics (TA, L). Note that the first step is the harmonization of intensity measures, and the definition of the typologies of structures and relative damage states. In practice, this means that it should be defined one common reference metric for the whole application, hereinafter referred to as [IM, TA, DS, L], in order to prepare a common base for merging different vulnerability/losses models. This process may imply a drastic reduction of classes of exposed elements with respect to single-model applications. This fact is in common to all processes of comparison among models. For example, the comparison among FMs performed in the EU-FP7 project SYNER-G1 (SYNER-G 2010-2012) implied the use of one reference IM (for example Peak Ground Acceleration, PGA, for seismic hazard), the conversion of all models to this IM, the definition of only two damage states (“yielding” and “collapse”), and the definition of correlations between the original and target (common) taxonomy (Crowley et al. 2011). The Bayesian assessment of damages, as described above in general terms, may be undertaken using an Event Tree logic. For each asset in the target area, we define the following nodes (for details, see D2.2):  Node D0 “Inventory”: classification of assets according to the TA.  Node D1 “Fragility”: assessment of DS, for each typology in the TA and each IM.  Node D2 “Losses”: assessment of L, for each DS of each asset. The treatment of uncertainties based on a Bayesian event tree approach is applied to quantify the ash-fall volcanic hazard and risk in one of the target areas defined for the Naples test case in MATRIX (Arenella, for details see, e.g., MATRIX Deliverables D2.2 for the methodology and D3.3 “Scenarios of cascade events” for the test case definition). The resulting risk curves are computed through Monte-Carlo simulations, combining the 1

http://www.vce.at/SYNER-G/ 17

conditional probabilities at all the nodes of the Event Tree, separately for each IM level (scenario). The exercise developed in D2.2 is restricted to the assessment of tangible economic losses in residential structures, in particular, to direct losses associated with physical damages to the roofs due to the accumulation of volcanic ash. The results obtained for this example are summarized as follows. Regarding the hazard assessment phase, volcanic ash-fall hazard assessment is performed, for the sake of simplicity, considering as the source just a single vent opening (node 4) in the Campi Flegrei caldera (about 6 km eastward of the Arenella district). All the other nodes are set according to Selva et al. 2012b (nodes 1 to 3), Orsi et al. 2009 (node 5), and Selva et al. 2010 (nodes 6 to 8). Note that an extension for a full hazard assessment considering any possible vent opening within the volcanic caldera does not pose any particular conceptual problem. The results of the hazard assessment, in the form of a hazard curve representing exceedance probabilities, are presented in Figure 4, where the intensity measure IM is ash loading, in kPa.

Figure 4. Hazard function for the Arenella district. The red area represents a 80% confidence interval (within the 10th and 90th th percentiles). The red line represents the ‘best guess’ values (50 percentile)

The inventory is based on the last available census data (ISTAT, 2001). The 1202 residential buildings in the study area were classified in terms of 3 categories of structures (R.C., Masonry and Mixed) and 7 categories of age (pre 1919,1919-1945,19461961,1962-1971,1972-1981,1982-1991, after 1991). Therefore, we have a total of 21 ISTAT typologies. The target taxonomy TA is defined as in Table 2 (Spence et al., 2005). We consider 4 classes (WE, ME, MS, ST), distinguished by their roof typology. Only one damage state DS is considered, that is, the collapse of the roof. Table 2. Application’s taxonomy TA, from Spence et al. (2005)

Roof typology (TAi) WE (weak) MW (medium weak) MS (medium strong)

ST (strong)

Description Sheet roofs, old or in poor condition. Tiled roof, old or in poor condition. Masonry vaulted roof. Sheet roof on timber; average quality; average or good quality tiled roof on timber rafters or trusses. Flat reinforced concrete roof; sloping reinforced concrete roof. Sheet roof on timber rafters or trusses, good quality and condition, designed for cyclone areas. Flat reinforced concrete roof designed for access; recent, good quality construction, younger than 20 years. 18

At node D0, each of the ISTAT typologies should be mapped in the target taxonomy TA, through the probability p(TA|TA(ISTAT)). The best estimate of this probability, for each ISTAT typology, is based on a historical consideration of the building practices in Naples, and is reported in Table 3. The maximum allowed epistemic uncertainty is considered, setting the “equivalent number of data” λ = 1 (e.g., Marzocchi et al., 2008). No independent data (e.g., survey of actual classification) are available, so that this prior distribution represents our best model. At node D1, we consider two different fragility models presented in Spence et al. 2005. The first model (FM1) is based on 5 typologies (RT1, RT2, RT3, RT4 and RT5), and it is specific for buildings in the Vesuvius area. The taxonomy of this model is presented in Table 4. The second model (FM2) is more generic, and it is the basis for the application’s TA (WE, ME, MS, ST). Given the more specific nature of FM1 to the Naples area, we assign to the FM1 model a larger credibility/applicability, that is, its weight w1 = 0.7, with respect to w2 = 0.3 assigned to FM2. (ISTAT)

Table 3. Node D0, best estimate of p(TAi|TA

)

R.C.

PRE 1919 1919-1945 1946-1961 1962-1971 1972-1981 1982-1991 1991-2001

WE

1

1

0.34

0

0

0

0

MW

0

0

0.33

0.5

0.2

0

0

MS

0

0

0.33

0.5

0.8

1

0

ST

0

0

0

0

0

0

1

Masonry

PRE 1919 1919-1945 1946-1961 1962-1971 1972-1981 1982-1991 1991-2001

WE MW

0.34 0.33

0.34 0.33

0.25 0.25

0.15 0.3

0.15 0.3

0 0.2

0 0

MS

0.33

0.33

0.25

0.3

0.55

0.8

0

ST

0

0

0.25

0.25

0

0

1

Mixed

PRE 1919 1919-1945 1946-1961 1962-1971 1972-1981 1982-1991 1991-2001

WE MW

0.34 0.33

0.34 0.33

0.25 0.25

0.15 0.3

0.15 0.3

0 0.2

0 0

MS

0.33

0.33

0.25

0.3

0.55

0.8

0

ST

0

0

0.25

0.25

0

0

1

Table 4. Taxonomy relative to FM2, from Spence et al. (2005)

Roof typology (TAi) RT1 RT2 RT3 RT4 RT5

Description Wooden pitched roof and flat iron R. C. SAP R. C. older than 20 years R. C. younger than 20 years or steel beam with r.c. slab R. C. pitched roof

In Table 5 we report how each fragility curve for each model (FM1 and FM2) is mapped into the application’s TA: in each cell, it is reported how representative in the Arenella district each starting typology is, for each target typology (TAi). 19

Table 5. Node D1, association between input FMs and application’s TA

FM1

FM2

RT1 RT2 RT3 RT4 RT5 WE MW MS ST

Application’s TA WE MW 0 0.5 0 0.5 0 0 0 0 0 0 1 0 0 1 0 0 0 0

MS 0 0 1 0 0 0 0 1 0

ST 0 0 0 1 0 0 0 0 1

Note that, since the same taxonomy is adopted for the application and for FM2, the relative matrix is an identity matrix. For FM1, different values are reported: these values depend on both (i) engineering considerations, that is, the definition of typologies in both taxonomies, and (ii) knowledge of the general features of the building stock in the target area. For example, note that RT5 is, due to engineering considerations, in the class TAi=”ST”, but its weight is 0, since this typology of buildings is surely not present in the Arenella area. The final weight for each fragility for each TA is obtained as the product of the FMs weight (0.7 and 0.3) and the applicability weight of each curve. The prior model is finally obtained from the FM1 and FM2 curves through a weighted procedure (Selva et al. 2012a). As for D0, no further data are available, and thus the prior model represents our best estimation. The resulting fragility curves are reported in Fig. 5 for each TA.

20

Figure 5. Fragility functions for all building typologies (TAk). The red area represents a 80% confidence interval.

At node D2, the model relative to losses is assumed equal for all. Since for all typologies of buildings the eventual reconstruction of the roof would be of the same type (RC roof) for all typologies, the total loss depends only on the size of each roof. However, the total area of the roof for each building is not known, but only for a subset of buildings (~ 60). Thus, the probability p(loss|collapse) is estimated through Bayesian inference, with a noninformative uniform prior distribution (Dirichlet with parameter αk =1), and with 60 observations (modelled with a multinomial distribution). The parameters of the obtained posterior (Dirichlet distribution) are reported in Table 6 for the discrete number of bins centred in Lk. Since no links between the 60 observed buildings and the ISTAT database exist, this distribution is assumed to be common for all buildings. Table 6. Node D2, alpha values of (posterior) Dirichlet distribution.

Cost (Lk) (x106 €) Alphai

0.050

0.073

0.095

0.118

0.140

0.163

0.185

0.208

0.230

0.253

6

10

13

10

11

7

6

3

2

2

21

The resulting risk curves are computed through Monte-Carlo simulations, combining the conditional probabilities at all the nodes of the Event Tree, separately for each IM level (scenario). With regards to epistemic uncertainty, the probability relative to each node is sampled from each distribution, providing an explicit evaluation of the epistemic uncertainty at the node. More precisely, for each IM, it is assessed by the probability of exceedance, in the exposure time ΔT, of given threshold value of losses l. That is: [p(>l; IM, ΔT)]=[p(>l|IM)] [p(IM; ΔT)] where [x] indicates the distribution of x; p(IM; ΔT) is obtained from nodes 1 to 8 (hazard assessment), while p(>l|IM) from nodes D0 to D2 (damage assessment). The contributions of all IMs are then summed at several level of confidence (percentiles), providing a proxy of the total variability of risk assessment (Selva et al. 2012a). The resulting risk curve is assessed by sampling from the distribution at all nodes, set as described above, and it is reported in Figure 6. The best estimate is relative to the 50th percentiles, and the confidence interval is represented by the 10th and 90th percentiles.

Figure 6. (preliminary) risk curve, for an exposure time of 1 year. The red area represents a 80% confidence interval.

2.2 Seismic hazard and risk For the Naples test case, the seismic risk has been considered as part of the activities of one of the two interaction scenarios analysed as examples of the multi-hazard / multi-risk environment (for details see Section 3). For this reason, a detailed description of the process for the seismic risk assessment is presented in section 3.2.2, where the interaction scenario characterized by the simultaneous effect of earthquakes and volcanic ash-fall is presented. In this part, a brief description of the adopted process is presented. The seismic hazard information has been obtained from the Italian national seismic hazard map produced by the Istituto Nazionale di Geofisica e Vulcanologia - INGV (e.g., Gruppo di Lavoro, 2004; Meletti and Montaldo, 2007). The intensity measure for seismic hazard is available as both macro-seismic intensity, and in PGA (g). The method used by the INGV for the seismic hazard assessment is based on the standard method of Cornell (Cornell,

22

1968; Bender and Perkins, 1987), and it has been implemented using a logic tree in order to explore epistemic and aleatory uncertainties, which refer mainly to:  Completeness of the used catalogues;  Determinations of seismicity rates;  Attenuation relationships. The information can be freely obtained from the website: http://zonesismiche.mi.ingv.it An example of the hazard curve (in PGA) for the area of interest in Naples is shown in Figure 7.

Figure 7. Seismic hazard curve for the Neapolitan area, obtained from the National seismic hazard map (Gruppo di Lavoro, 2004). th th The solid line represents the median value, whereas the dashed lines are the 16 and 84 percentiles and serve as an uncertainty measure for the results.

For the damage functions, a taxonomy with four classes of structures (A, B, C, D) has been considered, following the structure classification defined in Zuccaro et al., 2008 (for details, see section 3.2.2). The building classes were selected considering the type of vertical and horizontal structure. For each class of building, five damage states (DS) -plus no damage- are considered, ranging from 0 (no damage) to 5 (collapse). The seismic fragilities P(DS|PGA), were calculated using the earthquake damage probability functions for the different structural typologies of Zuccaro et al., 2008, which are based on the SAVE methodology (developed in the SAVE project within the GNDT-DPC programme 20002002, Zuccaro, 2004). The SAVE model deals with the physical vulnerability of buildings and is based on an empirical method to be used with macroseismic intensity hazard maps. The fragility curves are in this case the functional representation of the probability of having a given damage level for different earthquake intensities, and for a given building class. After assigning a fraction of the total rebuilding costs to each damage state (0 to 5), an analysis of uncertainties and the setting up of a rebuild cost distribution was performed.

23

As mentioned before, a more detailed description of the process for the seismic risk assessment is presented in section 3.2.2, where the interaction scenario characterized by the simultaneous effect of earthquakes and volcanic ash-fall is presented. The quantification of the loss distribution and expected losses were performed following a risk quantification procedure in which uncertainties (if available) from the hazard, the fragilities, the inventory of exposed elements (i.e., for each building, the probability of belonging to each building class is assessed), and the rebuilding cost were propagated up to the final risk quantification. Figure 8(a) shows the results of the total loss distribution (in terms of loss exceedance probability) from the seismic risk assessment. The solid line represents the median (50th percentile) value of the Loss exceedance probabilities (as a “best guess” value), whereas the uncertainty is represented by the area bounded by the 16th and 84th percentiles. Likewise, it is possible to calculate the expected annual loss, where instead of a single value, in this case we obtain a probability distribution showing the variability of the expected loss (the mean) given the considered uncertainties. Figure 8(b) shows the expected loss distribution for the seismic case. (a)

(b)

Figure 8. Loss exceedance probability (a) and (mean) expected total loss (b) resulting from the seismic risk analysis for the Arenella area (Naples)

24

3 Multi-hazard assessment for cascading effects Naples: considering interactions The core of the probabilistic assessment of cascading effects in a multi-hazard problem consists of the identification of the possible interactions that are likely to occur and that may result in an amplification of the expected damages in a given area of interest. This concept is also the fundamental part of a holistic multi-risk analysis. Within the framework of Work Package 3 of the MATRIX project, two possible kinds of interactions have been described, namely: (1) interactions at the hazard level, in which the occurrence of a given initial ‘triggering’ event, entails a modification of the probability of occurrence of a secondary event, and (2) interaction at the vulnerability (or damage) level, in which the main interest is to assess the effects that the occurrence of one event (the first one occurring in time) may have on the response of the exposed elements to another event (which may be of the same kind as the former, but also a different type of hazard). Implicitly, a combination of both kinds of interactions is another possibility, in which case, in the discussion about interactions at the vulnerability level, both dependent and independent hazards have been considered as possible scenarios. A detailed description of the probabilistic framework to quantify both kinds of interaction can be found in deliverable D3.4 of the MATRIX project (Garcia-Aristizabal et al., 2013a). 3.1 Identification of scenarios for the Naples test case An important initial step towards the assessment of cascading effects is the identification of possible scenarios. The term “scenario” is used in a wide range of fields and so different interpretations can be found in practical applications. In general, a scenario may be considered as a synoptic, plausible and consistent representation of an event or series of actions and events (e.g., see MATRIX deliverable D3.3, Garcia-Aristizabal et al., 2013b). In particular, it must be plausible because it must fall within the limits of what might conceivably happen, and must be consistent in the sense that the combined logic used to construct a scenario must not have any built-in inconsistencies. The process of scenario identification is the first and fundamental step towards a quantitative multi-risk assessment that considers cascading effects. To achieve the required complete set of scenarios, different strategies can be adopted, ranging from event-tree-like to fault-tree-like strategies. In many applications, an adaptive method combining both kinds of approaches is applied in order to ensure an exhaustive exploration of scenarios. From the multi-risk assessment point of view, the cascading scenarios of main interest are those that produce an amplification of the effects arising from these hazards when compared to those produced by the single events. With a set of cascading scenarios of interest, their quantification can be achieved by adopting different strategies, for example, by analysing databases of past events, performing physical modelling for the propagation of the intensity measures of interest, and/or by performing expert elicitations in order to obtain information about extremely complex problems, or in these cases, with poor data or needing rapid analysis. Table 7 summarizes the main chains of cascade events identified for the different hazards affecting the Naples test city.

25

Table 7. Possible cascade scenarios for the Naples test case (from MATRIX D3.3).

Earthquakes

Triggered events Volcanic

Landslides

Extreme temperature

Heavy precipitation

Extreme wind

Wildfires

Tsunami

general) (in floods

Volcanic eruption

Landslides

Earthquakes

Meteorological events

1 1

1, 4

Volcanic eruption

1d

Tephra fall

2

1

Pyroclastic flows

2

1

Lava flows

2

1

Lahars

2,3

2

Volcanic earthquake s

2

1

3

3

1?

1

1

Floods Tsunami Wildfires

1

Two different scenarios were considered in the Naples test case for more quantitative analysis:  Scenario 1: Volcanic earthquakes and seismic swarms triggered by volcanic activity (interaction at the hazard level).  Scenario 2: Simultaneous occurrence of ash-fall and earthquake hazards (interaction at the vulnerability level). The following results are presented in greater detail in MATRIX deliverable D3.4 (GarciaAristizabal et al., 2013a).

26

3.2 Scenarios of cascading effects analysed in the Naples test case This section summarizes the quantitative analysis of interaction scenarios performed in the Naples test case and that have been presented in the MATRIX deliverable D3.4 (GarciaAristizabal et al., 2013a). 3.2.1 Volcanic earthquakes and seismic swarms triggered by volcanic unrest (scenario 1) Naples is located in the midst of two volcanic areas, the Phlegrean Fields (in the western side) and Mount Vesuvius (in the E-NE side). It is a case of a densely populated area located near active volcanoes, particularly in the Campi Flegrei area. The Campi Flegrei is dominated by a resurgent nested caldera resulting from two main collapse events related to the formation of the 37 ka Campanian Ignimbrite and the 12 ka Neapolitan Yellow Tuff (e.g., Rosi et al., 1983; Orsi et al., 1992). The Campi Flegrei caldera (CFc) comprises both submerged and continental parts on the western portion of the Bay of Naples, and appears to have undergone resurgence since 10-15 ky (Orsi et al., 1996). The latest eruption occurred in 1538 during which the Monte Nuovo scoria cone was generated (Di Vito et al., 1987). During the last four decades it has experienced a huge uplift phase, which reached about 3.5 m in 1985, after which a subsidence phase started (e.g., Troise et al., 2008). Unrest episodes (e.g., from 1969 to 1972 and from mid-1982 to December 1984, Barberi et al., 1984) have been observed, during which swarms of volcanic seismic events were recorded (Dvorak and Gasparini, 1991; De Natale et al., 1995; Troise et al., 1997). Scenario description Eruption forecasting refers, in general, to the assessment of the occurrence probability of a given eruptive event, whereas volcanic hazards are normally associated with the generation of superficial and evident phenomena that usually accompany eruptions (e.g., lava, pyroclastic flows, tephra fall, lahars, etc.). Nevertheless, during the evolution of a volcanic system from a quiescent state to an eruptive state, a large number of small- to moderate-sized earthquakes occur. In fact, periods of volcanic unrest are normally accompanied by (volcanic) seismic activity located directly below the volcano (Endo et al., 1981; Sawada and Aramaki, 1989; Matsumura et al., 1991; Oshima et al., 1991; Power et al., 1994; Garcia-Aristizabal et al., 2007), and/or by triggered seismic swarms (sometimes categorized as ‘distant volcanotectonic events’, e.g., White and Power, 2002) associated with fault systems nearby the volcano (e.g., Minakami, 1974; Shimizu et al., 1992; Wolfe, 1992; Legrand et al., 2002). Several of the world’s active volcanoes are located near densely populated areas (for example, in Italy (e.g., Naples), Japan (e.g., Tokyo), Mexico (e.g., Mexico D.F.), Colombia (e.g., Manizales), Ecuador (e.g., Quito)), and therefore the seismic hazard associated with pre-eruptive earthquake activity can be considered from a cascading effect perspective within a multi-risk framework. The problem to be analysed is to assess the possible effect that local seismicity (normally characterized by many shallow micro-earthquakes) triggered during volcanic unrest may have in the area where the volcano is located, and compare it with the general seismic hazard analysis. In fact, seismic swarms triggered by volcanic unrest activity are composed of a high number of shallow events (rates are higher than normal seismicity of 27

the area), with generally low magnitude. Nevertheless, from a seismic hazard point of view, they may be of interest when exposed elements (e.g., populated areas or facilities of interest) are located very close to the volcano, since they are generally very shallow (few kilometres depth). In the Neapolitan area, this case may be of particular interest given the highly urbanized area in the surroundings of Campi Flegrei. For example, two unrest episodes (called bradisisma, from the Greek word for “slow earthquake”) occurred from 1969 to 1972 and from mid-1982 to December 1984 (Barberi et al., 1984). The first episode produced a net uplift of about 1.7 m (Troise et al., 2008), with only weak seismic activity was recorded during this phase. The second event presented a similar uplift pattern, but compared to the previous one, this episode was accompanied by greater seismicity, both in terms of the number of earthquakes and their magnitude (Dvorak and Gasparini, 1991; De Natale et al., 1995; Troise et al., 1997). Shallow earthquake hypocentres (depths from 1 to 4 km) occurred from mid-1983 to December 1984 and were able to cause damage to buildings in the Pozzuoli area (e.g., Troise et al., 2008). A magnitude 4.2 earthquake with its epicentre in Pozzuoli occurred on October 4, 1983. In 1984, the town was almost completely evacuated and about 40,000 people were relocated. The total uplift recorded between January 1982 and December 1984 amounted to 1.8 m. Since the end of 1984, the ground has started to subside with a peak subsidence rate of about 7–8 cm/year in 1985 and at a progressively slower rate until 2001–2002, when the subsidence ended and about 0.9 m of the total uplift was recovered (De Natale et al., 2006). Superimposed on the subsidence phase were some small and fast uplift episodes which occurred in 1989, 1994, and 2000, i.e., with an apparent period of about 5–6 years (Gaeta et al., 2003; Lanari et al., 2004). These smaller uplift episodes reached a maximum uplift of about 1–8 cm, lasted 4–10 months, and were accompanied by seismic swarms of small earthquakes (ML