Energy-Based Seismic Risk Evaluation of Tall Reinforced Concrete Building in Vancouver, BC, Canada, under Mw9 Megathrust Subduction Earthquakes and Aftershocks

This article presents a seismic performance evaluation framework for reinforced concrete (RC) buildings, comprising shear walls and gravity frames. The evaluation is undertaken within a performance-based earthquake engineering framework by considering regional seismicity and site-specific ground motion selection. Different engineering demand parameters (EDPs), i.e., maximum interstory drift ratio (MaxISDR) and energy-based damage index, are considered as performance indicators. Various prediction models of EDPs are developed by considering four ground motion intensity measures (IMs), i.e., spectral acceleration at the fundamental period, Arias intensity, cumulative absolute velocity (CAV), and significant duration of ground motion. For this study, a 15-story RC building, located in Vancouver, BC, Canada, is considered as a case study. By using 50 mainshock and 50 mainshock–aftershock (MS-AS) earthquake records (2 horizontal components per record and bidirectional loading), non-linear dynamic analyses are performed. Subsequently, the calculated MaxISDRs and damage indices are correlated with suitable IMs using cloud analysis, and the most efficient IM-EDP prediction models are selected by comparing standard deviations (SDs) of the regression errors. The MaxISDR of the shear walls is less than 1% for the mainshock and MS-AS records. The energy-based damage index shows sensitivity to delineate impact of earthquake types and aftershocks. The CAV is showed to be the most efficient IM for the energy-based damage index.


INTRODUCTION Motivation
Seismic performance of reinforced concrete (RC) shear wall systems designed with Canadian design codes has been investigated by various researchers (e.g., Tremblay et al., 2001;Adebar et al., 2010;Boivin andPaultre, 2010, 2012;Luu et al., 2014).For RC core buildings designed with the CSA standard A23.3-04,Boivin and Paultre (2010) showed that the RC core performs satisfactorily for flexural demand, while potential deficiency under significant shear demand may be a concern.Koduru and Haukaas (2010) studied the seismic performance and economic loss of a 15-story RC building constructed in 1988 and located in Vancouver, BC, Canada.Their study was comprehensive, covering from regional seismic hazard, seismic vulnerability assessment, and economic impact estimation.Nevertheless, important improvements can be made with regard to use of ground motion records that are applicable to megathrust interface records from the Cascadia subduction zone, noting that the records used by Koduru and Haukaas (2010) were calibrated based on shallow crustal records.On the other hand, adopting PEER's performancebased earthquake engineering (PBEE) framework (Cornell and Krawinkler, 2000), Yang et al. (2012) carried out a seismic loss assessment for a 42-story RC dual-system building, i.e., a centrally located core wall building with perimeter special momentresisting frames.With a design earthquake intensity level, the maximum interstory drift ratio (MaxISDR) calculated was less than 2%.
In older Canadian codes, shear wall buildings are the primary seismic force resisting systems, and detailing of gravity frames are often neglected (Adebar et al., 2010).The poor detailing associated with the gravity frames can be (Tesfamariam and Saatcioglu, 2008) inadequate lap splice length, lap splice located in a potential plastic hinge zone, poor detailing of transverse reinforcement anchorage, welded detailing, and lack of support to longitudinal bars.Gravity frames, however, located in the plastic hinge zone of the shear wall can experience excessive deformation and, if not detailed properly, can sustain severe damage (Adebar et al., 2010).This type of damage, for example, was reported in the 27 February 2010 Maule Chile earthquake (Naeim et al., 2011) and the 22 February 2011 Christchurch earthquake (Stirrat et al., 2014).Furthermore, older buildings in Canada lack consideration of large interface events in seismic design procedures.(Note: the potential risk due to the Cascadia subduction earthquakes was only recognized in late 1990s.)The problem is further compounded with the prevalence of mainshock-aftershock (MS-AS) earthquake sequences.
The structural analysis, using an appropriate structural model for the calculation of engineering demand parameter (EDP) and collapse capacity/probability, is an essential component of the seismic vulnerability evaluation.To assess the probability of attaining a specific structural response level conditioned on seismic excitation, incremental dynamic analysis (Vamvatsikos and Cornell, 2002) and cloud and stripe analyses (Jalayer et al., 2007) can be used.Different damage indices are used as a surrogate measure for EDPs that are categorized as follows (Williams and Sexsmith, 1995): (a) non-cumulative, (b) deformation-based cumulative, (c) energy-based cumulative, and (d) combined (non-cumulative and energy-based) indices.The most common approach to relate seismic demands to structural performance limits (i.e., capacity) is based on non-cumulative drift-based EDP, such as MaxISDR and residual interstory drift ratio.For RC shear wall buildings, however, the drift-based damage indicator may show satisfactory seismic resistance performance while underestimating overall damage in the plastic region due to cyclic loading.In this article, the cumulative energy-based damage index proposed by Mehanny and Deierlein (2000), which takes into account both peak amplitude and duration of non-linear responses of structural members in quantifying the structural damage, will be used.Alternatively, other damage indices that were proposed in the literature (e.g., Park and Ang, 1985;DiPasquale and Cakmak, 1989;Reinhorn and Valles, 1995) can be adopted.

