Probabilistic seismic hazard analysis for Bangalore

This article presents the results of probabilistic seismic hazard analysis (PSHA) for Bangalore, South India. Analyses have been carried out considering the seismotectonic parameters of the region covering a radius of 350 km keeping Bangalore as the center. Seismic hazard parameter ‘b’ has been evaluated considering the available earthquake data using (1) Gutenberg–Richter (G–R) relationship and (2) Kijko and Sellevoll (1989, 1992) method utilizing extreme and complete catalogs. The ‘b’ parameter was estimated to be 0.62 to 0.98 from G–R relation and 0.87 ± 0.03 from Kijko and Sellevoll method. The results obtained are a little higher than the ‘b’ values published earlier for southern India. Further, probabilistic seismic hazard analysis for Bangalore region has been carried out considering six seismogenic sources. From the analysis, mean annual rate of exceedance and cumulative probability hazard curve for peak ground acceleration (PGA) and spectral acceleration (Sa) have been generated. The quantified hazard values in terms of the rock level peak ground acceleration (PGA) are mapped for 10% probability of exceedance in 50 years on a grid size of 0.5 km × 0.5 km. In addition, Uniform Hazard Response Spectrum (UHRS) at rock level is also developed for the 5% damping corresponding to 10% probability of exceedance in 50 years. The peak ground acceleration (PGA) value of 0.121 g obtained from the present investigation is slightly lower (but comparable) than the PGA values obtained from the deterministic seismic hazard analysis (DSHA) for the same area. However, the PGA value obtained in the current investigation is higher than PGA values reported in the global seismic hazard assessment program (GSHAP) maps of Bhatia et al. (1999) for the shield area.     

Soil liquefaction studies at Mumbai city

Mumbai city is the economical capital of India and is situated about midway on the western coast of stable continental region of Peninsular India. Major part of the city being of reclaimed land, the soil type is of alluvium, sand, and recent conglomerate. There are some bigger water bodies within the city range. In this study, an attempt has been made to study the susceptibility of soil liquefaction using simplified empirical procedure based on number of blow counts (N values) of the soil layers from standard penetration test. The liquefaction susceptibility is quantified in terms of factor of safety along the borehole depths at available borehole locations using earthquake-induced cyclic stress on the soil and the cyclic resistance of the soil to withstand the load. The factor of safety against liquefaction is evaluated at different sites for two peak ground acceleration (PGA) levels pertaining to 10 and 2?% probability of exceedance in 50?years corresponding to uniform hazard response spectra for Mumbai city with 475- and 2,475-year return period, respectively. Contour maps are prepared that display the factor of safety at different depths for earthquake magnitude of M _w 6.5. These contour maps show the liquefaction vulnerability at different sites in the city.     

Liquefaction hazard mapping at the town of Edessa,Northern Greece

The town of Edessa is located on Northern Greece at a region that is characterized as low seismicity zone due to the fact that few moderate events of M < 6 occurred during the last century. According to the Greek Seismic Code, the expected acceleration having a 10% probability of being exceeded in 50 years is equal to 0.16g. However, an amplification of ground motion is likely to occur due the local geology that is consisted of Holocene fluvio-torrential deposits. The basic aim of this paper is to evaluate the site amplification due to geological conditions and to assess the liquefaction hazard. In order to achieve this, 1-D site response analyses were performed. The data that were employed for the construction of the numerical models have been collected from borings with standard penetrations tests (SPT) that were drilled for construction purposes. Afterward, the liquefaction potential of the subsoil layers was evaluated taking into consideration two seismic scenarios. The first scenario was based on the seismic parameters, earthquake magnitude and PGA, assigned by the Greek Seismic Code. On the second seismic model, we employed the values of acceleration, resulted from the 1-D analyses and the earthquake magnitude as it was defined by the Greek Seismic Code. In order to compile the liquefaction hazard maps, we initially estimated the liquefaction potential index (LPI) of the soil columns using the parameters provided by SPT, for both seismic loadings, and afterward we correlated these values with the proposed classification of the severity of liquefaction-induced deformations. In addition, having computed the value of probability based on the LPI, liquefaction manifestations probability maps were compiled for both scenarios. The result of this study was that liquefaction-induced ground disruptions are likely to occur at the center of the city, among the branches of Voda River, only when the amplified values of acceleration are taken into account to the computation of liquefaction potential.     

