Frequency response analysis of layered ground in Ahmedabad during Bhuj earthquake

This paper presents the results of seismic response analysis of layered ground in Ahmedabad City during the earthquake in Bhuj on 26(th) January 2001. An attempt has been made to understand the reasons for the failure of multistoreyed buildings founded on soft alluvial deposits in Ahmedabad. Standard Penetration test at a site very close to the Sabarmati river belt was carried out for geotechnical investigations. The program SHAKE91, widely used in the field of earthquake engineering for computing the seismic response of horizontally layered soil deposits, was used to analyse the soil profile at the selected site considering the ground as one dimensional layered elastic system. The ground accelerations recorded at the ground floor of the Regional Passport Staff Quarters building, which is very close to the investigated site, was used as input motion. Also, Finite Element Analysis was carried out for different configurations of multistorey building frames for evaluating their natural frequencies and is compared with the predominant frequency of the layered soil system. The results reveal that the varying degree of damage to multistorey buildings in the close proximity of Sabarmati river area was essentially due to the large amplification of the ground and the near resonance condition.

Review of correlations between SPT N and shear modulus: A new correlation applicable to any region

A low strain shear modulus plays a fundamental role in earthquake geotechnical engineering to estimate the ground response parameters for seismic microzonation. A large number of site response studies are being carried out using the standard penetration test (SPT) data, considering the existing correlation between SPT N values and shear modulus. The purpose of this paper is to review the available empirical correlations between shear modulus and SPT N values and to generate a new correlation by combining the new data obtained by the author and the old available data. The review shows that only few authors have used measured density and shear wave velocity to estimate shear modulus, which were related to the SPT N values. Others have assumed a constant density for all the shear wave velocities to estimate the shear modulus. Many authors used the SPT N values of less than 1 and more than 100 to generate the correlation by extrapolation or assumption, but practically these N values have limited applications, as measuring of the SPT N values of less than 1 is not possible and more than 100 is not carried out. Most of the existing correlations were developed based on the studies carried out in Japan, where N values are measured with a hammer energy of 78%, which may not be directly applicable for other regions because of the variation in SPT hammer energy. A new correlation has been generated using the measured values in Japan and in India by eliminating the assumed and extrapolated data. This correlation has higher regression coefficient and lower standard error. Finally modification factors are suggested for other regions, where the hammer energy is different from 78%. Crown Copyright (C) 2012 Published by Elsevier Ltd. All rights reserved.

Recent microtremors near the Idukki Reservoir, Kerala, South India

The continuing low-level seismicity in the vicinity of the Idukki Reservoir, Kerala, is interesting from the perspective of hydrologically triggered earthquakes. While the frequency of triggered earthquakes in the vicinity of a reservoir usually reduces with time and the largest earthquake usually occurs within a few years on the initial filling, the triggered seismicity in the proximity of the Idukki Reservoir seems to be showing a second, delayed peak, as the 1977 (M 3.5) tremor was followed by a slightly larger event in 2011, 24 years after the first burst of activity. Quite unprecedented in the context of reservoir-triggered sequences, we consider this delayed sequence as the hydrologic response of a critically stressed hypocentral region, to monsoonal recharging. The sustained activity several decades after the impoundment and the temporal relation with the monsoon suggest that at least some parts of the reservoir region continue to retain the potential for low-level seismic activity in response to hydrologic cycles.

The hazard potential of the western segment of the Makran subduction zone, northern Arabian Sea

