Probabilistic seismic hazard analysis for Bangalore

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

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.     

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.     

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.     

Long-term prediction for seismic hazard for radioactive waste disposal

We consider possible approaches to the long-term prediction for seismic hazard in relation to the practical need for the safety of geological disposal of long-lived radioactive waste. The required period of prediction significantly exceeds the one reflected in the set of maps of General Seismic Zoning of the territory of the Russian Federation (GSZ-97). The first geological repository in Russia is planned to be set up in the Nizhnii Kan granite massif in the Krasnoyarsk Krai. This region is an intraplate territory with a relatively high seismic activity. We summarize the analysis of the known empirical generalizations and theoretical principles underlying the seismic hazard prediction. Real seismic events constantly violate forward-looking statements even for relatively short periods of time. These and other arguments suggest that the hypothesis of stationarity of the seismic regime, which is the basis of long-term prediction today, has limited and uncertain applicability in time. Intraplate earthquake prediction is especially uncertain because of the uncertainty in the factor responsible for generating tectonic stresses in these regions. The short horizon of the prediction, based on statistical methods, can be attributed to the nonlinearity of seismic geodynamic processes. Fundamental laws of tectonic processes should be used as the scientific basis for long-term predictions for seismic hazard at the sites chosen for geological disposal of long-lived radioactive waste. These processes can be reflected in models for the migration of the seismically active boundaries of lithospheric plates and the occurrence of seismic activity in intraplate regions.     

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.     

Bayesian techniques for seismic hazard assessment using imprecise data

Some Bayesian methods of dealing with inaccurate or vague data are introduced in the framework of seismic hazard assessment. Inaccurate data affected by heterogeneous errors are modeled by a probability distribution instead of the usual value plus a random error representation; these data are generically called imprecise. The earthquake size and the number of events in a certain time are modeled as imprecise data. Imprecise data allow us to introduce into the estimation procedures the uncertainty inherent in the inaccuracy and heterogeneity of the measuring systems from which the data were obtained. The problem of estimating the parameter of a Poisson process is shown to be feasible by the use of Bayesian techniques and imprecise data. This background technique can be applied to a general problem of seismic hazard estimation. Initially, data in a regional earthquake catalog are assumed imprecise both in size and location (i.e errors in the epicenter or spreading over a given source). By means of scattered attenuation laws, the regional catalog can be translated into a so-called site catalog of imprecise events. The site catalog is then used to estimate return periods or occurrence probabilities, taking into account all sources of uncertainty. Special attention is paid to priors in the Bayesian estimation. They can be used to introduce additional information as well as scattered frequency-size laws for local events. A simple example is presented to illustrate the capabilities of this methodology.     

An indicator kriging model for investigation of seismic hazard

Time domain probabilistic techniques most often are used for assessment of seismic hazard. Such techniques are based on the historic frequency of ground motion. Hazard is expressed as a probability of experiencing a particular level of seismic activity over a given length of time. One of these techniques utilizes frequency of extreme values for assessment of hazard. The major disadvantage of this technique, however, becomes evident when maximum seismic activity for two consecutive years occurs only a few weeks or months apart. In this case, the extreme value approach overestimates seismic hazard. A new approach for hazard assessment is founded on principles of indicator kriging. This technique evaluates seismic hazard as a simple frequency record, which is more realistic for regions of little to moderate seismicity.     

Attenuation laws and seismic hazard assessment

In the present study, the seismic hazard in the city of Patras, central Greece, is estimated. The computations are based on a slightly modified version of the method proposed by Cornell, allowing the introduction of individual attenuation laws for each seismic source.The obtained results emphasize the dependence of hazard on attenuation and the importance of the use of local attenuation laws in seismic hazard assessment.Paper presented at the 21st General Assembly of the European Seismological Commission, held in Sofia, 1988.     

Attenuation laws and seismic hazard assessment

In the present study, the seismic hazard in the city of Patras, central Greece, is estimated. The computations are based on a slightly modified version of the method proposed by Cornell, allowing the introduction of individual attenuation laws for each seismic source. The obtained results emphasize the dependence of hazard on attenuation and the importance of the use of local attenuation laws in seismic hazard assessment.     