A performance-based framework for assessing liquefaction potential based on CPT data

This article describes a new performance-based approach for evaluating the return period of seismic soil liquefaction based on standard penetration test (SPT) and cone penetration test (CPT) data. The conventional liquefaction evaluation methods consider a single acceleration level and magnitude and these approaches fail to take into account the uncertainty in earthquake loading. The seismic hazard analysis based on the probabilistic method clearly shows that a particular acceleration value is being contributed by different magnitudes with varying probability. In the new method presented in this article, the entire range of ground shaking and the entire range of earthquake magnitude are considered and the liquefaction return period is evaluated based on the SPT and CPT data. This article explains the performance-based methodology for the liquefaction analysis – starting from probabilistic seismic hazard analysis (PSHA) for the evaluation of seismic hazard and the performance-based method to evaluate the liquefaction return period. A case study has been done for Bangalore, India, based on SPT data and converted CPT values. The comparison of results obtained from both the methods have been presented. In an area of 220 km² in Bangalore city, the site class was assessed based on large number of borehole data and 58 Multi-channel analysis of surface wave survey. Using the site class and peak acceleration at rock depth from PSHA, the peak ground acceleration at the ground surface was estimated using probabilistic approach. The liquefaction analysis was done based on 450 borehole data obtained in the study area. The results of CPT match well with the results obtained from similar analysis with SPT data.     

Comprehensive seismic hazard assessment of Tripura and Mizoram states

Northeast India is one of the most highly seismically active regions in the world with more than seven earthquakes on an average per year of magnitude 5.0 and above. Reliable seismic hazard assessment could provide the necessary design inputs for earthquake resistant design of structures in this region. In this study, deterministic as well as probabilistic methods have been attempted for seismic hazard assessment of Tripura and Mizoram states at bedrock level condition. An updated earthquake catalogue was collected from various national and international seismological agencies for the period from 1731 to 2011. The homogenization, declustering and data completeness analysis of events have been carried out before hazard evaluation. Seismicity parameters have been estimated using G–R relationship for each source zone. Based on the seismicity, tectonic features and fault rupture mechanism, this region was divided into six major subzones. Region specific correlations were used for magnitude conversion for homogenization of earthquake size. Ground motion equations (Atkinson and Boore 2003; Gupta 2010) were validated with the observed PGA (peak ground acceleration) values before use in the hazard evaluation. In this study, the hazard is estimated using linear sources, identified in and around the study area. Results are presented in the form of PGA using both DSHA (deterministic seismic hazard analysis) and PSHA (probabilistic seismic hazard analysis) with 2 and 10% probability of exceedance in 50 years, and spectral acceleration (T = 0. 2 s, 1.0 s) for both the states (2% probability of exceedance in 50 years). The results are important to provide inputs for planning risk reduction strategies, for developing risk acceptance criteria and financial analysis for possible damages in the study area with a comprehensive analysis and higher resolution hazard mapping.     

基于地面运动强度及标准贯入试验的上海地区砂土地震液化评价

上海市地处长江三角洲前缘,黄浦江和苏州河交汇区域,特殊的地理环境与沉积环境导致浅部砂层广泛发育。随着城市建设的不断推进,上海城市区域范围的砂土地震液化风险评价成为亟待研究的课题。文章基于上海市工程钻孔数据,结合地震地面运动加速度分布与标准贯入试验,建立区域性地震液化危险性评价模型,对上海市进行了地震液化危险性评价。研究认为当发生50年超越概率10%的地震条件下,上海市陆域面积的66.0%将不会产生地震砂土液化灾害,21.8%的陆域面积仅发生轻微液化,只有崇明、横沙、长兴三岛,黄浦江及苏州河两岸地震液化等级达到中等甚至严重,占全市陆域面积12.3%;50年超越概率2%的地震条件下,随着峰值地面运动加速度整体升高,全市范围内轻微—严重液化区域明显增多,可能发生地震液化的总面积达到全市陆域面积46.25%。上海市存在砂土地震液化的危险性,但是发生概率较低。研究认为,目前的抗震设计规范中上海市的设防烈度偏高,可能导致不必要的建设成本。同时研究中的不同超越概率下的地震液化危险性评价结果为上海市工程建设相关标准的合理化改进的提供了建议和参考。      

