Assessment of seismic hazard and liquefaction potential of Gujarat based on probabilistic approaches

Gujarat is one of the fastest-growing states of India with high industrial activities coming up in major cities of the state. It is indispensable to analyse seismic hazard as the region is considered to be most seismically active in stable continental region of India. The Bhuj earthquake of 2001 has caused extensive damage in terms of causality and economic loss. In the present study, the seismic hazard of Gujarat evaluated using a probabilistic approach with the use of logic tree framework that minimizes the uncertainties in hazard assessment. The peak horizontal acceleration (PHA) and spectral acceleration (Sa) values were evaluated for 10 and 2?% probability of exceedance in 50?years. Two important geotechnical effects of earthquakes, site amplification and liquefaction, are also evaluated, considering site characterization based on site classes. The liquefaction return period for the entire state of Gujarat is evaluated using a performance-based approach. The maps of PHA and PGA values prepared in this study are very useful for seismic hazard mitigation of the region in future.     

白龙江引水工程代古寺水库及其周围地区地震危险性评估

白龙江引水工程是我国拟建的一项重大战略工程,而代古寺水库是该工程的水源枢纽。代古寺水库及其周围地区(本文研究区)活动断层发育、大地震频发,故亟需开展可靠的地震危险性评估,为该研究区内的工程建设和运营保驾护航。由于传统评估方法物理依据不足,难以正确评估研究区的地震危险性,故本文采用了基于地震物理预测的地震危险性评估新方法。研究结果表明,该研究区位于海原地震区,未来100年内该研究区的地震危险性主要源于海原地震区的下一次M_S8.5标志性地震。根据断层地震活动、发震潜力与展布特征,我们预判了该标志性地震的可能发震断层和震中位置;应用地震烈度衰减关系,考虑不同震中位置,分别计算了其产生的地震烈度。为确保“百年大计”的白龙江引水工程代古寺水库水资源枢纽安全,我们建议该研究区的抗震设防烈度不宜低于8度。     

Developing a seismic source model for the Arabian Plate

A seismic source model is developed for the entire Arabian Plate, which has been affected by a number of earthquakes in the past and in recent times. Delineation and characterization of the sources responsible for these seismic activities are crucial inputs for any seismic hazard study. Available earthquake data and installation of local seismic networks in most of the Arabian Plate countries made it feasible to delineate the seismic sources that have a hazardous potential on the region. Boundaries of the seismic zones are essentially identified based upon the seismicity, available data on active faults and their potential to generate effective earthquakes, prevailing focal mechanism, available geophysical maps, and the volcanic activity in the Arabian Shield. Variations in the characteristics given by the above datasets provide the bases for delineating individual seismic zones. The present model consists of 57 seismic zones extending along the Makran Subduction Zone, Zagros Fold-Thrust Belt, Eastern Anatolian Fault, Aqaba-Dead Sea Fault, Red Sea, Gulf of Aden, Owen Fracture Zone, Arabian Intraplate, and a background seismic zone, which models the floating seismicity that is unrelated to any of the distinctly identified seismic zones. The features of the newly developed model make the seismic hazard results likely be more realistic.     

Seismicity of Gujarat

Paper describes tectonics, earthquake monitoring, past and present seismicity, catalogue of earthquakes and estimated return periods of large earthquakes in Gujarat state, western India. The Gujarat region has three failed Mesozoic rifts of Kachchh, Cambay, and Narmada, with several active faults. Kachchh district of Gujarat is the only region outside Himalaya-Andaman belt that has high seismic hazard of magnitude 8 corresponding to zone V in the seismic zoning map of India. The other parts of Gujarat have seismic hazard of magnitude 6 or less. Kachchh region is considered seismically one of the most active intraplate regions of the World. It is known to have low seismicity but high hazard in view of occurrence of fewer smaller earthquakes of M????6 in a region having three devastating earthquakes that occurred during 1819 (M _w7.8), 1956 (M _w6.0) and 2001 (M _w7.7). The second in order of seismic status is Narmada rift zone that experienced a severely damaging 1970 Bharuch earthquake of M5.4 at its western end and M????6 earthquakes further east in 1927 (Son earthquake), 1938 (Satpura earthquake) and 1997 (Jabalpur earthquake). The Saurashtra Peninsula south of Kachchh has experienced seismicity of magnitude less than 6.     

