Konstantinos Trevlopoulos
email: konstantinos.trevlopoulos@cea.fr

Philippe Guéguen
email: philippe.gueguen@univ-grenoble-alpes.fr

Agnès Helmstetter

Fabrice Cotton

Earthquake risk in reinforced concrete buildings during aftershock sequences based on period elongation and operational earthquake forecasting

Keywords:

des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Introduction

Appropriate decision-making after a major earthquake is important for enabling short-term decisions, such as those regarding building evacuation or repair. Furthermore, seismic crisis management and safety measures can be improved by combining rapid building damage assessment with shortterm forecasting of the seismicity and associated risks [START_REF] Bindi | On-site structure-specific real-time risk assessment: perspectives from the REAKT project[END_REF]. On an urban scale or in the case of critical structures, performance levels related to immediate occupancy or damage grades should be considered in the short-term forecasting. Specifically, a rapid risk assessment must be extended over the entire duration of the aftershock sequence to consider potential damage accumulations owing to the aftershocks [START_REF] Parolai | Bridging the gap between seismology and engineering: towards real-time damage assessment[END_REF]. Operative procedures are then required to combine aftershock occurrence forecasts, models of damage accumulation, and building health characterizations with respect to targeted levels of damage. Earthquake risk forecasts during an aftershock sequence are often based on operational earthquake forecasting (OEF). Hermann et al. [START_REF] Herrmann | Communicating time-varying seismic risk during an earthquake sequence[END_REF] used OEF to forecast the time-varying seismic risk during an earthquake sequence scenario simulating the 1356 Basel earthquake. They estimated the human losses and associated risks in the short-term to support decisions such as evacuation or the interruption of critical operations in reinforced concrete (RC) buildings. Chioccarelli and Iervolino [START_REF] Chioccarelli | Operational earthquake loss forecasting: a retrospective analysis of some recent Italian seismic sequences[END_REF] used OEF to perform a retrospective analysis of loss assessments for aftershock sequences in Italy.

Decision-making in seismic emergencies must consider the increased vulnerability owing to mainshock damage, in addition to the time-variant vulnerability of buildings in view of probable aftershock damage accumulations. Probabilistic frameworks have been developed to assess the contributions of aftershock damage accumulation with respect to the damage caused by the mainshock [START_REF] Li | Performance evaluation and damage assessment of steel frame buildings under main shock-aftershock earthquake sequences[END_REF][START_REF] Burton | Quantifying the reduction in collapse safety of main shock-damaged reinforced concrete frames with infills[END_REF][START_REF] Iervolino | Closed-form aftershock reliability of damage-cumulating elastic-perfectly-plastic systems: aftershock reliability of elastic-perfectly-plastic systems[END_REF]. Furthermore, approaches based on the Markov Chain have been used to assess seismic performance and to develop models for damage accumulation [START_REF] Shokrabadi | Risk-based assessment of aftershock and mainshock-aftershock seismic performance of reinforced concrete frames[END_REF][START_REF] Jalayer | Seismic risk assessment considering cumulative damage due to aftershocks[END_REF][START_REF] Iervolino | Markovian modeling of seismic damage accumulation: Markovian Modeling of seismic damage accumulation[END_REF], in combination with a probabilistic assessment of aftershock occurrence. Most studies report a large contribution of aftershocks in consequence and loss forecasting. The short-term variability of buildings' vulnerabilities, as related to the damage accumulated over the entire seismic sequence, is a key element to be assessed for short-term decision-making.

Monitoring the elongation of the fundamental period of buildings can help assess earthquake damage in buildings whose stiffness gradually degrades before failure [START_REF] Reuland | A model-based data-interpretation framework for post-earthquake building assessment with scarce measurement data[END_REF]. The fundamental period (or frequency) is assumed to be a proxy for the apparent structural stiffness and structural health [START_REF] Grimaz | How could cumulative damage affect the macroseismic assessment?[END_REF][START_REF] Michel | Quantification of fundamental frequency drop for unreinforced masonry buildings from dynamic tests[END_REF][START_REF] Brun | Pseudo-dynamic tests on lowrise shear walls and simplified model based on the structural frequency drift[END_REF][START_REF] Katsanos | Prediction of inelastic response periods of buildings based on intensity measures and analytical model parameters[END_REF][START_REF] Katsanos | Inelastic spectra to predict period elongation of structures under earthquake loading: inelastic spectra for period elongation of structures[END_REF]. For instance, the residual stiffness of masonry buildings from period measurements has been used to study the effect of seismic damage accumulation on a macro-seismic intensity assessment [START_REF] Grimaz | How could cumulative damage affect the macroseismic assessment?[END_REF]. Michel et al. [START_REF] Michel | Quantification of fundamental frequency drop for unreinforced masonry buildings from dynamic tests[END_REF] used laboratory tests to quantify the fundamental frequency shift for unreinforced masonry specimens as a function of the structural drift and the level of damage. Experimental and numerical approaches have found empirical relationships between the frequency shift and damage index for RC buildings [START_REF] Brun | Pseudo-dynamic tests on lowrise shear walls and simplified model based on the structural frequency drift[END_REF]. Katsanos et al. [START_REF] Katsanos | Prediction of inelastic response periods of buildings based on intensity measures and analytical model parameters[END_REF] showed that the transient period elongation during a seismic response did not exceed 1.2 and 1.7 for a designed earthquake and twice the designed earthquake, respectively. Katsanos and Sextos [START_REF] Katsanos | Inelastic spectra to predict period elongation of structures under earthquake loading: inelastic spectra for period elongation of structures[END_REF] also showed the sufficiency (i.e., independence from magnitude and distance) of the structural period elongation in damage prediction using single-degree-of-freedom oscillators. Reuland et al. [START_REF] Reuland | A model-based data-interpretation framework for post-earthquake building assessment with scarce measurement data[END_REF] predicted fragility curves for subsequent earthquakes based on a measured postearthquake structural frequency and a visual inspection. Masi and Vona [START_REF] Masi | Experimental and numerical evaluation of the fundamental period of undamaged and damaged RC framed buildings[END_REF] reported a 40-50% period elongation in laboratory-scaled specimens of RC buildings without masonry infills under strong motion. They also reported a 15-30% period elongation for site-monitored RC buildings with masonry infills subjected to moderate seismic ground motion. Numerical studies indicate that the period corresponding to a collapse threshold is sensitive to p-delta effects [START_REF] Norouzi | The collapse period of degrading SDOF systems considering a broad range of structural parameters[END_REF].

Regarding the aftershock collapse assessment, the maximum inter-story, roof, and residual inter-story drifts appear to be suitable as damage measures [START_REF] Raghunandan | Aftershock collapse vulnerability assessment of reinforced concrete frame structures: aftershock collapse assessment of reinforced concrete structures[END_REF]. Nevertheless, building tagging procedures often assign a red tag to buildings that have sustained heavy damage without collapsing and those that are considered unsafe, unrepairable, or not worth repairing. In this respect, the procedures herein concern the probabilities of damage states up to heavy damage, based on the period of elongation observed before collapse. Clinton et al. [START_REF] Clinton | The observed wander of the natural frequencies in a structure[END_REF] and Astorga et al. [START_REF] Astorga | Nonlinear elasticity observed in buildings during a long sequence of earthquakes[END_REF] showed the cumulative effects of successive earthquakes on the frequency drop in relation to the amplitude of seismic loading in the context of long-term structural monitoring. Furthermore, Gallipoli et al. [START_REF] Gallipoli | Fundamental period elongation of a RC building during the Pollino seismic swarm sequence[END_REF] reported a co-seismic frequency drop in permanently monitored buildings during mainshocks and after-shocks. Post-earthquake experimental surveys for building damage classification have also highlighted residual period elongations. For instance, a permanent frequency (or period) shift was first associated with a level of seismic damage after the 2003 Boumerdes earthquake in Algeria [START_REF] Dunand | Utilisation du bruit de fond pour l'analyse des dommages des bâtiments de Boumerdes suite au séisme du 21 mai 2003[END_REF]. The damage was assessed by visual screening and was based on a redorange-green traffic light classification. The observed post-earthquake period elongation was less than 30% in the case of buildings tagged green, between 10 and 70% in the case of buildings tagged orange, and 50-70% in the case of buildings tagged red. These values were confirmed by Vidal et al. [START_REF] Vidal | Changes in dynamic characteristics of Lorca RC buildings from pre-and post-earthquake ambient vibration data[END_REF] after the 2011 Lorca earthquake: the period elongations were equal to 20%, 43%, and 65% for the EMS98 damage grades 1 (no structural damage), 2 (slight structural damage, i.e., cracks in structural members), and 3-4 (moderate to very heavy structural damage, i.e., structural members with signs of inelastic deformation or failure), respectively.

This study presents an underlying probabilistic framework for the operational forecasting of time-variant damage-state probabilities for assisting decision-making in seismic emergencies. The main elements of the developed framework are operational aftershock forecasting, measurement of the fundamental period of the structure, damage accumulation and time-variant vulnerability modeling based on period elongation, and building tagging based on a traffic light scheme as a decision-making tool. A novel aspect of this study is the use of fragility curves for damage state thresholds in terms of period elongation. Period elongation is a damage proxy and can be tracked and used for damage probability forecasting. Such fragility curves should be evaluated before the mainshock, so that they can be used in decision making during an earthquake crisis. Although building modeling is used here to evaluate fragility curves, it is not used in the decision-making. First, using the epidemic type aftershock sequence model (ETAS) [START_REF] Kagan | Statistical short-term earthquake prediction[END_REF][START_REF] Ogata | Statistical models for earthquake occurrences and residual analysis for point processes[END_REF][START_REF] Kagan | Probabilistic forecasting of earthquakes[END_REF], synthetic sequences are generated for the given mainshock scenarios.

