Multichannel analysis of surface waves (MASW) for seismic site characterization using 2D genetic algorithm at Bahrah area,Wadi Fatima,Saudi Arabia

The present study deals with the seismic site classification of Bahrah area, Wadi Fatima, to characterize the local site conditions. The dynamic behavior of sediments was studied by the application of surface wave inversion. The multichannel analysis of surface waves (MASW) shallow geophysical technique was utilized for site classification. MASW survey was carried out at 66 sites along with 13 seismic refraction profiles at suitable localities. MASW and seismic refraction profiles were processed and compared with the available borehole data. The integration of MASW techniques with seismic refraction and borehole data progressively enhanced the subsurface visualization and reliability of the shear wave velocity estimation in the subsurface in the study area. The subsurface shear-wave velocity model was achieved by the solution of an inverse problem-based dispersion of surface waves and propagation in a vertically heterogeneous medium. The 2D genetic algorithm was employed for the inversion of dispersion curves to obtain velocity and thickness of subsurface layers. The depth to engineering bedrock and velocity of shear waves in the first 30 m was deciphered and mapped. The depth of bedrock in study area varies from 4 to 30 m, and V _{S 30} ranges from 320 to 800 m/s. The most of study area falls in B and C class categories in addition to few sites of D class according to the NEHRP guidelines.     

GIS-based multi-criteria earthquake hazards evaluation using analytic hierarchy process for a nuclear power plant site,west Alexandria,Egypt

Nuclear power plants are designed to prevent the hazardous effects of the earthquakes and any external events to keep the safety of the plant. Ninety-one shallow seismic refraction profiles were performed to determine shear wave velocity of the engineering layers at the site of El Dabaa area that is situated to the northern coastline of Egypt for seismic hazard microzonation evaluation according to hazard index values. A microzonation is a procedure of delineating an area into individual zones having different ranks of numerous seismic hazards. This will aid in classifying areas of high seismic risk which is vigorous for industrial design of nuclear structures. The site response analysis requires the characterization of subsurface materials considering local subsurface profiles of the site. Site classification of the area under investigation was undertaken using P- and S-waves and available borehole data. The studied nuclear power plant site has been characterized as per NEHRP site classification using an average velocity of transverse wave (V _s ³⁰ ) of depth 30 m which acquired from seismic survey. This site was categorized into two site classes: the major one is “site class B,” and the minor one is “site class A.” The attenuation coefficient, the damping ratio and the liquefaction potential are geotechnical parameters which were derived from P- and S-waves, and have their major effects on the seismic hazard contribution. 1D ground response analysis was carried out in the places of seismic profiles inside the site for estimating the amount of ground quaking using peak ground acceleration (PGA), site amplification, predominant frequency and spectral accelerations on the surface of ground by the DEEPSOIL software package. Seven factors (criteria) deliberated to assess the earthquake hazard index map are: (1) the peak ground acceleration at the bedrock, (2) the amplification of the site, (3) the liquefaction potential, (4) the main frequency of the earthquake signal, (5) the average V _s of the first 30 m from the ground surface, (6) the depth to the groundwater and (7) the depth to the bedrock. These features were exemplified in normalized maps after uniting them to 0–1 scores according to some criteria by the minimum and maximum values as linear scaling points. Multi-criteria evaluation is an application of multi-criteria decision analysis theory that used for developing a seismic hazard index map for a nuclear power plant site at El Dabaa area in ArcGIS 10.1 software. Two models of decision making were used in this work for seismic hazard microzonation. The analytic hierarchy process model was applied to conduct the relative weights of the criteria by pairwise comparison using Expert Choice Software. An earthquake hazard index map was combined using Weighted Linear Combination model of the raster weighted overlay tool of ArcGIS 10.1. The results indicated that most of the study site of the nuclear power plant is a region of low to moderate hazard; its values are ranging between 0.2 and 0.4.     

A newly developed seismic microzonation model of Erbaa (Tokat, Turkey) located on seismically active eastern segment of the North Anatolian Fault Zone (NAFZ)

