Development of Empirical Fragility Curves in Earthquake Engineering considering Nonspecific Damage Information

As a function of fragility curves in earthquake engineering, the assessment of the probability of exceeding a specific damage state according to the magnitude of earthquake can be considered. Considering that the damage states for fragility curves are generally nested to each other, the possibility theory, a special form of the evidence theory for nested intervals, is applied to generate fragility information from seismic damage data. While the lognormal distributions are conventionally used to generate fragility curves due to their simplicity and applicability, themethodology to use the possibility theory does not require the assumption of distributions. Seismic damage data classified by four damage levels were used for a case study. ,e resulted possibility-based fragility information expressed by two monotone measures, “possibility” and “certainty,” are compared with the conventional fragility curves based on probability. ,e results showed that the conventional fragility curves provide a conservative estimation at the relatively high earthquake magnitude compared with the possibility-based fragility information.


Introduction
In earthquake engineering, fragility curves have been used to estimate damages of infrastructures according to the magnitude of earthquake.While fragility curves can provide reasonable estimation of earthquake damages with damage levels, those might neglect the possible slight damage occurrence at the relatively low earthquake magnitude due to the nature of probability density functions consists of two parameters, expected mean value and dispersion.
Traditionally, probability theory has been used to model uncertainties in structural engineering, especially when addressing reliability for structural safety [1,2].However, the types of uncertainties considered in probability theory are random, chance, and likelihood, and there are limitations to model other types of uncertainties such as nonspecificity, fuzziness, and strife, using probability theory [3,4].Random uncertainty known as aleatory uncertainty is from inherent randomness and therefore is irreducible.However, other types of uncertainties known as epistemic uncertainties arise from lack of knowledge and therefore are reducible and subjective.Research on generalized information theory (GIT) [5,6] showed that three types of epistemic uncertainties due to lack of knowledge and/or variability thrive when modeling complex environments [7].While nonspecificity represents the difficulty to choose from many modeling alternatives, fuzziness represents the uncertainty due to lack of sharpness (imprecise boundaries) of the modeling parameter.Strife expresses the uncertainty due to conflict among alternatives.Given are there various types of uncertainties, appropriate modeling of uncertainty has been an interesting and challenging topic in many areas during the last few decades [8][9][10].A number of theories to model uncertainties adequately have been introduced: evidence theory [11,12], possibility theory [3,13] and fuzzy set theory [14,15].
In earthquake engineering, empirical fragility curves were generally presented in the form of lognormal cumulative distribution function (CDF) with respect to peak ground acceleration (PGA) representing the ground motion intensity due to earthquake [16].To generate empirical fragility curve for a damage state, the damage reports by experts are used.e damage reports usually present the damage states in linguistic ways such as "no damage," "slight damage," "moderate damage," "extensive damage," and "collapse" for a structure experiencing earthquake of a PGA.As the fragility curve of a damage state represents the fragility of "at least" of the damage level, the evidence for a damage state includes possible higher damage states.For example, the evidence for "moderate damage" of a structure by an earthquake means that the structure is damaged at least moderately, and it might be possible for the structure to be damaged extensively or collapsed.As there exist ambiguous boundaries between the damage states, fuzzy logic and possibility theory were applied to resolve the ambiguity [17,18].
Recently, extensive earthquake damage data are used to generate fragility curves [19][20][21][22][23][24][25].Postearthquake surveys of approximately 340000 reinforced concrete structures were used to derive fragility curves for a European seismic risk assessment scenario [19].A database of 7597 reinforced concrete buildings located in the city and the province of L'Aquila in Italy was used in order to derive fragility curves [20].e observed damage to 9500 of low-rise residential buildings from earthquakes in South Iceland was studied by typological fragility curves [21].Moreover, fragility curves were developed from millions of data on the basis of 665,515 building damage cases by earthquake in Nepal [22,23].In Italy, the postearthquake damage surveys of approximately 90,000 buildings in order to derive fragility curves were considered [24,25].Even with the increase of damage data to generate fragility curves, there is still an uncertainty of nonspecificity, the difficulty to choose from many modeling functions of fragility curves such as lognormal, extreme type I, extreme type II functions, and so on.
In this study, the evidence of damage state is dealt with possibility theory.It is noticeable that the fragility curves from possibility distribution representing the certainty of damage state and those are generated without any assumption of distributions.erefore, there is no uncertainty of nonspecificity to choose functions of fragility curves.

