Free‐field rotations during earthquakes: Relevance on buildings

Soil rotations around horizontal axes, during an earthquake, are studied through records collected by closely spaced arrays of strong motion accelerometers. The cross power spectrum of accelerations at nearby stations has been generally utilized to describe the spatial distribution of the motion. A number of cross spectra have been obtained during the training of these arrays. To take profit of these elaborations, a mathematical relation is established between the cross power spectrum and the power spectrum of the rotation. Rotation data presented by Liu et al, concerning 52 earthquake records collected at a single station in Taiwan, are compared with rotation data computed according to our procedure. The two series of data are suitably normalized to the peak horizontal acceleration. The data are shown in function of the distance from epicentre. The same ratio, computed according to our procedure, is in good agreement with the average value of these data. Direct measurements and the present approach have lead to evaluations of rotation higher than those predicted by mathematical investigations on the basis of the wave propagation theory, for comparable circumstances. The relevance of this input motion for relatively tall structures is examined, with reference to the structural effects that the horizontal motion concurrently provides. Meaningful will be ranked those effects of the order of magnitude of 20% or higher than those implied by the horizontal excitation. For understanding the relevance on building structures, the procedure has two areas of concern: 1) the coherence implicit in the cross power spectra, which depends on the interpolation process of the original records, collected in the arrays of instruments, and 2) the relative importance of the vertical to the horizontal input motion. As to the second item, the relevance of the rotation component on structures largely depends on the relative importance of the vertical to the horizontal input motion. When the records in an area a few km from the epicenter are considered, the response spectrum of vertical motion is comparable and in some records even higher, than that of horizontal motion, over the entire range of frequencies. This has been observed as well for the 2009 earthquake event of L'Aquila, Italy, and that at the Christchurch (New Zealand) 2011. When the response spectrum of vertical motion is comparable to that of horizontal motion, the effects of rotational motions on most engineering structures can be meaningful. Copyright © 2011 John Wiley & Sons, Ltd.     

On the reliability of long‐period response spectral ordinates from digital accelerograms

Using records from co‐located broadband and digital strong motion (SM) instruments, it is first shown that the displacement waveforms obtained by double integration of the accelerogram need not be free of unrealistic baseline drift to yield reliable spectral ordinates up to at least 10 s. Secondly, to provide objective criteria for selecting reliable digital SM records for ground motion predictions at long periods, a set of synthetic accelerograms contaminated by random long‐period noise has been used, and the difference between the original accelerograms and the spurious ones in terms of response spectra has been quantified, by introducing a noise index that can be easily calculated based on the velocity waveform of the record. The results of this study suggest that high‐pass filtering the digital acceleration record from a cutoff period selected to suppress baseline drifts on the displacement waveform appears to be in most cases too conservative and unduly depletes reliable information on long‐period spectral ordinates. Copyright © 2007 John Wiley & Sons, Ltd.     

Design response spectra for Aqaba City–Jordan

Strong motion records for earthquakes epicentred in the Gulf of Aqaba during swarms in August 1993 and November 1995 are studied. These recordings are made at two sites in the City of Aqaba, one near the seashore and the other inland. The former site has deep alluvial profile with low shear-wave velocity, while the latter's soil profile is constituted of shallow stiff alluvium overlaying rock. The magnitudes of these events range from 4·8 to 6·2 on the Richter scale. Ground motion on the seashore site is much greater than the inland site. Response spectra are derived for each of the available accelerograms. Thirteen records at the inland site are statistically treated to arrive at mean response spectrum using 5 per cent damping. The derived spectra are found to be sharp with small plateau, this is associated with small peak displacement. Empirical site-dependent response spectra are also derived based on recent research and on recommendations of the UBC'97. Comparison is made between the empirical and measured spectra. These were found to be quite close for the inland site but data are not sufficient for the seashore site to make a conclusive comparison. Copyright © 1999 John Wiley & Sons, Ltd.     

Characteristics of observed strong motion accelerograms from the 2011 Lorca (Spain) Earthquake

