Assessment of earthquake hazard in Panama based on seismotectonic regionalization

Assessment of seismic hazard in Panama is made using a seismotectonic regionalization model. The coefficients of Gumbel's Type-I distribution are calculated and return periods for several magnitudes are found. From these coefficients intensities, peak ground acceleration and earthquake hazard for a set of return periods and epicentral distances are estimated and substantial variations in the probability of occurrence are noted. The Panama Fracture Zone (PFZ) and the Panama-South America Suture Zone (PSZ) provinces are the most active in producing earthquakes with a magnitude of about 7.0 in less than 16 yr. Magnitude 7.0 earthquakes in the Azuero province have a return period of about 160 yr, whereas in the Panama Deformed Belt (PDB) province the return period for magnitude 7.5 events is about 175 yr.     

Probabilistic seismic hazard maps for Panama

Probabilistic seismic hazard maps in term of Modified Mercalli (MM) intensity are derived by applying the Cornell-McGuire method to four earthquake source zones in Panama and adjacent areas. The maps contain estimates of the maximum MM intensity for return periods of 5, 25 and 100 yr. The earthquake phenomenon is based on the point source model. The probabilistic iso-intensity map for a return period of 50 yr indicates that the Panama Suture Zone (PSZ) could experience a maximum (MM) intensity IX, and the Panama Fracture Zone (PFZ) an MM intensity VIII, for the rest of the area this varies from IV up to VIII. The present study intends to serve as a reference for more advanced approaches, to stimulate discussions and suggestions on the data base, assumptions and inputs, and path for the risk based assessment of the seismic hazard in the site selection and in the design of common buildings and engineering.     

Probabilities for the occurrences of medium to large earthquakes in northeast India and adjoining region

The return periods and occurrence probabilities related to medium and large earthquakes (M _w 4.0–7.0) in four seismic zones in northeast India and adjoining region (20°–32°N and 87°–100°E) have been estimated with the help of well-known extreme value theory using three methods given by Gumbel (1958), Knopoff and Kagan (1977) and Bury (1999). In the present analysis, the return periods, the most probable maximum magnitude in a specified time period and probabilities of occurrences of earthquakes of magnitude M ≥ 4.0 have been computed using a homogeneous and complete earthquake catalogue prepared for the period between 1897 and 2007. The analysis indicates that the most probable largest annual earthquakes are close to 4.6, 5.1, 5.2, 5.5 and 5.8 in the four seismic zones, namely, the Shillong Plateau Zone, the Eastern Syntaxis Zone, the Himalayan Thrusts Zone, the Arakan-Yoma subduction zone and the whole region, respectively. The most probable largest earthquakes that may occur within different time periods have been also estimated and reported. The study reveals that the estimated mean return periods for the earthquake of magnitude M _w 6.5 are about 6–7 years, 9–10 years, 59–78 years, 72–115 years and 88–127 years in the whole region, the Arakan-Yoma subduction zone, the Himalayan Thrusts Zone, the Shillong Plateau Zone and the Eastern Syntaxis Zone, respectively. The study indicates that Arakan-Yoma subduction zone has the lowest mean return periods and high occurrence probability for the same earthquake magnitude in comparison to the other zones. The differences in the hazard parameters from zone to zone reveal the high crustal heterogeneity and seismotectonics complexity in northeast India and adjoining regions.     

Quantitative assessment of landslide hazard along transportation lines using historical records

