Probabilistic Assessment of Earthquake Recurrence in Northeast India and Adjoining Regions

Northeast India and adjoining regions (20°–32° N and 87°–100° E) are highly vulnerable to earthquake hazard in the Indian sub-continent, which fall under seismic zones V, IV and III in the seismic zoning map of India with magnitudes M exceeding 8, 7 and 6, respectively. It has experienced two devastating earthquakes, namely, the Shillong Plateau earthquake of June 12, 1897 (M _w 8.1) and the Assam earthquake of August 15, 1950 (M _w 8.5) that caused huge loss of lives and property in the Indian sub-continent. In the present study, the probabilities of the occurrences of earthquakes with magnitude M ≥ 7.0 during a specified interval of time has been estimated on the basis of three probabilistic models, namely, Weibull, Gamma and Lognormal, with the help of the earthquake catalogue spanning the period 1846 to 1995. The method of maximum likelihood has been used to estimate the earthquake hazard parameters. The logarithmic probability of likelihood function (ln L) is estimated and used to compare the suitability of models and it was found that the Gamma model fits best with the actual data. The sample mean interval of occurrence of such earthquakes is estimated as 7.82 years in the northeast India region and the expected mean values for Weibull, Gamma and Lognormal distributions are estimated as 7.837, 7.820 and 8.269 years, respectively. The estimated cumulative probability for an earthquake M ≥ 7.0 reaches 0.8 after about 15–16 (2010–2011) years and 0.9 after about 18–20 (2013–2015) years from the occurrence of the last earthquake (1995) in the region. The estimated conditional probability also reaches 0.8 to 0.9 after about 13–17 (2008–2012) years in the considered region for an earthquake M ≥ 7.0 when the elapsed time is zero years. However, the conditional probability reaches 0.8 to 0.9 after about 9–13 (2018–2022) years for earthquake M ≥ 7.0 when the elapsed time is 14 years (i.e. 2009).     

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.     

印度及其毗邻地区的概率地震危险性图

本文介绍了全球地震危险性评估计划 ( GSHAP)在 0°～ 40°N和 65°～1 0 0°E范围内的印度及其毗邻地区的实施情况。通过对包括美国国家海洋及大气管理局 ( NOAA)汇编的地震目录在内的地方目录进行筛选 ,删除重复的 ,余震和没有震级的地震 ,将主要地震汇编成地震目录。根据主要的构造特征和地震活动性趋势 ,勾画出 86个潜在震源区。运用GSHAP所采用的 Mc Guire概率性地震评估法 ,计算出今后 50年内由 0 .5°× 0 .5°的网格确定的位置上、超越概率为 1 0 %的峰值地面加速度 ( PGA)。由于目前还没有适用于印度地区对衰减值估算的可靠的方法 ,所以我们还是采用乔伊纳和布尔( Joyner and Boore,1 981 )的衰减关系式。将格点上的 PGA值制成等值图来得到地震危险性图。危险性图反映出印度板块边界地区和青藏高原地区的危险性水平是0 .2 5g量级 ,在地震活动比较活跃的地带 ,如印度东北部、兴都库什地区和缅甸地区 ,此危险性水平是 0 .3 5～ 0 .4g量级。在印度地盾 ,在大部分地区的区域地震危险性量级大约是 0 .0 5～ 0 .1 g,而有些地区 ,如科伊纳地区的地震危险性水平是0 .2 g量级。该图还可改绘成分别标有因子 0 .1 g、 0 .2 g、 0 .3 g和 0 .4g四个地带的常规地震区划图     

Assessment of the Relative Largest Earthquake Hazard Level in the NW Himalaya and its Adjacent Region

In the present study, the level of the largest earthquake hazard is assessed in 28 seismic zones of the NW Himalaya and its vicinity, which is a highly seismically active region of the world. Gumbel’s third asymptotic distribution (hereafter as GIII) is adopted for the evaluation of the largest earthquake magnitudes in these seismic zones. Instead of taking in account any type of M_max, in the present study we consider the ω value which is the largest earthquake magnitude that a region can experience according to the GIII statistics. A function of the form Θ(ω, RP_6.0) is providing in this way a relatively largest earthquake hazard scale defined by the letter K (K index). The return periods for the ω values (earthquake magnitudes) 6 or larger (RP_6.0) are also calculated. According to this index, the investigated seismic zones are classified into five groups and it is shown that seismic zones 3 (Quetta of Pakistan), 11 (Hindukush), 15 (northern Pamirs), and 23 (Kangra, Himachal Pradesh of India) correspond to a “very high” K index which is 6.     

