Deterministic seismic hazard assessment for Sultanate of Oman

The Sultanate of Oman forms the southeastern part of the Arabian plate, which is surrounded by relatively high active tectonic zones. Studies of seismic risk assessment in Oman have been an important on-going socioeconomic concern. Using the results of the seismic hazard assessment to improve building design and construction is an effective way to reduce the seismic risk. In the current study, seismic hazard assessment for the Sultanate of Oman is performed through the deterministic approach with particular attention on the uncertainty analysis applying a recently developed method. The input data set contains a defined seismotectonic model consisting of 26 seismic zones, maximum magnitudes, and 6 alternative ground motion prediction equations that were used in four different tectonic environments: obduction zone earthquake (Zagros fold thrust belt), subduction zone earthquakes (Makran subduction zones), normal and strike-slip transform earthquakes (Owen and Gulf of Aden zones), and stable craton seismicity (Arabian stable craton). This input data set yielded a total of 76 scenarios at each point of interest. A 10 % probability that any of the 76 scenarios may exceed the largest median ground acceleration is selected. The deterministic seismic hazards in terms of PGA, 5 % damped spectral acceleration at 0.1, 0.2, 1.0 and 2.0 s are performed at 254 selected points. The ground motion was calculated at the 50th and 84th percentile levels for selected probability of exceeding the median value. The largest ground motion in the Sultanate of Oman is observed in the northeastern part of the country.     

Seismic risk assessment at the proposed site of Gemsa wind power station,southwestern coast of Gulf of Suez,Egypt

Gemsa has been chosen as the site for one of a new generation of power stations along the south-western margin of the Gulf of Suez. This site has been affected by a number of destructive earthquakes (M_w> 5), in addition to large number of earthquakes with magnitudes of less than 5. In this study seismic activities in the region were collected and re-evaluated, and the main earthquake prone zones were identified. It is indicated that this site is affected by the southern Gulf of Suez, northern Red Sea and Gulf of Aqaba source zones. The southern Gulf of Suez source zone is the nearest to the proposed site. The stochastic simulation method has been applied to estimate the Peak GroundAcceleration at the site of the proposed Gemsa power plant. It was noticed that the pseudo-spectral acceleration (PSA) reaches 175 cm/sec2 resulting from the southern Gulf of Suez seismic source. In addition, the response spectrum was conducted with a damping value of 5% of the critical damping, and the predominant period reached 0.1sec at the site. These results should be taken into consideration by civil engineers and decision-makers for designing earthquake resistant structures.     

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.     

Present Kinematic Regime and Recent Seismicity of Gulf Suez,Egypt

In this study we computed recent seismicity and present kinematic regime in the northern and middle zones of Gulf of Suez as inferred from moment tensor settlings and focal mechanism of local earthquakes that happened in this region. On 18 and 22 of July, 2014 two moderate size earthquakes of local magnitudes 4.2 and 4.1 struck the northern zone of Gulf of Suez near Suez City. These events are instrumentally recorded by Egyptian National Seismic Network (ENSN). The earthquakes have been felt at Suez City and greater Cairo metropolitan zone while no losses were reported. The source mechanism and source parameters of the calculated events were considered by the near-source waveform data listed at very broadband stations of ENSN and supported by the P-wave polarity data of short period stations. The new settling method and software used deem the action of the source time function, which has been ignored in most of the program series of the moment tensor settling analysis with near source seismograms. The obtained results from settling technique indicate that the estimated seismic moments of both earthquakes are 0.6621E + 15 and 0.4447E + 15 Nm conforming to a moment magnitude Mw 3.8 and 3.7 respectively. The fault plan settlings obtained from both settling technique and polarity of first-arrival indicate the dominance of normal faulting. We also evaluated the stress field in north and middle zones of Gulf of Suez using a multiple inverse method. The prime strain axis shows that the deformation is taken up mainly as stretching in the E–W and NE–SW direction.     

Expected peak ground acceleration in Uttarakhand Himalaya,India region from a deterministic hazard model

A method of seismic zonation based on the deterministic modeling of rupture planes is presented. Finite rupture planes along identified lineaments are modeled in the Uttarakhand Himalaya based on the semi empirical technique of Midorikawa (Tectonophysics 218:287–295, 1993). The expected peak ground acceleration thus estimated from this technique is divided into different zones similar to zones proposed by the Bureau of Indian standard, BIS (Indian standards code of practice for earthquake-resistant design of structures, 2002). The proposed technique has been applied to Kumaon Himalaya area and the surrounding region for earthquakes of magnitude M > 6.0. Approximately 56000 km² study area is classified into the highest hazard zone V with peak accelerations of more than 400 cm/s². This zone V includes the cities of the Dharchula, Almora, Nainital, Haridwar, Okhimath, Uttarkashi, Pithorahargh, Lohaghat, Munsiari, Rudraprayag, and Karnprayag. The Sobla and Gopeshwar regions belong to zone IV, where peak ground accelerations of the order from 250 to 400 cm/s² can be expected. The prepared map shows that epicenters of many past earthquakes in this region lie in zone V, and hence indicating the utility of developed map in defining various seismic zones.     

