Use of Seismotectonic Information for the Seismic Hazard Analysis for Surat City, Gujarat, India: Deterministic and Probabilistic Approach

Surat, the financial capital of Gujarat, India, is a mega city with a population exceeding five millions. The city falls under Zone III of the Seismic Zoning Map of India. After the devastating 2001 Bhuj earthquake of Mw 7.7, much attention is paid towards the seismic microzonation activity in the state of Gujarat. In this work, an attempt has been made to evaluate the seismic hazard for Surat City (21.170?N, 72.830?E) based on the probabilistic and deterministic seismic hazard analysis. After collecting a catalogue of historical earthquakes in a 350?km radius around the city and after analyzing a database statistically, deterministic analysis has been carried out considering known tectonic sources; a further recurrence relationship for the control region is found out. Probabilistic seismic hazard analyses were then carried out for the Surat region considering five seismotectonic sources selected from a deterministic approach. The final results of the present investigations are presented in the form of peak ground acceleration and response spectra at bed rock level considering the local site conditions. Rock level Peak Ground Acceleration (PGA) and spectral acceleration values at 0.01?s and 1.0?s corresponding to 10% and 2% probability of exceedance in 50 years have been calculated. Further Uniform Hazard Response Spectrum (UHRS) at rock level for 5% damping, and 10% and 2% probability of exceedance in 50 years, were also developed for the city considering all site classes. These results can be directly used by engineers as basic inputs in earthquake-resistant design of structures in and around the city.     

Probabilistic seismic hazard assessment of the Pyrenean region

A unified probabilistic seismic hazard assessment (PSHA) for the Pyrenean region has been performed by an international team composed of experts from Spain and France during the Interreg IIIA ISARD project. It is motivated by incoherencies between the seismic hazard zonations of the design codes of France and Spain and by the need for input data to be used to define earthquake scenarios. A great effort was invested in the homogenisation of the input data. All existing seismic data are collected in a database and lead to a unified catalogue using a local magnitude scale. PSHA has been performed using logic trees combined with Monte Carlo simulations to account for both epistemic and aleatory uncertainties. As an alternative to hazard calculation based on seismic sources zone models, a zoneless method is also used to produce a hazard map less dependant on zone boundaries. Two seismogenic source models were defined to take into account the different interpretations existing among specialists. A new regional ground-motion prediction equation based on regional data has been proposed. It was used in combination with published ground-motion prediction equations derived using European and Mediterranean data. The application of this methodology leads to the definition of seismic hazard maps for 475- and 1,975-year return periods for spectral accelerations at periods of 0 (corresponding to peak ground acceleration), 0.1, 0.3, 0.6, 1 and 2 s. Median and percentiles 15% and 85% acceleration contour lines are represented. Finally, the seismic catalogue is used to produce a map of the maximum acceleration expected for comparison with the probabilistic hazard maps. The hazard maps are produced using a grid of 0.1°. The results obtained may be useful for civil protection and risk prevention purposes in France, Spain and Andorra.     

Simulations of Seismic Hazard for the Pacific Northwest of the United States from Earthquakes Associated with the Cascadia Subduction Zone

-- We investigate the impact of different rupture and attenuation models for the Cascadia subduction zone by simulating seismic hazard models for the Pacific Northwest of the U.S. at 2% probability of exceedance in 50 years. We calculate the sensitivity of hazard (probabilistic ground motions) to the source parameters and the attenuation relations for both intraslab and interface earthquakes and present these in the framework of the standard USGS hazard model that includes crustal earthquakes. Our results indicate that allowing the deep intraslab earthquakes to occur anywhere along the subduction zone increases the peak ground acceleration hazard near Portland, Oregon by about 20%. Alternative attenuation relations for deep earthquakes can result in ground motions that differ by a factor of two. The hazard uncertainty for the plate interface and intraslab earthquakes is analyzed through a Monte-Carlo logic tree approach and indicates a seismic hazard exceeding 1 g (0.2 s spectral acceleration) consistent with the U.S. National Seismic Hazard Maps in western Washington, Oregon, and California and an overall coefficient of variation that ranges from 0.1 to 0.4. Sensitivity studies indicate that the paleoseismic chronology and the magnitude of great plate interface earthquakes contribute significantly to the hazard uncertainty estimates for this region. Paleoseismic data indicate that the mean earthquake recurrence interval for great earthquakes is about 500 years and that it has been 300 years since the last great earthquake. We calculate the probability of such a great earthquake along the Cascadia plate interface to be about 14% when considering a time-dependent model and about 10% when considering a time-independent Poisson model during the next 50-year interval.     

