A newly developed seismic microzonation model of Erbaa (Tokat, Turkey) located on seismically active eastern segment of the North Anatolian Fault Zone (NAFZ)

A methodology to model seismic microzonation maps is required in the hazard mitigation decision plans of the earthquake prone areas. The stage of disaster preparedness for new residential places is of great importance for detailed seismic microzonation models. The effects of local geological and geotechnical site conditions were considered in order to establish site characterization as the initial stage of the models in this study. Dynamic soil properties based on the empirical correlations between shear wave velocity (V _s) and standard penetration test blow counts were taken into account in order to define representative soil profiles extending down to the engineering bedrock. One-dimensional site response analyses were performed to analyze earthquake characteristics on the ground surface. The layers for soil classification, geology, depth to groundwater level, amplification, distance to fault, slope and aspect, and liquefaction-induced ground deformation potential of the study area were prepared in seismic microzonation models. The study area, Erbaa, is placed along the seismically active North Anatolian Fault Zone. Final seismic microzonation map of the study area was evaluated applying different GIS-based Multi-Criteria Decision Analysis (MCDA) techniques. Two of the MCDA techniques, simple additive weighting and analytical hierarchical process (AHP), are considered during the evaluation step of the final seismic microzonation map. The comparison is made in order to distinguish two different maps based on these MCDA techniques. Eventually, AHP-based seismic microzonation map is more preferable for the seismic design purposes in this study.     

Ground motion studies for microzonation in Iran

An important step in effectively reducing seismic risk and the vulnerability of a city located in an earthquake prone area is to conduct a ground motion microzonation study for the desired return period. The International Institute of Earthquake Engineering and Seismology (IIEES) initiated a number of seismic microzonation projects for Iran. This paper presents the steps followed by IIEES in ground motion microzonation. IIEES performs both probabilistic and deterministic seismic hazard analysis. IIEES uses his own fault map for seismotectonic studies and develops modulus and damping curves for the soils in the study area. The experience of ground motion microzonation shows that in almost all cases, the estimated 475-year peak ground acceleration (PGA) values are higher than the PGA proposed by the Iranian seismic code. Although ground motion microzonation in Iran has some shortcomings, IIEES is making new improvement. This includes development in deterministic seismic hazard analysis, two-dimensional and three-dimensional modelling of basin and topographical effects, using microtremor measurements to find shear-wave velocity profiles in high-density urban areas and providing maps for spectral acceleration in the study area.     

Lessons learned from two case histories of seismic microzonation in Italy

The prediction of the variability of the seismic ground motion in a given built-up area is considered an effective tool to plan appropriate urban development, to undertake actions on seismic risk mitigation and to understand the damage pattern caused by a strong-motion event. The procedures for studying the seismic response and the seismic microzonation of an urban area are well established; nevertheless, some controversial points still exists and are discussed here. In this paper, the selection of a reference input motion, the construction of a subsoil model and the seismic response analysis procedures are discussed in detail, based on the authors’ experience in two Italian case histories: the seismic microzonation of the city of Benevento, which was a predictive study, and the simulation of seismic response and damage distribution in the village of San Giuliano di Puglia, which was a retrospective analysis.     

Ground motion amplification due to shallow cavities in nonlinear soils

Subsurface cavities or stiff inclusions represent mechanical discontinuities for seismic waves propagating in soils. They modify the propagation pattern of seismic waves and alter soil response in correspondence to the ground level or building foundations. In the literature, different analytical and numerical solutions have been proposed to account for the effect of underground cavities or inclusions on the motion generated by P, S or R waves. In these former studies, the subsoil was assimilated to a homogeneous, isotropic and linear elastic halfspace containing one or more cavities. In the present study, the effect of subsurface cavities on ground motion amplification has been analysed accounting for soil stiffness degradation and associated damping increase with increasing levels of shear strains, a fundamental aspect of soil behaviour under earthquakes. The analysed model was inspired to a real case represented by the village of Castelnuovo (Italy), which during the 2009 Abruzzo earthquake suffered huge damage. The main shock (6 April 2009) caused the collapse of 50 % of the whole built environment. The historical centre of Castelnuovo rises on a hill. In its subsoil, there are many cavities with roofs 2–3 m below the ground level. The longitudinal NW–SE section of the hill has been investigated by 2D nonlinear site response analyses. A preliminary site response analysis was performed without modelling cavities, to identify ground motion amplification due to mere stratigraphic and topographic factors. The numerical model was later refined inserting: (1) a single cavity below the hilltop, (2) multiple cavities placed below the ground surface of the hill and (3) multiple cavities filled with concrete (inclusions). The performed study highlights the important role exerted by underground cavities on the ground motion computed at the hill surface. This effect should be properly considered for both microzonation studies and the correct determination of the seismic actions on specific buildings.     