A methodology to model seismic microzonation maps is required in the hazard mitigation decision plans of the earthquake prone areas. The stage of disaster preparedness for new residential places is of great importance for detailed seismic microzonation models. The effects of local geological and geotechnical site conditions were considered in order to establish site characterization as the initial stage of the models in this study. Dynamic soil properties based on the empirical correlations between shear wave velocity (V _s) and standard penetration test blow counts were taken into account in order to define representative soil profiles extending down to the engineering bedrock. One-dimensional site response analyses were performed to analyze earthquake characteristics on the ground surface. The layers for soil classification, geology, depth to groundwater level, amplification, distance to fault, slope and aspect, and liquefaction-induced ground deformation potential of the study area were prepared in seismic microzonation models. The study area, Erbaa, is placed along the seismically active North Anatolian Fault Zone. Final seismic microzonation map of the study area was evaluated applying different GIS-based Multi-Criteria Decision Analysis (MCDA) techniques. Two of the MCDA techniques, simple additive weighting and analytical hierarchical process (AHP), are considered during the evaluation step of the final seismic microzonation map. The comparison is made in order to distinguish two different maps based on these MCDA techniques. Eventually, AHP-based seismic microzonation map is more preferable for the seismic design purposes in this study.     

Surface wave studies for shear wave velocity and bedrock depth estimation over basalts

Shear wave velocity (V _S) estimation is of paramount importance in earthquake hazard assessment and other geotechnical/geo engineering studies. In our study, the shear wave velocity was estimated from ground roll using multichannel analysis of surface wave (MASW) technique making use of dispersive characteristics of Rayleigh type surface waves followed by imaging the shallow subsurface basaltic layers in an earthquake-prone region near Jabalpur, India. The reliability of MASW depends on the accurate determination of phase velocities for horizontally traveling fundamental mode Rayleigh waves. Inversion of data from surface waves resulted in a shear wave velocity (V _S) in the range of 200–1,200 m/s covering the top soil to weathering and up to bedrock corresponding to a depth of 10–30 m. The P-wave velocity (V _P) obtained from refraction seismic studies at these locations found to be comparable with V _S at an assumed specific Poisson’s ratio. A pair of selected set of V _S profiles over basalt which did not result in a hazardous situation in an earthquake of moderate magnitude are presented here as a case study; in other words, the shear wave velocity range of more than 200 m/s indicate that the area is highly unlikely prone to liquefaction during a moderate or strong earthquake. The estimated depth to basalt is found to be 10–12 m in both the cases which is also supported by refraction studies.     

Preliminary results for a semi-automated quantification of site effects using geomorphometry and ASTER satellite data for Mozambique,Pakistan and Turkey

Estimation of the degree of local seismic wave amplification (site effects) requires precise information about the local site conditions. In many regions of the world, local geologic information is either sparse or is not readily available. Because of this, seismic hazard maps for countries such as Mozambique, Pakistan and Turkey are developed without consideration of site factors and, therefore, do not provide a complete assessment of future hazards. Where local geologic information is available, details on the traditional maps often lack the precision (better than 1:10,000 scale) or the level of information required for modern seismic microzonation requirements. We use high-resolution (1:50,000) satellite imagery and newly developed image analysis methods to begin addressing this problem. Our imagery, consisting of optical data and digital elevation models (DEMs), is recorded from the ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) sensor system. We apply a semi-automated, object-oriented, multi-resolution feature segmentation method to identify and extract local terrain features. Then we classify the terrain types into mountain, piedmont and basin units using geomorphometry (topographic slope) as our parameter. Next, on the basis of the site classification schemes from the Wills and Silva (1998) study and the Wills et al (2000) and Wills and Clahan (2006) maps of California, we assign the local terrain units with V _s 30 (the average seismic shear-wave velocity through the upper 30m of the subsurface) ranges for selected regions in Mozambique, Pakistan and Turkey. We find that the applicability of our site class assignments in each region is a good first-approximation for quantifying local site conditions and that additional work, such as the verification of the terrain’s compositional rigidity, is needed.     

