First Order Seismic Microzonation of Delhi,India Using Geographic Information System (GIS)

A first order seismic microzonation map of Delhi is prepared using five thematic layers viz., Peak Ground Acceleration (PGA) contour, different soil types at 6 m depth, geology, groundwater fluctuation and bedrock depth, integrated on GIS platform. The integration is performed following a pair-wise comparison of Analytical Hierarchy Process (AHP), wherein each thematic map is assigned weight in the 5-1 scale: depending on its contribution towards the seismic hazard. Following the AHP, the weightage assigned to each theme are: PGA (0.333), soil (0.266), geology (0.20), groundwater (0.133) and bedrock depth (0.066). The thematic vector layers are overlaid and integrated using GIS. On the microzonation theme, the Delhi region has been classified into four broad zones of vulnerability to the seismic hazard. They are very high (> 52%), high (38–52%), moderate (23–38%) and less ( < 23%) zones of seismic hazard. The “very high” seismic hazard zone is observed where the maximum PGA varies from 140 to 210 gal for a finite source model of M_w 8.5 in the central seismic gap. A site amplification study from local and regional earthquakes for Delhi region using Delhi Telemetry Network data shows a steeper site response gradient in the eastern side of the Yamuna fluvial deposits at 1.5 Hz. The ‘high’ seismic hazard zone occupies most of the study area where the PGA value ranges from 90 to 140 gal. The ‘moderate’ seismic hazard zone occurs on either side of the Delhi ridge with PGA value varying from 60 to 90 gal. The ‘less’ seismic hazard zone occurs in small patches distributed along the study area with the PGA value less than 60 gal. Site response studies, PGA distribution and destruction pattern of the Chamoli earthquake greatly corroborate the seismic hazard zones estimated through microzonation on GIS platform and also establishes the methodology incorporated in this study.     

地震危险性、地震危害性和地震易损性

做好城市防震减灾工作 ,关键是对未来可能遭遇地震灾害定量化预测。本文论述地震灾害定量化的三要素 :地震危险性 ,地震危害性和地震易损性 ,它们概念完全不同 ,且很容易混淆 ,但又存在因与果的关系。本文还简要介绍厦门市地震科技工作者开展的“闽南地区综合防震减灾示范工程”。     

Seismic hazard of Egypt

Earthquake hazard parameters such as maximum expected magnitude,M _max, annual activity rate,, andb value of the Gutenberg-Richter relation have been evaluated for two regions of Egypt. The applied maximum likelihood method permits the combination of both historical and instrumental data. The catalogue used covers earthquakes with magnitude 3 from the time interval 320–1987. The uncertainties in magnitude estimates and threshold of completeness were taken into account. The hazard parameter determination is performed for two study areas. The first area, Gulf of Suez, has higher seismicity level than the second, all other active zones in Egypt.b-values of 1.2 ± 0.1 and 1.0 ± 0.1 are obtained for the two areas, respectively. The number of annually expected earthquakes with magnitude 3 is much larger in the Gulf of Suez, 39 ± 2 than in the other areas, 6.1 ± 0.5. The maximum expected magnitude is calculated to be 6.5 ± 0.4 for a time span of 209 years for the Gulf of Suez and 6.1 ± 0.3 for a time span of 1667 years for the remaining active areas in Egypt. Respective periods of 10 and 20 years were reported for earthquakes of magnitude 5.0 for the two subareas.     

Seismic hazard studies for Gaziantep city in South Anatolia of Turkey

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

Site-specific earthquake response study for hazard assessment in Kolkata city,India

The assessment of local site effects on seismic ground motions is of great importance in earthquake engineering practice. Several destructive earthquakes in the past have demonstrated that the amplification of ground motion and associated damage to structures due to local site conditions is a significant consideration in earthquake hazard analysis. A recent paper published in this journal highlights the hazard posed by earthquakes in the megacity of Kolkata in India due to its seismic and geological settings. The seismic hazard assessment study speculates that the deep alluvial deposit in the city may increase the seismic hazard probably due to the amplification of the seismic energies. This paper focuses on the seismic response studies of the various soil strata (i.e. for local subsurface conditions) obtained from various construction sites in the city for predicted earthquake. It is very well recognized that site response studies (a part of seismic microhazard zonation for urban areas) are the first step towards performance-based foundation design or seismic risk analysis and mitigation strategy. One of the problems for carrying out site-specific study in Kolkata is the lack of recorded strong motion data in the city. Hence, this paper outlines a methodology to carry out site-specific study, where no strong motion data or seismic data are available. The methodology uses wavelet-based spectrum compatibility approach to generate synthetic earthquake motions and equivalent linear method for seismic site response analysis. The Mega City of Kolkata has been considered to explain the methodology. Seismic hazard zonation map by the Bureau of Indian Standards classifies the City of Kolkata as moderate seismic zone (Zone III) with a zone factor 0.16. On the other hand, GSHAP(Global Seismic Hazard Assessment Program) map which is based on 10% probability of exceedance in 50 years specifies a maximum peak ground acceleration (PGA) of 1.6 m/s² (0.163 g) for this region. In the present study, the seismic response has been carried out based on GSHAP. The results of the analysis indicate the amplification of ground motion in the range of 4.46–4.82 with the fundamental period ranging from 0.81 to 1.17 s. Furthermore, the maximum spectral accelerations vary in the range of 0.78–0.95 g.     