Application of Hidden Markov Model for avalanche danger simulations for road sectors in North-West Himalaya

Seismic hazard assessment is the key tool for rational planning, safety and design of infrastructures in seismically vulnerable regions. Gujarat is the only state in peninsular India with the maximum seismic hazard of large shallow earthquakes originating from intra-plate seismicity. Probabilistic seismic hazard assessment (PSHA) of Gujarat is carried out in this paper. Three seismogenic sources, namely Kutch, Saurashtra and Mainland Gujarat, are considered, and seismicity parameters are estimated separately for each region taking into account the completeness of the available earthquake data. Peak ground acceleration (PGA) of the horizontal component and spectral acceleration at specific periods are considered as the intensity measures. Ground motion predictive equation chosen was reported to be based on simulated ground motions and verified against the strong motion records in the study region. Results are reported for the 17 major cities at the bedrock and also for the soil sites. Apart from hazard curves, 2475 and 475 years of return periods are considered for the PGA and uniform hazard spectra (UHS). The results are compared with the present recommendations of Indian Standards. Key observations include (1) Indian Standards underpredict PGA in the entire Gujarat when the soil sites are considered and in a few cities even at the bedrock; (2) amplification of PGA (or short period hazard) on account of soil sites should be included in the Indian Standard, which is currently absent; (3) shape of the UHS indicates that a separate amplification is required at the hyperbolic portion; and (4) ratio of 2475–475 years of PGA, which is considered 2.0 in Indian Standard, should be reduced to 1.5. Time-dependent recurrence model is also included in this paper and compared with conventional PSHA. General observations include that (1) hazard may increase significantly on account of time dependency; (2) this also influences the disaggregation and in turn the selection of ground motions; and (3) time since last earthquake significantly influences the extent of the effect of time dependency.     

Assessment of seismic hazard and liquefaction potential of Gujarat based on probabilistic approaches

Gujarat is one of the fastest-growing states of India with high industrial activities coming up in major cities of the state. It is indispensable to analyse seismic hazard as the region is considered to be most seismically active in stable continental region of India. The Bhuj earthquake of 2001 has caused extensive damage in terms of causality and economic loss. In the present study, the seismic hazard of Gujarat evaluated using a probabilistic approach with the use of logic tree framework that minimizes the uncertainties in hazard assessment. The peak horizontal acceleration (PHA) and spectral acceleration (Sa) values were evaluated for 10 and 2?% probability of exceedance in 50?years. Two important geotechnical effects of earthquakes, site amplification and liquefaction, are also evaluated, considering site characterization based on site classes. The liquefaction return period for the entire state of Gujarat is evaluated using a performance-based approach. The maps of PHA and PGA values prepared in this study are very useful for seismic hazard mitigation of the region in future.     

Probabilistic evaluation of seismic soil liquefaction potential based on SPT data

The performance-based liquefaction potential analysis was carried out in the present study to estimate the liquefaction return period for Bangalore, India, through a probabilistic approach. In this approach, the entire range of peak ground acceleration (PGA) and earthquake magnitudes was used in the evaluation of liquefaction return period. The seismic hazard analysis for the study area was done using probabilistic approach to evaluate the peak horizontal acceleration at bed rock level. Based on the results of the multichannel analysis of surface wave, it was found that the study area belonged to site class D. The PGA values for the study area were evaluated for site class D by considering the local site effects. The soil resistance for the study area was characterized using the standard penetration test (SPT) values obtained from 450 boreholes. These SPT data along with the PGA values obtained from the probabilistic seismic hazard analysis were used to evaluate the liquefaction return period for the study area. The contour plot showing the spatial variation of factor of safety against liquefaction and the corrected SPT values required for preventing liquefaction for a return period of 475 years at depths of 3 and 6 m are presented in this paper. The entire process of liquefaction potential evaluation, starting from collection of earthquake data, identifying the seismic sources, evaluation of seismic hazard and the assessment of liquefaction return period were carried out, and the entire analysis was done based on the probabilistic approach.     