Analysis of soil effects and distribution of damage from the Pyrgos 1993 (Greece) earthquake

Summary On 26 March, 1993, a moderate magnitude earthquake (M _s=5.5) occurred at 3 km epicentral distance from the town of Pyrgos, in Southern Greece, causing extensive damage to masonry houses. To explain the variability of seismic intensity over the town and to propose measures against future seismic activity, a microzonation study was undertaken which combined geological, geophysical and geotechnical investigations, site specific analyses of seismic ground response and detailed recording of structural damage. The analytical predictions of ground response are correlated to soil conditions and then used to identify (micro-)zones of sites with similar seismic response. Furthermore, they are compared to quantitative estimates of damage distribution over the town. It is concluded that the peak ground acceleration, normalized against the input peak seismic acceleration, is a function of the local soil conditions as well as the seismic excitation characteristics. Hence, it cannot be defined uniquely at a site, without reference to the seismic excitation. However, the normalized peak ground velocity and the acceleration response spectra are mainly functions of the soil conditions and can be used as criteria for the practical definition of (micro-)zones. The distribution of damage in various parts of the town is at least partially attributable to local soil effects. The small epicentral distance of the earthquake, connected with the direction of the fault rupture, as well as the quality and techniques of construction, are additional factors that may have influenced the extent and distribution of damage.     

Correlation between ground motion parameters and lining damage indices for mountain tunnels

In order to study the interdependencies between the overall lining damage indices for mountain tunnels and ground motion parameters, a total of 89 ground motion records were selected for analysis. These records were divided into three groups: (1) near-field ground motion with velocity pulses, (2) near-field ground motion without velocity pulses, and (3) far-field ground motion. Calculations of Pearson’s linear correlation coefficient have been used to clarify the grade of interrelationships between ground motion parameters. Further, nonlinear dynamic analyses were carried out on mountain tunnels to determine overall lining damage indices under seismic actions. Based on numerical results, linear correlation coefficients between ground motion parameters and overall lining damage indices for mountain tunnels were then calculated and analyzed. Overall lining damage indices were found to be highly correlated to the velocity-related seismic parameters but not so well correlated to spectra-related parameters.     

巨厚沉积土夹火山岩场地非线性地震反应特性

地震基岩深度和土体动力本构模型的选取对核岛场地地震效应评价结果的合理性具有重要影响。以拟建某沿海核电厂深度470 m沉积土夹火山岩层场地的3个钻孔剖面为研究对象,采用一维等效线性波传分析（ELA）法、基于Matasovic本构模型和Davidenkov-Chen-Zhao（DCZ）本构模型的一维非线性分析（NLA）法,选取不同剪切波速的5个岩土层作为地震基岩,研究了输入地震动特性、地震基岩深度和土体动力本构模型的选取对巨厚沉积土夹火山岩层场地非线性地震反应特性的影响。结果表明：（1）以浅层硬岩夹层或深部土层作为地震基岩,NLA法计算的5%阻尼比的地表谱加速度SA的短周期部分较之ELA法的计算值大,但两者计算的地表SA谱的长周期部分几乎一致;（2）基于Matasovic模型和DCZ模型的NLA法计算的地表SA谱谱形和峰值加速度随深度的变化趋势基本一致;（3）从NLA法计算的地表峰值加速度和累积绝对速度而言,以剪切波速约2 500 m/s的浅层硬岩夹层作为地震基岩是适宜的。     

Nonlinear Analysis of Local Site Effects on Seismic Ground Response in the Bam Earthquake

The influence of local geologic and soil conditions on the intensity of ground shaking is addressed in this study. The amplification of the ground motion due to local site effects resulted in severe damage to dwellings in the Bam area during the 2003 Bam Earthquake. A unique set of strong motion acceleration recordings was obtained at the Bam accelerograph station. Although the highest peak ground acceleration recorded was the vertical component (nearly 1 g), the longitudinal component (fault-parallel motion) clearly had the largest maximum velocity as well as maximum ground displacement. Subsurface geotechnical and geophysical (down-hole) data in two different sites have been obtained and used to estimate the local site condition on earthquake ground motion in the area. The ground response analyses have been conducted considering the nonlinear behavior of the soil deposits using both equivalent linear and nonlinear approaches. The fully nonlinear method embodied in FLAC was used to evaluate the nonlinear soil properties on earthquake wave propagation through the soil layer, and compare with the response from the equivalent linear approach. It is shown that thick alluvium deposits amplified the ground motion and resulted in significant damage in residential buildings in the earthquake stricken region. The comparison of results indicated similar response spectra of the motions for both equivalent and nonlinear analyses, showing peaks in the period range of 0.3–1.5 s. However, the amplification levels of nonlinear analysis were less than the equivalent linear method especially in long periods. The observed response spectra are shown to be above the NEHRP building code design requirements, especially at high frequencies.     