Objectives
This article presents a seismic vulnerability evaluation of a 15story RC shear wall structure located in Vancouver, BC, Canada.The RC shear wall building includes gravity frames, which are more vulnerable to low-amplitude repeated ground motions and is modeled in OpenSees finite element software by accounting for the non-linearity in the model.The shear walls for the tall building act as a cantilever beam, and the plastic hinge is formed at the base (assumed to be the first four stories), whereas the gravity columns within this plastic region are modeled with non-linear material elements.The seismic risk assessment that is carried out in this study accounts for: • Consideration of MS-AS earthquake records and earthquake types, i.e., shallow crustal earthquakes, deep inslab earthquakes, and megathrust Cascadia subduction earthquakes, by selecting applicable records for subduction environments from extensive MS-AS ground motion data sets (Goda and Taylor, 2012;Goda et al., 2015), including the 2011 Tohoku earthquake records that can be regarded as closest proxy for the Cascadia subduction events.• Energy-based damage indices, which are computed from the hysteretic responses of the structural model.In this article, Mehanny-Deierlein damage index (Mehanny and Deierlein, 2000) is considered, which captures the responses from longduration earthquakes.The damage indices for different structural elements are primarily integrated based on a combination rule suggested by Bracci et al. (1989).• Impact of different intensity measures (IMs) on efficiency.The cumulative absolute velocity (CAV) is identified as the most efficient IM for characterizing the damage index, and subsequently seismic demand prediction models are developed using CAV.Finally, convoluting the seismic demand model with the seismic hazard of Vancouver, seismic risk is computed quantitatively.

PBEE FOR SEISMIC VULNERABILITY ASSESSMENT
Disaster risk reduction against future earthquakes requires decision support tools for cost-effective risk mitigation options.For seismic risk assessment and design, a PBEE methodology, originally advocated by Cornell and Krawinkler (2000) and later extended by various researchers (e.g., Goulet et al., 2007), can be adopted.In PBEE, the mean annual rate of exceedance of earthquake impact expressed in terms of damage measures (DM) and ν(DM), is quantified, involving seismic hazard, structural, and damage analyses.Mathematically, ν(DM) can be expressed as: where λ(IM) is the mean annual rate that a certain level of IM is exceeded, G(EDP|IM) is the complementary cumulative distribution function of EDP given IM, and G(DM|EDP) is the complementary cumulative distribution function that can be characterized through damage analysis by relating EDP to the physical extent of structural damage, represented by DM.The accuracy of the earthquake impact assessment depends on the available data and the choices of the relevant models and parameters.

Hazard Consideration and IMs
In western BC, Canada, three dominant earthquake sources are present: crustal, inslab, and interface events (Hyndman and Rogers, 2010).The latter two types are originated from the Cascadia subduction zone (Figure 1A), where the oceanic Juan de Fuca Plate sinks beneath the continental North American Plate.When the stored strain along the fault is released, a megathrust subduction earthquake, similar to the 2010 Maule Chile and 2011 Tohoku Japan earthquakes, can happen.The structural damage potential and consequences due to these three earthquake types can be significantly different because of their ground motion characteristics, depending on buildings and infrastructure of interest.Generally, in comparison with crustal and inslab earthquakes, large interface ground motions, originated from the Cascadia subduction zone, have much longer duration (Figure 1B).The spectral content of the ground motion records for three earthquake types can differ significantly due to different earthquake source characteristics in terms of magnitude and distance (Figure 1C).For instance, the effects of long-duration ground motions on tall buildings have been highlighted for the 2011 Tohoku earthquake (Takewaki et al., 2011).
The three earthquake types have different input characteristics.Thus, besides selecting appropriate EDPs, selecting the corresponding IM is important in the overall risk assessment.As such, consideration of spectral acceleration at the fundamental period S a (T 1 ), which is most commonly adopted in modern seismic hazard and risk studies, may not be the most suitable indicator of the energy input from the ground motion.In this article, Arias intensity (AI; Arias, 1970), CAV (Electrical Power Research Institute , EPRI), and significant duration of ground motion (D 5-95% , Trifunac and Brady, 1975) are considered, and their correlations with structural damage are quantified.The definitions of AI and CAV are given as follows:  where a(t) is the accelerogram.D 5-95% is defined as duration between times the AI of a ground motion record reaches 5 and 95% of its final value.