Possibility Information in Fragility Curves
eories for modeling uncertainties present different types of uncertainty assignment and monotone measures.As uncertainty assignment terms, the degree of belief, probability distribution, and possibility distribution are used for evidence theory, probability theory, and possibility theory, respectively.To quantify the assigned uncertainties, monotone measures are used such as dual monotone measures of plausibility and belief, dual monotone measures of possibility and certainty, and single monotone measure of probability for evidence theory, probability theory, and possibility theory, respectively [26][27][28].Considering the relationship between uncertainty assignment terms and monotone measures used for each theory, it can be known as probability theory and possibility theory are special forms of evidence theory [6].Consider a discrete universe D that consists of a set of damage levels, where d N , d S , d M , d E , and d C represent no damage, slight damage, moderate damage, extensive damage, and collapse of a structure due to a seismic force level respectively.
In evidence theory, which is also known as Dempster-Shafer theory [11,12], the degree of belief m based on evidence is assigned to all countable subsets A (e.g., Ø, {d N }, {d N , d S }, ..., {d N , ..., d C }) with the constraint of (2) Dual monotone measures, belief bel (A) and plausibility pl (A), for a subset A are calculated as bel(A) �  B∈A m(B), (3) While belief measure represents the degree of evidence for a subset A, plausibility measure is defined as "Complement of the belief of the complement of a subset A" as pl(A) � 1 − bel(A). ( As belief measure is based on the degree of belief with its evidence, belief measure of "Complement of a subset A" also needs its evidence.erefore, if there is no evidence for "Complement of a subset A," one cannot determine the belief of "Complement of a subset A" as 1-bel (A).e difference between these two measures can represent our ignorance (lack of knowledge) of a subset A (denoted ign) as In probability theory, probability distribution, which is equivalent to the degree of belief m in evidence theory, is assigned to a single variable (e.g.d N , . .., d C ) on universe D such as where d i denotes the damage state, d N , d S , d M , d E , and d C .
Only one monotone measure, probability prob(A), for a subset A is defined as and probability measure of "Complement of a subset A" is defined as with the excluded middle axioms [4].Unlike evidence theory, "Complement of a subset A" can be determined as 1 -prob (A).erefore, our lack of knowledge measured by ign (A) in equation ( 6) cannot be measured in probability theory.
In possibility theory, possibility distribution π is assigned to a single damage level in possibility theory such as 2 Advances in Civil Engineering e relationship between the uncertainty assignment π and the degree of belief m in equation (10) indicates that possibility theory is a special form of evidence theory when the collective body of evidence is consonant [26][27][28] (see Figure 1).
Dual monotone measures, certainty cert (A) and possibility pos (A), for a subset A are determined as pos(A) sup π d i : is relationship can be converted to a single measure that represents the degree of con rmation C (A) of a subset A, whose range is from −1 to 1 [5]: Negative value of the degree of con rmation expresses the degree of discon rmation of a subset A. It is noticeable that the assignment of possibility distribution of 1 to a single variable means that the occurrence of the variable is possibly 1 but certainly "no evidence," while the assignment of degree of possibility of zero con rms that the occurrence of the variable is possibly zero and certainly zero.For a consonant body of evidence, the following relationship for two di erent subsets, A and B, can be proven [6]: When the ignorance in equation ( 6) is zero, the excluded middle axioms are satis ed and the evidence can be described using one monotone measure called the probability measure.Moreover, when body of evidence is consonant, monotone measures in evidence theory, belief and plausibility, can be represented as those in possibility theory, certainty and possibility, respectively (see Appendix for examples).
As the fragility curve of a damage state represents the fragility of "at least" of the damage level, probable damage information for higher damage state is included in a damage state.For example, the evidence for "no damage" means that a structure seems undamaged, but there might be possible damages for the structure: slight, moderate, extensive, or collapsed damages.Considering the characteristics of evidence, "at least," we can de ne subsets N, S, M, E, and C for "at least no damage," "at least slight damage," "at least moderate damage," "at least extensive damage," and "at least collapse" as shown in equation (16).It can be recognized that a damage subset includes higher than and equal to the damage levels.
It is seen that the subsets in equation ( 16) are nested to each other such as C ⊂ E ⊂ M ⊂ S ⊂ N. Based on the degree of belief, m is assigned to a consonant set of damage levels, possibility distributions for predetermined PGA region are generated using the assignment of the degree of belief.e sequence of monotone measures of each damage level with respect to PGA is presented and compared with conventional fragility curves generated by maximum likelihood estimation.It is noticeable that the evidence for at least no damage, N, is considered as any damage level.

Case Study
A schematic representation of the proposed framework to generate fragility curves using possibility distributions is shown in Figure 2. At the rst step, the empirical damage data are rearranged in the shape of binomial damage data for damage states.In this step, the damage data are rearranged in the ascending order with respect to PGA.As this procedure is out of scope of this study, the arranged data by other researchers [16] are used for case study.Some part of the arranged data is presented in Table 1.Using the arranged data set, PGA intervals are determined based on the inclusion of damage levels at the second step.At the third step, the degree of belief m is assigned to the sets in equation ( 16) based on the occurrence of damage levels in each PGA interval.At the fourth step, possibility distribution is generated based on the degree of belief.At the nal step, monotone measures for PGA intervals are calculated and the sequence of monotone measures with respect to PGA is determined.
For the case study, the damage data of bridges from the 1994 Northridge earthquake are used to develop empirical fragility curves [16].
e PGA value at the location of bridges is interpolated and extrapolated from the PGA data [16].e binomial damage information of damage states "no damage," "at least minor," "at least moderate," "at least major," and "collapse" in original reference [16] are used to generate possibility-based empirical fragility curves of "no damage fragility curve," "slight damage fragility curve," "moderate damage fragility curve," "extensive damage fragility curve," and "collapse fragility curve," respectively, in this study.
For the comparison, the binomial information for the damage state at PGA is used to generate empirical fragility curves by maximum likelihood estimation with the assumption of lognormal and extreme type I and II distributions.e results of the four families of fragility curves are presented in Figure 3. Advances in Civil Engineering

Results and Discussions
From the binomial information of damage state data, PGA intervals are determined based on the inclusion of damage levels and presented in Table 2.For the last PGA interval, it is divided into two regions as the first collapse evidence at PGA of 0.385 g seems possible, but the next evidence is far from the PGA such as 0.682 g.During the 0.385 g to 0.682 g, as there is no evidence for collapse, the PGA interval is divided for engineering sense.
Based on the occurrence of damage levels in each PGA interval, the degree of belief m is assigned to the sets in equation ( 16) and presented in Table 3.
e possibility distributions for each PGA intervals are generated as shown in Figure 4. e possibility distribution π values for each damage level are calculated based on the degree of belief and presented in Table 4.For example, π for d M for the interval [0.385, 0.680] is calculated as the summation of the degree of belief for N, S, M in the respective interval of Table 3.It follows the definition of π as presented in equation (10).
Using the generated possibility distributions in Figure 4, monotone measures for N, S, M, E, and C can be measured.An example to measure the dual monotone measures in possibility theory, certainty cert(M) and possibility pos(M), for M of the interval [0.323, 0.384] is presented in Figure 5.
e maximum possibility distribution π inside of the interval will be the measure of pos as 100% as defined in equation (11) and the complement of the maximum possibility distribution π outside of the interval will be the measure of cert as 8% as defined in equation (12).e degree of confirmation is calculated using the dual monotone measures as 8% as shown in equation (13).e certainty and possibility measures are presented in Tables 5 and 6, respectively.It is noticeable that the belief and plausibility measures for the sets in equation ( 16) are same with the certainty and possibility measures, respectively, as possibility theory is a special form of evidence theory.e belief and plausibility measures in evidence theory are presented in Tables 7 and 8.
As the possibility measures for all intervals are 100% as presented in Table 6, the degree of confirmation is going to be the same with the certainty measure in Table 5. e empirical fragility curves can be constructed by plotting certainty or the degree of confirmation with respect to the corresponding intervals as shown in Figure 6.
Comparisons with the conventional fragility curves that generated by assuming a distribution are presented in Figures 7-10.It is noticeable that the caption for y-axis in the figures, "fragility" represents "probability" in probability theory and "certainty" in possibility theory.For the comparison of "Slight" damage state as shown in Figure 7, the conventional fragility curves underestimate the failure probability at the PGA interval between 0.323 g and 0.5 g.For the comparison of "Moderate" and "Extensive" damage states as shown in Figures 8 and 9, respectively, the conventional fragility curves underestimate the failure probability at the PGA interval between 0.385 g and 0.5 g.Finally,  Table 1: Some part of the arranged binomial data to generate empirical fragility curves for the case study [16].
No. PGA None ≥Min ≥Mod ≥Maj ≥Col First appearance 1 0.069 Advances in Civil Engineering for the comparison of "Collapse" damage state as shown in Figure 10, the conventional fragility curves underestimate the failure probability at the PGA interval between 0.682 g and 0.8 g.While there is no evidence for damage over the PGA of 0.889 g, the conventional fragility curves for all damage states clearly overestimate the failure probability over the PGA of 0.889 g. ese results are obvious as the conventional fragility curves are generated by assuming the distributions with mathematical formulations such as lognormal and extreme type I and type II distributions.With extensive earthquake damage database, the uniform PGA interval, so called as "PGA bin," can be used instead of the PGA intervals determined by experts in Table 2.In this case, the shape of fragility curve also can be   Advances in Civil Engineering     16) using equation (11) based on the possibility distribution in Figure 4. Table 7: Belief measure of damage subsets in equation ( 16) using equation (3) based on the degree of belief m in Table 3.  16) using equation (4) based on the degree of belief m in Table 3. Advances in Civil Engineering formulated by connecting the necessity measure at the center of bin instead of using step function that has the same value in a PGA interval (bin).An example using "PGA bin" of 0.1 g is presented in Figures 11-14 for the respective damage state of "Slight," "Moderate," "Extensive," and "Collapse".In the gures, the proposed fragility curves using expert's interval are compared with those using the uniform interval.e alternative formulation of the proposed fragility curves by connecting the certainty measure at the center of PGA interval is also presented in Figures 11-14.e comparisons in the gures showed that the uniform interval would be used to generate the proposed fragility curves in this study with a large number of damage data.

