Evaluation of the exact conditional spectrum and generalized conditional intensity measure methods for ground motion selection

Two existing, contemporary ground motion selection and modification procedures – (i) exact conditional spectrum (CS‐exact) and (ii) generalized conditional intensity measure (GCIM) – are evaluated in their ability to accurately estimate seismic demand hazard curves (SDHCs) of a given structure at a specified site. The amount of effort involved in implementing these procedures to compute a single SDHC is studied, and a case study is chosen where rigorous benchmark SDHCs can be determined for evaluation purposes. By comparing estimates from ground motion selection and modification procedures with the benchmark, we conclude that estimates from CS‐exact are unbiased in many of the cases considered. The estimates from GCIM are even more accurate, as they are unbiased for most – but not all – of the cases where estimates from CS‐exact are biased. We find that it is possible to obtain biased SDHCs from GCIM, even after employing a very diverse collection of intensity measures to select ground motions and implementing its bias‐checking feature, because it is usually difficult to identify intensity measures that are truly ‘sufficient’ for the response of a complex, multi‐degree‐of‐freedom system. Copyright © 2015 John Wiley & Sons, Ltd.     

Ground motion selection for simulation‐based seismic hazard and structural reliability assessment

This paper examines four methods by which ground motions can be selected for dynamic seismic response analyses of engineered systems when the underlying seismic hazard is quantified via ground motion simulation rather than empirical ground motion prediction equations. Even with simulation‐based seismic hazard, a ground motion selection process is still required in order to extract a small number of time series from the much larger set developed as part of the hazard calculation. Four specific methods are presented for ground motion selection from simulation‐based seismic hazard analyses, and pros and cons of each are discussed via a simple and reproducible illustrative example. One of the four methods (method 1 ‘direct analysis’) provides a ‘benchmark’ result (i.e., using all simulated ground motions), enabling the consistency of the other three more efficient selection methods to be addressed. Method 2 (‘stratified sampling’) is a relatively simple way to achieve a significant reduction in the number of ground motions required through selecting subsets of ground motions binned based on an intensity measure, IM. Method 3 (‘simple multiple stripes’) has the benefit of being consistent with conventional seismic assessment practice using as‐recorded ground motions, but both methods 2 and 3 are strongly dependent on the efficiency of the conditioning IM to predict the seismic responses of interest. Method 4 (‘generalized conditional intensity measure‐based selection’) is consistent with ‘advanced’ selection methods used for as‐recorded ground motions and selects subsets of ground motions based on multiple IMs, thus overcoming this limitation in methods 2 and 3. Copyright © 2015 John Wiley & Sons, Ltd.     

Evaluation of ground motion selection and modification procedures using synthetic ground motions

This study presents a novel approach for evaluating ground motion selection and modification (GMSM) procedures in the context of probabilistic seismic demand analysis. In essence, synthetic ground motions are employed to derive the benchmark seismic demand hazard curve (SDHC), for any structure and response quantity of interest, and to establish the causal relationship between a GMSM procedure and the bias in its resulting estimate of the SDHC. An example is presented to illustrate how GMSM procedures may be evaluated using synthetic motions. To demonstrate the robustness of the proposed approach, two significantly different stochastic models for simulating ground motions are considered. By quantifying the bias in any estimate of the SDHC, the proposed approach enables the analyst to rank GMSM procedures in their ability to accurately estimate the SDHC, examine the sufficiency of intensity measures employed in ground motion selection, and assess the significance of the conditioning intensity measure in probabilistic seismic demand analysis. Copyright © 2015 John Wiley & Sons, Ltd.     

A framework for the evaluation of ground motion selection and modification procedures