Probabilistic Assessment of Earthquake Hazards in the North-East Indian Peninsula and Hindukush Regions

—The Himalayan region is one of the most seismic prone areas of the world. The North-East (NE) Indian peninsula and the Hindukush regions mark the zone of collision of the Indian and Eurasian plates. The probability of the occurrence of great earthquakes with magnitude greater than 7.0 during a specified interval of time has been estimated on the basis of four probabilistic models, namely, Weibull, Gamma, Lognormal and Exponential for the NE Indian peninsula and Hindukush regions. The model parameters have been estimated by the method of Maximum Likelihood Estimates (MLE) and the Method of Moments (MOM). The cumulative probability is estimated for a period of 40 years from 1964 and is ranging between 0.881 to 0.995 by the year 1995, using all four models for the NE Indian peninsula. The conditional probability is also estimated and it is concluded that the NE Indian peninsula would expect a great earthquake at any time in the remaining years of the present century. For the Hindukush region, the cumulative probability has already crossed its highest value, but no earthquake of magnitude greater than 7.0 has occurred after 1974 in this area. It may attribute to the occurrence of frequent shocks of moderate size, as seventeen earthquakes of magnitude greater than 6.0, including four greater than 6.4, have been reported until 1994 from this region.     

Probabilistic Seismic Hazard Assessment for Iraq Using Complete Earthquake Catalogue Files

Probabilistic seismic hazard analysis (PSHA) has been carried out for Iraq. The earthquake catalogue used in the present study covers an area between latitude 29°–38.5° N and longitude 39°–50° E containing more than a thousand events for the period 1905–2000. The entire Iraq region has been divided into thirteen seismogenic sources based on their seismic characteristics, geological setting and tectonic framework. The completeness of the seismicity catalogue has been checked using the method proposed by Stepp (1972). The analysis of completeness shows that the earthquake catalogue is not complete below M_s=4.8 for all of Iraq and seismic source zones S1, S4, S5, and S8, while it varies for the other seismic zones. A statistical treatment of completeness of the data file was carried out in each of the magnitude classes. The Frequency Magnitude Distributions (FMD) for the study area including all seismic source zones were established and the minimum magnitude of complete reporting (M_c) were then estimated. For the entire Iraq the M_c was estimated to be about M_s=4.0 while S11 shows the lowest M_c to be about M_s=3.5 and the highest M_c of about M_s=4.2 was observed for S4. The earthquake activity parameters (activity rate , b value, maximum regional magnitude m_max) as well as the mean return period (R) with a certain lower magnitude m_min m along with their probability of occurrence have been determined for all thirteen seismic source zones of Iraq. The maximum regional magnitude m_max was estimated as 7.87 ± 0.86 for entire Iraq. The return period for magnitude 6.0 is largest for source zone S3 which is estimated to be 705 years while the smallest value is estimated as 9.9 years for all of Iraq.The large variation of the b parameter and the hazard level from zone to zone reflects crustal heterogeneity and the high seismotectonic complexity. The seismic hazard near the source boundaries is directly and strongly affected by the change in the delineation of these boundaries. The forces, through which the geological structure along the plate boundary in Eastern and Northeastern Iraq are evolved, are still active causing stress-strain accumulation, deformation and in turn producing higher probabilities of earthquake activity. Thus, relatively large destructive earthquakes are expected in this region. The study is intended to serve as a reference for more advanced approaches and to pave the path for the probabilistic assessment of seismic hazard in this region.     

A Probabilistic Tsunami Hazard Assessment for Western Australia

The occurrence of the Indian Ocean Tsunami on 26 December, 2004 has raised concern about the difficulty in determining appropriate tsunami mitigation measures in Australia, due to the lack of information on the tsunami threat. A first step in the development of such measures is a tsunami hazard assessment, which gives an indication of which areas of coastline are most likely to experience tsunamis, and how likely such events are. Here we present the results of a probabilistic tsunami hazard assessment for Western Australia (WA). Compared to other parts of Australia, the WA coastline experiences a relatively high frequency of tsunami occurrence. This hazard is due to earthquakes along the Sunda Arc, south of Indonesia. Our work shows that large earthquakes offshore of Java and Sumba are likely to be a greater threat to WA than those offshore of Sumatra or elsewhere in Indonesia. A magnitude 9 earthquake offshore of the Indonesian islands of Java or Sumba has the potential to significantly impact a large part of the West Australian coastline. The level of hazard varies along the coast, but is highest along the coast from Carnarvon to Dampier. Tsunamis generated by other sources (e.g., large intra-plate events, volcanoes, landslides and asteroids) were not considered in this study.     