地震动力学与板块构造:地震灾害评估

作为地震灾害评估的理论基础,地震动力学主要研究与地震活动有关的断裂机制、破裂过程、震源辐射和由此而引起的地震波的传播及地面运动规律。对地震力学、震源辐射和能量释放等经典理论问题进行了系统研究。在此基础上,应用最新的定量地震学研究方法,以逻辑树的形式综合地震、地质和大地测量资料,提供了不同构造环境和断裂机制条件下地震灾害评估的概率分析和确定性分析实例。用于震源分析的典型构造类型包括板内地壳震源层、地壳活动断层及其速率、板块俯冲界面和俯冲板片。由于输入模型中不确定因素的存在,如输入参数的随机性和科学分析方法本身的不确定性,对分析结果的不确定性需审慎对待。通常对不同的模型或参量,包括地面衰减模型,进行加权平均可较为合理地减小结果的偏差：概率分析和确定性分析方法的结合亦为可取之有效途径。     

Seismicity of Sinai Peninsula, Egypt

The Sinai Peninsula has a triangular shape between the African and Arabian Plates and is bounded from the western and eastern borders by the Gulf of Suez and Gulf of Aqaba–Dead Sea rift systems, respectively. It is affected by strong and destructive earthquakes (e.g., March 31, 1969 and November 22, 1995) and moderate earthquakes (m _b?>?5) throughout its history. After the installation of the Egyptian National Seismic Network (ENSN), a great number of earthquakes has been recorded within and around Sinai. Consequently, the seismogenic source zones and seismotectonic behavior can be clearly identified. Available data, including both historical and instrumental (1900–1997), have been collected from national and international data centers. While the data from 1998 till December 2007 are gathered from ENSN bulletins. The seismogenic source zones that might affect Sinai Peninsula are defined more precisely in this work depending on the distribution of earthquakes, seismicity rate (a value), b value, and fault plane solution of the major earthquakes. In addition, the type of faults prevailed and characterized these zones. It is concluded that the Gulf of Aqaba zone–Dead Sea transform zone, Gulf of Suez rift zone, Cairo–Suez District zone, and Eastern Mediterranean dislocation zone represent the major effective zones for Sinai. Furthermore, there are two local seismic zones passing through Sinai contributing to the earthquake activities of Sinai, these are the Negev shear zone and Central Sinai fault (Themed fault) zone. The source parameters, a and b values, and the maximum expected moment magnitude have been determined for each of these zones. These results will contribute to a great extent in the seismic hazard assessment and risk mitigation studies for Sinai Peninsula to protect the developmental projects.     

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

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

Sensitivity analysis of seismic hazard for the northwestern portion of the state of Gujarat, India

We test the sensitivity of seismic hazard to three fault source models for the northwestern portion of Gujarat, India. The models incorporate different characteristic earthquake magnitudes on three faults with individual recurrence intervals of either 800 or 1600 years. These recurrence intervals imply that large earthquakes occur on one of these faults every 266–533 years, similar to the rate of historic large earthquakes in this region during the past two centuries and for earthquakes in intraplate environments like the New Madrid region in the central United States. If one assumes a recurrence interval of 800 years for large earthquakes on each of three local faults, the peak ground accelerations (PGA; horizontal) and 1-Hz spectral acceleration ground motions (5% damping) are greater than 1 g over a broad region for a 2% probability of exceedance in 50 years' hazard level. These probabilistic PGAs at this hazard level are similar to median deterministic ground motions. The PGAs for 10% in 50 years' hazard level are considerably lower, generally ranging between 0.2 g and 0.7 g across northwestern Gujarat. Ground motions calculated from our models that consider fault interevent times of 800 years are considerably higher than other published models even though they imply similar recurrence intervals. These higher ground motions are mainly caused by the application of intraplate attenuation relations, which account for less severe attenuation of seismic waves when compared to the crustal interplate relations used in these previous studies. For sites in Bhuj and Ahmedabad, magnitude (M) 7 3/4 earthquakes contribute most to the PGA and the 0.2- and 1-s spectral acceleration ground motion maps at the two considered hazard levels.     

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.     

Seismic hazard analysis of Lucknow considering local and active seismic gaps

