The evaluation of the seismic fragility of buildings is one key task of earthquake safety and loss assessment. Many research reports and papers have been published over the past 4 decades that deal with the vulnerability of buildings to ground motion caused by earthquakes in China. We first scrutinized 69 papers and theses studying building damage for earthquakes that occurred in densely populated areas. They represent observations where macroseismic intensities have been determined according to the official Chinese Seismic Intensity Scale. From these many studies we derived the median fragility functions (dependent on intensity) for four damage limit states of the two most widely distributed building types: masonry and reinforced concrete. We also inspected 18 publications that provide analytical fragility functions (dependent on PGA, peak ground acceleration) for the same damage classes and building categories. Thus, a solid fragility database based on both intensity and PGA is established for seismicity-prone areas in mainland China. A comprehensive view of the problems posed by the evaluation of fragility for different building types is given. Based on the newly collected fragility database, we propose a new approach in deriving intensity–PGA relations by using fragility as the bridge, and reasonable intensity–PGA relations are developed. This novel approach may shed light on new thought in decreasing the scatter in traditional intensity–PGA relation development, i.e., by further classifying observed macroseismic intensities and instrumental ground motions based on differences in building seismic resistance capability.
Field surveys after major disastrous earthquakes have shown that poor performance of buildings in earthquake-affected areas is the leading cause of human fatalities and economic losses (Yuan, 2008). The evaluation of seismic fragility for existing building stocks has become a crucial issue due to the frequent occurrence of earthquakes in the last decades (Rota et al., 2010). Building fragility curves, defined as expected probability of exceeding a specific building damage state under given earthquake ground shaking, have been developed for different typologies of buildings. They are required for the estimation of fatalities and monetary losses due to building structural damage. The development of fragility curves can be divided mainly into two approaches: empirical methods and analytical methods. Empirical methods are based on post-earthquake surveys for groups of buildings and considered to be the most reliable source, because they are directly correlated to the actual seismic behavior of buildings (Maio and Tsionis, 2015). Numerous post-earthquake investigations have been conducted for groups of buildings to derive the empirical damage matrices. A damage matrix is a table of predefined damage states and percentages of specific building types at which each damage state is exceeded due to particular macroseismic intensity levels. However, as pointed out by Billah and Alam (2015), empirical investigations are usually limited to particular sites or seismotectonic/geotechnical conditions with abundant seismic hazard and lack generality. Moreover, they usually refer to the macroseismic intensity, which is not an instrumental measure but is based on a subjective evaluation (Maio and Tsionis, 2015). By contrast, analytical methods are based on static and dynamic nonlinear analyses of modeled buildings, which can produce slightly more detailed and relatively more transparent assessment algorithms with direct physical meaning (Calvi et al., 2006). Therefore, analytical methods are conceived to be more reliable than empirical results (Hariri-Ardebili and Saouma, 2016). Nevertheless, variations in the different practices of analytical fragility studies, such as selection of seismic demand inputs, use of analysis techniques, characterization of modeling structures, definition of damage state thresholds and usage of damage indicators by different authorities, can create discrepancies among various analytical results even for exactly the same building typology. In addition, analytical fragility studies for groups of buildings are computationally demanding and often technically difficult to perform.
Despite the limitations of each fragility analysis method, both empirical and analytical fragility curves are essential in conducting seismic risk assessment. However, the application of the existing fragility curves has been considered to be a challenging task, since different approaches and methodologies are spread across scientific journals, conference proceedings, technical reports and software manuals, hindering the creation of an integrated framework that could allow the visualization, acquisition and comparison between all the existing curves (Maio and Tsionis, 2015). In this regard, the first purpose of this study is to describe and examine available fragility curves, specially developed for Chinese buildings from 87 papers and theses using empirical and analytical methods. The median fragility functions from these previous research findings for the main building types in seismicity-prone areas in mainland China are then outlined.
Furthermore, based on the empirical and analytical fragility database collected, the second purpose of this work is to propose a new approach in deriving intensity–PGA (peak ground acceleration) relations by using fragility as the bridge. The main concern behind this attempt is that the intensity–PGA relation is quite essential in seismic hazard assessment, while traditional practices in deriving such a relation are generally region-dependent and have large scatter (Caprio et al., 2015). Traditionally, intensity–PGA relations are developed using instrumental PGA records and empirical intensity observations within the same geographical range. In this work, we try to establish the intensity–PGA relation using fragility as a conversion medium. Formally, this is achieved by the elimination of the fragility values from the fragility–intensity relation and from the fragility–PGA relation. Theoretically, reasonable results should emerge if the building types used in analytic fragility analyses and those investigated in the empirical field surveys are close enough.
This study is organized as follows. In Sect. 1, the necessity of fragility database construction and the pros and cons of the main fragility analysis methods are briefly introduced. In Sect. 2, a literature review of fragility studies in mainland China and related concepts is provided. Section 3 presents the discrete fragility database extracted from reviewed papers and theses. In Sect. 4, median empirical and analytical fragility curves and their scatter are derived for major building types in seismicity-prone areas in mainland China. In Sect. 5, we introduce in detail our new approach in developing intensity–PGA relations by using fragility as a bridge, which is quite comparable with relations developed by traditional practice. In Appendix and Code and data availability, access to supplementary documents mentioned in the text are provided.
As documented in Calvi et al. (2006), the first application of an empirical
method to investigate building fragility at a large geographical scale was
carried out in the early 1970s. In mainland China, since the 1975 Haicheng
Given the importance of building fragility in seismic risk assessment and
loss mitigation, in total we reviewed 87 existing fragility analyses from
papers and theses for the main building typologies in seismicity-prone areas in
mainland China. It is worth noting that, in Ding (2016), a very detailed
collection of empirical fragility data was provided for 112 Sichuan (Chen et al., 2017; Gao et al.,
2010; He et al., 2002; Li et al., 2015, 2013; Sun et al., 2013, 2014; Sun and Zhang, 2012; Ye et al., 2017; Yuan, 2008; Zhang et
al., 2016), Yunnan (He et al., 2016; Ming et al., 2017; Piao, 2013;
Shi et al., 2007; Wang et al., 2005; Yang et al., 2017; Zhou et al., 2007, 2011), Xinjiang (Chang et al., 2012; Ge et al., 2014;
Li et al., 2013; Meng et al., 2014; Song et al., 2001; Wen et al., 2017), Qinghai (Piao, 2013; Qiu and Gao, 2015), Fujian (Bie et
al., 2010; Zhang et al., 2011; Zhou and Wang, 2015) and other
seismic active zones (A, 2013; Chen, 2008; Chen et al., 1999; Cui and Zhai,
2010; Gan, 2009; Guo et al., 2011; Han et al., 2017; He and Kang, 1999; He
and Fu, 2009; He et al., 2017; Hu et al., 2007; Li, 2014; Liu, 1986; Lv et
al., 2017; Ma and Chang, 1999; Meng et al., 2012, 2013; Shi et
al., 2013; Sun and Chen, 2009; Sun, 2016; Wang et al., 2011; Wang, 2007; Wei
et al., 2008; Wu, 2015; Xia, 2009; Yang, 2014; Yin et al., 1990; Yin, 1996;
Zhang and Sun, 2010; Zhang et al., 2017, 2014; Zhou et al.,
2013).
The main outputs of these post-earthquake surveys are empirical
damage probability matrices (DPMs), which can be used to derive the discrete
conditional probability of exceeding predefined damage limit states under
different macroseismic intensity degrees. That is, for the DPMs,
macroseismic intensity degree is usually used as the ground motion
indicator.
As summarized in the Introduction section, the main drawback of empirical method lies in the subjectivity on allocating each building to a damage state and the lack of accuracy in the determination of the macroseismic intensity affecting the region (Maio and Tsionis, 2015). Furthermore, the interdependency between macroseismic intensity and damage as well as the limited or heterogeneous empirical data are commonly identified as the main difficulties to overcome in the calibration process of empirical approaches (Del Gaudio et al., 2015). By contrast, analytical methodologies produce more detailed and transparent algorithms with direct physical meaning that not only allow detailed sensitivity studies to be undertaken, but also allow for the straightforward calibration of the various characteristics of the building stock and seismic hazard (Calvi et al., 2006). Different from the empirical fragility that is directly collected from post-earthquake surveys, the derivation of an analytical fragility curve is often based on nonlinear fine-element analysis. Popular analytical methods include pushover analysis (Freeman, 1998, 2004), the adaptive pushover method (Antoniou and Pinho, 2004) and incremental dynamic analysis (IDA) (Vamvatsikos and Cornell, 2002; Vamvatsikos and Fragiadakis, 2010). Within these approaches, most of the methodologies available in literature lie on two main and distinct procedures: the correlation between acceleration or displacement capacity curves and spectral response curves, such as the well-known Hazus or N2 methods (FEMA, 2003; Fajfar, 2000), and the correlation between capacity curves and acceleration time histories, as proposed in Rossetto and Elnashai (2003).
The major steps in using analytical methods to study building fragility include the selection of seismic demand inputs, the construction of building models, the selection of damage indicator and the determination of damage limit state criteria (Dumova-Jovanoska, 2000). To combine empirical post-earthquake damage statistics from actual building groups with simulated and analytical damage statistics from modeled building types under consideration, we examined quite a few studies deriving analytical fragility curves for masonry and RC buildings in mainland China. The analysis techniques in these studies vary from static pushover analysis or the adaptive pushover method (Cui and Zhai, 2010; Liu, 2017), to dynamic history analysis or incremental dynamic analysis (Liu et al., 2010; J. Liu, 2014; Y. Liu, 2014; Sun, 2016; Wang, 2013; Yang, 2015; Yu et al., 2017; Zeng, 2012; Zheng et al., 2015; Zhu, 2010) as well as analysis based on necessary statistical assumptions (Fang et al., 2011; Gan, 2009; Guo et al., 2011; Hu et al., 2010; Zhang and Sun, 2010).
As predefined, building fragility describes the exceedance probability of a specific damage state given an ensemble of earthquake ground motion levels. To describe the susceptibility of building structure to a certain ground motion level, four damage limit states are used to discriminate between different strengths of ground shaking: slight damage (LS1), moderate damage (LS2), serious damage (LS3) and collapse (LS4). These four limit states divide the building into five structural damage states, namely negligible (D1), slight damage (D2), moderate damage (D3), serious damage (D4) and collapse (D5). The relation between limit states and structural damage states is illustrated by Fig. 1. Hereafter, fragility curves in this study specifically refer to the probability of exceeding four damage limit states (LS1, LS2, LS3, LS4) under different ground motion levels.
Corresponding relation between structural damage states (DS1, D2, D3, DS4, DS5) and limit states (LS1, LS2, LS3, LS4) (modified from Wenliuhan et al., 2015).
Example of major damage state classification methods (modified after Rossetto and Elnashai, 2003).
Standard definitions of building structural damage states have been issued in different countries and areas. In the European Macroseismic Scale 1998 (EMS1998, 1998) proposed by the European Seismological Commission (ESC), five grades of structural damage are defined: negligible to slight damage (Grade 1), moderate damage (Grade 2), substantial to heavy damage (Grade 3), very heavy damage (Grade 4) and destruction (Grade 5). In the Hazus 99 Earthquake Model Technical Manual, developed by the Department of Homeland Security, Federal Emergency Management Agency of the United States (FEMA) in 1999, generally four structural damage classes are used for all building types: slight damage, moderate damage, extensive damage and complete damage. Other damage state classifications like MSK1969 proposed by Medvedev and Sponheuer (1969) and AIJ1995 (Nakamura, 1995) in Japan issued by the Architectural Institute of Japan are summarized in Table 1. In mainland China, the latest standard GB/T 17742-2008 (2008) was issued in 2008 by the China Earthquake Administration (CEA), in which detailed damage to structural and nonstructural components is defined for each damage state (Table 2).
Detailed definition of building damage states in GB/T 17742-2008 (2008), China.
Notes about qualifiers in italics: very few:
In the empirical method, the fragility curve is derived from damage probability
matrices (DPMs) based on post-earthquake field surveys. DPMs give the
proportions of buildings in each structural damage state (D1, D2, D3, D4,
D5), and they can be used to derive the probability of exceeding each damage
limit state
In the analytical method, the fragility curve is derived by Eq. (2), with the
assumption that building response to seismic demand inputs follows the
lognormal distribution:
During the past 4 decades, more than 2000
The distribution of earthquakes that occurred in mainland
China and its neighboring area, for which field surveys were conducted.
Detailed earthquake catalogue can be found in the Supplement, which
is newly compiled based on Ding (2016) and Xu and Gao (2014). The map was created using Generic Mapping Tools (
Divisions of seismic design level for Chinese buildings (modified after Lin et al., 2011).
The seismic resistance level of masonry and RC buildings is further divided into two classes: level A and level B. The assignment of seismic resistance level in this study is mainly based on supplementary information given in each scrutinized paper, including building age, construction material, seismic resistance code at construction time, load-bearing structure, etc. Given the changes in building quality and corresponding code standard over the past 4 decades in China, buildings constructed in different ages, though with the same nominal resistance level of each period, are reassigned with different seismic resistance levels according to the latest standard. The referred-to grouping criteria are given in Table 3 (more building classification details can be found in the Supplement). Generally, “level A” includes buildings with seismic resistance level assigned as pre-code, low or moderate, and “level B” includes buildings assigned as high.
After grouping the empirical and analytical fragilities based on building
type (masonry and RC) and seismic resistance level (A and B) in Sect. 3.1,
the empirical fragility database based on macroseismic intensity (Fig. 3)
and analytical fragility database based on PGA (Fig. 4) for four damage
limit states (LS1, LS2, LS3, LS4) are thus constructed (data can be found in the Supplement). The
The distribution of empirical fragility data from post-earthquake field surveys, depicting the relation between the exceedance probability of each damage limit state (LS1, LS2, LS3, LS4) at given macroseismic intensity levels. The fragility datasets are grouped by building types (masonry and RC) and seismic resistance levels (A and B).
The distribution of analytical fragility data derived from nonlinear analyses, depicting the relation between the exceedance probability of each damage limit state (LS1, LS2, LS3, LS4) at given PGA levels. The fragility datasets are grouped by building types (masonry and RC) and seismic resistance levels (A and B).
To figure out the outliers in the originally collected fragility database,
the box-plot check method was applied. For each building type
(Masonry_A, Masonry_B, RC_A,
RC_B) and in each damage limit state (LS1, LS2, LS3, LS4),
the corresponding series of fragility data was sorted from the lowest to the
highest value. Three quantiles (
Outlier check using the box-plot method for empirical fragility data. Five macroseismic intensity levels are used to classify the original fragility datasets: VI, VII, VIII, IX, X. “A” and “B” represent the pre-code, low and moderate level and the high seismic resistance level, respectively (more classification details are available in the Supplement). LS1, LS2, LS3 and LS4 are the four damage limit states. Outliers are marked by red crosses, and the red line within each box indicates the 50 % quantile fragility value.
Outlier check using the box-plot method for analytical fragility data. Twelve PGA levels are used to group the discrete analytical fragility datasets: 0.1–1.2 g. “A” and “B” represent the pre-, low and moderate level and the high seismic resistance level, respectively (more classification details are available in the Supplement). LS1, LS2, LS3 and LS4 are the four damage limit states. Outliers are marked by red crosses, and the red line within each box indicates the 50 % quantile fragility value.
After removing outliers, details of the remaining fragility dataset (e.g., the number of data points, the median and the standard deviation of these data) for each damage state of each building type are summarized in Appendix Table B1. The change of standard deviation of each fragility series is shown in Figs. A3 and A4 for empirical and analytical data, respectively. It is worth iterating that, as mentioned in the Introduction section, the organization of this study is centered on two focuses. The first one is to construct a comprehensive fragility database for Chinese buildings from 87 papers and theses using empirical and analytical methods, which is one key component of seismic risk assessment. Based on the empirical and analytical fragility database collected, the second focus is to propose a new approach in deriving intensity–PGA relations by using fragility as the bridge. In this regard, a representative fragility curve should be first derived for each damage state of each building type, and we use the median fragility values to derive such a curve.
To derive the representative fragility curve for each damage limit state
(LS1, LS2, LS3, LS4) of each building type (Masonry_A,
Masonry_B, RC_A, RC_B) for
further study (to derive the intensity–PGA relation in Sect. 5), the median
values (50 % quantile) of each fragility series in Figs. 5 and 6 are
used. For consecutive median fragility curve derivation, cumulative normal
distribution is assumed to fit the discrete median empirical fragilities, and
lognormal distributions are assumed to fit the discrete median analytical
fragilities. For each damage limit state of each building type, the
parameters
Median fragility curve and error-bar analysis derived
from empirical fragility datasets, which depicts the relation between
macroseismic intensity and exceedance probability of each damage limit
state (LS1, LS2, LS3, LS4) for masonry and RC building types (note these
median fragility curves are of varying robustness; see Sects. 4 and 5.3
for more details). The circle within each bar represents the median
exceedance probability of each damage limit state; the length of each bar
indicates the value of the corresponding standard deviation. Only intensity
and PGA values with truncated exceedance probability
Fragility curve and error-bar analysis derived from
analytical fragility datasets, which depict the relation between PGAs
(unit: g) and exceedance probability of each damage limit state (LS1, LS2,
LS3, LS4) for masonry and RC building types (note these median fragility
curves are of varying robustness; see Sects. 4 and 5.3 for more
details). The circle within each bar represents the median exceedance
probability of each damage limit state; the length of each bar indicates the
value of the corresponding standard deviation. Only intensity and PGA values
with truncated exceedance probability
The median fragility curves derived from the discrete fragilities for
empirical data and for analytical data are plotted in Figs. 7 and 8,
respectively. To better illustrate the scatter of the originally collected
discrete fragility data, the error analysis is attached with each regressed
median fragility curve. As can be clearly seen from the regressed fragility
curves in Figs. 7 and 8, there are two obvious trends: (1) for the same
building type (masonry or RC), the higher the seismic resistance level
(A
The median fragility curve parameters regressed from empirical and analytical fragility data. The three values in bold font indicate the low quality of the correspondingly regressed fragility curves (see Sect. 4 for more details).
Note: “fort_level” A and B
represent the pre-, low and moderate level and the high seismic resistance level,
respectively; “damage_state” LS1, LS2, LS3 and LS4
represent the four damage limit states: slight, moderate, serious and
collapse, respectively; “
Mathematically, the goodness of fit of the consecutive median fragility
curve from discrete median fragilities can be measured by statistical
indicator
The intensity–PGA relation has an important application in seismic hazard assessment, since the use of macroseismic data can compensate for the lack of ground motion records and thus help in reconstructing the shaking distribution for historical events. Traditionally, intensity–PGA relations are developed using instrumental PGA records and macroseismic intensity observations within the same geographical range (Bilal and Askan, 2014; Caprio et al., 2015; Ding et al., 2014, 2017; Ding, 2016; Ogweno and Cramer, 2017; Worden et al., 2012). These relations are generally region-dependent and have large scatter (Caprio et al., 2015). In this section, we propose a new approach in deriving intensity–PGA relations based on the newly collected empirical and analytical fragility database. For each building type and each damage limit state, an empirical fragility curve (exceedance probability vs. macroseismic intensity) and an analytic fragility curve (exceedance probability vs. PGA) are available, as derived from the median fragilities in Sect. 4. By eliminating the same fragility value, we can derive the corresponding pair of macroseismic intensity and PGA. Thus, for a series of fragility values, we can further regress the corresponding intensity–PGA relation based on the paired intensities and PGAs. Ideally, we would expect the overlap of all these regressed intensity–PGA relations, regardless of the difference in building type, seismic resistance level and damage state.
Compared with this new approach in intensity–PGA relation development, previous practices directly regressed intensity and PGA datasets within the same geographical range, but no further classification of datasets, for example based on building type or damage state as in this study, was conducted. The lack of further classification of PGA and intensity datasets may explain why the previously derived intensity–PGA relations generally have high scatter. The reason is that, although macroseismic intensity is a direct macro indicator of building damage, higher instrumental ground motion (e.g., PGA) does not necessarily mean higher damage to all buildings. Instead, damage is more determined by the seismic resistance capacity of different building types. Thus, further division of intensity and instrumental ground motion records based on affected building types should promisingly help decrease the scatter of the regressed intensity–PGA relation.
Furthermore, the local site effect also contributes to the amplification of instrumental peak ground motions (PGA or SA), when combining intensity and PGA datasets from areas with different geological background together. This in turn increases the scatter of regressed intensity–PGA relation. In this regard, it is worth emphasizing that, in our PGA-related analytical fragility database, the PGA parameter is not the real instrumental records as used in regressing the traditional intensity–PGA relation, but rather the input PGA records used in experimental fragility analysis (pushover analysis, incremental dynamic analysis, dynamic history analysis, etc.). Therefore, the regional dependence (here we mainly refer to site condition), which contributes to the scatter of the traditional PGA–intensity relation, is not a source of uncertainty in our relation.
As a tentative approach, here we derive the relation between intensity and
PGA using median fragility as the bridge for each damage limit state of each
building type. We are deeply aware that uncertainty is inherent in every
single step in both empirical and analytical fragility analysis. However,
the trial of using the median fragility as the bridge to develop the
intensity–PGA relation proposed here, more importantly, aims at providing a
new approach in this regard compared with traditional practice, not to reduce the uncertainties background (due to differences in building
structure, seismic demand inputs, computation methods, etc.) in deriving
empirical and analytical fragility. By using Eq. (3) for PGA–fragility and
intensity–fragility, respectively, and eliminating fragility as a variable, we
find
Intensity–PGA relations grouped by building types. Only
intensity and PGA values with truncated exceedance probability
Intensity–PGA relations grouped by damage limit states.
Only intensity and PGA values with truncated exceedance probability
These intensity–PGA relations are plotted in Fig. 9 (grouped by building types) and Fig. 10 (grouped by damage limit states). Theoretically, higher damage states can occur only for higher intensities or PGA values. For instance, a LS4 damage state at intensity III would not happen, as reflected by the curves in Figs. 9 and 10: LS1 has the lowest PGA or intensity starting point, while LS4 has the highest. Thus, we plot the intensity–PGA curves for fragility values above 1 %. Ideally, we would expect the overlap of all relation curves between intensity and PGA, whether grouped by building type or by damage state. As a matter of fact, for building types Masonry_A and Masonry_B in Fig. 9, the four intensity–PGA curves of the four damage limit states coincide very well. Meanwhile, the discrepancy in intensity–PGA relations of RC_A for damage states LS1, LS2 and LS3 in Fig. 9 is not surprising, given the relatively high scatter in the original analytical fragility datasets of RC_A (as discussed in Sect. 4 and verified by Appendix Figs. A3–A4).
For building types RC_A and RC_B in Fig. 9, it is observed that for the same intensity levels, the corresponding PGA values of damage state LS4 are much higher than those of damage limit states LS1, LS2 and LS3. For a fixed fragility value, this may be due to the underestimation of intensity by the median empirical fragility curve in Fig. 7 or the overestimation of PGA by the median analytical fragility curve in Fig. 8 or a combination of both effects. In this regard, damage data scarcity at higher damage limit states may contribute to the abnormally high PGA values of LS4. When reviewing the fragility data collection process, it is clear that the construction of an empirical fragility database requires the combination of damage statistics from multiple earthquake events that cover a wide range of ground motion levels. Generally, large-magnitude earthquakes occur more infrequently in densely populated areas; thus damage data tend to cluster around the low damage states and ground motion levels. This limits the validation of high-damage states or ground motion levels (Calvi et al., 2006). According to Yuan (2008), those seriously damaged buildings in earthquake-affected areas are mainly masonry buildings. Therefore, the cause of the abnormally high PGA values of damage state LS4 for RC_A and RC_B can be attributed to the relative scarcity of damage data at higher intensity–PGA levels, especially for RC buildings.
As for the building types Masonry_A and Masonry_B in Fig. 9, for the same intensity level, the PGA values revealed by four damage states of Masonry_B are generally higher than those in Masonry_A. This can be more clearly seen from Fig. 10, in which the intensity–PGA relations are grouped by damage limit states and the PGA values revealed by Masonry_B are generally higher than by all the other building types. To better understand this abnormality, we need to refer to the building seismic resistance level assignment process in this study. In fact, compared with Masonry_A, buildings assigned as type Masonry_B generally have much higher seismic resistance capacity. As mentioned in Sect. 3.1, level A refers to buildings with pre-, low and moderate seismic resistance capacity, and level B refers to buildings with high seismic resistance capacity. According to the grouping criteria in Table 3, buildings assigned as Masonry_B mainly refer to those built after 2001 with seismic resistance level VIII and above. This is obviously a very high code standard (more building classification details can be found in the Supplement). Thus, for the same ground motion level, the damage posed on Masonry_B should be much slighter than on Masonry_A. Consequently, the corresponding intensity revealed by Masonry_B should be lower than by Masonry_A.
Currently in mainland China, the macroseismic intensity level in post-earthquake field surveys is determined by damage states of three reference buildings types, namely (1) Type A, wood structure, old soil, stone or brick building; (2) Type B, single-story or multistory brick masonry without seismic resistance; and (3) Type C, single-story or multistory brick masonry sustaining shaking of intensity degree VII. In this study, buildings assigned as Masonry_B mainly refer to those constructed after 2001 with seismic resistance level VIII and above, and their seismic resistance capability is obviously much higher than all three A–B–C building types. Therefore, intensity levels derived from damage to those less fragile Masonry_B buildings tend to be underdetermined. This may help explain why for the same intensity level the corresponding PGA revealed by the intensity–PGA relation of Masonry_B is higher than that of Masonry_A.
Based on the above discussion and the initial analysis in Sect. 4, it can be
summarized that (a) due to the high scatter in the originally collected
fragility database, the intensity–PGA relations derived for LS1, LS2 and LS3 of
building type RC_A are of low robustness (as validated by
the low
According to the analysis in Sect. 5.3, intensity–PGA curves derived for the
four damage limit states of Masonry_A are of the
highest robustness. Therefore, we first focus only on building type
Masonry_A and average its four curves for discrete
intensity values, to derive the corresponding averaged PGA values, as listed
in Table 5. If we match the data points in Table 5 with a linear relation
between intensity and ln(PGA), we find Eq. (5):
The mean PGA values derived from intensity–PGA relations of Masonry_A based on the newly proposed approach (Sect. 5.4).
By integrating the uncertainty in both original empirical and analytical
fragility data of Masonry_A (as shown in Appendix Figs. A3–A4 and Table B1) into the intensity–PGA relation, the averaged standard
deviation
Based on the summarization in Sect. 5.3, if we only remove those obviously unreliable intensity–PGA curves, namely LS1, LS2, LS3 and LS4 of RC_A and LS4 of RC_B, the range of median PGA values corresponding to each intensity degree can be derived from the remaining intensity–PGA relations, as shown in Table 6. For comparison, the recommended PGA range for each intensity degree in the Chinese Seismic Intensity Scale (GB/T 17742-2008, 2008) is listed in Table 7. The PGA values for intensities VI and VII in our results are higher than those in GB/T 17742-2008 (2008), while for intensities VII, IX and X, the PGA values are quite comparable. We also found that the recommended PGA ranges in GB/T 17742-2008 (2008) are indeed the same as those given in GB/T 17742-1980, which was issued in the 1980s around 4 decades ago. At that time, available damage information used to derive the intensity–PGA relation in China was quite scarce. Therefore, damaging earthquakes that occurred in the United States before 1971 were also largely used, which may not be representative of the situation in China today. Thus, one possible explanation for the relatively low PGAs for low intensity levels (VI, VII) in Table 7 (GB/T 17742-1980/2008) is that the buildings in the 1980s were more fragile than buildings today. Since macroseismic intensity is a direct macro indicator of building damage, today buildings generally have better seismic resistance capacity and thus require higher ground motion (PGA) than buildings in the 1980s to be equally damaged.
The PGA ranges derived from more intensity–PGA relations (Sect. 5.5).
The recommended intensity–PGA relations in China (GB/T 17742-1980/2008).
Since the recommended PGA ranges in GB/T 17742-2008 (2008) are not so
representative of the current building status in mainland China, comparisons
with the latest intensity–PGA relation developed in Ding et al. (2017) are
also conducted. Ding et al. (2017) adopted the traditional practice in
regressing the macroseismic intensities and instrumental PGA records within
the same geographical range, by using records for 28
The latest intensity–PGA relation derived by traditional practice for mainland China (Ding et al., 2017).
We established an empirical fragility database by evaluating 69 papers and theses, mostly from the Chinese literature, that document observations of macroseismic intensities reflecting earthquake damage that has occurred in densely populated areas in mainland China over the past 4 decades. These publications provide empirical fragilities dependent on macroseismic intensities for four damage limit states (LS1, LS2, LS3, LS4) of four building types (Masonry_A, Masonry_B, RC_A, RC_B). We also established an analytical fragility database by scrutinizing 18 papers and theses with results on modeling fragilities for the nominally same building types and the same damage states either by response spectral methods or by time-history response analysis. These analytic methods provide fragilities as functions of PGA. From this wealth of data, we derived the median fragility curves for these building types by removing outliers using the box-plot method.
We proposed a new approach by using fragility as the bridge and derived intensity–PGA relations independently for each building type and each damage state. The potential sources of abnormalities in these newly derived intensity–PGA relations were discussed in detail. Ideally the individual intensity–PGA curves should all coincide and allow us to derive an average relation between intensity and PGA. The coincidence is not 100 % perfect and deviations for the cases where they occur were discussed. Given the high-damage data abundance and wide distribution of masonry buildings in mainland China, for studies referring to historic earthquakes and their losses in seismic active regions, e.g., Sichuan and Yunnan, we recommend utilizing the intensity–PGA relation derived from Masonry_A buildings in Eq. (5).
However, for engineering applications, due to the scatter in original fragility datasets and the simplification in using median fragility to derive the intensity–PGA relation in our proposed new approach, the use of the preliminary intensity–PGA relations developed here should be with caution. It is also worth noting that buildings used for empirical intensity determination and for analytical studies do not coincide: a Masonry_A building in a post-event field survey may encompass a wider range than in an analytic study. Therefore, following the novel idea of using fragility as the bridge to develop an intensity–PGA relation in this study, possible extensions in the future can be performed with fragility analysis for more specifically designed building types that are more representative of those widely damaged building types in the field.
In Fig. A1, the comparison of the Chinese Seismic Intensity Scale with other internationally adopted scales is presented. Additionally, the correspondence relation between intensity and PGA–PGV range suggested by the current seismic intensity scale in China (GB/T 17742-2008, 2008) is also graphically presented in Fig. A2. To better illustrate the scatter of the original fragility datasets we collected, standard deviations of each fragility series are also plotted in Fig. A3 (empirical data) and Fig. A4 (analytical data).
Comparison of the Chinese Seismic Intensity Scale with other internationally used seismic scales (Daniell, 2014; after the work of Gorshkov and Shenkareva, 1960; Barosh, 1969; Musson et al., 2010). In this figure, “Liedu-1980/1999” represents the Chinese Seismic Intensity Scale, which has marginal change compared with the current intensity scale GB/T 17742-2008 (2008) used in China.
The suggested correspondence relation between intensity and the PGA–PGV range by the Chinese Seismic Intensity Scale (GB/T 17742-2008, 2008).
Standard deviation of empirical fragility, namely the
exceedance probability of each damage limit state (LS1, LS2, LS3, LS4)
derived based on empirical fragility datasets for each building type
(Masonry_A, Masonry_B, RC_A,
RC_B; detailed values are given in Table B1). Only intensity
and PGA values with truncated exceedance probability
Standard deviation of analytical fragility, namely the
exceedance probability of each damage limit state (LS1, LS2, LS3, LS4)
derived based on analytical fragility datasets for each building type
(Masonry_A, Masonry_B, RC_A,
RC_B; detailed values are given in Table B1). Only intensity
and PGA values with truncated exceedance probability
In Table B1, more statistical details about our newly constructed fragility datasets, including the number of fragility data before and after removing the outliers, median fragility values used in deriving the fragility curve, and the standard deviation of each fragility dataset for each building type and each damage state in Figs. 7 and 8 are listed. Table B2 provides an unofficial English translation of China seismic intensity scale: GB/T 17742-2008 (2008), which is modified after CSIS (2019).
Statistics of fragility database for each damage limit state and each building type.
Continued.
Note: “origin fragility number” refers to the number of original
fragilities collected for each damage limit state of each building type from
previous studies; “fragility number after removing outliers” refers to the
remaining fragilities after removing outliers using the box-plot check method.
Only intensity and PGA values with truncated exceedance probability
Official Chinese Seismic Intensity Scale: GB/T 17742-2008 (2008; modified after CSIS, 2019).
Continued.
Notes about Qualifiers: very few:
The standard deviation in the intensity–PGA relation for each damage limit state of each building type.
The estimation of the uncertainty of the intensity–PGA relation (Eq. 5) is
not a standard procedure like regression analysis. We have fragility as a
function of intensity with an error on the fragility so that fragility is a
random variable. It is also a random variable when derived as function of
More fragility extraction and building classification details are available in the Supplement: Supplementary_building_ classification_ details.pdf. The earthquake catalog in plotting Fig. 2 is in EQ_list_with_field_survey.xlsx. The empirical and analytical fragility data in Figs. 3 and 4 are available in data_Fig3-4.
The supplement related to this article is available online at:
JED proposed the idea to review the fragility literature for buildings in mainland China. DX conducted the review work and proposed the new approach in deriving intensity–PGA relation and wrote the manuscript. FW proposed the methodology of uncertainty transmission from the fragility to intensity–PGA relation. All authors contributed to the revision of the manuscript.
The authors declare that they have no conflict of interest.
The authors thank the Editor Maria Ana Baptista for actively monitoring the whole reviewing process. Furthermore, we thank the reviewer Mustafa Erdik and the other six anonymous reviewers for their constructively critical and helpful comments, which improved this manuscript substantially. We also acknowledge the thorough review of the language copy-editor and the typesetter of the journal
This research has been supported by the China Scholarship Council (CSC) and by the Karlsruhe House of Young Scientists (KHYS) from the Karlsruhe Institute of Technology (KIT). The article processing charges for this open-access publication were covered by a Research Centre of the Helmholtz Association.
This paper was edited by Maria Ana Baptista and reviewed by Mustafa Erdik and one anonymous referee.