Probabilistic Appraisal of Earthquake Hazard Parameters Deduced from a Bayesian Approach in the Northwest Frontier of the Himalayas

A straightforward Bayesian statistic is applied in five broad seismogenic source zones of the northwest frontier of the Himalayas to estimate the earthquake hazard parameters (maximum regional magnitude M _max, β value of G–R relationship and seismic activity rate or intensity λ). For this purpose, a reliable earthquake catalogue which is homogeneous for M _W ≥ 5.0 and complete during the period 1900 to 2010 is compiled. The Hindukush–Pamir Himalaya zone has been further divided into two seismic zones of shallow (h ≤ 70 km) and intermediate depth (h > 70 km) according to the variation of seismicity with depth in the subduction zone. The estimated earthquake hazard parameters by Bayesian approach are more stable and reliable with low standard deviations than other approaches, but the technique is more time consuming. In this study, quantiles of functions of distributions of true and apparent magnitudes for future time intervals of 5, 10, 20, 50 and 100 years are calculated with confidence limits for probability levels of 50, 70 and 90 % in all seismogenic source zones. The zones of estimated M _max greater than 8.0 are related to the Sulaiman–Kirthar ranges, Hindukush–Pamir Himalaya and Himalayan Frontal Thrusts belt; suggesting more seismically hazardous regions in the examined area. The lowest value of M _max (6.44) has been calculated in Northern-Pakistan and Hazara syntaxis zone which have estimated lowest activity rate 0.0023 events/day as compared to other zones. The Himalayan Frontal Thrusts belt exhibits higher earthquake magnitude (8.01) in next 100-years with 90 % probability level as compared to other zones, which reveals that this zone is more vulnerable to occurrence of a great earthquake. The obtained results in this study are directly useful for the probabilistic seismic hazard assessment in the examined region of Himalaya.     

Probabilistic Assessment of Earthquake Hazards in the Indian Subcontinent

—The Indian subcontinent is one of the most seismic prone areas of the world. The Himalayan mountains in the north, mid-oceanic ridges in the south and earthquake belts surrounding the Indian plate all show that the subcontinent has undergone extensive geological and tectonic processes in the past. The probability of the occurrence of earthquakes with magnitude 6<Mb<7 during a specified interval of time has been estimated on the basis of four probabilistic models namely Lognormal, Weibull, Gamma and Exponential distribution for the Indian subcontinent. The seismicity map has been prepared using the earthquake catalogue from the period 1963–1994, and six different zones have been identified on the basis of clustering of events. The model parameters have been estimated by the method of maximum likelihood estimates (MLE) and method of moments (MOM). A computer program package has been developed for all four models, which represents the distributions of time intervals fairly well. The logarithmic of likelihood (ln L) is estimated for testing the models and different models have been found to be plausible. The probability of different magnitude thresholds has been evaluated using the Gutenberg–Richter formula Log N = a - bM for magnitude distribution. The constants a and b have been computed for each region and found to be varying between 5.46–8.53 and 0.87–1.34, respectively.     

Earthquake Hazard Parameters Estimated in Crete Island and the Adjacent Area

—?Earthquake hazard parameters are estimated by the application of the maximum likelihood method. The technique is based on a procedure which utilizes data of different quality, e.g., those in which the uncertainty in the assessment of the magnitudes is great and those in which the magnitudes are computed with great precision. In other words the data were extracted from both historical (incomplete) and recorded (complete) files. The historical part of the catalogue contains only the strongest events, whereas the complete part can be divided into several sub-catalogues; each one assumed to be complete above a specified magnitude threshold. Uncertainty in the determination of magnitudes has also been taken into account. The method allows us to estimate the earthquake hazard parameters which are the maximum regional magnitude, M _max, the activity rate, λ, of the seismic events and the well known value β (b=β?log?e), which is the slope of the magnitude-frequency relationship. All these parameters are of physical significance. The mean return periods, RP, of earthquakes with a certain lower magnitude M?≥?m are also determined. The method is applied in the Island of Crete and the adjacent area, where catastrophic earthquakes are known from the historical era. The earthquake hazard of the whole area is divided in a cellular manner which allow the analysis of the localized hazard parameters and the representation of their regional variation. The seismic hazard analysis, which is expressed by: (a) The annual probability of exceedance of a specified value of magnitude and (b) the return periods (in years) that are expected for given magnitudes, for shallow events is finally performed for shallow events. This hazard analysis is useful for both theoretical and practical reasons and provides a tool for earthquake resistant design in both areas of low and high seismicity.     