Seismic hazard assessment in Jordan and its vicinity

The earthquake hazard in Jordan and its vicinity is assessed on the basis of probabilistic methods. For this purpose, an updated earthquake catalog is compiled which covers the period between AD 1–1989. The earthquakes lie between latitudes 27.0°-35.5° N and longitudes 32.0°-39.0° E. Thirteen seismic zones are defined on a regional seismic and tectonic map presented for the area. Point-source and line-source models are used. The seismic hazard parameters, namely, theb-parameter (of the Gutenberg-Richter relation),m ₁ (the upper bound magnitude), and ₄ (the annual rate of occurrence of earthquakes with local magnitudeM _L 4.0) are calculated for each zone. The results of the seismic hazard assessment are displayed as iso-acceleration contours expected to be exceeded during typical economic life times of structures, i.e. 50 and 100 years. For each model, two seismic hazard maps are derived. In order to determine the importance of the South-eastern Mediterranean zone and the north part of the Red Sea zone from a seismic hazard point of view for Jordan, one seismic hazard map which corresponds to 50 years' economic life for every model, excluding the seismicity of these zones, is derived.     

Seismic hazard in the Sannio-Matese area of Italy

The assessment of seismic hazard at five selected sites in the Sannio-Matese region is based on the computer program SRAMSC. Owing to the extensive historical data base for the output parameter, the MSK intensiy is chosen. The seismicity model is made up of five narrow area seismic sources. Circular or elliptical macroseismic fields are assigned to individual sources. A generalized Kövesligethy equation is used for this purpose as the attenuation relationship. The study reveals similar and a rather high hazard at the sites at Benevento, Boiano, and Melfi, which are located in the zone of highest seismic activity. At the Pomigliano and Lucera sites, the assessed hazard is much lower.     

Seismic hazard maps of Mexico, the Caribbean, and Central and South America

The growth of megacities in seismically active regions around the world often includes the construction of seismically unsafe buildings and infrastructures due to an insufficient knowledge of existing seismic hazard and/or economic constraints. Minimization of the loss of life, property damage, and social and economic disruption due to earthquakes depends on reliable estimates of seismic hazard. We have produced a suite of seismic hazard estimates for Mexico, the Caribbean, and Central and South America. One of the preliminary maps in this suite served as the basis for the Caribbean and Central and South America portion of the Global Seismic Hazard Map (GSHM) published in 1999, which depicted peak ground acceleration (pga) with a 10% chance of exceedance in 50 years for rock sites. Herein we present maps depicting pga and 0.2 and 1.0 s spectral accelerations (SA) with 50%, 10%, and 2% chances of exceedance in 50 years for rock sites. The seismicity catalog used in the generation of these maps adds 3 more years of data to those used to calculate the GSH Map. Different attenuation functions (consistent with those used to calculate the U.S. and Canadian maps) were used as well. These nine maps are designed to assist in global risk mitigation by providing a general seismic hazard framework and serving as a resource for any national or regional agency to help focus further detailed studies required for regional/local needs. The largest seismic hazard values in Mexico, the Caribbean, and Central and South America generally occur in areas that have been, or are likely to be, the sites of the largest plate boundary earthquakes. High hazard values occur in areas where shallow-to-intermediate seismicity occurs frequently.     

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.     

Seismic hazard in Northern Algeria using spatially smoothed seismicity. Results for peak ground acceleration

Seismic hazard in terms of peak ground acceleration (PGA) has been evaluated in northern Algeria using spatially smoothed seismicity data. We present here a preliminary seismic zoning in northern Algeria as derived from the obtained results.Initially, we have compiled an earthquake catalog of the region taking data from several agencies. Afterwards, we have delimited seismic areas where the b and m_max parameters are different. Finally, by applying the methodology proposed by Frankel [Seismol. Res. Lett. 66 (1995) 8], and using four complete and Poissonian seismicity models, we are able to compute the seismic hazard maps in terms of PGA with 39.3% and 10% probability of exceedance in 50 years.A significant result of this work is the observation of mean PGA values of the order of 0.20 and 0.45 g, for return periods of 100 and 475 years, respectively, in the central area of the Tell Atlas.     

Seismic hazard in East Africa: An example of the application of incomplete and uncertain data