In this paper, a quantitative landslide hazard model is presented for transportation lines, with an example for a road and railroad alignment, in parts of Nilgiri hills in southern India. The data required for the hazard assessment were obtained from historical records available for a 21-year period from 1987 to 2007. A total of 901 landslides from cut slopes along the railroad and road alignment were included in the inventory. The landslides were grouped into three magnitude classes based on the landslide type, volume, scar depth, and run-out distance. To calculate landslide hazard, we estimated the total number of individual landslides per kilometer of the (rail) road for different return periods, based on the relationship between past landslides (recorded in our database) and triggering events. These were multiplied by the probability that the landslides belong to a given magnitude class. This gives the hazard for a given return period expressed as the number of landslides of a given magnitude class per kilometer of (rail) road. The relationship between the total number of landslides and the return period was established using a Gumbel distribution model, and the probability of landslide magnitude was obtained from frequency–volume statistics. The results of the analysis indicate that the total number of landslides, from 1- to 50-year return period, varies from 56 to 197 along the railroad and from 14 to 82 along the road. In total, 18 hazard scenarios were generated using the three magnitude classes and six return periods (1, 3, 5, 15, 25, and 50 years). The hazard scenarios derived from the model form the basis for future direct and indirect landslide risk analysis along the transportation lines. The model was validated with landslides that occurred in the year 2009.     

Some aspects of the seismic hazard in Panama City

The definition of earthquake sources in the Panama region on the basis of both tectonics and average seismicity rates, have recently led to the concept of a microplate surrounded by seismically active areas. The effects of these earthquakes on the place where the most important concentration of investments and population is located, the capital city of Panama, are analyzed in this paper using statistical approaches.The parameters of Gumbel's Type-I distribution of extreme values for a continuous interval of 60 yr annual maximum magnitudes were used to make probabilistic estimations of the seismic hazard in Panama City. An earthquake with magnitude 7.5 is capable of producing a modified Mercalli intensity VII in Panama City, provided the source distance is of the order of 100 km. This earthquake has a probability of occurrence of 69% in 50 yr.     

A quantitative appraisal of earthquake hazard parameters computed from Gumbel I method for different regions in and around Turkey

Useful information concerning the earthquake hazard parameters distributed in Turkey and the adjacent areas are estimated in the present work. Based on Gumbel’s I distribution parameters we are able to estimate the hazard values of the investigated area which are the mean return periods, the most probable maximum magnitude in the time period of t-years and the probability for an earthquake occurrence of magnitude ≥M during a time span of t-years. Figures concerning the spatial distribution of probabilities and the return periods are plotted and we considered them of particular interest for mapping the earthquake hazard in Turkey and the surrounding areas. These figures effectively produce a brief earthquake hazard atlas. The quantitative appraisal of the hazard parameters is useful for engineers, planners, etc., because it provides a tool for earthquake resistant design.     

Tsunamis and Tsunami Hazards in Central America

A tsunami catalogue for Central America is compiledcontaining 49 tsunamis for the period 1539–1996,thirty seven of them are in the Pacific and twelve inthe Caribbean. The number of known tsunamis increaseddramatically after the middle of the nineteenth century,since 43 events occurred between 1850 and 1996. This isprobably a consequence of the lack of populationliving near the coast in earlier times.The preliminary regionalization of the earthquakessources related to reported tsunamis shows that, inthe Pacific, most events were generated by theCocos-Caribbean Subduction Zone (CO-CA). At theCaribbean side, 5 events are related with the NorthAmerican-Caribbean Plate Boundary (NA-CA) and 7 withthe North Panama Deformed Belt (NPDB).There are ten local tsunamis with a specific damagereport, seven in the Pacific and the rest in theCaribbean. The total number of casualties due to localtsunamis is less than 455 but this number could behigher. The damages reported range from coastal andship damage to destruction of small towns, and theredoes not exist a quantification of them.A preliminary empirical estimation of tsunami hazardindicates that 43% of the large earthquakes (Ms 7.0) along the Pacific Coast of Central America and100% along the Caribbean, generate tsunamis. On thePacific, the Guatemala–Nicaragua coastal segment hasa 32% probability of generating tsunamis after largeearthquakes while the probability is 67% for theCosta Rica–Panama segment. Sixty population centers onthe Pacific Coast and 44 on the Caribbean are exposedto the impact of tsunamis. This estimation alsosuggests that areas with higher tsunami potential inthe Pacific are the coasts from Nicaragua to Guatemalaand Central Costa Rica; on the Caribbean side, Golfode Honduras Zone and the coasts of Panama and CostaRica have major hazard. Earthquakes of magnitudelarger than 7 with epicenters offshore or onshore(close to the coastline) could trigger tsunamis thatwould impact those zones.     