Ground motion prediction equation considering combined dataset of recorded and simulated ground motions

Himalayan region is one of the most active seismic regions in the world and many researchers have highlighted the possibility of great seismic event in the near future due to seismic gap. Seismic hazard analysis and microzonation of highly populated places in the region are mandatory in a regional scale. Region specific Ground Motion Predictive Equation (GMPE) is an important input in the seismic hazard analysis for macro- and micro-zonation studies. Few GMPEs developed in India are based on the recorded data and are applicable for a particular range of magnitudes and distances. This paper focuses on the development of a new GMPE for the Himalayan region considering both the recorded and simulated earthquakes of moment magnitude 5.3–8.7. The Finite Fault simulation model has been used for the ground motion simulation considering region specific seismotectonic parameters from the past earthquakes and source models. Simulated acceleration time histories and response spectra are compared with available records. In the absence of a large number of recorded data, simulations have been performed at unavailable locations by adopting Apparent Stations concept. Earthquakes recorded up to 2007 have been used for the development of new GMPE and earthquakes records after 2007 are used to validate new GMPE. Proposed GMPE matched very well with recorded data and also with other highly ranked GMPEs developed elsewhere and applicable for the region. Comparison of response spectra also have shown good agreement with recorded earthquake data. Quantitative analysis of residuals for the proposed GMPE and region specific GMPEs to predict Nepal–India 2011 earthquake of Mw of 5.7 records values shows that the proposed GMPE predicts Peak ground acceleration and spectral acceleration for entire distance and period range with lower percent residual when compared to exiting region specific GMPEs.     

Site specific probabilistic seismic hazard analysis at Dubai Creek on the west coast of UAE

A probabilistic seismic hazard analysis (PSHA) was conducted to establish the hazard spectra for a site located at Dubai Creek on the west coast of the United Arab Emirates (UAE). The PSHA considered all the seismogenic sources that affect the site, including plate boundaries such as the Makran subduction zone, the Zagros fold-thrust region and the transition fault system between them; and local crustal faults in UAE. PSHA indicated that local faults dominate the hazard. The peak ground acceleration (PGA) for the 475-year return period spectrum is 0.17 g and 0.33 g for the 2,475-year return period spectrum. The hazard spectra are then employed to establish rock ground motions using the spectral matching technique.     

Application of the time-predictable model in Peninsular India for future seismic hazard assessment

It has been the belief among Earth scientists that the Peninsular Shield is aseismic, as the region attained stability long ago. However, the earthquake at Koyna (10 December 1967), Bhadrachalam (13 April 1969), Broach (23 March 1970), Hyderabad (30 June 1983), Khillari (30 September 1993), Jabalpur (22 May 1997), Gujarat (26 January 2001), and additional ones of smaller magnitudes, altered this concept. This area has experienced many widely distributed shallow earthquakes, some of them having large magnitudes. It is now widely accepted that seismic activity still continues with moderate events. Therefore, a need has arisen to take into consideration recent seismological data to assess the future seismic status of Peninsular India. Earthquake generation model has been studied to develop the statistical relations with surface wave magnitude (M _S ≥ 4.5). Five seismogenic sources showing clustering of earthquakes and including at least three main shocks of magnitude 4.5 ≤ M _S ≤ 6.5 giving two repeat times, have been identified. It is mainly based on the so-called “regional time-predictable model”. For the considered region it is observed that the time interval between two consecutive main shocks depends on the preceding main shock magnitude (M _p ) not on the following main shocks magnitude M _f suggesting the validity of time predictable model in the region.     

基于特征地震模型含时间的概率地震危险性分析方法及其应用研究