This study develops a framework to evaluate ground motion selection and modification (GMSM) procedures. The context is probabilistic seismic demand analysis, where response history analyses of a given structure, using ground motions determined by a GMSM procedure, are performed in order to estimate the seismic demand hazard curve (SDHC) for the structure at a given site. Currently, a GMSM procedure is evaluated in this context by comparing several resulting estimates of the SDHC, each derived from a different definition of the conditioning intensity measure (IM). Using a simple case study, we demonstrate that conclusions from such an approach are not always definitive; therefore, an alternative approach is desirable. In the alternative proposed herein, all estimates of the SDHC from GMSM procedures are compared against a benchmark SDHC, under a common set of ground motion information. This benchmark SDHC is determined by incorporating a prediction model for the seismic demand into the probabilistic seismic hazard analysis calculations. To develop an understanding of why one GMSM procedure may provide more accurate estimates of the SDHC than another procedure, we identify the role of ‘IM sufficiency’ in the relationship between (i) bias in the SDHC estimate and (ii) ‘hazard consistency’ of the corresponding ground motions obtained from a GMSM procedure. Finally, we provide examples of how misleading conclusions may potentially be obtained from erroneous implementations of the proposed framework. Copyright © 2014 John Wiley & Sons, Ltd.     

Intensity measures for probabilistic assessment of non‐structural components acceleration demand

A fundamental issue in the framework of seismic probabilistic risk analysis is the choice of ground motion intensity measures (IMs). Based on the floor response spectrum method, the present contribution focuses on the ability of IMs to predict non‐structural components (NSCs) horizontal acceleration demand. A large panel of IMs is examined and a new IM, namely equipment relative average spectral acceleration (E‐ASA_R), is proposed for the purpose of NSCs acceleration demand prediction. The IMs efficiency and sufficiency comparisons are based on (i) the use of a large dataset of recorded earthquake ground motions; (ii) numerical analyses performed on three‐dimensional numerical models, representing actual structural wall and frame buildings; and (iii) systematic statistical analysis of the results. From the comparative study, the herein introduced E‐ASA_R shows high efficiency with respect to the estimation of maximum floor response spectra ordinates. Such efficiency is particularly remarkable in the case of structural wall buildings. Besides, the sufficiency and the simple formulation allowing the use of existing ground motion prediction models make the E‐ASA_R a promising IMs for seismic probabilistic risk assessment. Copyright © 2015 John Wiley & Sons, Ltd.     

A new seismic hazard analysis using FOSM algorithms

From recent lessons, it is evident that earthquake prediction is immature and impractical as of now. Under the circumstances, seismic hazard analysis is considered a more practical approach for earthquake hazard mitigation, by estimating the annual rate of earthquake ground motions (or seismic hazard) based on seismicity and other geological evidences. Like other earthquake studies for the high-seismicity region around Taiwan, this study aims to conduct a new seismic hazard assessment for the region using the well-established FOSM (first-order second-moment) algorithm, on the record of 55,000 earthquakes observed in the past 110 years. The new seismic hazard analysis from a different perspective shows that the annual rate for earthquake-induced PGA to exceed the current design value (i.e., 0.23g) in two major cities in Taiwan should be relatively low, with it no greater than 0.0006 per year. Besides, the FOSM estimates were found very close to those with Monte Carlo Simulation (MCS), mainly because the skewness of the three random variables (i.e., earthquake magnitude, location, and model error) considered in the probabilistic analysis is not very large.     

Ground motion selection for scenario ruptures using the generalised conditional intensity measure (GCIM) method

In this paper, the generalised conditional intensity measure (GCIM) method is extended to ground motion selection for scenario earthquake ruptures. The selection algorithm is based on generating random realisations of the considered intensity measure (IM) distributions for a specific rupture scenario and then finding the prospective ground motions that best fit the realisations using an optimal amplitude scale factor. Using different rupture scenarios and site conditions, two important aspects of the GCIM methodology are scrutinised: (i) different weight vectors for the various IMs considered and (ii) quantifying the importance of replicate selections for ensembles with different numbers of desired ground motions. It is demonstrated that considering only spectral acceleration (SA) ordinates in the selection process, as is common in many conventional selection procedures, may result in selected motions with a biased representation for duration and cumulative ground motion effects. In contrast, considering IMs other than SA ordinates (in particular, significant duration, cumulative absolute velocity, and Arias intensity) results in ensembles with an appropriate representation of these IMs, without a practically significant effect on SA ordinates. The benefit of conducting replicate selections to obtain a suite of motions with an improved representation for the distribution of the considered IMs is demonstrated, and a minimum number of replicates are suggested for different ground motion ensemble sizes. Copyright © 2015 John Wiley & Sons, Ltd.     