Seismic hazard map of Coimbatore using subsurface fault rupture

This study presents the future seismic hazard map of Coimbatore city, India, by considering rupture phenomenon. Seismotectonic map for Coimbatore has been generated using past earthquakes and seismic sources within 300 km radius around the city. The region experienced a largest earthquake of moment magnitude 6.3 in 1900. Available earthquakes are divided into two categories: one includes events having moment magnitude of 5.0 and above, i.e., damaging earthquakes in the region and the other includes the remaining, i.e., minor earthquakes. Subsurface rupture character of the region has been established by considering the damaging earthquakes and total length of seismic source. Magnitudes of each source are estimated by assuming the subsurface rupture length in terms of percentage of total length of sources and matched with reported earthquake. Estimated magnitudes match well with the reported earthquakes for a RLD of 5.2% of the total length of source. Zone of influence circles is also marked in the seismotectonic map by considering subsurface rupture length of fault associated with these earthquakes. As earthquakes relive strain energy that builds up on faults, it is assumed that all the earthquakes close to damaging earthquake have released the entire strain energy and it would take some time for the rebuilding of strain energy to cause a similar earthquake in the same location/fault. Area free from influence circles has potential for future earthquake, if there is seismogenic source and minor earthquake in the last 20 years. Based on this rupture phenomenon, eight probable locations have been identified and these locations might have the potential for the future earthquakes. Characteristic earthquake moment magnitude (M _w ) of 6.4 is estimated for the seismic study area considering seismic sources close to probable zones and 15% increased regional rupture character. The city is divided into several grid points at spacing of 0.01° and the peak ground acceleration (PGA) due to each probable earthquake is calculated at every grid point in city by using the regional attenuation model. The maximum of all these eight PGAs is taken for each grid point and the final PGA map is arrived. This map is compared to the PGA map developed based on the conventional deterministic seismic hazard analysis (DSHA) approach. The probable future rupture earthquakes gave less PGA than that of DSHA approach. The occurrence of any earthquake may be expected in near future in these eight zones, as these eight places have been experiencing minor earthquakes and are located in well-defined seismogenic sources.     

An overview on the seismic zonation and microzonation studies in India

The present study presents a review on the progressive development of the seismic zonation map of India both from official agencies and also from independent individual studies. The zonation map have been modified and updated regularly with the occurrence of major destructive earthquakes over the years in the Indian subcontinent with the addition of new data. This study discusses the criteria chosen for the progressive zonation and the major earthquakes that were responsible for retrospection of the earlier published maps. The seismic zonation maps of India have also been prepared by various independent workers by adopting different approaches to achieve the purpose of the zonation. Despite the endeavors from various sources to provide a solution for the problem of earthquake hazards in India, there were many limitations on the zonation map as it gives the picture at a regional scale mostly on the bedrock level without addressing the local site conditions. But nevertheless, the seismic zonation map gives basic guidelines for any region to know the hazard scenario and if any city or urban population is under threat from seismic point of view, further site specific seismic microzonation may be carried out. In the International scenario, the Global Seismic Hazard Assessment Program (GSHAP) in 1999 prepared a hazard map for world in terms of peak ground acceleration (PGA) with a 10% probability of exceedance in 50 years, but it turned out to be an underestimation of the hazard parameter when compared with the observed PGA. To tackle the problem of seismic hazards, there was a need to have a detail study on the local site conditions in terms of its geological, geophysical and geotechnical properties. With the advent of better instrumentation and knowledge on the mechanics of earthquakes, it was possible to identify zones of hazards at a local level and this gives rise to the study of seismic microzonation. Seismic microzonation work has been carried out in India in some of the strategic important mega cities and industrial build up that has the potential of being damaged from future earthquakes, as has been shown in the past. Though the microzonation map is not the final output map, as it can still be updated at later stage with more input data, it does provide a more realistic picture on the site specific seismic hazard.     