Ground motion estimation at Guwahati city for an Mw 8.1 earthquake in the Shillong plateau

In this paper, the ground motion at Guwahati city for an 8.1 magnitude earthquake on Oldham fault in the Shillong plateau has been estimated by stochastic finite-fault simulation method. The corresponding acceleration time histories on rock level at several sites in the epicentral region have been computed. These results are validated by comparing them with the estimates obtained from Medvedev–Sponheuer–Karnik (MSK) intensity observations of 1897 Shillong earthquake. Using the local soil parameters, the simulated rock level acceleration time history at Guwahati city is further amplified up to the ground surface by nonlinear site response analysis. The results obtained are presented in the form of peak ground acceleration (PGA) contour map. The maximum amplification for PGA over Guwahati city is as high as 2.5. Based on the simulated PGA, the liquefaction susceptibility at several locations in the city has been estimated. The results are presented in the form of contours of factor of safety against liquefaction at different depths below the ground surface. It is observed that over a large part of the Guwahati city, the factor of safety against liquefaction is less than one, indicating that the city is highly vulnerable to liquefaction in the event of this earthquake. The contour maps obtained can be used in identifying vulnerable areas and disaster mitigation.     

Probabilistic seismic hazard at the archaeological site of Gol Gumbaz in Vijayapura,south India

Probabilistic seismic hazard analysis (PSHA) is carried out for the archaeological site of Vijayapura in south India in order to obtain hazard consistent seismic input ground-motions for seismic risk assessment and design of seismic protection measures for monuments, where warranted. For this purpose the standard Cornell-McGuire approach, based on seismogenic zones with uniformly distributed seismicity is employed. The main features of this study are the usage of an updated and unified seismic catalogue based on moment magnitude, new seismogenic source models and recent ground motion prediction equations (GMPEs) in logic tree framework. Seismic hazard at the site is evaluated for level and rock site condition with 10% and 2% probabilities of exceedance in 50 years, and the corresponding peak ground accelerations (PGAs) are 0.074 and 0.142 g, respectively. In addition, the uniform hazard spectra (UHS) of the site are compared to the Indian code-defined spectrum. Comparisons are also made with results from National Disaster Management Authority (NDMA 2010), in terms of PGA and pseudo spectral accelerations (PSAs) at T = 0.2, 0.5, 1.0 and 1.25 s for 475- and 2475-yr return periods. Results of the present study are in good agreement with the PGA calculated from isoseismal map of the Killari earthquake, \({\hbox {M}}_{\mathrm{w}} = 6.4\) (1993). Disaggregation of PSHA results for the PGA and spectral acceleration (\({\hbox {S}}_{\mathrm{a}}\)) at 0.5 s, displays the controlling scenario earthquake for the study region as low to moderate magnitude with the source being at a short distance from the study site. Deterministic seismic hazard (DSHA) is also carried out by taking into account three scenario earthquakes. The UHS corresponding to 475-yr return period (RP) is used to define the target spectrum and accordingly, the spectrum-compatible natural accelerograms are selected from the suite of recorded accelerograms.     

A new probabilistic seismic hazard study of Spain

In this paper, we present a probabilistic seismic hazard analysis (PSHA) for mainland Spain that takes into account recent new results in seismicity, seismic zoning, and strong ground attenuation not considered in the latest PSHA of the Spanish Building Code. Those new input data have been obtained as a three-step project carried out in order to improve the existing hazard map for mainland Spain. We have produced a new earthquake catalogue for the area, in which the earthquakes are given in moment magnitude through specific deduced relationships for our territory based on intensity data (Mezcua et al. in Seismol Res Lett 75:75–81, 2004). In addition, we included a new seismogenetic zoning based on the recent partial zoning studies performed by different authors. Finally, as we have developed a new strong ground motion model for the area García Blanco (2009), it was considered in the hazard calculation together with other attenuations gathered from different authors using data compatible with our region. With this new data, a logic tree process is defined to quantify the epistemic uncertainty related to those parts of the process. A sensitivity test has been included in order to analyze the different models of ground motion and seismotectonic zonation used in this work. Finally, after applying a weighting scheme, a mean hazard map for PGA, based on rock type condition for 10% exceedance probability in 50 years, is presented, including 15th and 85th percentile hazard maps. The main differences with the present official building code hazard map are analyzed.     