本文介绍了特征地震的对数正态分布模型、 正态分布模型和布朗过程时间模型, 提出了使用地震破裂面源模型的特征地震含时间的概率地震危险性分析理论和方法. 通过具体算例对不同的特征地震模型进行了比较, 并对特征地震危险性分析方法进行了系统探索. 研究结果表明, 特征地震含时间模型在复发周期早期的地震危险性低于不含时间模型, 而在后期其地震危险性则高于不含时间模型. 特征地震复发周期的对数正态分布模型与布朗过程时间模型计算得出的地震危险性差别不大. 在未到期望复发时间时, 正态分布模型与前两种模型计算的地震危险性差别不大; 而接近期望复发时间及之后时段, 正态分布模型计算的地震危险性则迅速增大.      

Probabilistic Assessment of Earthquake Hazard in Gujarat and Adjoining Region of India

The Gujarat and adjoining region falls under all four seismic zones V, IV, III and II of the seismic zoning map of India, and is one of the most seismically prone intracontinental regions of the world. It has experienced two large earthquakes of magnitude M _w 7.8 and 7.7 in 1819 and 2001, respectively and several moderate earthquakes during the past two centuries. In the present study, the probability of occurrence of earthquakes of M ≥ 5.0 has been estimated during a specified time interval for different elapsed times on the basis of observed time intervals between earthquakes using three stochastic models namely, Weibull, Gamma and Lognormal. A complete earthquake catalogue has been used covering the time interval of 1819 to 2006. The whole region has been divided into three major seismic regions (Saurashtra, Mainland Gujarat and Kachchh) on the basis of seismotectonics and geomorphology of the region. The earthquake hazard parameters have been estimated using the method of maximum likelihood. The logarithmic of likelihood function (ln L) is estimated and used to test the suitability of models in three different regions. It was found that the Weibull model fits well with the actual data in Saurashtra and Kachchh regions, whereas Lognormal model fits well in Mainland Gujarat. The mean intervals of occurrence of earthquakes are estimated as 40.455, 20.249 and 13.338 years in the Saurashtra, Mainland Gujarat and Kachchh region, respectively. The estimated cumulative probability (probability that the next earthquake will occur at a time later than some specific time from the last earthquake) for the earthquakes of M ≥ 5.0 reaches 0.9 after about 64 years from the last earthquake (1993) in Saurashtra, about 49 years from the last earthquake (1969) in Mainland Gujarat and about 29 years from the last earthquake (2006) in the Kachchh region. The conditional probability (probability that the next earthquake will occur during some specific time interval after a certain elapsed time from last earthquake) is also estimated and it reaches about 0.8 to 0.9 during the time interval of about 57 to 66 years from the last earthquake (1993) in Saurashtra region, 31 to 51 years from the last earthquake (1969) in Mainland Gujarat and about 21 to 28 years from the last earthquake (2006) in Kachchh region.     

Potential rupture surface model and its application on probabilistic seismic hazard analysis

Potential sources are simplified as point sources or linear sources in current probabilistic seismic hazard analysis (PSHA) methods. Focus size of large earthquakes is considerable, and fault rupture attitudes may have great influ-ence upon the seismic hazard of a site which is near the source. Under this circumstance, it is unreasonable to use the simplified potential source models in the PSHA, so a potential rupture surface model is proposed in this paper. Adopting this model, we analyze the seismic hazard near the Chelungpu fault that generated the Chi-Chi (Jiji) earthquake with magnitude 7.6 and the following conclusions are reached. 1 This model is reasonable on the base of focal mechanism, especially for sites near potential earthquakes with large magnitude; 2 The attitudes of poten-tial rupture surfaces have great influence on the results of probabilistic seismic hazard analysis and seismic zoning.     

Probabilistic seismic hazard analysis: Early history

Probabilistic seismic hazard analysis (PSHA) is the evaluation of annual frequencies of exceedence of ground motion levels (typically designated by peak ground acceleration or by spectral accelerations) at a site. The result of a PSHA is a seismic hazard curve (annual frequency of exceedence vs ground motion amplitude) or a uniform hazard spectrum (spectral amplitude vs structural period, for a fixed annual frequency of exceedence). Analyses of this type were first conceived in the 1960s and have become the basis for the seismic design of engineered facilities ranging from common buildings designed according to building codes to critical facilities such as nuclear power plants. This Historical Note traces the early history of PSHA. Copyright © 2007 John Wiley & Sons, Ltd.     

Application of a Monte‐Carlo simulation approach for the probabilistic assessment of seismic hazard for geographically distributed portfolio