青藏高原巴塘断裂带地震滑坡危险性预测研究

全新世以来青藏高原东部巴塘断裂带活动强烈,地形地貌和地质构造复杂,历史地震频发,并诱发大量滑坡灾害。基于巴塘断裂带地震滑坡长期防控的需要,在分析区域地质灾害成灾背景和发育分布特征的基础上,采用Newmark模型完成了巴塘断裂带50年超越概率10%的潜在地震滑坡危险性预测评价,并完成地震滑坡危险性区划。结果表明:巴塘断裂带及其临近的金沙江断裂带区域、金沙江及其支流沿岸具有较高的潜在地震滑坡危险性,地震滑坡危险区具有沿断裂带和大江大河等峡谷区分布的总体趋势,受活动断裂和地形地貌影响显著;距离断层越近、坡度越大的斜坡,地震滑坡危险性越高;规划建设中的川藏铁路经巴塘县德达乡、白玉县沙马乡,向西北延伸,跨越金沙江,可以穿越较少的地震滑坡危险区,金沙江水电工程规划建设需加强潜在地震滑坡危害研判及防控。巴塘断裂带潜在地震滑坡危险性评价结果可为区域城镇开发和重大工程规划建设的地震滑坡长期防控提供科学参考。     

Occurrence and impact of characteristic earthquakes in northern Algeria

Earthquakes constitute the natural hazard that is one of the main natural threats to the northern part of Algeria because of its geographical setting at the Eurasia–Africa plate boundary. Several active multi-segment reverse faults have been identified near urban areas that may rupture during characteristic earthquakes and produce earthquakes of magnitudes ≥7.0. Characteristic earthquakes are extreme seismic events characterized by long return periods, which can have great societal impact. Earthquakes in northern Algeria are destructive for two main reasons: firstly, the shallow character of the faults and secondly, the vulnerability of the building stock built essentially prior to the implementation of seismic design codes that take into account the level of the seismic hazard. That is why even moderate earthquakes are disastrous in this area.     

Spatial-temporal variability of seismic hazard in peninsular India

This paper examines the variability of seismic activity observed in the case of different geological zones of peninsular India (10°N–26°N; 68°E–90°E) based on earthquake catalog between the period 1842 and 2002 and estimates earthquake hazard for the region. With compilation of earthquake catalog in terms of moment magnitude and establishing broad completeness criteria, we derive the seismicity parameters for each geologic zone of peninsular India using maximum likelihood procedure. The estimated parameters provide the basis for understanding the historical seismicity associated with different geological zones of peninsular India and also provide important inputs for future seismic hazard estimation studies in the region. Based on present investigation, it is clear that earthquake recurrence activity in various geologic zones of peninsular India is distinct and varies considerably between its cratonic and rifting zones. The study identifies the likely hazards due to the possibility of moderate to large earthquakes in peninsular India and also presents the influence of spatial rate variation in the seismic activity of this region. This paper presents the influence of source zone characterization and recurrence rate variation pattern on the maximum earthquake magnitude estimation. The results presented in the paper provide a useful basis for probabilistic seismic hazard studies and microzonation studies in peninsular India.     

Earthquake hazard of dams along the Mekong mainstream