Energy-Based Damage Index
Energy-based damage indices are cumulative and are computed with consideration of hysteretic response (Gosain et al., 1977;Park and Ang, 1985;Kraetzig et al., 1989;Mehanny and Deierlein, 2000; see Table 1).Gosain et al. (1977) formulated a model to describe damage by using energy absorption normalized by yield force and displacement.Park and Ang (1985) proposed a widely used damage index D PA as a linear combination of deformation and absorbed energy under cyclic loading.The weighing factor used for the energy term in D PA is calibrated through experimental work.Kraetzig et al. (1989) developed an energybased damage index D K that accounts for the energy dissipated in primary half cycles (PHCs) and follower half cycles (FHCs) for positive and negative parts of the response.Mehanny and Deierlein (2000) extended the Kraetzig et al. 's damage index by introducing weights on the PHCs and FHCs and the positive and negative damage indices and associated it with extent of physical damage.In this article, following Koduru and Haukaas (2010), the Mehanny-Deierlein damage index D MD is used.The limit states for D MD are shown in Table 2.
To define relationships between different earthquake return periods and acceptable performance limit states, seismic performance matrices are often adopted (Table 3).It is highlighted that for frequent [50% probability of exceedance (PE) in 30 years], occasional (50% PE in 50 years), rare (10% PE in 50 years), and very rare (2% PE in 50 years) earthquake design levels, the corresponding design performance limit states are immediate occupancy (IO), life safety (LS), and collapse prevention (CP), respectively.For the energy-based earthquake damage evaluation, limit states indicated in Table 3 can be used with the corresponding definitions of the damage index (Table 2).

Structural Model
The reference structure considered in this study is a 15-story RC building constructed in 1988 and located in Vancouver, BC, Canada (Ventura et al., 2001;Koduru and Haukaas, 2010).The primary lateral load-resistant element of the building is shear walls (Figure 2).The building is fairly regular in plan, with minor where (positive deformation) Kraetzig damage index for positive and negative parts of the response, respectively; Mehanny and Deierlein ( 2000) , where Mehanny-Deierlein damage index for positive and negative parts of the response, respectively; N + , n + = number of PHC and FHC, respectively, for positive part of the response; N − , n − = number of PHC and FHC, respectively, for negative part of the response; α, β, and γ = calibration parameters setbacks at the fourth and fourteenth story levels (Figure 2A).The mixed-use building has commercial occupancy at the first floor and residential occupancy at the remaining floors.The shear walls in the staircase and elevator shafts are concentrated at the central core and form the main lateral load-resisting system (Figure 2B).
The first story height varies from 2.7 to 4.7 m, and subsequent stories are 2.7 m. Mass and stiffness of the building are used as a base model, and further simplifications are considered in developing a numerical model (Koduru, 2008).For example, four levels of underground parking below grade were not considered in the model, and the foundation was considered to be fixed at base.The numerical model for this building was developed by Koduru and Haukaas (2010) in OpenSees (McKenna et al., 2000).Finite element modeling of the structure was implemented using a fiber element, where the element cross-sections are discretized, and corresponding non-linear material properties of the core concrete, cover concrete, and reinforcing bars were assigned.The structural model consists of three components, i.e., gravity support columns,  shear walls, and header beams.The modeling assumptions made by Koduru (2008) are outlined as follows:

Collapse prevention
• The effect of the RC flat plate slab is accounted for by means of the "rigid diaphragm" option in OpenSees.• The gravity support columns, header beams, and shear walls are modeled as beam-column elements.• The shear walls are designed to yield and form plastic hinges between the first and fourth stories.In the plastic region, the shear walls (first to fourth stories) and gravity columns (first to third stories) are modeled with fiber-discretized sections by accounting for bending moments and axial forces interaction.• From the 4th to 10th story, the shear walls and gravity columns are modeled with non-linear elements.A hysteretic material model includes the force-deformation curve from section analysis and a stiffness degradation factor of 0.05.The shear force-deformation model is separately included in the section models of all non-linear elements.The upper stories, from the 10th to 15th, are modeled as elastic elements.• The header beams, which connect the shear walls, are modeled as non-linear elements with a hysteretic model for the moment-curvature relationship.
Modal analysis is carried out; the first three modal vibration periods and corresponding damping ratios are obtained as follows.The first mode corresponds to the sway motion in the short structural axis direction; its vibration period and damping ratio are 0.9 s and 3.0%, respectively.The second mode is related to the torsional motion, and its vibration period and damping ratio are 0.84 s and 3.0%, respectively.The third mode is the sway mode in the long structural axis direction, and the corresponding vibration period and damping ratio are 0.25 s and 5.7%, respectively.The calculated vibration periods are in agreement with the measured vibration periods of the building by Ventura et al. (2001), i.e., first mode (0.81 s, short structural axis), second mode (0.79 s, torsional), and third mode (0.69 s, long structural axis).These dynamic characteristics are important for selecting ground motion records for use in non-linear dynamic analyses.