Deterministic seismic microzonation of Kolkata city

This paper presents the deterministic seismic microzonation of densely populated Kolkata city situated on the world’s largest delta island with very soft and thick soil deposit in the surficial layers. A fourth-order accurate staggered-grid finite-difference algorithm for SH-wave propagation simulation in visco-elastic medium is used for the linear computation of ground motion amplifications in sedimentary deposit. Different maps such as for fundamental frequency (F ₀), peak ground acceleration (PGA), peak ground velocity, and peak ground displacement are developed for variety of end-users communities, including structural and geotechnical engineers for performance-based designs, building officials, emergency managers, land-use planners, private businesses, and the general public. The scenario of simulated amplification factors in the different frequency bands revealed that the Kolkata city is very much prone to severe damage even during a moderate earthquake and very selective damage may occur at some of the localities during local and distant earthquakes. The deterministically predicted PGA at bedrock level is 0.0844 g and the maximum PGA predicted at the free surface is 0.6 g in Kolkata city due to maximum credible earthquake (M _w = 5.4) associated with Eocene Hinge Zone at a depth of 36 km. The seismic microzonation of Kolkata city reveals that the Nager Bazar and Nimtala areas are the safest regions with earthquake point of view.     

Dynamic soil properties for microzonation of Delhi,India

Delhi, the capital of India, has experienced mild seismic shaking during several earthquakes in the past. The large variations of depth to bedrock and ground water table coupled with different soil types at different locations of Delhi necessitate a seismic microzonation study. Dynamic soil properties such as shear wave velocity, modulus reduction and damping characteristics of local soils are the basic and essential input parameters for conducting even a preliminary ground response analysis which is an essential input in microzonation studies. Shear wave velocity is not measured routinely due to its high cost and lack of the required expertise. Several researchers in the past developed correlations between shear wave velocity (V _s ) and routinely measured N values. In the present study, shear wave velocity profiles measured in the field at more than 80 borehole locations to a depth of about 20 to 32m using Spectral Analysis of Surface Waves (SASW) are presented and correlations between shear wave velocity and N values are also presented for use by engineers and designers. Results of strain and stress controlled cyclic triaxial tests on remoulded samples of sand-silt mixtures in the high strain range are used for generating the modulus reduction and damping curves and are compared with the well-known curves in the literature. The results presented in this article can be used for microzonation studies as well as site specific ground response analyses at Delhi.     

Multilinear Regression Analysis for Seismic Response and Engineering Properties of Liquefiable Coromandel Coastal Soil Deposits

This paper investigates the seismic characteristics and geotechnical properties with respect to the liquefaction potential of the deposits in the Coromandel coastal line of Nagapattinam town, Tamilnadu, India. A series of field tests were conducted using standard penetration test, cone penetration test and plate load test. Laboratory tests were conducted on the collected samples. From the results, a microzonation map was developed for the liquefaction potential and settlements. Some multilinear regression models between permeability, fines content, relative density, coefficient of curvature, coefficient of uniformity, mean particle size, factor of safety, settlement, standard penetration test values and cone penetration test values were developed. The shear wave velocity and shear modulus were calculated from the field penetration tests and correlations between the normalized values of peak ground accelerations, velocities and displacements, which were obtained from the equivalent linear ground motion analysis using SHAKE software, with other parameters of soil. From the results, it was found that at some of the areas are vulnerable to high amplification of waves even for small earthquake.     

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.     

Local site effect of a clay site in Shiraz based on seismic hazard of Shiraz Plain