Subsequently, time-variant damage-state probabilities for numerical models of RC buildings are estimated based on the elongation of the period of the building models. The synthetic aftershock sequences are used in Markov chain Monte Carlo simulations [START_REF] Chioccarelli | Part 2 (AMRA -Multiple shock capacity reduction for non-evolutionary structural systems). D41 Fragility Funct. -Impact Repeated Events Var. Intensities Fragility Funtions Build[END_REF][START_REF] Trevlopoulos | Period elongation-based framework for operative assessment of the variation of seismic vulnerability of reinforced concrete buildings during aftershock sequences[END_REF], in order to estimate the time-variant cumulative probabilities of the considered damage states during seismic sequence scenarios. The fundamental period of the structure measured at a given time can be used for two purposes. First, it can be used to determine the initial damage state of the structure, if the period was measured before the mainshock. Second, if the measurement is made during the aftershock sequence, it can be used to update the assessments for the remainder of the aftershock sequence. Finally, a red-orange-green traffic light framework is developed for classifying buildings according to performance level thresholds. The classification is done based on the increment of the time-variant cumulative probabilities of the damage states in a selected time window. In this framework, the time window considered after the mainshock depends on the time required for emergency decisions (e.g., whether to provide structural support or evacuation). The earthquake scenarios, aftershock sequences, and building models are presented in the second section of this article and are based on the test site of Thessaloniki (Greece). The third section presents the underlying probabilistic approach to the time-variant vulnerability and accumulated damage states. Finally, the traffic light framework for decision-making is presented before the conclusion.

Earthquake scenarios, aftershock sequences, and building models

Synthetic aftershock sequences for earthquake scenarios

For each mainshock, multiple synthetic aftershock sequences are generated in order to be used in a subsequent step to estimate the timevariant cumulative probabilities of the considered damage states. The synthetic aftershock sequences in this framework are generated using the ETAS model. This model is based on empirical laws describing the distributions of earthquakes in space, time, and magnitude. It predicts the rate of future earthquakes based on previous seismicity; it can also be used for OEF [START_REF] Marzocchi | The establishment of an operational earthquake forecasting system in Italy[END_REF], and to generate synthetic catalogs. In this study, we use the ETAS model to generate two synthetic catalogs of aftershocks, as caused by two mainshock scenarios and their aftershocks. Scenarios 1 and 2 (Fig. 1a-b) correspond to the Thessaloniki (Greece) case study: the mainshock (depth = 10 km) magnitudes are 6.5 and 7.0, with epicenters located 10 km and 30 km northeast of the city, respectively. The scenarios are selected based on seismic hazard disaggregation results for Thessaloniki (Fig. 1c). The hazard disaggregation was performed using the OpenQuake Engine [START_REF] Pagani | OpenQuake engine: an open hazard (and risk) software for the global earthquake model[END_REF] and the computation input [START_REF] Pitilakis | Seismic Hazard Harmonization in Europe (SHARE): Online Data Resource[END_REF] from the 2013 Euro-Mediterranean Seismic Hazard Model (ESHM13) [START_REF] The | The 2013 European Seismic Hazard Model: key components and results[END_REF]. The ETAS model considers the uncertainty of the strike and the dip of the seismic fault. In Scenarios 1 and 2, these parameter values are sampled by the Monte Carlo simulation as random independent values from normal distributions. Sampling is performed with a standard deviation of 10° and median values equal to 120° and 60° for strike and dip, respectively. Such fault geometries are found in the regions of the scenario epicenters [START_REF] Papazachos | Surface fault traces and fault plane solutions of the May-June 1978 major shocks in the Thessaloniki area, Greece[END_REF][START_REF] Mountrakis | Neotectonic and seismological data concerning major active faults, and the stress regimes of Northern Greece[END_REF]. In the synthetic output catalogs, earthquakes with a magnitude greater than 2.0 are simulated, but only earthquakes with magnitudes of 4.0 or more are included in the damage model. Moreover, the Guttenberg-Richter law is capped at a magnitude of 8 for the generation of synthetic earthquake sequences, as a conservative assumption compatible with the estimated maximum possible magnitude (M 7.1) for the Thessaloniki region in the last 100 years [START_REF] Cha | Time-independent and time-dependent seismic hazard in greece based on seismogenic sources[END_REF][START_REF] Koravos | Assessment of seismic hazard in Greece and surrounding regions using a design earthquake model[END_REF].

In this study, the parameters required by the ETAS model [START_REF] Helmstetter | Comparison of short-term and time-independent earthquake forecast models for Southern California[END_REF] are estimated based on the earthquake catalogs from the seismological station of the Aristotle University of Thessaloniki (the 550BCE-2010 catalog [START_REF] Papazachos | A catalogue of earthquakes in Greece and surrounding area for the period 550BC-1999[END_REF][START_REF] Papazachos | A catalogue of earthquakes in the Mediterranean and surrounding area for the period 1901 -2010[END_REF] and the 1995-2014 catalog [START_REF]Permanent Regional Seismological Network operated by the Aristotle University of Thessaloniki[END_REF]). A sub-catalog, which excludes earthquakes outside a 200 km radius around Thessaloniki, is used to consider events of similar tectonic settings. Fig. 2a gives the cumulative distribution of the magnitudes in the catalogs and sub-catalog.

The ETAS model assumes that the seismicity rate is the sum of a background rate and the rate of the aftershocks triggered by each previous event. This work is concerned with the seismicity following a large mainshock; the background rate is therefore ignored. The seismicity rate φ at a time t following an earthquake of magnitude M located at a distance r is given by Figure 1. a) Location of study site; b) earthquake scenarios in the area of Thessaloniki, Greece; c) seismic hazard disaggregation for the study site using the OpenQuake Engine [START_REF] Pagani | OpenQuake engine: an open hazard (and risk) software for the global earthquake model[END_REF].

𝜑𝜑(𝑡𝑡, 𝑟𝑟, 𝑀𝑀) = 𝐾𝐾 • 10 𝛼𝛼𝛼𝛼 • 1 (𝑡𝑡 + 𝑐𝑐) 𝑝𝑝 • 1 �𝑟𝑟 + 𝐾𝐾(𝑀𝑀)� 1+𝑞𝑞 (1)
here, 𝐿𝐿(𝑀𝑀) is the associated rupture length of the earthquake, and 𝐾𝐾, 𝑐𝑐, 𝑝𝑝, α, and 𝑞𝑞 are the model parameters. The rate of aftershocks increases exponentially with the magnitude M of the triggering event. In time, the seismicity rate follows Omori's law: the rate of aftershocks decreases as a power law as a function of time after the triggering event. In space, the density is constant for distances smaller than the rupture length 𝐿𝐿(𝑀𝑀) = 0.01 • 10 𝛼𝛼/2 km and decreases according to a power law at large distances 𝑟𝑟 ≫ 𝐿𝐿(𝑀𝑀). The magnitude distribution is usually modeled by a Gutenberg-Richter law. The total seismicity rate 𝑅𝑅(𝑡𝑡, 𝑟𝑟) is the sum of the aftershock sequences of each previous earthquake, and is given by: 𝑅𝑅(𝑡𝑡, 𝑟𝑟) = � 𝜑𝜑(𝑡𝑡 -𝑡𝑡 𝑖𝑖 , |𝒓𝒓 -𝒓𝒓 𝒊𝒊 |, 𝑀𝑀 𝑖𝑖 ) 𝑡𝑡 𝑖𝑖 [START_REF] Parolai | Bridging the gap between seismology and engineering: towards real-time damage assessment[END_REF] In the above, |𝒓𝒓 -𝒓𝒓 𝒊𝒊 | is the norm of the vector 𝒓𝒓 -𝒓𝒓 𝒊𝒊 , i.e., the 3D distance between the point corresponding to 𝒓𝒓 (where the seismicity rate is estimated) and the point corresponding to 𝒓𝒓 𝒊𝒊 (location of earthquake 𝑖𝑖). Owing to this cascading effect, the global rate of aftershocks is much larger than the rate of direct aftershocks as modeled by Eq. [START_REF] Bindi | On-site structure-specific real-time risk assessment: perspectives from the REAKT project[END_REF]. It also obeys Omori's law, but with a smaller 𝑝𝑝 exponent, and with much larger fluctuations than the rate of direct aftershocks given by Eq. ( 1) [START_REF] Helmstetter | Subcritical and supercritical regimes in epidemic models of earthquake aftershocks: epidemic models of earthquake aftershocks[END_REF]. The cascade of secondary aftershocks also slightly modifies the spatial distribution of the aftershocks, which decays as a power law with an apparent exponent smaller than 1 + 𝑞𝑞, and induces a slow diffusion of aftershocks as a function of time [START_REF] Helmstetter | Diffusion of epicenters of earthquake aftershocks, Omori's law, and generalized continuous-time random walk models[END_REF].

