Abstract
Recent earthquakes have exposed the vulnerability of existing buildings; this is demonstrated by damage incurred after moderate-to-high magnitude earthquakes. This stresses the need to exploit available data from different sources to develop reliable seismic risk components. As far as it regards empirical fragility assessment, accurate estimation of ground-shaking at the location of buildings of interest is as crucial as the accurate evaluation of observed damage for these buildings. This implies that explicit consideration of the uncertainties in the prediction of ground shaking leads to more robust empirical fragility curves. In such context, the simulation-based methods can be employed to provide fragility estimates that integrate over the space of plausible ground-shaking fields. These ground-shaking fields are generated according to the joint probability distribution of ground-shaking at the location of the buildings of interest considering the spatial correlation structure in the ground motion prediction residuals and updated based on the registered ground shaking data and observed damage. As an alternative to the embedded coefficients in the ground motion prediction equations accounting for subsoil categories, stratigraphic coefficients can be applied directly to the ground motion fields at the engineering bedrock level. Empirical fragility curves obtained using the observed damage in the aftermath of Amatrice Earthquake for residential masonry buildings show that explicit consideration of the uncertainty in the prediction of ground-shaking significantly affects the results.
Similar content being viewed by others
1 Introduction
Accurate assessment of seismic risk for buildings at a territorial scale depends to a large extent on the availability of reliable and accurate fragility curves. The process of fragility estimation needs to make the most out of background information, expertise, available data and modelling/analysis capacities (Shinozuka et al. 2000; Jaiswal et al. 2011). The seismic fragility curves can be distinguished, based on the type of data used for deriving them, into four categories (Rossetto and Elnashai 2003); namely, analytic (e.g., Porter et al. 2007; Baker 2015; Jalayer et al. 2017), empirical (e.g., Shinozuka et al. 2000; Jaiswal et al. 2011), based on expert opinion (e.g., Mosleh and Apostolakis 1986; Porter et al. 2007), and hybrid (e.g., Singhal and Kiremidjian 1998). Several advantages and shortcomings can be associated to the empirical seismic fragility assessment. The main advantage is that, being based on observed damage, the empirical fragility curves are—potentially—able to provide a realistic picture of post-earthquake damage. For instance, they can consider the possible soil-structure interactions, the cumulative damage due to aftershocks (or back-to-back seismic events), the degradation in structural behaviour due to aging, and the effect of flexible diaphragm and secondary members (e.g., infills); just to name a few. The shortcomings are mainly associated to the difficulties in creating a homogenous (e.g., in terms of building class characterization, soil conditions, …, etc) spatial database of observed damage, inaccurate estimation of seismic ground shaking, site effects and spatial correlations in recorded ground motions and observed damage (e.g., Crowley et al. 2008; Colombi et al. 2008; Rossetto et al. 2014; Ioannou et al. 2015; Lallemant et al. 2015), the errors associated to the assignments of the observed damage state (e.g., Colombi et al. 2008; Rossetto et al. 2014) and finally a potential bias that could be introduced by underestimating the number of undamaged buildings (post-earthquake surveys tend to identify the damaged buildings) and hence the total number of buildings. In fact, the use of an established and standardized damage scale, in which different damage data should be converted, is one essential requirement for the derivation of empirical fragility curves. This permits to update the existing fragility curves as new data is gathered and to compare the vulnerability of different building classes. The European Macroseismic Scale (EMS-98, Grünthal 1998) is an example of a standard damage scale, very often used to report the empirical fragility curves in Europe.
1.1 A brief history of empirical fragility assessment for buildings in Italy
Several works have focused on empirical assessment of the vulnerability of Italian building stock (mainly masonry and reinforced concrete). Braga et al. (1982) used the MSK-76 damage scale (Medvedev–Sponheuer–Karník) to derive damage probability matrices (DPM) showing the probability of being in each of the five damage levels of MSK-76 for three building vulnerability classes and a 12-level damage-dependent seismic intensity scale (MSK-76). Sabetta et al. (1998) derived both the vulnerability (i.e., mean normalized damage versus seismic intensity relations, see Calvi et al. 2006) and the fragility curves by adopting PGA and Mercalli Cancani Sieberg (MCS) macroseismic scale as measures of seismic intensity. This study adopted the five damage levels of MKS-76 scale as measure of damage for three building vulnerability classes. Orsini (1999) has derived continuous vulnerability curves expressed as cumulative distribution of the seismic intensity level (a parameter-less scale of intensity, PSI, Spence et al. 1992) related to the excursion of a given damage level. Dolce et al. (2003) have derived damage probability matrices and scenario-based mean damage maps based on MSK-76 scale considering the site effects for four building vulnerability classes for town of Potenza in Southern Italy. Di Pasquale et al. (2005) have reported the damage probability matrices for four vulnerability classes as a function of MCS intensity. Lagomarsino and Giovinazzi (2006) have derived vulnerability curves reporting the mean damage ratio versus seismic intensity expressed in the EMS-98 scale for masonry and RC buildings. Zuccaro and Cacace (2009) mapped the Italian building inventory in terms of damage probability matrices based on EMS-98 scale. Zuccaro and Cacace (2015), instead, proposed a methodology to reduce the variability in the vulnerability classification of EMS-98 through the application of vulnerability modifiers. These modifiers expressed the shift in vulnerability (expressed in terms of a synthetic vulnerability index) for various “sub-classes” identified within a given vulnerability class with respect to the overall vulnerability index for that class.
Rota et al. (2008) derived empirical fragility curves based on damage survey data for 150,000 buildings from past Italian earthquakes, occurred in the period spanning from Irpinia (1980) to Molise (2002) earthquakes. These empirical fragility curves are reported using EMS-98 as the damage scale and PGA as the ground shaking intensity. The seismic intensity is calculated as a single mean estimate evaluated for the corresponding municipality based on Sabetta and Pugliese (1996) ground motion prediction equation (GMPE) and assuming rock as the soil type. Liel and Lynch (2012) employed post-earthquake investigations following the L’Aquila earthquake to quantify the seismic fragility of RC buildings. The seismic intensity in this study is characterized by PGA derived from the ShakeMap (Michelini et al. 2008) for L’Aquila’s main event. De Luca et al. (2015) and Del Gaudio et al. (2016) derived empirical fragility curves based on a database of post-earthquake building surveys conducted on 131 RC buildings located in the town of Pettino after the 2009 L’Aquila Earthquake. The EMS-98 compatible damage levels were deducted from the standard Italian survey sheet AeDES (Baggio et al. 2007) and seismic intensity was expressed in terms of PGA contours of the L’Aquila main event ShakeMap. Del Gaudio et al. (2017) have derived empirical fragility curves based on the large database of post L’Aquila 2009 Earthquake damage surveys conducted by the Italian Civil Protection (Dolce and Goretti 2015) on RC buildings. De Luca et al. (2018) have used the observed damage for infilled RC structures in the aftermath of the 24th August 2016 Amatrice Earthquake for Bayesian updating of existing empirical fragility curves for central Italy. Rosti et al. (2018) have used the database of damage data collected after the 2009 L’Aquila (Italy) earthquake, to derive damage probability matrices for several building types representative of the Italian building stock. Ioannou et al. (2020) have employed a Bayesian framework to study the importance of considering the uncertainty in the ground motion intensity on the shape of empirical fragility curves obtained based on post-disaster data aggregated at the municipality level for 1980 Irpinia Earthquake.
1.2 The focus: empirical fragility estimation based on conditional GMPE-based ground-shaking fields
One important aspect that emerges from the available rich literature on empirical vulnerability and fragility assessment is that accurate estimation of the ground shaking intensity is crucial (as much as accurate damage estimation) towards accurate and reliable empirical fragility assessment (e.g., Hsieh et al. 2013; Yazgan 2015; Lallemant et al. 2015; Ioannou et al 2015). Gaussian random ground shaking fields, whose first- and second-moment statistics are described by the ground motion prediction equations (GMPE) and the compatible spatial correlation models (i.e., models describing the intra-event variability in the residuals of the same GMPE), are used quite often in the literature (e.g., Park et al. 2007; Goda and Hong 2008; Jayaram and Baker 2009; Goda and Atkinson 2010; Goda 2011; Sokolov and Wenzel 2011, 2013, Weatherill et al. 2015). Based on a procedure used originally for generation of ground shaking fields in the frequency domain (Vanmarcke et al. 1993), Park et al. (2007) and also Crowley et al (2008) proposed an updating procedure that yields updated or “conditional” Gaussian random fields based on available recordings for various stations (a.k.a., correlated ground motion fields with “anchor points”). The conditional Gaussian random ground shaking field has been used mostly for portfolio loss assessment (e.g., Park et al. 2007; Crowley et al. 2008; Miano et al. 2015, 2016); however, it has not been used extensively for empirical fragility assessment. Straub and Der Kiureghian (2008) discuss the differences between empirical fragility assessment based on correlated observations (such as ground shaking for a given event) and that based on independent observations. Hsieh et al. 2013 use the 444 free-field ground motion recordings of the Chi-Chi Taiwan Earthquake to generate the map of ground shaking. They use a kriging model to interpolate the ground-shaking for the points with no recording stations nearby. Lallemant et al. (2015) recommend updating the ground-shaking field based on the recorded ground shaking data. They perform the updating by selecting the ground shaking field realizations that match the recorded ground shaking data (two stations in the example). Ioannou et al. (2015) perform updating using a kriging model in the case where the ground motion recordings are available at the location of a subset of the surveyed buildings. Yazgan (2015) uses the conditional ground shaking fields in a Bayesian framework to update the parameters of prescribed fragility models based on the likelihood of observing damage in the location of the surveyed buildings for a given class.
It is interesting to note that the procedure employed in this work for updating the statistical properties of the Gaussian random fields does not rely on the assumption that recording stations need to be close to the buildings’ location, as discussed later on. In fact, in the case of Amatrice Earthquake, most of the ground motion recording stations are located far field (this is shown later in the paper). Nevertheless, the updating procedure can include all the stations. The updating procedure can be carried out automatically and with no significant computational effort (the procedure can be implemented using Matlab’s toolboxes).
This work presents a procedure for the generation of conditional GMPE-based ground shaking fields for derivation of empirical fragility curves. Section 2 presents the methodology employed in this work for deriving empirical fragility curves. Section 3 presents the application of the methodology in deriving empirical fragility curves based on damage data for masonry buildings available right after the Amatrice Earthquake of 24th of August 2016. Section 4 is dedicated to the final conclusions drawn from the results outlined in Sect. 3.
2 Methodology
2.1 Basic assumptions for empirical fragility assessment for a class of buildings
The concept of fragility curve is applicable to Poisson-type limit state excursion (see e.g., Der Kiureghian 2005; Jalayer and Ebrahimian 2020). In the context of seismic risk, the uncertainty related to earthquake occurrence can be classified as leading to Poisson-type limit state excursion. Using the damage survey data from spatially-distributed surveyed buildings belonging to the same class for fragility modelling implies that the building-to-building variability within the same class is inevitably going to be considered. In other words, part of the dispersion captured by the resulting fragility curve is due to the differences between the surveyed buildings, within the same class, and not just due to the uncertainty in the ground motion representation. It can be argued that using a single class fragility curve is equivalent to the assumption that each of the surveyed buildings is replaced by a nominal average building for that class. One consequence of such assumption is that the building-to-building variability is going to be considered in “disguise” and as a contribution to record-to-record variability related to the ground motion. On the other hand, the ground motion intensity level at each building location (estimated in this work based on the simulation procedure discussed in the following paragraphs) cannot be independent as they correspond to the same earthquake event and they can exhibit various levels of correlation based on the mutual distance between the buildings. In other words, the ground-shaking levels registered at the location of each surveyed building can show significant correlation for close-by buildings. Excluding such correlation may lead to unrealistic and un-conservative fragility curves (Straub and Der Kiureghian 2008). Unconservative because the fragility curve that is derived based on the assumption of independent intensity values is going to rely on a larger information content compared to the fragility curve derived with explicit consideration of possible correlations between the intensity values—of course if such correlation is significant (this is the case when the buildings are close by). The spatial correlation in the ground motion is modelled herein by considering the spatial correlation structure in the residuals of the ground motion prediction equation (GMPE). Although assigning a damage grade to buildings based on rapid visual survey techniques (e.g., visual inspections, use of satellite imagery) is certainly subjected to errors (e.g., human/device error of judgment and/or an inevitable degree of subjectivity), the uncertainty in assigning the damage grades to surveyed buildings and the fragility model parameter uncertainties are not considered here.
2.2 Empirical fragility assessment for a class of buildings
The empirical fragility is calculated within an updated robust reliability assessment framework (Papadimitriou et al. 2001; Beck and Au 2002; Jalayer et al. 2010). This framework is also known in literature as “predictive fragility” (e.g., Sasani and Kiureghian 2001; Straub and Der Kiureghian 2008). The concept of robust reliability is used for calculating “robust” empirical fragility curves. The term robust in this context implies that the uncertainty in the evaluation of the ground shaking is considered (see Jalayer et al. 2017; Miano et al. 2018; Jalayer and Ebrahimian 2020 for examples of robust analytical fragility assessment). Let the survey dataset DDamage define the set of N damage states Di (e.g., the EMS-98 damage grade, i can vary between 0 and 5) for the buildings surveyed. Let PGA = {PGAi, i= 1:N} denote the vector of ground shaking values (expressed in terms of peak ground acceleration) estimated at the position of each surveyed building. The robust empirical fragility is calculated by considering the uncertainty in generating a plausible field of PGA values conditioned on both the vector of PGA values registered at various accelerometric stations, denoted by DPGA, and the observed damage pattern DDamage for the portfolio of the buildings surveyed. The updated empirical fragility conditioned on available data DDamage and DPGA can be written as follows:
where Di denotes a given damage grade; pga = x is a given level of ground-shaking; DDamage (the damage survey data) is the survey data acquired for the portfolio of buildings (for all the building classes); f(PGA|DPGA, DDamage) denotes the joint probability density function (PDF) for the ground-shaking field; Ω(PGA) denotes the domain of all plausible PGA fields; the joint PDF represents the joint probability of observing the ground shaking levels for building survey locations f(PGA|DPGA, DDamage) given the damage survey data DDamage and the registered PGA data DPGA, where PGA = [PGAi, i= 1:N] (see Sect. 2.3 for the derivation of f(PGA|DPGA, DDamage)). In general, it is quite complicated to calculate the integral in Eq. (1) analytically. It is usually estimated using standard Monte Carlo Simulation (MCS). In the adopted MCS approach, Nsim realizations of the ground shaking field PGA are generated based on the joint PDF f(PGA|DPGA, DDamage). The simulation procedure yields estimators of the updated empirical fragility and its standard deviation (see also Jalayer et al. 2016, 2017):
where P(D > Di|x,DPGA,DDamage) is the updated empirical fragility curve (for a given building class) and σ[P(D > Di|x,DPGA,DDamage)] is the standard deviation (error in the estimation) of the updated empirical fragility curve based on empirical damage data DDamage and registered PGA data DPGA.
2.3 Generation of conditional ground shaking fields
2.3.1 Generation of GMPE-based ground shaking fields
The joint probability density function f(PGA) for the vector of PGA = [PGAi, i= 1:N] values at the location of each building of interest for a given earthquake scenario can be evaluated by employing a ground motion prediction equation (GMPE). A full probabilistic representation (based on the first two moments) of GMPE can be expressed in terms of multi-variate Normal distribution which is identified by its expected value vector M and covariance matrix Σ. That is, once the first two moments are given, several realizations of the ground shaking field can be generated.
Herein, the ground motion prediction model proposed by Bindi et al. (2011, also known as ITA10) for the peak ground acceleration (PGA, the geometric mean of two horizontal components) as the intensity measure is used. A recent work by Michelini et al. (2019) ranks ITA10 quite favourably as the GMPE to use for Italian earthquakes from shallow active crustal zones. It is worth noting that there are more recent pan European BND14 (Bindi et al. 2014a, b) and Italian GMPE’s (ITA18, Lanzano et al. 2019) available. The latter includes the recent destructive Italian earthquakes such as Emilia Sequence 2012 and the Central Italy Sequence 2016. To authors’ knowledge, no well-established spatial correlation models have been calibrated so far to the residuals of these more recent GMPE’s. For instance, a very recent Italian correlation law (Sgobba et al. 2019) seems to have been calibrated to the residuals of a recent regional GMPE for Northern Italy (Lanzano et al. 2017, Emilia 2012 recorded ground motions) and not to ITA18. Nevertheless, the methodological procedure presented in this work is quite versatile and can use more recent GMPE’s and their related correlation structures as they become available. However, it is not imperative for this procedure that the adopted GMPE has the registrations of the earthquake of interest (in this case Amatrice Earthquake) embedded. In fact, the updating procedure (a sort of statistical base-line correction) aims to adjust potential biases with respect to a GMPE that does not include the registrations of the earthquake of interest. The functional form of the adopted GMPE model is the following:
where E[log10PGA] is the expected value (first moment) for the (base 10) logarithm of peak ground acceleration (PGA, in cm/s2); e1 is a constant term, FD (Rjb,M), FM(M), Fs and Fsof represent the distance function, the magnitude scaling, the site amplification and the style of faulting correction, respectively. M is the moment magnitude, Rjb is the Joyner–Boore distance in km (or the epicentral distance when the fault geometry is unknown—generally when M<5.5). The proposed equation for the distance function FD (Rjb,M) is:
the magnitude function FM(M) is expressed as:
where Mref, Mh, and Rref are fixed as Rref= 1 km, Mref= 5, Mh= 6.75. The functional form FS in Eq. 3 represents the site amplification and it is given by FS= sjCj, for j= 1:5, where sj are site amplification coefficients provided by Bindi et al. (2011) and Cj are dummy binary variables corresponding to the five different Eurocode 8 Part 3 (2007, EC8) site classes. Finally, the functional form Fsof represents the faulting correction coefficient and it is given by Fsof = fjEj, for j= 1:4, where fj are style of faulting coefficients and Ej are dummy binary variables used for the different faulting styles (i.e., normal (N), reverse (R), strike slip (SS) and unknown (U)). The values E[log10PGAi] (i= 1:N) from Eq. 3 constitute the components of the mean vector M. The covariance matrix, Σ, is defined as the sum of two inter-event and intra-event components:
where σintra represents the intra-event variability and σinter represents the inter-event variability (both parameters are tabulated in Bindi et al. 2011); e is the all ones matrix and R is the matrix of correlation coefficients. R is composed of unit diagonal terms and off-diagonals equal to ρjk, j≠k (both varying from 1 to N; where N is the total number of buildings surveyed). The covariance matrix is obtained according to the following formulation of ρjk (Esposito and Iervolino 2012):
where hjk represents the distance between sites j and k and b(T) is a tabulated coefficient equal to 10.8 km for PGA. It is to note that the above expression for the correlation coefficient has been calibrated to the residuals of the Bindi et al. (2011) GMPE adopted herein and therefore is reasonably consistent for the purpose used in this study.
2.3.2 Considering stratigraphic and topographic factors
Two alternatives are considered for considering the stratigraphic site effects; namely, (a) the coefficients imbedded in the GMPE (here Bindi et al. 2011, ITA10); (b) application of stratigraphic amplification factors to ground shaking at bedrock. Landolfi et al. (2011) and later Tropeano et al. (2018) propose stratigraphic coefficients that consider non-linear soil column propagation effects. Landolfi et al. (2011) have merged empirical, semi-empirical and analytical datasets to compute non-linear relationships quantifying stratigraphic amplification for different classes of subsoil profiles. That is, site effects are evaluated through the stratigraphic amplification factor, which is directly multiplied by the reference (i.e., propagated to bedrock) peak ground acceleration from the GMPE by Bindi et al. (2011) to obtain the peak acceleration at surface (where the subscribes B, C, D represent the soil class):
It is also possible to apply topographic amplification/deamplification factors (ST) to the GMPE. ST depends on the shape of slopes, because irregular surface geometry affects the focusing, defocusing, diffraction and scattering of seismic waves. This can lead to a change in amplitude, frequency and duration of ground motion compared to flat ground conditions (Sanchez-Sesma 1990; Paolucci 2002). Following Garcìa-Rodriguèz et al. (2008), a geometrical parameter more suitable for small scale studies seemed to be the slope curvature, which can be obtained from the digital elevation model (DEM) of the area (at 20 × 20 m resolution for the case-study presented here). This index permits to mark the concave and the convex features of a landscape, with negative and positive values respectively, accounting for attenuation in valleys and the seismic waves focusing on ridges. The effectiveness of this parameter was also validated by the numerical study of Torgoev et al. (2013) and adopted in seismic slope stability analyses by Silvestri et al. (2016) and reported in Table 1.
It should be noted that the above-mentioned stratigraphic and topographic factors are going to be applied herein in a deterministic manner. It can be shown that this changes the mean vector M by a constant (i.e., the inner product of a vector of constant factors and the median in the arithmetic scale) and leaves the covariance matrix Σ unaltered.
2.3.3 Updating the generated ground shaking fields based on registered PGA data
One interesting feature of the method adopted for generating the ground shaking field realizations is that it can be updated based on the registered values. It was assumed in Sect. 2.3.1 that the PGA values at the location of each surveyed building are distributed as a joint multi-variate lognormal distribution. One of the specific characteristics of a joint Normal distribution for a vector of variables is that any given partition of the vector conditioned on the remaining components of the vector is still going to be a joint Normal distribution. With specific reference to the case of the vector of log10PGA values denoted as data DPGA, let the vector of mean values M and the covariance matrix Σ be partitioned as follows (Park et al. 2007):
where M1 is the mean (of the base 10 logarithm) vector of PGA = [PGAi, i= 1:Ncl] values according to the adopted GMPE; M2 is the mean vector of calculated log10PGA at the stations within the area of interest (according to the adopted GMPE); Σ11 is the covariance matrix for the calculated (from the GMPE) log10PGA for the surveyed buildings of class CL; Σ12= Σ21 is the cross-covariance matrix for the log10PGA values calculated (from the GMPE) at the location of the surveyed buildings and those calculated at the location of the stations; Σ22 is the covariance matrix for the log10PGA values calculated at the stations.
As described briefly above, the conditional distribution of the calculated log10PGA values given the registered log10PGA values at the stations is a joint Normal distribution with mean vector M1|2 and covariance matrix Σ11|22:
where DPGA is the vector of the registered log10PGA values for the stations. Note that M1|2 and Σ11|22 represent the first two moments of the updated joint lognormal pdf f(PGA|DPGA).
2.3.4 Updating of the generated ground shaking fields based on observed damage
Once the first two moments of the updated ground shaking field corresponding to the lognormal PDF f(PGA|DPGA) are obtained, various plausible realizations of the ground-shaking random field can be generated. As the formulation in Eq. 1 indicates, the ground-shaking field is also conditioned on the observed damage survey data for the portfolio of surveyed buildings as reflected by the joint PDF f(PGA|DDamage,DPGA). A rejection sampling approach is adopted to reflect the conditioning of the random field joint PDF on the damage data. That is, the random field realizations that lead to non-monotonically increasing fragility functions are discarded. This is achieved by checking whether the gradient of the fragility curves’ slope is negative or not for all the classes and for all the damage levels considered. As a result, those generated random fields (for a given portfolio of buildings) which do not lead to meaningful fragility curves for all the classes and all the considered damage levels are discarded. This essentially means that the random field realizations that are not compatible with the observed damage pattern are not plausible. It is to note that, this updating procedure assumes that the damage survey is carried out with no error.
2.4 The empirical fragility curves
This section describes how the fragility curve P(D > Di|x,PGAk,DDamage) can be obtained for a given realization of vector of ground shaking values at the location of surveyed buildings PGAk generated based on pdf f(PGA|DPGA, DDamage) and given the damage survey DDamage. The empirical fragility curves P(D > Di|x,PGAk,DDamage) (k= 1:Nsim) have been calculated according to a logistic regression probability model (Wiesberg 1985, for instances of application in fragility assessment see e.g. Charvet et al. 2014; Lallemant et al. 2015; Yazdi et al. 2016; Jalayer and Ebrahimian 2017; De Risi et al. 2017). This type of regression is suitable for cases where the dependent variable is binary (i.e., either 0 or 1). Thus, it is especially suitable for estimating the probability of exceeding a damage state Di. That is, for damage state Di, each surveyed building can either: a) exceed the designated damage state denoted as 1 or b) NOT exceed the designated damage state, denoted as 0. Denoting probability of exceeding damage state Di as a function of the ground-shaking level PGA= xj as πi(xj), the likelihood of having ri buildings that exceed damage state Di for the NCL buildings surveyed for class CL can be expressed as follows:
Ncl is the number of surveyed buildings for building class CL; πi(x) is expressed through the following relation:
It is to note that Dcl is the vector of observed damage grades for the buildings in class CL. Note that Dcl is a subset of damage data DDamage for the entire portfolio of buildings surveyed. Figure 1 shows an example of empirical fragility curve obtained by employing logistic regression. As it can be observed from Fig. 1, the fragility is expressed as the conditional probability of exceeding damage level D5 for the class of buildings CL= MBC4 (Masonry buildings with tie rods or ring beams with number of stories > 2, described later in Sect. 3.1.2). The points correspond to the NMBC4 buildings surveyed for the class MBC4. The buildings whose damage state is estimated to be greater than or equal to D5 are assigned the value of 1 and rest of the surveyed buildings are assigned the value of 0.
3 Case study: masonry buildings damaged by Amatrice Earthquake
3.1 The study area
Between August and October 2016, Central Italy was stricken by three damaging earthquakes. The first Mw 6.0 event occurred on August 24th at 01:36 UTC close to Accumoli village (referred to as Amatrice Earthquake); it was followed by a long seismic sequence, which 2 months later produced a Mw 5.9 aftershock on October 26th at 19:18 UTC at 3 km West of Visso and a Mw 6.5 event on October 30th at 06:40 UTC, 6 km North of Norcia (see Ebrahimian and Jalayer 2017 for more details about the Central Italy seismic sequence; see Sextos et al. (2018) for more details about observed damage).
This paper focuses on the 3 municipalities and 4 fractions, namely, Amatrice, Accumoli, Arquata del Tronto, Tina (Accumuli), Illica (Accumuli), and Pescara del Tronto (Arquata del Tronto),Trisungo (Arquata del Tronto), affected by the 24th August 2016 event as shown in Fig. 2. The figure shows a simplified geological map of the area of interest, overlaid with the PGA contour map provided by INGV (Italian Institute of Geophysics and Volcanology) official ShakeMaps and the fault projection.
According to the PGA distribution provided by INGV, ground acceleration values as high as 0.5 g affected several small towns in the vicinity of the epicentre. Among these towns, Amatrice is the one that has had the most widespread destruction and the highest number of fatalities. The two hamlets of Tino and Illica have been considered as examples of the stratigraphic and the topographic amplification of ground motion. The masonry buildings located in the above-mentioned towns were strongly damaged by the August 24th event and its aftershocks.
3.1.1 Damage survey data DDamage: satellite imagery (Copernicus) and visual survey
The European Macroseismic Scale EMS 1998 (Grünthal 1998) classification has been used in order to identify the damage to the portfolio of masonry buildings considered. The grades of damage are defined as follows: Grade 0 (D0): No damage; Grade 1 (D1): Negligible to slight damage; Grade 2 (D2): Moderate damage; Grade 3 (D3): Substantial to heavy damage; Grade 4 (D4): Very heavy damage; Grade 5 (D5): Destruction. Identification of the damage level for different buildings has been carried out through two alternative surveying techniques; namely satellite Imagery through Copernicus-EMS damage grading map and virtual survey. Copernicus-EMS: The Copernicus (https://www.copernicus.eu/en, last access 31/01/2019) EMS provides rapid assessment of the damages through generation of “damage grading” maps, made possible by comparing pre- and post-event satellite images. With specific reference to damage data observed for Amatrice Earthquake, Masi et al. (2017) indicate that satellite imagery is more reliable for damage grades higher than or equal to Grade 4 and less reliable for damage grades lower than Grade 4. Therefore, although it is quite efficient for performing damage survey for vast areas, its use for a complete damage classification should be ideally accompanied by other means of damage surveying such as field surveys and areal photos. Visual Survey: It is to note that the Copernicus-EMS maps provide the damage grades for the buildings. However, these maps do not provide information about the sub-division of buildings into different classes. Hence, with the objective of complementing the damage grading maps provided by Copernicus, the same portfolio of buildings is also subjected to visual survey. The visual survey was based on photography available from field trips done at the end of September 2016 (courtesy of G. Forte and A. Santo), videos provided by drone (Feliziani et al. 2017, taken days after the Amatrice Earthquake) for areas that were difficult to access and google street view (https://www.google.it/streetview, last access 01/06/2017, used only for verifying the building classes and not the damage grades). The visual survey done building by building is arguably one of the most reliable means for assessing the occurred damage. However, the visual survey carried out herein is limited by the number of available photos/videos.
3.1.2 Structural characterization of the buildings
The portfolio of surveyed buildings is limited to residential masonry buildings. Further breaking down was based on the recommendations in Rota et al. (2008). In Rota et al. (2008), the masonry buildings are classified according to four parameters: number of floors, the presence of tie rods or tie beams, the type of horizontal structure (flexible or rigid floors), regular or irregular masonry layout. However, due to limitations posed by visual survey, the breaking down into more detailed classifications within the portfolio of residential masonry buildings in this work is limited to two factors: the number of storeys and the presence of tie rods or ring beams. As highlighted in a study conducted by Sorrentino et al. (2019), around 1/5th (20%) of the portfolio of masonry buildings in Amatrice and Accumuli presented either ring beams or tie rods. This percentage is confirmed for Amatrice also in a pre-earthquake survey (Fumagalli et al. 2017). In fact, the percentage of buildings that present either tie rods or rings beams according to visual survey in the present study for Amatrice and Accumuli is 28% (for all the seven towns considered this percentage is around 24%). The pre- and post-earthquake surveys conducted in the literature (Fumagalli et al. 2017; Sorrentino et al. 2019) indicate that the two municipalities of Amatrice and Accumuli are for the most part consisted of masonry buildings with following characteristics: up to 5 storeys (with a prevalence of 3 storey buildings); multi-layered masonry walls in the ground floor; units approximately dressed or undressed; medium to low mortar quality. It seems that, for the most part, the surveyed buildings could be classified as irregular (due to insufficient connections between the horizontal and vertical structures, between the masonry layers and the low to medium strength mortar used). Therefore, it has been assumed that the surveyed portfolio consists of irregular masonry buildings. As a result, four distinct classes of masonry buildings have been defined:
-
(a)
Masonry buildings without tie rods and without ring beams with number of stories ≤ 2 (Masonry Buildings Class 1, MBC1);
-
(b)
Masonry buildings without tie rods and without ring beams with number of stories > 2 (Masonry Buildings Class 2, MBC2);
-
(c)
Masonry buildings with tie rods or ring beams with number of stories ≤ 2 (Masonry Buildings Class 3, MBC3);
-
(d)
Masonry buildings with tie rods or ring beams with number of stories > 2 (Masonry Buildings Class 4, MBC4).
MBC1 can be associated to the union of classes IMA2 (Masonry-irregular layout-flexible floors-without tie rods and ring beams-1 or 2 stories) and IMA4 (Masonry-irregular layout-rigid floors-without tie rods and ring beams-1 or 2 stories) of Rota et al. (2008); MBC2 to the union of classes IMA6 (Masonry-irregular layout-flexible floors-without tie rods and ring beams-more than 2 stories) and IMA8 (Masonry-irregular layout-rigid floors-without tie rods and ring beams-more than 2 stories); MBS3 to the union of classes IMA1 (Masonry-irregular layout-flexible floors-with tie rods or ring beams-1 or 2 stories) and IMA3 (Masonry-irregular layout-rigid floors-with tie rods or ring beams-1 or 2 stories); MBC4 to union of classes IMA5 (Masonry-irregular layout-flexible floors-with tie rods or ring beams-more than 2 stories) and IMA7 (Masonry-irregular layout-rigid floors-with tie rods or ring beams-more than 2 stories). Figure 3 reports the numeric break-down of the portfolio of buildings into four classes, classified by the town to which they belong.
It is to note (Fig. 3) that the three main municipalities Amatrice, Accumuli and Arquata del Tronto show similar proportionality between the different masonry buildings classes considered. This confirms a rather homogenous structural characterization of the buildings in the largest centres of the studied area. Tables 2 and 3 report the number of buildings for each class and for each damage grade obtained based on Copernicus EMS and the visual survey, respectively. Figure 4 shows the damage grades provided by Copernicus for the surveyed buildings for the towns of (a) Amatrice, (b) Accumoli, (c) Arquata del Tronto, (d) Pescara del Tronto, (e) Illica, (f) Tino, and (g) Trisungo (the buildings belonging to damage levels D0 and D1 are grouped with label D0 + D1 in Fig. 4). It is worth noting that damage levels D0 and D1 could not be distinguished using the survey techniques adopted in this work. Consequently, fragility curves for D1 are not reported hereafter. It is to note that the surveyed building stock does not include all the masonry buildings of the seven towns under investigation. The buildings included in the survey are those for which both the damage level and the relative building class could be determined with confidence. From Figs. 3 and 4, it is possible to see that the surveyed buildings manage to provide a significant coverage of the entire portfolio. Moreover, comparisons of the survey results for the portfolio of buildings in this study with the damage patterns demonstrated in the recent literature (Fiorentino et al. 2018; Sorrentino et al. 2019) confirm compatibility in terms of the proportion of the buildings surveyed and the distribution of the different damage levels. More specifically, damage distribution for masonry buildings following the August 2016 earthquakes for Amatrice according to the visual survey demonstrates a fine agreement with Fiorentino et al. (2018) for D3 and higher. As far as it regards the Copernicus-EMS survey results for Amatrice only, the agreement is quite good for D4 and higher (see electronic supplementary material, Figure S1). In the work of Fiorentino et al. 2018, it is stated that 48.6% of the masonry buildings of the historical centre of Amatrice had collapsed; while 19.2% experienced partial collapse. Associating collapse to D5 and partial collapse to D4, the collapse and partial collapse percentages are 37.9% and 20.7% based on visual survey and 40.0% and 22.1% based on Copernicus. Figure S2 (in the electronic supplementary file) shows the histograms of number of cases > Di (i= 2:5) for visual survey and Copernicus for all four classes of masonry buildings considered. An overall fine agreement can be observed (especially for D3, D4 and D5).
3.2 Site effects for the surveyed buildings
Post-earthquake field recognition identified, as the most affected area, the valley of Tronto river. This valley is host to several municipalities and hamlets. The local geological and geomorphological setting can be sketched as shown in Fig. 5. Indeed, the villages occupy either the valley close to the river (e.g., Trisungo) or the cliffs overlooking it; with the latter being located usually at the top of small ridges and ancient erosional terraces (e.g. Amatrice, Accumoli, Arquata del Tronto, Fig. 5a) or located on the slopes (e.g. Pescara del Tronto, Illica, Tino, Fig. 5b).
The villages located on cliff-type morphology (Fig. 5a) are almost always bordered by steep slopes (25°–35°) with heights varying from 20 to 80 m. For these areas, buildings labelled with D4 and D5 damage grades are widespread and mainly localized near the steep escarpments and in the narrower part of the ridges as shown in Figs. 4b (Accumuli), 4c (Arquata del Tronto) and Fig. 5 part 1 (Amatrice). These buildings were affected by seismic waves’ focalization due to topographic shape effects (e.g., Grelle et al. 2018). These topographic effects were not present in lowland areas of the valley which suffered less damage, relatively speaking (see Trisungo, Fig. 4g). On the other hand, other towns shown in Figs. 4a (Amatrice), 4d (Pescara del Tronto), 4e (Illica), 4f (Tino), suffered widespread damage due to both topographic and stratigraphic effects (see Fig. 5 part 3). Some of these towns lie on slopes characterized by few meters of soft soils resting on a stiffer material (see Fig. 5b, Accumuli); where stiff arenaceous formation of the Laga Flysch is buried by few meters of weathered deposits and colluvium mainly made of silty sands. In the case of Pescara del Tronto (Fig. 4d), the hamlet lies on debris and travertine sands resting above a limestone bedrock. Figures S3 (in the electronic supplementary file) shows the overall damage break-down with respect to geologic units for MBC 1–2 and MBC 3–4, respectively. To account for stratigraphic amplification effects at a small scale, Eurocode 8 seismic soil classes are attributed to lithologic units. This is done according to the suggestions provided by Forte et al. (2017) for a comparable geo-lithological setting and is reported in Table 4.
Both adopted approaches (a) and (b) for considering the stratigraphic amplification are based on the soil classes. Therefore, they can be classified as 2-Grade level of microzonation following the guidelines proposed by ISSMGE (1999). In this study, the soil class was derived by the most recent soil class map of Italy published in Forte et al. (2019). This map was developed accounting for the shear-waves measurements distribution for each geological formation identified on a reinterpretation of 1:100.000 geological maps. The results are provided in a 500 × 500 m grid resolution, which is higher than the other available maps such as those in Michelini et al. (2008, Italian ShakeMap).
3.3 Generation of GMPE-based ground-shaking fields for the Central Italy Earthquake (24 August 2016) scenario
According to the procedure outlined in Sects. 2.2 and 2.3, Nsim= 1,000,000 realizations of the GMPE-based ground shaking fields are generated for the Amatrice Earthquake scenario (Mw= 6.0) providing the peak ground acceleration (PGA) at each building location. It is worth noting that the realizations are generated for the entire portfolio of buildings surveyed (i.e., all the classes together). Depending on the approach followed for considering the stratigraphic amplification, following steps are taken: for approach (a) outlined in Sect. 2.3.2, the ground-shaking field realizations are generated directly by employing the GMPE and considering the specific coefficients embedded in the GMPE for EC8 soil categories for stratigraphic amplification. These realizations are further modified in order to consider the topographic effect; for approach (b) outlined in Sect. 2.3.2, the generated ground-shaking fields on bedrock (i.e., assuming soil type A) are multiplied by both the stratigraphic and topographic amplification/deamplification factors calculated according to Eq. 9 and Table 1, respectively. It should be emphasized that the topographic effects are considered in the same manner for both approaches. Since the uncertainty in the evaluation of these amplification factors is not considered, their application affects only the median M (multiplied linearly by amplification factors) and leaves the covariance matrix Σ unaffected. Finally, in order to obtain the conditional GMPE-based fields, PGA registrations for eighty one accelerometric stations within a 100 km distance of the epicentre (see Fig. 6a and b) are employed in order to update the ground motion fields according to the procedure described in Sect. 2.3.3. Figure 6a shows the location of the 81 stations used for this work. Figure 6b shows the 16th, median, and 84th percentiles of the GMPE (Bindi et al. 2011), considering the stratigraphic (soil B, Landolfi et al. 2011) and topographic (curvature class 3, Silvestri et al. 2016) factors applied to the GMPE, before and after the updating with the available registrations. It should be noted that the results can be sensitive both to the position and the number of accelerometric stations considered. It is worth noting that: 1-theoretically-speaking, the recorded ground shaking information from all the accelerometric stations can be included in the updating; 2-the results are most sensitive to the accelerometric recordings inside the polygon delineating the surveyed buildings; 3- in the case of Amatrice earthquake, only one station is located within the polygon that delineates the surveyed buildings (see Fig. 6a). Therefore, in lieu of enough accelerometric stations within the polygon of the surveyed buildings, a 100 km distance limit from the epicentre was set. Roughly speaking, this marks the distance beyond which the ground shaking intensity is not going to be significant for shallow crustal earthquakes. Figure 6c maps the median PGA values for 25,000 realizations (note that smaller number of realizations are necessary for generating the map in Fig. 6c with respect to those necessary for estimating empirical fragility) rendered by a mesh-grid of 500 m × 500 m resolution. Figure 6d, demonstrates the median of the conditional GMPE-based ground-shaking field for the surveyed buildings in Amatrice (the largest of the seven towns considered).
3.4 The empirical fragility
3.4.1 Empirical fragility assessment for masonry building classes MBC1-MBC4
The empirical fragility curves and the corresponding plus/minus one standard deviation interval for each class of masonry buildings defined in Sect. 3.1.2 (MBC1–4) and for the different damage grades defined in Sect. 3.1.1 (EMS-98, D2–D5) have been calculated following the procedure described in Sects. 2.3 and 2.4 and plotted in Figs. 7 and 8 based on Copernicus EMS and visual survey (Sect. 3.1.1), respectively. In each sub-figure of Figs. 7 and 8, the fragility curves (from Eqs. 2 and 14) are plotted in thick solid lines (using established colour codes, i.e., green for D2, yellow for D3, orange for D4 and red for D5). It should be noted that the buildings marked as D0 and D1 are grouped together in this study. Therefore, no empirical fragility curves are reported for D1. The fragility curves are calculated based on the conditional ground-shaking fields generated as described in Sect. 2.3, corresponding to a total set of 1,000,000 simulations. Moreover, the expected value plus/minus one standard deviation fragility curves (Eq. 2) are plotted in dashed lines; the binned damage fractiles are plotted as circles. The observed damage fractiles, which are reported only for the sake of comparison with the resulting fragility curves, are obtained by dividing the PGA domain (ranging from smallest value to the maximum PGA value generated in the ground-shaking fields) into equidistant bins—making sure that each bin contains at least one surveyed building. For instance, the 0.50 damage fractile (shown in full circles) for damage grade Di and a given bin of PGA is calculated as the 50th percentile of the ratio of number of buildings in that bin whose damage grade is equal to or exceeds Di to the total number of buildings in the same bin over all the realizations. The damage fractiles for each bin are positioned at the centre of each bin. It can be observed that the fragility curves obtained through the logistic regression procedure (Sect. 2.4) show a very good agreement with the observed damage fractiles. As mentioned before, those ground-shaking fields that fail to produce monotonically increasing fragility curves are discarded as not plausible (see Sect. 2.3.4). Tables S1 and S2 (reported in the electronic supplement to the paper) summarize the equivalent lognormal median and logarithmic standard deviation (dispersion) for the fragility curves, respectively for Copernicus EMS and visual survey damage surveys.
3.4.2 Comparing the results based on ShakeMap (INGV) and on the proposed GMPE-based conditional ground-shaking fields
In this section, the empirical fragility curves obtained using the proposed conditional GMPE-based ground-shaking fields and the fragility curves obtained based on the INGV ShakeMap (http://shakemap.rm.ingv.it/shake/index.html, last access 15/12/2018) are compared for the building classes MBC1–4. It is to note that approach (b) as defined in Sect. 2.3.2 is used for consideration of the stratigraphic amplification. Figure 9 shows the fragility curves obtained using the INGV ShakeMap based on both Copernicus and visual survey and plotted with solid and dashed grey lines, respectively. The figure shows the empirical fragility curves (for D2–D5 and based on the standard colour code) obtained by employing the proposed conditional ground-shaking fields based on Copernicus and visual survey (plotted as solid and dashed lines). The consideration of ShakeMap for estimation of ground shaking (as it is common in the literature) shows significant discrepancies between the fragility curves obtained based on Copernicus and visual survey for MBC1 (D2 to D4), MBC2 (D2 to D3), MBC4(D2). This difference can be attributed to the fact that areal ground shaking values according to ShakeMap consisted of a maximum of six distinct levels of ground shaking. In some cases, differences in damage level attribution has led to significant changes in the number of buildings associated to each of the six ground shaking levels. Consequently, this can lead to significant difference in the resulting fragility curves. On the other hand, the fragility curves calculated based on the conditional GMPE-based random fields are based on point-wise estimates of ground shaking. Therefore, they are going to be less sensitive to variations in damage grade assignments (see Figure S4 in the electronic supplementary material for more detailed analysis of MBC1, D2). In fact, the fragility curves obtained by employing the proposed conditional GMPE-based ground-shaking fields based on Copernicus and visual survey show better agreement. The empirical fragility curves based on visual survey (as DDamage) and ShakeMap (as PGA) seem in closer agreement with those obtained through the proposed procedure (MBC1: D2–D5, MC4: D2). The empirical fragility curves obtained based on Copernicus (as damage survey) and the ShakeMap (for ground shaking estimation) seem to overestimate, with respect to the proposed procedure, the structural capacity for all the damage levels and for all the structural classes studied. In other words, the fragility curves obtained by employing the generated conditional ground-shaking field demonstrate higher vulnerability with respect to the fragility curves obtained based on the ShakeMap (paired up with Copernicus) for all four classes (the equivalent lognormal statistics of the ShakeMap-based fragility curves are reported in the electronic supplement to the paper, Table S3).
3.4.3 Stratigraphic amplification: Bindi et al. (2011) GMPE versus Landolfi et al. (2011)
In this section, the empirical fragility curves obtained based on ground shaking fields PGA that consider the effects of stratigraphic amplification as in Bindi et al. (2011) GMPE (approach (a), Sect. 2.3.2) and as proposed in Landolfi et al. (2011) model (approach (b), Sect. 2.3.2) are compared. The Copernicus-EMS damage grading maps are used as observed damage data (DDamage). As far as it regards approach (b), Bindi et al. (2011, ITA10) attenuation law is used to propagate the ground-shaking up to bedrock and the coefficients provided in Landolfi et al. (2011) are used to consider the stratigraphic amplification. Figure 10 shows the comparison between the fragility curves obtained using the stratigraphic coefficients from Bindi et al. (2011) (continuous grey lines) and Landolfi et al. (2011) (continuous lines, standard damage grade colours) for masonry buildings MBC1-MBC4 and for the damage levels D2–D5. The curves are obtained based on the proposed procedure in Sects. 2.3 and 2.4 in both cases; the only difference regards the factor accounting for stratigraphic amplification. It is to note that, for the fragility curves obtained based on Landolfi et al. (2011) stratigraphic factors are the same as those plotted in Fig. 7. It can be observed that the two approaches adopted for considering the stratigraphic amplification lead to a good agreement. No significant improvement is observed, in terms of the relative predictive capacity of the empirical fragility curves, when the Landolfi et al. (2011) model with the ground shaking-dependent non-linear soil column effect is adopted instead of the GMPE-embedded constant coefficients. More specifically, the use of Landolfi et al. (2011) model leads to a slight reduction both in the median and standard deviation with respect to Bindi et al. (2011) model for all the four building classes considered (with the exception of MBC3 where approach (a) fragility curves demonstrate slightly lower dispersion). Therefore, in lieu of site-specific information on stratigraphic amplification, the use of approach (a) based on ITA10-embedded coefficients is recommended. The equivalent lognormal statistics of the approach (a)-based fragility curves are reported in the electronic supplement to the paper, Table S4.
3.4.4 Comparison with literature
Figure 11 shows a comparison with empirical fragility curves proposed in Rota et al. (2008). The figure reports fragility curves (solid lines with the standard colour code for various damage levels) together with the corresponding plus/minus one-standard deviation confidence intervals (dashed lines), calculated according to the methodology proposed in Sect. 2 based on Copernicus-EMS damage grading maps (stratigraphic coefficients based on approach (b) Sect. 2.3.2). The fragility curves reported in Rota et al. (2008) (shade and dotted black lines) are plotted for analogous irregular masonry building classes. MBC1 can be seen as analogous to the union of classes IMA2 (dashed) and IMA4 (dotted); MBC2 as the union of classes IMA6 (dashed) and IMA8 (dotted); MBC3 as the union of IMA1 (dashed) and IMA3 (dotted); and MBC4 as the union of classes IMA5 (dashed) and IMA7 (dotted). The hazard curve (dotted light grey curve) for the town of Amatrice (Meletti et al. 2007, http://esse1.mi.ingv.it/d2.html, last access 15/01/2019) is also plotted in Fig. 11. It is to note that, the Bindi et al. (2011) GMPE furnishes estimates of the geometric mean between maxima of the two horizontal components of the ground shaking IMg.m.. This is while, the national hazard curves (Meletti et al. 2007; Meletti and Montaldo 2007, Project S1) are expressed in terms of the maximum of one arbitrary horizontal component IMarb. Since the empirical fragility curves obtained in this work are based on ground shaking estimates from the GMPE (i.e., IMg.m.), it is important to ensure they are consistent with the hazard curves in terms of intensity measure (Baker and Cornell 2006). That is, the risk convolution requires that hazard and fragility curves are expressed in terms of the same intensity measure. Based on a set of working assumptions, Baker and Cornell (2006) demonstrate that the expected value of the logarithm of IMg.m. is the same as that of IMarb. However, the logarithmic standard deviation of IMg.m. is lower than that or equal to that of IMarb, depending on the correlation between the maxima of ground shaking for the two horizonal components. A constant correction factor (estimated in Ebrahimian et al. 2019) has been applied to both inter-event and intra-event components of logarithmic standard deviation of the GMPE to make the fragility and hazard curves consistent (for the purpose of risk convolution).Table 5 reports the annual probabilities of exceeding the different damage levels (risk R). RF denotes the risk calculated based on fragility curve F; RF±σ denote the risk calculated based on fragility plus/minus one standard deviation curves, respectively. Except for MBC2 (D5), MBC3 (D3–D5) and MBC4 (D4–D5), the overall agreement observed between the fragility curves is fine (when compared to Rota et al. 2008 fragility curves with flexible floors). The quantified comparison in Table 5 reveals that the annual probability of damage level exceedance estimates rendered by fragility curves of Rota et al. 2008 and the ones proposed in study are quite different (although some agreement in ballpark can be observed). The following factors can be responsible for the lack of a good agreement: differences in building portfolio characteristics, differences in the definitions of the building classes, the methods used for estimating the ground-shaking intensity, and the methods used for fitting the damage data. Furthermore, the results in Table 5 show an overall agreement between visual survey and Copernicus damage gradings paired up the conditional GMPE-based ground shaking fields (the first two rows of the table). Moreover, the results in Table 5 indicate that the fragility curves rendered by pairing up of visual survey and INGV ShakeMap provide un-conservative estimates with respect to both Rota et al. 2008 and conditional GMPE-based ground shaking fields. This provides evidence for the remarkable sensitivity of the empirical fragility curves to the methods used for the estimation of ground shaking level.
4 Conclusions
Empirical fragility assessment can reveal significant sensitivity to the definition of the building class, building-to-building variability within a given class, observed damage and ground shaking levels for the considered portfolio of buildings for a given class. In such context, accurate estimation of the ground shaking level at the location of buildings and considering the possible correlations between these ground shaking levels is as important as accurate estimation of the damage. The (Bayesian) concept of robust updated fragilities (a.k.a. predictive fragilities) can be applied to consider the uncertainty in the estimation of the ground shaking level at the position of the buildings of the interest. The joint probability distribution of the ground shaking levels—and hence the spatial correlation structure—can be modelled through simulation of plausible realizations of a random ground motion field. This random field is defined based on the (updated) statistics of the GMPE considering the ground shaking levels recorded at the nearby accelerometric stations. The set of plausible ground motion fields is further filtered as those fields which lead to physically meaningful fragility curves for all damage levels. This work uses damage data obtained based on Copernicus-EMS damage grading maps and complemented by visual survey for buildings damaged due to the Amatrice Earthquake of 24th of August 2016 and the aftershocks immediately following it. Using Copernicus EMS and visual survey after the Amatrice Earthquake of 24 August 2016 as data basis, the paper explores several factors that may severely affect empirical fragility assessment.
Effect of stratigraphic amplification It can be observed that explicit consideration of the non-linear soil column amplification by using (deterministic) class-specific coefficients applied to ground shaking at bedrock level (approach b, Sect. 2.3.2) shows very slight improvement with respect to embedded coefficients of ITA10 GMPE (approach a, Sect. 2.3.2). This could be somewhat expected because, technically speaking, both approaches considered in this work can be classified by the same level of micro-zonation (since they are both based on soil class). Therefore, in lieu of site-specific information on stratigraphic amplification, the use of the embedded coefficients of ITA10 is recommended.
Comparison between Copernicus and Visual Survey The two approaches seem to be in good agreement in terms of the cumulative probability of exceeding a damage state for damage states of D3 and higher. It can be observed that the two approaches for damage estimation lead to comparable results in ballpark in terms of the annual probability of exceeding a specific damage level. Nevertheless, the Copernicus-EMS database does not provide the building classification; therefore, it needs to rely on complementary information regarding the building classification.
Comparison with ShakeMap While Copernicus and visual survey demonstrate overall fine agreement with each other when the proposed procedure is followed for the estimation of the ground shaking levels, the two damage grading datasets lead in some cases to totally different results when paired up with ShakeMap. This can be attributed to the fact that ShakeMap relies on areal estimates of ground shaking and does not consider the uncertainty in the estimation of ground-shaking at the position of buildings of interest. In the context of the survey performed in this work, pairing up Copernicus damage grades with ShakeMap leads to significant underestimation of vulnerability and risk with respect to the conditional GMPE-based ground shaking fields.
Comparison with literature Comparison with the fragility curves reported in Rota et al. 2008 for analogous classes of masonry buildings, shows agreement in the ballparks of the probability of exceeding various damage levels for buildings without tie rods or ring beams (MBC1–2).
The results demonstrated are subjected to the following limitations:
-
In some cases, the adopted survey techniques might not have assigned the correct building class. For example, the ring beams may not be visible from the building façade.
-
There might be some bias in the ratio of un-damaged buildings to the total number of buildings surveyed. This can lead to a bias in the survey sample size. The potential bias can be attributed to the nature of the survey that tended to identify damaged buildings.
-
The survey coverage is not full. It is assumed that the surveyed portfolio of buildings provides a sufficient and un-bias coverage of the portfolio of buildings in the seven towns considered.
-
There was only one accelerometric stations inside the polygon defining the portfolio of buildings considered (see Fig. 6b).
-
It is assumed that the damage survey is not subjected to error. However, if there is information available for calibrating the accuracy of damage results, the predictive fragilities can be extended to also consider the error in damage survey.
-
The choice of the GMPE and the corresponding spatial correlation structure between intra-events residuals can affect the results. The results could potentially improve, if site-specific data was available for the stratigraphic amplification. In other words, the fragility curves would be more accurate if based on higher levels of seismic microzonation.
-
The empirical fragility parameters are estimated without enforcing the ordinal hierarchy in the damage scale.
As a final word, the proposed method provides practical means for considering the uncertainties in the estimated ground shaking level and correlations in the ground motion at the location of the damaged buildings. The statistical updating considering recorded ground motion constitutes a fundamental step in incorporating the local site conditions and the ground motion levels effectively experienced. After being subjected to sufficient validation and testing, the proposed methods can be considered as a potentially valid alternative to the use of ShakeMap for estimation of the ground shaking level, for the purpose of empirical fragility estimation. Moreover, the proposed method can incorporate damage estimation methods using remote sensing techniques for fast fragility estimation.
References
Baggio C, Bernardini A, Colozza R, Corazza L, Della Bella M, Di Pasquale G, et al (2007) Field manual for post-earthquake damage and safety assessment and short term countermeasures (AeDES). European Commission—Joint Research Centre—Institute for the Protection and Security of the Citizen, EUR, 22868
Baker JW (2015) Efficient analytical fragility function fitting using dynamic structural analysis. Earthq Spectra 31(1):579–599
Baker JW, Cornell CA (2006) Which spectral acceleration are you using? Earthq Spectra 22(2):293–312
Beck JL, Au SK (2002) Bayesian updating of structural models and reliability using Markov chain Monte Carlo simulation. J Eng Mech 128(4):380–391
Bindi D, Pacor F, Luzi L, Puglia R, Massa M, Ameri G, Paolucci R (2011) Ground motion prediction equations derived from the Italian strong motion database. Bull Earthq Eng 9(6):1899–1920
Bindi D, Massa M, Luzi L, Ameri G, Pacor F, Puglia R, Augliera P (2014a) Pan-European ground-motion prediction equations for the average horizontal component of PGA, PGV, and 5%-damped PSA at spectral periods up to 3.0 s using the RESORCE dataset. Bull Earthq Eng 12(1):391–430
Bindi D, Massa M, Luzi L, Ameri G, Pacor F, Puglia R, Augliera P (2014b) Erratum to: Pan-European ground-motion prediction equations for the average horizontal component of PGA, PGV, and 5%-damped PSA at spectral periods up to 3.0 s using the RESORCE dataset. Bull Earthq Eng 12(1):431–448
Braga F, Dolce M, Liberatore D (1982) A statistical study on damaged buildings and an ensuing review of the MSK-76 scale. In: 7th European conference on earthquake engineering, Athens, Greece, pp 431–450
Calvi GM, Pinho R, Magenes G, Bommer JJ, Restrepo-Vélez LF, Crowley H (2006) Development of seismic vulnerability assessment methodologies over the past 30 years. ISET J Earthq Technol 43(3):75–104
Charvet I, Ioannou I, Rossetto T, Suppasri A, Imamura F (2014) Empirical fragility assessment of buildings affected by the 2011 Great East Japan tsunami using improved statistical models. Nat Hazards 73(2):951–973
Colombi M, Borzi B, Crowley H, Onida M, Meroni F, Pinho R (2008) Deriving vulnerability curves using Italian earthquake damage data. Bull Earthq Eng 6(3):485–504
Crowley H, Bommer JJ, Stafford PJ (2008) Recent developments in the treatment of ground-motion variability in earthquake loss models. J Earthq Eng 12(S2):71–80
De Luca F, Verderame GM, Manfredi G (2015) Analytical versus observational fragilities: the case of Pettino (L’Aquila) damage data database. Bull Earthq Eng 13(4):1161–1181
De Luca F, Woods GE, Galasso C, D’Ayala D (2018) RC infilled building performance against the evidence of the 2016 EEFIT Central Italy post-earthquake reconnaissance mission: empirical fragilities and comparison with the FAST method. Bull Earthq Eng 16(7):2943–2969
De Risi R, Goda K, Mori N, Yasuda T (2017) Bayesian tsunami fragility modeling considering input data uncertainty. Stoch Environ Res Risk Ass 31(5):1253–1269
Del Gaudio C, Ricci P, Verderame GM, Manfredi G (2016) Observed and predicted earthquake damage scenarios: the case study of Pettino (L’Aquila) after the 6th April 2009 event. Bull Earthq Eng 14(10):2643–2678
Del Gaudio C, De Martino G, Di Ludovico M, Manfredi G, Prota A, Ricci P, Verderame GM (2017) Empirical fragility curves from damage data on RC buildings after the 2009 L’Aquila earthquake. Bull Earthq Eng 15(4):1425–1450
Der Kiureghian A (2005) Non-ergodicity and PEER’s framework formula. Earthq Eng Struct Dyn 34(13):1643–1652
Di Pasquale G, Orsini G, Romeo RW (2005) New developments in seismic risk assessment in Italy. Bull Earthq Eng 3(1):101–128
Dolce M, Goretti A (2015) Building damage assessment after the 2009 Abruzzi earthquake. Bull Earthq Eng 13(8):2241–2264
Dolce M, Masi A, Marino M, Vona M (2003) Earthquake damage scenarios of the building stock of Potenza (Southern Italy) including site effects. Bull Earthq Eng 1(1):115–140
Ebrahimian H, Jalayer F (2017) Robust seismicity forecasting based on Bayesian parameter estimation for epidemiological spatio-temporal aftershock clustering models. Sci Rep 7(1):9803
Ebrahimian H, Jalayer F, Forte G, Convertito V, Licata V, d’Onofrio A, Santo A, Silvestri F, Manfredi G (2019) Site-specific probabilistic seismic hazard analysis for the Western area of Naples, Italy. Bull Earthq Eng 17(9):4743–4796
Esposito S, Iervolino I (2012) Spatial correlation of spectral acceleration in European data. Bull Seismol Soc Am 102(6):2781–2788
Eurocode 8 Part 3 (2007) Design of structures for earthquake resistance-Assessment and retrofitting of buildings
Feliziani F, Lorusso O, Ricci A, Massabò A, Di Lolli A, Colangeli A, Fiorini M, (2017) 2D and 3D models in emergency scenarios: UAVs for planning search and rescue operations and for preliminary assessment of high buildings. In: UAV and SAR 2017-International workshop, Rome, 29th March, 2017. http://www.vigilfuoco.it/aspx/isaViewDoc.aspx?id=222&t=1; youtube link: https://www.youtube.com/watch?v=aekrS7MOcCQ
Fiorentino G, Forte A, Pagano E, Sabetta F, Baggio C, Lavorato D et al (2018) Damage patterns in the town of Amatrice after August 24th 2016 Central Italy earthquakes. Bull Earthq Eng 16(3):1399–1423
Forte G, Fabbrocino S, Fabbrocino G, Lanzano G, Santucci de Magistris F, Silvestri F (2017) A geolithological approach to seismic site classification: an application to the Molise Region (Italy). Bull Earthq Eng 15(1):175–198
Forte G, Chioccarelli E, De Falco M, Cito P, Santo A, Iervolino I (2019) Seismic soil classification of Italy based on surface geology and shear-wave velocity measurements. Soil Dyn Earthq Eng 122:79–93
Fumagalli F, Liberatore D, Monti G, Sorrentino L (2017). Building features of Accumoli and Amatrice in a pre-earthquake survey. Building features of Accumoli and Amatrice in a pre-earthquake survey, pp 44–54
Garcìa-Rodriguèz MJ, Havenith HB, Benito B (2008) Evaluation of Earthquake- triggered landslides in El Salvador using a GIS based Newmark model. In: 14th conference on earthquake engineering, Beijing, China
Goda K (2011) Interevent variability of spatial correlation of peak ground motions and response spectra. Bull Seismol Soc Am 101(5):2522–2531
Goda K, Atkinson GM (2010) Intraevent spatial correlation of ground-motion parameters using SK-net data. Bull Seismol Soc Am 100(6):3055–3067
Goda K, Hong HP (2008) Spatial correlation of peak ground motions and response spectra. Bull Seismol Soc Am 98(1):354–365
Grelle G, Bonito L, Maresca R, Maufroy E, Revellino P, Sappa G, Guadagno F (2018) Topographic effects in Amatrice suggested from the SISERHMAP predictive model, seismic data and damage. In: 16th European conference on earthquake engineering, Thessaloniki, 18–21 June 2018
Grünthal G (1998) Cahiers du Centre Europe´en de Ge´odynamique et de Se´ismologie: volume 15—European Macroseismic Scale 1998. European Center for Geodynamics and Seismology, Luxembourg
Hsieh MH, Lee BJ, Lei TC, Lin JY (2013) Development of medium-and low-rise reinforced concrete building fragility curves based on Chi–Chi Earthquake data. Nat Hazards 69(1):695–728
Ioannou I, Douglas J, Rossetto T (2015) Assessing the impact of ground-motion variability and uncertainty on empirical fragility curves. Soil Dyn Earthq Eng 69:83–92
Ioannou I, Chandler RE, Rossetto T (2020) Empirical fragility curves: the effect of uncertainty in ground motion intensity. Soil Dyn Earthq Eng 129:105908
ISSMGE (1999) Manual for zonation on seismic geotechnical hazards. In: Technical committee for earthquake geotechnical engineering, TC4, international society for soil mechanics and geotechnical engineering. The Japanese Geot Soc, Tokyo
Jaiswal K, Wald D, D’Ayala D (2011) Developing empirical collapse fragility functions for global building types. Earthq Spectra 27(3):775–795
Jalayer F, Ebrahimian H (2017) Seismic risk assessment considering cumulative damage due to aftershocks. Earthq Eng Struct Dyn 46(3):369–389
Jalayer F, Ebrahimian H (2020) Seismic reliability assessment and the non-ergodicity in the modelling parameter uncertainties. Earthq Eng Struct Dyn 49(5):434–457
Jalayer F, Iervolino I, Manfredi G (2010) Structural modeling uncertainties and their influence on seismic assessment of existing RC structures. Struct Saf 32(3):220–228
Jalayer F, Carozza S, De Risi R, Manfredi G, Mbuya E (2016) Performance-based flood safety-checking for non-engineered masonry structures. Eng Struct 106:109–123
Jalayer F, Ebrahimian H, Miano A, Manfredi G, Sezen H (2017) Analytical fragility assessment using un-scaled ground motion records. Earthq Eng Struct Dyn 46(15):2639–2663
Jayaram N, Baker JW (2009) Correlation model for spatially distributed ground-motion intensities. Earthq Eng Struct Dyn 38(15):1687–1708
Lagomarsino G, Giovinazzi S (2006) Macroseismic and mechanical models for the vulnerability and damage assessment of current buildings. Bull Earthq Eng 4:415–443
Lallemant D, Kiremidjian A, Burton H (2015) Statistical procedures for developing earthquake damage fragility curves. Earthq Eng Struct Dyn 44(9):1373–1389
Landolfi L, Caccavale M, d’Onofrio A, Silvestri F, Tropeano G (2011) Preliminary assessment of site stratigraphic amplification for ShakeMap processing. In: V international conference on earthquake geotechnical engineering, Santiago, Chile, p 5.3
Lanzano G, Pacor F, Luzi L, D’Amico M, Puglia R, Felicetta C (2017) Systematic source, path and site effects on ground motion variability: the case study of Northern Italy. Bull Earthq Eng 15(11):4563–4583
Lanzano G, Luzi L, Pacor F, Felicetta C, Puglia R, Sgobba S, D’Amico M (2019) A revised ground-motion prediction model for shallow crustal earthquakes in Italy. Bull Seism Soc Am 109(2):525–540
Liel A, Lynch K (2012) Vulnerability of reinforced-concrete-frame buildings and their occupants in the 2009 L’Aquila, Italy, earthquake. Nat Hazards Rev 13(1):11–23
Masi A, Chiauzzi L, Santarsiero G, Liuzzi M, Tramutoli V (2017) Seismic damage recognition based on field survey and remote sensing: general remarks and examples from the 2016 Central Italy earthquake. Nat Hazards 86(1):193–195
Meletti C, Montaldo V (2007) Stime di pericolosità sismica per diverse probabilità di superamento in 50 anni: valori di ag. Project DPC-INGV S1, deliverable D2 (in Italian). Istituto Nazionale di Geofisica e Vulcanologia—Sezione di Milano-Pavia. http://esse1.mi.ingv.it/d2.html. Accessed 15 Jan 2019
Meletti C, Calvi GM, Stucchi M (2007) Project S1—continuation of assistance to DPC for improving and using the seismic hazard map compiled according to the Prime Minister “Ordinanza” 3274/2003 and planning future initiatives—final report, INGV (in Italian). Interactive maps of seismic hazard (WebGis). http://esse1.mi.ingv.it/. Accessed 15 Jan 2019
Miano A, Jalayer F De Risi R, Prota A, Manfredi G (2015) A case-study on scenario-based probabilistic seismic loss assessment for a portfolio of bridges. In: 12th international conference on applications of statistics and probability in civil engineering, Vancouver, Canada, July 12–15, 2015
Miano A, Jalayer F, De Risi R, Prota A, Manfredi G (2016) Model updating and seismic loss assessment for a portfolio of bridges. Bull Earthq Eng 14(3):699–719
Miano A, Jalayer F, Ebrahimian H, Prota A (2018) Cloud to IDA: efficient fragility assessment with limited scaling. Earthq Eng Struct Dyn 47(5):1124–1147
Michelini A, Faenza L, Lauciani V, Malagnini L (2008) ShakeMap implementation in Italy. Seismol Res Lett 79(5):688–697
Michelini A, Faenza L, Lanzano G, Lauciani V, Jozinović D, Puglia R, Luzi L (2019) The new shakemap in Italy: progress and advances in the last 10 Yr. Seismol Res Lett. https://doi.org/10.1785/0220190130
Mosleh A, Apostolakis G (1986) The assessment of probability distributions from expert opinions with an application to seismic fragility curves. Risk Anal 6(4):447–461
Orsini G (1999) A model for buildings’ vulnerability assessment using the Parameterless Scale of Seismic Intensity (PSI). Earthq Spectra 15(3):463–483
Paolucci R (2002) Amplification of earthquake ground motion by steep topographic irregularities. Earthq Eng Struct Dyn 31(10):1831–1853
Papadimitriou C, Beck JL, Katafygiotis LS (2001) Updating robust reliability using structural test data. Probab Eng Mech 16(2):103–113
Park J, Bazzurro P, Baker JW (2007) Modeling spatial correlation of ground motion intensity measures for regional seismic hazard and portfolio loss estimation. Applications of statistics and probability in civil engineering. Taylor & Francis Group, London, pp 1–8
Porter K, Kennedy R, Bachman R (2007) Creating fragility functions for performance-based earthquake engineering. Earthq Spectra 23(2):471–489
Rossetto T, Elnashai A (2003) Derivation of vulnerability functions for European-type RC structures based on observational data. Eng Struct 25(10):1241–1263
Rossetto T, Ioannou I, Grant DN, Maqsood T (2014) Guidelines for empirical vulnerability assessment, Technical Report 2014-08 V1.0.0, 140 pp, GEM Foundation, Pavia, Italy
Rosti A, Rota M, Penna A (2018) Damage classification and derivation of damage probability matrices from L’Aquila (2009) post-earthquake survey data. Bull Earthq Eng 16(9):3687–3720
Rota M, Penna A, Strobbia CL (2008) Processing Italian damage data to derive typological fragility curves. Soil Dyn Earth Eng 28(10–11):933–947
Sabetta F, Pugliese A (1996) Estimation of response spectra and simulation of nonstationary earthquake ground motions. Bull Seismol Soc America 86(2):337–352
Sabetta F, Goretti A, Lucantoni A (1998) Empirical fragility curves from damage surveys and estimated strong ground motion. In: 11th European conference on earthquake engineering, Paris, France, September 6–11
Sanchez-Sesma FJ (1990) Elementary solutions for response of a wedge-shaped medium to incident SH and SV waves. Bull Seismol Soc Am 80:737–742
Sasani M, Kiureghian AD (2001) Seismic fragility of RC structural walls: displacement approach. J Struct Eng 127(2):219–228
Sextos A, De Risi R, Pagliaroli A et al (2018) Local site effects and incremental damage of buildings during the 2016 Central Italy earthquake sequence. Earthq Spectra 34(4):1639–1669
Sgobba S, Lanzano G, Pacor F, Puglia R, D’Amico M, Felicetta C, Luzi L (2019) Spatial correlation model of systematic site and path effects for ground-motion fields in Northern Italy. Bull Seism Soc Am 109(4):1419–1434
Shinozuka M, Feng MQ, Lee J, Naganuma T (2000) Statistical analysis of fragility curves. J Eng Mech 126(12):1224–1231
Silvestri F, Forte G, Calvello M (2016) Multi-level approach for zonation of seismic slope stability: experiences and perspectives in Italy. Landslides and engineered slopes. Experience, theory and practice. In: Aversa et al (eds) 2016 Associazione Geotecni
Singhal A, Kiremidjian AS (1998) Bayesian updating of fragilities with application to RC frames. J Struct Eng 124(8):922–929
Sokolov V, Wenzel F (2011) Influence of spatial correlation of strong ground motion on uncertainty in earthquake loss estimation. Earthq Eng Struct Dyn 40(9):993–1009
Sokolov V, Wenzel F (2013) Spatial correlation of ground motions in estimating seismic hazards to civil infrastructure. In: Handbook of seismic risk analysis and management of civil infrastructure systems. Woodhead Publishing Series in Civil and Structural Engineering 2013, pp 57–78
Sorrentino L, Cattari S, da Porto F, Magenes G, Penna A (2019) Seismic behaviour of ordinary masonry buildings during the 2016 central Italy earthquakes. Bull Earthq Eng 17(10):5583–5607
Spence RJS, Coburn AW, Pomonis A, Sakai S (1992) Correlation of ground motion with building damage: the definition of a new damage-based seismic intensity scale. In: 10th world conference on earthquake engineering, Madrid, Spain, vol 1, pp 551–556
Straub D, Der Kiureghian A (2008) Improved seismic fragility modeling from empirical data. Struct Safety 30(4):320–336
Torgoev A, Havenith HB, Lamair L (2013) Improvement of seismic landslide susceptibility assessment through consideration of geological and topographic amplification factors. JAG 2013, 17–18/09 Grenoble, France
Tropeano G, Soccodato FM, Silvestri F (2018) Re-evaluation of code-specified stratigraphic amplification factors based on Italian experimental records and numerical seismic response analyses. Soil Dyn Earthq Eng 110:262–275
Vanmarcke EH, Heredia-Zavoni E, Fenton GA (1993) Conditional simulation of spatially correlated earthquake ground motion. J Eng Mech 119(11):2333–2352
Weatherill GA, Silva V, Crowley H, Bazzurro P (2015) Exploring the impact of spatial correlations and uncertainties for portfolio analysis in probabilistic seismic loss estimation. Bull Earthq Eng 13(4):957–981
Wiesberg S (1985) Applied linear regression. John Wiley, New York
Yazdi AJ, Haukaas T, Yang T, Gardoni P (2016) Multivariate fragility models for earthquake engineering. Earthq Spectra 32(1):441–461
Yazgan U (2015) Empirical seismic fragility assessment with explicit modeling of spatial ground motion variability. Eng Struct 100:479–489
Zuccaro G, Cacace F (2009) Revisione dell’inventario a scala nazionale delle classi tipologiche di vulnerabilita` ed aggiornamento delle mappe nazionali di rischio sismico. XIII Convegno ANIDIS ‘‘L’ingegneria sismica in Italia’’, June 28–July 2, Bologna, Italy. Paper S14.39 (in Italian)
Zuccaro G, Cacace F (2015) Seismic vulnerability assessment based on typological characteristics. The first level procedure ‘‘SAVE’’. Soil Dyn Earthq Eng 69:262–269
Acknowledgements
This work is based on visual damage survey results described in the B.Sc. thesis of Filomena Castagna. This support is gratefully acknowledged. Moreover, the visit of the third and fourth authors on the stricken area was funded by the Italian Civil Protection (DPC) Progetto RELUIS 2016 WP4 - Approcci a larga scala per la previsione della risposta locale e dell’instabilità del sottosuolo. (Scientific Coordinator Prof. Antonio Santo). The authors would also like to acknowledge the anonymous reviewers who have contributed significantly to improving and enriching the paper.
Funding
Open access funding provided by Università degli Studi di Napoli Federico II within the CRUI-CARE Agreement.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Miano, A., Jalayer, F., Forte, G. et al. Empirical fragility assessment using conditional GMPE-based ground shaking fields: application to damage data for 2016 Amatrice Earthquake. Bull Earthquake Eng 18, 6629–6659 (2020). https://doi.org/10.1007/s10518-020-00945-6
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10518-020-00945-6