The Himalayas are one of very active seismic regions in the world where devastating earthquakes of 1803 Bihar–Nepal, 1897 Shillong, 1905 Kangra, 1934 Bihar–Nepal, 1950 Assam and 2011 Sikkim were reported. Several researchers highlighted central seismic gap based on the stress accumulation in central part of Himalaya and the non-occurrence of earthquake between 1905 Kangra and 1934 Bihar–Nepal. The region has potential of producing great seismic event in the near future. As a result of this seismic gap, all regions which fall adjacent to the active Himalayan region are under high possible seismic hazard due to future earthquakes in the Himalayan region. In this study, the study area of the Lucknow urban centre which lies within 350 km from the central seismic gap has been considered for detailed assessment of seismic hazard. The city of Lucknow also lies close to Lucknow–Faizabad fault having a seismic gap of 350 years. Considering the possible seismic gap in the Himalayan region and also the seismic gap in Lucknow–Faizabad fault, the seismic hazard of Lucknow has been studied based on deterministic and the probabilistic seismic hazard analysis. Results obtained show that the northern and western parts of Lucknow are found to have a peak ground acceleration of 0.11–0.13 g, which is 1.6- to 2.0-fold higher than the seismic hazard compared to the other parts of Lucknow.     

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 hazard of Egypt

Earthquake hazard parameters such as maximum expected magnitude,M _max, annual activity rate,, andb value of the Gutenberg-Richter relation have been evaluated for two regions of Egypt. The applied maximum likelihood method permits the combination of both historical and instrumental data. The catalogue used covers earthquakes with magnitude 3 from the time interval 320–1987. The uncertainties in magnitude estimates and threshold of completeness were taken into account. The hazard parameter determination is performed for two study areas. The first area, Gulf of Suez, has higher seismicity level than the second, all other active zones in Egypt.b-values of 1.2 ± 0.1 and 1.0 ± 0.1 are obtained for the two areas, respectively. The number of annually expected earthquakes with magnitude 3 is much larger in the Gulf of Suez, 39 ± 2 than in the other areas, 6.1 ± 0.5. The maximum expected magnitude is calculated to be 6.5 ± 0.4 for a time span of 209 years for the Gulf of Suez and 6.1 ± 0.3 for a time span of 1667 years for the remaining active areas in Egypt. Respective periods of 10 and 20 years were reported for earthquakes of magnitude 5.0 for the two subareas.     

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.     

Fuzzy probabilistic seismic hazard analysis with applications to Kunming city,China

China is prone to highly frequent earthquakes due to specific geographical location, which could cause significant losses to society and economy. The task of seismic hazard analysis is to estimate the potential level of ground motion parameters that would be produced by future earthquakes. In this paper, a novel method based on fuzzy logic techniques and probabilistic approach is proposed for seismic hazard analysis (FPSHA). In FPSHA, we employ fuzzy sets for quantification of earthquake magnitude and source-to-site distance, and fuzzy inference rules for ground motion attenuation relationships. The membership functions for earthquake magnitude and source-to-site distance are provided based on expert judgments, and the construction of fuzzy rules for peak ground acceleration relationships is also based on expert judgment. This methodology enables to include aleatory and epistemic uncertainty in the process of seismic hazard analysis. The advantage of the proposed method is in its efficiency, reliability, practicability, and precision. A case study is investigated for seismic hazard analysis of Kunming city in Yunnan Province, People’s Republic of China. The results of the proposed fuzzy logic-based model are compared to other models, which confirms the accuracy in predicting the probability of exceeding a certain level of the peak ground acceleration. Further, the results can provide a sound basis for decision making of disaster reduction and prevention in Yunnan province.     

A deterministic approach for preparation of seismic hazard maps in North East India

A method of seismic zonation based on deterministic modeling of rupture plane is presented in this work. This method is based on the modeling of finite rupture plane along identified lineaments in the region using the semi-empirical technique, of Midorikawa [(1993) Tectonophysics 218:287–295]. The modeling procedure follows ω² scaling law, directivity effects, and other strong motion parameters. The technique of zonation is applied for technoeconomically important NE part of Brahmaputra valley that falls in the seismic gap region of Himalaya. Zonation map prepared for Brahmaputra valley for earthquakes of magnitude M > 6.0 show that approximately 90,000 km² area fall in the highly hazardous zone IV, which covers region that can have peak ground accelerations of order more than 250 cm/s². The zone IV covers the Tezu, Tinsukia, Dibrugarh, Ziro, North Lakhimpur, Itanagar, Sibsagar, Jorhat, Golaghat, Wokha, Senapati, Imphal, and Kohima regions. The Pasighat, Daring, Basar, and Seppa region belong to zone III with peak ground accelerations of the order 200–250 cm/s². The seismic zonation map obtained from deterministic modeling of the rupture is consistent with the historical seismicity map and it has been found that the epicenter of many moderate and major earthquakes fall in the identified zones.     

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.     

Seismological hazard assessment for a site in Northern Germany, an area of low seismicity