The model parameters are usually estimated using a maximum likelihood method [START_REF] Ogata | Statistical models for earthquake occurrences and residual analysis for point processes[END_REF]. However, this method often yields biased and unrealistic values for some parameters. For example, the exponent α (describing the increase in aftershock productivity with the mainshock magnitude) can be underestimated if the spatial distribution of aftershocks is not isotropic [START_REF] Hainzl | Impact of earthquake rupture extensions on parameter estimations of point-process models[END_REF]. In this study, we use a different approach. We start by extracting aftershocks from the catalogs, using the solution proposed by Helmstetter et al. [START_REF] Helmstetter | Importance of small earthquakes for stress transfers and earthquake triggering[END_REF]. We then estimate the model parameters of Eq. (1) (𝑝𝑝 = 1.1, 𝑐𝑐 = 0.001 day, 𝛼𝛼 = 0.8, 𝐾𝐾 = 0.008 and 𝑞𝑞 = 2.0) by comparing the synthetic aftershock sequences generated by ETAS and the aftershock sequences extracted from the regional subcatalog, including events within 200 km from the site. Specifically, 𝛼𝛼 is estimated by fitting the aftershock productivity as a function of the mainshock magnitude (Fig. 2b). The parameters 𝐾𝐾, 𝑝𝑝, and 𝑐𝑐 of Omori's law are found by fitting the aftershock decay as a function of the time since the mainshock (Fig. 2c). The total aftershock rate is characterized by an Omori exponent 𝑝𝑝 = 0.95, i.e., smaller than the value 𝑝𝑝 = 1.1 that characterizes the rate of direct aftershocks as defined by Eq. (1). In space, the exponent 𝑞𝑞 of Eq. ( 1) is obtained by analyzing the distribution of distances between the mainshock and the aftershock hypocenters. Finally, for seismic risk analysis, the duration of the simulated aftershock activity should ideally enable a return to a safe situation, i.e., after all cumulative damage to buildings. In our case, the maximum simulation time estimated by the earthquake catalog analysis is 1,000 days.

Building models

A series of two-dimensional low-rise and mid-rise RC building models is used (Fig. 3), and the buildings are selected as generic buildings, i.e., the focus of this study. These models are used to evaluate the fragility curves for the considered damage states, and the evaluation of the fragility curves must be made before the mainshock. The evaluated fragility curves are used in a subsequent step in order to compute the time-variant cumulative probabilities of the damage states. The models are considered representative of three-dimensional structural systems in Greece and are used to estimate losses [START_REF] Kappos | A hybrid method for the vulnerability assessment of R/C and URM buildings[END_REF]. The "HAZUS" [START_REF]HAZUS-MH 2.1 Technical Manual[END_REF] model-naming convention is used, and the model characteristics are provided in Table 1.

These models include (1) low-code and low-rise models of buildings with a structural system of bare (C1L) and regularly infilled (C3L) moment-resisting frames; (2) low-(C2L) and mid-rise (C2M) low-code models with frames coupled with shear walls; and (3) mid-rise frame models (C1M) with low, moderate, and high seismic code design. These low, moderate, and high code levels refer to seismic provisions for earthquake design. The low-code design level refers to seismic design according to the first Greek seismic code in 1959; the moderate design level refers to the supplementary clauses added in 1985; and the high-code level refers to the 2001 Greek seismic code [START_REF]EAK 2000 -Greek Antiseismic Code 2000[END_REF]. Numerical modeling is performed with OpenSees [START_REF] Mazzoni | Open system for earthquake engineering simulation -user manual[END_REF], and details regarding the building modeling approach can be found in [START_REF] Trevlopoulos | Period elongation-based framework for operative assessment of the variation of seismic vulnerability of reinforced concrete buildings during aftershock sequences[END_REF]. Here, the inelastic behavior of the structural elements is modeled using a distributed plasticity model of hysteretic response with stiffness degradation. The bilinear moment-curvature hysteretic models for the columns and beams are calibrated according to fiber section analysis results. In these analyses, reinforcement steel is modeled using the Giuffre-Menegotto-Pinto model [START_REF] Filippou | Effects of bond deterioration on hysteretic behavior of reinforced concrete joints[END_REF] and the Steel02 material in Open-Sees. The strength degradation of confined concrete is modeled according to [START_REF] Mander | Theoretical stress-strain model for confined concrete[END_REF] using the material Concrete02 in OpenSees. The shear-flexural interaction in the shear walls is not modeled. In that regard, modeling the shear-flexure interaction using models such as the SFI_MVLEM element in OpenSees [START_REF] Kolozvari | Modeling of cyclic shear-flexure interaction in reinforced concrete structural walls. I: theory[END_REF] is incompatible with the scheme herein for evaluating the fragility curves of degraded buildings in damage states. If such models were used for the shear-flexure interaction, one would have to use other procedures, such as sequential non-linear time-history analyses (e.g. [START_REF] Shokrabadi | Risk-based assessment of aftershock and mainshock-aftershock seismic performance of reinforced concrete frames[END_REF]), to compute the fragility curves of the buildings in the damage states. Pinching in the moment-curvature loops is used only for shear walls with a factor of 0.3, assuming a strong pinching effect. The infill walls are modeled with double struts [START_REF] Karapetrou | Time-building specific" seismic vulnerability assessment of a hospital RC building using field monitoring data[END_REF] with a compressive strength of 3.2 MPa and an elastic modulus of 3.2 GPa. The inelastic behavior of the struts is modeled using the Hysteretic material in OpenSees. Moreover, the p-delta effects are not simulated. According to the Greek code for interventions in existing buildings [START_REF]Code for interventions[END_REF], the p-delta effects can be ignored if the inter-story drift sensitivity coefficient is less than a prescribed value. In the cases considered in this study, this value is not exceeded, or it is exceeded for an inter-story drift higher than 0.02. This value of the inter-story drift corresponds to the collapse threshold for RC buildings in Greece [START_REF] Kappos | Seismic Vulnerability and Loss Assessment for Buildings in Greece[END_REF]. We observe that the masonry in the C3L model leads to a considerably lower first eigenperiod with respect to C1L (Table 1) and that the collapse mechanism of the C3L model is the same as that in the C1L model, i.e., plastic hinge formation at the bottom and top of the ground floor columns.

A detailed analysis of the contribution of masonry to period elongation or structural performance is beyond the scope of this study, and readers can refer to Refs. [START_REF] Burton | Quantifying the reduction in collapse safety of main shock-damaged reinforced concrete frames with infills[END_REF][START_REF] Ricci | Influence of infill distribution and design typology on seismic performance of low-and mid-rise RC buildings[END_REF] for further information.

Framework for time-variant vulnerability and cumulative damage state probability

In order to estimate the probabilities of the damage states during an aftershock sequence, we consider the time-variant vulnerability of the buildings through Markov chain Monte Carlo simulations. The simulations rely on the synthetic aftershock sequences, on the fragility curves in terms of period elongation, and on measurements of the period elongation. We consider two period elongation thresholds (Δ𝑇𝑇 1 = 40% and Δ𝑇𝑇 1 = 60%) as the most critical values related to the damage state, according to post-earthquake observations [START_REF] Dunand | Utilisation du bruit de fond pour l'analyse des dommages des bâtiments de Boumerdes suite au séisme du 21 mai 2003[END_REF][START_REF] Vidal | Changes in dynamic characteristics of Lorca RC buildings from pre-and post-earthquake ambient vibration data[END_REF]. These values represent moderate and substantial-to-very-heavy damage, respectively, and are used as performance levels for short-term decision-making in earthquake emergencies, e.g., decisions regarding immediate occupancy. At the core of the framework is the estimation of the cumulative probabilities of exceeding the period elongation thresholds at a given time 𝑡𝑡 ( 𝑃𝑃 [40% < Δ𝑇𝑇 1 ≤ 60%](𝑡𝑡), 𝑃𝑃[60% < Δ𝑇𝑇 1 ](𝑡𝑡)). Based on these estimations, the time-variant probabilities of the damage states are computed as probability increments during a predefined time-window. Fig. 4 shows a flowchart of the framework developed in this study, which is summarized below.

Step 1. A probabilistic seismic demand analysis (PSDA) of the undamaged building model is performed using inelastic finite element time-history analyses. Herein, the probabilistic model is arbitrarily formulated using a ground motion bin approach [START_REF] Mackie | Seismic Demands of Performance-Based Design of Bridges[END_REF]. Each bin consists of 15 synthetic ground motions, which are considered to be hazard-consistent, generated with SeismoArtif [START_REF][END_REF]. For each bin, a different magnitude-distance pair is considered, and the ground motions are simulated while considering a physical model (i.e, the specific barrier model) [START_REF] Halldorsson | Calibration of the specific barrier model to earthquakes of different tectonic regions[END_REF]. Subsequently, they are adjusted in the frequency domain so that the bin average is compatible with the median spectrum for the bin's magnitude-distance pair, as given by a ground motion prediction equation (GMPE). The GMPE is a model proposed by Ambraseys et al. [START_REF] Ambraseys | Rrediction of horizontal response spectra in Europe[END_REF], which has been used to assess the seismic vulnerability of RC buildings in Greece [START_REF] Pitilakis | Consideration of aging and SSI effects on seismic vulnerability assessment of RC buildings[END_REF]. Periods in the [0.02 s-2.0 s] range are selected for spectrum-adjustment, i.e., more than twice the longest eigenperiod of the models herein. Moreover, the procedure for adjusting the records leads to ground motions with reasonable variability. In particular, we consider the dispersions of the fragility curves, which reflect the variability of the selected ground motions (among other uncertainties). Specifically, the computed dispersions are similar to those reported by Kappos and Panagopoulos [START_REF] Kappos | Fragility curves for reinforced concrete buildings in Greece[END_REF], who used the building models used herein.

