A performance‐based framework for adaptive seismic aftershock risk assessment

Operative seismic aftershock risk forecasting can be particularly useful for rapid decision‐making in the presence of an ongoing sequence. In such a context, limit state first‐excursion probabilities (risk) for the forecasting interval (a day) can represent the potential for progressive state of damage in a structure. This work lays out a performance‐based framework for adaptive aftershock risk assessment in the immediate post‐mainshock environment. A time‐dependent structural performance variable is adopted in order to measure the cumulative damage in a structure. A set of event‐dependent fragility curves as a function of the first‐mode spectral acceleration for a prescribed limit state is calculated by employing back‐to‐back nonlinear dynamic analyses. An epidemic‐type aftershock sequence model is employed for estimating the spatio‐temporal evolution of aftershocks. The event‐dependent fragility curves for a given limit state are then integrated together with the probability distribution of aftershock spectral acceleration based on the epidemic‐type aftershock sequence aftershock hazard. The daily probability of limit state first‐excursion is finally calculated as a weighted combination of the sequence of limit state probabilities conditioned on the number of aftershocks. As a numerical example, daily aftershock risk is calculated for the L'Aquila 2009 aftershock sequence (central Italy). A representative three‐story reinforced concrete frame with infill panels, which has cyclic strength and stiffness degradation, is used in order to evaluate the progressive damage. It is observed that the proposed framework leads to a sound forecasting of limit state first‐excursion in the structure for two limit states of significant damage and near collapse. Copyright © 2014 John Wiley & Sons, Ltd.     

Effects of aftershocks on peak ductility demand due to strong ground motion records from shallow crustal earthquakes

Aftershocks induced by a large mainshock can cause additional damage to structures and infrastructure, hampering building reoccupation and restoration activities in a post‐disaster situation. To assess the nonlinear damage potential due to aftershocks, this study investigates the effects of aftershocks by using real as well as artificially generated mainshock–aftershock sequences. The real mainshock–aftershock sequences are constructed from the Pacific Earthquake Engineering Research Center—Next Generation Attenuation database for worldwide shallow crustal earthquakes; however, they are deemed to be incomplete because of missing records. To supplement incomplete real dataset, artificial sequences are generated on the basis of the generalized Omori's law, and a suitable aftershock record selection procedure is then devised to simulate time‐series data for mainshock–aftershock sequences. The results from nonlinear dynamic analysis of inelastic single‐degree‐of‐freedom systems using real and artificial sequences indicate that the incremental effects of aftershocks on peak ductility demand using the real sequences are relatively minor and that peak ductility demand estimates based on the generalized Omori's law are greater, particularly in the upper tail, than those for the real sequences. The results based on the generalized Omori's law also suggest that the aftershock effects based on the real sequences might underestimate the aftershock impact because of the incompleteness of the real dataset. Copyright © 2012 John Wiley & Sons, Ltd.     

Aftershock risk assessment and the decision to open traffic on bridges

Critical issues in emergency management after a seismic event are assessing the functionality of the main infrastructures (hospitals, road network, etc.) and deciding on their usability just after the mainshock. The use of a pure analytical tool to assess the aftershock risk of a structure can be contrasted with the limited time available to make a decision about the usability of a structure. For this reason, this paper presents a method for evaluating post‐earthquake bridge practicability based on a rational combination of information derived from numerical analyses and in situ inspections. In particular, we propose an effective tool to speed up the decision‐making process involved in evaluating the seismic risk of mainshock‐damaged bridges in the context of aftershocks. The risk is calculated by combining the aftershock hazard using the Omori law and the fragility curves of the structure, which are calculated using the regression analysis of a sample of results obtained with data randomly generated by the Latin Hypercube Sampling technique and updated based on the results of in situ inspection. Different decision criteria regarding the practicability of bridges are discussed, and a new criterion is proposed. This tool was applied to an old highway RC viaduct. There are two main findings, including the high sensitivity to Bayesian updating (especially when the damage predicted by numerical analysis does not match the real damage) and the criteria used to decide when re‐open bridges to traffic. Copyright © 2013 John Wiley & Sons, Ltd.     

Closed‐form aftershock reliability of damage‐cumulating elastic‐perfectly‐plastic systems