Earthquake hazard parameters maximum regional magnitude M _max and annual activity rate , and the b parameter of the Gutenberg-Richter relation have been evaluated for parts of East Africa. The applied maximum likelihood method permits the combination of both historical and instrumental data available (the catalog used here covers the interval 1880–1979). In addition, the uncertainty involved in magnitude determination and thresholds of completeness were taken into account. The hazard-parameter determination was performed for two study areas corresponding to segments of the eastern and western branches of the East African rift system. The results for these areas show differences that can be partly explained through the characteristics of the data.     

Distribution of the energy release, b-values and seismic hazard in Egypt

A review of the seismicity and seismic history of Egypt indicates areas of high activity concentrated along Oligocene-Miocene faults. This supports the idea of recent activation of the Oligocene-Miocene stress cycle. There are similarities in the spatial distribution of recent and historical epicenters. Destructive earthquakes in Egypt are mostly concentrated in the highly populated areas of the Nile Valley and Nile Delta. Some big earthquakes located near the plate boundary as far away as Turkey and Crete were strongly felt in Egypt. The distribution of the energy release shows a possible tectonic connection between active zones in Egypt and the complicated tectonic zones in Turkey and Crete through geologically verified fault systems. The distribution of intensity shows a strong directivity along the Nile Valley. This is due to the presence of a thick layer of loose sediments on top of the hard rock in the Nile Valley graben. The distribution of b-values indicates two different zones, comparable with stable and unstable shelf areas. Stress loads in the northern Red Sea and northern Egypt are similar. Geologically, northern Egypt is a part of the Unstable Shelf area. The probability to have an earthquake with intensity V or larger within 94 years is more than 80% in the Nile Valley and Nile Delta areas, Egypt-Mediterranean coastal area, Aswan High Dam area, Gulf of Aqaba-Levant Fault zone and in the oil fields of the Gulf of Suez. The maximum expected intensity in these areas and within the same period is V–VI for a 80% probability and VII–VIII+ for a 10% probability. Intensity VIII–IX has been reported for several earthquakes in both historical and recent time.     

Seismic vulnerability assessment of earthquake-prone mega-city Shillong,India using geophysical mapping and remote sensing

ABSTRACT

The concept of seismic vulnerability is a yard-stick of damage estimation from a probable earthquake considering physical cum social dimension and enables a basis for decision-makers to develop preparedness and mitigation strategies. We aim at vulnerability assessment of the typical urban system of capital city Shillong situated on hilly terrain. High-resolution satellite imagery of Shillong facilitates analysis of building footprints, communication network, and open ground. Different building typologies are identified taking into account the building’s structural configuration assessed through a rapid visual survey of more than 15% of total residential households. Slope map demarcates the landslide-prone area through discrete elevation modelling. A methodology incorporating several parameters e.g. building typology, slope angle, shear wave velocity characteristics, geomorphology, and the number of occupants in correlation with a physical measurement of vulnerability is presented and is applied to estimate the dimension of vulnerability. Additionally, MASW survey indicates lithology up to 30?m deep along with the existence of stiff soil and rocks at different depths whereas resonant frequency is identified to be in the range of 6–8?Hz through H/V ratio. Integrating all, it is observed that more than 60% of Shillong city falls under moderate to higher vulnerability and the rest is less vulnerable.     

Ground zoning against seismic hazard in Athens, Greece

The 1999, Ms=5.9, Athens earthquake caused serious structural damage to buildings in the western part of Athens, Greece. This paper presents the ground zoning against seismic hazard proposed shortly after the earthquake in order to aid reconstruction of the area. Existing engineering geological and geotechnical data were combined with local observations to provide a unified set of classification criteria, consistent with provisions of the Greek Seismic Code EAK. The accuracy and the possible limitations of this zoning procedure are addressed through comparison with observed damage distribution as well as results from seismic ground response analyses performed at sites with well established soil profiles. There is clear evidence that the proposed zones correspond to geological formations exhibiting grossly different seismic response with regard to the design of common engineering structures. However, the mostly qualitative nature of the guidelines for ground categorisation provided by EAK and the general lack of systematic, site-specific geotechnical data for the whole area induce uncertainties in the definition of the seismic design actions for the different zones. These objective uncertainties certainly demand increased conservatism but do not limit application of the proposed methodology for first aid, preliminary planning in the event of destructive earthquakes.     

Tmax of asphaltenes: a parameter for oil maturity assessment

A new maturity parameter determined on both oil and bitumen samples, the asphaltene T_max, is proposed and discussed. This parameter could be very useful to address the maturity of the source rock. The asphaltene T_max is measured by programmed Rock-Eval pyrolysis, using a modified temperature program. Some phases of the experimental procedure, such as the asphaltene preparation and the Rock-Eval measurement substratum choice, are crucial in order to achieve reliable data. Laboratory simulations were carried out in order to assess the possible effects of both primary and secondary migration on asphaltene T_maxin the expelled oil: the original value of the asphaltene T_max in the bitumen is not substantially modified and it is very close to that measured on kerogen. Examples of the determination of asphaltene T_max on many samples, collected from different areas and with different organic matter composition, are given. Results show that T_max values from oil asphaltenes are reasonable indicators of source rock maturity.     

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.     

