Non-Poisson probabilistic seismic hazard assessment

The development of the new seismic hazard map of metropolitan Tehran is based on probabilistic seismic hazard computation using the non-Poisson recurrence time model. For this model, two maps have been prepared to indicate the earthquake hazard of the region in the form of iso-acceleration contour lines. They display the non-Poisson probabilistic estimates of peak ground accelerations over bedrock for 10 and 63 % probability of exceedance in 50 years. To carry out the non-Poisson seismic hazard analysis, appropriate distributions of interoccurrence times of earthquakes were used for the seismotectonic provinces which the study region is located and then the renewal process was applied. In order to calculate the seismic hazard for different return periods in the probabilistic procedure, the study area encompassed by the 49.5–54.5°E longitudes and 34–37°N latitudes was divided into 0.1° intervals generating 1,350 grid points. PGA values for this region are estimated to be 0.30–0.32 and 0.16–0.17 g for 10 and 63 % probability of exceedance, respectively, in 50 years for bedrock condition.     

Probabilistic seismic hazard assessment in the northeastern part of Algeria

This work involves updating the evaluation of seismic hazard in Northeast Algeria by a probabilistic approach. This reassessment attempts to resolve inconsistencies between seismic zoning in regional building codes and is further motivated by the need to refine the input data that are used to evaluate seismic hazard scenarios. We adopted a seismotectonic model that accounts for differences in interpretations of regional seismicity. We then performed a probabilistic assessment of regional seismic hazard in Northeast Algeria. Based on a homogeneous earthquake catalog and geological and seismotectonic data gathered in the first part of the study, a seismotectonic zoning map was created and seven risk areas were identified. For each area, peak ground acceleration hazard maps were produced. Details of the calculations are provided, including hazard curves at periods of 0.1, 0.2, 0.33, 0.5, 1.0, and 2.0 s and uniform hazard spectra at urban locations in the area, including Sétif, Constantine, Kherrata, Bejaia, and Jijel.     

Deterministic and probabilistic seismic hazard analysis for Gwadar City, Pakistan

Gwadar City is located at the coastline of Pakistan. The city is currently in a phase of development, which is expected to become a future economic hub for Pakistan. This has led us to choose Gwadar for seismic hazard evaluation. Seismic hazard analysis for Gwadar is carried out using deterministic and probabilistic seismic hazard analysis techniques. The present study will help in sustainable development of a future large city and economic hub for Pakistan on ways of coping from a major threat of earthquake hazard. In deterministic seismic hazard analysis, line sources were identified close to Gwadar. Based on the analysis of maximum magnitude and closest distance (worse conditions), Makran subduction zone was identified out of all the line sources with earthquake potential of 8.2 at a distance of 30 km. This yielded a peak ground acceleration value of 0.38 g for Gwadar City. In second phase, probabilistic seismic hazard analysis technique with the area source modeling was adopted to acquire results at different return periods. For this purpose, seismic data were collected from the Pakistan Meteorological Department and International Seismological Center (2010) databases for development of a comprehensive data catalog. The a and b values were obtained using regression analysis for each source zone, and probabilistic analysis yielded the results of 0.34 g for a return period of 500 years. As per building codes of Pakistan, areas or cities with ground acceleration greater than 0.32 g are considered in seismic zone 4, and both deterministic and probabilistic hazard analysis place the city in seismic zone 4. These values correspond to rock site with shear wave velocity of 760 m/s.     

Near-field probabilistic seismic hazard analysis with characteristic earthquake effects

In conventional seismic hazard analysis, uniform distribution over area and magnitude range is assumed for the evaluation of source seismicity which is not able to capture peculiar characteristic of near-fault ground motion well. For near-field hazard analysis, two important factors need to be considered: (1) rupture directivity effects and (2) occurrence of scenario characteristic ruptures in the nearby sources. This study proposed a simple framework to consider these two effects by modifying the predictions from the conventional ground motion model based on pulse occurrence probability and adjustment of the magnitude frequency distribution to account for the rupture characteristic of the fault. The results of proposed approach are compared with those of deterministic and probabilistic seismic hazard analyses. The results indicate that characteristic earthquake and directivity consideration both have significant effects on seismic hazard analysis estimates. The implemented approach leads to results close to deterministic seismic hazard analysis in the short period ranges (T < 1.0 s) and follows probabilistic seismic hazard analysis results in the long period ranges (T > 1.0 s). Finally, seismic hazard maps based on the proposed method could be developed and compared with other methods.     

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.     

Probabilistic seismic hazard analysis for Kakrapar atomic power station, Gujarat, India