Major earthquakes (i.e., mainshocks) typically trigger a sequence of lower magnitude events clustered both in time and space. Recent advances of seismic hazard analysis stochastically model aftershock occurrence (given the main event) as a nonhomogeneous Poisson process with rate that decays in time as a negative power law. Risk management in the post‐event emergency phase has to deal with this short‐term seismicity. In fact, because the structural systems of interest might have suffered some damage in the mainshock, possibly worsened by damaging aftershocks, the failure risk may be large until the intensity of the sequence reduces or the structure is repaired. At the state‐of‐the‐art, the quantitative assessment of aftershock risk is aimed at building tagging, that is, to regulate occupancy. The study, on the basis of age‐dependent stochastic processes, derived closed‐form approximations for the aftershock reliability of simple nonevolutionary elastic‐perfectly‐plastic damage‐cumulating systems, conditional on different information about the structure. Results show that, in the case hypotheses apply, the developed models may represent a basis for handy tools enabling risk‐informed tagging by stakeholders and decision makers. Copyright © 2013 John Wiley & Sons, Ltd.     

Seismic risk assessment considering cumulative damage due to aftershocks

Calculating the limit state (LS) exceedance probability for a structure considering the main seismic event and the triggered aftershocks (AS) is complicated both by the time‐dependent rate of aftershock occurrence and also by the cumulative damage caused by the sequence of events. Taking advantage of a methodology developed previously by the authors for post‐mainshock (MS) risk assessment, the LS probability due to a sequence of mainshock and the triggered aftershocks is calculated for a given aftershock forecasting time window. The proposed formulation takes into account both the time‐dependent rate of aftershock occurrence and also the damage accumulation due to the triggered aftershocks. It is demonstrated that an existing reinforced concrete moment‐resisting frame with infills subjected to the main event and the triggered sequence exceeds the near‐collapse LS. On the other hand, the structure does not reach the onset of near‐collapse LS when the effect of triggered aftershocks is not considered. It is shown, based on simplifying assumptions, that the derived formulation yields asymptotically to the same Poisson‐type functional form used when the cumulative damage is not being considered. This leads to a range of approximate solutions by substituting the fragilities calculated for intact, MS‐damaged, and MS‐plus‐one‐AS‐damaged structures in the asymptotic simplified formulation. The latter two approximate solutions provide good agreement with the derived formulation. Even when the fragility of intact structure is employed, the approximate solution (considering only the time‐dependent rate of aftershock occurrence) leads to higher risk estimates compared with those obtained based on only the mainshock. Copyright © 2016 John Wiley & Sons, Ltd.     

Stochastic procedure for the simulation of synthetic main shock‐aftershock ground motion sequences

According to the current seismic codes, structures are designed to resist the first damaging earthquake during their service life. However, after a strong main shock, a structure may still face damaging aftershocks. The main shock‐aftershock sequence may result in major damage and eventually the collapse of a structure. Current studies on seismic hazard mainly focus on the modeling and simulation of main shocks. This paper proposes a 3‐step procedure to generate main shock‐aftershock sequences of pairs of horizontal components of a ground motion at a site of interest. The first step generates ground motions for the main shock using either a source‐based or site‐based model. The second step generates sequences of aftershocks' magnitudes, locations, and times of occurrence using either a fault‐based or seismicity‐based model. The third step simulates pairs of ground motion components using a new empirical model proposed in this paper. We develop prediction equations for the controlling parameters of a ground motion model, where the predictors are the site condition and the aftershock characteristics from the second step. The coefficients in the prediction equations and the correlation between the model parameters (of the 2 horizontal components of 1 record and of several records in 1 sequence) are estimated using a database of aftershock accelerograms. A backward stepwise deletion method is used to simplify the initial candidate prediction equations and avoid overfitting the data. The procedure, based on easily identifiable engineering parameters, is a useful tool to incorporate effects of aftershocks into seismic analysis and design.     

Building service life economic loss assessment under sequential seismic events

Aftershocks have been shown to exacerbate earthquake‐induced financial losses by causing further damage to structural and nonstructural components in buildings that have already been affected by a mainshock event and increasing the duration of disrupted functionality. Whereas seismic loss assessment under isolated events has been addressed thoroughly in previous studies, comparatively less has been accomplished in the area of loss assessment under sequences of mainshock‐aftershock ground motions. The main objective of the current study is to formulate a comprehensive framework for quantifying financial losses under sequential seismic events. The proposed framework is capable of accounting for the uncertainties in the state of structure due to accumulation of earthquake‐induced damage, the time‐dependent nature of seismic hazard in the post‐mainshock environment, and the uncertainties in the occurrence of mainshock and aftershock events. Application of the proposed framework to a 4‐story reinforced concrete moment frame shows that consideration of aftershocks could increase lifecycle earthquake‐induced losses by up to 30% compared with mainshock‐only assessments.     

Aftershock collapse vulnerability assessment of reinforced concrete frame structures