Seismic Hazard for Vancouver and Ground Motion Selection
The development of seismic damage prediction models requires a series of non-linear dynamic analyses of a structural model subjected to a set of ground motion records, which reflect the regional seismic hazard of interest.In this article, record selection based on multiple conditional mean spectra (CMS) for different earthquake types is carried out by following the same procedures described in the studies by Tesfamariam et al. (2015) and Tesfamariam and Goda (2015).The target CMS are developed for crustal, interface, and inslab earthquakes, based on full probabilistic seismic hazard analysis (PSHA) results by Atkinson and Goda (2011).The site of interest is Vancouver, and its surface soil is classified as site class C (average top 30 m shear-wave velocity ranges from 360 to 760 m/s).Figures 3A,B show the uniform hazard spectrum and seismic disaggregation result, respectively, at the return period of 2,500 years.The seismic disaggregation is based on spectral acceleration at 0.9 s (i.e., same as the fundamental vibration period of the building model).To develop CMS for different earthquake types, mean record characteristics for individual earthquake types are obtained from the PSHA results.Three CMS are included in Figure 3A, illustrating different spectral characteristics for the three earthquake types.For the considered case (i.e., Vancouver, site class C, vibration period of 0.9 s, and return period of 2,500 years), crustal, interface, and inslab events contribute equally to the overall seismic hazard.
To select a set of suitable ground motion records that match with the constructed CMS, an extended data set of real MS-AS sequences is used, which was developed by combing the worldwide (NGA) database (Goda and Taylor, 2012) with the Japanese (K-KiK-SK) database (Goda et al., 2015).The number of available MS-AS sequences is 606; among them, there are 197 crustal earthquakes, 340 interface earthquakes, and 69 inslab earthquakes.The interface events are from the M w 8.3 2003 Tokachi-oki earthquake and the M w 9.0 2011 Tohoku earthquake (which have similar event characteristics as the expected Cascadia subduction earthquake).
A set of ground motion records is selected by comparing response spectra of candidate records with the target spectra (i.e., CMS).The total number of records is set to 50.(Note: each record has two horizontal components.)The contributions from the 3 earthquake types are equal; as a result, 17, 17, and 16 records are selected for the crustal, interface, and inslab earthquakes, respectively.Response spectral matching is conducted in a least squares sense by considering the geometric mean of the response spectra of two horizontal components.(Note: spectral matching is performed for mainshock records of MS-AS sequences.)The vibration period range for spectral matching is from 0.1 to 2.0 s, which is inclusive of major vibration periods of the tall building model.Figures 3C-E show the statistics of the response spectra of the selected ground motion records (i.e., median as well as 16th and 84th percentile curves).For comparison, Figures 3C-E include the target CMS as well as the CMS ± 1 conditional standard deviation (SD) (Jayaram et al., 2011).The response spectra of the selected records and the target CMS are similar for the crustal and interface records; for inslab records, the selected records contain richer short-period spectral content than the target spectra.Given the availability of ground motion records and the dataset size of ground motion records (i.e., 16-17 for each earthquake type), matching of the candidate response spectra with the target is judged as adequate.Figure 3F shows the response spectra of the unscaled mainshock ground motions that are selected based on the preceding method.It is noted that the mean spectral acceleration at 0.9 s of the 50 records is about 0.36 g, which corresponds to the return period of 1,300 years, while 8 of 50 records exceed the spectral acceleration at 0.9 s that corresponds to the return period of 2,500 years (i.e., 0.5 g).
Figure 4A shows the magnitude-distance distribution of the selected earthquake records; in the figure, record characteristics for mainshocks and major aftershocks (i.e., events having the second largest magnitude within individual sequences) are included.Finally, Figures 4B,C compare the D 5-95% -AI plot and the D 5-95% -CAV plot for different earthquake types (mainshocks only).Figures 4B,C indicate that the interface records are associated with longer duration and larger CAV values than the crustal and inslab records.
The selected ground motion records are used for seismic performance assessment of the tall building in Vancouver.The records reflect regional seismic hazard and dominant record characteristics.In particular, consideration of the 2011 Tohoku records is relevant to the seismic performance assessment in Vancouver because of the anticipated macro-level similarity between the 2011 Tohoku earthquake and possible Cascadia events.The relative contributions from the crustal, interface, and inslab events are equal (for the considered scenario), and thus, these records can also be used for evaluating the effects of ground motion records having different record characteristics (i.e., spectral content and duration) on non-linear seismic demand and earthquake damage potential.The records can be employed in cloud analysis to develop probabilistic seismic demand models.For this purpose, target spectral acceleration levels need to be defined.