Conclusions
A framework to consider the "nonspeci city" of damage data is presented.Considering that the damage states for fragility curves are generally nested to each other, the possibility theory, a special form of the evidence theory for nested intervals, is applied to generate fragility information from seismic damage data.Seismic damage data classi ed by four damage levels were used for a case study.Based on the damage state evidence, the degree of belief is assigned to a consonant set of damage levels.Possibility distributions for predetermined PGA region are generated using the assignment of the degree of belief.e sequence of monotone measures of each damage level  8 Advances in Civil Engineering with respect to PGA is presented and compared with conventional fragility curves generated by maximum likelihood estimation.It was shown that the generated sequence of certainty measure could be used as fragility curves in alternative perspective.Noticeably, the fragility curves from possibility distribution represent the certainty of damage state, and those are generated without any assumption of distributions.e empirical fragility curves using possibility theory were compared with those generated by maximum likelihood estimation.e results showed that the conventional fragility curves generated by assuming the distributions overestimate the failure probability at the relatively high PGA, while those underestimate the failure probability at the medium range of PGA.erefore, the proposed empirical fragility curves using possibility theory can be used as an alternative methodology for earthquake engineering.
with the summation of m as unity as shown in equation ( 2) such as e monotone measures for a subset are calculated by equation ( 3) for belief and equation ( 4) for plausibility as presented in Table 9.
Example 2. Special case of evidence theory (same with probability theory).
For the same universe X in example 2, the degree of belief (evidence) m is assigned to the corresponding power set such as e monotone measures, belief and plausibility, are calculated as presented in Table 10.
As shown in results, belief and plausibility measures have the same values for all subsets, and ignorance measures are   Advances in Civil Engineering going to be zero.erefore, belief and plausibility measures can be represented as one monotone measure such as probability.
Example 3. Special case of evidence theory (same with possibility theory).
As shown in Figure 1, if a body of evidence is nested within each other for a discrete universe {a, b, c}, this body of evidence is called "consonant body of evidence".For a consonant body of evidence, belief and plausibility are equal to necessity and possibility, respectively.
For the same universe X in example 2, when the degree of belief (evidence) m is assigned to the corresponding power set such as e monotone measures, belief and plausibility, are calculated as presented in Table 11.Moreover, possibility distributions can be generated using equation (10)
Empirical damage dataBinomial damage data for damage states Determine PGA intervals Generate possibility distribution for each PGA intervalAssign the degree of belief m for damage state subsets in equation (16) at each PGA interval Generate fragility curves by measuring certainty measure for each PGA interval

