Nonlinear soil response to ground motions during different earthquakes in Nepal,to arrive at surface response spectra

Catastrophic damages reported during an earthquake include building damages, excessive ground shaking, uneven settlements and liquefaction. While most of the seismic hazard studies map the probable level of ground shaking at the bedrock level, their use in assessing the above damages is very limited until the response of the local soil is also taken into account. Determination of the local soil response needs regionally recorded ground motions, dynamic soil properties, in situ geotechnical details, etc., which most of the time are not readily available for the region under study. In the present work, the response of local soil for Nepal has been studied indirectly taking into account the surface level of ground shaking during various past as well as recent EQs observed at various locations. Based on the present analysis, a low value of amplification factor for high peak horizontal acceleration and vice versa is observed in central, western as well as southern parts of Nepal. These observations suggest nonlinear soil behavior and are in accordance with the available literature. Further, the ground motion records during 2015 Nepal EQ show maximum soil response at 0.3 s which is exactly matching with the site class C obtained from in situ data for the above locations. Based on the above observations, various correlations between the high peak horizontal acceleration and the surface spectral acceleration are proposed to obtained site specific surface response spectrum for Nepal.     

Seismic hazard and site-specific ground motion for typical ports of Gujarat

Economic importance of major ports is well known, and if ports are located in seismically active regions, then site-specific seismic hazard studies are essential to mitigate the seismic risk of the ports. Seismic design of port sites and related structures can be accomplished in three steps that include assessment of regional seismicity, geotechnical hazards, and soil structure interaction analysis. In the present study, site-specific probabilistic seismic hazard analysis is performed to identify the seismic hazard associated with four typical port sites of Gujarat state (bounded by 20°–25.5°N and 68°–75°E) of India viz. Kandla, Mundra, Hazira, and Dahej ports. The primary aim of the study is to develop consistent seismic ground motion for the structures within the four port sites for different three levels of ground shaking, i.e., operating level earthquake (72 years return period), contingency level earthquake (CLE) (475 year return period), and maximum considered earthquake (2,475 year return period). The geotechnical characterization for each port site is carried out using available geotechnical data. Shear wave velocities of the soil profile are estimated from SPT blow counts using various empirical formulae. Seismicity of the Gujarat region is modeled through delineating the 40 fault sources based on the seismotectonic setting. The Gujarat state is divided into three regions, i.e., Kachchh, Saurashtra, and Mainland Gujarat, and regional recurrence relations are assigned in the form of Gutenberg-Richter parameters in order to calculate seismic hazard associated with each port site. The horizontal component of ground acceleration for three levels of ground shaking is estimated by using different ground motion attenuation relations (GMAR) including one country-specific GMAR for Peninsular India. Uncertainty in seismic hazard computations is handled by using logic tree approach to develop uniform hazard spectra for 5% damping which are consistent with the specified three levels of ground shaking. Using recorded acceleration time history of Bhuj 2001 earthquake as the input time motion, synthetic time histories are generated to match the developed designed response spectra to study site-specific responses of port sites during different levels of ground shaking. It is observed that the Mundra and Kandla port sites are most vulnerable sites for seismic hazard as estimated CLE ground motion is in order of 0.79 and 0.48 g for Mundra and Kandla port sites, respectively. Hazira and Dahej port sites have comparatively less hazard with estimated CLE ground motion of 0.17 and 0.11 g, respectively. The ground amplification factor is observed at all sites which ranges from 1.3 to 2.0 for the frequency range of 1.0–2.7 Hz. The obtained spectral accelerations for the three levels of ground motions and obtained transfer functions for each port sites are compared with provisions made in Indian seismic code IS:1893-Part 1 (2002). The outcome of present study is recommended for further performance-based design to evaluate the seismic response of the port structures with respect to various performance levels.     

Generation of stochastic earthquake ground motion in western Saudi Arabia as a first step in development of regional ground motion prediction model