Prediction of spatially distributed seismic demands in specific structures: Ground motion and structural response

The efficacy of various ground motion intensity measures (IMs) in the prediction of spatially distributed seismic demands (engineering demand parameters, (EDPs)) within a structure is investigated. This has direct implications to building‐specific seismic loss estimation, where the seismic demand on different components is dependent on the location of the component in the structure. Several common IMs are investigated in terms of their ability to predict the spatially distributed demands in a 10‐storey office building, which is measured in terms of maximum interstorey drift ratios and maximum floor accelerations. It is found that the ability of an IM to efficiently predict a specific EDP depends on the similarity between the frequency range of the ground motion that controls the IM and that of the EDP. An IMs predictability has a direct effect on the median response demands for ground motions scaled to a specified probability of exceedance from a ground motion hazard curve. All of the IMs investigated were found to be insufficient with respect to at least one of magnitude, source‐to‐site distance, or epsilon when predicting all peak interstorey drifts and peak floor accelerations in a 10‐storey reinforced concrete frame structure. Careful ground motion selection and/or seismic demand modification is therefore required to predict such a spatially distributed demands without significant bias. Copyright © 2009 John Wiley & Sons, Ltd.     

Impact of hazard‐consistent ground motion duration in structural collapse risk assessment

This study evaluates the effect of considering ground motion duration when selecting hazard‐consistent ground motions for structural collapse risk assessment. A procedure to compute source‐specific probability distributions of the durations of ground motions anticipated at a site, based on the generalized conditional intensity measure framework, is developed. Targets are computed for three sites in Western USA, located in distinct tectonic settings: Seattle, Eugene, and San Francisco. The effect of considering duration when estimating the collapse risk of a ductile reinforced concrete moment frame building, designed for a site in Seattle, is quantified by conducting multiple stripe analyses using groups of ground motions selected using different procedures. The mean annual frequency of collapse (λ_collapse) in Seattle is found to be underestimated by 29% when using typical‐duration ground motions from the PEER NGA‐West2 database. The effect of duration is even more important in sites like Eugene (λ_collapse underestimated by 59%), where the seismic hazard is dominated by large magnitude interface earthquakes, and less important in sites like San Francisco (λ_collapse underestimated by 7%), where the seismic hazard is dominated by crustal earthquakes. Ground motion selection procedures that employ causal parameters like magnitude, distance, and Vs₃₀ as surrogates for ground motion duration are also evaluated. These procedures are found to produce poor fits to the duration and response spectrum targets because of the limited number of records that satisfy typical constraints imposed on the ranges of the causal parameters. As a consequence, ground motions selected based on causal parameters are found to overestimate λ_collapse by 53%. Copyright © 2016 John Wiley & Sons, Ltd.     

Modification of stochastic ground motion models for matching target intensity measures

Stochastic ground motion models produce synthetic time‐histories by modulating a white noise sequence through functions that address spectral and temporal properties of the excitation. The resultant ground motions can be then used in simulation‐based seismic risk assessment applications. This is established by relating the parameters of the aforementioned functions to earthquake and site characteristics through predictive relationships. An important concern related to the use of these models is the fact that through current approaches in selecting these predictive relationships, compatibility to the seismic hazard is not guaranteed. This work offers a computationally efficient framework for the modification of stochastic ground motion models to match target intensity measures (IMs) for a specific site and structure of interest. This is set as an optimization problem with a dual objective. The first objective minimizes the discrepancy between the target IMs and the predictions established through the stochastic ground motion model for a chosen earthquake scenario. The second objective constraints the deviation from the model characteristics suggested by existing predictive relationships, guaranteeing that the resultant ground motions not only match the target IMs but are also compatible with regional trends. A framework leveraging kriging surrogate modeling is formulated for performing the resultant multi‐objective optimization, and different computational aspects related to this optimization are discussed in detail. The illustrative implementation shows that the proposed framework can provide ground motions with high compatibility to target IMs with small only deviation from existing predictive relationships and discusses approaches for selecting a final compromise between these two competing objectives.     