Evaluating the hazard potential of the Makran subduction zone requires understanding the previous records of the large earthquakes and tsunamis. We address this problem by searching for earthquake and tectonic proxies along the Makran Coast and linking those observations with the available constraints on historical seismicity and the tell-tale characteristics of sea floor morphology. The earthquake of Mw 8.1 of 1945 and the consequent tsunami that originated on the eastern part of the Makran are the only historically known hazardous events in this region. The seismic status of the western part of the subduction zone outside the rupture area of the 1945 earthquake remains an enigma. The near-shore shallow stratigraphy of the central part of Makran near Chabahar shows evidence of seismically induced liquefaction that we attribute to the distant effects of the 1945 earthquake. The coastal sites further westward around Jask are remarkable for the absence of liquefaction features, at least at the shallow level. Although a negative evidence, this possibly implies that the western part of Makran Coast region may not have been impacted by near-field large earthquakes in the recent past-a fact also supported by the analysis of historical data. On the other hand, the elevated marine terraces on the western Makran and their uplift rates are indicative of comparable degree of long-term tectonic activity, at least around Chabahar. The offshore data suggest occurrences of recently active submarine slumps on the eastern part of the Makran, reflective of shaking events, owing to the great 1945 earthquake. The ocean floor morphologic features on the western segment, on the contrary, are much subdued and the prograding delta lobes on the shelf edge also remain intact. The coast on the western Makran, in general, shows indications of progradation and uplift. The various lines of evidence thus suggest that although the western segment is potentially seismogenic, large earthquakes have not occurred there in the recent past, at least during the last 600 years. The recurrence period of earthquakes may range up to 1,000 years or more, an assessment based on the age of the youngest dated coastal ridge. The long elapsed time points to the fact that the western segment may have accumulated sufficient slip to produce a major earthquake.

On the Compatibility of a Flux Transport Dynamo with a Fast Tachocline Scenario

The compatibility of the fast-tachocline scenario with a flux-transport dynamo model is explored. We employ a flux-transport dynamo model coupled with simple feedback formulae relating the thickness of the tachocline to the amplitude of the magnetic field or to the Maxwell stress. The dynamo model is found to be robust against the nonlinearity introduced by this simplified fast-tachocline mechanism. Solar-like butterfly diagrams are found to persist and, even without any parameter fitting, the overall thickness of the tachocline is well within the range admitted by helioseismic constraints. In the most realistic case of a time-and latitude-dependent tachocline thickness linked to the value of the Maxwell stress, both the thickness and its latitudinal dependence are in excellent agreement with seismic results. In nonparametric models, cycle-related temporal variations in tachocline thickness are somewhat larger than admitted by helioseismic constraints; we find, however, that introducing a further parameter into our feedback formula readily allows further fine tuning of the thickness variations.

Strain localization effect on system reliability based design of bridge abutments under earthquake loading

The component and system reliability based design of bridge abutments under earthquake loading is presented in the paper. Planar failure surface has been used in conjunction with pseudo-dynamic approach to compute seismic active earth pressures on an abutment. The pseudo-dynamic method, considers the effect of phase difference in shear waves, soil amplification along with the horizontal seismic accelerations, strain localization in backfill soil and associated post-peak reduction in the shear resistance from peak to residual values along a previously formed failure plane. Four modes of stability viz. sliding, overturning, eccentricity and bearing capacity of the foundation soil are considered in the analysis. The series system reliability is computed with an assumption of independent failure modes. The lower and upper bounds of system reliability are also computed by taking into account the correlations between four failure modes, which is evaluated using the direction cosines of the tangent planes at the most probable points of failure.

Influence of strain localization on reliability based design of bridge abutments using pseudo-dynamic method

The paper focuses on reliability based design of bridge abutments when subjected to earthquake loading. Planar failure surface has been used in conjunction with pseudo-dynamic approach to compute the seismic active earth pressures on the bridge abutment. The proposed pseudo dynamic method, considers the effects of strain localization in the backfill soil and associated post-peak reduction in the shear resistance from peak to residual values along a previously formed failure plane, phase difference in shear waves and soil amplification along with the horizontal seismic accelerations. Four modes of stability viz. sliding, overturning, eccentricity and bearing capacity of the foundation soil are considered for the reliability analysis. The influence of various design parameters on the seismic reliability indices against four modes of failure is presented, following the suggestions of Japan Road Association, Caltrans Bridge Design Specifications and U.S Department of the Army.

Estimation of design basis earthquake using region-specific M-max, for the NPP site at Kalpakkam, Tamil Nadu, India