Step 2. The stiffness degradation of the structural elements is modeled as a function of inter-story drift. First, the stiffness degradation ratio is defined as the ratio of the secant modulus (𝐾𝐾 𝑠𝑠𝑠𝑠𝑠𝑠 ) of the momentcurvature loops resulting at each section along the length of the structural elements to the elastic stiffness (𝐾𝐾 0 = 𝐸𝐸𝐸𝐸) of these sections. The stiffness degradation ratio used herein is equivalent to the stiffness reduction factors for concentrated plasticity hinges [START_REF] Marasco | Cascading hazard analysis of a hospital building[END_REF][START_REF] Polese | Damage-dependent vulnerability curves for existing buildings: damage-dependent vulnerability curves for existing buildings[END_REF], which are based on the secant modulus in moment-chord rotation loops. Specifically, the hinge stiffness reduction factors are estimated as a function of inter-story drift, as normalized to the yield chord rotation [START_REF] Ludovico | A proposal for plastic hinges modification factors for damaged RC columns[END_REF]. Fig. 5 shows the numerical results of the stiffness degradation ratio obtained for the bottom section of an outer column on the ground floor of the C1M low-code model. Several assumptions are made: (1) similarly to [START_REF] Carrillo | Damage index based on stiffness degradation of low-rise RC walls: damage index based on stiffness degradation of low-rise RC walls[END_REF][START_REF] Carrillo | Stiffness degradation model of thin and lightly reinforced concrete walls for housing[END_REF], power functions are computed for sections of the model's beam-column finite elements, to estimate the stiffness variations along the structural elements; (2) the location of the sections along the length of the elements is defined by the Gauss-Legendre quadrature rule; (3) prismatic beam elements are used to model the damaged buildings, owing to the fact that the stiffness degradation varies along the length of each element (e.g., more damage is expected to occur at the header and foot of a column than at the middle); (4) the cross-section dimensions remain constant from one section to another for each element, regardless of variations in bending stiffness; and (5) the degraded stiffness in each section is computed with analytical relationships estimated from the moment-curvature results (as given for the C1M low-code model in Fig. 5).

Finite element models of the building with increasing stiffness degradations are built using the power functions for the stiffness degradation ratio. For a given drift, the degraded stiffness values at the sections along each structural element are computed, and the corresponding model is built. One model is built for each drift in a series of drift values, e.g., from 10 -4 to 5•10 -2 . These models are used in modal analyses, to compute the period of each building at increasing levels of stiffness degradation. The elongation (Δ𝑇𝑇 1 ) of the first eigenperiod of the degraded building (𝑇𝑇 1,𝑑𝑑𝑠𝑠𝑑𝑑. ), relative to the period of the original building (𝑇𝑇 1,𝑜𝑜𝑜𝑜𝑖𝑖𝑑𝑑. ), is computed as follows:

Δ𝑇𝑇 1 = 𝑇𝑇 1,𝑑𝑑𝑠𝑠𝑑𝑑. -𝑇𝑇 1,𝑜𝑜𝑜𝑜𝑖𝑖𝑑𝑑. 𝑇𝑇 1,𝑜𝑜𝑜𝑜𝑖𝑖𝑑𝑑. (3)
Fig. 6 shows Δ𝑇𝑇 1 with respect to the maximum inter-story drift, in the case of the C1L and C1M low-code building models configured according to the procedure described. We observe that, for the same inter-story drift, the period of elongation for the four-story C1M low-code model is lower than that for the two-story C1L low-code model. This makes sense based on engineering intuition, as the columns on the ground floor, where the most stiffness degradation is likely to occur, have bigger sections and additional reinforcement in the case of a four-story building. Therefore, for the same drift, we expect the columns of the four-story building to develop less damage. The maximum inter-story drift values corresponding to Δ𝑇𝑇 1 at 20%, 40%, and 60% are then obtained from Fig. 6 and are used as the damage state thresholds. Here, we consider that period elongations of 20%, 40%, and 60% correspond to the thresholds of slight, moderate, and heavy damage, respectively, based on [START_REF] Dunand | Utilisation du bruit de fond pour l'analyse des dommages des bâtiments de Boumerdes suite au séisme du 21 mai 2003[END_REF][START_REF] Vidal | Changes in dynamic characteristics of Lorca RC buildings from pre-and post-earthquake ambient vibration data[END_REF]. By applying the analytical relationships for the stiffness degradation at the element sections (Fig. 5), the stiffness of the element sections is computed for Δ𝑇𝑇 1 at 20% and 40%. The buildings in the considered damage states are then modeled using the degraded stiffness values corresponding to Δ𝑇𝑇 1 at 20% and 40%. Fig. 7 shows the degraded stiffness ratio in the case of the C1M lowcode model with 20% period elongation. For the sake of simplicity, the uncertainty in the stiffness reduction functions is not propagated analytically in the computation of Δ𝑇𝑇 1 as a function of the drift or in the process of estimating the stiffness of the structural elements for each damage state. However, this uncertainty and the damage state threshold uncertainty are considered by adding empirical components to the dispersion of the estimated fragility curves.

Step 3. A PSDA is performed on the degraded building models for Δ𝑇𝑇 1 = 20% and 40%. The seismic response analyses for the PSDA for each degraded model are run successively and in a single-thread mode: on an Intel i7-4710HQ, this takes approximately 30 min, roughly equal to the time required for every undamaged model. The fragility curves are then computed, giving the probability of each degraded building model to exceed the selected period elongation thresholds as a function of an intensity measure (IM). Although other IMs may correlate better with the structural response [START_REF] Luco | Structure-specific scalar intensity measures for nearsource and ordinary earthquake ground motions[END_REF][START_REF] Perrault | Correlation between ground motion and building response using California earthquake records[END_REF], peak ground acceleration (PGA) has been used in operational post-event response systems for structural evaluation [START_REF] Bindi | On-site structure-specific real-time risk assessment: perspectives from the REAKT project[END_REF]. Therefore, the PGA is adopted as the IM in this study. The fragility curves are evaluated with the "Cornell Method" [START_REF] Iervolino | Assessing uncertainty in estimation of seismic response for PBEE[END_REF][START_REF] Zentner | A general framework for the estimation of analytical fragility functions based on multivariate probability distributions[END_REF], i.e., by fitting a linear function of the PGA in log-scale to the inter-story drifts. As an example, Fig. 8 shows the fragility curves corresponding to the probability to exceed the threshold of moderate damage (Δ𝑇𝑇 1 = 40%) for the C1M lowcode undamaged building (Δ𝑇𝑇 1 = 0%), and for the building model with slight damage (Δ𝑇𝑇 1 = 20%). We observe that the probability of moderate damage is higher in the case of the model which has slight damage, as expected. This is generally observed, i.e., the fragilities of the building models for any period elongation threshold increase with increasing existing period elongation.

Step 4. For each seismic event in the mainshock-aftershocks sequence, a random IM value is computed. Specifically, the IM value is a sample from a random variable with a median and standard deviation given by the GMPE for the seismic event. The fragility curves of the damaged and undamaged building models are then used to compute the elements of Markov chain transition matrices [START_REF] Chioccarelli | Part 2 (AMRA -Multiple shock capacity reduction for non-evolutionary structural systems). D41 Fragility Funct. -Impact Repeated Events Var. Intensities Fragility Funtions Build[END_REF][START_REF] Trevlopoulos | Period elongation-based framework for operative assessment of the variation of seismic vulnerability of reinforced concrete buildings during aftershock sequences[END_REF]. By using the Markov chain, we model a building during an aftershock sequence as a system going through a chain of states, with the hypothesis that every subsequent state depends only on the current state, as follows:

𝑇𝑇 (𝑡𝑡𝑖𝑖) = � 𝑃𝑃 11 (𝑡𝑡 𝑖𝑖 ) 0 𝑃𝑃 12 (𝑡𝑡 𝑖𝑖 ) 𝑃𝑃 22 (𝑡𝑡 𝑖𝑖 ) … 0 … 0 … … 𝑃𝑃 1𝑚𝑚 (𝑡𝑡 𝑖𝑖 ) 𝑃𝑃 2𝑚𝑚 (𝑡𝑡 𝑖𝑖 ) … … … 1 � (4)
In the above, 𝑇𝑇 (𝑡𝑡𝑖𝑖) are the transition matrices, and 𝑃𝑃 is the probability of transition from the original state (first index) to the next state (second index). The elements in the last column are computed by:

𝑃𝑃 𝑗𝑗𝑚𝑚 (𝑡𝑡 𝑖𝑖 ) = 𝑃𝑃[𝐷𝐷𝐷𝐷(𝑡𝑡 𝑖𝑖 ) = 𝑚𝑚|𝐷𝐷𝐷𝐷(𝑡𝑡 𝑖𝑖-1 ) = 𝑗𝑗, 𝐸𝐸𝑀𝑀(𝑡𝑡 𝑖𝑖 ) = 𝛼𝛼 𝑖𝑖 ] (5)
here, 𝑃𝑃 𝑗𝑗𝑚𝑚 (𝑡𝑡 𝑖𝑖 ) is the probability of exceeding the threshold of the highest considered damage state (𝐷𝐷𝐷𝐷 = 𝑚𝑚) at the time of a seismic event (𝑡𝑡 𝑖𝑖 ), given that the building is in a lower damage state (𝐷𝐷𝐷𝐷 = 𝑗𝑗, 𝑗𝑗 < 𝑚𝑚) at the previous time-step (𝑡𝑡 𝑖𝑖-1 ). This value is equal to the value of the fragility curve of the building model in damage state 𝑗𝑗 for damage state 𝑚𝑚, considering the intensity measurement (𝛼𝛼 𝑖𝑖 ). In other words, 𝑃𝑃 𝑗𝑗𝑚𝑚 is equal to the probability of transition from damage state 𝑗𝑗 to damage state 𝑚𝑚 for a given 𝛼𝛼 𝑖𝑖 . Subsequently, the elements of each matrix between the diagonal and the last column are computed by:

𝑃𝑃 𝑗𝑗𝑗𝑗 (𝑡𝑡 𝑖𝑖 ) = 𝑃𝑃[𝐷𝐷𝐷𝐷(𝑡𝑡 𝑖𝑖 ) = 𝑘𝑘|𝐷𝐷𝐷𝐷(𝑡𝑡 𝑖𝑖-1 ) = 𝑗𝑗, 𝐸𝐸𝑀𝑀(𝑡𝑡 𝑖𝑖 ) = 𝛼𝛼 𝑖𝑖 ] -𝑃𝑃[𝐷𝐷𝐷𝐷(𝑡𝑡 𝑖𝑖 ) = 𝑘𝑘 + 1|𝐷𝐷𝐷𝐷(𝑡𝑡 𝑖𝑖-1 ) = 𝑗𝑗, 𝐸𝐸𝑀𝑀(𝑡𝑡 𝑖𝑖 ) = 𝛼𝛼 𝑖𝑖 ] (6)
In Eq. ( 6), 𝑃𝑃 𝑗𝑗𝑗𝑗 (𝑡𝑡 𝑖𝑖 ) is the probability of transition from damage state 𝑗𝑗 to damage state 𝑘𝑘, i.e., the probability of exceeding the threshold of damage state 𝑘𝑘 but not the threshold of state 𝑘𝑘 + 1, given that the building is in damage state 𝑗𝑗 (𝑗𝑗 < 𝑘𝑘). Both terms on the right-hand side of Eq. ( 6) are the fragilities for a given intensity measurement (𝛼𝛼 𝑖𝑖 ), i.e., equal to the value of the fragility curves of the building model in damage state j for damage states 𝑘𝑘 and 𝑘𝑘 + 1, respectively. Finally, the diagonal elements are computed by:

𝑃𝑃 𝑗𝑗𝑗𝑗 (𝑡𝑡 𝑖𝑖 ) = 1 -1[𝐷𝐷𝐷𝐷(𝑡𝑡 𝑖𝑖 ) = 𝑗𝑗 + 1|𝐷𝐷𝐷𝐷(𝑡𝑡 𝑖𝑖-1 ) = 𝑗𝑗, 𝐸𝐸𝑀𝑀(𝑡𝑡 𝑖𝑖 ) = 𝛼𝛼 𝑖𝑖 ], 𝑖𝑖 = {1, … , 𝑛𝑛}, 𝑗𝑗 = {1, … , 𝑚𝑚} (7)
here, 𝑃𝑃 𝑗𝑗𝑗𝑗 (𝑡𝑡 𝑖𝑖 ) is the probability of remaining in damage state 𝑗𝑗 and not exceeding the threshold of damage state 𝑗𝑗 + 1 for a given intensity measurement (𝛼𝛼 𝑖𝑖 ), given that the building is in state 𝑗𝑗. As with Eq. ( 6), the second term of the right-hand side corresponds to the fragility. The elements below the diagonal of the transition matrices are assumed to be equal to zero, i.e., it is assumed that transitions to lower damage states are impossible. In other words, it is assumed that once the structural period has been elongated, its reduction is impossible. One transition matrix is computed for each seismic event of the mainshock-aftershock sequence, and one series of matrices is computed for each building model. As the value of the 𝐸𝐸𝑀𝑀 is different for each event, the Markov chain is time-inhomogeneous, i.e., the values of the elements of the transition matrix are time-dependent. The probabilities of the damage states at the times of the seismic events are computed with a Markov chain with discrete time and space elements, as follows:

�𝑃𝑃[𝐷𝐷𝐷𝐷(𝑡𝑡 𝑛𝑛 ) = 1] … 𝑃𝑃[𝐷𝐷𝐷𝐷(𝑡𝑡 𝑛𝑛 ) = 𝑗𝑗] … 𝑃𝑃[𝐷𝐷𝐷𝐷(𝑡𝑡 𝑛𝑛 ) = 𝑚𝑚]� = [1 0 … 0] ∏ 𝑇𝑇(𝑡𝑡 𝑖𝑖 ) 𝑛𝑛 1 , 𝑖𝑖 = {1, … , 𝑛𝑛}, 𝑗𝑗 = {1, … , 𝑚𝑚} (8)
In our case studies, zero values are assigned to all m elements of the initial vector, which is on the right-hand side of Eq. ( 8), except for the first element of the vector, which is equal to 1. It is thus assumed that the building is undamaged before the mainshock. In a real application, it would be possible to define the initial vector based on measurements of the fundamental period of the building, i.e., considering the actual structural health. Moreover, it would be possible to update the assessments at a given time during the aftershock sequence. To do so, in addition to using the measured fundamental period to define the initial vector, the transition matrices for the aftershocks after the given time should be used.

Step 5. The probabilities of the damage states computed with Eq. ( 8) are computed at the times of the seismic events. Therefore, their timestep is variable, and depends on the time of the earthquakes within a sequence. Furthermore, the number and times of the seismic events in each synthetic catalog vary according to the ETAS algorithm. The computed probabilities are re-sampled with a constant time-step of 0.1 day. This is performed to compute the percentiles of the probability of each damage state as a function of time. Figs. 9 and 10 show the 5th, 50th, and 95th percentiles of the computed probabilities 𝑃𝑃 [40% < Δ𝑇𝑇 1 ≤ 60%](𝑡𝑡) and 𝑃𝑃[60% < Δ𝑇𝑇 1 ](𝑡𝑡) in scenarios 1 and 2, respectively, for the C1L and C1M low-code models. The time of the mainshock is 𝑡𝑡 = 1.0 day. As expected, the highest accumulation of 𝑃𝑃[Δ𝑇𝑇 1 > 60%] is during the first day, given that the seismicity rate decreases constantly with time (Fig. 2).

Step 6. The final step of the procedure is described in Section 4. The slopes of the curves in Figs. 9 and 10 depend on the rate of aftershocks and the rate at which the probability of period elongation is accumulated in each Monte Carlo simulation. The increasing trends of 𝑃𝑃 [60% < Δ𝑇𝑇 1 ](𝑡𝑡) (Figs. 9 and 10) are explained by the fact that the buildings may accumulate damage during the aftershock sequence. In addition, the curves become gradually horizontal over the course of time. This is because the rate of aftershocks reduces with time.

There are cases where 𝑃𝑃 [60% < Δ𝑇𝑇 1 ](𝑡𝑡) is higher than 𝑃𝑃[40% < Δ𝑇𝑇 1 ≤ 60%](𝑡𝑡) (e.g., in the case of the C1M low-code model), i.e., the probability of extensive damage is higher than the probability of moderate damage. Moreover, there are cases in which the probability of damage states corresponding to light or moderate damage may decrease over time. Such cases include the 𝑃𝑃[40% < Δ𝑇𝑇 1 ≤ 60%](𝑡𝑡) result for the C1L low-code mode in Scenario 1 (Fig. 9a). Such results may appear confusing, but they are coherent with the fact that damage states are mutually exclusive, i.e., a high probability of one damage state is accompanied by lower probabilities of the other damage states.

The cumulative probability of exceeding Δ𝑇𝑇 1 = 60% also provides information on the probability of observing the damage level (based on the period elongation threshold) at a given time after the mainshock, as compared with the damage expected owing to the mainshock alone. The ratio of the 50 th percentile of 𝑃𝑃[60% < Δ𝑇𝑇 1 ](𝑡𝑡) at the end of the scenarios, divided by the value of the 50 th percentile owing to the mainshock, is provided in Table 2. This ratio varies according to building model, from 1.92 to 3.41 and from 2.10 to 2.34 in Scenarios 1 and 2, respectively. For both scenarios, this suggests that the aftershock sequence contributes significantly to the probabilities of the final damage states.