Seismic ground motion caused by earthquakes mainly affects the constructions and structures around its area of influence. In this context, the probabilistic seismic hazard analysis (PSHA) is a scientific step towards the safety analysis of any major construction such as nuclear power plant. Thus, the present study focused to estimate seismic hazard level at different probabilities for Kakrapar nuclear power plant located in the Western India. The hazard curves for the study area are developed following the procedure of PSHA suggested by Cornell–McGuire. Three source zones, Narmada-Tapti zone (NTZ), Rann of Kuchchh (ROK), and west passive margin (WPM), are classified on the basis of seismicity and tectonic setting of the study area. The estimated maximum magnitude (m _max) for NTZ, ROK, and WPM are 6.9 ± 0.57, 6.5 ± 0.64, and 6.1 ± 0.64, respectively. Logic tree approach has been used for the development of hazard curves to account the epistemic uncertainties associated with the analysis. For maximum credible earthquake [MCE, i.e., the probability of exceedance of 2 % in 50 years (return period of ~2,500 years)], the peak spectral acceleration (i.e., PSA at 0.2 s) expected around 5 km of the Kakrapar nuclear power plant (site) is 0.23 g from all source zones; however, at exact site location, it is 0.18 g. The PSA values due to NTZ, ROK, and WPM based on MCE are 0.22, 0.065, and 0.052 g, respectively. In case of design-based earthquake (DBE, i.e., 50 % probability in 50 years (return period of ~110 years)), the calculated maximum spectral acceleration (SA) from all source zones is about 0.045 g. The PSA distribution for the DBE from the NTZ has reached a maximum value of 0.042 g; however, PSA for ROK and WPM is considerably low with a maximum value of 0.022 and 0.021 g, respectively. Considering the MCE and DBE, the estimated PSA at 0.2 s has a highest value of ~0.23 g from all source zones. Spectral accelerations (SAs) correspond to different periods are presented, and SA plots for NTZ zone can be considered as response spectra for the KAPS site. Deaggregation of PSHA in the present study is also discussed. PGA values reported in seismic zonation map and global seismic hazard analysis program around the present study area range from 0.05 to 0.2 g which is slightly lower than the peak acceleration obtained in this study. The results of this study would facilitate in the performance of the site-specific seismic probabilistic safety analysis.     

Seismic hazard analysis of India using areal sources

In view of the major advancement made in understanding the seismicity and seismotectonics of the Indian region in recent times, an updated probabilistic seismic hazard map of India covering 6–38°N and 68–98°E is prepared. This paper presents the results of probabilistic seismic hazard analysis of India done using regional seismic source zones and four well recognized attenuation relations considering varied tectonic provinces in the region. The study area was divided into small grids of size 0.1° × 0.1°. Peak Horizontal Acceleration (PHA) and spectral accelerations for periods 0.1 s and 1 s have been estimated and contour maps showing the spatial variation of the same are presented in the paper. The present study shows that the seismic hazard is moderate in peninsular shield, but the hazard in most parts of North and Northeast India is high.     

Assessment of drought hazard impact on wheat cultivation using standardized precipitation index in Iran

The assessment of drought hazard impacts on wheat cultivation as a strategic crop in Iran is essential for making mitigation plans to reduce the impact of drought. Standardized precipitation index has gained importance in recent years as a potential drought indicator and is being used more frequently for assessment of drought hazard in many countries. In the present study, the calculated standardized precipitation index for 48 stations dataset in the 30-year time scale fulfilled 30 statistical matrices. The drought hazard index map was produced by sum overlaying the spatial representations of 30 statistical matrices and categorized into four levels of low, moderate, high, and very high, which demonstrated probability of drought occurrences of 10–20 %, 20–30 %, 30–40 %, and 40–50 %, respectively. Finally, after the general division of zonal statistics in drought hazard index map of Iran, major drought hazard zones were geographically classified into five zones. The statistical analysis showed a significant correlation (R ²?=?0.701 to 0.648) between drought occurrences and wheat cultivation including surface area and total production for these drought hazard zones.     

The New Seismic Zonation Map for the Territory of Armenia

In the present work a new, deterministic, seismic zonation map for the territory of Armenia is presented. It has been compiled on the basis of the newest seismological, geological and geophisical data. Creation of the map included: determination of a seismic hazard area for the territory of Armenia; compilation of a map of active faults; identification of active blocks; identification of seismic source zones; calculation of a seismic effect from the seismic source zones (SSZs). Seismic effect on the Earth's surface from the singled out SSZs is calculated in the form of an intensity of seismic influences expressed in units of MSK-64 scale, and horizontal soil accelerations expressed in fractions of gravity force acceleration (g). The map compiled will allow comparison between the deterministic model with a probabilistic seismic zonation model for the territory of Armenia (in preparation), in order to analyze the western and eastern standards of seismic hazard assessment on the basis of the same complete data base.     