A comparison of intensity‐based demand distributions and the seismic demand hazard for seismic performance assessment

This paper compares the seismic demands obtained from an intensity‐based assessment, as conventionally considered in seismic design guidelines, with the seismic demand hazard. Intensity‐based assessments utilize the distribution of seismic demand from ground motions that have a specific value of some conditioning intensity measure, and the mean of this distribution is conventionally used in design verification. The seismic demand hazard provides the rate of exceedance of various seismic demand values and is obtained by integrating the distribution of seismic demand at multiple intensity levels with the seismic hazard curve. The seismic demand hazard is a more robust metric for quantifying seismic performance, because seismic demands from an intensity‐based assessment: (i) are not unique, with different values obtained using different conditioning intensity measures; and (ii) do not consider the possibility that demand values could be exceeded from different intensity ground motions. Empirical results, for a bridge‐foundation‐soil system, illustrate that the mean seismic demand from an intensity‐based assessment almost always underestimates the demand hazard value for the exceedance rate considered, on average by 17% and with a large variability. Furthermore, modification factors based on approximate theory are found to be unreliable. Adopting the maximum of the mean values from multiple intensity‐based assessments, with different conditional intensity measures, provides a less biased prediction of the seismic demand hazard value, but with still a large variability, and a proportional increase the required number of analyses. For an equivalent number of analyses, direct computation of the seismic demand hazard is a more logical choice and provides additional performance insight. Copyright © 2013 John Wiley & Sons, Ltd.     

Evaluation of bias on the probability of collapse from amplitude scaling using spectral-shape-matched records

Although for many years it was thought that amplitude scaling of acceleration time series to reach a target intensity did not introduce any bias in the results of nonlinear response history analyses, recent studies have showed that scaling can lead to an overestimation of deformation demands with increasing scale factors. Some studies have suggested that the bias can be explained by differences in spectral shape between the response spectra of unscaled and scaled records. On the basis of these studies, some record selection procedures assume that if records are selected using spectral-shape-matching procedures, amplitude scaling does not induce any bias on the structural response. This study evaluates if bias is introduced on lateral displacement demands and seismic collapse risk estimates even when spectral shape is carefully taken into consideration when selecting ground motions. Several single-degree-of-freedom and multiple-degree-of-freedom systems are analyzed when subjected to unscaled and scaled ground motions selected to approximately match the mean and the variance of the conditional spectrum at the target level of intensity. Results show that an explicit consideration of spectral shape is not enough to avoid a systematic overestimation of lateral displacement demands and collapse probabilities as the scale factor increases. Moreover, the bias is observed in practically all cases for systems with strength degradation and it increases with decreasing period and decreasing lateral strength relative to the strength required to remain elastic. Key reasons behind the bias are presented by evaluating input energy, causal parameters, and damaging pulse distributions in unscaled and scaled ground motion sets.     

Effects of near-fault ground motions on the nonlinear behaviour of reinforced concrete framed buildings

The design provisions of current seismic codes are generally not very accurate for assessing effects of near-fault ground motions on reinforced concrete(r.c.)spatial frames,because only far-fault ground motions are considered in the seismic codes.Strong near-fault earthquakes are characterized by long-duration(horizontal)pulses and high values of the ratio α_(PGA)of the peak value of the vertical acceleration,PGA_V,to the analogous value of the horizontal acceleration,PGA_H,which can become critical for girders and columns.In this work,six- and twelve-storey r.c.spatial frames are designed according to the provisions of the Italian seismic code,considering the horizontal seismic loads acting(besides the gravity loads)alone or in combination with the vertical ones.The nonlinear seismic analysis of the test structures is performed using a step-by-step procedure based on a two-parameter implicit integration scheme and an initial stress-like iterative procedure.A lumped plasticity model based on the Haar-Karman principle is adopted to model the inelastic behaviour of the frame members.For the numerical investigation,five near-fault ground motions with high values of the acceleration ratio α_(PGA) are considered.Moreover,following recent seismological studies,which allow the extraction of the largest(horizontal) pulse from a near-fault ground motion,five pulse-type(horizontal)ground motions are selected by comparing the original ground motion with the residual motion after the pulse has been extracted.The results of the nonlinear dynamic analysis carried out on the test structures highlighted that horizontal and vertical components of near-fault ground motions may require additional consideration in the seismic codes.     