Potential sites of nuclear waste deposits in the Federal Republic of Germany are situated in areas of low seismicity. Nevertheless, seismic hazard assessment has to be performed for a very long time period in order to prove the facilities of the repositories able to withstand seismic induced loads; even though there is a considerable debate whether or not it is possible to quantify the seismic risk in such an area.

A combination of deterministic and probabilistic methods is used to assess the seismic hazard for a site in Northern Germany, fulfilling the standards of the German building code for nuclear power plants. As an example, the site of the former iron ore mine Konrad is investigated. The deterministic method is based on the assumption that the strongest earthquakes inside a tectonic region can happen everywhere there, also near the selected site. For the probabilistic method, several models describing the seismicity in an area of 200 km around the site are used to show the influence of the variability in input parameters, like the maximum intensity of each source region on the exceeding probability of the site intensity. It can be shown that the seismic hazard of a site in an area of low seismicity is mainly caused by the effects of distant but strong source regions and the background seismicity for very low probabilities.

Probabilistic evaluation has the advantage of quantifying the seismic risk. But deterministic and probabilistic methods together seem a practical tool for mutual control of the results and to overcome the weakness of each approach alone. The historical German earthquake catalog with an observation period of about 1200 years is the basis for the input data for a probabilistic model. From a deep knowledge of geological development and structural geology, the time history of the surrounding faults is developed. Indications were found that the nearest and most important fault was active at least 5 Ma ago. The combination of both seismicity and tectonics provides the basis for a long term prognostic with probabilities of exceedance in the order of 10⁻⁵ per year.

For the investigated site the following parameters were derived: site intensity as a function of exceeding probability; site acceleration; strong motion duration; site dependent response spectra for the surface and the underground inside the mine.     

Seismic Hazard Assessment in the Area of Mystras-Sparta,South Peloponnesus,Greece, Based on Local Seismotectonic,Seismic, Geologic Information and on Different Models of Rupture Propagation

Strong seismic events once again confirm the view that great destructive earthquakes are produced by the reactivation of pre-existing faults although they have usually remained inactive for many, perhaps thousands of years. It is evident that such active seismogenic zones, with little or no seismicity, have presumably been ignored in the determination of the region's seismic hazard. At south Peloponnesus, Greece, is situated at Taygetos mountain. At its eastern front lies a large normal fault system, the southern segment being the Sparta fault. This area has been characterized by low seismicity for the last 25 centuries. However, during the 6th and 5th centuries B.C. several destructive earthquakes have been reported. That of 464 B.C., was the most destructive and devastated the city of Sparta. Detailed morphotectonic observations of this area, suggest that the earthquake of 464 B.C. could be related to the most recent reactivation of this fault. The ground accelerations that would be produced by a future activation of the Sparta fault, were calculated, by applying a method which takes into account information mainly from the seismotectonic parameters of the Sparta fault, the rupture pattern, the properties of the propagation medium and the local ground conditions. Moreover, these results were compared with those of other independent studies based mainly on the seismic data of the area. This method estimated greater expected values of ground acceleration than those computed by the conventional seismic hazard methods. The highest values correspond to the activation of the Sparta fault either in a unilateral rupture, which would start from the southernmost point of the fault, or in a circular one. Furthermore, an increase is observed of the order of 50% in the ground acceleration values in unconsolidated soft ground in relation to the corresponding values of hard ground. This revised version was published online in June 2006 with corrections to the Cover Date.     

Seismic hazard assessment of northern Pakistan

This article presents probabilistic seismic hazard analyses of northern Pakistan region carried out to produce macro-seismic hazard maps for the region that define new regional ground motion design parameters for 95-, 475-, 975- and 2475-year return period earthquakes as regional contour maps and horizontal uniform hazard at important cities. The Cornell–McGuire approach (Cornell in Bull Seismol Soc Am 58(05):1583–1606, 1968; McGuire in FORTRAN computer program for seismic risk analysis. US Geological Survey, Open file Report, 76-6768, 1976) is used to carry out the analyses at 0.1° rectangular grid. The seismotectonic model of the region used in analysis consists of shallow and deep area zones differentiated based on the focal depths of the earthquakes. Earthquake catalogue compiled and used in the analysis is a composite catalogue composed of 19,373 events. Ground motion prediction equations (GMPEs) used are calibrated using goodness-of-fitness measures and visual inspection with local strong motion data. Epistemic uncertainty in the GMPEs is taken into account through the logic tree approach. Comparison of ground motions due to deep earthquakes is made for the first time for the region. The comparison between ground motion due to shallow and deep earthquakes indicates that the seismic hazard would be underestimated if the deep earthquakes are excluded. Ground motion values obtained in this study considering all the earthquakes suggest ground motions are dominant towards the north east of the region. The proposed study indicates that the ground motion hazard values suggested by the current Building Code of Pakistan underestimate the seismic hazard. Final results of this study are in close agreement with the recent studies on the region.