Dynamic Analysis and Cloud Analysis
The structural analysis is carried for the 50 MS and 50 MS-AS earthquake records discussed in the previous section.The simulations are carried out using unscaled ground motions in bidirectional horizontal excitations of the 3D model (shown in Figure 2A).Various structural responses are stored for postprocessing, including time history data for interstory drift ratios and floor accelerations at all 16 stories, and Mehanny-Deierlein damage indices in the plastic hinge zone.The bidirectional interstory drift ratios were combined through geometric mean, and the corresponding drift values are used in the subsequent analysis.Figure 5 illustrates calculated time histories of structural responses subjected to the three crustal/interface/inslab ground motion records, which are the same as those shown in Figure 1B. Figure 5A shows results for the interstory drift ratio at the ninth story, while Figure 5B shows results for the Mehanny-Deierlein damage index for gravity column at the third story.The results of the third-story column is selected because the damage index of this structural element is in the middle of all other structural elements of the gravity frame system and is thus suitable to show the average trends.It can be inspected from Figure 5A that under the three input records, the structure behaves elastically as the MaxISDR is only 0.3-0.4%, and there is no significant residual drift after the earthquake sequences.On the other hand, accumulation of cumulative seismic demands to the structure can be observed in Figure 5B; in particular, a major aftershock following the interface mainshock record (middle panel in Figure 5B) shows noticeable increase of the damage index.The maximum values of D MD for the gravity column are still less than 0.3 for the three records; therefore, the structural damage to this element is considered to be negligible for these cases (see Table 2).Figures 6A,B show storywise profiles of interstory drift ratios (left) and relationships between MaxISDR and spectral acceleration (right), S a (T 1 )-MaxISDR, for MS and MS-AS records, respectively.The results are obtained from the cloud analysis, and the responses due to different earthquake types are color coded in the figures.For both MS and MS-AS records, MaxISDR is less than 1% (minimal damage) and tends to take the largest values at the 8th to 10th stories.This is a result of higher mode effects on the response of the structure.The drift is minimal in the plastic hinge zone (the first four stories), as a combined effect of the shear walls and gravity columns.The S a (T 1 )-MaxISDR plots display that MaxISDR is well correlated with S a (T 1 ).In the figure panels, fitted prediction models, having a form of log 10 EDP = a + blog 10 IM, where a and b are the regression coefficients, are also included.As a quantitative measure of efficiency (Luco and Cornell, 2007) Maximum inter-story drift ratio 10 -2 10 -1 10 0 10 1 Spectral acceleration at 0.9 s (g) Mainshock Standard deviation: 0.091

FIGURE 6 | Peak drift demands from cloud analysis displaying storywise profiles of interstory drift ratios and maximum interstory drift ratio and spectral acceleration relationships: (A) MS records and (B) mainshock-aftershock records.
SD of the regression residuals is indicated in the figure panels.It is noteworthy that although detailed results are not shown, similar regression analyses were performed for different IM variables, such as spectral accelerations at different vibration periods, AI, CAV, and D 5-95% .It was found that for MaxISDR, S a (T 1 ) is the most efficient; this conclusion may be due to the elastic responses of the shear wall core under the ground motion records considered.
Overall, the results shown in Figure 6 indicate the effects of aftershocks are not significant in terms of interstory drift.It is noteworthy that the average increase due to major aftershocks is about 8%.However, this increase is mainly caused by only 6 sequences of 50; the majority of the aftershocks do not increase the overall interstory drift.This observation is in agreement with Goda and Taylor (2012) and Goda et al. (2015).Under the considered ground motions, the shear wall core structure remains to be mainly in the linear elastic range, and this indeed reiterates that the shear walls RC cores are not vulnerable (e.g., Yang et al., 2012).For this reason, the subsequent investigations focus upon the seismic damage evaluation of the plastic zone, gravity columns (first to third stories), header beams (second to fourth stories), and shear walls (first to fourth stories) based on the energy-based damage index; see Figure 2.This focus is justified because the members in the plastic zone area are susceptible to severe damage due to cyclic loading (Koduru and Haukaas, 2010).

Efficient IM for Energy-Based Damage Index
Effectiveness of each IM on the estimation of EDPs is assessed using the concept of efficiency.An efficient IM results in relatively small variability of EDP given IM (Luco and Cornell, 2007).This property can be quantified by the SD of the regression residuals of predicted EDP values for a given IM.In this section, the Mehanny-Deierlein damage indices D MD for the structural elements in the plastic zone are considered as EDP.More specifically, in total, 10 damage indices are focused upon; 3 are for the gravity columns (first/second/third story denoted as D MD -C 1 , D MD -C 2 , and D MD -C 3 , respectively), 3 are for the header beams (second/third/fourth story denoted as D MD -HB 2 , D MD -HB 3 , and D MD -HB 4 , respectively), and 4 are for the shear walls (first/second/third/fourth story denoted as D MD -SW 1 , D MD -SW 2 , D MD -SW 3 , and D MD -SW 4 , respectively).On the other hand, four IMs are considered to identify the efficient IM parameters: S a (T 1 ), AI, CAV, and D 5-95% .
To compute efficiency of each IM, the log-linear model, i.e., log 10 EDP = a + blog 10 IM, is fitted using a least squares method.Table 4 summarizes the logarithmic SD of the regression residuals for all combinations of EDP and IM for both MS records and MS-AS records.To show the results graphically, scatter plots of the damage index (D MD -C 3 ) versus four IMs for MS-AS records are presented in Figures 7A-D.The results shown in Table 4 and Figure 7 indicate that the SDs of the residuals are smallest for CAV, followed by AI, S a (T 1 ), and D 5-95% .The results are consistent for all EDPs and for MS/MS-AS records.It can be concluded that CAV is the most efficient IM for D MD .Moreover, the results shown in Figures 7A-D suggest that interface records tend to result in greater damage index values in comparison with crustal and inslab records.This is because interface records are long-duration ground motions (Figure 4C), and thus, their cumulative damage potential is higher than other short-duration ground motions.The consideration of D MD as EDP facilitates the incorporation of cumulative damage modes into the seismic performance evaluation.Figure 7E compares the scatter plots of CAV and D MD -C 3 for MS and MS-AS records, respectively.The results clearly show that the effects of major aftershocks for D MD -C 3 are significant, resulting in increased earthquake damage.Based on the two fitted curves, the average increase of D MD can be quantified as 54%, which is in sharp contrast with the increase of MaxISDR shown in Figure 6.It is important to emphasize the differences of the aftershock effects on D MD -C 3 and MaxISDR.For D MD -C 3 , the effects due to major aftershocks are noticeable for the majority of the cases, rather than a small fraction of the cases (which was applicable to MaxISDR).To demonstrate this clearly, a histogram of the ratios of D MD -C 3 between MS-AS records and MS records is shown in Figure 7F.The results highlight the widespread influence of the major aftershocks on the damage index.