A new approach of selecting real input ground motions for seismic design: The most unfavourable real seismic design ground motions

This paper presents a new way of selecting real input ground motions for seismic design and analysis of structures based on a comprehensive method for estimating the damage potential of ground motions, which takes into consideration of various ground motion parameters and structural seismic damage criteria in terms of strength, deformation, hysteretic energy and dual damage of Park & Ang damage index. The proposed comprehensive method fully involves the effects of the intensity, frequency content and duration of ground motions and the dynamic characteristics of structures. Then, the concept of the most unfavourable real seismic design ground motion is introduced. Based on the concept, the most unfavourable real seismic design ground motions for rock, stiff soil, medium soil and soft soil site conditions are selected in terms of three typical period ranges of structures. The selected real strong motion records are suitable for seismic analysis of important structures whose failure or collapse will be avoided at a higher level of confidence during the strong earthquake, as they can cause the greatest damage to structures and thereby result in the highest damage potential from an extended real ground motion database for a given site. In addition, this paper also presents the real input design ground motions with medium damage potential, which can be used for the seismic analysis of structures located at the area with low and moderate seismicity. The most unfavourable real seismic design ground motions are verified by analysing the seismic response of structures. It is concluded that the most unfavourable real seismic design ground motion approach can select the real ground motions that can result in the highest damage potential for a given structure and site condition, and the real ground motions can be mainly used for structures whose failure or collapse will be avoided at a higher level of confidence during the strong earthquake. Copyright © 2007 John Wiley & Sons, Ltd.     

Probabilistic seismic hazard analysis: Early history

Probabilistic seismic hazard analysis (PSHA) is the evaluation of annual frequencies of exceedence of ground motion levels (typically designated by peak ground acceleration or by spectral accelerations) at a site. The result of a PSHA is a seismic hazard curve (annual frequency of exceedence vs ground motion amplitude) or a uniform hazard spectrum (spectral amplitude vs structural period, for a fixed annual frequency of exceedence). Analyses of this type were first conceived in the 1960s and have become the basis for the seismic design of engineered facilities ranging from common buildings designed according to building codes to critical facilities such as nuclear power plants. This Historical Note traces the early history of PSHA. Copyright © 2007 John Wiley & Sons, Ltd.     

Identification of critical ground motions for seismic performance assessment of structures

A method is established to identify critical earthquake ground motions that are to be used in physical testing or subsequent advanced computational studies to enable seismic performance to be assessed. The ground motion identification procedure consists of: choosing a suitable suite of ground motions and an appropriate intensity measure; selecting a computational tool and modelling the structure accordingly; performing Incremental Dynamic Analysis on a non‐linear model of the structure; interpreting these results into 50th (median) and 90th percentile performance bounds; and identifying the critical ground motions that are close to these defining probabilistic curves at ground motion intensities corresponding to the design basis earthquake and the maximum considered earthquake. An illustrative example of the procedure is given for a reinforced concrete highway bridge pier designed to New Zealand specifications. Pseudodynamic tests and finite element based time history analyses are performed on the pier using three earthquake ground motions identified as: (i) a Design Basis Earthquake (10% probability in 50 years) with 90 percent confidence of non‐exceedance; (ii) a Maximum Considered Event (2% probability in 50 years) representing a median response; and (iii) a Maximum Considered Event representing 90 percent confidence of non‐exceedance. Copyright © 2006 John Wiley & Sons, Ltd.     

Seismic demand models for probabilistic risk analysis of near fault vertical ground motion effects on ordinary highway bridges

