International Journal of Modern Trends in Engineering and Research www.ijmter.com e-issn No.:2349-9745, Date: 28-30 April, 2016 A REVIEW ON DEVELOPMENT OF SEISMIC HAZARD ANALYSIS OF INDIA Shivamanth A 1, Shama P S 2, Mayank K. Desai 3, Chandresh H. Solanki 4 ( 1 Applied Mechanics Department, SVNIT, Surat-395007, Email: shiva05cv@gmail.com ) ( 2 Department of Computer Science & Engineering, CUK, Kalaburagi 585367, Emai : shm.san1@gmail.com) ( 3 Applied Mechanics Department, SVNIT, Surat-395007, Email : mkd@amd.svnit.ac.in) ( 4 Applied Mechanics Department, SVNIT, Surat-395007, Email : chs@amd.svnit.ac.in) ABSTRACT The present report highlights seismic vulnerability is closely dependent on the social and economic condition of the population. A brief review of past Probabilistic seismic hazard analysis (PSHA), Deterministic Seismic Hazard Analysis (DSHA) and site characterization of efforts for estimating seismic hazard in India Seismic hazard analysis (SHA) involves the quantitative estimation of ground shaking hazards at a particular site. Seismic hazard studies are needed for the preparation of earthquake loading regulations and for various earthquake risk management purposes. Probabilistic seismic hazard analysis (PSHA) is a technique for estimating the annual rate of exceedance of a specified ground motion at a site due to known and suspected earthquake sources. Probabilistic analysis provides a framework in which uncertainties in the size, location, and rate of recurrence and effects of earthquakes to be explicitly considered in the evaluation of seismic hazard. Well established probabilistic analysis procedure is adopted to compute the prevalent hazard in terms of peak ground acceleration (PGA), short period and long period spectral accelerations for different return periods. Keywords Earthquake PSHA, DSHA, PGA and spectral accelerations. I. GENERAL INTRODUCTION The earthquake is well-known fact that occurs without warning and involves violent shaking of the ground and everything over it. It s sudden release of accumulated strain energy of the moving lithospheric or crustal plates. The earthquake hazard can be reduced by earthquake-resistant design which may require many precautions be taken to prevent the structures from most possible types of failures. The main goal of the earthquake-resistant design is to produce a structure or facility that can withstand a certain level of shaking without excessive damage [12]. Sufficient knowledge for ground conditions is important for analyses, design and construction of projects and also delays of the project, failures etc., to the soil investigations. Hence, includes soil investigation for the part of the design process. Site characterization is one of the processes of a collection of information, consideration of data, assessment and depiction through maps without which the hazards in the below ground the site cannot be known. Deterministic Seismic Hazard Analysis (DSHA) and another type of seismic hazard analysis is known as Probabilistic Seismic Hazard Analysis (PSHA). In this case, we considered the uncertainties involved in the size of the earthquake, the location of the earthquake, and the rate of recurrence that is in the interval through which earthquake repeats or reoccur that of the earthquake. Considering uncertainties involved in all these parameters, then the Probabilistic Seismic Hazard Analysis is done, it was initially developed by [9]. @IJMTER-2016, All rights Reserved
II. LITERATURE REVIEW Chandra (1977) attributed the interpolate seismicity by using different processes and also observed that some of the past earthquakes in Peninsular India. In this case, during earthquakes they were felt over a much larger area than one would expect earthquakes of equivalent magnitude to be felt in most other parts of the world. Researcher has been studied the seismicity of Peninsular India from a detailed consideration of the historical as well as recent earthquake data, and a catalogue of earthquakes. The average focal depths for these earthquakes are reported to be within the uppercrustal (0-12 km) layers complying with the depth of the Moho varying from 34 to 41 km beneath the region of southern India [18]; [13]. S. Akkar and J. J. Bommer (2010) in this article derive new equations using the strong-motion database for the seismically active areas of Europe and the Middle East. A total of 532 strong-motion accelerograms recorded at distances of up to 100 km from 131 earthquakes with moment magnitudes ranging from M 5 to 7.6 are used to obtain equations for both the larger and the geometric mean of the horizontal components Figure. 1 compares the predictions for earthquakes of surface-wave magnitude (Ms) 5.5 and 7.0 obtained from the Greek equation of [19] and the European equation of [27]. Das et al (2006) had prepared Seismic hazard maps for Northeast India based on uniform hazard response spectra for stiff sites. Through the use 261 recorded accelerograms had developed a new attenuation model for pseudo-spectral velocity. In this approach, the entire region is divided into a fine grid of around 2,500 nodes, and hazard level is predicted for each of the nodes by assuming it to be dependent on the earthquake events surrounded by 300 km radius area. Peterson et al. (2004) Gujarat state (surrounded by Latitude 20 0 25.5 0 N and Longitude 68 0 75 0 E) is seismically one of the most active regions in India. There are three fault source models were adapted to the seismic hazard sensitivity test for the Kutch region of Gujarat, India. Developed Seismic hazard maps were for each of the 3 models for 2% and 10% probability of exceedance in 50 years. The PGA s for 2% probability in 50 years was observed to be greater than 1g and for 10% probability in 50 years the was found to be between 0.2g and 0.7g. Anbazhagan et al. (2009) carried out probabilistic Seismic Hazard Analysis of Bangalore city taking into consideration the seismotectonic features within radius of 350 km around Bangalore city. The quantified hazard values in terms of the rock level peak ground acceleration (PGA) are mapped and 5% damping for 10% probability of exceedance in 50 years. Seismic hazard parameter b was estimated to be 0.62 to 0.98 using G R relation. Probabilistic seismic hazard analysis for Bangalore region was carried out considering six seismogenic sources. Shukla and Choudhury (2012) carried out site-specific probabilistic seismic hazard analysis for four typical port sites of Gujarat namely, Kandla, Mundra, Hazira and Dahej ports. PSHA was performed using seven ground motion prediction equations (GMPEs) including attenuation relation given by [24]. It is observed that the Mundra and Kandla port sites are most vulnerable sites for seismic hazard as estimated contingency level earthquake (CLE) ground motion is in order of 0.79 and 0.48 g for Mundra and Kandla port sites, respectively. And other hand comparatively less hazard with estimated CLE ground motion of 0.17 and 0.11 g for Hazira and Dahej port sites, respectively. Ramanna and Dodagaodar (2012) Probabilistic Seismic Hazard Analysis (PSHA) has been widely used to determine the peak ground acceleration in this paper the researcher have been estimated peak ground acceleration for 10% probability of exceedance in 50 years using three approaches namely, Cornell-McGuire approach, fixed kernel and adaptive kernel techniques. Raghukanth and Iyenger (2006) present seismic hazard and obtained a recurrence relation for Mumbai region using a catalog of 1274 past earthquakes from AD 1594 2002 specific to the region around Mumbai. The parameters of the Gutenberg and Richter equation were found to be a = 0.77 ± @IJMTER-2016, All rights Reserved 451
0.04 and b = 0.86 ± 0.02. These were used to compute the probability of ground motion that can be induced by each of the twenty three known faults that exist around the city. Ghobadi M H and Fereidooni D (2012) conduct a seismic hazard analysis for the city of Hamedan at the west of Iran. In this case the both of DSHA and PSHA approaches have been used for the assessment of seismic hazards and earthquake risk evaluation. Gutenberg Richter relationship method was used for available earthquake data for getting the a and b parameters were estimated 5.53 and 0.68, respectively. III. STEP BY STEP PROCEDURE FOR SEISMIC HAZARD ANALYSIS 3.1. Introduction A seismic hazard analysis is a physical phenomenon, which is ground shaking or ground failure. Seismic hazard is associated with an earthquake and that can produce adverse effects on human activities. To estimate the seismic hazards for a particular region, all possible sources of seismic activity have to be identified and their possible for generating future strong ground motion. Recognition of seismic sources requires some detective work like natures clues, some of which are clear and others quite obscure necessity to observed and interpreted. Seismic hazards may be analyzed deterministically, at particular earthquake scenario is assumed. Otherwise analyzed probabilistically, in which uncertainties in earthquake location, size and time of occurrence are explicitly considered. Although seismic hazard analysis is a critical part of the development of attenuation relations for design ground motions and site condition consists of four primary steps: 1) Identification and characterization of all sources 2) Selection of source-site distance parameter 3) Selection of controlling earthquake 4) Definition of hazard using controlling earthquake 3.1. Declustering process Declustering is the process of removal of dependent earthquake events (foreshocks and aftershocks) from an earthquake catalogue. The Poisoning assumption of earthquake recurrence implies that the earthquake event occurs randomly with no memory of time, size or location of any preceding event. Therefore, foreshocks and aftershocks are removed from the earthquake catalogue when a PSHA based on the Cornell-McGuire method is employed. 3.2. Earthquake catalogue completeness An important problem in defining a recurrence relation is ensuring completeness. For the time period under consideration, it is crucial to ensure that all earthquakes have been recorded for each magnitude range of interest. For example, there is a threshold magnitude below which earthquake occurrences have not been recorded completely in a given time frame for a given layout of seismographs and a given distribution of the population. The completeness interval for the higher magnitude classes would be relatively difficult to determine. The graph would exhibit a stepped behaviour due to the fact that stronger events tend to be separated by relatively long time intervals and sometimes occur within a short period, both owing to the physical nature of earthquakes in a seismogenic zone. IV. CASE STUDY Seismic hazard for Gujarat region has been studied, its located central western part of Gujarat state and which is surrounded by Arabian Sea. The area of Gujarat region has been subdivided into three major regions are: one is Kachchh region; second is Saurashtra region; and third one is Mainland Gujarat region. These cover the area of the 25 cities, which are selected for case study of present analysis. Further later on, four port sites: Kandla port, Mundra port, Hazira port and Dahej port are also considered for site specific ground estimation analysis [16]. @IJMTER-2016, All rights Reserved 452
Figure 4. Seismic Zones of Gujarat, IS: 1893-Part I (2002),[12],[16] Figure 4 shows the seismic zonation map of Gujarat state of India, These are the four colour codes like red, orange, yellow and blue. These are showing four different zones of the Gujarat state: zone 5, zone 4, zone 3 and zone 2. Among these, zone 5 is the most vulnerable one; zone 4 is lesser than that; zone 3 is lesser than that; and, zone 2 is least vulnerable in the Gujarat region. So, Gujarat is having seismic hazard for all the four zones as per our IS seismic design code 1893-part 1. Another reason is 2001 Bhuj earthquake is very well-known, very damaging earthquake, which occurred in Bhuj region of Gujarat state. And, it affected several other parts of the Gujarat as well [16]. V. ANALYSIS AND RESULTS 5.1. Deterministic seismic hazard analysis (DSHA) Here carried out the DSHA, it describes the potential for dangerous earthquake related natural phenomenon, is that the maximum realistic earthquake or maximum considered earthquake. The earthquake hazard for the site is peak ground acceleration is 0.57 g, which is resulting from an earthquake magnitude of 5.7 on a particular fault. It has been obtained Narmada Son fault at a distance of 11.42 kilometer from the site. This is just for one particular site for the DSHA result for any other site in the entire Gujarat region [16]. 5.2. Ground motion prediction equations (GMPEs) Ground Motion Prediction Equations are needed for the select to do the seismic hazard analysis. So that, these are six ground motion prediction equations, are mentioned over here [2]; [6]; [7]; [25]; [26]., 1997; and [23]. So, among these seven ground motion prediction equations or attenuation relationships, which they were used for this seismic hazard analysis, only the last one is for peninsular India; remaining all other from worldwide data, not from Indian data 5.3. Probabilistic seismic hazard analysis (PSHA) After study of all standard attenuation relationships, had calculated PSHA also, above study requires some parameters. 1) 40 number of major faults has been identified and those are considered only for @IJMTER-2016, All rights Reserved 453
further calculation, 2) number of magnitude recurrence relation which they were used, 3) 7 numbers of GMPE s were used. So, to create the PSHA map with the grid points of 8430 that will mean that you have to get the value of 8430 times of 4 and 7, Because already mentioned earlier that magnitude recurrence relation and the attenuation relation are gives the independency to these two. There are several computer software are available which can perform the probabilistic seismic hazard analysis like Seisrisk 3 crisis 2007 e z-frisk 88 etcetera. And using the M S excel of course, after getting all these grid point values have to assign them and prepare a two dimensional chart [16]. VI. CONCLUSION AND DISCUSSIONS Several; study s of seismic hazards assessments have been carried out whole country of India have been revived, in this conclusions. From the sismotectonic map city area has been developed. Comparison between instrumented PGA values from the Koyna-Warna region [15]. [22] independently carried out PSHA for the city of Mumbai. Seismicity of India has been addressed by many researchers in particular [11], [8], [24], [17], [19], [15] and [4]. As per IS 1893 (BIS 2002), seismic study area falls in the zones II and III in the seismic zonation map of India. Gujarat region has carryout the seismic hazard analysis for using regional seismicity parameters which are consistent with present state of seismicity in the Gujarat. The seismicity parameter (bvalue) is calculated by two different approaches namely Least Square Fit method and Maximum Likelihood method. These parameters constitute the basic framework required for seismic hazard assessment for Gujarat region. The b-values obtained using Least Square Fit are 0.41, 0.64 and 0.62 for Kachchh, Saurashtra and Mainland Gujarat, respectively. Though the peak ground acceleration (PGA) values are not in close agreement with those recommended by IS:1893 Part 1 (2002) for some cities but overall seismic ground motion distribution across the Gujarat region is more or less in agreement with regional distribution recommended by IS: 1893 Part 1 (2002). The present PSHA study results in close agreement with the seismic ground motions recommended by IS: 1893 Part 1 (2002). REFERENCES [1] Anbazhagan, P., Vinod, J. S., and Sitharam, T. G, Probabilistic seismic hazard Analysis for Bangalore, Nat. hazards, Vol.8, pp.145 166, 2009. [2] Abrahamson N A, Silva W J., Empirical response spectral attenuation relations for shallow crustal earthquakes, Seismol Res Lett, Vol.68, pp.94 127, (1997). [3] Akkar, S and J. J Bommer, Empirical equations for the prediction of PGA, PGV, and spectral accelerations in Europe, the Mediterranean region, and the Middle East, Seismological Research Letters, Vol.81, pp.195 206, (2010). [4] Bilham, R., Earthquakes in India and the Himalaya: Tectonics, Geodesy and History, Annals of Geophysics, Vol. 47, pp. 839-858, (2004). [5] Boore D M, Joyner W B, Fumal T E., Equations from estimating horizontal response spectra and peak acceleration from Western North American earthquakes: a summary of recent work, Seismol Res Lett 68:128 153, (1997). [6] Campbell K W., Empirical near-source attenuation relationships for horizontal and vertical components of peak ground acceleration, peak ground velocity, and pseudo-absolute acceleration response spectra. Seismol Res Lett 68:154 179, (1997). [7] Chandra, U., Earthquakes of Peninsular India, A Seismotectonic Study, Bulletin of the Seismological Society of America, Vol. 65, pp. 1387-1413, (1977). [8] Cornell, C.A., Engineering seismic risk analysis, Bulletin of Seismological Society of America, Vol. 58, pp.1583-1606, (1968). [9] Das, S., Gupta, I. D. and Gupta, V. K., A Probabilistic Seismic Hazard Analysis of Northeast India, Earthquake Spectra, Vol. 22, pp. 1-27, (2006). [10] Kaila, K.L., Gaur, V.K. and Narain Hari., Quantitative seismicity maps of India. Bull Seismol Soc Am., Vol.62 pp.1119 1131, (1972). @IJMTER-2016, All rights Reserved 454
[11] Kramer, S.L., Geotechnical Earthquake Engineering, Prentice-Hall, Upper Saddle River, New Jersey, (1996). [12] Gupta, I. D. Delineation of Probable Seismic Sources in India and Neighbourhood by a Comprehensive Analysis of Seismotectonic Characteristics of the Region, Soil Dynamics and Earthquake Engineering, Vol. 26, pp. 766 790, (2006). [13] Indian Standard 1893 (Part I), Criteria for Earthquake Resistance Design of Structures, Part-1: General Provisions and Buildings, Fifth revision, Bureau of Indian Standards, New Delhi, (2002). [14] Iyengar, R.N. and Raghukanth, S.T.G., Attenuation of strong ground motion in peninsular India, Seismological Research Letters, vol. 75(4), 530-540, (2004). [15] J. Shukla and D. Choudhury., Estimation of seismic ground motions using deterministic approach for major cities of Gujarat, Nat. Hazards Earth Syst. Sci., Vol. 12, pp.2019 2037, (2012). [16] Khattri, K. N., Seismic Hazard in Indian Region, Current Science, Vol. 62, pp. 109-116, (1992).. [17] Mandal, P., Intraplate Stress Distribution Induced by Topography and Crustal Density Heterogeneities Beneath the South Indian Shield, India, Tectonophysics, Vol. 302, pp. 159-172, (1999). [18] Margaris, B. N., C. Papazachos, C. Papaioannou, N. Theodulidis, I. Kalogeras, and A. Skarlatoudis., Ground motion attenuation relations for shallow earthquakes in Greece, in Proc. of the Twelfth European Conf. on Earthquake Engineering, 9 13 September, London, paper no. 385, (2002). [19] Mohammad, Ghobadi ; Davood, Fereidooni., Seismic hazard assessment of the city of Hamedan and its vicinity, west of Iran. Natural Hazards, Vol. 63 pp.1025-1038, (2012). [20] Parvez, A. I., Vaccari, F. and Panza, G. F., A Deterministic Seismic Hazard Map of India and Adjacent Areas, International Journal of Geophysics, Vol. 155, pp. 489-508, (2003). [21] Peterson, M. D., Rastogi, B. K., Schweig, G.E.S., Harmsen,S.C. and. Gomberg J.S., Sensitivity analysis of seismic hazard for the northwestern portion of the state of Gujarat, India, Technophysics, Vol. 390, pp.105-115, (2004). [22] Raghu Kanth S T G, Iyengar R N., Seismic hazard estimation for Mumbai city, Curr Sci 11:1486 1494, (2006). [23] Raghu Kanth S T G, Iyengar R N., Estimation of seismic spectral acceleration in peninsular India, J Earth Syst Sci Vol.3, pp.199-214, (2007). [24] Ramanna, C. K. and Dodagoudar, G. R., Seismic Hazard Analysis using the adaptive kernel density estimation technique for Chennai City, Pure and Applied Geophysics, (2012). [25] Sadigh K, Chang C Y, Egan J A, Makdisi F, Youngs R R., Attenuation relations for shallow crustal earthquakes based on California strong motion data, Seismol Res Lett, Vol. 68, pp.180-18, (1997). [26] Toro G R, Abrahamson N A, Schneider J F., Model of strong ground motions from earthquakes in central and eastern North America: best estimates and uncertainties, Seismol Res Lett, Vol.68, pp.41-58, (1997). [27] Tromans, I. J., and J. J. Bommer., The attenuation of strong-motion peaks in Europe, in Proc. of the Twelfth European Conf. On Earthquake Engineering, London, paper no. 394, (2002). @IJMTER-2016, All rights Reserved 455