A Probabilistic Tsunami Hazard Study of the Auckland Region, Part II: Inundation Modelling and Hazard Assessment

Regional source tsunamis pose a potentially devastating hazard to communities and infrastructure on the New Zealand coast. But major events are very uncommon. This dichotomy of infrequent but potentially devastating hazards makes realistic assessment of the risk challenging. Here, we describe a method to determine a probabilistic assessment of the tsunami hazard by regional source tsunamis with an “Average Recurrence Interval” of 2,500-years. The method is applied to the east Auckland region of New Zealand. From an assessment of potential regional tsunamigenic events over 100,000 years, the inundation of the Auckland region from the worst 100 events is modelled using a hydrodynamic model and probabilistic inundation depths on a 2,500-year time scale were determined. Tidal effects on the potential inundation were included by coupling the predicted wave heights with the probability density function of tidal heights at the inundation site. Results show that the more exposed northern section of the east coast and outer islands in the Hauraki Gulf face the greatest hazard from regional tsunamis in the Auckland region. Incorporating tidal effects into predictions of inundation reduced the predicted hazard compared to modelling all the tsunamis arriving at high tide giving a more accurate hazard assessment on the specified time scale. This study presents the first probabilistic analysis of dynamic modelling of tsunami inundation for the New Zealand coast and as such provides the most comprehensive assessment of tsunami inundation of the Auckland region from regional source tsunamis available to date.     

地震损失概率预测研究中的地震危险性分析方法探讨

本文基于假定未来某段时间内研究地区或场地上地震烈度影响事件遵从贝努利(Bernoulli)随机独立试验过程，利用历史地震烈度-额度统计关系，提出了一种适用于地震损失概率预测的危险性分析方法。     

Earthquake Hazard Assessment for Dehradun, Uttarakhand, India, Including a Characteristic Earthquake Recurrence Model for the Himalaya Frontal Fault (HFF)

The seismic hazard has been computed for the city of Dehradun, Uttarakhand, India. The city lies in the Himalayan foothills between two faults: the Main Boundary Thrust (MBT) and the Himalayan Frontal Fault (HFF). The contributions from these two faults have been modelled differently in a probabilistic model. While the MBT has been modelled with a Poissonian earthquake distribution, the HFF has been modelled both with a characteristic earthquake recurrence model and a Poissonian model. The hazard scenarios reveal different patterns depending on the classical approaches and the characteristic models applied, and the obtained results indicate that Dehradun may experience PGA shaking around 2.2?m/s² for 225 years return period and around 4.6?m/s² for a 2,500?years return period.     

Estimation of Earthquake Hazard Parameters Based on the Unequal-Accuracy Catalog

In this paper,the maximum likelihood estimation of earthquake hazard parameters(β,Mmax)is extended to those data including incomplete and uncertain parts.We calculated earthquake hazard parameters of the Fenwei seismic belt,and the border area between Shaanxi,Sichuan,Hubei and Henan Provinces.The result is comparatively stable and gives an objective evaluation of earthquake hazard.     

A Probabilistic Displacement-based Vulnerability Assessment Procedure for Earthquake Loss Estimation

Earthquake loss estimation studies require predictions to be made of the proportion of a building class falling within discrete damage bands from a specified earthquake demand. These predictions should be made using methods that incorporate both computational efficiency and accuracy such that studies on regional or national levels can be effectively carried out, even when the triggering of multiple earthquake scenarios, as opposed to the use of probabilistic hazard maps and uniform hazard spectra, is employed to realistically assess seismic demand and its consequences on the built environment. Earthquake actions should be represented by a parameter that shows good correlation to damage and that accounts for the relationship between the frequency content of the ground motion and the fundamental period of the building; hence recent proposals to use displacement response spectra. A rational method is proposed herein that defines the capacity of a building class by relating its deformation potential to its fundamental period of vibration at different limit states and comparing this with a displacement response spectrum. The uncertainty in the geometrical, material and limit state properties of a building class is considered and the first-order reliability method, FORM, is used to produce an approximate joint probability density function (JPDF) of displacement capacity and period. The JPDF of capacity may be used in conjunction with the lognormal cumulative distribution function of demand in the classical reliability formula to calculate the probability of failing a given limit state. Vulnerability curves may be produced which, although not directly used in the methodology, serve to illustrate the conceptual soundness of the method and make comparisons with other methods.