Dynamic soil properties for microzonation of Delhi,India

Delhi, the capital of India, has experienced mild seismic shaking during several earthquakes in the past. The large variations of depth to bedrock and ground water table coupled with different soil types at different locations of Delhi necessitate a seismic microzonation study. Dynamic soil properties such as shear wave velocity, modulus reduction and damping characteristics of local soils are the basic and essential input parameters for conducting even a preliminary ground response analysis which is an essential input in microzonation studies. Shear wave velocity is not measured routinely due to its high cost and lack of the required expertise. Several researchers in the past developed correlations between shear wave velocity (V _s ) and routinely measured N values. In the present study, shear wave velocity profiles measured in the field at more than 80 borehole locations to a depth of about 20 to 32m using Spectral Analysis of Surface Waves (SASW) are presented and correlations between shear wave velocity and N values are also presented for use by engineers and designers. Results of strain and stress controlled cyclic triaxial tests on remoulded samples of sand-silt mixtures in the high strain range are used for generating the modulus reduction and damping curves and are compared with the well-known curves in the literature. The results presented in this article can be used for microzonation studies as well as site specific ground response analyses at Delhi.     

Assessment of potential seismic hazard and site effect in Antakya (Hatay Province), SE Turkey

Antakya city is at risk because of strong earthquakes occurring in the area, and different soil conditions that can produce variation of the ground motion amplification. Microzonation of cities provides a basis for site-specific hazard analysis in urban settlements. In particular, seismic microzonation can be provided by means of detailed seismic assessment of the area, including earthquake recordings and geological studies. In this paper, we propose a preliminary microzonation map for the city of Antakya, based on the variation of the dominant periods and shear velocities of the sediments covering the area. The periods are retrieved from microtremor measurements conducted at 69 sites, using the horizontal-to-vertical spectral ratio technique. The results of microtremor analysis were compared with data obtained from refraction microtremor (ReMi) measurements at four profiles crossing the studied area. According to the classification of dominant periods, Antakya city can be divided into five zones, probably prone to different levels of seismic hazard. The shorter natural periods are in inner Antakya and both the sides of Asi River (i.e., northern and southern parts). The eastern and western parts of Antakya have maximum dominant periods. The V _s ³⁰ values were calculated by using the ReMi method along the profiles. Antakya city has V _s ³⁰ values in the range of category C of the national earthquake hazard reduction programme site classification.     

Seismic microzonation of Bangalore,India

In the present study, an attempt has been made to evaluate the seismic hazard considering local site effects by carrying out detailed geotechnical and geophysical site characterization in Bangalore, India to develop microzonation maps. An area of 220 km², encompassing Bangalore Mahanagara Palike (BMP) has been chosen as the study area. Seismic hazard analysis and microzonation of Bangalore are addressed in three parts: in the first part, estimation of seismic hazard is done using seismotectonic and geological information. Second part deals with site characterization using geotechnical and shallow geophysical techniques. In the last part, local site effects are assessed by carrying out one-dimensional (1-D) ground response analysis (using the program SHAKE2000) using both standard penetration test (SPT) data and shear wave velocity data from multichannel analysis of surface wave (MASW) survey. Further, field experiments using microtremor studies have also been carried out for evaluation of predominant frequency of the soil columns. The same has been assessed using 1-D ground response analysis and compared with microtremor results. Further, the Seed and Idriss simplified approach has been adopted to evaluate the soil liquefaction susceptibility and liquefaction resistance assessment. Microzonation maps have been prepared with a scale of 1:20,000. The detailed methodology, along with experimental details, collated data, results and maps are presented in this paper.     

Examination of the structural characteristics arising in gypsums by the GPR and MASW methods (Sivas,Turkey)