Moreover, the computation of 𝑃𝑃[60% < Δ𝑇𝑇 1 ](𝑡𝑡) allows for the assessment of the comparative seismic risk of the building models in the selected scenarios. In the context of an operational framework for postearthquake risk management, such rankings are intended to help decision-makers prioritize actions for building classes. The building models are ranked in Table 3 according to their 𝑃𝑃[60% < Δ𝑇𝑇 1 ](𝑡𝑡) results at the ends of Scenarios 1 and 2. A higher rank corresponds to a higher risk relative to the other building models. In both scenarios, the low-rise bare moment-resisting frame model (C1L) and the mid-rise shear wall model (C2M) are at the top and bottom of the ranking, respectively. The only difference between the scenarios is the position of the low-code C3L and C1M models, which switch between the 4 th and 5 th ranks.

Interestingly, in Table 3, the C1M high-code model has a worse rank than the C1M low-code model. Both models are mid-rise models with bare frames; their only difference is the code under which they were designed. The high-code model has a worse performance according to this ranking. However, this should not lead to the conclusion that it is less safe. This result may be explained by the differences between the philosophies of modern and old codes. In Greece and other countries, a behavior factor is used in building design with recent codes. This behavior factor leads to reduced design spectra with respect to the elastic code spectra. This implies that, under the design loads, the buildings will most likely have an inelastic response. In the old code, an inelastic response was undesirable, and the design relied on strength. As an example, it is noted that the yield moment at the foot of the outer columns of the high-code model is approximately 2/3 of the yield moment at the same section in the low-code model. Recent codes aim to prevent collapse under the design loads by designing structural members with sufficient ductility to develop damage and thus absorb energy, thereby avoiding brittle collapse mechanisms.

More ductility means more stiffness degradation and a higher period elongation before collapse. Therefore, the same amount of period elongation appears to correspond to different levels of damage in old and modern code buildings. In that respect, the ranking in Table 3 may be misleading with respect to the high-code model. A better ranking would be possible if the damage state thresholds were calibrated to reflect the same level of damage across all models. However, this is out of the scope of this article, and further investigation is required in this direction.

Time-variant damage state probability and traffic light concept for decision-making

The so-called "traffic light" model is often used for decision-making. It evaluates the safety of structures according to criteria related to their operability and occupancy, assigning each building to one of three categories: the tolerable zone (green), the intermediate zone (orange), and the intolerable zone (red). This strategy has already been applied for post-earthquake crisis management, which consists in visually inspecting buildings and classifying them according to occupancy safety based on an empirical expert appraisal [START_REF] Dunand | Utilisation du bruit de fond pour l'analyse des dommages des bâtiments de Boumerdes suite au séisme du 21 mai 2003[END_REF][START_REF] Vidal | Changes in dynamic characteristics of Lorca RC buildings from pre-and post-earthquake ambient vibration data[END_REF]. It has also been used for continuous and online condition-based solutions for structural aging [START_REF] Guéguen | Consideration of the effects of air temperature on structural health monitoring through traffic light-based decision-making tools[END_REF]. Typically, the "traffic light" defines a "green" level as indicating that the buildings are safe, have no apparent damage, and can be occupied again, an "orange" level as indicating that the buildings present all of the characteristics of damaged buildings, but with a high level of uncertainty, and a "red" level as indicating that occupancy should be suspended immediately and the building demolished. Previous studies [START_REF] Dunand | Utilisation du bruit de fond pour l'analyse des dommages des bâtiments de Boumerdes suite au séisme du 21 mai 2003[END_REF][START_REF] Vidal | Changes in dynamic characteristics of Lorca RC buildings from pre-and post-earthquake ambient vibration data[END_REF] have shown the benefits of experiment-based assessment of the fundamental period of structures in reducing classification uncertainties, particularly for the "orange" level.

In our study, a period elongation-based red/orange/green "traffic light" scheme is proposed to classify buildings according to the estimated time-variant probabilities of damage states. A probability increment between 𝑡𝑡 and 𝑡𝑡 + 7 days [START_REF] Bindi | On-site structure-specific real-time risk assessment: perspectives from the REAKT project[END_REF][START_REF] Iervolino | Closed-form aftershock reliability of damage-cumulating elastic-perfectly-plastic systems: aftershock reliability of elastic-perfectly-plastic systems[END_REF] (𝑃𝑃[40% < Δ𝑇𝑇 1 ≤ 60%] (𝑡𝑡, 𝑡𝑡 + 7), 𝑃𝑃[60% < Δ𝑇𝑇 1 ](𝑡𝑡, 𝑡𝑡 + 7)) is used to tag buildings, based on tolerable damage state probabilities. The 5 th , 50 th , and 95 th percentiles of the computed probability are considered for four tolerable annual rates (𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡 ) of damage state probabilities (Table 4), based on the performance levels selected for seismic risk analysis. Silva et al. [START_REF] Silva | Exploring risk-targeted hazard maps for Europe[END_REF] proposes 5 • 10 -5 and 10 -4 as tolerable annual rates of collapse leading to loss of human life for low and moderate earthquake-prone countries, such as Greece ( 10 -4 ). These values are adopted herein (𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡.1 and 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡.2 ) for low and moderate seismicity regions, respectively, as "human life safety" performance levels.

The thresholds 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡.3 and 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡.4 correspond to the thresholds of the probabilities of structural failure for normal buildings and buildings that are part of critical infrastructure, respectively. Herein, these thresholds correspond to the Eurocode 0 [START_REF]Eurocode -Basis of structural design[END_REF] reliability classes RC2 and RC3. EC0 does not specify tolerable rates of structural failure and exceedance probabilities, but it does indicate an annual probability of failure, as follows:

𝑃𝑃(𝑇𝑇 = 1𝑦𝑦) = Φ(-𝛽𝛽) (9)
In the above, Φ is a cumulative distribution function of the standard normal distribution (mean 𝜇𝜇 = 0 and standard deviation 𝜎𝜎 = 1). The EC0 recommended minimum values of the reliability index (𝛽𝛽) for the RC2 and RC3 structures are 4.7 and 5.2, respectively. Using Eq. ( 9), the tolerable annual rates of failure (𝜆𝜆 𝑓𝑓𝑜𝑜𝑡𝑡,3 and 𝜆𝜆 𝑓𝑓𝑜𝑜𝑡𝑡,4 ) are computed as equal to 10 -6 and 10 -7 , respectively. Finally, the four thresholds of tolerable risk, according to the standard values for earthquake design and performancebased assessment, are computed using a Poisson model as follows:

𝑃𝑃(𝑇𝑇) = 1 -𝑒𝑒 -𝜆𝜆𝜆𝜆 (10)
here, 𝜆𝜆 = 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡 , and a time window 𝑇𝑇 of 475 years is arbitrarily assumed.

This arbitrary selection is similar to the arbitrary selection of a tolerable level for the probability of collapse in 7 days in [START_REF] Iervolino | Closed-form aftershock reliability of damage-cumulating elastic-perfectly-plastic systems: aftershock reliability of elastic-perfectly-plastic systems[END_REF].

In the context of the proposed tagging scheme, tags are given based on the levels of the probability of damage states, with respect to tolerable probability thresholds. If there were no uncertainties in the estimated damage state probabilities, one would simply compare the probability of a damage state with its tolerable threshold. If the probability of the damage state were higher than the tolerable threshold, then the red tag would be assigned. Otherwise, the green would be assigned. However, given the uncertainty in the estimated damage-state probabilities, tagging is considered possible after defining a level of confidence. Here, we arbitrarily select a 95% level of confidence. Therefore, when the 5 th percentile of the probability of a damage state is higher than its tolerable threshold, the red tag is assigned. The green tag is assigned when the 95% percentile of the probability of a damage state is lower than its tolerable threshold. This leads to a period between the red and green tag periods, during which neither the red nor the green tag can be assigned according to this scheme. Throughout the period between the red and green tag periods, the orange tag is assigned, as proposed in [START_REF] Bindi | On-site structure-specific real-time risk assessment: perspectives from the REAKT project[END_REF]. The orange tag marks a period during which we cannot say if the probability of the damage state is higher or lower than its tolerable threshold, based on the defined level of confidence.

At this point, it is worth elaborating on the tagging scheme. A red tag based on 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡1 and 𝑃𝑃[60% < Δ𝑇𝑇 1 ](𝑡𝑡, 𝑡𝑡 + 7) signifies that the accumulated probability of human life loss owing to extensive structural damage between 𝑡𝑡 and 𝑡𝑡 + 7 days is intolerable. This assumption is justified if we consider that 𝑃𝑃 [60% < Δ𝑇𝑇 1 ](𝑡𝑡, 𝑡𝑡 + 7 ) corresponds to the accumulated probability of period elongation corresponding to extensive damage between 𝑡𝑡 and 𝑡𝑡 + 7 days, and that the level of this probability exceeds, with 95% confidence, the corresponding tolerable threshold. Similarly, a red tag based on 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡4 and 𝑃𝑃[40% < Δ𝑇𝑇 1 ≤ 60%] (𝑡𝑡, 𝑡𝑡 + 7) means that the accumulated probability of moderate structural damage between 𝑡𝑡 and 𝑡𝑡 + 7 days is intolerable. In general, all building models are tagged green, regardless of the tolerable damage-state probability thresholds, once a certain amount of time has elapsed since the mainshock; this time depends on the building model and the scenario. This happens because the rate of aftershocks reduces with time, and the accumulation of the probability of damage practically stops at some point.