Earthquake ground motion model is an essential part of seismic hazard assessment. The model consists in several empirical ground motion prediction equations (GMPEs) that are considered to be applicable to the given region. When the recorded ground motion data are scarce, numerical modeling of ground motion based on available seismological information is widely used. We describe results of stochastic simulation of ground motion acceleration records for western Saudi Arabia. The simulation was performed using the finite fault model and considering peak ground acceleration and amplitudes of spectral acceleration at natural frequencies 0.2 and 1.0 s. Based on the parameters of the input seismological model that were accepted in similar previous studies, we analyze influence of variations in the source factor (stress drop) and in the local attenuation and amplification factors (kappa value, crustal amplification). These characteristics of the model are considered as the major contributors to the ground motion variability. The results of our work show that distribution of simulated ground motion parameters versus magnitude and distance reveals an agreement with the GMPEs recently used in seismic hazard assessment for the region. Collection of credible information about seismic source, propagation path, and site attenuation parameters using the regional ground motion database would allow constraining the seismological model and developing regional GMPEs. The stochastic simulation based on regional seismological model may be applied for generation of ground motion time histories used for development of analytical fragility curves for typical constructions in the region.     

Site-specific earthquake response study for hazard assessment in Kolkata city,India

The assessment of local site effects on seismic ground motions is of great importance in earthquake engineering practice. Several destructive earthquakes in the past have demonstrated that the amplification of ground motion and associated damage to structures due to local site conditions is a significant consideration in earthquake hazard analysis. A recent paper published in this journal highlights the hazard posed by earthquakes in the megacity of Kolkata in India due to its seismic and geological settings. The seismic hazard assessment study speculates that the deep alluvial deposit in the city may increase the seismic hazard probably due to the amplification of the seismic energies. This paper focuses on the seismic response studies of the various soil strata (i.e. for local subsurface conditions) obtained from various construction sites in the city for predicted earthquake. It is very well recognized that site response studies (a part of seismic microhazard zonation for urban areas) are the first step towards performance-based foundation design or seismic risk analysis and mitigation strategy. One of the problems for carrying out site-specific study in Kolkata is the lack of recorded strong motion data in the city. Hence, this paper outlines a methodology to carry out site-specific study, where no strong motion data or seismic data are available. The methodology uses wavelet-based spectrum compatibility approach to generate synthetic earthquake motions and equivalent linear method for seismic site response analysis. The Mega City of Kolkata has been considered to explain the methodology. Seismic hazard zonation map by the Bureau of Indian Standards classifies the City of Kolkata as moderate seismic zone (Zone III) with a zone factor 0.16. On the other hand, GSHAP(Global Seismic Hazard Assessment Program) map which is based on 10% probability of exceedance in 50 years specifies a maximum peak ground acceleration (PGA) of 1.6 m/s² (0.163 g) for this region. In the present study, the seismic response has been carried out based on GSHAP. The results of the analysis indicate the amplification of ground motion in the range of 4.46–4.82 with the fundamental period ranging from 0.81 to 1.17 s. Furthermore, the maximum spectral accelerations vary in the range of 0.78–0.95 g.     

Estimation of seismic spectral acceleration in Peninsular India

Peninsular India (PI), which lies south of 24°N latitude, has experienced several devastating earthquakes in the past. However, very few strong motion records are available for developing attenuation relations for ground acceleration, required by engineers to arrive at rational design response spectra for construction sites and cities in PI. Based on a well-known seismological model, the present paper statistically simulates ground motion in PI to arrive at an empirical relation for estimating 5% damped response spectra, as a function of magnitude and source to site distance, covering bedrock and soil conditions. The standard error in the proposed relationship is reported as a function of the frequency, for further use of the results in probabilistic seismic hazard analysis.     

Site-Specific Analysis of Strong Motion Data from the September 7, 1999 Athens,Greece Earthquake

The strong ground motion from Athens, Greece 07/09/1999 earthquake has been recorded by eighteen (18) stations, fourteen (14) within the central Athens area and four (4) at the centers of nearby towns. The ground conditions for most of the recording sites were identified, based on previous geotechnical investigations carried out in the wider area of the sites, and consequently correlated to the seismic motion characteristics. Hence, it has been possible to evaluate the accuracy of different seismological methods for site characterization and also estimate soil effects on peak ground acceleration and elastic response spectra. In addition, preliminary estimates are drawn for the seismic motion characteristics at the epicentral area, where no strong motion recordings are available. The detailed soil profiles at the recordingsites are placed in the Appendix.     

汶川地震紫坪铺面板堆石坝地震波输入研究

