The Tesistan,Mexico earthquake (Mw 4.9) of 11 May 2016: seismic-tectonic environment and resonance vulnerability on buildings

To highlight the importance of small earthquakes in seismic hazard, a study of the 11 May 2016, M_w 4.9, Tesistan, Mexico earthquake is presented. Due to the close proximity of the event to the city, accelerations were considerably higher than those caused by historical severe earthquakes(6.0 M_w 8.2). This paper addresses two objectives related to the Tesistan event: the first is to estimate the focal mechanism solution in order to place the event in the context of the tectonic environment of this area. The second is focused on a vulnerability evaluation of buildings that suffered resonance. Several building′s typologies with variations in construction system and height are assessed in terms of resonance with the structural and soil periods. The results show that around the Zapopan station, strong damage is expected in intermediate to high-rise buildings(12-30 m) with moment resistant frame systems and in reinforced concrete shear walls. Masonry structures around this station may not present resonance. In contrast, in the surroundings of the Guadalajara station, all intermediate height buildings from 9 to 21 m may present resonance.     

Wavelet-Hilbert transform-based simulation of pulse-like ground motion

In nonlinear dynamic structural analysis, a suite of pulse-like ground motions is required for the performance-based design of structures near active faults. The dissimilarity in the amplitude and frequency content of the earthquake time series referred to nonstationary properties in temporal and spectral, respectively. An approach is proposed based on the nonstationary properties of the far-field records and the seismological information in an event for simulating pulse-like records. The pulse-like earthquake time history is estimated via the superposition of the residual part of the earthquake with the estimated pulse. The wavelet-based Hilbert transform is utilized to characterize the nonstationary properties, the instantaneous amplitude, and frequencies of far-field records to model residual part. The effects of near-fault and pulse are estimated based on the seismological properties of the region. The validation of the procedure is indicated by comparing simulated time-series, response spectra, and Arias intensity with recorded pulse-like records in two different earthquakes in California; the M_w 6.7 1994 Northridge and the M_w 6.5 1979 Imperial valley.

     

Prediction of the maximum credible ground motion in Singapore due to a great Sumatran subduction earthquake: the worst‐case scenario

Although Singapore is located in a low‐seismicity region, huge but infrequent Sumatran subduction earthquakes might pose structural problems to medium‐ and high‐rise buildings in the city. Based on a series of ground motion simulations of potential earthquakes that may affect Singapore, the 1833 Sumatran subduction earthquake (M_w=9.0) has been identified to be the worst‐case scenario earthquake. Bedrock motions in Singapore due to the hypothesized earthquake are simulated using an extended reflectivity method, taking into account uncertainties in source rupture process. Random rupture models, considering the uncertainties in rupture directivity, slip distribution, presence of asperities, rupture velocity and dislocation rise time, are made based on a range of seismologically possible models. The simulated bedrock motions have a very long duration of about 250 s with a predominant period between 1.8 and 2.5 s, which coincides with the natural periods of medium‐ and high‐rise buildings widely found in Singapore. The 90‐percentile horizontal peak ground acceleration is estimated to be 33 gal and the 90‐percentile horizontal spectral acceleration with 5% damping ratio is 100 gal within the predominant period range. The 90‐percentile bedrock motion would generate base shear force higher than that required by the current design code, where seismic design has yet to be considered. This has not taken into account effects of local soil response that might further amplify the bedrock motion. Copyright © 2002 John Wiley & Sons, Ltd.     

Seismic response of high plasticity clays during extreme events

Mexico City high plasticity clays exhibit a small degree of nonlinearity for shear strains as large as 0.1%, which leads to both moderate shear stiffness degradation and small to medium damping increment, even for long duration subduction strong ground motions, such as the 8.1M_w 1985Michoacan earthquake. Nonetheless, current seismic design criteria of strategic infrastructure used worldwide have striven for having larger return periods for establishing the seismic environment, considering recent large magnitude (M>8.5M_w) events. This paper presents the study of the seismic response of typical high plasticity clays found in the so-called Texcoco Lake, in the surrounding of Mexico City valley, for larger to extreme earthquakes. The shear wave velocity profile was characterized using a down-hole test. The seismic environment was established from a set of uniform hazard response spectra developed for a nearby rock outcrop for return periods of 125, 250, 475 and 2475 years. A time-domain spectral matching was used to develop acceleration time histories compatible with each uniform hazard response spectrum. Both frequency and time domain site response analyses were carried out considering each seismic scenario. Ground nonlinearities were clearly observed in the soil response during extreme ground shaken, which increases rapidly with the return period. This fact must be taken into account to avoid costly and potentially unsafe seismic designs.     