Current building tagging procedures typically give a tag which remains unchanged. In the procedure herein, a tag may change from red to orange to green. It should be noted that current tagging operations are made with respect to safety during future earthquakes. The proposed tagging procedure assigns tags are based on the tolerable damage state probabilities and the accumulated damage state probability in 7 days. The selection of a 7-day interval is arbitrary, and was adopted by [START_REF] Iervolino | Closed-form aftershock reliability of damage-cumulating elastic-perfectly-plastic systems: aftershock reliability of elastic-perfectly-plastic systems[END_REF] [START_REF] Iervolino | Closed-form aftershock reliability of damage-cumulating elastic-perfectly-plastic systems: aftershock reliability of elastic-perfectly-plastic systems[END_REF], where the accumulated probability of failure in 7 days was computed. This interval can be tuned according to the decision-making context. For example, if decisions should be made based on the cumulative earthquake risk over 3 days, then a 3-day interval should be used. The 7-day interval means that tags are assigned with respect to the safety during a 7-day window during the ongoing aftershock sequence. This also means that even buildings tagged red after the mainshock according to Table 5. Estimated time in days after the mainshock until the red (𝑡𝑡 𝑜𝑜𝑠𝑠𝑑𝑑 ) and orange (𝑡𝑡 𝑜𝑜𝑜𝑜𝑜𝑜𝑛𝑛𝑑𝑑𝑠𝑠 ) tag periods end in Scenario 1 (𝑀𝑀 6.5, 𝑅𝑅 10 km), according to the four tolerable risk thresholds considered in this study (see Table 4). the proposed tagging procedure will eventually be tagged green, when the rate of aftershocks drops to a level at which practically no damage state probability is accumulated. Moreover, a tag from the proposed procedure can coexist with a tag from current building tagging operations. For example, a severely damaged building would be tagged red by current building tagging operations. The same building would be tagged green by the proposed tagging procedure at the end of the aftershock sequence. Figs. 11 and12 show the results of the operational damage state probability forecasting for different building classes and different performance levels. As 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡1 and 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡2 refer to life safety and 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡3 and 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡4 refer to the probability of observing any structural failure, 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡1 and 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡2 are represented for 𝑃𝑃 [60% < Δ𝑇𝑇 1 ]( 𝑡𝑡, 𝑡𝑡 + 7 ), and 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡3 and 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡4 are represented for 𝑃𝑃[40% < Δ𝑇𝑇 1 ≤ 60%](𝑡𝑡, 𝑡𝑡 + 7). The time-variant probabilities of the damage states depend on the scenarios and building models. As expected, the highest probability of the two damage states is observed at the time of the mainshock. Subsequently, a rapid decrease of the 5 th , 50 th , and 95 th percentiles of the computed probability is observed on the first day after the mainshock, followed by a period of slower decrease. The end of the risk period corresponds to the moment when the percentiles of the computed probability drop to very small values (e.g., less than 10 -5 ). The computed duration of this period is also highly variable, and depends on the building typology and the tolerable annual rate of damage state thresholds 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡 . With these figures, decisions can be made based on the structure's health as estimated on the basis of the period elongation, as measured using experimental data.

Fig. 11 shows the computed 𝑃𝑃[40%< Δ𝑇𝑇 1 ≤ 60%](𝑡𝑡, 𝑡𝑡 + 7), corresponding to moderate damage, in the cases of the C1L and C1M lowcode building models in Scenarios 1 and 2. In all cases, this probability exceeds the tolerable risk thresholds after the mainshock, with 95% confidence. These buildings are given a red tag with respect to moderate damage, as long as the 5th percentile of 𝑃𝑃[40% < Δ𝑇𝑇 1 ≤ 60%] (𝑡𝑡, 𝑡𝑡 + 7) exceeds the considered tolerable threshold. Of all of the cases in Fig. 11, the C1L low-code model in Scenario 2 (Fig. 11b) is tagged red the longest with respect to moderate damage. This seems reasonable, as the mainshock in Scenario 2 has a larger magnitude and more aftershocks are expected and given that the C1L lowcode model has fragility curves with higher values [START_REF] Trevlopoulos | Period elongation-based framework for operative assessment of the variation of seismic vulnerability of reinforced concrete buildings during aftershock sequences[END_REF] than the C1M lowcode model. According to Fig. 11b, the C1L low-code building model is tagged red with respect to 𝑃𝑃[40% < Δ𝑇𝑇 1 ≤ 60%](𝑡𝑡, 𝑡𝑡 + 7) for 2.9 days after the mainshock. This means that, over this period, moderate damage is forecasted for this building model, with 95% confidence. Fig. 12 presents 𝑃𝑃[40% < Δ𝑇𝑇 1 ≤ 60%](𝑡𝑡, 𝑡𝑡 + 7) for the same models and scenarios included in Fig. 11. As in Fig. 11, the estimated probabilities for the C1L low-code model are higher than those for the C1M low-code model. Here, we focus on the duration of the orange tag, whose end signifies that the aftershock sequence is no longer a safety concern. Moreover, the duration of the orange tag reflects the amount of uncertainties considered. For both building models in Fig. 12, the duration of the orange tag is longer in Scenario 1. This may seem counterintuitive, as in Scenario 1, the mainshock magnitude is smaller, and less aftershocks are expected. However, the mainshock in Scenario 1 is closer to the building, leading to higher uncertainty with respect to aftershock distances, the intensities of the aftershock ground motions, the estimated 𝑃𝑃[40% < Δ𝑇𝑇 1 ≤ 60%](𝑡𝑡, 𝑡𝑡 + 7), and, ultimately, the duration of the orange tag.

Finally, Tables 5 and6 summarize the time after the mainshock for the buildings to switch from red to orange (𝑡𝑡 𝑜𝑜𝑠𝑠𝑑𝑑 ) and from orange to green (𝑡𝑡 𝑜𝑜𝑜𝑜𝑜𝑜𝑛𝑛𝑑𝑑𝑠𝑠 ) for the four tolerable risk thresholds (λtoli, i = 1-4) for Scenarios 1 and 2, respectively. With respect to 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡1 , only the C1L and C2L models are tagged red for 0.1 days after the mainshock in Scenario 1. In Scenario 2, the C2M model is tagged green after the mainshock based on 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡1 owing to the lower vulnerability of this model. For performance levels corresponding to 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡3 and 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡4 , all building models are tagged red immediately after the mainshock, i.e., the probability of moderate structural damage is intolerable at the beginning of Scenarios 1 and 2. In Scenario 2, for performance levels corresponding to 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡3 and 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡4 , the buildings remain red-tagged for longer after the mainshock, with the length of time increasing with fragility (e.g., 2.9 and 4.4 days for the C1L and C2L models for 𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡4 , respectively). The most striking observation from the results is the length of the orange tag period in most cases for this tolerable risk threshold: for Scenarios 1 and 2, the orange tag period regarding moderate structural damage may exceed one year in certain cases.

Conclusion

In this study, we presented an operational framework to assess damage risk and its variation over time during a seismic crisis. Our objective was to develop an underlying probabilistic framework to be used as an emergency decision support tool with respect to the integrity of structures damaged during a mainshock-aftershocks sequence. The originality of this underlying framework is in the use of fragility curves for damage state thresholds in terms of period elongation as a proxy for damage. This allows the structural period measurements during seismic sequences [START_REF] Dunand | Utilisation du bruit de fond pour l'analyse des dommages des bâtiments de Boumerdes suite au séisme du 21 mai 2003[END_REF][START_REF] Vidal | Changes in dynamic characteristics of Lorca RC buildings from pre-and post-earthquake ambient vibration data[END_REF] to be used to define the initial state vector in the Markov chain computation and to update the forecasts.

The underlying probabilistic framework for the operational forecast of damage states is based on aftershock sequences for two mainshock scenarios in the vicinity of Thessaloniki. The operational earthquake forecasting during the aftershock sequences is provided by an ETAS model. The risk assessment is then completed by associating the spatiotemporal variation of the sequence's events with the probable structural damage and levels of performance considered in this study with respect to moderate damage or human life safety. This study confirms the major contribution of the aftershocks to the risk associated with each scenario. In Scenario 2, the median probability of Δ𝑇𝑇 1 = 60% at the end of the aftershock sequence is more than twice the median probability for the mainshock for all considered building models. In Scenario 1, this ratio is equal to 3.41 in the case of the C2M low-code building model. An original framework is also proposed for assessing risk for decision-making processes during seismic crises based on a physical parameter for characterizing the level of damage in the structures, rather than an empirical visual screening approach. It is worth noting that there is room for further reduction of uncertainty in the results, i.e., mainly in the duration of the orange classification. For example, in the case of the C1L low-code building model, the duration of the orange tag is 17.5 days Table 6. Estimated time in days after the mainshock until the red (𝑡𝑡 𝑜𝑜𝑠𝑠𝑑𝑑 ) and orange (𝑡𝑡 𝑜𝑜𝑜𝑜𝑜𝑜𝑛𝑛𝑑𝑑𝑠𝑠 ) tag periods end in Scenario 2 (𝑀𝑀 7.0, 𝑅𝑅 30 km), according to the four tolerable risk thresholds considered in this study (see Table 4). with respect to the heavy damage state based on the tolerable risk threshold 2. This uncertainty can be reduced by integrating more information on the seismic response of the structure, for example, by using higher modes or damping, by improving the levels of damage proportional to the period elongation, and particularly by selecting intensity measures as efficient and sufficient as possible in predicting the structure's response. In particular, spectral values (e.g., spectral velocity) for integrating the period elongation would significantly decrease the uncertainties [START_REF] Perrault | Correlation between ground motion and building response using California earthquake records[END_REF].