In this study, the earthquake hazard was evaluated for all of 19 of the proposed or built dams along the Mekong River. All values representing a potential indication of hazardous earthquakes, such as the closest earthquake and seismogenic faults and including the seismic parameters required for a seismic safety evaluation, were clarified. The results will be useful in reviewing the safety of existing dams and for the design of suitable earthquake resistant specifications for any currently or future planned dam construction in this area. Seismotectonically, 14 of the 19 proposed Mekong River dams are located within an earthquake source zone. Most of faults are potentially still active, according to both seismicity and paleoseismological evidence. In addition, the maximum credible earthquakes were estimated to be in the range of 7–8 Mw for the closest fault zone of each dam. Previous isoseismal maps indicated a risk of shaking intensities of around scale III–IV (Modified Mercalli Intensity Scale) for the dams. According to the preliminary ranging of the International Commission on Large Dams, 9 of these 19 dams are classified as in an extreme hazard class and so need careful observation and monitoring of hazardous earthquakes. An effective mitigation plan should also be prepared for each operating dam.     

中国活动构造研究的进展与展望

本文试图在世界活动构造研究背景下概括中国活动构造研究的进展，并对最近20年来中国的活动构造研究，包括几何学、滑动速率、古地震和大震重复间隔、分段理论和应用、活动褶皱、深浅构造耦合关系、城市和海域活动构造探测、块体划分和块体运动、地震危险性的未来断裂滑动量评价及中国现代构造活动动力学等进行简要评述，对今后在这些方面的发展进行概略的讨论。     

Study of active faults: Theoretical and applied implications

Theoretical and applied implications of the study of active faults and their identification and parametrization are discussed. The term active fault is defined as a fault with displacements that occurred in the late Pleistocene and Holocene and are expected to resume in the future. The displaced young landforms, late Quaternary sediments, and man-made constructions are the main reference marks to recognize active faults and estimate their kinematics and intensity. Since the structural pattern and parameters of all active faults are referred to the same, geologically short time interval, they are important for the study of recent geodynamics and young tectogenesis on the global, regional, and local scales. The opportunities that are opened for such investigations are illustrated by verification of the real existence of the Okhotsk and Bering minor plates. With allowance for active faults, it is possible to make tectonic and geodynamic reconstructions of the past events more plausible. Natural hazards, primarily, earthquakes, are related to active faults. The geological and geomorphic assessment of the seismic potential of active fault zones is discussed with emphasis on the maximum possible magnitude of earthquakes (M _max) estimated from the segmentation of the active zones, the length of particular segments (L), the empirical relationships M _max/L for strong recent earthquakes, as well as the resumption periods of strong earthquakes and measurement of particular seismogenic offsets by trenching and other techniques, including archeoseismological methods. A question is posed about possible perennial variations in the stress-and-strain state of active zones, which are expressed on the scale of large seismoactive regions in oscillations of released seismic energy and should be taken into account by assessment of the seismic hazard.     

Potential source zones for Himalayan earthquakes: constraints from spatial�Ctemporal clusters

The Himalayan fold-thrust belt has been visited by many disastrous earthquakes (magnitude > 6) time and again. This active collisional orogen bordering Indian subcontinent in the north remains a potential seismic threat of similar magnitude in the adjoining countries like India, Pakistan, Nepal, Bhutan and China. Though earthquake forecasting is riddled with all conjectures and still not a proven presumption, identifying likely source zones of such disastrous earthquakes would be an important contribution to seismic hazard assessment. In this study, we have worked out spatio-temporal clustering of earthquakes (Mb ?? 4.5; 1964?C2006) in the Himalayas. ??Point density?? spatial statistics has helped in detecting 22 spatial seismicity clusters. Earthquake catalog is then treated with a moving time-distance window technique (inter-event time 35 days and distance 100 ± 20 km) to bring out temporal clusters by recognizing several foreshock-main shock-aftershock (FMA) sequences. A total of 53 such temporal sequences identified in the process are confined within the 22 spatial clusters. Though each of these spatio-temporal clusters deserves in-depth analysis, we short-listed only eight such clusters that are dissected by active tectonic discontinuities like MBT/MCT for detail study. Spatio-temporal clusters have been used to constrain the potential source zones. These eight well-defined spatio-temporal clusters demonstrate recurrent moderate to large earthquakes. We assumed that the length of these clusters are indicating the possible maximum rupture lengths and thus empirically estimated the maximum possible magnitudes of eight clusters that can be generated from them (from west to east) as 8.0, 8.3, 8.2, 8.3, 8.2, 8.4, 8.0 and 7.7. Based on comparative study of the eight cluster zones contemplating with their temporal recurrences, historical seismic records, presence of intersecting faults and estimated magnitudes, we have guessed the possibility that Kangra, East Nepal, Garhwal and Kumaun?CWest Nepal clusters, in decreasing order of earthquake threat, are potential source zones for large earthquakes (??7.7 M) in future.     