Geotechnical engineering aspects of the catastrophic Tohoku (Japan) earthquake of March 11, 2011: the review of first results

Geotechnical engineering aspects of the catastrophic earthquake, which occurred in Japan on March, 11, 2011 and called Tohoku earthquake are discussed. A review is presented of the first results obtained mainly by Japanese scientists based on records of seismic networks of Japan K-NET, Kik-net and on GPS data. The basic concepts of seismic zoning in Japan and the location of the Tohoku-oki earthquake on the seismic zoning maps are described, as well as models of the source process obtained by various authors based on teleseismic data, strong motion data, GPS data, and tsunami observations. The recorded peak accelerations and velocities and their correspondence to the current empirical attenuation curves are discussed. The records of the Tohoku earthquake made by Japanese seismic networks K-NET, Kik-net and some others represent unique seismological material and the most complete seismic database (including vertical array records) in the near-source zone of a strongest earthquake with magnitude M _w = 9. These data will be studied by seismologists all over the world for many years and, probably, they will answer many questions of geotechnical seismic engineering.     

The role of alternating outcrops of sediments and basaltic lavas on seismic urban scenario: the study case of Catania, Italy

Experimental data and numerical modelling were used to study the effect of local geology on the seismic response of the Catania area. The town extends on a marly clays bedrock and terraced deposits made up by coastal sands and alluvial conglomerates. This sedimentary substratum is deeply entrenched by paleo-valleys filled by lava flows and pyroclastics. Available borehole data and elastic parameters were used to reconstruct a geotechnical model in order to perfome 1D numerical modeling. Seismic urban scenarios were simulated considering destructive (M _w = 7.0), strong (M _w = 6.2) and moderate (M _w = 5.7) earthquakes to assess the shaking level of the different outcropping formations. For each scenario seven real accelerograms were selected from the European Strong Motion Database to assess the expected seismic input at the bedrock. PGA and spectral acceleration at different periods were obtained in the urban area through the equivalent linear numerical code EERA, and contour maps of different levels of shaking were drawn. Standard and horizontal-to-vertical spectral ratios were achieved making use of a dataset of 172 seismic events recorded at ten sites located on the main outcropping lithotypes. Spectral ratios inferred from earthquake data were compared with theoretical transfer functions. Both experimental and numerical results confirm the role of the geological and morphologic setting of Catania. Amplification of seismic motion mainly occurs in three different stratigraphic conditions: (a) sedimentary deposits mainly diffused in the south of the study area; (b) spots of soft sediments surrounded by lava flows; (c) intensely fractured and scoriaceous basaltic lavas.     

云南地区GNSS应变场时空演化与地震事件关联分析

采用GAMIT/GLOBK软件对云南境内及邻区近400个GNSS测点1999～2018年的观测数据进行解算,在各个测点时间序列和速度场的基础上,采用克里金插值方法分时段估计该区域在1999～2004年、2004～2007年、2009～2013年、2013～2015年、2015～2018年共计5个时间区域应变率场;根据区域地壳面应变率和最大剪应变率的空间变化以及相应时段之后3年内的MS≥5. 0地震事件分布特征,分析发现:绝大部分震例发生在面应变高梯度带的张压转换区和最大剪应变高值区,可见研究区各个观测时段GNSS应变率场对后期1～3年内的中强震发生区域有一定的指示意义;以2014年盈江6. 1级、鲁甸6. 5级和景谷6. 6级地震为样本,建立监视块体获取应变时间序列,分析发现:地震前三个月左右均出现震中附近短期应变趋势改变、快速增强、转折的现象,这些形变异常变化或许反映了发震区应力-应变积累在接近临界破裂状态时的非线性调整,为地震短临预测尤其是时间要素的判断提供参考。     

The moment magnitude <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> and the energy magnitude <Emphasis Type="Italic">M</Emphasis><Subscript>e</Subscript>: common roots and differences