In principle, the proposed framework is empirical. It uses measurements of period elongation before and after seismic events, rather than building modeling. The use of fragility curves by building type is common in decision-making during an earthquake crisis. In this context, generic building models were chosen to evaluate the fragility curves in terms of period elongation. Based on the measured period elongations at a given time, updated assessments for the remaining aftershock sequence can be performed. In that case, the assessments should use an initial state vector defined based on the measured elongation and the events in the synthetic aftershock catalogs after the time of the determination of the period elongation. The intent is to identify buildings that are considered unsafe once they exceed tolerable risk thresholds owing to all seismic events in a selected time window. The developed framework is different and complementary to existing building tagging operations, as discussed above.

The need for measurements is a sine qua non condition. This implies obtaining the periods of buildings before damage, or the more expensive option of using permanent instrumentation, which could be installed in certain critical structures. The continuous monitoring of the fundamental period during a seismic sequence can help assess the risk associated with vulnerability changes and the spatiotemporal distribution of the aftershocks. Based on the variations in the vulnerability and estimated probability of damage, characterizing buildings with respect to performance levels is possible (e.g., immediate occupancy, collapse prevention). In this case, a real-time online algorithm could be developed using a solution like that proposed by Guéguen and Tiganescu [START_REF] Guéguen | Consideration of the effects of air temperature on structural health monitoring through traffic light-based decision-making tools[END_REF] based on the traffic light concept. False alarms would have to be considered by integrating natural period variations and the slow recovery of the elastic properties after a seismic event, as reported by Guéguen et al. [START_REF] Guéguen | Nonlinear dynamics induced in a structure by seismic and environmental loading[END_REF] and Astorga et al. [START_REF] Astorga | Nonlinear elasticity observed in buildings during a long sequence of earthquakes[END_REF]. Such elements could be incorporated into the developed framework and used to support decision-making with respect to recovery time.

Figure 2 .

 2 Figure 2. a) Cumulative magnitude distribution in the two historical earthquake catalogs for Greece (550 BCE-2010, 1995-2014) and in the catalog including earthquakes up to 200 km from Thessaloniki; b) number of aftershocks as a function of the magnitude of the triggering earthquake in the catalog of Greece (1995-2014 earthquakes with 𝑀𝑀 > 2 and distance less than 200 km from Thessaloniki). The solid line is a linear fit with an exponent of 𝛼𝛼 = 0.8; c) aftershock decay as a function of time after mainshock magnitude with 0.5 intervals (color scale). The straight lines represent Omori's law with an exponent of 𝑝𝑝 = 0.95.

Figure 3 .

 3 Figure 3. a) 2D building models representing reinforced concrete buildings in Greece (for details, see Table 1); b) diagram of the finite element model for the C1L lowcode building.

Figure 4 .

 4 Figure 4. Flowchart of the proposed framework. Steps 1-4 are executed "offline", i.e. before the mainshock, and steps 5-6 are executed immediately after the mainshock.

Figure 5 .

 5 Figure 5. Example of the degradation of stiffness using the moment-curvature hysteretic loops (left insert) at the cross-section at the bottom of the outer columns of the C1M low-code model for each input motion represented by crosses. The analytical equation is obtained by considering a bi-linear degradation model (continuous line) to fit the data.

Figure 7 .

 7 Figure 7. Example of the ratio of the stiffness of the structural elements with varying cross-section properties to the stiffness of the corresponding elements in the initial C1M low-code building model (K/K0) for Δ𝑇𝑇 1 equal to 20%.

Figure 6 .

 6 Figure 6. Period elongation (Δ𝑇𝑇 1 ) of the building models as a function of maximum inter-story drift.

Figure 8 .

 8 Figure 8. Examples of fragility curves for the C1M low-code model in terms of period elongation according to[START_REF] Ogata | Statistical models for earthquake occurrences and residual analysis for point processes[END_REF]. Probability of period elongation corresponding to Δ𝑇𝑇 1 > 40% at 𝑡𝑡 𝑖𝑖 given that Δ𝑇𝑇 1 = 0% (for the undamaged building model) or 20% (for a damaged building model) at 𝑡𝑡 𝑖𝑖-1 .

Figure 9 .

 9 Figure 9. Examples of cumulative probability of period elongation corresponding to 40% < Δ𝑇𝑇 1 ≤ 60% after the time of the mainshock (𝑡𝑡 = 1) during synthetic aftershock sequences in the case of the C1L low-code model in a) Scenario 1 and b) Scenario 2, and the C1M low-code model in c) Scenario 1 and d) Scenario 2.

Figure 10 .

 10 Figure 10. Examples of cumulative probability of period elongation corresponding to 60% < Δ𝑇𝑇 1 after the time of the mainshock (𝑡𝑡 = 1) during synthetic aftershock sequences in the case of the C1L low-code model in a) Scenario 1 and b) Scenario 2, and the C1M low-code model in c) Scenario 1 and d) Scenario 2.

Figure 11 .

 11 Figure 11. Examples of probability of period elongation corresponding to 40% < Δ𝑇𝑇 1 ≤ 60% between t and 𝑡𝑡 + 7 days during synthetic aftershock sequences in the case of the C1L low-code model in a) Scenario 1 and b) Scenario 2 and the C1M low-code model in c) Scenario 1 and d) Scenario 2.

Figure 12 .

 12 Figure 12. Examples of probability of period elongation corresponding to 60% < Δ𝑇𝑇 1 between t and 𝑡𝑡 + 7 days during synthetic aftershock sequences in the case of the C1L low-code model in a) Scenario 1 and b) Scenario 2 and the C2L low-code model in c) Scenario 1 and d) Scenario 2.

Table 1

 1 Description of the building models.

	Type Structural system	Height	Design level	T 1,orig. (s)
	C1L	Moment resisting frames	Low-rise Low code	0.48
	C2L	Shear walls equivalent	Low-rise Low code	0.18
	C3L	Regularly infilled frames	Low-rise Low code	0.20
	C1M Moment resisting frames	Mid-rise	Low code	0.61
	C1M Moment resisting frames	Mid-rise	Moderate code	0.62
	C1M Moment resisting frames	Mid-rise	High code	0.64
	C2M Shear walls equivalent	Mid-rise	Low code	0.38

Table 3

 3 Comparative risk ranking of the building models in Scenarios 1 and 2.

	Type	Design level	Rank
			Scenario 1	Scenario 2
	C1L	Low code	1	1
	C2L	Low code	2	2
	C3L	Low code	5	4
	C1M	Low code	4	5
	C1M	Moderate code	6	6
	C1M	High code	3	3
	C2M	Low code	7	7

Table 2

 2 Ratio of the 50th percentile of the cumulative probability of period elongation 60% (𝑃𝑃[60% < 𝛥𝛥𝑇𝑇 1 ](𝑡𝑡)) at the end of the aftershock sequences to the 50th percentile of 𝑃𝑃[60% < 𝛥𝛥𝑇𝑇 1 ] owing to the mainshock in Scenarios 1 and 2.

	Type	Design level	Scenario 1
			50 th percentile of 𝑃𝑃[60% < Δ𝑇𝑇 1 ](𝑡𝑡)	Ratio	50 th percentile of 𝑃𝑃[60% < Δ𝑇𝑇 1 ](𝑡𝑡)	Ratio
			tmain	tend		tmain	tend
	C1L	Low code	0.3582	0.6912	1.93	0.0854	0.1796	2.10
	C2L	Low code	0.3002	0.5751	1.92	0.0836	0.1754	2.10
	C3L	Low code	0.0356	0.0954	2.68	0.0034	0.0076	2.21
	C1M	Low code	0.0384	0.1059	2.76	0.0028	0.0064	2.26
	C1M	Moderate code	0.0250	0.0748	2.99	0.0015	0.0034	2.29
	C1M	High code	0.1367	0.3226	2.36	0.0173	0.0377	2.18
	C2M	Low code	0.0152	0.0520	3.41	0.0011	0.0026	2.34

Table 4

 4 Tolerable risk thresholds used in the building tagging framework.

	Threshold	Tolerable annual rate	Tolerable 𝑃𝑃(𝑇𝑇 = 475 y)	Element at risk
	1	𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡.1 = 1 • 10 -4	4.6 • 10 -2	Structure and human life (moderate seismicity region)
	2	𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡.2 = 5 • 10 -5	2.4 • 10 -2	Structure and human life (low seismicity region)
	3	𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡.3 = 10 -6	4.7 • 10 -4	Normal structure (failure with moderate damage)
	4	𝜆𝜆 𝑡𝑡𝑜𝑜𝑡𝑡.3 = 10 -7	4.7 • 10 -5	Critical structure (failure with moderate damage)

Acknowledgements

The authors thank Prof. Kyriazis Pitilakis of the Aristotle University of Thessaloniki for providing the details of the building models. Part of this work was supported by the harmonized approach to stress tests for critical infrastructures against natural hazards (STREST) project, funded by the EU Seventh Framework programme (Grant agreement No. 603389), and the Seismology and Earthquake Engineering Research Infrastructure Alliance for Europe (SERA) project, funded by the EU Horizon 2020 programme under Grant Agreement Number 730900.