This study was performed at an area of 50?×?48 m² being defined as a new settlement in the northeast of Sivas. In the study, the discontinuities that are not deep and their geophysical characteristics were examined by the GPR and MASW methods. For interpretation, GPR cross sections were prepared as 2D–3D, and MASW cross sections were prepared as 2D. As for geophysical cross sections, about 10 m depth was examined. It was understood that the reflections observed in the form of hyperbolas in GPR cross sections correspond to areas having low S wave velocity (V_s) in MASW cross sections. It was understood that the S wave velocities are lower than 653 m/s, that the seismic velocities in between 653 and 275 m/s indicate partially deteriorated areas and that the S wave velocities of unweathered gypsums are higher than 1275 m/s at these low-velocity zones. Thus, it was thought that the fill material that may arise in the fracture, crack and deterioration areas arises from intercalation and clastic gypsum units, and that it plays a role in having low value S wave velocities. In all the geophysical cross sections, it was understood that the structures with gypsum are intense at the initial 5 m. And a fracture at the south of the study area, that it was estimated might be longer than 40 m, was determined as the largest gypsum structure. It was understood that this fracture starts from a depth of about 5 m in the west and that it slopes down to 7 m depth in the east. According to these results, it was understood that the damage amount arising in time in the gypsum structures from the effect of water may increase, the study area was defined as risky, and the required importance should be attached to these structures especially in foundation engineering.     

Estimation of shearwave velocity in drifts using multichannel analysis of surface wave (MASW) technique — A case study from Jammu & Kashmir,India

Multichannel analysis of surface waves (MASW) is a non-destructive seismic prospecting method utilizing Rayleigh waves for imaging and characterizing shallow sub-surface structure. Multichannel analysis of surface waves (MASW) studies were conducted in drift areas of two bridge sites in the hilly terrain of J&K for imaging and characterizing shallow sub-surface structure. The purpose of the present study is to estimate the shear wave velocity (V_S) and subsurface structure in four drifts made in a hilly terrain for construction of two bridges. Rayleigh waves are having dispersive properties, travelling along or near the ground surface and are usually characterized by relatively low velocity, low frequency, and high amplitude. The study area comprises of Tertiary group of rocks which are underlain by Siwalik group. The main rock type in the study area is dolomite which has undergone various geological processes like weathering, jointing, fracturing and shearing. MASW data was collected inside four drifts in the mountainous terrain of J&K state which are located on either sides of Chenab river. The data was analyzed by relevant processing software using dispersion and inversion technique. Shear wave velocities were estimated up to 30 m depth. Average shear wave velocity (V_S ³⁰) up to top 30m was also computed. It is observed that, V_S in the range 400–800 m/s upto 10–15 m corresponding to weathered rock, followed by compact dolomite rock up to the depth of about 30 m with V_S in the range 1200–1600 m/s. Some low velocity zones are also identified from these sections which represent shear zones.     

Assessing shear wave velocity profiles using multiple passive techniques of Shillong region of northeast India

Shillong region of northeast India falls under seismic zone V. Being a commercial hub, urbanization in this study region is at its peak. In order to qualitatively assess the subsurface velocity profiling of this area, we have blended two robust techniques, namely spatial autocorrelation (SPAC) and frequency wavenumber (FK) method. Corresponding to array noise data collected at five sites, situated in the close proximity of boreholes, we have evaluated V_S and V_P as well. The shear wave velocity estimates yielded by these techniques are found to be consistent with each other. The computed V_s values up to depth of 30 m are observed to be in the range of 275–375 m/s, in most of the sites which implies prevalence of low-velocity zone at some pocket areas. The estimates so found are systematically analyzed and implications are outlined. The results were corroborated by substantial evidence of site geology as well as geotechnical information.     

Applications of a GIS-based geotechnical tool to assess spatial earthquake hazards in an urban area