In a seismically active region, structures may be subjected to multiple earthquakes, due to mainshock–aftershock phenomena or other sequences, leaving no time for repair or retrofit between the events. This study quantifies the aftershock vulnerability of four modern ductile reinforced concrete (RC) framed buildings in California by conducting incremental dynamic analysis of nonlinear MDOF analytical models. Based on the nonlinear dynamic analysis results, collapse and damage fragility curves are generated for intact and damaged buildings. If the building is not severely damaged in the mainshock, its collapse capacity is unaffected in the aftershock. However, if the building is extensively damaged in the mainshock, there is a significant reduction in its collapse capacity in the aftershock. For example, if an RC frame experiences 4% or more interstory drift in the mainshock, the median capacity to resist aftershock shaking is reduced by about 40%. The study also evaluates the effectiveness of different measures of physical damage observed in the mainshock‐damaged buildings for predicting the reduction in collapse capacity of the damaged building in subsequent aftershocks. These physical damage indicators for the building are chosen such that they quantify the qualitative red tagging (unsafe for occupation) criteria employed in post‐earthquake evaluation of RC frames. The results indicated that damage indicators related to the drift experienced by the damaged building best predicted the reduced aftershock collapse capacities for these ductile structures. Copyright © 2014 John Wiley & Sons, Ltd.     

Post‐repair effect of column jackets on aftershock fragilities of damaged RC bridges subjected to successive earthquakes

In light of recent earthquakes, structures damaged during an initial seismic event (mainshock) may be more vulnerable to severe damage and collapse during a subsequent event (aftershock). In this paper, a framework for the development of aftershock fragilities is presented; these aftershock fragilities define the likelihood that a bridge damaged during an initial event will exhibit a given damage state following one or more subsequent events. The framework is capable of (i) quantifying the cumulative damage of unrepaired bridges subjected to mainshock–aftershock sequences (effect of multiple earthquakes) and (ii) evaluating the effectiveness of column repair schemes such as steel and fiber‐reinforced‐polymer jackets (post‐repair effect of jackets). To achieve this aim, the numerical model of repaired columns is validated using existing experimental results. A non‐seismically designed bridge is chosen as a case study and is modeled for three numerical bridge models: a damaged (but unrepaired) bridge model, and two bridge models with columns repaired with steel and fiber‐reinforced polymer jackets. A series of back‐to‐back dynamic analyses under successive earthquakes are performed for each level of existing damage. Using simulated results, failure probabilities of components for multiple limit states are computed for each bridge model and then are used to evaluate the relative vulnerability of components associated with cumulative damage and column repair. Copyright © 2016 John Wiley & Sons, Ltd.     

Wavelet‐based simulation of spectrum‐compatible aftershock accelerograms

In damage‐based seismic design it is desirable to account for the ability of aftershocks to cause further damage to an already damaged structure due to the main shock. Availability of recorded or simulated aftershock accelerograms is a critical component in the non‐linear time‐history analyses required for this purpose, and simulation of realistic accelerograms is therefore going to be the need of the profession for a long time to come. This paper attempts wavelet‐based simulation of aftershock accelerograms for two scenarios. In the first scenario, recorded main shock and aftershock accelerograms are available along with the pseudo‐spectral acceleration (PSA) spectrum of the anticipated main shock motion, and an accelerogram has been simulated for the anticipated aftershock motion such that it incorporates temporal features of the recorded aftershock accelerogram. In the second scenario, a recorded main shock accelerogram is available along with the PSA spectrum of the anticipated main shock motion and PSA spectrum and strong motion duration of the anticipated aftershock motion. Here, the accelerogram for the anticipated aftershock motion has been simulated assuming that temporal features of the main shock accelerogram are replicated in the aftershock accelerograms at the same site. The proposed algorithms have been illustrated with the help of the main shock and aftershock accelerograms recorded for the 1999 Chi–Chi earthquake. It has been shown that the proposed algorithm for the second scenario leads to useful results even when the main shock and aftershock accelerograms do not share the same temporal features, as long as strong motion duration of the anticipated aftershock motion is properly estimated. Copyright © 2008 John Wiley & Sons, Ltd.     

Spatiotemporal seismic hazard and risk assessment of M9.0 megathrust earthquake sequences of wood‐frame houses in Victoria,British Columbia,Canada