Combined Damage Index
The abovementioned results clearly indicate that the CAV is the efficient IM for all D MD , and the effects of earthquake types (long-duration interface events versus other earthquake types) and aftershocks have major influence on the earthquake damage evaluation.To perform probabilistic seismic risk analysis of the system in the plastic zone of the tall building, a combined measure of earthquake damage needs to be defined.It is noteworthy that the damage index computed for each component of the system represents local damage.
By assigning relative importance or weight to each local damage index, a global damage index can be computed.Park et al. (1987) proposed a story damage index D story as: where D i is the local damage index at location i and E i is the corresponding energy absorbed at location i.The energy dissipated, however, is also incorporated in the damage computation.This index can potentially misrepresent the overall damage state (Williams and Sexsmith, 1995).A more general combination rule for computing a global damage index was proposed by Bracci et al. (1989).They proposed weighing each local damage index by importance weight w i , and corresponding D story is computed as follows: Parameter c is used to give higher importance to the most severely damaged elements.Bracci et al. (1989) suggested c = 1 and an equal weight to structural elements within the same story level.The global damage index is computed for each structural element, i.e., gravity columns, header beams, shear walls, within the plastic region.For simplicity, in this article, in assigning the weight, the number of structural elements at each floor is taken into consideration.For example, at the first and fourth stories, where only two types of structural members are considered, an equal weight of 0.5 can be assigned to each.On the other hand, for the second and third stories, three structural member types are present, and the weight can be assigned as 1/3.Furthermore, the damage index computed for each floor level is extended for the global damage index over the first four stories of the plastic regions, by assigning an equal weight (=0.25) to each floor.
Alternatively, more uniform weighting schemes may be considered.For instance, c = 0 in Eq. 5 corresponds to the arithmetic mean.Another popular choice for a uniform combination rule is the geometric mean of all contributing elements.In the following, these combination rules, in addition to Eq. 5 with c = 1, as proposed by Bracci et al. (1989), will be considered as a part of epistemic uncertainty associated with the structural damage assessment.
Figure 8A compares the 10 local damage indices of the gravity frame system for MS records with the combined damage index based on Eq. 5 with c = 1.It can be observed that the coupling beam at the second story shows the lowest damage index.The shear walls and gravity columns show higher damage indices, whereas the fourth story shear walls exhibit the highest damage potential.2 and 3. D MD,C is calculated based on Eq. 5 with c = 1.

Fragility Curves for Energy-Based Damage Index
The developed CAV-D MD,C models (e.g., Eqs 6 and 7) can be used to derive fragility curves.The distribution of seismic demand about its median is often assumed to follow a two-parameter lognormal probability distribution.After estimating dispersion β D'CAV of the demand about its median, the fragility, i.e., probability that D MD > D MD,C at a given CAV, can be computed as: where Ĉ is the median structural capacity associated with the limit state.By taking the damage limit states shown in Tables 2 and  3 and the combination rule for the global damage index based on Eq. 5 with c = 1 (i.e., Eqs 6 and 7), three fragility curves are developed for the LS threshold (i.e., D MD,C = 0.30), near collapse (NC) threshold (i.e., D MD,C = 0.60), and collapse (C) threshold (i.e., D MD,C = 0.95).The fragility curves derived using Eq. 12 for MS records and MS-AS records are compared in Figure 9.The comparison of the fragility curves indicates that the effects of major aftershocks can be significant.
It is noteworthy that Eq. 12 accounts for statistical uncertainty associated with seismic demand predictions only.On the other hand, there are other important uncertain elements in assessing seismic fragility, such as aleatory uncertainty of capacity C and epistemic modeling uncertainty, denoted by β C and β M , respectively (Ellingwood et al., 2007).To include these effects in evaluating the seismic fragility, the SD β D|CAV in Eq. 12 can be replaced by: β M is assumed to be 0.20, by considering that the modeling process yields an estimate of building frame response that, with 90% confidence, is within ±30% of the actual value (Ellingwood et al.,

2007
). β C is assumed to be 0.25 for IO and LS, following Celik and Ellingwood (2009).For CP, however, β C can be considered to be 0.17 for the four-story gravity frame and 0.13 for the shear walls.Note that the mentioned values of β M and β C are obtained from the literature and are applicable to structural models that were considered therein.Therefore, caution should be exercised before adopting these recommended values of logarithmic SDs.