A geotechnical information system (GTIS) was constructed within a spatial geographic information system (GIS) framework to reliably predict geotechnical information and accurately estimate site effects at Gyeongju, an urban area in South Korea. The system was built based on both collected and performed site investigation data in addition to acquired geo-knowledge data. Seismic zoning maps were constructed using the site period (T _G) and mean shear wave velocity to a depth of 30 m (V _S30), and these maps were presented as a regional strategy to mitigate earthquake-induced risks in the study area. In particular, the T _G distribution map indicated the susceptibility to ground motion resonance in periods ranging from 0.2 to 0.5 s and the corresponding seismic vulnerability of buildings with two to five stories. Seismic zonation of site classification according to V _S30 values was also performed to determine the site amplification coefficients for seismic design and seismic performance evaluation at any site and administrative subunit in the study area. In addition, we investigated the site effects according to subsurface and surface ground irregularities at Gyeongju by seismic response analyses in time domains based on both two- and three-dimensional spatial finite element models, which were generated using spatial interface coordinates between geotechnical subsurface layers predicted by the GTIS. This practical study verified that spatial GIS-based geotechnical information can be a very useful resource in determining how to best mitigate seismic hazards, particularly in urban areas.     

Geospatial contour mapping of shear wave velocity for Mumbai city

Shear wave velocity is one of the most important input parameter in the analysis of geotechnical earthquake engineering problems, particularly to estimate site-specific amplification factor and ground response study. Dynamic in situ tests such as spectral analysis of surface waves (SASW) or multichannel analysis of surface waves (MASW) are very expensive. Also due to lack of specialized personnel, these tests are generally avoided in many soil investigation programs. Worldwide, several researchers have developed correlations between the SPT ‘N’ value and shear wave velocity ‘V _s’, which are useful for determining the dynamic soil properties. In the present study, more than 400 numbers of soil borehole data were collected from various geotechnical investigation agencies, government engineering institutes and geotechnical laboratories from different parts of Mumbai city, which is financial capital of India with highest population density. In this paper, an attempt has been made to develop the correlation between the SPT ‘N’ value and shear wave velocity ‘V _s’ for various soil profile of Mumbai city and compared with other existing correlations for different cities in India. Using Geographical Information System (GIS), a geospatial contour map of shear wave velocity profile for Mumbai city is prepared with contour intervals of 25 and 50 m/s. The scarcity of database or maps of shear wave velocity profile for Mumbai city will make the present geospatial contour maps extremely useful and beneficial to the designer, practitioners for seismic hazard study involved in geotechnical earthquake engineering.     

Use of Surface Waves in Statistical Correlations of Shear Wave Velocity and Penetration Resistance of Chennai Soils

Shear wave velocity (V _s) is one of the most important input parameter to represent the stiffness of the soil layers. It is preferable to measure V _s by in situ wave propagation tests, however it is often not economically feasible to perform the tests at all locations. Hence, a reliable correlation between V _s and standard penetration test blow counts (SPT-N) would be a considerable advantage. This paper presents the development of empirical correlations between V _s and SPT-N value for different categories of soil in Chennai city characterized by complex variation of soil conditions. The extensive shear wave velocity measurement was carried out using Multichannel Analysis of Surface Waves (MASW) technique at the sites where the SPT-N values are available. The bender element test is performed to compare the field MASW test results for clayey soils. The correlations between shear wave velocity and SPT-N with and without energy corrections were developed for three categories of soil: all soils, sand and clay. The proposed correlations between uncorrected and energy corrected SPT-N were compared with regression equations proposed by various other investigators and found that the developed correlations exhibit good prediction performance. The proposed uncorrected and energy corrected SPT-N relationships show a slight variation in the statistical analysis indicating that both the uncorrected and energy corrected correlations can predict shear wave velocity with equal accuracy. It is also found that the soil type has a little effect on these correlations below SPT-N value of about 10.     

Spatial distribution of soil shear-wave velocity and the fundamental period of vibration – a case study of the Saguenay region,Canada

The spatial distribution of soil shear-wave velocity and the fundamental period of vibration were selected as input parameters for the determination of potential seismic site effects in the Saguenay region, Canada. The methodology used in this study involved three clear steps. First, a 3D geological model of the surficial deposits was built taking into consideration the type, spatial distribution and thickness of the deposits. Second, representative average V_s values were determined for each of the major soil units. Finally, the average shear-wave velocity from the ground surface to bedrock (V_sav), the shear-wave velocity of the upper 30 m (V_s30) and the fundamental site resonance period (T₀) were calculated over a regular grid for the study area. The results include the spatial distribution of the fundamental site resonance period, the average shear-wave velocity in the first 30 m of the ground and the spatial distribution of National Building Code of Canada seismic soil classes for the Saguenay region.     