Re-evaluations of seismic hazard of Iraq

In recent years, Iraq has experienced an increase in seismic activity, especially, near the east boundary with Iran. Previous studies present their results in terms of PGA and for return periods of 500 years and less, and other studies not continued to include the whole PSHA process whereas some recent studies continued to include the whole PSHA process using earthquakes data till 2009 including dependent events. This study includes two main stages, the first is collecting the earthquakes records including the recent events till the end of March 2016 and applying data processing to get the net catalog to independent events. The second stage is applying the steps of PSHA method. Matlab programs have been built to execute these two stages and to convert the results of PSHA computations into contours of 5% damping PGA and spectral accelerations at 0.2 and 1.0 s for a return period of 2475 years, and for rock sites. Also, spectral acceleration against period has been presented for main cities. Also, the PGA map, for a return period of 475 years, has been plotted and then prepared together with similar maps of neighbor countries in one map for comparison. In general, this comparison indicates the similarity in behavior but, the values reveal a relative agreement and they are between Turkish and Iranian values.     

Regionalization of seismic hazard in Greece based on seismic sources

A semi-probabilistic approach to the seismic hazard assessment of Greece is presented. For this reason, a recent seismotectonic model for shallow and intermediate depth earthquake sources, based on historical as well as on instrumental data, was used. Different attenuation formulae were proposed for the macroseismic intensity and the strong ground motion parameters for the shallow and the intermediate focal depth shocks. The data were elaborated in terms of McGuire's computer program, which is based on the Cornell's method.A grid of equally spaced points at 20 km distance was made and the seismic hazard recurrence curves for various parameters of the seismic intensity was estimated for each point. Finally, seismic hazard maps for the area of Greece were compiled utilizing the entire range of recurrence curves. These maps depict areas of equal seismic hazard and for every area the analytical relations of the typeSI =f(Tm), whereSI is a seismic intensity parameter andTm is the mean return period, were determined.     

Preliminary seismic hazard maps of Slovenia

The preparation of the preliminary seismic hazard maps of the territory of Slovenia has been based on an expansion of the basic approach laid out by Cornell in 1968. Three seismic source models were prepared. Two of them are based mainly on the earthquake catalogue using the Poissonian probability model. A map of seismic energy release and a map of earthquake epicenter density are used to delineate seismic sources in these models. The geometry of the third model which is based on a rough estimate of seismotectonic setting is taken from the probabilistic seismic hazard analysis of a nuclear power plant in Slovenia. Published ground motion attenuation models based on strong motion records of recent strong earthquakes in Italy are used. Test maps for variable and uniform b-values are presented. The computer program, Seisrisk III, developed by the U.S. Geological Survey is used.     

Probabilistic seismic hazard zonation of Syria

Earthquake hazard maps for Syria are presented in this paper. The Peak Ground Acceleration (PGA) and the Modified Mercalli Intensity (MMI) on bedrock, both with 90% probability of not being exceeded during a life time of 50, 100 and 200 years, respectively are developed. The probabilistic PGA and MMI values are evaluated assuming linear sources (faults) as potential sources of future earthquakes. A new attenuation relationship for this region is developed. Ten distinctive faults of potential earthquakes are identified in and around Syria. The pertinent parameters of each fault, such as theb-parameter in the Gutenberg-Richter formula, the annual rate ₄ and the upper bound magnitudem ₁ are determined from two sets of seismic data: the historical earthquakes and the instrumentally recorded earthquake data (AD 1900–1992). The seismic hazard maps developed are intended for preliminary analysis of new designs and seismic check of existing civil engineering structures.     

A seismic source zone model for the seismic hazard assessment of the Italian territory

We designed a new seismic source model for Italy to be used as an input for country-wide probabilistic seismic hazard assessment (PSHA) in the frame of the compilation of a new national reference map.

We started off by reviewing existing models available for Italy and for other European countries, then discussed the main open issues in the current practice of seismogenic zoning.