The Lorca earthquake (southeast Spain) on May 11, 2011; Mw 5.1, and its aftershocks, have provided an important set of accelerograms recorded by the strong motion network of the Instituto Geográfico Nacional of Spain. It is particularly noticeable that the record obtained in Lorca town, very near of the fault rupture, presents a PGA value close to 0.37 g in the N30W component. This paper provides an overview of the strong motion data recorded during the Lorca seismic series, with particular attention to the accelerograms from the mainshock and foreshock and its characteristics. Due to the special circumstances of these two records, a more detailed processing has been required, in which various alternatives to adjust the baseline have been considered and analyzed. Based on this special processing, small residual displacements were obtained and reported in some of these cases. However, given the sensitivity of the process and the small obtained values, these findings should be taken with caution. Besides, response spectra have been analyzed and compared with design spectra proposed by the Spanish Seismic Code (NCSE-02) for the towns of Lorca, Alhama de Murcia and Mula. Large amplitude differences were observed in these spectra if compared to those recorded in Lorca. Also noteworthy is that the design spectra proposed for this town were exceeded by some horizontal directions of the response spectra.     

Evaluation of near-source seismic records based on damage potential parameters: Case study: Greece

The study presents the results of an investigation that identifies among the available records in Greece those that reveal near-source characteristics, using a procedure based on damage potential parameters. The findings reaffirm the opinion that earthquakes with magnitude less than M=6.0 present near-source phenomena that can cause severe damage to relatively stiff structures, common in urban areas. Spectra from selected accelerograms of near-source records are compared with the corresponding elastic spectra of the current Greek Seismic Code, (EAK 2000), the impulsive character of the Greek near-source records is ascertained and shortcomings of EAK 2000 to account for near-source effects are demonstrated.Based on records from seismic events in Greece and records from international earthquakes of small-to-moderate magnitude, the study demonstrates that there exists a near-source magnitude-distance region, where the velocity pulses have smaller amplitude and period in comparison with earthquakes of great magnitude. A simplified representation of three pulse types is, also, adopted for near-source events. It is found that the type-A pulse related to permanent displacement phenomena does not characterize Greek records. In addition a simple criterion is developed to identify the most appropriate simplified pulse type for near-source seismic events independently of magnitude.     

The study of the radiation characteristics and propagation of seismic waves in the North Caucasus by modeling the accelerograms of the recorded earthquakes

For evaluating the parameters of the vibrations of the Earth’s surface in the case of strong earthquakes, which are possible in the future, the regular patterns of the emission and propagation of seismic waves in the North Caucasus regions are investigated. The regional parameters of emission and propagation of seismic waves are evaluated by solution of the inverse problems of stochastic modeling of the accelerograms of the earthquakes, recorded by the seismic station in Sochi. The horizontal components of the strongest earthquakes (M _w ～ 3.9?5.6), that occurred in 2002–2006 within a radius of ～300 km from the seismic station, with source depths up to 60 km are modeled. For calculations of accelerograms, estimates of the quality are used, obtained earlier for this region in the form: Q(f) ～ 80 ～ f ^0.9. The parameter settings are carried out, which determine the shapes of the source spectra, the amplification of the seismic waves in the Earth’s crust, the weakening of the waves at high frequencies (κ), the parameters that determine the shape and duration of accelerograms, etc. Sufficiently good agreement of the calculated and recorded accelerograms is obtained, the regional characteristics of emission and propagation of seismic waves, which can be used for prediction of the parameters of strong motions in the North Caucasus, are evaluated; however, in the future these characteristics should be studied in more detail.     

Site-dependent design spectra and strength modification factors, based on records from Greece

Elastic and inelastic spectra are derived, based on a representative sample of acceleration records from Greece, carefully selected based on magnitude, distance and peak ground acceleration criteria, and grouped into three ground condition categories according to the 2004 Eurocode 8 (EC8) provisions. Using software developed in-house, elastic (pseudoacceleration, pseudovelocity and displacement), as well as inelastic (strength and displacement) spectra are computed for various critical damping ratios and ductility levels. After appropriate scaling, mean spectra are computed both irrespective of, as well as for each different, ground condition, and comparisons with EC8 provisions are made. As a further evaluation of the code spectra, three additional earthquake scenarios are considered representing ground-motion characteristics not reflected in the compiled dataset of records. Subsequently, modification factors for strength (q_μ) are derived from statistical analysis of constant ductility spectra, and corresponding empirical relationships, suitable for design purposes, are proposed.     

Displacement controlled rocking behaviour of rigid objects