Seismic activity changes progressing simultaneously with slow-slip in the Tokai area

Anomalous movements were detected simultaneously in both the seismic and the GPS observations in the Tokai area, the central part of the Japanese islands from the late 1990s to 2000. The anomalies are of great concern since the pending risk of a large megathrust earthquake in this area has been predicted for more than 20 years. The GPS data revealed that a slow-slip on the plate interface had commenced beneath Lake Hamana, the center of which is positioned around the edge of the assumed focal zone. On the other hand, the seismic data indicated that a delicate but clear quiescence appeared over a wide area that spreads into the main focal zone. Analyses of the seismicity changes in space and time confirmed that the contrast in the seismicity rate is distinct inside the focal zone. While the integrated seismicity indicated quiescence, some locations were distinguished as activated zones, possibly indicating the appearance of asperities. The rise of the seismicity rate in a quasi-stationary manner suggests an increase in the stress rate at that location. The following hypothesis is proposed based on the simultaneously detected evidences. The slow-slip progressing beneath Lake Hamana will induce a stress shift that invades the interior of the main locked zone, which will increase the contrast of the seismicity rate, possibly reflecting inhomogeneity in the locking strength. Even in this stage, the activated zones still maintain a locked state to prevent overall breakage. Investigations of the b-value changes and of tidal dependence in seismicity that reveal the stress-concentrated state also support the hypothesis. If this is the case, the observed change in seismicity would indicate the process of stress redistribution in the locking state, which represents the preparatory process toward final breakage. Tracking such seismicity changes would yield valid information for predictions of the next Tokai earthquake.     

Seismic investigation of near-surface geological structure in the Paducah, Kentucky, area: application to earthquake hazard evaluation

Seismics method were used to evaluate shallow geological conditions at 33 sites in the vicinity of Paducah, Kentucky. A combined set of P- and S-wave seismic refraction and reflection soundings were used, in addition to local borehole information, to produce structure maps of (1) a shallow (< 30 m deep) horizon believed to represent an unconformity surface at the top of the Eocene, and (2) the Paleozoic bedrock surface (< 85 to > 160 m deep). Shear-wave velocity contrasts across the shallow unconformity were generally 2-to-1 while the contrast at the top of the Paleozoic bedrock exceeds 5-to-1. These seismic boundaries have been determined to be very important in modelling and interpreting earthquake ground motion amplification in the Paducah area. The quality and accuracy of the data, and the cost effective nature of the methods, suggest that other communities in areas at risk to damage from seismic activity, with foundation conditions comparable to Paducah, might benefit from similar characterization in order (1) to identify seismically hazardous, near-surface, geological conditions, and (2) to develop geological models that could be used in computer simulations of site response.     

The Olmsted fault zone, southernmost Illinois: A key to understanding seismic hazard in the northern new Madrid seismic zone

Geological deformation in the northern New Madrid seismic zone, near Olmsted, Illinois (USA), is analyzed using integrated compressional-wave (P) and horizontally polarized-wave (SH) seismic reflection and regional and dedicated borehole information. Seismic hazards are of special concern because of strategic facilities (e.g., lock and dam sites and chemical plants on the Ohio River near its confluence with the Mississippi River) and because of alluvial soils subject to high amplification of earthquake shock. We use an integrated approach starting with lower resolution, but deeper penetration, P-wave reflection profiles to identify displacement of Paleozoic bedrock. Higher resolution, but shallower penetration, SH-wave images show deformation that has propagated upward from bedrock faults into Pleistocene loess. We have mapped an intricate zone more than 8 km wide of high-angle faults in Mississippi embayment sediments localized over Paleozoic bedrock faults that trend north to northeast, parallel to the Ohio River. These faults align with the pattern of epicenters in the New Madrid seismic zone. Normal and reverse offsets along with positive flower structures imply a component of strike-slip; the current stress regime favors right-lateral slip on northeast-trending faults. The largest fault, the Olmsted fault, underwent principal displacement near the end of the Cretaceous Period 65 to 70 million years ago. Strata of this age (dated via fossil pollen) thicken greatly on the downthrown side of the Olmsted fault into a locally subsiding basin. Small offsets of Tertiary and Quaternary strata are evident on high-resolution SH-wave seismic profiles. Our results imply recent reactivation and possible future seismic activity in a critical area of the New Madrid seismic zone. This integrated approach provides a strategy for evaluating shallow seismic hazard-related targets for engineering concerns.