对国家强震动台网中心紫坪铺面板堆石坝区域台站实测主震记录以及大坝台网实测余震记录进行分析,研究主震与强余震地震动的基本特征。分别选取茂县地办、郫县走石山、成都中和这3组基岩台站实测主震地震动,紫坪铺台站2008年11月6日实测余震地震动以及按水工抗震规范人工生成地震动作为数值计算的地震动输入,对紫坪铺大坝进行三维动力有限元分析,并与实测结果进行对比。研究表明,郫县走石山与成都中和2个远场台站位于断层下盘,其实测地震动的加速度反应谱长周期（0.65 s以后）分量过于丰富,不宜作为断裂带附近紫坪铺大坝的地震动输入;紫坪铺大坝台站实测的余震地震动1 Hz附近（大坝基频）的频率成分相对较少,且持续时间较短,以至于难以激发大坝响应;对比坝顶实测地震动加速度反应谱和数值计算反应谱,建议汶川地震中紫坪铺大坝动力计算时可采用茂县地办台站实测地震动或按抗震设计规范反应谱人工生成地震动。     

Multiple source and attenuation relationships for evaluation of deterministic seismic hazard: logic tree approach considering local site effects

The most important seismic hazard parameters required to demarcate seismic zones are the peak horizontal acceleration (PHA) and spectral acceleration (SA). The two approaches for evaluation of seismic hazard are the probabilistic seismic hazard analysis and the deterministic seismic hazard analysis (DSHA). The present study evaluates the seismic hazard of the South Indian Peninsular region based on the DSHA methodology. In order to consider the epistemic uncertainties in a better manner, a logic tree approach was adopted in the evaluation of seismic hazard. Two types of seismic sources and three different attenuation relations were used in the analysis. The spatial variation of PHA (mean and 84th percentile values) and SA values for 1 Hz and 10 Hz at bedrock level (84th percentile values) for the entire study area were evaluated and the results are presented here. The surface level peak ground acceleration (PGA) values will be different from that of the bedrock level values due to the local site conditions. The PGA values at ground surface level were evaluated for four different National Earthquake Hazard Reduction Program site classes by considering the non-linear site response of different soil types. The response spectra for important cities in South India were also prepared using the deterministic approach and the results are presented in this paper.     

A review of strong motion studies in Gujarat State of western India

The present work reviews the strong motion studies done in Gujarat State of western India. Prior to the 2001 Bhuj earthquake, no strong motion instrument was in operation in Gujarat. After the earthquake, number of research institutes/universities from India and abroad deployed strong motion instruments to study aftershock activity, source dynamics, path and site effects. The strong motion recordings have enhanced the general understanding of the physics of earthquakes in the region. An attempt has been made to develop attenuation relationship for the Gujarat region from the actual ground motions recorded by the strong motion networks. The Government of Gujarat with the help from Asian Development Bank, World Bank (WB), Ministry of Science and Technology and Ministry of Earth Sciences, Government of India, has established a permanent dense network of strong motion accelerograph (SMA) all over Gujarat. In addition, the Institute of Seismological Research has been established in Gandhinagar, Gujarat, with the help of WB for carrying out seismological research. Recently, many important studies have been carried out using actual acceleration data obtained from a dense network of 54 SMA, as well as synthetic data generated using region-specific ground motion parameters. The recorded data are used to obtain region-specific ground motion parameters and ground motion prediction equation. A deterministic hazard analysis for the entire state of Gujarat has been carried out using site-specific ground motion parameters. The estimated peak ground acceleration and modified Mercalli intensity values have been used to estimate the vulnerability of the different types of buildings in 31 cities of Gujarat. As Gujarat has three distinct regions having varied geological conditions, the recorded strong motion data gave an opportunity to study the effect of geological and local-site conditions on the response spectra. This study for an intra-plate region like Gujarat is a pioneer work. Still, lots of research work need to be carried out as more and more data are available, such as development of more robust ground motion prediction equations and a 3D-velocity structure of Gujarat. Generation of shake maps in real time and a credible early earthquake warning system is need of the hour for disaster mitigation and management.     

Modeling and synthesis of strong ground motion