Local site effect of a clay site in Shiraz based on seismic hazard of Shiraz Plain

Propagation of seismic waves through soil layers would drastically change the frequency content and amplitude-based features of ground motions at the surface. These alterations are known as seismic site effects. Computation of site effects of high-populated areas such as large cities is of great importance (e.g., it is used in development of seismic microzonation of a region). Shiraz is one of the most populous cities of Iran and is located in a high seismic hazardous region. A representative clay site in this city is selected to assess local site effects. The time series and random vibration theory procedure in the frequency domain are implemented to analyze the aforementioned site. Furthermore, the nonlinear dynamic soil behavior is simulated by the equivalent linear method and the nonlinear method via DEEPSOIL program. Three types of soil column uncertainties such as shear wave velocity, modulus reduction, and damping ratio of soil layers as well as depth of underlying rock half-space (D _bed) are considered herein. The mean amplification and standard deviation of natural logarithm of amplification factors are computed for a variety of analysis types. The results of the current study show that the computed mean and standard deviation of amplification factor in ln units by considering only V _S uncertainty are in good agreement with the corresponding ones by considering V _S and modulus reduction and damping ratio variabilities simultaneously for the studied site. Furthermore, it seems that the effect of bedrock depth in definition of spectral shapes of the Iranian seismic building code should be taken into account.     

Variations in shear wave velocity and soil site class in Kolkata city using regression and sensitivity analysis

The detrimental effects of an earthquake are strongly influenced by the response of soils subjected to dynamic loading. The behavior of soils under dynamic loading is governed by the dynamic soil properties such as shear wave velocity, damping characteristics and shear modulus. Worldwide, it is a common practice to obtain shear wave velocity (V _s in m/s) using the correlation with field standard penetration test (SPT) N values in the absence of sophisticated dynamic field test data. In this paper, a similar but modified advanced approach has been proposed for a major metro city of eastern India, i.e., Kolkata city (latitudes 22°20′N–23°00′N and longitudes 88°04′E–88°33′E), to obtain shear wave velocity profile and soil site classification using regression and sensitivity analyses. Extensive geotechnical borehole data from 434 boreholes located across 75 sites in the city area of 185 km² and laboratory test data providing information on the thickness of subsoil strata, SPT N values, consistency indices and percentage of fines are collected and analyzed thoroughly. A correlation between shear wave velocity (V _s) and SPT N value for various soil profiles of Kolkata city has been established by using power model of nonlinear regression analysis and compared with existing correlations for other Indian cities. The present correlations, having regression coefficients (R ²) in excess of 0.96, indicated good prediction capability. Sensitivity analysis predicts that significant influence of soil type exists in determining V _s values, for example, typical silty sand shows 30.4 % increase in magnitude of V _s as compared to silt of Kolkata city. Moreover, the soil site classification shows Class D and Class E category of soil that exists typically in Kolkata city as per NEHRP (Recommended provisions for seismic regulations for new buildings and other structures—Part 1: Provisions. Prepared by the Building Seismic Safety Council for the Federal Emergency Management Agency (Report FEMA 450), Washington, DC, 2003) guidelines and thereby highlighting the seismic vulnerability of the city. The results presented in this study can be utilized for seismic microzonation, ground response analysis and hazard assessment for Kolkata city.

     

Earthquake hazard in Northeast India — A seismic microzonation approach with typical case studies from Sikkim Himalaya and Guwahati city

