1. Introduction
Earthquakes often trigger the closure of the road network, causing obstacles in rescue and evacuation procedures. Efficient mitigation strategies and appropriate emergency plans should consider the probability of road closures in order to ensure rapid emergency management, and to minimize the impact of disasters in terms of casualties, the disruption of activities, and economic losses. Although there are growing studies on the accessibility modeling of road networks for the essential elements of emergency management in earthquake cases, such as operational centers, hospitals, fire stations and emergency areas, the literature that incorporates road closure probabilities into the accessibility analysis is quite limited. An interesting integration can be found in Ertugay et al. [
1].
In the aftermath of a disastrous seismic event, road blockages can be direct or indirect (
Figure 1). The former directly affects the road pavement due to ground failures from liquefaction, landslides, and fault ruptures, or through the structural damage of its additional components (bridges, tunnels, retaining walls, etc.), while the latter are due to the debris distribution generated by the collapse of the surrounding environment.
Generally, outside inhabited centers, this distributed debris is due to the collapse of elements in the natural environment (rockfalls and landslides), while inside inhabited centers, it is caused by the collapse of buildings adjacent to road edges (interfering buildings). In urban areas, road blockages due to collapsed buildings are the most important mobility limitation for the road network [
2,
3]. Therefore, the vulnerability of interfering buildings plays a key role in assessing the seismic risk of an urban road network.
In Argyroud et al. [
4], a summary is provided of the different methodologies to develop fragility curves for transportation infrastructures subject to seismic activity; however, the authors focused on the geotechnical effects, and did not consider road blockages from damage to interfering buildings.
Analyzing the collapsed buildings triggered by the earthquakes that occurred between 1988 and 2003 in Algeria, Armenia, India, Japan, Mexico, Romania, Taiwan and Turkey, researchers [
5] considered ten different damage types categorized into the following five groups (
Figure 2): (i) inclined layers; (ii) pancake collapses; (iii) debris heaps; (iv) overturn collapses; and (v) overhanging elements.
In the group “inclined layers”, the structure supporting a floor slab or a flat roof collapses on one side, generating an inclined plane that can affect one floor, several floors, or the whole building. “Pancake collapses” are due to the uniform collapse of the structural elements of a floor or of all or some floors (so-called soft storeys). These collapses are typical of reinforced concrete (RC) buildings. They generate modest debris areas on the road pavement. In the group “debris heaps”, all or some of the structural elements are completely destroyed, generating significant debris areas. Commonly, such collapses are typical of masonry buildings. “Overturn collapses” are characterized by two types of damage, namely “overturn collapse–separated” and “overturn collapse”. In the first damage type, the building footprint area is still recognizable, but a new structure (the upper part of the building) can be found nearby. In the second damage type, the building maintains a single body, but the inclination of its vertical axis is changed until the building is completely overturned (in urban areas, complete overturning is often prevented by the surrounding buildings or debris). For this damage type, soil conditions play an important role. Finally, in the “overhanging elements” group, the load-bearing external walls are destroyed, but the floor slab or the roof remains in its initial position. Depending upon the building type (RC, masonry, etc.), the five groups highlight the need to consider certain appropriate predictive models capable of assessing the debris area of the collapsed buildings and their possible distributions around the damaged building.
This work aims to evaluate the probability of the opening or closure of the urban road network in the emergency management phase through a methodology that allows fragility curves to be generated based on the debris distribution of interfering buildings. To the authors’ knowledge, there are no fragility curves that interpret this phenomenon. In particular, these curves should express, for a given seismic intensity measure (IM) level, the probability that the width of the road pavement is at least 3.50 m, in order to allow the passage of emergency vehicles after a debris fall.
The proposed methodology is specifically tailored for a minimum amount of information provided by the analysis of the limiting conditions of the emergency (LCE). This analysis was formulated in 2012 by the Italian Civil Protection Department [
6,
7]. Based on a limited amount of information, the analysis identifies those physical elements (strategic buildings = ESs, emergency areas = AEs, road network segments = RAs, structural units = USs, structural aggregates = ASs) strictly necessary for the emergency management of an urban settlement. Together with LCE, these minimal physical elements are assimilated within a network of arcs and nodes, where arcs are the road network segments.
2. Methodology
The proposed methodology is valid for RC and masonry buildings, and has been validated for post-earthquake observations and experimental literature data. It can be developed into four key steps:
- Step 1.
Definition and characterization of the road type;
- Step 2.
Assignment of vulnerability functions;
- Step 3.
Construction of debris graphs;
- Step 4.
Construction of a fragility curve for each road segment.
- Step 1.
Definition and characterization of the road type.
In order to define the type of road related to the road network segments of the urban settlement, it is necessary to identify the following:
(i) The geometric characteristics of the analyzed road segment (the average width of the road pavement (wr), the average width of the sidewalks (wbr), and the average distance between the opposite interfering buildings (wbb) along the road sides);
(ii) the types of interfering buildings constituting the structural aggregates or isolated structural units (structural material, number of floors, age of the building);
(iii) the geometric characteristics of the structural aggregates or of the interfering isolated structural unit (average height of the roof covering (hb) and the footprint dimensions in transversal (a) and longitudinal (b) directions).
This information can be easily found in the LCE analysis of RAs, USs and Ass, as well as through satellite images.
- Step 2.
Assignment of vulnerability functions.
Once the types of interfering buildings constituting the structural aggregates or the isolated structural units are identified, two or more vulnerability functions should be assigned for each side of the road. In this way, the typological differences of the structural units constituting the aggregates and all those characteristics (structural and architectural systems, building use, and structural regularity/irregularities), not easily identifiable during the field survey necessary for the LCE analysis, can be taken into account.
Vulnerability functions express the relationship between the selected
IM (e.g., peak ground acceleration, PGA) and the damage factor (
DF), where the letter represents the expected average ratio of the repair cost to the replacement cost of the building. Such vulnerability functions can be defined using the methodologies and results of scientific literature (interesting approaches for structural aggregates and/or structural units can be found in [
8,
9,
10]), or utilizing datasets and catalogs that are available on web platforms as resources for open risk modeling. Alternatively, they can be generated by sets of fragility curves using, for example, the following expression:
where
n is the number of damage states considered;
is the mean damage factor value for a given damage state
dsi (using, for example, the collected values from the vulnerability literature [
11]);
is the probability of a building of sustaining the damage state
dsi for a specific level of
IM. Another interesting approach to define vulnerability functions from fragility curves can be found in Vona et al. [
12]. Obviously, the vulnerability functions or the fragility curves used should represent the structural types and the real conditions of the analyzed buildings.
- Step 3.
Construction of debris graphs.
For each assigned vulnerability function, it is necessary to construct different graphs of debris width (wd) extending further than the initial width of the structure (a), in order to consider the possible debris distributions around the damaged building. These graphs express the relationship between the selected IM and the debris heap width resulting from building damage.
As shown in
Figure 3, using the assigned vulnerability functions, it is possible to define the
IM value in which a certain performance level is reached, thanks to the
DF values provided by the Italian guidelines for the classification of the seismic risk of buildings [
13,
14], the so-called SISMABONUS. These
DF values are shown in
Table 1. They are valid for typical RC and masonry buildings with residential use, and can conventionally be used for buildings with different uses or structural materials.
According to the current Italian seismic code [
15] (NTC 2018 hereafter), in the event of an earthquake, any building meeting the operative limit state (OLS) is considered to have no damage to its structural elements and only minor damage to its non-structural elements. Therefore, for this damage state, the building should not suffer significant damage, so as to produce appreciable debris on the road. Consequently, consistent with
Table 1, for each
IM in which
DF ≤ 7%, the structure retains its footprint area (
Af), and the debris width can be assumed equal to zero (
wd = 0). For damage beyond the OLS, the structure instead has a different behavior.
For masonry buildings (
Figure 3a), it is assumed that the maximum debris area (
Ad,max) is reached in the collapse limit state (CLS). Such debris area encompasses the building footprint area and the relative debris heap on the road. Meeting this limit state,
Ad,max remains constant, while the angle of inclination of the debris heap (therefore the debris volume
Vb) increases up to the maximum
DF value (100%). In the transition from OLS to CLS, it is assumed that the
Ad grows linearly with the increasing
IM. In order to evaluate the maximum debris area in the SLC, it is possible to use the methodology developed by Domaneschi et al. [
16]. Based on the geometric and typological characteristics of the masonry buildings of central Italy (building size, wall thickness, arrangement of openings, indoor spaces distribution, etc.), in [
16], the debris area of collapsed buildings was determined by using an amplification factor (
ε) of the building footprint area, which depended on the geometric characteristics of the structure. In particular,
Ad,max and
ε can be calculated as follows:
where
a,
b,
hb,
Af =
a·
b, and
Vb =
a·
b·
hb can be assumed to be the geometric characteristics of the structural aggregates or of the interfering isolated structural unit (see step 1). The debris area can have different distributions around the damaged building. To take into account the uncertainty of the possible distributions, it is possible to use the simplified collapse models proposed by some authors [
17,
18]. These models correspond to a collapsed building in one or two directions. Consequently, an equation that estimates the debris width (
wd) can be defined for each debris distribution. For the models corresponding to a collapsed building in one direction, the corresponding
wd can be evaluated as follows:
where
wd,2 describes the case in which the debris distribution, which does not affect the analyzed road network segment. For the collapsed building in two directions, the relevant
wd can be calculated as follows:
Based on observations from past earthquakes, the main debris distribution of RC buildings is due to the collapse of nonstructural elements. For these buildings (
Figure 3b), it is therefore assumed that the maximum debris area is reached in the life safety limit state (SLS), beyond which the debris area remains constant, while debris volume increases up to the maximum
DF value. This assumption is justified by the fact that, according to NTC 2018, in meeting the SLS, the failure of nonstructural components is expected. In the transition from OLS to SLS, it is assumed that the debris area grows linearly with increasing
IM. The infill panels may not collapse on the analyzed road network segment, or they may collapse with complete or partial overturning. To take into consideration the randomness of these collapse mechanisms, in this case it is also possible to assign different graphs of the debris widths. In order to calculate the debris widths expected in SLS, the following expressions are suggested:
where Δ
h is the average height of the infill panels or approximately the average floor height of the structure. Equations (7) and (9) correspond to collapse mechanisms with complete or partial overturning, respectively, while Equation (8) considers the case in which the overturning of the infill panels does not affect the analyzed road network segment. For example, this case can also be found in pancake collapses in which there is a minimal debris area.
For each assigned vulnerability function, three graphs of debris width should be constructed in order to take into account the possible debris distributions around the damaged building. For each debris distribution, the graph is built based on the debris widths and
IM values in which the considered building performance levels are reached. These performance levels can be identified with the
DF threshold values shown in
Table 1 on the
DF values axis of the vulnerability function. According to the assumptions of the proposed methodology, for each
IM in which
DF ≤ 7%, the debris width can be assumed equal to zero (
wd = 0), while for each
IM in which
DF ≥ 80% (for a masonry building) or
DF ≥ 50% (for an RC building), the debris width can be defined using Equations (4)–(6) for the masonry building, or Equations (7)–(9) for the RC building. Based on these boundary conditions, a reasonable assumption can be made that for the other
IM values, the debris width varies with a linear law, and the generic graph of debris width can be defined.
- Step 4.
Construction of the fragility curve for each road segment.
Considering the geometric characteristics of the analyzed road network segment, and comparing the debris graphs related to the interfering buildings constituting the considered structural aggregate along the road edges, for each
IM it is possible to determine different width values of the road pavement that remains free after the debris fall (
wfr), using the following expression (
Figure 4):
where
wbb is the average distance between the opposite interfering buildings along the road sides,
wbr1 and
wbr2 are the average widths of the sidewalk on side 1 and 2 of the road, respectively,
wd1,i is the debris width of the
i-th debris graph related to the structural aggregates on side 1, and
wd2,j is the debris width of the
j-th debris graph pertaining to the structural aggregates on side 2.
As shown in
Figure 5, by applying Equation (10) for each
IM considered, it is possible to identify the total number of analysis points and the number of points where the width of the road pavement is reduced beyond the 3.50 m (
wfr ≤ 3.50 m), required for the passage of emergency vehicles (number of closures or collapses). In particular, for these closures or collapse points, the required performance level can be defined through the width of the debris heap on the road pavement (
wdr), which can assume the following values:
where
wr is the average width of the road pavement expressed in meters. The proposed methodology, therefore, concludes by applying a statistically appropriate method for fitting the fragility curve to all these points.
Fragility curves are usually described by a lognormal probability distribution function according to the following equation:
where
is the probability of exceeding a particular performance level,
PL, for a given
IM, Φ is the standard cumulative probability function,
β is the standard deviation of the natural logarithm of
IM for the performance level, and
θ is the median threshold value of
IM required to cause the performance level. According to Equation (12), the development of fragility curves requires the estimation of two parameters, namely
θ and
β.
At each intensity level
IM =
xj, the proposed procedure provides a certain number of closures or collapses out of a total number of analyses. In order to estimate the fragility curve parameters, the appropriate fitting technique for this type of data is to use the method of maximum likelihood estimation [
19]. More specifically, the estimate of the fragility function parameters is obtained by maximizing the logarithm of the likelihood function as follows:
where
m is the number of
IM levels and
zj is the number of collapses out of
nj ground motions with
IM =
xj. To solve Equation (13), Baker [
19] implemented simple software tools for the calculations. The user only needs to provide observed data on the
IM levels considered, and the relevant numbers of analyses and collapses, and the software tools estimate the fragility curve parameters. An application example of this analysis is provided by the same author in order to obtain the fragility curves of structures with a minimal number of structural analyses. Obviously,
IM levels should focus the fitting on the region of the fragility curve most important for risk assessments, rather than considering analysis at high or low
IM levels. They can be defined according to the seismic hazard of the study area.
Figure 6a shows the fragility curve construction for a generic road network segment using the aforementioned software tools and the plotted data in
Figure 5.
Based on the fragility curve in terms of the probability of closure
, for each
IM, the probability that the width of the road pavement is at least 3.50 m to allow the passage of emergency vehicles after debris fall
, can be evaluated as follows:
Figure 6b shows the relevant fragility curve. Consistent with LCE analysis, this type of fragility curve can be used directly to assess the operational efficiency of the road network segments (RAs), providing a quantitative measure of the probability to preserve their efficiency after a seismic event with a given
IM level.
3. Case Study
The proposed methodology was applied to a case study of the “Corso Umberto I” in Amatrice, Italy. Amatrice extends along the East–West direction with parallel roads and is crossed by Corso Umberto I, the main street of the city. Corso Umberto I is located in the historic center of the city. Although it plays a key role in the LCE analysis of the municipality of Amatrice, the earthquakes that hit Central Italy on 24 August, 26 October, and 30 October 2016 triggered major damage, mostly to buildings and architectural heritage structures along the road, causing closure of the road already with the seismic events of 24 August. More specifically, two shocks were felt on 24 August: the main shock at 3:36 a.m. (local time) with magnitude Mw = 6.2 and epicenter close to Accumoli (Rieti province), and the aftershock at 4.33 a.m. with a magnitude Mw = 5.5 and epicenter close to Norcia (Perugia province) [
20]. According to the Italian Civil Protection Department [
21], these seismic events caused a total of 299 fatalities, and 386 people were injured and about 4800 rendered homeless.
Based on LCE analysis, Corso Umberto I consisted of two road network segments (AC1 and AC2). In addition, there were eight interfering structural aggregates and one isolated structural unit along the road. Using the LCE information available on the web-GIS platform [
22], all data were found for each road segment and structural aggregate in order to identify their geometric and typological characteristics. More specifically, AC1 and AC2 had an average width of road pavement of 5.50 m (
wr = 5.50 m), an average width of sidewalks of 1.00 m on each side (
wbr1 =
wbr2 = 1.00 m), and an average distance between the opposite interfering buildings along the road sides of 7.50 m (
wbb = 7.50 m). In terms of structural aggregates, there were 100 masonry buildings and only three RC buildings (the debris of masonry buildings was clearly prevalent). The buildings had a covered area of approximately 15,600 square meters. They varied from one to four storeys. Depending on the geometric and typological characteristics of the aggregates, the critical areas of AC1 and AC2 are circled in red in
Figure 7 (road segment between aggregate 1 and 2 for AC1 and road segment between aggregate 4 and 6 for AC2). According to Equations (2)–(4), in these areas the structural aggregates reached their maximum debris widths along the infrastructures considered. These critical areas were the road portions with higher risk of loss of connectivity. The maximum area and width of the debris (
wd,max) for each structural aggregate in the critical areas are shown in
Table 2.
In order to describe the different seismic behaviors of the buildings constituting the structural aggregates present in the critical area of AC1 and AC2, six types of fragility curves were selected, and consequently, using Equation (1), six vulnerability functions were generated and assigned for both sides of the two road segments (
Figure 8). In particular, using the OpenQuake platform [
23], the fragility curves developed by Ahmad et al. [
24] were selected, which were analytical fragility functions for masonry buildings and building aggregates expressed in terms of PGA. According to these authors, the prototypical buildings represent Euro–Mediterranean buildings in general, and Italian and Slovenian buildings in particular. The fragility curves used were valid for unreinforced masonry buildings with low and high void percentages and for stone buildings of two and four storeys, respectively.
For each assigned vulnerability curve, based on Equations (4)–(6), three graphs of debris width were constructed in order to take into account the different debris distributions along the critical areas of AC1 and AC2. These graphs illustrate the debris distribution generated by the interfering buildings constituting the structural aggregates along the road edges of the critical areas. In fact, the debris distribution can affect the following:
- (a)
Corso Umberto I (wd = wd,max on Corso Umberto I);
- (b)
the roadway opposite Corso Umberto I (wd = 0 on Corso Umberto I);
- (c)
both roads (wd = wd,max/2 on Corso Umberto I).
As indicated in step 3 of
Section 2, the graphs of debris width depend on the geometric and typological characteristics of buildings, the collapse model considered, and the
IM values for which the OLS (
DF = 7%) and SLC (
DF = 80%) are reached.
Taking into account the geometric characteristics of AC1 and AC2 and comparing the debris graphs (18 diagrams for each side of the road), using Equation (10), 324
wfr values (number of analyses) were calculated for each
IM. Considering the number of points for which
wfr ≤ 3.50 m (or
wdr ≥ 2.00 m), the number of closures or collapses was identified for each
IM. Based on the number of analyses and closures or collapses identified for 0.1 g ≤
IM ≤ 0.5 g (region of the fragility curve most important for risk assessments), using a simple software tool [
19] to solve Equation (13), the fragility curve parameters were estimated for AC1 and AC2.
Figure 9 shows the fragility curves of the analyzed infrastructure with the relevant construction points. Finally, using Equation (14), the probability that the road was open for the passage of emergency vehicles after debris fall could be defined for each
IM (
Figure 10).