Success of earthquake resistant design practices critically depends on how accurately the future ground motion can be determined at a desired site. But very limited recorded data are available about ground motion in India for engineers to rely upon. To identify the needs of engineers, under such circumstances, in estimating ground motion time histories, this article presents a detailed review of literature on modeling and synthesis of strong ground motion data. In particular, modeling of seismic sources and earth medium, analytical and empirical Green’s functions approaches for ground motion simulation, stochastic models for strong motion and ground motion relations are covered. These models can be used to generate realistic near-field and far-field ground motion in regions lacking strong motion data. Numerical examples are shown for illustration by taking Kutch earthquake-2001 as a case study.     

Tectonic stress,seismicity, and seismic hazard in the southeastern Carpathians

Intermediate-depth earthquakes in the Vrancea region occur in response to stress generation due to descending lithosphere beneath the southeastern Carpathians. In this article, tectonic stress and seismicity are analyzed in the region on the basis of a vast body of observations. We show a correlation between the location of intermediate-depth earthquakes and the predicted localization of maximum shear stress in the lithosphere. A probabilistic seismic hazard assessment (PSHA) for the region is presented in terms of various ground motion parameters on the utilization of Fourier amplitude spectra used in engineering practice and risk assessment (peak ground acceleration, response spectra amplitude, and seismic intensity). We review the PSHA carried out in the region, and present new PSHA results for the eastern and southern parts of Romania. Our seismic hazard assessment is based on the information about the features of earthquake ground motion excitation, seismic wave propagation (attenuation), and site effect in the region. Spectral models and characteristics of site-response on earthquake ground motions are obtained from the regional ground motion data including several hundred records of small and large earthquakes. Results of the probabilistic seismic hazard assessment are consistent with the features of observed earthquake effects in the southeastern Carpathians and show that geological factors play an important part in the distribution of the earthquake ground motion parameters.     

Seismicity, seismotectonics and seismic hazard of Italy

A first generation of probabilistic seismic hazard maps of the Italian country are presented. They are based on seismogenic zoning deriving from a kinematic model of the structural tectonic units and on an earthquake catalogue with the foreshock and aftershock events filtered out. The following ground motion parameters have been investigated and mapped using attenuation equations based on strong-motion recordings of Italian earthquakes: peak ground acceleration and velocity; Arias intensity; strong motion duration; and the pseudovelocity and pseudoacceleration spectral values at 14 fixed frequencies both for the vertical and the largest horizontal component. A Poissonian model of earthquake occurrence is assumed as a default and the hazard maps are presented in terms of ground motion values expected to be exceeded at a 10% probability level in 50 years (return period 475 years) according to the requirement of Eurocode 8 for the seismic classification of national territories, as well as in terms of exceedance probabilities of selected ground motion values. Finally, as a tentative study, the use of hybrid methods (implementing both seismogenic zones and structures), renewal processes (including earthquake forecasting) and the influence of site effects (as the basis for the planning of earthquake scenarios) were explored.     

沙牌坝址基岩场地地震动输入参数研究

重大水利水电工程地震动输入参数必须根据专门的地震危险性分析结果来确定。目前由地震危险性分析得到的一致概率反应谱具有包络的意义,不能反映实际地震的频谱特性,输入“一致概率反应谱”可能导致地震作用偏大;拟合设计反应谱人工生成地震动加速度时程的频率非平稳性也没有得到很好解决。为了解决这些问题,得到与坝址地震危险性一致、具体地震的输入参数,结合沙牌大坝提出了一套适用于重大水利水电工程基岩场地地震动输入参数确定方法：通过以有效峰值加速度为参数的概率地震危险性计算分析,确定坝址不同超越概率下的有效峰值加速度及对坝址贡献最大的潜在震源区;在最大贡献潜在震源内利用震级空间联合分布概率最大法确定坝址设定地震,依据加速度反应谱衰减关系确定与坝址设定地震对应的设计反应谱;根据设定地震结果和时变功率谱模型参数衰减关系确定时变功率谱,将时变功率谱和最小相位谱按三角级数叠加法进行强度和频率非平稳地震加速度时程合成。在对沙牌坝址区域的地震活动性及地震构造环境分析评价的基础上,采用上述方法,得到了坝址基岩场地不同超越概率下的有效峰值加速度、设计反应谱、强度和频率非平稳地震加速度时程等地震动输入参数。     