A comprehensive analytical as well as numerical treatment of seismological, geological, geomorphological and geotechnical concepts has been implemented through microzonation projects in the northeast Indian provinces of Sikkim Himalaya and Guwahati city, representing cases of contrasting geological backgrounds — a hilly terrain and a predominantly alluvial basin respectively. The estimated maximum earthquakes in the underlying seismic source zones, demarcated in the broad northeast Indian region, implicates scenario earthquakes of M _W 8.3 and 8.7 to the respective study regions for deterministic seismic hazard assessments. The microzonation approach as undertaken in the present analyses involves multi-criteria seismic hazard evaluation through thematic integration of contributing factors. The geomorphological themes for Sikkim Himalaya include surface geology, soil cover, slope, rock outcrop and landslide integrated to achieve geological hazard distribution. Seismological themes, namely surface consistent peak ground acceleration and predominant frequency were, thereafter, overlaid on and added with the geological hazard distribution to obtain the seismic hazard microzonation map of the Sikkim Himalaya. On the other hand, the microzonation study of Guwahati city accounts for eight themes — geological and geomorphological, basement or bedrock, landuse, landslide, factor of safety for soil stability, shear wave velocity, predominant frequency, and surface consistent peak ground acceleration. The five broad qualitative hazard classifications — ‘low’, ‘moderate’, ‘high’, ‘moderate high’ and ‘very high’ could be applied in both the cases, albeit with different implications to peak ground acceleration variations. These developed hazard maps offer better representation of the local specific seismic hazard variation in the terrain.     

Microtremor study for evaluating the site response characteristics in the Surat City of western India

The Surat City, which is the second most populated city in the state of Gujarat in western India, warrants site-specific seismic hazard assessment due to its rapid urbanization and proximity to major seismogenic zones. This study reports results of microtremor investigations at 72 single stations and 4 arrays in an area of 325 km² spanning the city. The resonant frequencies, associated peak amplification values and liquefaction vulnerability indices were deduced from the horizontal to vertical spectral ratios. Ground amplification (A_HVSR) in the range of 3.0–5.0 was observed in the 2.0–4.0-Hz frequency band at most of the sites. A secondary A_HVSR between 2.0 and 3.0 is also observed in the 6.0–7.0-Hz frequency band at a few sites. Locales that are most susceptible to liquefaction are identified based on their vulnerability index (K _g) exceeding the value of 10. The shear wave velocities (V _s) ≥ 500 m/s inferred from array measurements occur at 38 m depth in the western part and ~16 m depth in the eastern part of city. The response spectra estimated from strong motion data recorded at an accelerograph site in Surat from three earthquakes of M _w ≥ 3.2 that occurred in Kachchh, Saurashtra and Narmada regions are in accordance with our inferences of characteristic site frequencies and amplification. Our results, in agreement with the damage scenario during the 2001 Bhuj earthquake, provide valuable inputs for site-specific seismic hazard evaluation of the Surat City.     

Seismic hazard assessment and mitigation in India: an overview

The Indian subcontinent is characterized by various tectonic units viz., Himalayan collision zone in North, Indo-Burmese arc in north-east, failed rift zones in its interior in Peninsular Indian shield and Andaman Sumatra trench in south-east Indian Territory. During the last about 100 years, the country has witnessed four great and several major earthquakes. Soon after the occurrence of the first great earthquake, the Shillong earthquake (M _w: 8.1) in 1897, efforts were started to assess the seismic hazard in the country. The first such attempt was made by Geological Survey of India in 1898 and since then considerable progress has been made. The current seismic zonation map prepared and published by Bureau of Indian Standards, broadly places seismic risk in different parts of the country in four major zones. However, this map is not sufficient for the assessment of area-specific seismic risks, necessitating detailed seismic zoning, that is, microzonation for earthquake disaster mitigation and management. Recently, seismic microzonation studies are being introduced in India, and the first level seismic microzonation has already been completed for selected urban centres including, Jabalpur, Guwahati, Delhi, Bangalore, Ahmadabad, Dehradun, etc. The maps prepared for these cities are being further refined on larger scales as per the requirements, and a plan has also been firmed up for taking up microzonation of 30 selected cities, which lie in seismic zones V and IV and have a population density of half a million. The paper highlights the efforts made in India so far towards seismic hazard assessment as well as the future road map for such studies.