The comparison of the NDSHA, PSHA seismic hazard maps and real seismicity for the Italian territory

Rigorous and objective testing of seismic hazard assessments against the real seismic activity must become the necessary precondition for any responsible seismic risk estimation. Because seismic hazard maps seek to predict the shaking that would actually occur, the reference hazard maps for the Italian seismic code, obtained by probabilistic seismic hazard assessment (PSHA), and the alternative ground shaking maps based on the neo-deterministic approach (NDSHA), are cross-compared and tested against the real seismicity for the territory of Italy. The comparison between predicted intensities and those reported for past earthquakes shows that models generally provide rather conservative estimates, except for PGA with 10 % probability of being exceeded in 50 years, which underestimates the largest earthquakes. In terms of efficiency in predicting ground shaking, measured accounting for the rate of underestimated events and for the territorial extent of areas characterized by high seismic hazard, the NDSHA maps appear to outscore the PSHA ones.     

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.     

Site-specific probabilistic seismic hazard analysis for northern part of the Qeshm Island,Iran

This paper presents results of a site-specific probabilistic seismic hazard analysis for northern part of the Qeshm Island, one the most seismic prone areas of Iran. Seismotectonic and seismicity properties of seismic sources in the study area were characterized and used for evaluation of various strong ground motion parameters implementing the classical Cornell’s PSHA approach. The results show that peak rock accelerations for 475-year return period are 0.4 and 0.27 g, respectively, for 84th and 50th percentiles while being about 0.37 and 0.61 g for 2475-year return period. These values are slightly smaller than those read from national seismic zonation maps which can be attributed to the considered conservatism for development of such design maps. In order to incorporate local site conditions, a series of dynamic site response analyses based on the equivalent linear approach were also employed. The results indicate that the presence of soft subsurface deposits at the site significantly alters the fundamental characteristics of the response spectra. The obtained median (50th percentile) peak ground accelerations for 975-year return period range between 0.49 and 0.54 g at different locations in the study site showing minor amplifications relative to their corresponding bedrock acceleration of 0.48 g. Finally, the obtained site-specific spectrum was compared with the standard spectrum mandated by the design codes. In this regard, the agreement was found to be reasonable at period ranges shorter than about 0.5 s, while the differences were more obvious at longer periods. This reveals the need for implementation of site-specific design spectrum to avoid underestimation or overestimation of seismic forces for designing critically important structures especially when softer deposits are encountered.     

Seismic hazard studies for Gaziantep city in South Anatolia of Turkey

Seismic hazard studies were conducted for Gaziantep city in the South Anatolia of Turkey. For this purpose, a new attenuation relationship was developed using the data of Zaré and Bard and accelerations were predicted employing this new equation. Deterministic approach, total probability theorem and GIS methodology were all together utilized for the seismic assessments. Seismic hazard maps with 0.25° grid intervals considering the site conditions were produced by the GIS technique. The results indicated that the acceleration values by the GIS hazard modelings were matched with the ones from the deterministic approach, however, they were underestimated comparing with the total probability theorem. In addition, the GIS based seismic hazard maps showed that the current seismic map of Turkey fairly yields conservative acceleration values for the Gaziantep region. Therefore, the constructed GIS hazard models are offered as a base map for a further modification of the current seismic hazard map.     

Regional seismic hazard study for the eastern Mediterranean (Trans-Jordan,Levant and Antakia) and Sinai region

The aim of this study is to assess the seismic hazard in the eastern Mediterranean and Sinai region using a probabilistic approach. An updated earthquake catalogue for the period 1 to 1993 AD that covers the area between latitude 27°–37°N and longitude 32°–39°E, has been used. Using the new seismic-tectonic map for the area, 10 line-sources are delineated. These lines or fault zones are thought to represent the main sources for the seismic potential in the area. The results are demonstrated as iso-contour lines of the peak-ground acceleration. The iso-acceleration contours represent 90% probability that these peak values will not be exceeded over periods of 50, 100 and 200 years, respectively. This study concludes that the seismic hazard severity is highest along the Jordan Dead Sea transform fault system, namely from south of the Gulf of Aqaba, Dead Sea-Jordan River, Tiberia Lake, Rachaya, Ed Damur, Yammuneh Fault, and Ghab Fault in the north. For the 50 year iso-contour map, the major cities of Amman, Damascus, and Beirut lay around the 2 m s⁻² contour line, while Jerusalem lies along the 3 m s⁻² line. Antakia in Turkey has the highest seismic potential severity (around 5 m s⁻²) while in Cyprus the maximum hazard is expected to reach 4 m s⁻² for the coming 50 years.     

Probabilistic seismic hazard map of NW Himalaya and its adjoining area,India