Megathrust earthquake sequences, comprising mainshocks and triggered aftershocks along the subduction interface and in the overriding crust, can impact multiple buildings and infrastructure in a city. The time between the mainshocks and aftershocks usually is too short to retrofit the structures; therefore, moderate‐size aftershocks can cause additional damage. To have a better understanding of the impact of aftershocks on city‐wide seismic risk assessment, a new simulation framework of spatiotemporal seismic hazard and risk assessment of future M9.0 sequences in the Cascadia subduction zone is developed. The simulation framework consists of an epidemic‐type aftershock sequence (ETAS) model, ground‐motion model, and state‐dependent seismic fragility model. The spatiotemporal ETAS model is modified to characterise aftershocks of large and anisotropic M9.0 mainshock ruptures. To account for damage accumulation of wood‐frame houses due to aftershocks in Victoria, British Columbia, Canada, state‐dependent fragility curves are implemented. The new simulation framework can be used for quasi‐real‐time aftershock hazard and risk assessments and city‐wide post‐event risk management.     

Framework of aftershock fragility assessment–case studies: older California reinforced concrete building frames

Current seismic design codes and damage estimation tools neglect the influence of successive events on structures. However, recent earthquakes have demonstrated that structures damaged during an initial event (mainshock) are more vulnerable to severe damage and collapse during a subsequent event (aftershock). This increased vulnerability to damage translates to increased likelihood of loss of use, property, and life. Thus, a reliable risk assessment tool is required that characterizes the risk of the undamaged structure subjected to an initial event and the risk of the damaged structure under subsequent events. In this paper, a framework for development of aftershock fragilities is presented; these aftershock fragilities define the likelihood that a building damaged during a mainshock will exhibit a given damage state following one or more aftershocks. Thus, the framework provides a method for characterizing the risk associated with damage accumulation in the structure. The framework includes the following: (i) creation of a numerical model of the structure; (ii) characterization of building damage states; (iii) generation of a suite of mainshock–aftershocks; (iv) mainshock–aftershock analyses; and (v) development of aftershock fragility curves using probabilistic aftershock demand models, defined as a linear regression of aftershock demand–intensity pairs in a logarithmic space, and damage‐state prediction models. The framework is not limited to a specific structure type but requires numerical models defining structural response and linking structural response with damage. In the current study, non‐ductile RC frames (low‐rise, mid‐rise, and high‐rise) are selected as case studies for the application of the framework. Copyright © 2015 John Wiley & Sons, Ltd.     

Aftershock collapse fragility curves for non‐ductile RC buildings: a scenario‐based assessment

Recent studies have addressed the computation of fragility curves for mainshock (MS)‐damaged buildings. However, aftershock (AS) fragilities are generally conditioned on a range of potential post‐MS damage states that are simulated via static or dynamic analyses performed on an intact building. Moreover, there are very few cases where the behavior of non‐ductile reinforced concrete buildings is analyzed. This paper presents an evaluation of AS collapse fragility conditioned on various return periods of MSs, allowing for the rapid assessment of post‐earthquake safety variations based solely on the intensity of the damaging earthquake event. A refined multi‐degree‐of‐freedom model of a seven‐storey non‐ductile building, which includes brittle failure simulations and the evaluation of a system level collapse, is adopted. Aftershock fragilities are obtained by performing an incremental dynamic analysis for a number of MS–AS ground motion sequences and a variety of MS intensities. The AS fragilities show that the probability of collapse significantly increases for higher return periods for the MS. However, this result is mainly ascribable to collapses occurred during MSs. When collapse cases that occur during a MS are not considered in the assessment of AS collapse probability, a smaller shift in the fragility curves is observed as the MS intensity increases. This result is justified considering the type of model and collapse modes introduced, which strongly depend on the brittle behavior of columns failing in shear or due to axial loads. The analysis of damage that is due to MSs when varying the return period confirms this observation. Copyright © 2017 John Wiley & Sons, Ltd.     

唐山东矿区地震易损性分析

地震易损性分析是小区域震害预测的一种新方法,是以震害经验为基础,综合分析地基土、第四系厚度、地下水埋深、砂土液化、边坡、断层和地形等地质因素对震害的影响,并将这些影响以数字形式表示。以唐山东矿区为例,介绍地震易损性分析在我国的应用     

Evaluation of a recently proposed record selection and scaling procedure for low‐rise to mid‐rise reinforced concrete buildings and its use for probabilistic risk assessment studies

This paper evaluates a recent record selection and scaling procedure of the authors that can determine the probabilistic structural response of buildings behaving either in the elastic or post‐elastic range. This feature marks a significant strength on the procedure as the probabilistic structural response distribution conveys important information on probability‐based damage assessment. The paper presents case studies that show the utilization of the proposed record selection and scaling procedure as a tool for the estimation of damage states and derivation of site‐specific and region‐specific fragility functions. The method can be used to describe exceedance probabilities of damage limits under a certain target hazard level with known annual exceedance rate (via probabilistic seismic hazard assessment). Thus, the resulting fragility models can relate the seismicity of the region (or a site) with the resulting building performance in a more accurate manner. Under this context, this simple and computationally efficient record selection and scaling procedure can be benefitted significantly by probability‐based risk assessment methods that have started to be considered as indispensable for developing robust earthquake loss models. Copyright © 2013 John Wiley & Sons, Ltd.     