Probabilistic seismic hazard of Pakistan, Azad-Jammu and Kashmir

The seismic hazard study for Pakistan and Azad Jammu and Kashmir has been conducted by using probabilistic approach in terms of peak ground acceleration (PGA) in m/s² and also seismic hazard response spectra for different cities. A new version of Ambraseys et al. (Bull Earthq Eng 3:1–53, 2005) ground acceleration model is used, and parameterization is based on most recent updated earthquake catalogs that consisted of 14,000 events. The threshold magnitude was fixed at M _w 4.8, but seismic zones like northern Pakistan–Tajikistan, Hindukush and northern Afghanistan–Tajikistan border had M _w 5.2. The average normalized ‘a’ and ‘b’ values for all zones are 6.15 and 0.95, respectively. Seismicity of study area was modeled, and ground motion was computed for eight frequencies (0.025, 0.1, 0.2, 0.5, 1.0, 1.5, 2.0, 2.5 s) for different annual exceedance rates of 0.02, 0.01, 0.005, 0.002 and 0.001 (return periods 50, 100, 200, 500 and 1,000 years) for stiff rocks at the gridding of 0.1° × 0.1°. Seismic hazard maps based on computed PGA for 0.02, 0.01 and 0.002 annual exceedance are prepared. These maps indicate the earthquake hazard of Pakistan and surrounding areas in the form of acceleration contour lines, which are in agreement with geological and seismotectonic characteristics of the study area. The maximum seismic hazard values are found at Muzaffarabad, Gilgit and Quetta areas.     

Seismic hazard assessment of the city of Hamedan and its vicinity, west of Iran

This article presents the results of deterministic and probabilistic seismic hazard analyses (DSHA and PSHA) of the city of Hamedan and its neighboring regions. This historical city is one of the developing cities located in the west of Iran. For this reason, the DSHA and PSHA approaches have been used for the assessment of seismic hazards and earthquake risk evaluation. To this purpose, analyses have been carried out considering the historic and instrumented earthquakes, geologic and seismotectonic parameters of the region covering a radius of 100?km, keeping Hamedan as the center. Therefore, in this research, we studied the main faults and fault zones in the study area and calculated the length and distance of faults from the center of Hamedan. In the next step, we measured the maximum credible earthquake (MCE) and peak ground acceleration (PGA) using both DSHA and PSHA approaches and utilized the various equations introduced by different researchers for this purpose. The results of DSHA approach show that the MCE-evaluated value is 7.2 Richter, which might be created by Nahavand fault activities in this region. The PGA value of 0.56?g will be obtained from Keshin fault. The results of PSHA approach show that the MCE-evaluated value is 7.6 Richter for a 0.64 probability in a 50-year period. The PGA value of 0.45?g will be obtained from Keshin fault. Seismic hazard parameters have been evaluated considering the available earthquake data using Gutenberg?CRichter relationship method. The ??a?? and ??b?? parameters were estimated 5.53 and 0.68, respectively.     

PGA distributions and seismic hazard evaluations in three cities in Taiwan

This study first presents the series of peak ground acceleration (PGA) in the three major cities in Taiwan. The PGAs are back-calculated from an earthquake catalog with the use of ground motion models. The maximums of the 84th percentile (mean?+?one standard deviation) PGA since 1900 are 1.03, 0.36, and 0.10?g, in Taipei, Taichung, and Kaohsiung, respectively. Statistical goodness-of-fit testing shows that the series of PGA follow a double-lognormal distribution. Using the verified probability distribution, a probabilistic analysis was developed in this paper, and used to evaluate probability-based seismic hazard. Accordingly, given a PGA equal to 0.5?g, the annual exceedance probabilities are 0.56, 0.46, and 0.23?% in Taipei, Taichung, and Kaohsiung, respectively; for PGA equal to 1.0?g, the probabilities become 0.18, 0.14, and 0.09?%. As a result, this analysis indicates the city in South Taiwan is associated with relatively lower seismic hazard, compared with those in Central and North Taiwan.     