The capacity of a gravity structure to counter seismically induced overturning can only be estimated with good accuracy using a dynamic analysis of the rotational (rocking) motion involving large displacement theory. Seismic assessment employing quasi‐static analysis can be overly conservative if the reserve capacity of this type of rocking structure to displace without overturning is not taken into account. It was revealed through dynamic testing on a shaking table that the overturning hazards of ground shaking are best represented by the peak displacement demand (PDD) parameter and that the vulnerability to overturning instability decreases with the increasing size of the object when the aspect ratio is held constant. This finding has important implications on the engineering of structures for countering moderate ground shaking in regions of low and moderate seismicity. Experimental data were validated and supplemented by computer simulations that involved generating artificial accelerograms of designated earthquake scenarios and non‐linear time‐history analyses of the overturning motions. Based on these simulations, fragility curves were constructed for estimating the probability of overturning for given levels of PDD and for different specimen dimensions. An expression was developed for estimating the level of PDD required to overturn rectangular objects of given dimensions for 5% probability of exceedance. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical simulation of ground acceleration spectra and accelerograms for engineering application

A numerical approach to the earthquake ground motion analysis is proposed for regions where no accelerograms are available. Using Haskell matrix techniques, the response spectra of a layered substratum for SV waves were calculated and then multiplied by the spectra corresponding to Brune's type pulses. The ground acceleration spectra were obtained for different angles of pulse incidence at the substratum base. The spectrum shape depends upon the substratum response and the pulse shape, while its level was related to the maximum ground acceleration corresponding to the expected maximum intensity. Transformation of the ground spectra into the time domain produced numerical accelerograms for horizontal and vertical components and for different angles of pulse incidence. Finally, a standard statistical procedure was applied to obtain the design response spectra used in engineering applications.     

Stress Release and Source Scaling of the 2010–2011 Canterbury,New Zealand Earthquake Sequence from Spectral Inversion of Ground Motion Data

The Canterbury earthquake sequence beginning with the 2010 M _W 7.2 Darfield earthquake is one of the most notable and well-recorded crustal earthquake sequences in a low-strain-rate region worldwide and as such provides a unique opportunity to better understand earthquake source physics and ground motion generation in such a tectonic setting. Ground motions during this sequence ranged up to extreme values of 2.2 g, recorded during the February 2011 M _W 6.2 event beneath the city of Christchurch. A better understanding of the seismic source signature of this sequence, in particular the stress release and its scaling with earthquake size, is crucial for future ground motion prediction and hazard assessment in Canterbury, but also of high interest for other low-to-moderate seismicity regions where high-quality records of large earthquakes are lacking. Here we present a source parameter study of more than 200 events of the Canterbury sequence, covering the magnitude range M _W 3–7.2. Source spectra were derived using a generalized spectral inversion technique and found to be well characterized by the ω ^?2 source model. We find that stress drops range between 1 and 20 MPa with a median value of 5 MPa, which is a factor of 5 larger than the median stress drop previously estimated with the same method for crustal earthquakes in much more seismically active Japan. Stress drop scaling with earthquake size is nearly self-similar, and we identify lateral variations throughout Canterbury, in particular high stress drops at the fault edges of the two major events, the M _W 7.2 Darfield and M _W 6.2 Christchurch earthquakes.     

What is a Poor Quality Strong-Motion Record?