Propagation of seismic waves through soil layers would drastically change the frequency content and amplitude-based features of ground motions at the surface. These alterations are known as seismic site effects. Computation of site effects of high-populated areas such as large cities is of great importance (e.g., it is used in development of seismic microzonation of a region). Shiraz is one of the most populous cities of Iran and is located in a high seismic hazardous region. A representative clay site in this city is selected to assess local site effects. The time series and random vibration theory procedure in the frequency domain are implemented to analyze the aforementioned site. Furthermore, the nonlinear dynamic soil behavior is simulated by the equivalent linear method and the nonlinear method via DEEPSOIL program. Three types of soil column uncertainties such as shear wave velocity, modulus reduction, and damping ratio of soil layers as well as depth of underlying rock half-space (D _bed) are considered herein. The mean amplification and standard deviation of natural logarithm of amplification factors are computed for a variety of analysis types. The results of the current study show that the computed mean and standard deviation of amplification factor in ln units by considering only V _S uncertainty are in good agreement with the corresponding ones by considering V _S and modulus reduction and damping ratio variabilities simultaneously for the studied site. Furthermore, it seems that the effect of bedrock depth in definition of spectral shapes of the Iranian seismic building code should be taken into account.     

Earthquake hazard in Northeast India — A seismic microzonation approach with typical case studies from Sikkim Himalaya and Guwahati city

A comprehensive analytical as well as numerical treatment of seismological, geological, geomorphological and geotechnical concepts has been implemented through microzonation projects in the northeast Indian provinces of Sikkim Himalaya and Guwahati city, representing cases of contrasting geological backgrounds — a hilly terrain and a predominantly alluvial basin respectively. The estimated maximum earthquakes in the underlying seismic source zones, demarcated in the broad northeast Indian region, implicates scenario earthquakes of M _W 8.3 and 8.7 to the respective study regions for deterministic seismic hazard assessments. The microzonation approach as undertaken in the present analyses involves multi-criteria seismic hazard evaluation through thematic integration of contributing factors. The geomorphological themes for Sikkim Himalaya include surface geology, soil cover, slope, rock outcrop and landslide integrated to achieve geological hazard distribution. Seismological themes, namely surface consistent peak ground acceleration and predominant frequency were, thereafter, overlaid on and added with the geological hazard distribution to obtain the seismic hazard microzonation map of the Sikkim Himalaya. On the other hand, the microzonation study of Guwahati city accounts for eight themes — geological and geomorphological, basement or bedrock, landuse, landslide, factor of safety for soil stability, shear wave velocity, predominant frequency, and surface consistent peak ground acceleration. The five broad qualitative hazard classifications — ‘low’, ‘moderate’, ‘high’, ‘moderate high’ and ‘very high’ could be applied in both the cases, albeit with different implications to peak ground acceleration variations. These developed hazard maps offer better representation of the local specific seismic hazard variation in the terrain.     

Seismic microzonation of a nuclear power plant site with detailed geotechnical,geophysical and site effect studies

This paper highlights the seismic microzonation carried out for a nuclear power plant site. Nuclear power plants are considered to be one of the most important and critical structures designed to withstand all natural disasters. Seismic microzonation is a process of demarcating a region into individual areas having different levels of various seismic hazards. This will help in identifying regions having high seismic hazard which is vital for engineering design and land-use planning. The main objective of this paper is to carry out the seismic microzonation of a nuclear power plant site situated in the east coast of South India, based on the spatial distribution of the hazard index value. The hazard index represents the consolidated effect of all major earthquake hazards and hazard influencing parameters. The present work will provide new directions for assessing the seismic hazards of new power plant sites in the country. Major seismic hazards considered for the evaluation of the hazard index are (1) intensity of ground shaking at bedrock, (2) site amplification, (3) liquefaction potential and (4) the predominant frequency of the earthquake motion at the surface. The intensity of ground shaking in terms of peak horizontal acceleration (PHA) was estimated for the study area using both deterministic and probabilistic approaches with logic tree methodology. The site characterization of the study area has been carried out using the multichannel analysis of surface waves test and available borehole data. One-dimensional ground response analysis was carried out at major locations within the study area for evaluating PHA and spectral accelerations at the ground surface. Based on the standard penetration test data, deterministic as well as probabilistic liquefaction hazard analysis has been carried out for the entire study area. Finally, all the major earthquake hazards estimated above, and other significant parameters representing local geology were integrated using the analytic hierarchy process and hazard index map for the study area was prepared. Maps showing the spatial variation of seismic hazards (intensity of ground shaking, liquefaction potential and predominant frequency) and hazard index are presented in this work.     

Estimating the vulnerability and loss functions of residential buildings

This paper is concerned with an investigation of the damage to residential buildings in two areas within Gilan and Zanjan provinces, Iran, caused by the Manjil-Rudbar earthquake of 20 June 1990. A statistical correlation between the observed ground motion and the damage to the residential buildings is derived for overall damaged buildings and expressed as the vulnerability function. The loss function is calculated by combining the seismic hazard with the vulnerability function.The study of vulnerability and annual seismic hazard shows that the specific annual risk for the range of motion of 0.18 to 0.5 g is equal to 0.02. This indicates that the specific risk for semi-engineered residential buildings with a lifetime of 20 years is about 33%. This study also shows that in large cities, such as Tehran, located in seismic areas, the extent of damage according to the vulnerability function will be 45 and 70% for expected maximum accelerations of 0.3 and 0.4 g, respectively.     