The influence of vertical ground motions on the seismic response of highway bridges is not very well understood. Recent studies suggest that vertical ground motions can substantially increase force and moment demands on bridge columns and girders and cannot be overlooked in seismic design of bridge structures. For an evaluation of vertical ground motion effects on the response of single‐bent two‐span highway bridges, a systematic study combining the critical engineering demand parameters (EDPs) and ground motion intensity measures (IMs) is required. Results of a parametric study examining a range of highway bridge configurations subjected to selected sets of horizontal and vertical ground motions are used to determine the structural parameters that are significantly amplified by the vertical excitations. The amplification in these parameters is modeled using simple equations that are functions of horizontal and vertical spectral accelerations at the corresponding horizontal and vertical fundamental periods of the bridge. This paper describes the derivation of seismic demand models developed for typical highway overcrossings by incorporating critical EDPs and combined effects of horizontal and vertical ground motion IMs depending on the type of the parameter and the period of the structure. These models may be used individually as risk‐based design tools to determine the probability of exceeding the critical levels of EDP for pre‐determined levels of ground shaking or may be included explicitly in probabilistic seismic risk assessments. Copyright © 2011 John Wiley & Sons, Ltd.     

Spectral shape,epsilon and record selection

Selection of earthquake ground motions is considered with the goal of accurately estimating the response of a structure at a specified ground motion intensity, as measured by spectral acceleration at the first‐mode period of the structure, Sa(T₁). Consideration is given to the magnitude, distance and epsilon (ε) values of ground motions. First, it is seen that selecting records based on their ε values is more effective than selecting records based on magnitude and distance. Second, a method is discussed for finding the conditional response spectrum of a ground motion, given a level of Sa(T₁) and its associated mean (disaggregation‐based) causal magnitude, distance and ε value. Records can then be selected to match the mean of this target spectrum, and the same benefits are achieved as when records are selected based on ε. This mean target spectrum differs from a Uniform Hazard Spectrum, and it is argued that this new spectrum is a more appropriate target for record selection. When properly selecting records based on either spectral shape or ε, the reductions in bias and variance of resulting structural response estimates are comparable to the reductions achieved by using a vector‐valued measure of earthquake intensity. Copyright © 2006 John Wiley & Sons, Ltd.     

Fully probabilistic seismic displacement analysis of spatially distributed slopes using spatially correlated vector intensity measures

Earthquake‐induced slope displacement is an important parameter for safety evaluation and earthquake design of slope systems. Traditional probabilistic seismic hazard analysis usually focuses on evaluating slope displacement at a particular location, and it is not suitable for spatially distributed slopes over a large region. This study proposes a computationally efficient framework for fully probabilistic seismic displacement analysis of spatially distributed slope systems using spatially correlated vector intensity measures (IMs). First, a spatial cross‐correlation model for three key ground motion IMs, that is, peak ground acceleration (PGA), Arias intensity, and peak ground velocity, is developed using 2686 ground motion recordings from 11 recent earthquakes. To reduce the computational cost, Monte Carlo simulation and data reduction techniques are utilized to generate spatially correlated random fields for the vector IMs. The slope displacement hazards over the region are further quantified using empirical predictive equations. Finally, an illustrative example is presented to highlight the importance of the spatial correlation and the advantage of using spatially correlated vector IMs in seismic hazard analysis of spatially distributed slopes. Copyright © 2013 John Wiley & Sons, Ltd.     

Information theory measures for the engineering validation of ground‐motion simulations

This short communication introduces a quantitative approach for the engineering validation of ground‐motion simulations based on information theory concepts and statistical hypothesis testing. Specifically, we use the Kullback‐Leibler divergence to measure the similarity of the probability distributions of recorded and simulated ground‐motion intensity measures (IMs). We demonstrate the application of the proposed validation approach to ground‐motion simulations computed by using a variety of methods, including Graves and Pitarka hybrid broadband, the deterministic composite source model, and a stochastic white noise finite‐fault model. Ground‐motion IMs, acting as proxies for the (nonlinear) seismic response of more complex engineered systems, are considered herein to validate the considered ground‐motion simulation methods. The list of considered IMs includes both spectral‐shape and duration‐related proxies, shown to be the optimal IMs in several probabilistic seismic demand models of different structural types, within the framework of performance‐based earthquake engineering. The proposed validation exercise (1) can highlight the similarities and differences between simulated and recorded ground motions for a given simulation method and/or (2) allow the ranking of the performance of alternative simulation methods. The similarities between records and simulations should provide confidence in using the simulation method for engineering applications, while the discrepancies should help in improving the tested method for the generation of synthetic records.