The conventional integral approach is very well established in probabilistic seismic hazard assessment (PSHA). However, Monte‐Carlo (MC) simulations can become an efficient and flexible alternative against conventional PSHA when more complicated factors (e.g. spatial correlation of ground shaking) are involved. This study aims at showing the implementation of MC simulation techniques for computing the annual exceedance rates of dynamic ground‐motion intensity measures (GMIMs) (e.g. peak ground acceleration and spectral acceleration). We use multi‐scale random field technique to incorporate spatial correlation and near‐fault directivity while generating MC simulations to assess the probabilistic seismic hazard of dynamic GMIMs. Our approach is capable of producing conditional hazard curves as well. We show various examples to illustrate the potential use of the proposed procedures in the hazard and risk assessment of geographically distributed structural systems. Copyright © 2015 John Wiley & Sons, Ltd.     

Distant‐earthquake simulations considering source rupture propagation: refining the seismic hazard of Hong Kong

Deterministic and probabilistic seismic hazard analyses should be complementary, in the sense that probabilistic analysis may be used to identify the controlling deterministic design‐level earthquake events, and more sophisticated models of these events may then be developed to account for uncertainties that could not have been included directly in the probabilistic analysis. De‐aggregation of the tentative uniform hazard spectra (UHS) in Hong Kong resulting from a probabilistic seismic hazard assessment (PSHA) indicates that strong and major distant earthquakes, rather than moderate local earthquakes, make the largest contribution to the seismic hazard level within the natural‐period range longer than 0.3 s. Ground‐motion simulations of controlling events located 90 and 340 km from Hong Kong, taking into account uncertainties in the rupture process, reveal that the tentative UHS resulting from the PSHA may have significantly underestimated the mid‐to‐long period components. This is attributed mainly to the adoption of double‐corner source‐spectrum models in the attenuation relationships employed in the PSHA. The results of the simulations indicate clearly that rupture directivity and rupture velocity can significantly affect the characteristics of ground motions, even from such distant earthquakes. The rupture‐directivity effects have profound implications in elongating the second corner period where the constant velocity intersects the constant displacement, thus increasing the associated displacement demand. However, demands for acceleration and velocity are found to be not sensitive to the presence of the directivity pulses. Ground pulses resulting from forward rupture directivity of distant earthquakes have longer predominant periods than the usual near‐fault directivity pulses. These long‐period pulses may have profound implications for metropolises, such as Hong Kong and others in Southeast Asia, having large concentration of high‐rise buildings. Copyright © 2005 John Wiley & Sons, Ltd.     

Earthquake hazard in Marmara Region, Turkey

Earthquake hazard in the Marmara Region, Turkey has been investigated using time-independent probabilistic (simple Poissonian) and time-dependent probabilistic (renewal) models. The study culminated in hazard maps of the Marmara Region depicting peak ground acceleration (PGA) and spectral accelerations (SA)'s at 0.2 and 1 s periods corresponding to 10 and 2% probabilities of exceedance in 50 yrs. The historical seismicity, the tectonic models and the known slip rates along the faults constitute the main data used in the assignment. Based on recent findings it has been possible to provide a fault segmentation model for the Marmara Sea. For the main Marmara Fault this model essentially identifies fault segments for different structural, tectonic and geometrical features and historical earthquake occurrences. The damage distribution and pattern of the historical earthquakes have been carefully correlated with this fault segmentation model. The inter-event time period between characteristic earthquakes in these segments is consistently estimated by dividing the seismic slip estimated from the earthquake catalog by the GPS-derived slip rate of 22±3 mm/yr. The remaining segments in the eastern and southern Marmara region are also identified using recent geological, geophysical studies and historical earthquakes. The model assumes that seismic energy along the segments is released by characteristic earthquakes. For the probabilistic studies characteristic earthquake based recurrence relationships are used. Assuming normal distribution of inter-arrival times of characteristic earthquakes, the ‘mean recurrence time’, ‘covariance’ and the ‘time since last earthquake’ are developed for each segment. For the renewal model, the conditional probability for each fault segment is calculated from the mean recurrence interval of the characteristic earthquake, the elapsed time since the last major earthquake and the exposure period. The probabilities are conditional since they change as a function of the time elapsed since the last earthquake. For the background earthquake activity, a spatially smoothed seismicity is determined for each cell of a grid composed of cells of size 0.005°×0.005°. The ground motions are determined for soft rock (NEHRP B/C boundary) conditions. Western US-based attenuation relationships are utilized, since they show a good correlation with the attenuation characteristics of ground motion in the Marmara region. The possibility, that an event ruptures several fault segments (i.e. cascading), is also taken into account and investigated by two possible models of cascading. Differences between Poissonian and renewal models, and also the effect of cascading have been discussed with the help of PGA ratio maps.     