Probabilistic Seismic Risk Analysis
It is important to assess the seismic performance of the gravity frame system of the tall RC building by taking into account uncertainties associated with regional seismic hazards, ground motions, and seismic vulnerability.The PBEE-based risk analysis methodology, as formulated in Eq. 1, is suitable for such assessments.In carrying out such risk assessments for the tall building in Vancouver, the regional seismic hazard model by Atkinson and Goda (2011) can be used as a starting point.Changes to the ground motion models are necessary because the adopted IM for the gravity frame system is CAV, instead of S a (T 1 ).In this study, two ground motion models for CAV are considered.A model by Campbell and Bozorgnia (2012) is applicable to crustal earthquakes because this model was developed by using strong motion data from worldwide crustal earthquakes compiled in the PEER-NGA database.For the interface and inslab earthquakes, a model by Foulser-Piggott and Goda (2015) can be used because it was developed based on the extensive strong motion database of Japanese earthquakes, including the 2011 Tohoku data.For the seismic vulnerability assessment, the CAV-D MD,C models for MS records and MS-AS records developed in the study (see Eqs 6 and 7) are adopted.The numerical evaluation of Eq. 1 is conducted using Monte Carlo simulations.More specifically, a synthetic earthquake catalog of regional seismicity having five million years is generated from the regional seismicity model of Atkinson and Goda (2011).For a given seismic event in the catalog, a value of CAV is simulated by taking into account earthquake types (i.e., crustal, interface, and inslab).Subsequently, a value of D MD,C is sampled from the developed prediction models (i.e., Eqs 6-11).This is repeated for all seismic events in the catalog.In the CAVbased seismic hazard analysis, the average shear-wave velocity in the uppermost 30 m is set to 555 m/s (i.e., site class C).Once all values of CAV and D MD,C are evaluated, annual maximum hazard and risk values can be extracted to develop a seismic hazard curve (i.e., exceedance probability curve) for CAV [i.e., λ(IM) in Eq. 1] as well as a seismic risk curve for D MD,C [i.e., ν(DM) in Eq. 1].
In addition to the hazard and risk curves, seismic disaggregation plots can be obtained through the postprocessing of the results.
Figure 10 depicts the CAV-based seismic hazard results (i.e., hazard curve and corresponding disaggregation) at the return period of 2,500 years.The seismic hazard curve shown in Figure 10A indicates that the CAV values at the return periods of 500 and 2,500 years correspond to 1,230 and 2990 cm/s, respectively.The disaggregation plot at the 2,500-year return period level highlights significant contribution from the interface events (i.e., events having magnitudes greater than 8.0).The increase is in sharp contrast with the counterpart for S a (T 1 ), shown in Figure 3B.This is because long-duration characteristics of the interface events result in greater values of CAV.To relate the estimated damage potential to the limit states, four ranges of the damage limit states for D MD (i.e., Tables 2 and 3) are indicated along the upper boundary of the figure panel.The results shown in Figures 11A-C indicate that the influence of major aftershocks on the damage potential is significant, increasing the damage index values by approximately 40% for a given probability level.When the mainshock effects only are considered, the return periods that correspond to incipient of LS and NC damage states (i.e., D MD,C = 0.3 and 0.6, respectively) are 650 years and 2050 years.These return period levels are decreased to 350 years and 1,050 years (i.e., greater risks) when the aftershock effects are taken into account in addition to those due to the mainshocks.Comparison of the disaggregation plots for CAV and D MD,C suggests that they are very similar; these are because the differences of the earthquake types effectively capture the long-duration effects, as shown in Figure 7C.
Finally, the effects of the combination rule for defining the global damage index on seismic risk assessment are investigated.For this purpose, seismic risk assessments are carried out by considering three sets of seismic demand prediction models of D MD,C , i.e., Eqs 6 and 7 versus Eqs 8 and 9 versus Eqs 10 and 11.The results are shown in Figures 11D,E.The base case (i.e., Eqs 6 and 7) leads to greater damage index values, in comparison with the two other cases (i.e., arithmetic mean and geometric mean) because more weights are given to severely damaged structural elements.The differences of the seismic risk curves for the three cases are significant, highlighting the importance of capturing this uncertainty in seismic risk assessments.