Recent Advances in Prediction and Processing of Strong Ground Motions

We present an overview of our recent results on utilizing small earthquakes in the earthquake engineering practice. Site-specific ground motion time-histories of large earthquakes can be successfully simulated using recordings of small earthquakes which are often referred to as 'empirical Green's functions' in seismology. Another important practical problem is whether and how these observations can be used in seismic risk studies which are based on empirical attenuation relations for ground motion parameters. We study a possibility of extrapolating attenuation relations for small earthquakes, to larger magnitudes using the data from the Garner Valley downhole array in Southern California. Finally we introduce efficient ground motion processing techniques in frequency- and time-domains and apply them to site response estimation.     

Seismic performance evaluation of dam-reservoir-foundation systems to near-fault ground motions

Ground motion records obtained in recent major strong earthquakes have provided evidence that ground motions recorded near the near-fault regions differ in many cases from those observed further away from the seismic source. As the forward directivity and fling effect characteristics of the near-fault ground motions, they have the potential to cause more considerable damage to structures during an earthquake. Therefore, understanding the influence of near-fault ground motions on the performance of structures is critical to mitigate damage and perform effective response. This paper presents results of a study aimed at evaluating the effects of near-fault and far-fault ground motions on seismic performance of concrete gravity dams including dam-reservoir-foundation interaction. Koyna gravity dam is selected as a numerical application. Four different near-fault ground motion records with an apparent velocity pulse are used in the analyses. The earthquake ground motions recorded at the same site from other events that the epicenter far away from the site are employed as the far-fault ground motions. The seismic performance evaluation method based on the demand-capacity ratio, the cumulative overstress duration and the spatial extent of overstressed regions is presented. The concrete damaged plasticity model including the strain hardening or softening behavior is employed in nonlinear analyses. Nonlinear seismic damage analyses of the selected concrete dam subjected to both near-fault and far-fault ground motions are performed. The results obtained from the analyses show the effects of near-fault ground motions on seismic performance of concrete gravity dams and demonstrate the importance of considering the near-fault ground excitations.     

PSHA of Van province for performance assessment using spectrally matched strong ground motion records

Within the framework of the performance based earthquake engineering, site specific earthquake spectra for Van province has been obtained. It is noteworthy that, in probabilistic seismic hazard assessment, as a first stage data from geological studies and records from the instrumental period were compiled to make a seismic source characterization for the study region. The probabilistic seismic hazard curves were developed based on selected appropriate attenuation relationships, at rock sites, with a probability of exceedance 2, 10 and 50% in 50 yrs period. The obtained results are compared with the spectral responses proposed for seismic evaluation and retrofit of building structure in Turkish Earthquake Code (2007), section 7. The acceleration response spectrums obtained from probabilistic seismic hazard analysis are matched to adjust earthquake accelerograms recorded during the 2011 Van earthquakes by using SeismoMatch v2.0 software. The aim of this procedure is to obtain a set of reasonable earthquake input motions for the seismic evaluation of existing buildings.     

香港地区地震风险评价和设防区划

香港地区隶属于中国板内地震区中的东南沿海块缘地震带。港-九块体为晚中生代至早第三纪以来以持续稳定上升为主的块体, 块体活动性明显地低于其周边相对下沉的中新生代盆地。为此, 赋予港-九块体为最大可信震级M=5.5级潜在震源区, 而其周边中新生代断陷盆地则为最大可信震级M6.0潜在震源。从大陆地震构造成因的环境、潜在震源可信震级范围M=5.0～7.5和近源地震动饱和等三个方面的可比性, 结合中国大陆地震震源破裂尺度和地震烈度影响场, 所建立的反映中国地震构造和震源破裂及地震烈度影响场特点的PGA和反应谱地震动衰减预测关系式可用于香港地区。以50%概率时的中值对中国强震记录的对比, 本研究所提供的经验期望预测关系式, 能包络实际的资料。因此, 用此地震动衰减预测模式对香港地区地震危险性进行评估, 将会得到相对保守地震动预测值。通过香港地区基岩地震动危险性分析和计算, 参照中国大陆地震设防标准, 进行了基岩设计参数确定和区划。以年概率P=0.02、P=0.002、P=0.0004三个概率标准, 对应的地震动重复周期大约分别约为50a, 500a, 2500a的基岩PGA和反应谱, 作为香港地区基岩上构筑物和建筑物可选的基     