Seismic hazard in Northern Algeria using spatially smoothed seismicity. Results for peak ground acceleration

Seismic hazard in terms of peak ground acceleration (PGA) has been evaluated in northern Algeria using spatially smoothed seismicity data. We present here a preliminary seismic zoning in northern Algeria as derived from the obtained results.Initially, we have compiled an earthquake catalog of the region taking data from several agencies. Afterwards, we have delimited seismic areas where the b and m_max parameters are different. Finally, by applying the methodology proposed by Frankel [Seismol. Res. Lett. 66 (1995) 8], and using four complete and Poissonian seismicity models, we are able to compute the seismic hazard maps in terms of PGA with 39.3% and 10% probability of exceedance in 50 years.A significant result of this work is the observation of mean PGA values of the order of 0.20 and 0.45 g, for return periods of 100 and 475 years, respectively, in the central area of the Tell Atlas.     

Seismic Hazard and Risk Assessment for Tulbagh,South Africa: Part I – Assessment of Seismic Hazard

In this part of our study the probabilistic seismic hazard analysis (PSHA) for Tulbagh was performed. The applied procedure is parametric and consists essentially of two steps. The first step is applicable to the area in the vicinity of Tulbagh and requires an estimation of the area-specific parameters, which, in this case, is the mean seismic activity rate, , the Gutenberg-Richter parameter, b, and the maximum regional magnitude, m_max. The second step is applicable to the Tulbagh site, and consists of parameters of distribution of amplitude of the selected ground motion parameter. The current application of the procedure provides an assessment of the PSHA in terms of peak ground acceleration (PGA) and spectral acceleration (SA). The procedure permits the combination of both historical and instrumental data. The historical part of the catalogue only contains the strongest events, whereas the complete part can be divided into several subcatalogues, each assumed complete above a specified threshold of magnitude. In the analysis, the uncertainty in the determination of the earthquake was taken into account by incorporation of the concept of `apparent magnitude'. The PSHA technique has been developed specifically for the estimation of seismic hazard at individual sites without the subjective judgement involved in the definition of seismic source zones, when the specific active faults have not been mapped or identified, and where the causes of seismicity are not well understood. The results of the hazard assessment are expressed as probabilities that specified values of PGA will be exceeded during the chosen time intervals, and similarly for the spectral accelerations. A worst case scenario sketches the possibility of a maximum PGA of 0.30g. The results of the hazard assessment can be used as input to a seismic risk assessment.     

Seismic hazard analysis of Laoag City,Northern Philippines for liquefaction potential assessment

A seismic hazard analysis was conducted in Laoag City, Northern Philippines to determine the design ground motion for liquefaction potential assessment of the area. Because the hazard analysis was done within the framework of liquefaction potential assessment, only those earthquakes with magnitude–distance combinations that are capable of generating liquefaction were considered in the study. Both probabilistic and deterministic approaches were used in the analysis. From the results of the probabilistic analysis, seismic hazard curves were generated from which the ground motion with a 10% probability of exceedance in 50years was obtained. This was then modified in consideration of the soft soil condition in the study area. Deaggregation was performed to determine the most likely earthquake to generate the said level of ground shaking.     

Seismic Hazard Assessment at Esfaraen‒Bojnurd Railway,North‒East of Iran

The objective of this study is to evaluate the seismic hazard at the Esfarayen-Bojnurd railway using the probabilistic seismic hazard assessment (PSHA) method. This method was carried out based on a recent data set to take into account the historic seismicity and updated instrumental seismicity. A homogenous earthquake catalogue was compiled and a proposed seismic sources model was presented. Attenuation equations that recently recommended by experts and developed based upon earthquake data obtained from tectonic environments similar to those in and around the studied area were weighted and used for assessment of seismic hazard in the frame of logic tree approach. Considering a grid of 1.2 × 1.2 km covering the study area, ground acceleration for every node was calculated. Hazard maps at bedrock conditions were produced for peak ground acceleration, in addition to return periods of 74, 475 and 2475 years.