The new model, termed ZS9, is largely based on data collected in the past 10 years, including historical earthquakes and instrumental seismicity, active faults and their seismogenic potential, and seismotectonic evidence from recent earthquakes. This information allowed us to propose new interpretations for poorly understood areas where the new data are in conflict with assumptions made in designing the previous and widely used model ZS4.

ZS9 is made out of 36 zones where earthquakes with M_w > = 5 are expected. It also assumes that earthquakes with M_w up to 5 may occur anywhere outside the seismogenic zones, although the associated probability is rather low. Special care was taken to ensure that each zone sampled a large enough number of earthquakes so that we could compute reliable earthquake production rates.

Although it was drawn following criteria that are standard practice in PSHA, ZS9 is also innovative in that every zone is characterised also by its mean seismogenic depth (the depth of the crustal volume that will presumably release future earthquakes) and predominant focal mechanism (their most likely rupture mechanism). These properties were determined using instrumental data, and only in a limited number of cases we resorted to geologic constraints and expert judgment to cope with lack of data or conflicting indications. These attributes allow ZS9 to be used with more accurate regionalized depth-dependent attenuation relations, and are ultimately expected to increase significantly the reliability of seismic hazard estimates.     

Integrated seismic hazard evaluation and disaster management approach for Turkey

A significant proportion of the urban areas in Turkey is subject to high seismic risk. An important step for seismic risk mitigation is to define the hazard and damage after an earthquake. This paper proposes an integrated seismic hazard assessment and disaster management processes for Turkey. The proposed methodology utilizes information technologies in its seismic assessment component that provides fast results for assessment. First, image process methodology by using satellite images was implemented in the seismic assessment process for fast evaluation right after an earthquake. Second, the seismic assessment process was integrated with disaster management process. As a result, through integrated seismic hazard evaluation and disaster management procedure, an effective, fast and dependable estimation of loss for Turkey was planned.     

Modeling and assessment of bridge structure for seismic hazard prevention

This study uses data from Mao-Luo-Hsi Bridge to model the bridge structure and a set of developed alarm and action values to formulate guidelines for bridge maintenance and seismic hazard prevention. The bridge model is improved by incorporating on-site ambient vibration measurement to perform modal analyses. Dynamic analyses of the bridge are implemented using the established 3D model subjected to uniform loading and seismic force, with or without consideration of soil interaction with the structure. The maximum displacements for different sections of the bridge are compared, and statistical regression analyses are used to explore their correlation. Information for bridge safety assessment is proposed, which can mitigate loss of property and lives due to bridge failure. Regression analyses of the maximum displacements between abutments D and E of the considered bridge in the axial, horizontal, and vertical directions under various seismic intensities are conducted, giving R ² values of 0.9462, 0.9352, and 0.9010, respectively. The developed maintenance guidelines are reliable since all parameters from regression analyses have a 95% confidence interval excluding the zero value. The bridge alarm value and action value are determined for this bridge site at earthquake intensity scales of 4 and 5, respectively.     

A method for producing digital probabilistic seismic landslide hazard maps

The 1994 Northridge, California, earthquake is the first earthquake for which we have all of the data sets needed to conduct a rigorous regional analysis of seismic slope instability. These data sets include: (1) a comprehensive inventory of triggered landslides, (2) about 200 strong-motion records of the mainshock, (3) 1:24 000-scale geologic mapping of the region, (4) extensive data on engineering properties of geologic units, and (5) high-resolution digital elevation models of the topography. All of these data sets have been digitized and rasterized at 10 m grid spacing using ARC/INFO GIS software on a UNIX computer. Combining these data sets in a dynamic model based on Newmark's permanent-deformation (sliding-block) analysis yields estimates of coseismic landslide displacement in each grid cell from the Northridge earthquake. The modeled displacements are then compared with the digital inventory of landslides triggered by the Northridge earthquake to construct a probability curve relating predicted displacement to probability of failure. This probability function can be applied to predict and map the spatial variability in failure probability in any ground-shaking conditions of interest. We anticipate that this mapping procedure will be used to construct seismic landslide hazard maps that will assist in emergency preparedness planning and in making rational decisions regarding development and construction in areas susceptible to seismic slope failure.