Compliance-based estimation of seismic collapse risk of an existing reinforced concrete frame building

Seismic evaluation of existing structures is based on determining the damage likely to occur during the lifetime of these structures due to earthquake ground motion excitation. However, there is not a consensus about the acceptable level of seismic damage, the expected lifetime of these structures, and the seismic hazard level(s) to evaluate the structures at. This paper presents a methodology for a parametric calculation of the seismic collapse risk of an existing Reinforced Concrete (RC) frame building based on its seismic code compliance, quantified by a dimensionless metric. This metric, defined as compliance factor, compares the seismic capacity of an existing structure with the seismic demand for a new structure at a predetermined hazard level. The inelastic seismic behavior of four models of the RC frame building of varying compliance was analytically investigated in this study to demonstrate the novel methodology. The four models of the RC building were chosen to represent existing RC frame structures designed and constructed before the introduction of modern seismic code provisions. These four building models were excited by a group of earthquake ground motion excitations using Incremental Dynamic Analysis. The collapse probabilities of the four building models, representing varying values of seismic code compliance, were determined for two different locations corresponding to regions of moderate and high seismic hazard, thus laying the basis for the compliance-based estimation of the seismic collapse risk of existing structures.


Introduction
A large percentage of the existing housing inventory worldwide consists of buildings designed before the adoption and enforcement of the modern seismic design codes. Therefore, most of these buildings do not comply with the current seismic code provisions and are exposed to a higher seismic damage risk compared to a new structure. The seismic evaluation and rehabilitation of these buildings is of utmost importance for the maintenance of public safety as well as for the preservation of the cultural heritage of our communities. However, the seismic evaluation and the prioritization of the seismic retrofitting of a portfolio of existing buildings (Grant et al. 2007;Crowley et al. 2008) necessitate the determination of their seismic vulnerability and seismic risk.
The assessment of the seismic vulnerability, as an essential component of the determination of the seismic risk of existing buildings, can be performed using three different approaches: An approach based on seismic damage observed in past earthquake events (Anagnos et al. 1995), an approach based on analytical simulation of the inelastic seismic response of buildings (Silva et al. 2019) and a third approach, which is defined as hybrid vulnerability assessment and is based on an efficient combination of the advantages of the first two approaches (Kappos et al. 2006).
The lack of adequate seismic damage data from past earthquakes inhibits the sufficiency of the implementation of a pure empirical approach for the assessment of the seismic vulnerability of RC buildings. Hence, a wide range of analytical methodologies have been proposed during the last decades for the determination of fragility curves towards the assessment of the seismic vulnerability of existing RC buildings. Singhal and Kiremidjian (1996) performed nonlinear dynamic analyses of RC frame buildings subjected to artificial ground motion records to propose analytical fragility curves based on the Park-Ang global damage indices (1985). The application of a similar methodology was presented by Masi (2003) for the determination of fragility curves of post-1970 RC buildings designed only for vertical loads, including two additional RC building configurations: RC buildings with masonry infills and Pilotis frames. Zeris et al. (2002Zeris et al. ( , 2006 highlighted the importance of conducting nonlinear dynamic analysis for an accurate estimation of the local inelastic deformation demand of RC frame buildings, designed before the adoption of the modern seismic code provisions. The use of a predetermined, finite number of natural or artificial ground motion records for the derivation of the aforementioned fragility curves leads to a discrete representation of the seismic vulnerability of RC buildings. Nevertheless, the increase of the accuracy of this representation necessitates the use of a large set of ground motion records for the conduction of nonlinear dynamic analysis, which can lead to high computational cost. In light of the aforementioned limitations, the use of the Capacity Spectrum Method (Fajfar 1999) and the Static Pushover Analysis were proposed for a simplified and highly efficient estimation of the seismic vulnerability of RC buildings (FEMA HAZUS 2003;Borzi et al. 2008;Fardis et al. 2012;Milosevic et al. 2020). Pinho et al. (2006Pinho et al. ( , 2013 presented a comprehensive analysis of the displacement dispersion obtained by five different nonlinear static procedures for the seismic assessment of existing frame buildings. These methodologies provide a continuous representation of the seismic vulnerability of existing buildings and can be easily applied by the civil engineering community due to their simplicity (D'Ayala et al. 2015;Beyer et al. 2015). However, these methods do not consider the dynamic effect of earthquake ground motion excitation on the inelastic response and seismic damage of existing buildings.
The use of Incremental Dynamic Analysis-IDA (Vamvatsikos and Cornell 2002) emerged from the necessity for a fine balance between the consideration of ground motion variability and the continuous representation of the seismic vulnerability of RC buildings. A software tool (SPO2IDA) was proposed by Vamvatsikos and Cornell (2006) to simplify the use of IDA and reduce its computational cost. Chaulagain et al. (2016) employed SPO2IDA for the assessment of the seismic vulnerability of RC frame buildings in Nepal. Dolšek (2008Dolšek ( , 2010 developed a software tool based on Matlab (MathWorks 2007) and OpenSees (McKenna et al. 2000), which automates the application of IDA for the seismic performance assessment of RC frames. Zeris and Repapis (2018) and Carvalho et al. (2012) compared the use of Static Pushover Analysis and Dynamic Analysis for the assessment of the seismic performance of existing RC buildings.
A common theme across the aforementioned studies is the significant dependence of the fragility curves on the compliance of the RC structure to the seismic code provisions, which is usually defined on a qualitative basis. However, the use of a normalized metric for the quantification of the compliance of these buildings to seismic code provisions would shed more light in the correlation between the code compliance of these buildings and their seismic vulnerability. The quantification of this correlation would make these curves comparable with other fragility curves developed for similar building configurations subjected to the same seismic code provisions. Furthermore, the quantification of the code compliance of these buildings can lead to the determination of their seismic collapse risk (Elwood et al. 2012;Liel et al. 2017;Duvernay et al. 2018;Žižmond andDolšek 2019, Shahnazaryan andO'Reilly 2021) and can justify the application of seismic retrofitting towards the mitigation of this risk (Galanis et al. 2018). However, the assessment of seismic code compliance of existing buildings on a qualitative basis inhibits the accurate estimation of the implications of their compliance on their seismic collapse risk.
In view of this challenge, the use of ratios that express the capacity of an existing structure over its seismic demand has been proposed by several researchers in the past to provide a quantitative probabilistic framework for the performance-based seismic evaluation of existing structures. Cornell et al. (2002) and Jalayer and Cornell (2002) suggested the use of the ratio of the factored demand over the factored capacity of an existing structure, considering aleatory and epistemic uncertainties. Jalayer et al. (2007) proposed the use of this ratio for the structural element that corresponds to the most critical failure mechanism for the system. A simplified form of this ratio that utilizes a seismic design spectrum for the determination of the seismic demand is defined in the Swiss Seismic Code Provisions (SIA 2017) to quantify the seismic code compliance of existing structures. This ratio α s denotes the strength-based compliance factor of an existing structure, defined as (SIA 2017): where C e represents the seismic strength capacity of an existing structure and C c denotes the seismic strength demand for a new, code-compliant structure designed at the same location. The corresponding displacement-based compliance factor of an existing structure is defined as (SIA 2017): where D e represents the seismic displacement capacity of an existing structure and D c denotes the seismic displacement demand for a new, code-compliant structure designed at the same location. Within this frame, the determination of the strength-and displacementbased compliance of an existing structure can be demonstrated using a base shear coefficient C b vs. displacement S d spectrum in a format similar with the Capacity Spectrum method (Freeman et al. 1975;Freeman 1998) and the N2 method (Fajfar 1999): Figs. 1, 2 and 3 illustrate the determination of the compliance of three types of seismically deficient existing structures, assuming in the general case that the strength of the existing and the corresponding new, code-compliant structures (expressed in terms of base shear coefficient C b = S a /g) are lower than the strength demand given by an elastic design spectrum C c,elastic .
(1) a s = C e C c (2) a D = D e D c Figure 1 shows this determination for a strength-and displacement-deficient existing structure, whose strength capacity C e and displacement capacity D e are both lower than the corresponding demand values C c and D c determined for a new structure at the same location, so that α s < 1 and α D < 1 [Eqs.
(1), (2)]. The aforementioned capacity values C e and D e are obtained from the idealized pushover curve of the equivalent Fig. 1 Strength-based compliance factor α s and displacement-based compliance factor α D for a strengthand displacement-deficient existing structure (α s < 1, α D < 1)

Fig. 2
Strength-based compliance factor α s and displacement-based compliance factor α D for a strengthdeficient existing structure (α s < 1, α D > 1) SDOF system of the existing structure (EC8 2004;Fajfar 1999). The corresponding displacement demand D c is derived by the intersection (Demand Point DP c ) of the inelastic demand spectrum for the new, code-compliant structure with the idealized pushover curve of the new structure with strength C c and displacement capacity D c,ult . Figure 2 elucidates the determination of the compliance of a strength-deficient existing structure (α s < 1) whose displacement capacity is greater than the demand for a new structure (α D > 1). The compliance of a displacement-deficient existing structure (α D < 1) is shown in Fig. 3 for the case of α s = 1, representing cases in which the existing structure has the same strength with the new structure, but lower displacement capacity. Tsiavos et al. (2018) and Bender et al. (2017Bender et al. ( , 2019 proposed the use of compliance factors for the estimation of seismic risk of existing structures using constantperiod and constant-yield-displacement approaches . Nevertheless, the simplified fragility curves used in these studies limited their application spectrum to non-ductile RC buildings, such as the one investigated by Fardis and Negro (2005), Bento et al. (2010) and Belejo et al. (2017), which collapse immediately after the exceedance of a predetermined ground motion intensity measure. In an extension, this study aims to integrate the compliance factor-based metric in a comprehensive approach for the estimation of the seismic collapse risk of existing RC frame buildings.
Along these lines, the goal of this study is to parametrize the probability of collapse of existing buildings using strength-based compliance factors and to derive compliancebased fragility curves that represent the behavior of a typical example of RC frame buildings, constructed before the adoption of modern seismic code provisions. The illustration of this correlation between compliance and seismic risk will lay the basis for the risk-based seismic evaluation of existing structures and the decision making about the seismic retrofit of these structures, based on a targeted seismic risk level that is accepted by our communities.

Methodology
The methodology proposed in this paper includes three steps: 1. Estimation of seismic hazard in the location of the RC building. 2. Determination of the compliance of the RC building to the existing seismic code provisions. 3. Quantification of the seismic risk of the RC building, expressed by its probability of collapse.

Estimation of seismic hazard in the location of the RC building
The first step of the proposed methodology comprises the estimation of the seismic hazard in the location of the building through the use of seismic hazard curves: These curves show the probability that an earthquake intensity level is exceeded during a predetermined time period. The Peak Ground Acceleration (PGA) was chosen as the earthquake intensity measure (IM) in this study to facilitate a comparison of the analytically derived results with data obtained from experimental campaigns conducted using the chosen RC frame building that are based on the excitation of the structure with different PGA levels (Fardis 1996;Dolšek 2008Dolšek , 2010. Furthermore, the choice of PGA as an IM facilitates a comparison of this study with the analytically derived, PGA-based fragility curves developed for RC frame buildings of varying compliance to seismic code provisions, such as the ones determined by Kappos et al. (2006) and the PGA-based fragility curves considering ageing effects on RC frame buildings (Couto et al. 2020). Kostinakis et al. (2015), Cantagallo et al. (2012), Yakut and Yılmaz (2008) and Iervolino et al. (2018), among others, demonstrated the high correlation of the spectral acceleration at the fundamental period of the structure with the maximum and average interstorey drifts of RC buildings, indicating that such period-dependent intensity measure could be more appropriate. However, Cantagallo et al. (2012) suggest the use of the cracked (secant) period of the structure to select a spectral acceleration as an intensity measure: determining such a period for an existing building is not easy in practice, considering the variations from the as-designed state due to construction, use and ageing. Within this frame, the use of a vibration-period-independent IM in this study is a good compromise between simplicity, computational efficiency and accuracy. Two different locations have been chosen to represent the effect of two different seismic hazard levels on the seismic risk of existing RC buildings: Sion, Switzerland and Athens, Greece. The seismic hazard curves for both regions are determined using the seismic hazard platform of the European Facilities for Earthquake Hazard and Risk (EFEHR 2017) and are shown in Fig. 4.

Numerical simulation of the RC building
The second step of the methodology proposed in this study consists of the determination of the strength-based compliance factor α s of a RC building through the numerical simulation of the inelastic behavior of the building. A four-storey RC frame building with masonry infills (Dolšek 2010) was chosen to simulate the design standards for RC frame buildings constructed between 1970 and 1980, thus illustrating the methodology proposed in this study. The geometry and the reinforcement of the building, defined as Model A, are shown in Fig. 5. The slab thickness is 15 cm. The concrete quality used for the construction of the building is C25/30, while the quality of the steel reinforcement is B500. The floor masses of the building are: 87 t (first floor), 86 t (second and third floor) and 83 t (fourth floor). The design base shear of the building is 529 kN (Fardis 1996). The transverse reinforcement of the beams and columns was capacity-designed to promote flexural plastic hinge formation. Stirrups Ø10/100 are used in the critical regions of columns (Negro et al. 1994). The vibration periods of the first three vibration modes of the building and the corresponding mass participation ratios are shown in Table 1.
The numerical model of the building was generated using OpenSees (McKenna et al. 2000) and the PBEE toolbox (Dolšek 2010) and consists of beam and column elements. The floor slabs are assumed to be rigid in their plane.

Moment-rotation relation
The use of the PBEE toolbox (Dolšek 2010) enabled the simulation of the inelastic behavior of the columns and beams of the existing RC frame building using plastic hinges with a trilinear moment-rotation relation and material softening (Dolšek 2008). The plastic hinges, modelled as zero-length elements are introduced in the beginning and the end of each structural element of the model (Amrein 2020;Dolšek 2008), using a uniaxial hysteretic material model (Mazzoni et al. 2006). An elastoplastic strain-stress relationship is used for the steel bars, whereas in the case of the concrete the stress-strain relationship prescribed in EC 2 (2004) for non-linear analysis is assumed for calculation of the moment-curvature relationship. The moment-curvature analysis is performed up until the ultimate deformation of the concrete ε cu = − 3.50‰ or the ultimate deformation of the reinforcing steel ε su = 10%. For the beams, the plastic  Fardis (1996) according to EC8 and investigated by Dolšek (2008Dolšek ( , 2010  hinge was used for major axis bending only. For the columns, two independent plastic hinges for bending about the two principal axes were used.

Pushover curve
The force distribution for the application of the horizontal, incrementally increasing static load to the RC frame building is determined by the multiplication of the fundamental mode shape with the floor masses. The nonlinear static analysis includes the P-Delta effects. The pushover curve for the building presented in Fig. 5 is shown in Fig. 6 and is similar to the one obtained by Dolšek (2008) for the same building. A bilinear approximation of the pushover curve is presented on the same Figure. This bilinear approximation is based on the principle of equal areas under the curves (Fajfar 1999), considered up to the point where the first significant drop in strength (usually about 20%) occurs in the nonlinear, actual pushover curve, as proposed by Kappos et al. (2006) and Kappos and Panagopoulos (2010).

Determination of the strength-based compliance of the RC building
The strength-based compliance factor α s of the RC building [Eq.
(1)] was used to define and parametrize its seismic compliance. The elastic seismic design spectra for a new building located in Athens, Greece (Soil Class A, ZII, α gd = 0.24 g) and a new building located in Sion, Switzerland (Soil Class B, Z3b, α gd = 0.16 g) are used to define the base shear coefficient demand C c and are presented in Fig. 7, using the spectrum format presented in Fig. 1, 2 and 3. The EC8 (2004) elastic design spectrum for α gd = 0.24 g is presented in the same Figure. The use of elastic design spectra was chosen for the definition of the code-compliant new building as a special case of the inelastic design approach shown in Figs. 1, 2 and 3 to facilitate a consistent comparison between the seismic design spectra corresponding to Sion, Athens and EC8 (2004). The base shear coefficient-displacement (C b -S d ) capacity of the building is derived by the transformation of the MDOF system to a SDOF system, prescribed by EC8 (2004), using the values derived from the eigenvalue analysis for the numerically simulated

Quantification of the seismic risk of the RC building
The last step of the proposed methodology comprises the compliance-based determination of the seismic risk of the building, expressed by its probability of collapse.

Incremental dynamic analysis
Incremental Dynamic Analysis (Vamvatsikos and Cornell 2002) is used in this study to determine the maximum interstorey drift of the RC frame building presented in Fig. 5 for varying values of the chosen Intensity Measure (PGA). Each IDA curve is derived from the nonlinear dynamic analysis of the RC building subjected to an incrementally scaled ground motion record and comprises 20 points. These points correspond to 20 increments of the PGA of the scaled record, as performed by Tsiavos and Stojadinovic (2016). The scaling is automated by the PBEE toolbox using a hunt and fill algorithm developed by Vamvatsikos and Cornell (2004). The ground motion ensemble chosen for the conduction of the IDA comprises 12 ground motions, obtained from the PEER ground motion database (2018). The ground motion ensemble is shown in Table 2. The minimum and maximum scaling factors for the motions used in the IDA analysis are shown in the same Table.

Collapse limit state criterion
The interstorey drift ratio of 4% is determined as the collapse limit state criterion in this study, as presented by Kappos et al. (2006). For each ground motion, the collapse intensity is determined as the intensity, for which the maximum interstorey drift exceeds the predetermined interstorey drift. ratio of 4%. The IDA curves corresponding to the 12 ground motions used in this study and a vertical line denoting the selected drift-based collapse criterion are shown in the following Fig. 8.

Fragility curve
The collapse fragility curve of the RC building for each value of the selected intensity measure (PGA) is derived by the probability of exceedance of the drift-based collapse criterion for this intensity, based on the results of the IDA curves (Fig. 8).
Assuming a lognormal distribution, the corresponding parameters can be estimated from the obtained collapse intensities (Baker 2015). The cumulative distribution function and the corresponding parameters (θ, β) are determined using Eqs. 3, 4 and 5: The fragility curve of the building presented in Fig. 5, derived through the aforementioned procedure (Amrein 2020), is shown in Fig. 9.

Probability of collapse
The mean annual rate of collapse can be calculated through the integration of the fragility curve of the building (Fig. 9) over the seismic hazard curve (Fig. 1) using Eq. (6) (Eads et al. 2013): Fig. 8 IDA curves of the RC frame building Fig. 9 Fragility curve of the RC frame building 1 3 The probability of collapse during the service life of the building T L can be then calculated using Eq. (7), assuming that the occurrence of earthquakes in time follows a Poisson distribution (Eads et al. 2013): The derived probability of collapse of the building (Model A) for varying values of its service life and two different locations (Sion and Athens) is shown in Fig. 10. As shown in the Figure, the compliance factor of the building located in Sion α s = 1.40 leads to a probability of collapse of 1% in 50 years. However, the location of the same building in a region of higher seismic hazard, such as Athens, with a compliance α s = 0.83 corresponds to a five times higher probability of collapse: 5% in 50 years. This result is in accordance with the investigation of Kappos et al. (2007), who determined the economic losses attributed to RC frame structures of low code-compliance subjected to the 1999 Athens earthquake.

Parametric investigation
A parametric investigation was performed in this study following the proposed methodology to extend the results obtained for one building to four RC frame buildings of varying compliance.
Thus, the longitudinal reinforcement of the RC frame building shown in Fig. 5 (Model A) was reduced to represent three additional RC frame building classes of lower seismic code compliance, denoted as Model B, Model C and Model D. The longitudinal reinforcement diameters and the strength-based compliance factors α s of the four models according to the Swiss and the Greek seismic code provisions are shown in Table 3. The pushover curves of the buildings for varying vaues of their compliance, defined according to Greek seismic code provisions, are presented in Fig. 11.

IDA curves of the four building models of varying compliance
The IDA curves for the four buildings types of varying seismic code compliance, defined according to the Greek seismic code provisions, are shown in Fig. 12. As presented in the figure, the dynamic instabilities manifested by increased values of drift for small intensity increments are observed for more ground motions when the seismic code compliance decreases, as observed by Tsiavos et al. (2017Tsiavos et al. ( , 2020. These instabilities indicate more extensive structural collapse (Vamvatsikos and Cornell 2002;D'Ayala et al. 2015) for decreased seismic code compliance and can be more clearly illustrated in the fragility curves shown in Fig. 13 for the four selected building models of varying code compliance. 1 3

Fragility curves of the four building models of varying compliance
As shown in Fig. 13, the fragility curves become steeper due to a decrease of the code compliance from 0.83 to 0.67 (Model A to Model B) indicating higher probability of exceedance of the collapse limit state for these buildings at all intensity levels. This increase of the probability of exceedance manifests itself only for intensities (PGA) exceeding 0.25 g when the code compliance is further reduced from 0.67 to 0.49 (Model B to Model C). However, this collapse probability increase due to a reduction of the compliance factor from 0.49 to 0.34 (Model C to Model D) is more uniform and intensity-independent.

Probability of collapse of the four building models of varying compliance
The mean annual rate of collapse of each of the four buildings λ c is calculated through the integration of the fragility curve of the building (Fig. 13) over the seismic hazard curve using Eq. (6) (Eads et al. 2013). The emerging deaggregation of the collapse rate λ c is shown in Fig. 14 for Sion, Switzerland and in Fig. 15 for Athens, Greece. These graphs indicate the contributions of different levels of ground motion intensity to the total collapse risk (Eads et al. 2013).
As presented in both Figures, the ground motion intensity, at which the maximum annual rate of collapse occurs, shifts from the initial value of PGA = 0.35 g to smaller intensities when the code compliance of the building decreases, leading to a PGA = 0.25 g for Model D (α s = 0.41, according to the Greek seismic code provisions).
Along these lines, the low seismic code compliance of RC buildings increases the spectrum of ground motions that can can lead to collapse of these buildings, thus increasing their seismic vulnerability and seismic risk, as shown by Iervolino et al. (2007). This increase of seismic risk for lower values of compliance is expressed through the estimation of the corresponding probabilities of collapse during the service life for Sion and Athens locations shown in Fig. 16 and Fig. 17, respectively. The decrease of the compliance factor from α s = 1.40 to α s = 0.62 leads to an increase of the collapse probability for a 50 year service life from 1 to 3% for Sion. This increase of collapse probability for a decrease of the corresponding compliance factor from α s = 0.83 to α s = 0.34 is much more pronounced for regions of high seismic hazard, such as Athens, where the collapse probability increases from 5 to 15% for a 50 year service life.
The American Seismic Code Provisions (ASCE 7 2016; ASCE 41 2017; Luco et al. 2007) prescribe a risk-based ground motion selection procedure for the design of new structures that corresponds to an accepted level of seismic risk, expressed by a probability Fig. 14 Deaggregation of mean annual collapse rate λ c for Sion, Switzerland of collapse of 1% in 50 years. Interestingly, the investigated RC frame building with compliance α s = 1.40 located in Sion corresponds roughly to this allowable seismic risk level with a probability of collapse of 0.97% in 50 years (Fig. 16). However, locating the same RC frame building in Athens leads to reduced code compliance expressed by strengthbased compliance factor values below 1, which are associated with probabilities of collapse that are significantly higher than the aforementioned value of 1% in 50 years (Fig. 17).
Within this frame, the presented compliance-based, seismic risk assessment of existing structures should not be building-specific, but also location-dependent: Different communities will be led to different seismic assessments and different optimal seismic retrofit  decisions about their existing non-compliant infrastructure, depending on their seismic hazard level and an accepted level of seismic risk. Determining the code compliance that corresponds to the accepted level of seismic risk is of paramount importance in this process: the additional risk posed by existing structures that do not satisfy the compliance associated with the accepted seismic risk level will have to be mitigated either by seismic retrofitting or demolishion and replacement, depending on their remaining service life and a holistic cost-benefit analysis.

Conclusions
This paper presents a methodology for the seismic evaluation of existing structures, based on their seismic code compliance. The normalized metric that is chosen to quantify the code compliance of existing structures is defined as the compliance factor of the structure. This dimensionless metric normalizes the seismic capacity of an existing structure with the seismic demand for a new structure designed at the same location for a predetermined hazard level.
The inelastic seismic behavior of RC frame buildings of varying compliance was analytically investigated in this study to illustrate the application of the proposed methodology for the seismic evaluation of existing structures using the PBEE toolbox (Dolšek 2010). Four RC frame buildings of varying strength were modelled to represent typical examples of RC frame buildings, designed and constructred before the introduction of the modern seismic code provisions. The strength and displacement capacity of these buildings was derived based on static pushover analysis. The strength-based compliance factors of the four buildings were calculated based on bilinear idealizations of the obtained pushover curves and elastic seismic design spectra.
The four RC buildings of varying compliance were subjected to a ground motion excitation ensemble using Incremental Dynamic Analysis (IDA), implemented through the PBEE toolbox. The comparison of the obtained IDA curves with a collapse criterion established for RC frame buildings led to the development of compliance-based fragility curves for the four buildings. The probability of collapse of the four buildings was derived through the integration of the fragility curves over the seismic hazard curves for two locations representing regions of moderate and high seismic hazard: Sion (Switzerland) and Athens (Greece).
The proposed methodology is easy to use and comprises three steps: The estimation of seismic hazard in the location of the structure, the determination of the compliance of the structure to the existing seismic code provisions and the quantification of the seismic risk of the structure, expressed by its probability of collapse.
This study focused on the application of this methodology to an RC frame building of a predetermined geometry and strength for two seismic hazard levels. However, the presented methodology can be applied for the determination of the relation between the code compliance and the seismic risk of a wide range of existing RC, masonry, steel or timber structures of varying geometries and can be performed for more regions, representing different seismic hazard levels. This determination can provide a location-dependent benchmark database that leads to a direct, compliance-based estimation of the seismic risk of existing structures of a known geometry and building typology. The determination of seismic risk targeted-compliance values in this database could provide a valuable tool to local authorities towards a fast, compliance-based seismic assessment of the existing building inventories at a community level.
Along these lines, this methodology can lay the basis for the efficient determination of the seismic risk of existing structures of varying code compliance, expressed by their compliance factors. This seismic risk assessment can yield valuable insight into the decisionmaking process for the prioritization of the seismic retrofitting of existing structures: The optimal seismic retrofitting strategy for the seismic upgrade of our existing infrastructure can be chosen as the one leading to a seismic risk level accepted by our communities. Within this frame, this seismic risk-based assessment and decision-making process can lead to the development of resilient communities.
Funding Open Access funding provided by ETH Zurich. The research presented in this study was funded by ETH Zurich.

Declarations
Conflict of interest The authors of this paper certify that they have no affiliations with or involvement in any organization or entity with any financial interest (such as honoraria; educational grants; participation in speakers' bureaus; membership, employment, consultancies, stock ownership, or other equity interest; and expert testimony or patent-licensing arrangements), or non-financial interest (such as personal or professional relationships, affiliations, knowledge or beliefs) in the subject matter or materials discussed in this manuscript.
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:// creat iveco mmons. org/ licen ses/ by/4. 0/.