Proposed seismic hazard maps of Sumatra and Java islands and microzonation study of Jakarta city,Indonesia

This paper presents the development of spectral hazard maps for Sumatra and Java islands, Indonesia and microzonation study for Jakarta city. The purpose of this study is to propose a revision of the seismic hazard map in Indonesian Seismic Code SNI 03-1726-2002. Some improvements in seismic hazard analysis were implemented in the analysis by considering the recent seismic activities around Java and Sumatra. The seismic hazard analysis was carried out using 3-dimension (3-D) seismic source models (fault source model) using the latest research works regarding the tectonic setting of Sumatra and Java. Two hazard levels were analysed for representing 10% and 2% probability of exceedance (PE) in 50 years ground motions for Sumatra and Java. Peak ground acceleration contour maps for those two hazard levels and two additional macrozonation maps for 10% PE in 50 years were produced during this research. These two additional maps represent short period (0.2 s) and long-period (1.0 s) spectra values at the bedrock. Microzonation study is performed in order to obtain ground motion parameters such as acceleration, amplification factor and response spectra at the surface of Jakarta. The analyses were carried out using nonlinear approach. The results were used to develop contour of acceleration at the surface of Jakarta. Finally, the design response spectra for structural design purposes are proposed in this study.     

Liquefaction potential evaluation for subsurface soil layers of Delhi region

The study area Delhi is second most populous city and third largest urban area in the world. Though the area lies in seismic high damage risk zone, number of high rise building and construction of mega structure at several sites of the city increase rapidly. In this study field Standard Penetration Test (SPT) values of soil collected from 750 boreholes data were analyzed to identify liquefiable sub-surface soil layers. Finally, liquefaction susceptible sub-surface maps of the region at various depth (20 m, 15 m, 12 m, 9 m, 6 m and 3 m) from ground level is prepared. The outcome of this study will be useful input for preliminary foundation and designing of earthquake resistant high rise building and seismic microzonation studies of Delhi.     

Seismic Hazard Analysis for the Bangalore Region

Indian peninsular shield, which was once considered to be seismically stable, is experiencing many earthquakes recently. As part of the national level microzonation programme, Department of Science and Technology, Govt. of India has initiated microzonation of greater Bangalore region. The seismic hazard analysis of Bangalore region is carried out as part of this project. The paper presents the determination of maximum credible earthquake (MCE) and generation of synthetic acceleration time history plot for the Bangalore region. MCE has been determined by considering the regional seismotectonic activity in about 350 km radius around Bangalore city. The seismotectonic map has been prepared by considering the faults, lineaments, shear zones in the area and historic earthquake events of more than 150 events. Shortest distance from the Bangalore to the different sources is measured and then peak ground acceleration (PGA) is calculated for the different source and moment magnitude. Maximum credible earthquake found in terms of moment magnitude is 5.1 with PGA value of 0.146 g at city centre with assuming the hypo central distance of 15.88 km from the focal point. Also, correlations for the fault length with historic earthquake in terms of moment magnitude, yields (taking the rupture fault length as 5% of the total fault length) a PGA value of 0.159 g. Acceleration time history (ground motion) and a response acceleration spectrum for the corresponding magnitude has been generated using synthetic earthquake model considering the regional seismotectonic parameters. The maximum spectral acceleration obtained is 0.332 g for predominant period of 0.06 s. The PGA value and synthetic earthquake ground motion data from the identified vulnerable source using seismotectonic map will be useful for the PGA mapping and microzonation of the area.     

场地条件对地震动和震害影响的研究进展与建议

场地条件对地震动和震害有重要影响。本文简要回顾了场地条件对地震动和震害影响的研究历史,总结和评述了这一领域的研究进展。重点对地形地貌、岩土类型、覆盖层厚度、土层结构、地下水、岩土动力性质及物理地质现象等对地震动和震害的影响进行了总结和评述。提炼和概括了目前场地条件对地震动和震害影响研究中比较一致的几点认识。在此基础上,提出了这一领域当前应进一步研究的问题,主要包括地震现场资料调查、局部地形几何尺度、场地分类、覆盖层安全厚度界限、近断层地震动参数分布规律、活断层避让距离,以及强震资料积累等关键问题。