Influence of slow active faults in probabilistic seismic hazard assessment: the northwestern margin of the València trough

Areas of low strain rate are typically characterized by low to moderate seismicity. The earthquake catalogs for these regions do not usually include large earthquakes because of their long recurrence periods. In cases where the recurrence period of large earthquakes is much longer than the catalog time span, probabilistic seismic hazard is underestimated. The information provided by geological and paleo-seismological studies can potentially improve seismic hazard estimation through renewal models, which assume characteristic earthquakes. In this work, we compare the differences produced when active faults in the northwestern margin of the València trough are introduced in hazard analysis. The differences between the models demonstrate that the introduction of faults in zones characterized by low seismic activity can give rise to significant changes in the hazard values and location. The earthquake and fault seismic parameters (recurrence interval, segmentation or fault length that controls the maximum magnitude earthquake and time elapsed since the last event or T_e) were studied to ascertain their effect on the final hazard results. The most critical parameter is the recurrence interval, where shorter recurrences produce higher hazard values. The next most important parameter is the fault segmentation. Higher hazard values are obtained when the fault has segments capable of producing big earthquakes. Finally, the least critical parameter is the time elapsed since the last event (T_e), when longer T_e produces higher hazard values.     

Probabilistic seismic hazard analysis for Sumatra, Indonesia and across the Southern Malaysian Peninsula

The ground motion hazard for Sumatra and the Malaysian peninsula is calculated in a probabilistic framework, using procedures developed for the US National Seismic Hazard Maps. We constructed regional earthquake source models and used standard published and modified attenuation equations to calculate peak ground acceleration at 2% and 10% probability of exceedance in 50 years for rock site conditions. We developed or modified earthquake catalogs and declustered these catalogs to include only independent earthquakes. The resulting catalogs were used to define four source zones that characterize earthquakes in four tectonic environments: subduction zone interface earthquakes, subduction zone deep intraslab earthquakes, strike-slip transform earthquakes, and intraplate earthquakes. The recurrence rates and sizes of historical earthquakes on known faults and across zones were also determined from this modified catalog. In addition to the source zones, our seismic source model considers two major faults that are known historically to generate large earthquakes: the Sumatran subduction zone and the Sumatran transform fault. Several published studies were used to describe earthquakes along these faults during historical and pre-historical time, as well as to identify segmentation models of faults. Peak horizontal ground accelerations were calculated using ground motion prediction relations that were developed from seismic data obtained from the crustal interplate environment, crustal intraplate environment, along the subduction zone interface, and from deep intraslab earthquakes. Most of these relations, however, have not been developed for large distances that are needed for calculating the hazard across the Malaysian peninsula, and none were developed for earthquake ground motions generated in an interplate tectonic environment that are propagated into an intraplate tectonic environment. For the interplate and intraplate crustal earthquakes, we have applied ground-motion prediction relations that are consistent with California (interplate) and India (intraplate) strong motion data that we collected for distances beyond 200 km. For the subduction zone equations, we recognized that the published relationships at large distances were not consistent with global earthquake data that we collected and modified the relations to be compatible with the global subduction zone ground motions. In this analysis, we have used alternative source and attenuation models and weighted them to account for our uncertainty in which model is most appropriate for Sumatra or for the Malaysian peninsula. The resulting peak horizontal ground accelerations for 2% probability of exceedance in 50 years range from over 100% g to about 10% g across Sumatra and generally less than 20% g across most of the Malaysian peninsula. The ground motions at 10% probability of exceedance in 50 years are typically about 60% of the ground motions derived for a hazard level at 2% probability of exceedance in 50 years. The largest contributors to hazard are from the Sumatran faults.     