Starting from the classical empirical magnitude-energy relationships, in this article, the derivation of the modern scales for moment magnitude M _w and energy magnitude M _e is outlined and critically discussed. The formulas for M _w and M _e calculation are presented in a way that reveals, besides the contributions of the physically defined measurement parameters seismic moment M ₀ and radiated seismic energy E _S, the role of the constants in the classical Gutenberg–Richter magnitude–energy relationship. Further, it is shown that M _w and M _e are linked via the parameter Θ = log(E _S/M ₀), and the formula for M _e can be written as M _e = M _w + (Θ + 4.7)/1.5. This relationship directly links M _e with M _w via their common scaling to classical magnitudes and, at the same time, highlights the reason why M _w and M _e can significantly differ. In fact, Θ is assumed to be constant when calculating M _w. However, variations over three to four orders of magnitude in stress drop Δσ (as well as related variations in rupture velocity V _R and seismic wave radiation efficiency η _R) are responsible for the large variability of actual Θ values of earthquakes. As a result, for the same earthquake, M _e may sometimes differ by more than one magnitude unit from M _w. Such a difference is highly relevant when assessing the actual damage potential associated with a given earthquake, because it expresses rather different static and dynamic source properties. While M _w is most appropriate for estimating the earthquake size (i.e., the product of rupture area times average displacement) and thus the potential tsunami hazard posed by strong and great earthquakes in marine environs, M _e is more suitable than M _w for assessing the potential hazard of damage due to strong ground shaking, i.e., the earthquake strength. Therefore, whenever possible, these two magnitudes should be both independently determined and jointly considered. Usually, only M _w is taken as a unified magnitude in many seismological applications (ShakeMap, seismic hazard studies, etc.) since procedures to calculate it are well developed and accepted to be stable with small uncertainty. For many reasons, procedures for E _S and M _e calculation are affected by a larger uncertainty and are currently not yet available for all global earthquakes. Thus, despite the physical importance of E _S in characterizing the seismic source, the use of M _e has been limited so far to the detriment of quicker and more complete rough estimates of both earthquake size and strength and their causal relationships. Further studies are needed to improve E _S estimations in order to allow M _e to be extensively used as an important complement to M _w in common seismological practice and its applications.     

A numerical study on response modification,overstrength, and displacement amplification factors for steel plate shear wall systems

Design recommendations for steel plate shear wall (SPSW) systems have recently been introduced into seismic provisions for steel buildings. Response modification (R), overstrength (Ω_o), and displacement amplification (C_d) factors for SPSW systems presented in design codes were based on professional experience and judgment. A numerical study has been undertaken to evaluate these factors for SPSW systems. Forty‐four unstiffened SPSW possessing different geometrical characteristics were designed based on the recommendations given in the AISC Seismic Provisions. Bay width, number of stories, story mass, and steel plate thickness were considered as the prime variables that influence the response. Twenty records were selected to include the variability in ground motion characteristics. In order to provide a detailed analysis of the post‐buckling response, three‐dimensional finite element analyses were conducted for the 44 structures subjected to the selected suite of earthquake records. For each structure and earthquake record, two analyses were conducted in which the first includes geometrical nonlinearities and the other includes both geometrical and material nonlinearities, resulting in a total of 1760 time history analyses. In this paper, the details of the design and analysis methodology are given. Based on the analysis results, response modification (R), overstrength (Ω_o), and displacement amplification (C_d) factors for SPSW systems are evaluated. Copyright © 2008 John Wiley & Sons, Ltd.     

A spatial correlation model of peak ground acceleration and response spectra based on data of the Istanbul Earthquake Rapid Response and Early Warning System

Ground motion intensity measures such as the peak ground acceleration (PGA) and the pseudo-spectral acceleration (PSA) at two sites due to the same seismic event are correlated. The spatial correlation needs to be considered when modeling ground-motion fields for seismic loss assessments, since it can have a significant influence on the statistical moments and probability distribution of aggregated seismic loss of a building portfolio.Empirical models of spatial correlation of ground motion intensity measures exist only for a few seismic regions in the world such as Japan, Taiwan and California, since for this purpose a dense observation network of earthquake ground motion is required. The Istanbul Earthquake Rapid Response and Early Warning System (IERREWS) provides one such dense array with station spacing of typically 2 km in the urban area of Istanbul. Based on the records of eight small to moderate (M_w3.5–M_w5.1) events, which occurred since 2003 in the Marmara region, we establish a model of intra-event spatial correlation for PGA and PSA up to the natural period of 1.0 s.The results indicate that the correlation coefficients of PGA and short-period PSA decay rapidly with increasing interstation distance, resulting in correlation lengths of approximately 3–4 km, while correlation lengths at longer natural periods (above 0.5 s) exceed 6 km. Finally, we implement the correlation model in a Monte Carlo simulation to evaluate economic loss in Istanbul's district Zeytinburnu due to a M_w7.2 scenario earthquake.     