Post‐earthquake functionality of highway overpass bridges

Bridges are crucial to the transportation network in a region struck by an earthquake. Collapse of a bridge determines if a road is passable. Ability of a bridge to carry traffic load after an earthquake determines the weight and speed of vehicles that can cross it. Extent of system and component structural damage in bridges determines the cost and time required for repair. Today, post‐earthquake bridge evaluation is qualitative rather than quantitative. The research presented in this paper aims to provide a quantitative engineering basis for quick and reliable evaluation of the ability of a typical highway overpass bridge to function after an earthquake. The Pacific Earthquake Engineering Research (PEER) Center's probabilistic performance‐based evaluation approach provides the framework for post‐earthquake bridge evaluation. An analytical study was performed that linked engineering demand parameters to earthquake intensity measures. The PEER structural performance database and reliability analysis tools were then used to link demand parameters to damage measures. Finally, decision variables were developed to describe three limit states, repair cost, traffic function, and collapse, in terms of induced damage. This paper presents the analytical models used to evaluate post‐earthquake bridge function, decision variables and their correlation to the considered limit states, and fragility curves that represent the probability of exceeding a bridge function limit state given an earthquake intensity. Copyright © 2005 John Wiley & Sons, Ltd.     

Damage spectra for the mainshock–aftershock sequence-type ground motions

Strong aftershocks have the potential to increase the damage state of the structures due to the damage accumulation. This paper investigates the damage spectra for the mainshock–aftershock sequence-type ground motions with Park–Ang damage index. A method of simulating the mainshock–aftershock sequence-type ground motions is proposed based on the modified form of Bath's law and NGA ground motion prediction equation. The damage spectra are computed using the recorded and simulated sequence-type ground motions, and the effects of period of vibration, strength reduction factor, site condition, seismic sequence, damping ratio and post-yield stiffness on damage spectra are studied statistically. The results indicate that the effect of aftershock on structural damage is significant and recorded sequence-type ground motions may underestimate the damage of long-period structures due to the incompleteness of dataset. A simplified equation is proposed to facilitate the application of damage spectra in the seismic practice for mainshock–aftershock sequence-type ground motions.     

基于速率-状态依从摩擦定律的前郭震群余震活动率及模型参数相关性研究

基于速率-状态依从摩擦定律的地震活动率时空预测模型,以同震库伦应力变化作为模型初始应力扰动,模拟了2013年吉林前郭MS5.8震群的余震活动率变化。考虑模型参数相关性,在模拟中采用2种不同的拟合方案,一是余震持续时间t_a不固定条件下的拟合,二是余震持续时间t_a固定条件下的拟合。结果显示,t_a不固定条件下的拟合方式可获得较好的AIC评价,适用于震后早期的趋势判定; t_a固定条件下的拟合计算耗时更短,拟合误差更小,理论模拟结果与前郭震群实际地震时序特征更为吻合。采用该方案对截至2016年10月24日的余震活动率变化进行了回溯性预测检验,结果显示模型预期的余震日频次与实际记录呈较好的正相关关系。研究还发现,主震破裂面附近的同震应力影区导致震后早期模型预测值相对于实际偏低,说明前郭序列余震活动可能还存在其他触发机制。     

人工主余震作用下桥墩结构损伤预测评估分析

主余震作用下,断层距（R_JB）和余震次数的不同对结构造成的损伤也大不相同,进行主余震序列下桥梁损伤评估时需要考虑断层距、余震次数等影响因素。基于OpenSees平台,以一座连续梁桥为例,根据人工主余震构造时不同的影响因素,选用余震衰减和PGA调幅的人工主余震构造法构造主余震序列进行结构损伤分析,并与实际主余震事件的损伤结果作对比,研究人工构造主余震序列预测结构损伤的可行性。结果表明：基于PGA调幅构造主余震的方法,当考虑余震衰减时,能够较好模拟出结构在实际主余震地震作用下造成的最终损伤,能够较合理的预测结构在实际主余震序列中的损伤指标;对区域主余震作用下的桥梁进行损伤评估时,需要考虑断层距和余震次数影响因素,桥梁距断层距R_JB越近时结构的损伤状态越严重,余震发生次数越多则结构的损伤增量就越大,并且结构在人工主余震作用下R_JB越大的地方,结构损伤增量误差程度越小,预测的损伤状态就越接近实际的损伤状态。