Potential of future seismogenesis in Hebei Province (NE China) due to stress interactions between strong earthquakes

Northeast China, a densely populated area, is affected by intense seismic activity, which includes large events that caused extensive disaster and tremendous loss of life. For contributing to the continuous efforts for seismic hazard assessment, the earthquake potential from the active faults near the cities of Zhangjiakou and Langfang in Hebei Province is examined. We estimate the effect of the coseismic stress changes of strong (M ⩾ 5.0) earthquakes on the major regional active faults, and mapped Coulomb stress change onto these target faults. More importantly our calculations reveal that positive stress changes caused by the largest events of the 1976 Tangshan sequence make the Xiadian and part of Daxing fault, thus considered the most likely sites of the next strong earthquake in the study area. The accumulated static stress changes that reached a value of up to 0.4 bar onto these faults, were subsequently incorporated in earthquake probability estimates for the next 30 years.     

Earthquake prediction activities and Damavand earthquake precursor test site in Iran

Iran has long been known as one of the most seismically active areas of the world, and it frequently suffers destructive and catastrophic earthquakes that cause heavy loss of human life and widespread damage. The Alborz region in the northern part of Iran is an active EW trending mountain belt of 100 km wide and 600 km long. The Alborz range is bounded by the Talesh Mountains to the west and the Kopet Dagh Mountains to the east and consists of several sedimentary and volcanic layers of Cambrian to Eocene ages that were deformed during the late Cenozoic collision. Several active faults affect the central Alborz. The main active faults are the North Tehran and Mosha faults. The Mosha fault is one of the major active faults in the central Alborz as shown by its strong historical seismicity and its clear morphological signature. Situated in the vicinity of Tehran city, this 150-km-long N100° E trending fault represents an important potential seismic source. For earthquake monitoring and possible future prediction/precursory purposes, a test site has been established in the Alborz mountain region. The proximity to the capital of Iran with its high population density, low frequency but high magnitude earthquake occurrence, and active faults with their historical earthquake events have been considered as the main criteria for this selection. In addition, within the test site, there are hot springs and deep water wells that can be used for physico-chemical and radon gas analysis for earthquake precursory studies. The present activities include magnetic measurements; application of methodology for identification of seismogenic nodes for earthquakes of M ≥ 6.0 in the Alborz region developed by International Institute of Earthquake Prediction Theory and Mathematical Geophysics, IIEPT RAS, Russian Academy of Science, Moscow (IIEPT&MG RAS); a feasibility study using a dense seismic network for identification of future locations of seismic monitoring stations and application of short-term prediction of medium- and large-size earthquakes is based on Markov and extended self-similarity analysis of seismic data. The establishment of the test site is ongoing, and the methodology has been selected based on the IASPEI evaluation report on the most important precursors with installation of (i) a local dense seismic network consisting of 25 short-period seismometers, (ii) a GPS network consisting of eight instruments with 70 stations, (iii) magnetic network with four instruments, and (iv) radon gas and a physico-chemical study on the springs and deep water wells.     

Buckle-controlled seismogenic faulting in peninsular India

As intraplate earthquakes are often not associated with major known faults their location as well as their timing is unpredictable. In peninsular India the larger (M5.0) events occur mainly on reverse faults in a series of belts 400–800 km apart which are aligned roughly normal to the azimuth of convergence between the Indian and Eurasian plates. The location of the belts is controlled largely by the buckling wavelength of the lithosphere, and the seismogenic faults do not generate folding and sometimes result from it. There is consequently no need to postulate the creation of regularly spaced normal faults in an antecedent extensional phase, and the deformation is consistent with a plate-driving force such as gravity glide which is unlikely to reverse its polarity and which creates structures that are influenced by plate geometry at the leading edge. The thesis is potentially of value to seismic hazard mitigation as it identifies the zones that are most at risk.     

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

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