CONCLUSION
Seismic performance of an RC shear wall system designed with Canadian design codes has shown acceptable performance in terms of drift limits.Recent damaging earthquakes have highlighted that MS-AS earthquake records are important factors in the overall risk assessment.Furthermore, gravity columns, which are not seismically detailed and thus have exhibited severe damage, can potentially lead to localized collapse.However, driftbased limit states showed little sensitivity to MS-AS earthquake records and impact of earthquake types.In this article, an energybased damage index is considered to capture the effects of longduration earthquake ground motions.It is important to assess the seismic performance of the gravity frame system of the tall RC building by taking into account uncertainties associated with regional seismic hazards, ground motions, and seismic vulnerability.The PBEE-based risk analysis methodology is suitable for such assessments.Thus, in this article, seismic performance of the 15-story RC shear wall building located in Vancouver, BC, Canada, was investigated.For the seismicity of Vancouver, BC, Canada, scope of the study and conclusions are summarized below.

Frequent
PE in 50 years) ♢ × Very rare (2% PE in 50 years) ♢ , basic objective-proposed National Building Code of Canada (NBCC) normal importance; , essential service objective-proposed NBCC high importance; ♢, safety critical objective-not proposed NBCC category; ×, unacceptable performance for new construction.
FIGURE 2 | (A) Schematic view of the Heritage Court Tower building.(B) Schematic of a typical floor plan.

FIGURE 5 |
FIGURE 3 | (A)Uniform hazard spectrum and conditional mean spectra for crustal, interface, and inslab events in Vancouver (site class C; return period of 2,500 years).(B) Seismic disaggregation for spectral acceleration at 0.9 s in Vancouver (site class C; return period of 2,500 years).Comparison of response spectra of the selected (mainshock) ground motion records (50th/16th/84th curve) and conditional mean spectra (mean and mean ± 1 SD): (C) crustal earthquakes, (D) interface earthquakes, and (E) inslab earthquakes.(F) Geometric mean response spectra of the selected mainshock ground motions.

FIGURE 7 |
FIGURE 7 | Scatter plots of damage index (D MD -C 3 ) versus four intensity measures by considering mainshock-aftershock (MS-AS) records: (A) spectral acceleration at 0.9 s, (B) Arias intensity, (C) cumulative absolute velocity, and (D) significant duration.(E) Scatter plot of damage index (D MD -C 3 ) versus cumulative absolute velocity by considering MS records and MS-AS records.(F) Histogram of the damage index ratios (D MD -C 3 ) between MS-AS records and MS records.

FIGURE 8 |
FIGURE 8 | Comparison of local damage indices with the combined damage index in the plastic hinge zone for MS records.(B) Comparison of three combined damage indices in the plastic hinge zone for MS records.
for MS records, and log 10 D MD,C = −3.741+ 1.011log 10 CAV (9) for MS-AS records.The SDs for Eqs 8 and 9 are 0.120 and 0.115, respectively.Moreover, by adopting the geometric mean combination rule, seismic demand prediction models for D MD,C can be obtained as: log 10 D MD,C = −4.159+ 1.032log 10 CAV (10) for MS records, and log 10 D MD,C = −3.911+ 1.008log 10 CAV (11) for MS-AS records.The SDs for Eqs 10 and 11 are 0.123 and 0.110, respectively.

FIGURE 9 |
FIGURE 9 | CAV-D MD,C fragility models for MS and mainshockaftershock records based on the damage index thresholds listed in Tables2 and 3. D MD,C is calculated based on Eq. 5 with c = 1.

Figure 11
Figure11shows the D MD,C -based seismic risk results.Two risk curves for MS records and MS-AS records are presented in Figure11A, whereas the seismic risk disaggregation plots for MS records and MS-AS records at the 2,500-year return period are shown in Figure11Band Figure11C, respectively.The seismic demand prediction models used for obtaining the results shown in Figures11A-Care Eqs 6 and 7. To relate the estimated damage potential to the limit states, four ranges of the damage limit states for D MD (i.e., Tables2 and 3) are indicated along the upper boundary of the figure panel.The results shown in Figures11A-Cindicate that the influence of major aftershocks on the damage potential is significant, increasing the damage index values by approximately 40% for a given probability level.When the mainshock effects only are considered, the return periods that correspond to incipient of LS and NC damage states (i.e., D MD,C = 0.3 and 0.6, respectively) are 650 years and 2050 years.These return period levels are decreased to 350 years and 1,050 years (i.e., greater risks) when the aftershock effects are taken into account in addition to those due to the mainshocks.Comparison of the disaggregation plots for CAV and D MD,C suggests that they are very similar; these are because the differences of the earthquake types effectively capture the long-duration effects, as shown in Figure7C.Finally, the effects of the combination rule for defining the global damage index on seismic risk assessment are investigated.For this purpose, seismic risk assessments are carried out by considering three sets of seismic demand prediction models of D MD,C , i.e., Eqs 6 and 7 versus Eqs 8 and 9 versus Eqs 10 and 11.The results are shown in Figures11D,E.The base case (i.e., Eqs 6 and 7) leads to greater damage index values, in comparison with the two other cases (i.e., arithmetic mean and geometric mean) because more weights are given to severely damaged structural elements.The differences of the seismic risk curves for the three cases are significant, highlighting the importance of capturing this uncertainty in seismic risk assessments.