Figure 2 :
Figure 2: Procedure to generate empirical fragility curves using possibility distributions.

Figure 3 :
Figure 3: Empirical fragility curves determined by maximum likelihood estimation using different distribution functions.At each figure, curves from top to bottom are fragility curves for damage state subsets of S, M, E, and C, respectively.(a) Lognormal CDF.(b) Extreme type I (Gumbel) CDF.(c) Extreme type II (Frechet) CDF.

Figure 6 :
Figure 6: Empirical fragility curves of damage state subsets determined by measuring certainty of possibility distributions for PGA (g) intervals.

Figure 7 :
Figure 7: Comparison of empirical fragility curves of damage state with Slight.

Figure 11 :
Figure 11: Comparison of the proposed fragility curves for Slight about the di erent interval selections with alternative formulation of the proposed fragility curves by using step function and connecting the necessity measure at the center of PGA interval.

Figure 13 :
Figure 13: Comparison of the proposed fragility curves for Extensive about the different interval selections with alternative formulation of the proposed fragility curves by using step function and connecting the necessity measure at the center of PGA interval.

Figure 14 :
Figure 14: Comparison of the proposed fragility curves for Collapse about the different interval selections with alternative formulation of the proposed fragility curves by using step function and connecting the measure at the center of PGA interval.

Table 2 :
Number of damage information at each damage state in PGA intervals.

Table 3 :
Degree of belief m based on the fraction of each damage state in PGA intervals.

Table 4 :
Assigning π based on the degree of belief m in Table3.

Table 6 :
Possibility measure of damage sets in equation (

Table 8 :
Plausibility measure of damage subsets in equation ( Figure 10: Comparison of empirical fragility curves of damage state with Collapse.
as presented in Table12.Using these possibility distributions in Figure 12: Comparison of the proposed fragility curves for Moderate about the different interval selections with alternative formulation of the proposed fragility curves by using step function and connecting the necessity measure at the center of PGA interval.

Table 12 ,
certainty and possibility measures in possibility theory are calculated using equations (11) and (12) as presented in Table13.As expected, the belief and plausibility measures in Table11are equal to the certainty and possibility measures in Table13respectively.