Critical evaluation of strong motion in Kocaeli and Düzce (Turkey) earthquakes

Turkey was struck by two major events on August 17th and November 12th, 1999. Named Kocaeli (M_w=7.4) and Düzce (M_w=7.2) earthquakes, respectively, the two earthquakes provided the most extensive strong ground motion data set ever recorded in Turkey. The strong motion stations operated by the General Directorate of Disaster Affairs, the Kandilli Observatory and Earthquake Research Institute of Bogazici University and Istanbul Technical University have produced at least 27 strong motion records for the Kocaeli earthquake within 200 km of the fault. Kocaeli earthquake has generated six motions within 20 km of the fault adding significantly to the near-field database of ground motions for M_w>=7.0 strike–slip earthquakes. The paper discusses available strong motion data, studies their attenuation characteristics, analyses time domain, as well as spectral properties such as spectral accelerations with special emphasis on fault normal and fault parallel components and the elastic attenuation parameter, kappa. A simulation of the Kocaeli earthquake using code FINSIM is also presented.     

Water-level changes in the wells of Kamchatka at the time of the <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> 7.6, April 20, 2006 Olyutorskii earthquake

The seismic waves excited by the M _w 7.6 Olyutorskii earthquake that occurred on April 20, 2006 in the Koryak Upland gave rise to water-level changes in five wells situated in continental areas of Kamchatka at hypocentral distances of 750–1150 km. We describe the effects due to seismic waves, as well as the water-level anomalies for February–April 2006 before the earthquake. We used an original technique for the processing of water-level records based on the study of barometric and tidal water-level responses in order to estimate the volume strain in water-saturated rocks during synchronous level variations at two wells. We discuss possible mechanisms for producing anomalous water-level changes due to elastic deformation of monitored groundwater reservoirs and to crack dilatancy in the water-saturated rocks.     

The August 27, 2008, M w = 6.3 Kultuk earthquake (South Baikal): The stress-strain state of the source area from the aftershock data

We consider the results of reconstructing the stress-strain state of the Earth’s crust in South Baikal from the focal mechanism data for the Kultuk earthquake of August 27, 2008 (M _w = 6.3) and its aftershocks. The source parameters of the main shock were determined by calculating the seismic moment tensor. The focal mechanism solutions of 32 aftershocks (M _w ≥ 2.3) were obtained through the deployment of a local seismic network at South Baikal. It is found that the main shock and first aftershocks (August–September) gave rise to the activation of latitudinal fragments of the segmented near-edge fault, and the sources of the consequent aftershocks were dominated by the NW-striking planes related to the small intrabasin structures. The calculations of seismotectonic deformations based on the data on the focal mechanisms of the earthquakes show that the area of activation is dominated by the transtension regime (with deformation in the form of extension with shear). The epicentral and hypocentral fields of the aftershocks and the mechanisms of their sources reflect the complex tectonic structure of the source zone of the Kultuk earthquake, which exhibits a clear subvertical zonality of the local seismically active volume and a wedge-shaped area of crustal destruction.     

Simulation of ground motion in the Moscow region using the empirical Green’s function

Historically, the Moscow region regularly experienced rather weak but quite perceptible seismic vibrations produced by intermediate-depth earthquakes of the Vrancea zone (Romania), located at a distance of 1400 km from Moscow. The coincidence of a number of unique factors such as a slowly varying focal depth, predominant source mechanisms, weak attenuation of seismic radiation in the north-northeast direction provide favorable conditions for application of the empirical Green’s function method. Using the digital seismogram of the Vrancea Mw-5.8 earthquake recorded at the Moscow seismic station, we simulated synthetic seismograms of a scenario (expected maximum) earthquake with Mw = 8.0, by application of the empirical Green’s function method adjusted for the given conditions. The calculation procedure was verified using analog records of strong earthquakes available at the Moscow seismic station. Digital records of the Obninsk seismic station included in the Incorporated Research Institutions for Seismology (IRIS) system were used for additional control. Here, the scenario earthquake was modeled using the data on a much stronger earthquake of 1990 (MW = 6.9). It is shown that, despite a certain scatter (quite adequately assessed in the scope of the method), the ultimate estimates of expected seismic impacts are quite reliable and can be recommended for practical use.     

四川荣县MS4.7、MS4.3、MS4.9地震房屋震害特征和人员伤亡分析