QLg tomography in Gujarat,Western India

We propose a novel Lg attenuation tomography model (Q_Lg tomography) for the state of Gujarat, Western India, using earthquake data recorded by the Gujarat Seismic Network, operated by the Institute of Seismological Research in Gandhinagar. The waveform dataset consist of 400 3-component recordings, produced by 60 earthquakes with magnitude (M_L) spanning from 3.6 to 5.1, recorded at 60 seismic stations having epicentral distances spanning between 200 and 500 km. Spectral amplitude decays for Lg wave displacement were obtained by generalized inversion at 17 frequencies spanning between 0.9 and 9 Hz. Lg wave propagation efficiency was measured by Lg/Pn spectral ratio categorizing as efficient ratio ≥6 for 86%, intermediate ratio of 3–6 for 10% and inefficient ratio <3 for 4% paths of total 400 ray paths. The earthquake size and quality of waveform recorded at dense network found sufficient to resolve lateral variation of Q_Lg in Gujarat.Average power-law attenuation relationship obtained for Gujarat as Q_Lg(f) = 234f^0.64, which corresponds to high attenuation in comparison to peninsular India shield region and other several regions around the world. Q_Lg tomography resolves the highly attenuating crust of extremely fractured Saurashtra region and tectonically active Kachchh region. The Gujarat average attenuation is also lying in between them. The low attenuation in Cambay and Narmada rift basins and extremely low attenuation in patch of Surendranagar area is identified. This study is the first attempt and can be utilized as pivotal criteria for scenario hazard assessment, as maximum hazard has been reported in highly attenuating tectonically active Kachchh region and in low attenuating Cambay, Narmada and Surendranagar regions. The site and source terms are also obtained along with the Q_Lg inversion. The estimated site responses are comparable with observed local geological condition and agree with the previously reported site amplifications at the same sites. The source terms are comparable with local magnitude estimated from Network. The M_w (Lg) is nearly equivalent to M_L (GSN) and the slight differences are noted for larger magnitude events.     

Uncertainty Analysis and Expert Judgment in Seismic Hazard Analysis

The large uncertainty associated with the prediction of future earthquakes is usually regarded as the main reason for increased hazard estimates which have resulted from some recent large scale probabilistic seismic hazard analysis studies (e.g. the PEGASOS study in Switzerland and the Yucca Mountain study in the USA). It is frequently overlooked that such increased hazard estimates are characteristic for a single specific method of probabilistic seismic hazard analysis (PSHA): the traditional (Cornell?CMcGuire) PSHA method which has found its highest level of sophistication in the SSHAC probability method. Based on a review of the SSHAC probability model and its application in the PEGASOS project, it is shown that the surprising results of recent PSHA studies can be explained to a large extent by the uncertainty model used in traditional PSHA, which deviates from the state of the art in mathematics and risk analysis. This uncertainty model, the Ang?CTang uncertainty model, mixes concepts of decision theory with probabilistic hazard assessment methods leading to an overestimation of uncertainty in comparison to empirical evidence. Although expert knowledge can be a valuable source of scientific information, its incorporation into the SSHAC probability method does not resolve the issue of inflating uncertainties in PSHA results. Other, more data driven, PSHA approaches in use in some European countries are less vulnerable to this effect. The most valuable alternative to traditional PSHA is the direct probabilistic scenario-based approach, which is closely linked with emerging neo-deterministic methods based on waveform modelling.     

Response spectral attenuation relationships for Sumatran-subduction earthquakes and the seismic hazard implications to Singapore and Kuala Lumpur