Soft sediments and damage pattern: a few case studies from large Indian earthquakes vis-a-vis seismic risk evaluation

India is prone to earthquake hazard; almost 65 % area falls in high to very high seismic zones, as per the seismic zoning map of the country. The Himalaya and the Indo-Gangetic plains are particularly vulnerable to high seismic hazard. Any major earthquake in Himalaya can cause severe destruction and multiple fatalities in urban centers located in the vicinity. Seismically induced ground motion amplification and soil liquefaction are the two main factors responsible for severe damage to the structures, especially, built on soft sedimentary environment. These are essentially governed by the size of earthquake, epicentral distance and geology of the area. Besides, lithology of the strata, i.e., sediment type, grain size and their distribution, thickness, lateral discontinuity and ground water depth, play an important role in determining the nature and degree of destruction. There has been significant advancement in our understanding and assessment of these two phenomena. However, data from past earthquakes provide valuable information which help in better estimation of ground motion amplification and soil liquefaction for evaluation of seismic risk in future and planning the mitigation strategies. In this paper, we present the case studies of past three large Indian earthquakes, i.e., 1803 Uttaranchal earthquake (Mw 7.5); 1934 Bihar–Nepal earthquake (Mw 8.1) and 2001 Bhuj earthquake (Mw 7.7) and discuss the role of soft sediments particularly, alluvial deposits in relation to the damage pattern due to amplified ground motions and soil liquefaction induced by the events. The results presented in the paper are mainly focused around the sites located on the river banks and experienced major destruction during these events. It is observed that the soft sedimentary sites located even far from earthquake epicenter, with low water saturation, experienced high ground motion amplification; while the sites with high saturation level have undergone soil liquefaction. We also discuss the need of intensifying studies related to ground motion amplification and soil liquefaction in India as these are the important inputs for detailed seismic hazard estimation.     

Seismic Hazard for Selected Sites in Greece: A Bayesian Estimate of Seismic Peak Ground Acceleration

A procedure for estimating maximum values of seismic peak ground accelerationat the examined site and quantiles of its probabilistic distribution in a future timeinterval of a given length is considered. The input information for the method areseismic catalog and regression relation between peak seismic acceleration at a givenpoint and magnitude and distance from the site to epicenter (seismic attenuation law).The method is based on Bayesian approach, which simply accounts for influenceof uncertainties of seismic acceleration values. The main assumptions for the method are Poissonian character of seismic events flow and distribution law of Gutenberg-Richter's type. The method is applied to seismic hazard estimation in six selected sitesin Greece.     

Spatial variation of probabilistic seismic hazard for Mumbai and surrounding region

Mumbai city, the economical capital of India, is located on the west coast of stable intra-plate continental region of Peninsular India which has an experience of significant historical earthquakes in the past. The city stood as the fourth most populous city in the world. Recent seismo-tectonic studies of this city highlighted the presence of active West coast fault and Chiplun fault beneath the Deccan basalt. In the present study, spatial variability of probabilistic seismic hazard for Mumbai region (latitudes of 18.85–19.35°N and longitudes of 72.80–73.15°E at a grid spacing of 0.05°) which includes Mumbai city, Suburban, part of Thane district and Navi Mumbai, in terms of ground motion parameters; peak horizontal acceleration and spectral acceleration at 1.0-s period for 2 and 10 % probability of exceedance in 50 years are generated. The epistemic uncertainty in hazard estimation is accounted by employing seven different ground motion prediction equations developed for worldwide shallow crustal intra-plate environments. Further, the seismic hazard results are deaggregated for Mumbai (latitude 18.94°N, longitude 72.84°E) to understand the relative contributions of earthquake sources in terms of magnitude and distance. The generated hazard maps are compared with the zoning specified by Indian seismic code (IS1893: Part 1 in Indian standard criteria for earthquake-resistant design of structures, Part 1—General provisions and buildings. Bureau of Indian Standards, New Delhi, India, 2002) for rocky site. Present results show an underestimation of potential seismic hazard in the entire study region by non-probabilistic zoning prescribed by IS1893: Part 1 with significantly higher seismic hazard values in the southern part of Navi Mumbai.