根据四川荣县M_S4.7、M_S4.3、M_S4.9地震现场灾害调查资料,分析房屋震害特征和人员伤亡情况,结果表明震区房屋破坏类型主要为砖混结构、砖木结构和土木结构,其中砖木结构和土木结构受损比较严重,人员伤亡主要由房屋损坏导致;造成此现象的原因主要是该地区房屋建造年代久远,房屋结构不合理、抗震性能差,短时间地震频发造成震害累积,再加上民众防震减灾意识薄弱等。     

Source Parameters of the Deadly M<Subscript>w</Subscript> 7.6 Kashmir Earthquake of 8 October, 2005

During the last six years, National Geophysical Research Institute, Hyderabad has established a semi-permanent seismological network of 5–8 broadband seismographs and 10–20 accelerographs in the Kachchh seismic zone, Gujarat with a prime objective to monitor the continued aftershock activity of the 2001 M_w 7.7 Bhuj mainshock. The reliable and accurate broadband data for the 8 October M_w 7.6 2005 Kashmir earthquake and its aftershocks from this network as well as Hyderabad Geoscope station enabled us to estimate the group velocity dispersion characteristics and one-dimensional regional shear velocity structure of the Peninsular India. Firstly, we measure Rayleigh-and Love-wave group velocity dispersion curves in the period range of 8 to 35 sec and invert these curves to estimate the crustal and upper mantle structure below the western part of Peninsular India. Our best model suggests a two-layered crust: The upper crust is 13.8 km thick with a shear velocity (Vs) of 3.2 km/s; the corresponding values for the lower crust are 24.9 km and 3.7 km/sec. The shear velocity for the upper mantle is found to be 4.65 km/sec. Based on this structure, we perform a moment tensor (MT) inversion of the bandpass (0.05–0.02 Hz) filtered seismograms of the Kashmir earthquake. The best fit is obtained for a source located at a depth of 30 km, with a seismic moment, Mo, of 1.6 × 10²⁷ dyne-cm, and a focal mechanism with strike 19.5°, dip 42°, and rake 167°. The long-period magnitude (M_A ~ M_w) of this earthquake is estimated to be 7.31. An analysis of well-developed sPn and sSn regional crustal phases from the bandpassed (0.02–0.25 Hz) seismograms of this earthquake at four stations in Kachchh suggests a focal depth of 30.8 km.     

Characterization of the seismogenic process in the Aleutian island arc: I. Source relations of the large earthquakes of 1957, 1986, and 1996 in the Andreanof Islands

An interpretation of the type, size, and interrelations of sources is proposed for the three large Aleutian earthquakes of March 9, 1957, May 7, 1986, and June 10, 1996, which occurred in structures of the Andreanof Islands. According to our interpretation, the earthquakes were caused by steep reverse faults confined to different structural units of the southern slope of the Andreanof Islands and oriented along the strike of these structures. An E-W reverse fault that generated the largest earthquake of 1957 is located within the Aleutian Terrace and genetically appears to be associated with the development of the submarine Hawley Ridge. The western and eastern boundaries of this source are structurally well expressed by the Adak Canyon in the west (～177°W) and an abrupt change in isobaths in the east (～173°W). The character of the boundaries is reflected in the focal mechanisms. The source of the earthquake of 1957 extends for about 300 km, which agrees well with modern estimates of its magnitude (M _w = 8.6). Because the earthquake of 1957 caused, due to its high strength, seismic activation of adjacent areas of the Aleutian island arc, its aftershock zone appreciably exceeded in size the earthquake source. Reverse faults that activated the seismic sources of the earthquakes of 1986 and 1996 were located within the southern slope of the Andreanof Islands, higher than the Aleutian Terrace, outside the seismic source of the 1957 earthquake. The boundaries of these sources are also well expressed in structures and focal mechanisms. According to our estimate, the length of the 1986 earthquake source does not exceed 130–140 km, which does not contradict its magnitude (M _w = 8). The length of the 1996 earthquake source is ～100 km, which also agrees with the magnitude of the earthquake (M _w = 7.8).     

漾濞MS6.4地震中土木结构房屋抗震加固效果研究

对漾濞6.4级地震灾区开展土木结构房屋加固情况与破坏形式调查。从云南相近震级地震破坏情况的对比、加固与未加固房屋的对比、强震动记录、与烈度衰减关系对比等方面进一步分析,认为抗震加固能减轻土木结构建筑的破坏。甚至可降低其宏观烈度1～2度,对抗震加固的防震减灾效果进行评估;就抗震加固工作及烈度评估工作提出建议。