Singapore and Kuala Lumpur, the capital of Malaysia, may well represent the classic examples of area with low seismic hazard but with high consequence. Both cities are located in a low-seismicity region of Southeast Asia, where active seismic sources are located more than 300 km away. Seismic designs have not been implemented in this seemingly low-hazard region though distant earthquakes in Sumatra had frequently shaken high-rise structures in the two cities. Several studies have been conducted to systematically assess the seismic hazards of Singapore and the Malay Peninsula. The present research particularly addresses issues in deriving a new set of attenuation relationships of peak ground acceleration (PGA), peak ground velocity (PGV) and response spectral acceleration (RSA) for the Sumatran-subduction earthquakes. To be relevant for the seismic hazard assessment of the remote metropolises, the derived attenuation relationships cover a long distance range from 150 to 1500 km. The attenuation relationships are derived using synthetic seismograms that account for source and path effects. The uncertainties in rupture parameters, such as stress drop, strike, dip and rake angles, have been defined according to the regional geological and tectonic settings as well as the ruptures of previous earthquakes. The seismic potential of the Sumatran subduction zone are high in the region from 2°N to 5°S as there has been no recurrence of great thrust events since 1861. A large event with M_w greater than 7.8 in this particular subduction zone may be capable of generating destructive ground motions in Singapore and Kuala Lumpur, even at a distance of 700 km.     

Neo-Deterministic and Probabilistic Seismic Hazard Assessments: a Comparison over the Italian Territory

Estimates of seismic hazard obtained using the neo-deterministic approach (NDSHA) and the probabilistic approach (PSHA) are compared for the Italian territory. The NDSHA provides values larger than those given by the PSHA in areas where large earthquakes are observed and in areas identified as prone to large earthquakes, but lower values in low-seismicity areas. These differences suggest the adoption of the flexible, robust and physically sound NDSHA approach to overcome the proven shortcomings of PSHA, thus allowing for a reliable seismic hazard estimation, especially for those areas characterized by a prolonged quiescence, i.e. in tectonically active sites where events of only moderate size have occurred in historical times.     

Seismic hazard assessment in continental Ecuador

A seismic hazard assessment study of continental Ecuador is presented in this paper. The study begins with a revision of the available information on seismic events and the elaboration of a seismic catalog homogenized to magnitude Mw. Different seismic source definitions are revised and a new area-source model, based on geological and seismic data, is proposed. The available ground motion prediction equations for crustal and subduction sources are analyzed and selected for the tectonic environments observed in Ecuador. A probabilistic seismic hazard assessment approach is carried out to evaluate the exceedance probability of several levels of peak ground acceleration PGA and spectral accelerations SA (T) for periods (T) of 0.1, 0.2, 0.5, 1 and 2s. The resulting hazard maps for continental Ecuador are presented, together with the uniform hazard spectra of four province capital cities. Hazard disaggregation is carried out for target motions defined by the PGA values and SA (1s) expected for return periods of 475 and 2475 years, providing estimates for short-period and long-period controlling earthquakes.     

累积绝对速度在核电厂地震危险性分析中的应用研究

指出了运用我国现行的考虑时空非均匀性的地震危险性分析计算方法对核电厂等设计精良的设施进行地震危险性分析时所存在的问题.介绍了累积绝对速度(CAV)的概念,并将其引入到我国现行的考虑时空非均匀性的地震危险性分析计算方法之中,用以排除厂址周围小震对核电厂地震危险性分析的影响,并选取实际工程场点进行了试算.试算结果表明,此方法能明显排除厂址周围小震对地震危险性分析结果的影响.     

新版地震区划图地震活动性模型与参数确定

地震活动性模型和地震动预测模型是概率地震危险性分析的两个核心。在新版地震区划图中,依据板内地震活动空间不均匀性分布的特点,在概率地震危险性分析方法(CPSHA)中采用了由地震统计区、背景潜在震源区和构造潜在震源区构成的三级层次性潜在震源区模型,并构建了相应的地震活动性模型。本文在论述CPSHA方法及其地震活动性模型基本概念的基础上,重点介绍了新版地震区划图地震活动性模型的三级潜在震源区模型的构成、地震活动性假定和基本特点,同时,也对新版地震区划图地震活动性模型的重要参数确定思路、方法与结果进行了介绍。本文将为更好地认识与理解我国新版地震动参数区划图提供有益的参考。