The objective of the paper is to estimate Safe Shutdown Earthquake (SSE) and Operating/Design Basis Earthquake (OBE/DBE) for the Nuclear Power Plant (NPP) site located at Kalpakkam, Tamil Nadu, India. The NPP is located at 12.558 degrees N, 80.175 degrees E and a 500 km circular area around NPP site is considered as `seismic study area' based on past regional earthquake damage distribution. The geology, seismicity and seismotectonics of the study area are studied and the seismotectonic map is prepared showing the seismic sources and the past earthquakes. Earthquake data gathered from many literatures are homogenized and declustered to form a complete earthquake catalogue for the seismic study area. The conventional maximum magnitude of each source is estimated considering the maximum observed magnitude (M-max(obs)) and/or the addition of 0.3 to 0.5 to M-max(obs). In this study maximum earthquake magnitude has been estimated by establishing a region's rupture character based on source length and associated M-max(obs). A final source-specific M-max is selected from the three M-max values by following the logical criteria. To estimate hazard at the NPP site, ten Ground-Motion Prediction Equations (GMPEs) valid for the study area are considered. These GMPEs are ranked based on Log-Likelihood (LLH) values. Top five GMPEs are considered to estimate the peak ground acceleration (PGA) for the site. Maximum PGA is obtained from three faults and named as vulnerable sources to decide the magnitudes of OBE and SSE. The average and normalized site specific response spectrum is prepared considering three vulnerable sources and further used to establish site-specific design spectrum at NPP site.

Constraining large earthquakes along the Andaman trench using deepwater turbidites: prospects and challenges

The 2004 Sumatra-Andaman earthquake was unprecedented in terms of its magnitude (M-w 9.2), rupture length along the plate boundary (1300 km) and size of the resultant tsunami. Since 2004, efforts are being made to improve the understanding of the seismic hazard in the Sumatra-Andaman subduction zone in terms of recurrence patterns of major earthquakes and tsunamis. It is reasonable to assume that previous earthquake events in the Myanmar Andaman segment must be preserved in the geological record in the form of seismo-turbidite sequences. Here we present the prospects of conducting deep ocean palaeoseismicity investigations in order to refine the quantification of the recurrence pattern of large subduction-zone earthquakes along the Andaman-Myanmar arc. Our participation in the Sagar Kanya cruise SK-273 (in June 2010) was to test the efficacy of such a survey. The primary mission of the cruise, along a short length (300 km) of the Sumatra Andaman subduction front was to collect bathymetric data of the ocean floor trenchward of the Andaman Islands. The agenda of our piggyback survey was to fix potential coring sites that might preserve seismo-turbidite deposits. In this article we present the possibilities and challenges of such an exercise and our first-hand experience of such a preliminary survey. This account will help future researchers with similar scientific objectives who would want to survey the deep ocean archives of this region for evidence of extreme events like major earthquakes.

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

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

Assessment of Seismically Induced Landslide Hazard for the State of Karnataka Using GIS Technique

Landslide hazards are a major natural disaster that affects most of the hilly regions around the world. In India, significant damages due to earthquake induced landslides have been reported in the Himalayan region and also in the Western Ghat region. Thus there is a requirement of a quantitative macro-level landslide hazard assessment within the Indian subcontinent in order to identify the regions with high hazard. In the present study, the seismic landslide hazard for the entire state of Karnataka, India was assessed using topographic slope map, derived from the Digital Elevation Model (DEM) data. The available ASTER DEM data, resampled to 50 m resolution, was used for deriving the slope map of the entire state. Considering linear source model, deterministic seismic hazard analysis was carried out to estimate peak horizontal acceleration (PHA) at bedrock, for each of the grid points having terrain angle 10A degrees and above. The surface level PHA was estimated using nonlinear site amplification technique, considering B-type NEHRP site class. Based on the surface level PHA and slope angle, the seismic landslide hazard for each grid point was estimated in terms of the static factor of safety required to resist landslide, using Newmark's analysis. The analysis was carried out at the district level and the landslide hazard map for all the districts in the Karnataka state was developed first. These were then merged together to obtain a quantitative seismic landslide hazard map of the entire state of Karnataka. Spatial variations in the landslide hazard for all districts as well as for the entire state Karnataka is presented in this paper. The present study shows that the Western Ghat region of the Karnataka state is found to have high landslide hazard where the static factor of safety required to resist landslide is very high.

Recent Seismicity in India and Adjoining Regions

In the present work, historical and instrumental seismicity data of India and its adjoining areas (within 300km from Indian political boundary) are compiled to form the earthquake catalog for the country covering the period from 1505 to 2009. The initial catalogue consisted of about 139563 earthquake events and after declustering,the total number of events obtained was 61315. Region specific earthquake magnitude scaling relations correlating different magnitude scales were achieved and a homogenous earthquake catalogue in moment magnitude (MW) scale was developed for the region. This paper also presents the results of the use of Geographic Information Systems (GIS) to prepare a digitized seismic source map of India. The latest earthquake data were superimposed on the digitized source map to get a final Seismotectonic map of India. The study area has been divided into 1225 grid points (approximately 110km×110km) and the seismicity analysis has been done to get the spatial variation of seismicity parameters ‘a’ and ‘b’ across the country. The homogenized earthquake catalogue with the event details is listed in the website http://civil.iisc.ernet.in/~sreevals/resource.htm

Spatial variation of seismicity parameters across India and adjoining areas

An attempt has been made to quantify the variability in the seismic activity rate across the whole of India and adjoining areas (0–45°N and 60–105°E) using earthquake database compiled from various sources. Both historical and instrumental data were compiled and the complete catalog of Indian earthquakes till 2010 has been prepared. Region-specific earthquake magnitude scaling relations correlating different magnitude scales were achieved to develop a homogenous earthquake catalog for the region in unified moment magnitude scale. The dependent events (75.3%) in the raw catalog have been removed and the effect of aftershocks on the variation of b value has been quantified. The study area was divided into 2,025 grid points (1°91°) and the spatial variation of the seismicity across the region have been analyzed considering all the events within 300 km radius from each grid point. A significant decrease in seismic b value was seen when declustered catalog was used which illustrates that a larger proportion of dependent events in the earthquake catalog are related to lower magnitude events. A list of 203,448 earth- quakes (including aftershocks and foreshocks) occurred in the region covering the period from 250 B.C. to 2010 A.D. with all available details is uploaded in the website http://www.civil.iisc.ernet.in/*sreevals/resource.htm.

A performance-based framework for assessing liquefaction potential based on CPT data

This article describes a new performance-based approach for evaluating the return period of seismic soil liquefaction based on standard penetration test (SPT) and cone penetration test (CPT) data. The conventional liquefaction evaluation methods consider a single acceleration level and magnitude and these approaches fail to take into account the uncertainty in earthquake loading. The seismic hazard analysis based on the probabilistic method clearly shows that a particular acceleration value is being contributed by different magnitudes with varying probability. In the new method presented in this article, the entire range of ground shaking and the entire range of earthquake magnitude are considered and the liquefaction return period is evaluated based on the SPT and CPT data. This article explains the performance-based methodology for the liquefaction analysis – starting from probabilistic seismic hazard analysis (PSHA) for the evaluation of seismic hazard and the performance-based method to evaluate the liquefaction return period. A case study has been done for Bangalore, India, based on SPT data and converted CPT values. The comparison of results obtained from both the methods have been presented. In an area of 220 km2 in Bangalore city, the site class was assessed based on large number of borehole data and 58 Multi-channel analysis of surface wave survey. Using the site class and peak acceleration at rock depth from PSHA, the peak ground acceleration at the ground surface was estimated using probabilistic approach. The liquefaction analysis was done based on 450 borehole data obtained in the study area. The results of CPT match well with the results obtained from similar analysis with SPT data.

Slow Slip Acceleration beneath Andaman Islands Triggered by the 11 April 2012 Indian Ocean Earthquakes

The M-w 8.6 and 8.2 strike-slip earthquakes that struck the northeast Indian Ocean on 11 April 2012 resulted in coseismic deformation both at near and distant sites. The slip distribution, deduced using seismic-wave analysis for the orthogonal faults that ruptured during these earthquakes, is sufficient to predict the coseismic displacements at the Global Positioning System (GPS) sites, such as NTUS, PALK, and CUSV, but fall short at four continuous sites in the Andaman Islands region. Slip modeling, for times prior to the events, suggests that the lower portion of the thrust fault beneath the Andaman Islands has been slipping at least at the rate of 40 cm/yr, in response to the 2004 Sumatra-Andaman coseismic stress change. Modeling of GPS displacements suggests that the en echelon and orthogonal fault ruptures of the 2012 intraplate oceanic earthquakes could have possibly accelerated the ongoing slow slip, along the lower portion of the thrust fault beneath the islands with a month-long slip of 4-10 cm. The misfit to the coseismic GPS displacements along the Andaman Islands could be improved with a better source model, assuming that no local process contributed to this anomaly.