The seismically active Northwest (NW) Himalaya falls within Seismic Zone IV and V of the hazard zonation map of India. The region has suffered several moderate (~25), large-to-great earthquakes (~4) since Assam earthquake of 1897. In view of the major advancement made in understanding the seismicity and seismotectonics of this region during the last two decades, an updated probabilistic seismic hazard map of NW Himalaya and its adjoining areas covering 28–34°N and 74–82°E is prepared. The northwest Himalaya and its adjoining area is divided into nineteen different seismogenic source zones; and two different region-specific attenuation relationships have been used for seismic hazard assessment. The peak ground acceleration (PGA) estimated for 10% probability of exceedance in 50 and 10 years at locations defined in the grid of 0.25 × 0.25°. The computed seismic hazard map reveals longitudinal variation in hazard level along the NW Himalayan arc. The high hazard potential zones are centred around Kashmir region (0.70 g/0.35 g), Kangra region (0.50 g/0.020 g), Kaurik-Spitti region (0.45 g/0.20 g), Garhwal region (0.50 g/0.20 g) and Darchula region (0.50 g/0.20 g) with intervening low hazard area of the order of 0.25 g/0.02 g for 10% probability in 50 and 10 years in each region respectively.     

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.     

Estimation of Seismic Landslide Hazard in the Eastern Himalayan Syntaxis Region of Tibetan Plateau

The eastern Himalayan syntaxis is one of the most tectonically active and earthquake-prone regions on Earth where earthquake-induced geological disasters occur frequently and caused great damages. With the planning and construction of Sichuan-Tibet highway, Sichuan-Tibet railway and hydropower development on the Yarlung Zangbo River etc. in recent years, it is very important to evaluate the seismic landslide hazard of this region. In this paper, a seismic landslide hazard map is produced based on seismic geological background analysis and field investigation using Newmark method with 10% PGA exceedance probabilities in future 50 years by considering the influence of river erosion, active faults and seismic amplification for the first time. The results show that the areas prone to seismic landslides are distributed on steep slopes along the drainages and the glacier horns as well as ridges on the mountains. The seismic landslide hazard map produced in this study not only predicts the most prone seismic landslide areas in the future 50 years but also provides a reference for mitigation strategies to reduce the exposure of the new building and planning projects to seismic landslides.     

PSHA of Van province for performance assessment using spectrally matched strong ground motion records

Within the framework of the performance based earthquake engineering, site specific earthquake spectra for Van province has been obtained. It is noteworthy that, in probabilistic seismic hazard assessment, as a first stage data from geological studies and records from the instrumental period were compiled to make a seismic source characterization for the study region. The probabilistic seismic hazard curves were developed based on selected appropriate attenuation relationships, at rock sites, with a probability of exceedance 2, 10 and 50% in 50 yrs period. The obtained results are compared with the spectral responses proposed for seismic evaluation and retrofit of building structure in Turkish Earthquake Code (2007), section 7. The acceleration response spectrums obtained from probabilistic seismic hazard analysis are matched to adjust earthquake accelerograms recorded during the 2011 Van earthquakes by using SeismoMatch v2.0 software. The aim of this procedure is to obtain a set of reasonable earthquake input motions for the seismic evaluation of existing buildings.     

Spatial variation of probabilistic seismic hazard for Mumbai and surrounding region

Mumbai city, the economical capital of India, is located on the west coast of stable intra-plate continental region of Peninsular India which has an experience of significant historical earthquakes in the past. The city stood as the fourth most populous city in the world. Recent seismo-tectonic studies of this city highlighted the presence of active West coast fault and Chiplun fault beneath the Deccan basalt. In the present study, spatial variability of probabilistic seismic hazard for Mumbai region (latitudes of 18.85–19.35°N and longitudes of 72.80–73.15°E at a grid spacing of 0.05°) which includes Mumbai city, Suburban, part of Thane district and Navi Mumbai, in terms of ground motion parameters; peak horizontal acceleration and spectral acceleration at 1.0-s period for 2 and 10 % probability of exceedance in 50 years are generated. The epistemic uncertainty in hazard estimation is accounted by employing seven different ground motion prediction equations developed for worldwide shallow crustal intra-plate environments. Further, the seismic hazard results are deaggregated for Mumbai (latitude 18.94°N, longitude 72.84°E) to understand the relative contributions of earthquake sources in terms of magnitude and distance. The generated hazard maps are compared with the zoning specified by Indian seismic code (IS1893: Part 1 in Indian standard criteria for earthquake-resistant design of structures, Part 1—General provisions and buildings. Bureau of Indian Standards, New Delhi, India, 2002) for rocky site. Present results show an underestimation of potential seismic hazard in the entire study region by non-probabilistic zoning prescribed by IS1893: Part 1 with significantly higher seismic hazard values in the southern part of Navi Mumbai.     

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.