Effect of Aftershocks on Earthquake Hazard Estimation: An Example from the North Anatolian Fault Zone

In order to investigate the effect of aftershocks on earthquake hazard estimation, earthquake hazard parameters (_m, b and M_max) have been estimated by the maximum likelihood method from the main shocks catalogue and the raw earthquakes catalogue for the North Anatolian Fault Zone (NAFZ). The main shocks catalogue has been compiled from the raw earthquake catalogue by eliminating the aftershocks using the window method. The raw earthquake catalogue consisted of instrumentally detected earthquakes between 1900 and 1992, and historical earthquakes that occurred between 1000–1900. For the events of the mainshock catalogue the Poisson process is valid and for the raw earthquake catalogue it does not fit. The paper demonstrates differences in the hazard outputs if on one hand the main catalogues and on the other hand the raw catalogue is used. The maximum likelihood method which allows the use of the mixed earthquake catalogue containing incomplete (historical) and complete (instrumental) earthquake data is used to determine the earthquake hazard parameters. The maximum regional magnitude (M_max, the seismic activity rate (_m), the mean return period (R) and the b value of the magnitude-frequency relation have been estimated for the 24°–31° E, 31°–41° E, 41°–45° E sections of the North Anatolian Fault Zone from the raw earthquake catalogue and the main shocks catalogue. Our results indicate that inclusion of aftershocks changes the b value and the seismic activity rate _m depending on the proportion of aftershocks in a region while it does not significantly effect the value of the maximum regional magnitude since it is related to the maximum observed magnitude. These changes in the earthquake hazard parameters caused the return periods to be over- and underestimated for smaller and larger events, respectively.     

Seismic hazard in mega city Kolkata, India

The damages caused by recent earthquakes in India have been a wake up call for people to take proper mitigation measures, especially the major cities that lie in the high seismic hazard zones. Kolkata City, with thick sediment deposit (∼12 km), one of the earliest cities of India, is an area of great concern as it lies over the Bengal Basin and lies at the boundary of the seismic zones III and IV of the zonation map of India. Kolkata has been affected by the 1897 Shillong earthquake, the 1906 Calcutta earthquake, and the 1964 Calcutta earthquake. An analysis on the maximum magnitude and b-value for Kolkata City region is carried out after the preparation of earthquake catalog from various sources. Based on the tectonic set-up and seismicity of the region, five seismic zones are delineated, which can pose a threat to Kolkata in the event of an earthquake. They are broadly classified as Zone 1: Arakan-Yoma Zone (AYZ), Zone 2: Himalayan Zone (HZ), Zone 3: Shillong Plateau Zone (SPZ), Zone 4: Bay of Bengal Zone (BBZ), and Zone 5: Shield Zone (SZ). The maximum magnitude (m _max) for Zones 1, 2, 3, 4, and 5 are 8.30 ± 0.51, 9.09 ± 0.58, 9.20 ± 0.51, 6.62 ± 0.43 and 6.61 ± 0.43, respectively. A probability of 10% exceedance value in 50 years is used for each zone. The probabilities of occurrences of earthquakes of different magnitudes for return periods of 50 and 100 years are computed for the five seismic zones. The Peak Ground Acceleration (PGA) obtained for Kolkata City varies from 0.34 to 0.10 g.     

Seismic hazard across Bulgaria and neighbouring areas: extreme magnitude recurrence and strong ground shaking

A contemporary probabilistic seismic hazard assessment (PSHA) study for Bulgaria and the surrounding Balkan area is performed under constraints of a newly developed, fit-for-purpose historical earthquake catalogue and the theory of extreme values. Sensitivity analyses are first adopted as preparatory reviews on subsets of the adopted data to determine suitable values for the constraints of cut-off magnitude threshold, sample extreme interval and start year of catalogue data to impose on the parent database for both the full region considered as well as significant urban centres within it. Maximum estimates are then determined for magnitude recurrence hazard using Gumbel’s third asymptotic extreme values distribution for return periods of 50 and 100 years, and also these time intervals at 90 % probability of not being exceeded (PNBE). Gumbel’s first asymptotic extreme values distribution is also used with carefully selected, geographically relevant ground motion models for peak horizontal ground acceleration, PGA(h), and peak horizontal ground velocity, PGV(h), for the same return periods. The former provides direct comparison with the current EUROCODE 8 anti-seismic building code standard promoted across Europe, the previous GSHAP and SESAME hazard mapping projects as well as a number of recent studies. Sofia is forecast an upper bound magnitude of 7.33 M _w (±0.78) compared with 7.31 M _w (±0.55) for the full Balkan extent and 7.24 M _w (±0.70) for the political triple junction area of southwest Bulgaria, viz., Bulgaria, Greece and The Former Yugoslav Republic of Macedonia. Sofia is also forecast a 475-year return period (equivalent to a 50-year return period at 90 % PNBE) magnitude of 7.27 M _w, with an equivalent PGA (the standard EUROCODE 8 metric) of 156 cm s^?2 and PGV of 13 cm s^?1.     

Flood frequency analysis of river swat using Log Pearson type 3, Generalized Extreme Value,Normal, and Gumbel Max distribution methods

In Pakistan, floods are among the most devastating and recurring natural hazards. Flood hazard assessment requires flood event magnitude and probability of occurrence. Flood frequency analysis is the most common technique used for the at-site estimation of flood recurrence magnitude. This paper evaluates four most commonly used distribution methods, i.e., Generalized Extreme Value (GEV), Log Pearson 3 (LP3), Gumbel Max, and Normal for the flood frequency and estimation of flood recurrence. Different hydrological stations data namely Khwazakhela, Chakdarra, Panjkora, and Munda Headwork located at Swat river was taken from Provincial Irrigation Department, Khyber Pakhtunkhwa. The analysis is done for 5-, 10-, 25-, 50-, and 100-year return periods by using annual maximum discharge data from 1980 to 2016 (37 yr). Three goodness-of-fit tests were applied to the fitted distributions, i.e., Kolmogorov–Smirnov, Anderson–Darling, and Chi-squared at 5% significance level. Results indicate that LP3 and GEV were ranked top two distributions at all locations while Gumbel Max and Normal were the least fitted having rank 3 and 4, respectively. Based on the goodness-of-fit ranking, LP3 was selected for the estimation of flood magnitude and return periods at Khwazakhela. Designed hydrographs based on probabilistic approach and flood 2010 hydrograph are presented for flood simulation.     

Crustal deformation in the northern Andes – A new GPS velocity field

We present a velocity field for northwestern South America and the southwest Caribbean based on GPS Continuously Operating Reference Stations in Colombia, Panama, Ecuador and Venezuela. This paper presents the first comprehensive model of North Andean block (NAB) motion. We estimate that the NAB is moving to the northeast (060°) at a rate of 8.6 mm/yr relative to the South America plate. The NAB vector can be resolved into a margin-parallel (035°) component of 8.1 mm/yr rigid block motion and a margin-normal (125°) component of 4.3 mm/yr. This present-day margin-normal shortening rate across the Eastern Cordillera (EC) of Colombia is surprising in view of paleobotanical, fission-track, and seismic reflection data that suggest rapid uplift (7 km) and shortening (120 km) in the last 10 Ma. We propose a “broken indenter” model for the Panama-Choco arc, in which the Choco arc has been recently accreted to the NAB, resulting in a rapid decrease in shortening in the EC. The Panama arc is colliding eastward with the NAB at approximately 15–18 mm/yr, and the Panama-Choco collision may have been responsible for much of the uplift of the EC. The present on-going collision poses a major earthquake hazard in northwestern Colombia from the Panama border to Medellin area. Since the northeastward margin-parallel motion of the NAB is now greater than the rate of shortening in the EC, northeast trending right-lateral strike-slip faulting is the primary seismic hazard for the 8 million inhabitants of Bogota, the capital city of Colombia. There continues to be a high risk of a great megathrust earthquakes in southern Colombia along the Ecuador-Colombia trench. Trench earthquakes have only released a fraction of the energy accumulated in the Ecuador-Colombia trench since the 1906 Ecuador earthquake, and interseismic strain is accumulating rapidly at least as far north as Tumaco, the rupture area of the 1979 earthquake.     

Hydrological impacts of land use–land cover change and detention basins on urban flood hazard: a case study of Poisar River basin,Mumbai, India

Flooding in urban area is a major natural hazard causing loss of life and damage to property and infrastructure. The major causes of urban floods include increase in precipitation due to climate change effect, drastic change in land use–land cover (LULC) and related hydrological impacts. In this study, the change in LULC between the years 1966 and 2009 is estimated from the toposheets and satellite images for the catchment of Poisar River in Mumbai, India. The delineated catchment area of the Poisar River is 20.19 km². For the study area, there is an increase in built-up area from 16.64 to 44.08% and reduction in open space from 43.09 to 7.38% with reference to total catchment area between the years 1966 and 2009. For the flood assessment, an integrated approach of Hydrological Engineering Centre-Hydrological Modeling System (HEC-HMS), HEC-GeoHMS and HEC-River analysis system (HEC-RAS) with HEC-GeoRAS has been used. These models are integrated with geographic information system (GIS) and remote sensing data to develop a regional model for the estimation of flood plain extent and flood hazard analysis. The impact of LULC change and effects of detention ponds on surface runoff as well as flood plain extent for different return periods have been analyzed, and flood plain maps are developed. From the analysis, it is observed that there is an increase in peak discharge from 2.6 to 20.9% for LULC change between the years 1966 and 2009 for the return periods of 200, 100, 50, 25, 10 and 2 years. For the LULC of year 2009, there is a decrease in peak discharge from 10.7% for 2-year return period to 34.5% for 200-year return period due to provision of detention ponds. There is also an increase in flood plain extent from 14.22 to 42.5% for return periods of 10, 25, 50 and 100 years for LULC change between the year 1966 and year 2009. There is decrease in flood extent from 4.5% for 25-year return period to 7.7% for 100-year return period and decrease in total flood hazard area by 14.9% due to provisions of detention pond for LULC of year 2009. The results indicate that for low return period rainfall events, the hydrological impacts are higher due to geographic characteristics of the region. The provision of detention ponds reduces the peak discharge as well as the extent of the flooded area, flood depth and flood hazard considerably. The flood plain maps and flood hazard maps generated in this study can be used by the Municipal Corporation for flood disaster and mitigation planning. The integration of available software models with GIS and remote sensing proves to be very effective for flood disaster and mitigation management planning and measures.     

Probabilities of earthquake occurrences in Mainland Southeast Asia

The frequency–magnitude distributions of earthquakes are used in this study to estimate the earthquake hazard parameters for individual earthquake source zones within the Mainland Southeast Asia. For this purpose, 13 earthquake source zones are newly defined based on the most recent geological, tectonic, and seismicity data. A homogeneous and complete seismicity database covering the period from 1964 to 2010 is prepared for this region and then used for the estimation of the constants, a and b, of the frequency–magnitude distributions. These constants are then applied to evaluate the most probable largest magnitude, the mean return period, and the probability of earthquake of different magnitudes in different time spans. The results clearly show that zones A, B, and E have the high probability for the earthquake occurrence comparing with the other seismic zones. All seismic source zones have 100 % probability that the earthquake with magnitude ≤6.0 generates in the next 25 years. For the Sagaing Fault Zone (zones C), the next Mw 7.2–7.5 earthquake may generate in this zone within the next two decades and should be aware of the prospective Mw 8.0 earthquake. Meanwhile, in Sumatra-Andaman Interplate (zone A), an earthquake with a magnitude of Mw 9.0 can possibly occur in every 50 years. Since an earthquake of magnitude Mw 9.0 was recorded in this region in 2004, there is a possibility of another Mw 9.0 earthquake within the next 50 years.     

Seismic hazard assessment for low-seismicity areas — Case study: Northern Germany

The problem of assessing seismic hazard in low-seismicity areas becomes obvious in many practical applications. A typical low-seismicity area, which experienced damaging earthquakes in historical times, is the North German Plain, for which a case study is presented. It is shown how seismic hazard assessments are influenced by different interpretations of historical key earthquakes, changes in b-value as well as variations of the upper bound magnitude assumed for the seismic source regions. The latter strongly influences the hazard results in the case of very low b-values for long return periods.     

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.     

Probabilistic seismic hazard maps for the sultanate of Oman

This study presents the results of the first probabilistic seismic hazard assessment (PSHA) in the framework of logic tree for Oman. The earthquake catalogue was homogenized, declustered, and used to define seismotectonic source model that characterizes the seismicity of Oman. Two seismic source models were used in the current study; the first consists of 26 seismic source zones, while the second is expressing the alternative view that seismicity is uniform along the entire Makran and Zagros zones. The recurrence parameters for all the seismogenic zones were determined using the doubly bounded exponential distribution except the zones of Makran, which were modelled using the characteristic distribution. Maximum earthquakes were determined and the horizontal ground accelerations in terms of geometric mean were calculated using ground-motion prediction relationships developed based upon seismic data obtained from active tectonic environments similar to those surrounding Oman. The alternative seismotectonic source models, maximum magnitude, and ground-motion prediction relationships were weighted and used to account for the epistemic uncertainty. Hazard maps at rock sites were produced for 5?% damped spectral acceleration (SA) values at 0.1, 0.2, 0.3, 1.0 and 2.0?s spectral periods as well as peak ground acceleration (PGA) for return periods of 475 and 2,475?years. The highest hazard is found in Khasab City with maximum SA at 0.2?s spectral period reaching 243 and 397?cm/s² for return periods 475 and 2,475 years, respectively. The sensitivity analysis reveals that the choice of seismic source model and the ground-motion prediction equation influences the results most.     

2017年8月8日四川九寨沟地震诱发地质灾害空间分布规律及易发性评价初步研究

2017年8月8日四川省九寨沟发生7.0级地震,是继2008年汶川大地震和2013年芦山地震之后,四川省发生的又一次7.0级以上的强震。为了给现场救援和震后地质灾害防治提供科学依据,本文作者第一时间赶赴现场,并基于震前和震后高精度遥感影像,完成了震区地质灾害解译和复核工作,共解译地质灾害1883处,主要以中小型浅层滑坡和崩塌为主。基于解译结果,对同震灾害的空间分布规律和控制因素（距断层距离、地面峰值加速度PGA、高程、坡度和坡向等）进行了分析,研究表明地质灾害主要沿北西-南东向呈带状分布,且沿公路、沟谷较为发育,在野外推测发震断层2km范围内高度集中,呈现明显的断层效应,但与塔藏断层之间的空间相关性相对较弱。在上述分析的基础上,采用逻辑回归模型,利用地震、地形和地质3大因素（8个因子）对地质灾害易发性进行了快速评价和分区,经统计校验证明该模型的准确率达0.851,模型精度较高。     

Seismic Zoning of Croatia

The territory of Croatia and neighboring regions is divided into 17 seismic source zones, considering available seismological and geological data. On this basis, seismic hazard elements (seismicity rate, maximum magnitude, b-value, probabilities of exceedance and return periods for a predefined set of magnitudes) are computed using the maximum likelihood method appropriate for treating data-sets with variable completeness thresholds. The values of long term expected peak horizontal acceleration obtained by using a combination of the deterministic and the probabilistic procedure are the highest in the Dubrovnik zone, while the Zagreb zone has the highest earthquake hazard in the continental part of the country. This revised version was published online in June 2006 with corrections to the Cover Date.     

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.