Some accelerograms are affected by non-standard recording and digitization problems that mean they are often not used in strong-motion studies. These non-standard problems cannot be corrected by the standard processing techniques that remove low and high-frequency noise from the time-history. Records from analogue instruments are more prone to these problems but even records from digital instruments, which are becoming increasingly common, can be affected by such errors. Since all strong-motion data is valuable it is important to know whether any useful information can be obtained from accelerograms that are affected by such problems. This article examines whether strong-motion records from analogue instruments that are missing their initial part due to late triggering of the instrument and also strong-motion records from digital instruments with low A/D converter resolution can be used for response spectral studies. It is found, by simulating such errors on high-quality strong-motion records, that good response spectral ordinates can be obtained from such `poor-quality' records within the period range of most engineering interest. This revised version was published online in July 2006 with corrections to the Cover Date.     

Damage Potential of the 1999 Athens, Greece, Accelerograms

Athens experienced a devastating earthquake on September 7, 1999, centred about 10km beneath the outskirts of the city, with peak accelerations in the meizoseismal area estimated above 0.5 g. While the magnitude of the earthquake was moderate (M_S=5.9), damage and loss of life were extensive, numbering over 143 fatalities, 90 collapsed buildings, and over 100,000 rendered homeless. The most severe losses took place within 8 km of the fault, and peak rock accelerations in excess of 0.35 g were recorded 12 km from fault in downtown Athens. Local geologic effects played an important role, as demonstrated by concentrations of building damage on pockets of soft soil and near river canyons in the northwest part of the city. The earthquake was the first to severely damage Athens in over 2500 years. This study examines the engineering characteristics of 12 triaxial strong motion accelerograms recorded during the main shock. Properties investigated include peak ground acceleration and velocity, bracketed duration, Housner intensity, and response spectra. Spectral values at different orientations in the horizontal plane (“planar spectra”) are computed to infer the interaction of source directivity effects and building period, on damage potential. The role of rupture directivity close to the fault is investigated by means of idealized triangular pulses. Inelastic effects are examined using ductility spectra and sliding block analyses. Evidence is presented that low-rise structures in the area may have higher strength and longer natural period than those anticipated by building codes. The implications of the damage potential calculations for earthquake hazard assessment in Eastern North America and Northern Europe are also addressed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Near-source ground motion attenuation relationship for PGA and PSA of the vertical and horizontal components

In this paper, empirical ground-motion models for the vertical and average horizontal components of peak ground-motion and acceleration response spectra from shallow crustal earthquakes are derived using near-source database. These attenuation relationships were derived using a worldwide dataset consisted of corrected and processed accelerograms of 678 strong-motion records recorded with 60 km of the rupture plane of earthquakes between M_w 5.2 and 7.9. Ground motion models are functions of earthquake mechanism, distance from source to site, local average shear wave velocity, nonlinear soil response, sediment depth, depth-to-top of the rupture, hanging wall effects and faulting mechanism.     

Assessment of liquefaction evaluation procedures and severity index frameworks at Christchurch strong motion stations

The objective of the study presented herein is to assess three commonly used CPT-based liquefaction evaluation procedures and three liquefaction severity index frameworks using data from the 2010–2011 Canterbury earthquake sequence. Specifically, post-event field observations, ground motion recordings, and results from a recently completed extensive geotechnical site investigation programme at selected strong motion stations (SMSs) in the city of Christchurch and surrounding towns are used herein. Unlike similar studies that used data from free-field sites, accelerogram characteristics at the SMS locations can be used to assess the performance of liquefaction evaluation procedures prior to their use in the computation of surficial manifestation severity indices. Results from this study indicate that for cases with evidence of liquefaction triggering in the accelerograms, the majority of liquefaction evaluation procedures yielded correct predictions, regardless of whether surficial manifestation of liquefaction was evident or not. For cases with no evidence of liquefaction in the accelerograms (and no observed surficial evidence of liquefaction triggering), the majority of liquefaction evaluation procedures predicted liquefaction was triggered. When all cases are used to assess the performance of liquefaction severity index frameworks, a poor correlation is shown between the observed severity of liquefaction surface manifestation and the calculated severity indices. However, only using those cases where the liquefaction evaluation procedures yielded correct predictions, there is an improvement in the correlation, with the Liquefaction Severity Number (LSN) being the best performing of the frameworks investigated herein. However scatter in the relationship between the observed and calculated surficial manifestation still remains for all liquefaction severity index frameworks.     

Numbers of scaled and matched accelerograms required for inelastic dynamic analyses

Selecting, scaling and matching accelerograms are critically important to engineering design and assessment, enabling structural response to be determined with greater confidence and through fewer analyses than if unscaled accelerograms are employed. This paper considers the response of an 8‐storey multiple‐degree‐of‐freedom reinforced concrete structure to accelerograms selected, linearly scaled or spectrally matched using five different techniques. The first method consists of selecting real records on the basis of seismological characteristics, while the remaining methods make an initial selection on the basis of magnitude and spectral shape before (1) scaling to the target spectral acceleration at the initial period; (2) scaling to the target spectrum over a range of periods; (3) using wavelet adjustments to match the target spectrum and (4) using wavelet adjustments to match multiple target spectra for multiple damping ratios. The analyses indicate that the number of records required to obtain a stable estimate of the response decreases drastically as one moves through these methods. The exact number varies among damage measures and is related to the predictability of the damage measure. For measures such as peak roof and inter‐storey drift, member end rotation and the Park and Ang damage index, as few as one or two records are required to estimate the response to within ±5% (for a 64% confidence level) if matching to multiple damping ratios is conducted. Bias checks are made using predictive equations of the expected response derived from the results of 1656 nonlinear time‐domain analyses of the structure under the action of unscaled accelerograms. Copyright © 2008 John Wiley & Sons, Ltd.     

根据渐进谱的统计规律生成地震加速度时程

基于美国西部80条基岩上的近场强震记录, 采用Nakayama方法生成记录的渐进谱, 并参照Kameda方式,用统计方法建立了根据震级、距离等地震参数预测渐进谱的统计模型. 提出一种以渐进谱为目标谱的生成幅值和频率非平稳地震加速度时程的迭代方法. 由于考虑了渐进谱幅值和相位的相互影响,所生成的时程的相位也是时频非平稳的,并在相位调整中识别了相位谱增量符号以加速迭代收敛进程. 最后根据统计回归的目标渐进谱模型和本文提出的拟合目标渐进谱的方法,可生成不同震级、距离条件下的幅值和频率均非平稳的地震加速度时程.      

Earthquake accelerogram simulation with statistical law of evolutionary power spectrum

By using the technique for evolutionary power spectrum proposed by Nakayama and with reference to the Kameda formula, an evolutionary spectrum prediction model for given earthquake magnitude and distance is established based on the 80 near-source acceleration records at rock surface with large magnitude from the ground motion database of western U.S.. Then a new iteration method is developed for generation of random accelerograms non-stationary both in amplitude and frequency which are compatible with target evolutionary spectrum. The phase spectra of those simulated accelerograms are also non-stationary in time and frequency domains since the interaction between amplitude and phase angle has been considered during the generation. Furthermore, the sign of the phase spectrum increment is identified to accelerate the iteration. With the proposed statistical model for predicting evolutionary power spectra and the new method for generating compatible time history, the artificial random earthquake accelerograms non-stationary both in amplitude and frequency for certain magnitude and distance can be provided.     

Evaluating nonlinear effective stress site response analyses using records from the Canterbury earthquake sequence

A widely used one-dimensional nonlinear effective stress site response analysis program is used to model the response of potentially liquefiable soils during strong shaking. Ground motion records from six events of the 2010–2011 Canterbury earthquake sequence and the extensive site investigation data that have been obtained for the Christchurch area provide the basis for the analyses. The results of the analyses depend significantly on the input motions and soil profile characterization, so these important aspects are examined. Deconvolved Riccarton Gravel input motions were generated, because recorded rock or firm layer motions were not available. Nonlinear effective stress seismic site response analyses are shown to capture key aspects of the observed soil response through the comparison of acceleration response spectra of calculated surface motions to those of recorded surface motions; however, equivalent-linear and total stress nonlinear analyses capture these aspects as well. Biases in the computed motions compared to recorded motions were realized for some cases but they can be attributed primarily to the uncertainty in the development of the input motions used in the analyses.     

2010年新西兰MW7.0主震与2011年MW6.1余震近场强地面运动特征比较分析

2010年9月3日16时35分46秒新西兰南岛Greendale附近发生了M_W7.0地震, 震源深度约10.0 km. 2011年2月21日新西兰南岛又发生了M_W6.1地震, 为2010年M_W7.0主震后最大的一次余震, 震源深度约5.0 km, 发震断层为Christchurch南约9 km一条近东西走向逆冲的隐伏断层, 该地震造成Christchurch城内多处建筑物严重损毁. 本文分析了2010年新西兰地震事件M_W7.0主震与M_W6.1余震强地面运动的特征. 新西兰M_W6.1余震近场强地面运动整体高于M_W7.0主震. 将主震和余震的强震观测记录分别与新一代衰减关系(NGA)进行对比, 发现余震强震观测数据整体高于其震级对应的NGA. 分别选取距离主震和余震震中最近且强震观测记录最高的两个台站(GDLC台站和HVSC台站)作为参照台站, 建立动态复合震源模型(DCSM)及有限断层随机振动模型(SFFM)进行强地面运动的模拟计算, 分析两种模型的模拟结果并对比二者的优势及局限, 以便在未来工作中更好地通过模型计算强地面运动特征, 实现区域化特征快速、 实时分析及局部重点、 细致分析相结合的目标.