Seismic - Microzonation Guwahati PDF
Seismic - Microzonation Guwahati PDF
On request and initiative of the Hon’ble Chief Minister of Assam, Sri Tarun
Gogoi, the Department of Science & Technology, Government of India
constituted a national level Expert Group vide DST Office Order No. DST/
Exp Group/Guwahati-microzonation/2002 dated August 19, 2002
(Annexure-I) to carry out the Seismic Microzonation of Guwahati Region.
The Expert Group met for the first time at AMTRON in Guwahati on 12th
and 13th of September, 2002 and thereby set the plan of action and the
broad framework for the participating institutions. It was acknowledged
that the Seismic Microzonation of Guwahati Region posed a challenge
hitherto unattempted, when compared with the similar work being carried
out by DST for the city of Jabalpur and New Delhi. In case of Guwahati,
most of the data would have to be collected afresh, meaning it would require
extensive ground survey, including geophysical, geomorphological, land
use, land cover, ground noise response spectra, basement configurations
and landslide hazard mapping, deployment of strong motion
accelerographs, studies on Peak Ground Acceleration etc.
Accordingly, the Expert Group set out the tasks and priorities for the
participating members and institutions. The Groups so far held eight sittings
(two at AMTRON, two at IIT Guwahati & four at IIT Kharagpur Salt Lake
City Campus, Kolkata), and has monitored the progress of work, held critical
review of the outcome, and thus, ensured quality output at every stage.
vii
surveys carried out by GSI during 1986-87 and 1987-88, supplemented
by similar data generated by GSI during 2002-03 field seasons. These
data were compared and contrasted with similar set of data obtained
by DGM from about 30 bore holes in the basin area. Finally a
synthesized basement contour map was prepared by examining both
the sets of data and superimposed the same on the geological map of
the area.
4. Preparation of landslide hazard zonation map of the hilly tract of the
area on 1:25,000 scale by GSI.
5. A seismotectonic map, prepared by GSI of the area of 200km radius
around Guwahati based mostly on the data published by GSI in the
Seismotectonic Atlas of India (2000).
6. Recent satellite imageries of the area studied by AMTRON to verify
and update the change in landforms and land use pattern from that
given in the SOI toposheets surveyed during 1986-87 to achieve an
updated landuse map of the study area.
7. Instrumental site response studies by using ambient noise (Nakamura
ratio technique) were carried out by GSI, IMD and RRL-J to determine
the maximum site amplification factors and corresponding peak
frequency covering almost the entire area at 141 sites.
8. 200 Borehole geotechnical studies by Assam Engineering College.
9. Compilation of results of macroseismic surveys of the past major
earthquakes affecting Guwahati town by D.R. Nandy, Retd. Director,
GSI.
10. All the data digitized and stored at AMTRON, Guwahati.
11. Map analysis (raster & vector) for Microzonation by IIT Kharagpur.
12. Empirical site response synthesis by IIT Kharagpur.
13. Estimation of scenario earthquake magnitude from the earthquake
catalogue of 140 years by IIT Kharagpur.
14. Site amplification, classification, strong motion spectral acceleration
analysis and synthesis by IIT Kharagpur, IIT Guwahati, Assam
Engineering College, Geological Survey of India, IMD and RRL Jorhat.
15. Site classification map on GIS platform using shear wave velocity (VS30),
viii
site response, predominant frequency and factor of safety by IIT
Kharagpur.
16. Demography and Preliminary Seismic Population Risk
Assessment by Jorhat Engineering College and IIT Kharagpur.
17. Integration of all the thematic layers and preparation of Seismic
Microzonation map on GIS platform in the scale of 1:25,000 by
IIT Kharagpur.
18. Integration of Seismic Microzonation and demographic distribution
maps on GIS platform for the preparation of Preliminary Seismic
Population Risk map in the scale of 1:25,000 by IIT Kharagpur.
19. Data repository, web based GIS and hosting on a special website
by AMTRON.
Due to untiring effort of the Group and initiatives of DST in sorting out
issues arising out of inter institutional coordination; the basic groundwork
has been completed, collated and compiled for Guwahati Region. Finally
following maps have been prepared and included in the Atlas
ix
Map XIII Ambient noise survey location map (Predominant
frequency observation sites)
Map XIV Boreholes, Ambient noise survey and Strong motion
station map
Map XV Shear wave velocity(Vs30) contour map of Guwahati
Region
Map XVI Shear wave velocity(Vs30) zonation map of Guwahati
Region
Map XVII Factor of Safety contour map of Guwahati Region
Map XVIII Factor of Safety zonation map of Guwahati Region
Map XIX Bulk density contour map of Guwahati Region
Map XX Predominant frequency contour map of Guwahati
Region
Map XXI Predominant frequency distribution map of Guwahati
Region
Map XXII Site response contour map of Guwahati Region
Map XXIII Site response distribution map of Guwahati Region
Map XXIV Site classification of Guwahati Region based on
integration of Vs30 , predominant frequency, site
response and factor of safety
Map XXV Peak Ground Acceleration (PGA) map of Guwahati
Region. PGA is estimated through Spectral Strong
Motion Synthesis by Brune ω2-circular Crack source
model for an SEM of Mw 8.7 nucleating from the 1897
Great Shillong Earthquake of Mw 8.7
Map XXVI Peak Ground Acceleration (PGA) map of Guwahati
Region. PGA is estimated through Spectral Strong
Motion Synthesis by Synthetic Green’s Function
Simulation for an SEM of Mw 8.7 nucleating from the
1897 Great Shillong Earthquake of Mw 8.7
Map XXVII Seismic Microzonation map of Guwahati Region
(PGA from Green’s Function Simulation)
x
Map XXVIII Seismic Microzonation map of Guwahati Region with
PGA contours overlaid
Map XXIX Demographic distribution map of Guwahati Region
Map XXX Preliminary Seismic Population Risk map of Guwahati
Region
This report gives indepth analysis of the data collated and compiled for the
above mentioned themes. The data has been given in the Data volume, 27
Map Plates and 3 Section profiles.
The Site response was computed empirically at all 141 predominant frequency
observation sites and as well as at five strong motion stations using waveform
data. Site response contour map was drawn taking all the sites into
consideration. This map matches very well with the surface geology and
shear wave velocity map (Vs30). Earthquake catalogue derived from ISC and
USGS from 1866 - 2006 were used for computing the Scenario Earthquake
Magnitude (SEM) for Guwahati region. It is to be noted that the historical
earthquakes 1897 Shillong, 1918 Srimangal, and 1950 Assam constrained
this catalogued estimation. Peak Ground Acceleration (PGA) is estimated
2
through Spectral Strong Motion Synthesis by both Brune ω -circular Crack
source model and impulsive source function using Empirical Green’s Function
Simulation for an SEM of Mw 8.7 nucleating from hypocenter of the 1897
Great Shillong Earthquake. PGA is contoured for the region with the SEM at
the focus of the 1897 Shillong Earthquake to create a seismic hazard scenario
in the region. It is found that PGA varies monotonically between 0.15g - 0.93g
(1g = 1000gal) in the region. Through the overlaying and subsequent
integration of Geology and Geomorphology (GG), Basement (B), Landslide
Zones (L) and Seismological themes such as Shear Wave Velocity (Vs30),
Site Response (SR), Peak Ground Acceleration (PGA) and Predominant
Frequency (PF) on GIS platform, the final seismic hazard Microzonation map
is prepared.
Five zones are mapped where the average PSHI index is 0.55, 0.44, 0.35,
0.26 and 0.16. We termed these zones as very high, high, moderate, low
and very low hazard regions. An integration of this microzonation vector layer
with the demographic distribution of the urban center yielded a preliminary
seismic population risk map of the region.
Prof. Sankar Kumar Nath
IIT Kharagpur
xi
CHAPTER 1
1.1 INTRODUCTION
The northeast India and its adjoining regions are characterized by high seismic
activity. This region encompasses the northern part of the Assam-Arakan geological
province and includes the eastern Himalayas, the Indo-Myanmar arc, the Mishmi
Massif, the Shillong Plateau and adjoining parts, the Tripura folded belt, the Assam
intermountain depression, and the northern part of the Bengal basin. The earlier
seismicity studies of the region were mostly those of investigations of the major
earthquakes originating from the region. One of the first scientific study could be
associated with the classic work of Oldham (Oldham,1899) on the 1897 Great Assam
Earthquake. Pioneering works on seismicity studies of the region also include that
of Bullore, 1904; Pendse, 1948; Tandon, 1954, and Banerjee, 1957. Moreover several
investigators have also studied the seismic phenomena in the region (Dutta, 1964;
Dutta, 1967; Santo, 1969; Than, 1975; Tandon and Srivastava, 1975; Gyi, 1973;
Verma et al., 1976; Khattri and Wyss, 1978). Dutta (1964) studied the annual
frequency relationship. In another paper, Dutta (1967) introduced a further division
of the region into four zones, which coincide with the main geotectonic units of the
region comprising the eastern Himalayas, the Shillong Plateau, the Mishmi Massif
and the Burmese arc. Later, Saikia et al. (1977) introduced another zone
encompassing the Upper Assam depression, parts of the Surma Valley and the
Tripura folded belt. Dutta and Saikia (1976) examined the factors responsible for
structural instability of the fold resistant Assam wedge and the pattern of seismicity,
and attributed the cause to the transformation force systems originating from the
two adjoining arcuate systems - the Himalayan and the Indo-Myanmar. Santo (1969)
related the seismicity to the interaction between India and the Asian continent during
2 Seismic Microzonation Atlas
Though most part of the northeast India is prone to earthquakes, the intensity is
quite different, ranging from 5 to 8 and above, in Richter scale. Earthquakes of low
magnitude of less than 5 on Richter scale are scattered all over the region. However,
most of the earthquakes having a magnitude of more than 5 have been observed in
northern most part of Arunachal Pradesh, some part of lower Brahmaputra Valley,
Central Manipur and West Tripura. Earthquakes ranging between 5 and 6 have
been experienced mainly in northern rim of Assam, West Tripura, Southern Manipur
and eastern rim of Mizoram. Most of the earthquakes, ranging from 6 to 7 magnitudes
of Richter scale, were common in the Brahmaputra Valley, northern Manipur, eastern
Nagaland, and east-central Mizoram. Maximum concentrations of severe
earthquakes, ranging from 7 to 8 magnitudes, have been observed along eastern
Arunachal Pradesh, Central and lower Brahmaputra Valley, and South Eastern
Manipur. Earthquakes having highest magnitude of above 8.0 have been observed
only in the eastern most part, beyond the borders of Arunachal Pradesh and east
Khasi Hills.
The whole of northeast India falls in zone V of the seismic hazard zonation map
(BIS, 2002) of India as shown in Figure 1.1, the highest vulnerable zone in the
country. However, it has been observed during the recent past that earthquake
shaking affects nearby areas differently. Severity of shaking is closely related to
local site conditions.
Chapter 1 3
The Guwahati city being located almost at the center of the region, acts as the
gateway and transit point for communication and transportation for the seven sister
states of northeastern region of India. Shifting of the capital of Assam from Shillong
to Guwahati in 1972 has increased its importance manifold. People from all over
Assam and from the neighboring states have been migrating to Guwahati for job,
business and education. This has resulted in very fast and unplanned growth of the
city which, otherwise, geographically has very limited space in between the relatively
steep sided granite hillocks as seen in Figure 1.2. In early 70s only a few multistoried
buildings existed in the city except in the Downtown area. Most of the houses were
kutchcha-pucca Assam type with corrugated (G.I.) sheet or thatched roofs. But
during the recent past many tall buildings have come up which are juxtaposed to
each other in the main commercial and business hubs and markets, perhaps due
4 Seismic Microzonation Atlas
to acute shortage of space. Many natural water bodies have been filled up for
construction of houses making them more vulnerable to earthquake hazard. In
some cases steep slope of the hillocks has been occupied making life and property
vulnerable to landslides, especially during earthquake and heavy rains. During the
heavy rains of 6-7 October, 2004, 17 people were killed by landslides right within
the city from places like Dhiren Para, Jyotinagar, Chandmari, Jorabat, 9-mile and
Lalganesh.
Figure 1.1 Seismic Zonation Map of India (BIS, 2002) with an inset showing the area
occupied by different seismic hazard zones.
Chapter 1 5
Table 1.2 Decadal Growth of Population of Guwahati City since 1901 to 2001
(Under Guwahati Municipal Corporation Area)
Source: District Census Handbook (1991) Part XII-A & B, Census of India, Assam
Keeping in view the importance of Guwahati city which is located in seismic hazard
zone V, and on request of Sri Tarun Gogoi, Hon’ble Chief Minister of Assam in
2002, the Department of Science and Technology (DST), Government of India
took the initiative for the preparation of Seismic Microzonation map of Guwahati
Region of about 600 sq km bounded by latitudes 26005’ and 26015’, and longitudes
91055’ and 91055’. The Government of Assam designated Assam Electronics
Development Corporation Ltd. (AMTRON), a Govt. of Assam undertaking and Nodal
IT agency of the state, to coordinate on behalf of the State Government with DST
on matters concerning the Microzonation work of the city. An Expert Group was
promptly constituted by DST for the preparation of comprehensive proposal on
Microzonation of Guwahati region vide DST’s Office Memorandum dated
19.08.2002. From the deliberation in the meetings the Expert Group identified
different institutions viz. Geological Survey of India (GSI), India Meteorological
Department (IMD), Survey of India (SOI), Regional Research Laboratory, Jorhat
(RRL-J), IIT Kharagpur (IITKGP), IIT Guwahati (IITG), IIT Roorkee (IITR), Assam
6 Seismic Microzonation Atlas
Number of items of work relevant to the project was identified including collation of
already available data and were assigned to different organizations. Of them, GSI
took the major load of work. AMTRON acted as the nodal agency to compile, digitize
and store all the multidisciplinary data sets supplied/generated by the participating
organizations. Under the main project the DST sanctioned a few sub-projects to
IITKGP, IITG, AEC, JEC and AMTRON.
A regional seismicity map for northeast India (Nandy, 2001) and the adjoining region,
covering the area between 88°E and 98°E longitudes and 20°N and 31°N latitudes
is given in Figure 1.3. Most of the events falling in the tectonic domains of the
eastern Himalaya, Mishmi block, Assam shelf, Meghalaya Plateau and Mikir Hills,
Surma and Bengal basins have shallow focal depth i.e. <70km except for a few
events in the N-W trending wedge shaped block lying between the Kopilli and Bomdila
faults (Figures 1.3 and 1.4).The earthquake events in these tectonic domains occur
in diffused pattern having post-collisional intracratonic characteristics. On the other
hand, most of the earthquake events falling in the Indo-Myanmar (Burmese) tectonic
domain have focal depths varying from 70-200km where seismicity is more intense
and defines the westerly convex broadly N-S subduction zone of the Indian plate.
Seismicity in the Eastern Himalaya, Tethyan Himalaya and over the Yarlung Zangbo
or Tsangpo suture zone is relatively sparse. Events mostly locate in-between the
MBT and MCT and are more concentrated in the areas traversed by the transverse
faults/lineaments running across the Himalayan trend. There is a narrow N-S cluster
of events in the Yarlung-Zangbo suture zone area falling in the line of 90°E meridian,
which are associated with the Yadong-Gulu garben. Most of the events are of shallow
focus (0-40km depth) in these domains with a few having a depth range of 41 to
Chapter 1 7
70km, and a very few having depth of focus varying from 71 to 150km and occur in
the zone where E-W thrust sheets of the Himalaya take NE trend and are intersected
by cross structures. Magnitudes of the earthquakes were from Mw 4.0 to Mw 5.9
except one event of 1947 having magnitude Mw=7.9 with its epicenter at 28°30´,
94°00´. So far no great earthquakes have been reported from these tectonic domains
(Nandy, 2001).
The Mishmi block and the transverse mountain range are seismically more active
than the adjoining Eastern Himalaya. This tectonic domain generated the 1950
Great Assam earthquake (Mw=8.7), perhaps due to strike-slip movement along the
Po Chu fault. This tectonic domain is traversed by many NW-SE thrusts and faults
(Nandy, 2001).
Figure 1.3 Seismicity map of Northeast India and adjoining region (after Nandy, 2001)
Chapter 1 9
Two great and three major earthquakes occurred around Guwahati during the period
of last 137 years. They are
Of these the 1897 event had severe damaging effects in and around Guwahati city.
During 1869 earthquake severe shock was felt at Guwahati (Oldham, 1882). Almost
all the brick buildings suffered damage. East and west wall of the jail developed
horizontal cracks and were tilted; wicket damaged with slight cracks in the arches.
One side of the roof of graveyard porch slipped down. The native infantry hospital
Chapter 1 11
had developed bad cracks either due to faulty construction or settlement of the
wall. The spire of the church got badly cracked. None of the native house in the
bazaar area suffered any damage and no loss of life was reported.
During 1897-Great Assam Earthquake (Oldham, 1899) one of the gate pillars of
Loki Rani’s house fell to N600E, the other gate pillar was removed; a pair of gate
pillars in the compound east of Telegraph Club was severely broken and overthrown;
in the cemetery a marble cross was broken across the socket; the base of stone
pillar of Robert Beecher monument had been projected to a distance of 84cm
horizontally and 1.15m vertically; the clock of the telegraph office stopped at 5.15
pm, local time, the pendulum was broken up the shock; one of the 3m high pillars of
the Deputy Commissioner’s Bungalow had broken through 84cm from the ground
and got twisted by 6 degrees, other pillar had broken at the same level but was not
twisted, low compound wall had fallen; out of the 4 gate pillars made of masonry at
two entrance of the old Dak Bungalow three were broken and twisted and the fourth
one was only broken and fallen; the Commissioner’s and Deputy Commissioner’s
offices were wrecked with the collapse of the walls; capping of the small gate pillar
was shoot off to a distance of 1.30m from the center of the pillar; the obelisk of
almost all the tombs in the old cemetery situated over alluvial ground south of the
station had broken and fallen, one of them was twisted by 480, one recently built
(1875) marble tomb was badly damaged whereas the new cemetery situated on
the flanks of the low hill of gneissic rock east of the station, had suffered little or no
damage; newly constructed brick building of the Railway Station was greatly damaged
with fallen center gable and broken walls and diagonal cracks; E-W wall of the old
temple at the streamer ghat had fallen completely.
The bridge made of three girders over a small stream near the western end of the
bazaar along the G.T. Road had shortened by 45cm due to fissuring of the banks on
both side of the river-abutment having been carried forward, one of the piers was
tilted. Ground in the country surrounding Guwahati had suffered fissuring along the
bank of the river and consequent subsidence of surface, and ejection of sands had
filled up rivers and water bodies; Strand Road from Sukleswar ghat to Bhorolumukh
furnished striking illustration of fissures which developed parallel to the river bank
having 60cm to 120cm width. Near Chouki, 8km north of Guwahati, a small landslip
12 Seismic Microzonation Atlas
occurred from a hill on the roadside. Further north, near Chutiapara, a huge rock
fell down the Kohra and Deodual hill.
During 1931 Dhubri Earthquake, both the Revenue and Judicial Record rooms
were badly damaged. The walls of the buildings were separated at the corners and
cracked from the post-plates through door and window lintels. The wall of the Jail
was damaged for a length of 12m from the height of 1.5m to the top. The walls were
cracked horizontally and moved out about 2.50cm. The Commissioner’s Record
room was damaged over two windows.
From the foregoing description it is evident that the alluvial tract of the city is more
vulnerable to damage than the hilly tracts by earthquake shaking excepting the
chance of land/rock slides; and alluvial tracts by the side of the river are highly
vulnerable to liquefaction phenomenon.
i. Topography,
iii. Depth of sedimentary basin and bed rock (basement) topography, basin edge
conditions and lateral inhomogeneities, vertical/inclined layering of sediments
etc.
Amplitude and shape of the seismic waves are modified by surficial soil through
which they are transmitted in at least three ways:
(iii) By causing wave to increase in amplitude in soil that are less tightly packed.
Basin effect generally influences long-period wave energy and cause amplified
motion in high-rise building and multi-span bridges. Scattering alters the short-period
wave energy that affects shorter structures. Impedance contrast affects the amplitude
and phase of both short and long period wave. Detailed analysis of macroseismic
data of past earthquakes, if available, with reference to topographic and geotechnical
maps may lead to qualitative appraisal of most hazardous zone.
Beside considering the above factors for the purpose of seismic microzonation, in
the absence of well spaced strong motion data, site condition may directly be obtained
by direct measurement of shear wave velocity at each site by down hole instrumental
probing up to 30m depth or indirectly through noise survey and microtremor
measurement over the area. It has been observed that the ratio between Fourier
spectra of horizontal and vertical components can be used to identify the resonance
fundamental frequency of soft soil (Nakamura, 1989), but it fails for higher harmonics
and that peak amplitude is somewhat different from amplification measured in the
spectral ratio.
CHAPTER 2
2.1 INTRODUCTION
The northeast Indian region falls into the category of high seismic zone (BIS, 2002).
There have been incidents of 20 large earthquakes (8.0>Mw>7.0), in addition to the
two great earthquakes (Mw 8.7) - the Great Shillong Earthquake of 1897 and the
Great Assam Earthquake of 1950 in the region since 1897 (Kayal et al., 2006). The
region as such, is placed in Zone V, the highest level of seismic hazard potential,
according to the seismic zonation map of India (BIS, 2002). The Global Seismic
Hazard Assessment Programme (GSHAP) also classifies the region in the zone of
highest risk with Peak Ground Acceleration (PGA) values in the order of 0.35-0.4g
(Bhatia et al., 1999). Moreover, in the region, there is a higher level of man-made
constructions and a significantly larger population than at the time of great earth-
quakes implicating more vulnerability to earthquake disasters.
Data sources considered are USGS catalogue, ISC catalogue (ISC, 2006), IRIS
catalogue, and ISET catalogue (Bapat et al., 1983). Data for the study pertains to a
spatial scope of latitudes 21oN to 31oN and longitudes 87oE to 99oE that include an
extra of 1o on each side of the rectangular area as a buffer for the spatial windowing.
Lower limit of magnitude 4.0 is taken since magnitude = 4.0 is considered well
recorded. We look into the aspect of heterogeneity of the catalogue, which is due to
different types of magnitude scales employed. However, magnitude scales: surface-
wave magnitude (Ms), body-wave magnitude (Mb) and Local magnitude (ML) are
related to the largest amplitude that is recorded on a seismogram (Shearer, 1999).
Ms and Mb are generally adopted and calibrated to agree with the ML for the small
events. Mb and ML are found to be saturated for the large/great events such that Ms
or moment magnitude (Mw) is used for the purpose. Mw and ML are approximated as
equal below 5.5. Mw and Ms are identical up to magnitude 7.5 (Chen and Scawthorn,
2003). So, Mw and Ms can be taken as equivalent within their measuring limits. We
can convert Mb into Ms using the empirical relation constructed with record’s entry
of both type of magnitude scales.
A simple method is adopted in order to enrich the database yet avoid duplications.
The data sources are compared with each year’s record. The record of the particular
year with the maximum number of earthquakes is selected to construct the data-
base for the specific year. This is done for each year, till all the years are covered and
the database is completed. Though events for the period from 1764 to 1897 are found
in the ISET and USGS catalogues, there are large gaps in the recorded earthquake
history prior to 1866. Therefore, the period accounted for the present study is from
May 23, 1866 to May 15, 2006. The USGS catalogue refers to the merged NEIC
database of PDE, India and significant earthquakes worldwide for the mentioned
period and the study region. For the period prior to 1964, USGS provides an overall
higher count. However ISET catalogue shows higher count for the years 1880 and
1956, and ISC catalogue has more events in the years 1914, 1931-35, 1938-39 and
1941. From 1964 onwards IRIS catalogue is found to have higher counts for most of
the years. ISET catalogue has more events in the year 1964. Higher counts are
found in the ISC catalogue for the years 1965, 1967-70, 1973-75 and 2004-06. Higher
Chapter 2 17
number of events for the year 1976 is found with the USGS catalogue. After merger,
the catalogue is scanned through to remove any trivial entries.
The database obtained contains a good number of entries in Mb. The entire cata-
logue histogram is given in Figure 2.1(a). To establish an empirical relation between
Ms and Mb, a linear regression is performed using records with both magnitudes
which are found in the ISC and IRIS catalogue as shown in the Figure 2.1(b). The
relation obtained is
The linear regression is compared with that of Richter (1958) given below,
The former relation gives an average weight to the higher and the lower magnitudes
which is rational to the data. The Richter (1958) relation pulls down the lower
magnitudes of Mb to further lower values and elevates the higher value of Mb to
further higher value of Ms.
2
values of stress drop, ó = 50bars and rigidity, ì = 5x1011dyn/cm for earthquakes in
the crust and upper mantle is
Also we have,
where Mo is in Nm.
Figure 2.1 (a) The number of events per year in the merged catalogue from ISET,
IRIS, ISC and USGS for the period May 23, 1866 - May 15, 2006, consisting only of
events with M ³ 4 earthquakes. The difference in the clustering of earthquakes can
be seen in the data for period prior to 1964 and that of 1964 onwards. A difference
owed to the deployment of the World Wide Seismograph Station Network (WWSSN)
in 1964. (b) Linear regression on Ms against Mb of the total 1201 records of events
found with both of the magnitudes in the ISC and IRIS catalogues merged, which
are most likely duplicated (obviously not determinant in this workout), yields an
empirical relation Ms=(0.89 ± 0.049) Mb + (0.18 ± 0.024) as given by thick line. The
dash lines represents upper and lower bounds with 95% confidence. The normal
line plot of Richter (1958) Ms=1.59 Mb -3.97 is given for comparison purpose.
20 Seismic Microzonation Atlas
Owing to the improvement of quality of recordings from 1964 onwards, a sub catalogue
derived from the ISC catalogue for the period 1964-2006 is employed for the
assessment of the b-value distribution. The detail of the analysis is discussed by
Thingbaijam et al. (2007).
The b-value is found to vary both spatially and temporally and is often employed as
one parameter approach for earthquake studies (Gibowicz and Lasocki, 2001; Nuannin
et al., 2005). The different interpretations of b-value variations are discussed in (Mogi,
1962; Scholz, 1968; Schorlemmer et al., 2005; Warren and Latham, 1970; Wyss,
1973; Wesnouski et al., 1983). A low value implies that majority of earthquakes are
of higher magnitude and a high value implies that the majority of earthquakes are of
lower magnitude. The variations of b-values are seen to be inversely related to stress
(Mogi, 1962; Schorlemmer et al., 2005; Wesnouski et al., 1983). Large material
heterogeneities are reported with higher b-values (Scholz, 1968). High b-values are
associated with aftershocks and low b-values with foreshocks (Suyehiro et al., 1964;
Nuannin et al., 2005). The maximum likelihood method for the evaluation of b-values
is widely accepted (Aki, 1965; Bender, 1983; Utsu, 1965). The b-value is estimated
as
log10 (e )
b=
Äm ..................................................................... (2.7)
m Mean - m t - 2
where mMean is the average magnitude, mt is the threshold or the minimum magnitude
and Dm is the magnitude bin size. The standard deviation, db, of the b-value in
equation (2.7) is estimated as
∑ (m - m ) /N (N - 1) ............................................ (2.8)
2
δb = 2.30b2 i Mean
i=1
where N is the total number of events in the sample (Shi and Bolt, 1982).
The temporal seismicity pattern involving the number of events indicates period of
commencement of seismicity buildup from 1964 onwards. The period is also
accounted for improvement in quality of records due to the advent of digital
22 Seismic Microzonation Atlas
instrumentation. Hence, the sub catalogue for the period 1964 - 2006 is preferred
and used for the estimation of spatial distribution of the seismicity parameters.
The estimated b-value ranges from 0.59 to 1.42. The variation of b-values is shown in
Figure 2.3(a). An examination of the spatial distribution reveals relatively higher
b-values on the eastern Himalayan area. The b-values are low in most of the central
Eastern Boundary Thrust (EBT) encompassing the Arakan Yoma range with an
increase towards both the south and the north. However eastern end of the Lohit
thrust and southern most end of the EBT zone also have lower b-values. Southeast
part of the map is observed to have low b-values. The standard deviations of the
b-value, db, is found to vary from as low as 0.033 to as high as 0.165 and are depicted
in Figure 2.3(b). Though comparative quite high b-values are observed in Dhansiri
Kopili fault and Bomdila lineament zone and adjoining region, large standard deviations
are also seen. High deviation is also marked on the northwest end of the Lohit thrust.
However, across the EBT and the upper Himalayan zone, the standard deviations
associated with the b-values are low to moderate.
framework of classification. The tentative precincts are setup with the historical great
earthquakes and thereafter, the seismotectonic settings are segregated into four
zones according to indicative trends observed in the spatial distributions of b-value.
The trends indicate high values on the eastern Himalayan region. While low values
are seen in the thrust zones especially in the central EBT dominated by intermediate
depth earthquakes. However higher values are observed on the north and the south
of the EBT zone. High values are again seen on the regions of Mishmi block and Po
Chu fault. The Shillong plateau region has dominantly moderate values.
(a)
Figure 2.3 cont’d...
24 Seismic Microzonation Atlas
(b)
Figure 2.3 The spatial distributions over the backdrop of seismicity with
Mw ³ 4.0 from the sub catalogue covering a period 1964-2006:
(a) b-value, and (b) standard deviation of the b-value (db)
Figure 2.4 The classified source zones demarcated by the polygons on the
seismotectonic map of Figure 2.2.
Chapter 2 25
Finally the boundaries are demarcated with the consideration of tectonic regimes.
Four source zones are classified as given in Figure 2.4. Zone 1 is Himalayan Source,
Zone 2 is figured on Mishmi block. Zone 3 outlines the Indo Myanmar arc and Zone
4 encompasses the Bhramaputra basin, Kopili fault, Naga thrust, Shillong plateau
and Duaki fault. Zone 1 is delineated in the coordinates 87°13.2’E, 31°N; 93°19.8’E,
31°N; 93°19.8’E, 27°30’N; 92°12’E, 26°48’N; and 87°13.2’E, 26°48’N. The coordinates
demarcating Zone 2 are 93°19.8’E, 30°15’ N; 98°06.6’E, 30°15’N; 98°06.6’E, 27°30’N;
and 93°19.8’E, 27°30’N. Zone 3 is bound within 95°48’E, 27°30’N; 98°06.6’E, 27°3’N;
94°24’E, 21°N; and 91°48’E, 23°06’N and the coordinates of Zone 4 are 88°07.8’E,
26°48’N; 92°12’E, 26°4’2N; 93°19.8’E, 27°30’N; 95°48’E, 27°30’N; and 91°48’E,
23°06’N.
The approach for the estimation of the maximum earthquake, Mmax, is similar to one
employed for the Aegean area by Koravos et al.,(2003). With the idea that there is a
bound on the maximum energy release, the predicted maximum earthquake
magnitudes are generally constrained with finite seismic moment release or tectonic
moment release or both. In the present analysis we employ the seismic moment
release computed from the data catalogue. We looked into the four seismic zones
separately and in addition, the whole region as singular source regime using five data
catalogs, one for each zone and the main catalogue for the entire region. The data
catalogs are extracted from the main data catalogue according to the coordinates
defined by the source zones. It is noted that a potential source of uncertainty is the
number of year gaps in the preceding years of the catalogue duration. The gaps are
made less promiscuous by considering that there is neither significant number of
events, nor events of significantly big magnitude during the gaps. The seismic moment
release rate for each catalogue, is estimated by summing up the moment release for
each earthquake in the data and dividing with the catalogue duration.
26 Seismic Microzonation Atlas
The results for all four source zone are shown in Table 2.1. The source zone 4
encompassing Shillong plateau which is the potential source of earthquake hazard
to the Guwahati region, is associated with Mmax of Mw 9.1 ± 0.14 implicating the
historical maximum of Mw 8.7 as a more realistic choice for the SEM.
Figure 2.5 Incremental frequency magnitude distribution plots for Source Zone 4. The
thick lines in the plot indicate the best fit line. The dashed lines illustrate the 95 per cent
intervals on this fit. The error bars given for each magnitude are estimated as F/ÖN where
F is the frequency and N is the number of observations used to compute F (Koravos et al.,
2003).
Chapter 2 27
Table 2.1: The estimated maximum earthquake is compared with historical maximum
earthquake. The completeness thresholds mc estimated from the frequency magnitude
distribution which conforms to of the fitting model at the lower magnitude bound are also
listed. The thresholds are employed for the determination of annual number of events and
annual seismic moment release for each source zone. # indicates the entire region.
3.1 INTRODUCTION
It has long been known that each soil type responds differently when subjected to
ground motion from earthquakes. Usually the younger softer soil amplifies the ground
motion relative to older more competent soils or bedrock. The potentially severe
consequences of this phenomenon were recently demonstrated in the damage
patterns of the 1985 Michoacan, Mexico earthquake (Singh et al., 1988), the 1988
Armenian earthquake (Borcherdt et al., 1989), the 1989 Loma Prieta earthquake
(Hough et al., 1990; Borcherdt and Glassmoyer, 1992) and the Northridge earthquake
in Los Angeles, California (EERI, 1994). Numerous other studies have also
demonstrated the ability of surface geologic conditions to alter seismic motions
(Borcherdt, 1970; King and Tucker, 1984; Aki, 1988; Field et al., 1992).
There are many factors that influence the way a site will respond to earthquake
ground motion (Aki, 1988; Aki and Irikura, 1991; Bard, 1995). These include: (i) the
source location, (ii) the prevalence of energy focused or scattered from lateral
heterogeneity, and (iii) the degree to which sediments behave nonlinearly, which
causes the response to depend on the level of input motion. However, the site
effects in the assessment of seismic hazard follow a simple approach wherein, for
the potential sources of earthquake ground motion in a region the unique behavior
of one site in relation to others is calculated.
30 Seismic Microzonation Atlas
After the occurrence of large destructive earthquakes during the last 20 years, such
as the Mexico 1985 (Bard and Chávez-García, 1993), the Armenia 1988 (Borcherdt
et al., 1989), the Loma Prieta 1989 (Hough et al., 1990), the Northridge 1994 (EERI,
1994), and the Kobe 1995 (EERI, 1995), both the seismologists and earthquake
engineers have focused their attention on the importance of local site response on
seismic ground motion. The observed unequal distribution of damage from these
earthquakes prompted a series of studies on scenarios and parameters concerning
local geological conditions responsible for the differentiation of the seismic ground
motion. Indeed, many factors affect the seismic motion near the recording site:
seismic–wave velocities and geometry of the stratigraphy (the thickness and position
of the discontinuities of the geological formations), topography (Boore, 1972, 1973;
Bard, 1995), and intensity of excitation (weak or strong) that can induce nonlinear
phenomena in correlation with local geology (Aki, 1993; Field et al., 1997).
The term ‘site effect’ means different things to different audiences. We take the term
to represent local ground response, basin effects, and the influence of surface
topography. The definition of surface topography is obvious. ‘Local ground response’
refers to the influence of relatively shallow geologic formations on (nearly) vertically
propagating body waves. These effects are ideally modeled using the full soil profile,
but for deep alluvial basins the modeling domain generally does not extend beyond
depths of 100-200m.
Site effects play a very important role in characterizing seismic ground motions
because they may strongly amplify (or de-amplify) seismic motions at the last moment
just before reaching the surface of the ground or the basement of man-made
structures.
The greatest challenge in estimating site response from earthquake data is removing
the source and path effects. Borcherdt (1970) introduced a simple procedure to divide
the spectrum observed at the site in question by the same observed at a nearby
reference site, preferably on competent bedrock. The resulting spectral ratio
constitutes an estimate of the site response if the reference site has a negligible site
response. Andrews (1986) introduced a generalized inverse technique to compute
site response by solving data of a number of recorded events for all source/path
effects and site effects simultaneously. These techniques for computing site response
Chapter 3 31
depend on the availability of an adequate reference site (on competent bedrock) with
negligible site response. Since such a site may not always be available, it is desirable
to develop alternative methods that do not depend on a reference site. Boatwright et
al. (1991b) suggested a generalized inversion scheme where shearwave spectra are
represented with a parameterized source- and path-effect model and a frequency-
dependent site response term for each station. Another non-reference-site-dependent
technique involves dividing the horizontal-component shear wave spectra at each
site by the vertical-component spectrum observed at that site (Lermo and Chávez-
García, 1993). This method analogous to the so-called receiver function technique
(Langston, 1979) used to study the upper mantle and crust from tele-seismic records,
assumes that the local site conditions are relatively transparent to the motion that
appears on the vertical component. Nakamura (1989) introduced another technique
for analyzing the ambient seismic noise. He hypothesized that site response could
be estimated by dividing horizontal-component noise spectra by vertical-component
noise spectra. Several studies have since shown that Nakamura’s procedure can be
successful in identifying the fundamental resonant frequency of sedimentary deposits
(Omachi et al., 1991; Lermo and Chávez-García, 1992; Field and Jacob, 1993; Field
et al., 1995).
The receiver function analysis exploits the fact that tele-seismic P-waves that are
incident upon the crustal section below a station produce P to S conversions at
crustal boundaries as well as multiple reverberations in the shallow layers. By
deconvolving the vertical-component signal from the horizontal-components, the
obscuring effects of source function and instrument response can be removed,
leaving a signal composed of primarily S-wave conversions below the station. The
deconvolved horizontal component called receiver function trace is a best
representative of the site response as the local site conditions are relatively
transparent to the motion that appears on the vertical component.
To evaluate site response several techniques have been utilized and compared in
recent studies (Field and Jacob, 1995). Two of the proposed methods - the standard
spectral ratio (SSR) (Borcherdt, 1970) and the receiver function technique (HVSR)
32 Seismic Microzonation Atlas
(Langston, 1979; Lermo and Chávez-García, 1993; Nath et al., 2000, 2002a, 2002b)
are based on a spectral ratio scheme. In both these techniques, the source and path
contributions are removed from the seismic recordings by means of a deconvolution
operation using a function free of site effects.
The seismograms of the selected events were first corrected for the system response.
Next the S-wave packets recorded by the seismographs were windowed with a window
width containing the maximum amplitude. The window length was selected following
the results of Seekins et al. (1996). A Hanning taper is applied to the time windowed
data and then butterworth bandpass filtered before the amplitude spectra were
computed.
Let the S-wave spectral amplitude and that of the background noise be O(rij, fk) and
B(rij, fk) respectively at the hypocentral distance rij. Then the signal amplitude spectrum
at the frequency fk can be expressed as,
The corrected spectra are smoothed in order to reduce the data variance using a five-
point smoothing window and a spline interpolator at 0.1Hz interval.
The receiver function HVSRij(fk) can be computed at each j site for the ith event at the
central frequency fk from the root mean square average of the amplitude spectra as,
1
absHij ( fk ) NS
2
+absHij ( fk ) EW
2
Finally, the event average receiver function HVSRij j ( fk ) (Field and Jacob, 1995) is
ave
computed at each j site for the kth frequency to consider the contribution of all the
seismic events recorded at that station.
ρβ
Amp ( f ) = exp ( −π ko f ) .......................................................
ρs ( f ) β s ( f ) (3.4)
where b and r are the shearwave velocity and density of the crust beneath the site.
The effective velocity bs(f), effective density rs(f ), and site damping factor Ko are calculated
from the expressions given by Boore (2003).
In this study two hundred boreholes were used for the site response calculation. All
the geotechnical parameters as well as estimated site response at 141 noise survey
locations is given in Table 3.1.
34 Seismic Microzonation Atlas
The tectonic setting of the Himalaya and its relation with the seismicity of the region
has drawn world attention from earth scientists. The Guwahati region, which is a part
of the northeast India, is placed in Zone V (PGA > 0.4g), the highest region of the
seismic zonation map of India (BIS, 2002, Figure 1.1). Northeast India lies at the
junction of the Himalayan Arc to the north and the Burmese Arc to the east. The high
seismicity in the region is attributed to the collision tectonics between the Indian
plate and the Eurasian plate to the north and Indo-Myanmar range to the east. Global
Positioning System (GPS) measurements show that India and southern Tibet
converge at 20±3mm/year (Bilham et al., 2001). Bilham et al. (2001) divided the
Himalaya into 10 imaginary sections, each around 220km in length. At the observed
convergence rate of 20mm/year, at least 6 of these 10 regions have an accumulated
slip potential of 4m. This is equivalent to the slip believed to be associated with the
1934 Nepal-Bihar earthquake that killed 10,700 people. Site amplification is one of
the important factors which contribute strongly to the hazard. Strong motion data is
record of short period ground motion which strongly accumulate the soil effect due to
which amplification takes place. Amplification of shear-wave at free surface is easily
estimated from strong motion data by deconvolving the vertical component of
accelerogram from radial and transverse component.
Qs = (134.15 ± 29.92 ) f (
0.52 ± 0.09 ) ......................................................... (3.5)
The QS obtained in the present study represents the overall attenuation of the seismic
wave energy, which includes the direct S-wave, early coda, and possibly Lg phase of
the recorded data from events with focal depths more than 10km, except for a few
shallow foci earthquakes. Since in our case there is no deeper event with a focal
depth beyond 35km, we restrict our observations to only one seismic wave energy
attenuation relation with a trade off between the source and the attenuation factors
built in the convolution model.
36 Seismic Microzonation Atlas
Figure 3.1 The plots of QS verses frequency. The dashed lines represent the
zone of scatter of the data points with respect to the power law relation between
QS and frequency with + one standard deviation in the constant and the exponent
Station site amplification has been computed from 2 strong motion events with signal-
to-background noise ratio greater than 3. The site response has been calculated by
using HVSR technique for all the events at each site for different source azimuths.
Figure 3.2 (d) and (e) represents the site response due to events recorded at 302.34°N
and 335.17°N azimuth at the station .
(a) (b)
Chapter 3 37
(c) (d)
(e)
Figure 3.2 (a) Radial component, (b) Transverse component, (c) Vertical
component of accelerogram recorded at an azimuth 335.17oN, (d) RMS site
response by HVSR for source azimuth 302.34°N and (e) RMS site response by
HVSR for source azimuth 335.17°N
The US Army Waterways Experiment Station (WES) has been using the computer
program SHAKE to calculate site response for level-ground soil sites for more than
15 years, including use on a number of USACE projects. WES has continually
made adaptations to SHAKE as the use for each new project required. The original
version for use on a personal computer was obtained from the University of California
at Berkeley (UCB) around 1985. This program at WES is now called WESHAKE91
38 Seismic Microzonation Atlas
to reflect the numerous changes that have been made to keep pace with state-of-the-
art technology, to provide for needs of USACE users, and to provide a user-friendly
interface. These adaptations facilitate transfer technology to, and wide-spread use
among USACE personnel.
past studies have shown that reasonable results are obtained for a much broader
spectrum of in situ conditions. The equivalent-linear constitutive model, assumption
(d), is described later in this section. The last assumption (e) narrows the focus to a
simple class of problems, but, is a common assumption for this type of problem. It is
important to realize that the formulation of SHAKE for wave propagation is based on
a total stress analysis. The materials are considered to be continua and pore water
pressures are non-existent. The calculation of shear modulus using values of K2
does involve the determination of mean effective stress using the depth of the water
table and the unit weight of water.
Resonance frequency of each soil type also differs depending upon the physical
property, and depth to bedrock. Quantification of this amplification of ground motion
and determination of natural resonance frequency is the main objective of site
response study. This frequency dependent amplification forms an important factor
for seismic hazard analysis and microzonation studies. The site response parameters
are also used to distinguish regions where seismic hazard is high due to amplification
from the surface geology and match of natural frequency of the soil with the
construction. The mapping of the soil behavior before a seismic wave - field
propagates also provides an overview of the possible damage to individual structure
or a set of buildings.
The total area of the Guwahati city is about 600 sq km lying between longitude 910
30¢-910 50¢E and latitude 260 05¢- 260 12¢N, covering almost the entire urban part of
the city. The Guwahati city is covered with recent alluvium with some Archean hillocks
at places. The area falls in the Lower Assam valley bounded by the Eastern Himalaya
to the north and Shillong Plateau to the south. The Lower Assam valley consists of
crystalline rocks that are covered by gently dipping Tertiary and younger sediments.
The sediment thickness varies from ten to few hundred meters in the study area.
Chapter 3 41
Several methods for soil characterization are described in literature and used for site
response study of several cities (Bard, 1997, 2000; Kudo, 1995). Broadly these
methods can be categorized as follows:
i) Numerical,
iii) Experimental.
This method considers that spectral amplification of a surface layer could be obtained
by evaluating the horizontal to vertical spectral ratio of the microtremors recorded at
the site. The main challenge to determine site amplification characteristics from
microtremors is to remove source and path effect. Nakamura (1989, 1996) proposed
an approximate procedure for removing source effects from microtremor records based
on a modification of the conventional transfer function of the site. The following
hypothesis are adopted:
(ii) The effect of Rayleigh waves remarkably appears in the vertical tremor.
Use of microtremors appears to be very convenient because of its cost and time
effective nature. The studies when conducted in conjunction with other
complementary techniques, such as shear wave velocity measurements (Aki’s
technique), receiver-function type technique based on weak motion data and
especially in reference site dependent mode, would be an effective tool for site
characterization ascertaining peak/resonance frequency and amplification.
Hence, the site response studies at Guwahati were conducted resorting to Nakamura
type studies based on microtremor data and geological attributes. 141 sites (Annexure
II) were selected for response studies to characterize the ground condition defining
their response parameters.
The response studies were conducted in long term array and rapid modes for
ascertaining response parameters for different ground conditions. The data
acquisition was made following the recommendations of Bard (2000) and Mucciarelli
(1998) using velocity sensor of natural period 1sec.
To avoid errors due to day and night variations in noise pattern at Guwahati, enough
data of about 24 hours to several days have been collected at each representative
site covering day, night and local variation and each hour data is used for spectrum
analysis. At most of the stations arrangements have been made to ensure perfect
coupling between sensor and the natural soil of the site.
44 Seismic Microzonation Atlas
For processing of acquired data, the following steps have been adopted
a) Continuous waveform data collected for several days have been divided in one
hour waveform data files and thus 80 to 100 waveform data files have been
created for each representative site.
b) Continuous waveform data collected for one hour have been divided in 8 to
10min waveform data files and thus about 6 to 8 waveform data files have been
created for each representative site.
c) For the data collected in trigger mode for 24 to 48 hours, each waveform file
created through trigger algorithm has been examined and approximately 80 to
100 waveform data files have been created for each site.
d) It is seen that H/V ratio, gets affected by “heavy” traffic; similarly, walking near
the sensor create spurious transients. To avoid spurious inconsistent data, a
portion of uniform wave train has been selected for the analysis. Each waveform
data file generated has been plotted and a portion of smooth common wave
train of about 120 seconds has been separated from these and new waveform
data file has been created for selected portion of data.
e) Appropriate data portion of each file then has been used for spectrum analysis
as per the criteria of time window length selection described below.
In such a study, selection of time window for spectral analysis is an important aspect.
The application of the spectral ratio technique to micro tremor records present a
basic problem; it is very difficult to identify a common wave train for a longer duration
for the station involved, particularly in urban areas due to the movement of heavy
transport, adjacent to the observational sites, which cannot be avoided.
Chapter 3 45
Generally accepted rule of thumb in site response studies is that the window length
should have at least 10 cycles for the lowest frequency analyzed (Bard, 1997). The
frequency of interest can be ascertained on the basis of frequency of structures in
the area of interest, which can be determined using the relationship between the
height of the building and its fundamental period of vibration T=Number of storey/10.
In the Jabalpur area frequency of built environment ranging from 0.33 sec (3Hz) to
0.16sec (6Hz), therefore, about 30 sec time window length is found to be appropriate.
The SPEC program computes H/V of individual component. Then, root mean square
is calculated for all the individual plots to get the final H/V ratio. To avoid spurious
peak links with sharp troughs on spectrum and for clear identification of Peak frequency,
spectra need to be appropriately smoothened. In the present study, smoothing has
been performed 300 times.
Following the Ibs-von Seht and Wohlenberg (1999) relationship which is given by
equation (3.6), predominant frequency is empirically derived from the basement depth
(Table 3.1). Good match has been observed between predominant frequencies derived
from empirical relation and those estimated from ambient noise survey.
46 Seismic Microzonation Atlas
m = 96 ( f −1.388 ) .........................................................................(3.6)
‘m’ is the basement depth in metre and ‘f’ in Hz is the predominant frequency of the
site under consideration.
3.6 REMARKS
Due to highly variable geological features of the Guwahati urban agglomeration, the
peak resonance frequency showed high variation from site to site. The frequency
for most part of the urban area lies between 0.5Hz to 5.0Hz. However, there are
isolated sites, where peak frequency is more than 5 Hz and lies between 5 to 10Hz.
Figure 3.5 illustrates the corresponding peak frequency contour map. It shows that
Geological domain of Active Flood plain and Levee have the peak frequency >0.5Hz,
Digaru Surface have peak frequency between 0.5-0.8Hz, Bordang Surface have
peak frequency between 0.8-2.0Hz; Sonapur Surface between 2.0–4.0Hz; Rocky
areas of Pediment have peak frequency between 4-6Hz and Denuded Hill have peak
frequency greater than 6.0Hz.
Most important parameter for site response study is the shear wave velocity of each
site, determination of which is very expensive and time consuming. Some empirical
relations between SPT-N values and S-wave velocity have been established that can
be used to obtain ideas about the expected shear wave velocity of the surfacial soft
soil. Fumal and Tinslay (1985), after measuring N-values and corresponding shear
wave velocities of Holocene-Pleistocene sediments of the Los Angeles region in
California, USA, have developed three separate correlations between N-values and
shear wave velocities for three textural group of sediments, viz. clay and silty clay,
silt loam and sandy clay, and gravelly sand respectively. The three correlations are
as follows:
Vs = 4.3N + 218 for silt loam and sandy clay ................................... (3.8)
Tonouchi, et al. (1983) has given another empirical relation between N-values and
S-wave velocities as under –
SPT data of 200 borehole sites (Figure 3.6) were collected having depth range varying
from 6m to 30m. Majority of locations of these sites fall in the Downtown area and
along a zone on either side of the G-S road. Depth-wise N-value data have been
tabulated for each site; then each N-value has been converted in to S-wave velocity
from the Equations 3.10 and 3.7. The values thus obtained for each layer were averaged
to arrive at average Vs30 according to the individual relations as stated. Ultimately the
S-wave velocities obtained through each relationship were again averaged to arrive at
expected Vs30 at the site. Details of geotechnical/lithological parameter used in these
analyses have been given in the Annexure: Annexure III for Location, Depth of Ground
48 Seismic Microzonation Atlas
Water Table, and Date of drilling and Litho-Log of each SPT borehole; Annexure IV for
Physical and Shear Parameters of sediment as obtained from boreholes; Annexure
V for Shear wave velocity data at different depth; Annexure VI for Shear wave velocity
(Vs30) at each borehole; Annexure VII for Soil density at different depth and Annexure
VIII for Factor of safety.
The values of shear wave velocity (Vs30), thus obtained, vary from 200m/s to 360m/s.
High shear wave velocities have been found at a few spots along the G-S road,
northern portion of the Cantonment area and over the rocky ground by the side of
Brahmaputra River in the Ujan Bazar area. While plotting on the base map the
shear wave velocity data have been grouped into 4 classes viz. 200-240, 240-280,
280-320 and >320 meter per second. The first group indicates soft soil and the rest
indicate stiff soil as per the classification of UBC, 1977.
Chapter 3 49
Contour map of bulk density and effective shear wave velocity (Vs30) calculated from
Boore’s (2003) relation is shown in Figures 3.7 and 3.8. Site response for Guwahati
region was calculated at 131 predominant frequency observation sites out of 141.
Rest of the observations were taken on hillocks where borehole information was not
available. Site response calculated from Equation (3.4), along with geotechnical
parameters used in site amplification is given in Table 3.1. Table 3.2 represents a
comparison between the site response calculated from HVSR (strong motion data)
and empirical method at predominant frequencies. Value obtained by both the method
depicts good match. We have used the empirical site response at borehole locations
for our analysis to avoid extrapolation in the contouring because boreholes are well
distributed in the region. Figure 3.9 represents site response contour map of Guwahati
region.
Table 3.1: Site Response and Geotechnical Parameters in Guwahati Region at 141
Noise Survey Locations
Sl. Lat Long Predominant Predominant Effective Effective Q-factor Damping Site
No. (°N) (°E) frequency frequency Shear Wave Density Factor Response
estimated estimated by Velocity (gm/cc) ( K0*10-4)
by noise empirical Vs30 (m/sec)
survey (Hz) relation (Hz)
Sl. Lat Long Predominant Predominant Effective Effective Q-factor Damping Site
No. (°N) (°E) frequency frequency Shear Wave Density Factor Response
estimated estimated by Velocity (gm/cc) ( K0*10-4)
by noise empirical Vs30 (m/sec)
survey (Hz) relation (Hz)
Sl. Lat Long Predominant Predominant Effective Effective Q-factor Damping Site
No. (°N) (°E) frequency frequency Shear Wave Density Factor Response
estimated estimated by Velocity (gm/cc) ( K0*10-4)
by noise empirical Vs30 (m/sec)
survey (Hz) relation (Hz)
Sl. Lat Long Predominant Predominant Effective Effective Q-factor Damping Site
No. (°N) (°E) frequency frequency Shear Wave Density Factor Response
estimated estimated by Velocity (gm/cc) ( K0*10-4)
by noise empirical Vs30 (m/sec)
survey (Hz) relation (Hz)
100 26.137 91.782 1.5 1.40 251.57 1.9769 202.89 5.88 4.65
101 26.148 91.811 2.2 1.19 271.74 1.9105 242.91 4.54 4.55
102 26.152 91.793 2.4 1.40 271.74 2.0005 253.05 4.36 4.45
103 26.159 91.788 1.6 2.63 304.88 2.066 209.14 4.70 4.14
104 26.133 91.737 3.6 ——— 294.12 2.0375 306.17 3.33 4.24
107 26.123 91.829 2.4 1.10 268.82 1.8582 253.05 4.41 4.64
108 26.163 91.707 3.6 2.63 304.88 2.0660 209.14 4.70 4.14
110 26.144 91.771 3.4 ——— 264.90 1.8907 298.06 ——— ———
111 26.135 91.758 7.2 2.63 284.09 2.1160 424.09 2.49 4.23
112 26.160 91.742 0.6 1.14 270.27 1.9162 131.90 8.42 4.56
113 26.121 91.771 0.9 0.91 222.22 1.9525 159.59 7.03 4.37
115 26.147 91.752 4.8 ——— 265.96 1.9195 350.50 2.57 4.59
116 26.102 91.751 2.2 0.91 234.74 1.8980 242.91 5.26 4.91
117 26.113 91.738 0.9 1.40 267.38 2.1150 159.59 7.03 4.37
118 26.178 91.830 4.9 2.63 243.90 1.9205 353.92 3.48 4.79
119 26.203 91.792 2.4 ——— 255.10 2.0740 253.05 4.65 4.51
120 26.182 91.814 4.7 1.88 255.10 1.8165 347.05 3.39 4.82
Chapter 3 53
Sl. Lat Long Predominant Predominant Effective Effective Q-factor Damping Site
No. (°N) (°E) frequency frequency Shear Wave Density Factor Response
estimated estimated by Velocity (gm/cc) ( K0*10-4)
by noise empirical Vs30 (m/sec)
survey (Hz) relation (Hz)
121 26.190 91.801 0.9 ——— 270.27 1.9150 159.59 5.56 4.56
122 26.130 91.770 1.2 1.40 246.31 1.9615 182.69 6.67 4.72
123 26.191 91.767 1.6 ——— 233.92 1.7560 209.14 4.91 5.12
124 26.181 91.813 2.9 1.4 259.07 1.8665 276.59 4.19 4.72
125 26.130 91.822 2.1 3.09 219.78 1.7147 237.66 4.59 5.35
126 26.171 91.830 3.5 ——— 265.49 1.6054 302.15 2.24 5.03
129 26.102 91.810 1.1 3.09 263.16 2.0140 175.37 6.50 4.51
130 26.173 91.748 1.5 0.91 260.42 2.1005 202.89 5.68 4.44
131 26.122 91.749 3.1 2.63 277.78 1.9620 285.40 3.78 4.44
133 26.111 91.769 1.1 1.14 250.00 1.7012 175.37 5.47 5.03
134 26.132 91.808 1.1 0.97 292.40 1.9035 175.37 5.85 4.40
135 26.154 91.837 1.1 0.91 257.73 1.8435 175.37 6.64 4.76
136 26.144 91.830 1.1 1.14 273.22 1.9915 175.37 6.26 4.45
137 26.135 91.808 4.2 ——— 304.88 2.0480 329.18 2.99 4.15
138 26.135 91.808 4.6 ——— 304.88 2.0480 329.18 2.99 4.15
139 26.114 91.817 1.6 1.32 282.49 1.8520 209.14 5.08 4.54
140 26.192 91.752 1.6 ——— 264.55 1.9700 209.14 5.42 4.55
141 26.178 91.758 2.2 ——— 284.09 1.9279 242.91 4.35 4.43
Table 3.2: Comparison of site response calculated from empirical relation and strong
ground motion data
Station (Latitude, Longitude) Predominant Empirically Site Response
Frequency Calculated from strong
Site ground motion
Response data
AEC (26.141910 N, 91.66107oE) 7.2 4.97 5.60
AMTRON (26.185730 N,
91.7861oE) 6.2 5.12 4.38
Cotton College Guwahati
(26.185860 N, 91.74524oE) 0.8 4.41 3.98
0
IIT Guwahati (26.18748 N,
91.69056oE) 0.6 6.63 6.73
0
Irrigation (26.184713 N,
91.77269oE), 5.4 4.88 5.45
54 Seismic Microzonation Atlas
Figure 3.8 Shear wave velocity (Vs30) contour map of Guwahati Region
Chapter 3 55
Site classification in Guwahati Region is done on the basis of shear wave velocity,
site response, predominant frequency and factor of safety. These themes are
integrated using multi criteria decision making method (Saaty, 1980). Figure 3.10
depicts the site classification of Guwahati Region. Range of effective shear wave
velocity and predominant frequency are given in Table 3.3.
Table 3.3: Shear wave velocity and Predominant frequency of each site class
4.1 INTRODUCTION
The introduction of performance based earthquake resistant design for buildings and
other civil engineering structure has increased the need for simulating realistic ground
motions. Combined with recent developments in software tools and structural modeling
techniques for time domain transient nonlinear dynamic analysis, the use of simulated
time histories of ground motion has gained major importance. Although the use of
recorded ground motion under condition similar to the design earthquake is appealing,
there may never be an adequate suite of such data in terms of tectonic structures,
earthquake size, local geology, and near field conditions. The variability of ground
acceleration traces depends on the source, propagation path, and site characteristics.
A number of methods have been developed and are in use for the adjustment of
selected recorded time histories to provide conformity to site condition (so-called
accelerogram scaling) and other spectrum technique that aim to generate time
histories of ground motion whose response spectra match the design response
spectrum. There are several approaches to the modeling and simulation of strong
ground motion taking into account the physics of source and propagation process.
Each approach accommodates, to varying degrees, the seismic wave radiation from
a fault rupture propagation through the crust, and modification by local site conditions.
we perform the simulation and modeling of ground motion by two approaches, first by
recorded strong motion data of two earthquakes considering omega (ω k = 2π f k )
squared circular crack source model (Brune,1970) and second by generating synthetic
seismogram using continuous wave number integration method (Wang and Herrmann,
1980).
58 Seismic Microzonation Atlas
In the previous chapter, it has been estimated that propagation path term P(rij,fk), and
site effect term SIj(fk) follow the equation (3.2). Brune (1970) gives the source term
as,
RθφF (2πfk )2
⋅ Moi
SO j ( fk ) = ................................................ (4.1)
(
2 4πρβ3 ) fk
1 + γi
fci
where Moi , fci and gi control the Brune’s source model (Brune, 1970). We initialize
Moi by computing the value of M0 in dyne-cm for a given Mw using Kanamori’s relation
(Hanks and Kanamori, 1979),
2
Mw = log (M0 ) − 10.73 ............................................................... (4.2)
3
The initial values of fci and gi are assumed to be 0.1Hz (Dutta et al., 2003) and 2
(Hwang and Huo, 1997) respectively. Estimation of these values is carried out in an
iterative fashion. Iteration continues until the difference between the observed and
the simulated source spectra tends to approach a minimum value. This exercise
was carried out for 2 events at five stations as well as for borehole locations those
coming in azimuthal range of these events. The location of all the five stations and
200 boreholes are shown in Figure 4.1.The simulated moment rate spectra and the
spectral acceleration have also been compared with those extracted from observed
data at the respective stations. The site response, source spectra, moment rate
spectra and the spectral acceleration at station IIT Guwahati (26.19° N, 91.69°E) are
furnished in Figure 4.2.
At IIT Guwahati, the site response as given in Figure 4.2(a) has been computed for
the event recorded in the azimuthal direction of 306.15°N. The observed and simulated
source amplitude spectra in Figure 4.2(b) find good match, the same has been
.
observed for moment rate spectra, Moi in Figure 4.2(c) and the spectral acceleration
in Figure 4.2(d). A corner frequency 5.4Hz for this event of Mw 3.1 with a scalar
moment of 2.5 × 1021dyne-cm seems to have replicated the observed spectra to a
great extent with the resulting stress drop of 1.8 bars.
Chapter 4 59
Figure 4.1 Borehole location, strong motion station and ambient noise survey sites
in Guwahati Region
Figure 4.2 (a) Site Response for the source azimuth 306.15°N, (b) Observed and simulated
source spectra, (c) Computed and best fit moment rate spectra, and (d) Observed and
simulated spectral acceleration with source parameters at IIT Guwahati (26.190 N, 91.69oE).
60 Seismic Microzonation Atlas
Table 4.1: List of source parameters: Moment Magnitude (Mw), Corner frequency
(fc), Scalar moment (M0) and Stress Drop (Ds) determined for representative events
in this study.
Sl. No. Mw fc M0 Ds
(Hz) (dyne-cm) (bars)
1 3.1 2.20 2.5E+21 1.8
2 5.1 5.04 4.8E+22 0.30
The value of corner frequency fci can also be computed using the relation
(Hwang and Huo, 1997),
1
∆σ 3
fci = 4.9 ∗10000 ∗β∗ ............................................................. (4.3)
M0
Equation (4.2) given by Kanamori can be used to compute M0 from ML. To estimate
the value of stress drop (Ds) we used the relation (Kanamori and Anderson, 1975),
3 16∆σ
log (M0 ) = log (S ) + log
2 3
.................................................. (4.4)
7π 2
where S represents the surface area, that can further be calculated using the relation
given by Bath and Duda (1964).
The corresponding spectral acceleration computed at AMTRON along with the PGA
values is presented in Figure 4.3. Figure 4.4 shows the site response and spectral
acceleration at borehole location 1 (26.141670N, 91.66125oE). Figure 4.5 represents
PGA distribution map of Guwahati Region considering omega squared circular crack
model. The highest PGA of 0.53g being observed at IIT Guwahati and 0.49g at AEC,
diminishing towards the east where the value reached 0.26g at Cotton College
Guwahati. This aspect is further analyzed through Green’s function simulation and
attenuation study from strong motion accelerometric data.
Figure 4.3 Simulated spectral acceleration for the Scenario Earthquake Magnitude
(SEM) Mw 8.7 nucleating at the focus of 1897 Shillong earthquake at AMTRON
Figure 4.4 (a) Site response, (b) Spectral acceleration at the bore hole with
the spatial coordinate 26.141670N, 91.66125oE
62 Seismic Microzonation Atlas
Figure 4.5 Spatial distribution of PGA in the Guwahati Region estimated from the
Simulation of scenario earthquake of Mw 8.7 nucleating from 1897 Shillong earthquake
Wave number integration method or Green’s function approach is useful if the range
of possible fault rupture history is narrow enough to functionally constrain the predicted
strong ground motion as is the case here. For the computation of synthetic
accelerogram impulsive source has been used as a first approximation for the near-
field effect. The wave number integration method of Herrmann and Mandal (1986) is
then followed.
The generation of synthetic seismograms for point sources in simply layered structures
has made rapid advances in the past decade. Two approaches, involving Laplace
transform and Fourier transform techniques, are actively being pursued. The Laplace
transform or Cagniard-de Hoop technique, usually referred to as the generalized ray
method (Helmberger, 1968), constructs the solution by tracking the individual seismic
arrivals ray by ray from the source to receiver. This method is valid at high frequencies
Chapter 4 63
and works well at predicting particular phases, but is poorly suited to models with
many layers and larger distances when a complete seismogram is desired. The
other approach involves expressing the solutions in terms of a double integral
transformation over wave number and frequency (Hudson, 1969). The complete solution
rather than individual rays, is considered in such a full wave theory approach. This
method can handle a larger number of plane layers, but requires considerable
computational effort, especially at high frequencies.
ZDD, ZDS, ZSS, RSS, RDS, RDD, TSS and TDS are referred as Green’s Functions.
It is apparent, that these three components are necessary to represent the P-SV
motion for any shear-dislocation source (Harkrider, 1976). RSS and RDS in equation
(4.6) also include the near-field terms. These terms decrease faster than the others
and therefore, are important only at short distances.
64 Seismic Microzonation Atlas
where S(ù) is the Fourier spectra of the impulse source function as described by
Herrmann (1979).
The S-wave part of the accelerogram is convolved with site response of each station
to obtain the response on engineering bed rock.
We have used the crustal model derived by inverting Rayleigh wave dispersion curve,
(Mitra et al., 2006). The model and dispersion curve, which has been taken for this
study are shown in Figure 4.6(a-b). This model is derived by inverting Rayleigh wave
dispersion curve. In each layer of the model the seismic energy attenuation due to
both absorption by intrinsic anelasticity and scattering by heterogeneities is
parameterized by the quality factor of S-wave Qβ and of P-wave Qα . They are assumed
to be frequency independent. Nath et al. (2005) has derived Qβ for Sikkim. The
relation between the P-wave quality factor Qα and the S-wave quality factor Qβ
proposed in the literature (Anderson et al., 1965; Kijko and Mitchell, 1983) implicates
Qα = k Qβ with 1 ≤ k ≤ 2.5 . We have chosen k =2.5, so that in our study the
effect of attenuation for a given Qβ can be considered to be minimal.
The SEM as considered in the foregoing analysis is assigned with the focal mechanism
of 1897 Great Shillong Earthquake for the synthesis of strong motion acceleration at
all the borehole locations in the study region. Fault plane parameters that have been
used for the Green’s Function synthesis (after Bilham et al., 2003) corresponds to a
Chapter 4 65
slip of 16m on a fault plane striking ESE for 110km and dipping SSW at 57°beneath
the northern edge of the plateau, slip on the plane extends from 9 to 45km beneath
the surface, with a rake of 76°. The dip of the fault and the rake of the slip are
constrained within 15°by the distribution of shear strains.
Our next step is to validate the PGA calculated by simulation using Brune’s source
with the PGA estimated by Synthetic Accelerogram Simulation. Hence the PGA is
calculated using synthetic seismogram simulation at 200 borehole locations as well
as at five strong motion instrument locations. Figure 4.7 shows the calculated Green’s
Function at borehole location 141 (26.11113oN, 91.79567oE). This computed synthetic
accelerogram was band pass filtered between the frequency range 0.5Hz to 25Hz.
Figure 4.8(a)-(e) show synthetic accelerogram convolved with site response and
simulated spectral acceleration at borehole location 141(26.11113oN, 91.79567oE ).
Figures 4.9 and 4.10 represent Green’s function and simulated spectral acceleration
at IIT Guwahati station. Figure 4.11 represents spatial distribution of simulated PGA
using Green’s function simulation in Guwahati Region. A comparative plot of the PGA
estimated by point source approximation w 2 - circular crack source model, and by
Wave Number Integration is presented in Figure 4.12 with a good 1:1 correspondence
observed in the cross plot with + 1 standard deviation bounds. However, in case of a
homogeneous, elastic 2D or 3D horizontally stratified half space the wave fields
radiated by an arbitrary internal point source may be described as a superposition of
cylindrical waves (Green’s function) characterized by horizontal wave numbers
pertaining to different directions.
Response at any seismometer is the convolution of source, site and path terms.
Path effect is the one which is dependent on the media in which the wave traverses.
In Green’s function simulation, superposition of impulsive layered earth responses
are considered with the fault attitude and Haskel source approximation, thereby
estimating realistic near - field effect having better accuracy at lower frequencies.
The radiation pattern in case of Brune source approximation of w 2 - circular crack
source model does not incorporate layered earth responses rather a simplified
reciprocal term for the path effect in an otherwise homogeneous half space for free-
field modeling especially at higher frequencies. Therefore, simulation by Brune’s
source approximation provides over-estimated ground motions at seismic stations.
66 Seismic Microzonation Atlas
Figure 4.6 (a) Dispersion curve for the velocity model (after Mitra et al., 2006),
(b) Inverted shear wave crustal model used for continuous wave number integration
Chapter 4 67
Figure 4.7 Green’s functions at the bore hole no. 141(26.11113oN, 91.76oE)
68 Seismic Microzonation Atlas
Figure 4.8 (a) Radial component of synthetic accelerogram, (b) Transverse com-
ponent of synthetic accelerogram, (c) Radial component of synthetic accelerogram
convolved with site response, (d) Transverse component of synthetic accelerogram
convolved with site response, (e) Simulated spectral acceleration for Mw 8.7 at
borehole no. 141(26.11113oN, 91.79567oE).
Chapter 4 69
Figure 4.11 Spatial distribution of PGA in the Guwahati Region estimated through
simulation of the scenario earthquake of Mw 8.7 by continuous Wave Number
(F-K) integration method (Green’s Function Approach)
Figure 4.12 Scatter plot between PGA (by Brune source approximation) and PGA
(by F-K Integration) depicting data clustering around the 45°, 1:1 correspondence line.
72 Seismic Microzonation Atlas
A ground response analysis consists of studying the behavior of a soil deposit subjected
to an acceleration time history applied to a layer of the profile. Ground response
analyses are used to predict ground surface motions for the evaluation of amplification
potential and for the development of design response spectrum. In the present study,
ground response spectra have been computed using WESHAKE91 software (from
the Waterways Experiment Station of US Navy) for the analysis of non-linear behavior
of soil very often approximated by a ‘linear-equivalent’ method that uses an iterative
procedure to adapt the soil parameters (i.e., rigidity and damping) to the actual strain
it undergoes. The behavior of soil under irregular cyclic loading is modeled by using
modulus reduction (G/Gmax) and damping (â) vs. strain curves, by calibrating with a
standard database input. The degradation curves for sand and rock used for the
present work are those proposed by Seed and Idriss (1970) and Schnabel (1973)
respectively.
Figure 4.13 (a) Absolute and (b) Normalized spectral acceleration at site class IIIA
Chapter 4 75
Figure 4.14 (a) Absolute and (b) Normalized spectral acceleration at site class IIIB
Chapter 4 77
Figure 4.15: (a) Absolute and (b) Normalized spectral acceleration at site
class IIIC
We have already observed in our previous analysis that the spectral acceleration
depends on site amplification, topography, source azimuth and the local site
conditions and hence it became necessary to work out a site specific attenuation
relation.
We started with the general form of equation given by Campbell’s attenuation law
(Campbell, 1997) for spectral acceleration as,
where SAĤ is the horizontal spectral acceleration, AĤ is the PGA, SSR and SHR are
variables representing local site conditions for soft rock and hard rock respectively, D
is the depth to the basement rock and fSA is a function of D. We introduced a term for
site amplification to take into account local site conditions. Our established second
order attenuation relation, therefore, takes the following shape,
30
ln (PGA) = ln(SA) – a1 – (a2 + a3 m)r –a4Vs -a5 ln(SR) .....................(4.10)
with SEM 8.7 .The median PGA is tabulated in Table 4.3 along with PGA values
computed by the simulation as well as those estimated using various global
relationships. The PGA estimated through simulation and spectral attenuation relations
closely follow each other, both being consistent with the occurrence of a scenario
earthquake of great magnitude in the Guwahati Region.
Figure 4.17 Comparative plots of simulated spectral acceleration for Mw 6.0 using
the Brune’s source model and by the spectral attenuation relation at AMTRON station
Table 4.3: Comparison of PGA values at all the stations computed by Campbell’s
Relation, Simulation of spectral acceleration, and the 2nd order local spectral
attenuation laws derived in this study
Cotton
College 0.21 0.26 0.23 0.21 0.26 0.26
Guwahati
5.1 INTRODUCTION
The GIS framework allowed us to account for added levels of details and complexity.
It is very important that the relevant data layers be consistent in their level of detail,
in order to successfully combine them and cross analyze them in the pair-wise
comparison process. A schematic process flow is shown in Figure 5.1. Prior to
presenting the methodologies used to produce the integrated seismic hazard maps,
a review of all the thematic coverage with their analytical detail is presented here.
84 Seismic Microzonation Atlas
5.2 LANDUSE
Guwahati, the capital of Assam and gateway to the east has continued to hold a
place of prominence throughout the history due to its very strategic location. Nestled
on the southern bank of the mighty Brahmaputra, where the river gets probably the
narrowest in its westward journey, and surrounded all around by hills, Guwahati has
been growing in prominence all along, and today is one of the fastest growing city of
India. It has established itself as the business and commercial capital of the northeast
India. Especially after shifting of capital of Assam from Shillong to Guwahati in 1972,
there has been spurt of growth in the city. Population of Guwahati has crossed 15
lakhs as per the 2001 census.
For preparation of Landuse map of the city, the 1990 landuse map prepared by the
Assam State Remote Sensing Application Center, Guwahati was used as primary
guide to carry out the baseline classification and ground checking to arrive at an
appropriate landuse. Keeping in mind the microzonation, old Survey of India topographic
sheets pertaining to 1927, and 1968 editions were used to obtain the possible land
fill areas. Satellite image, as described in the Para on Base Map was used for
detailed classification of different landuse classes. A detailed road to road GPS based
survey was conducted to arrive at cultural parameters such as residential, commercial
and institutional areas in order to further sub-classify the broad habitation areas
shown in greenish hue on the image. Ground checking was also conducted to collect
data on other landuse classes such as water bodies, swamps, and agricultural field.
A survey of the industrial areas was also conducted.
The landuse map was developed in-house at AMTRON and is shown in Figure 5.2.
86 Seismic Microzonation Atlas
In order to facilitate compilation and collation of data and maps from varied sources
on a single platform, the first task was to bring out a Base map of the area surrounding
Guwahati where microzonation was contemplated. For this purpose, DST took up
the matter with the Survey of India who provided the 1:25,000 scale topographic
sheets. The required information, coverage wise, such as roads, contour, railways,
streams and certain administrative boundaries were traced on a transparency and
digitized using onscreen digitization techniques. Further, IRS PAN and LISS III satellite
images were acquired in digital format for the Greater Guwahati region. The images
were corrected and georeferenced with respect to the topographic map points and
GPS readings taken of the ground control points. Thereafter, the LISS and the PAN
scenes were merged to produce a sharpened high resolution FCC image having a
ground resolution of 5.8m. The stream network, road network, river and water bodies
Chapter 5 87
were, thereafter updated using the satellite image. Further, an extensive GPS based
point survey was done in the greater Guwahati to capture all possible road networks,
important places, and other requisite features on the ground. Based on the above
exercise, a comprehensive map of Guwahati was prepared. This map has been used
all throughout the microzonation exercise. All other maps, data and GPS point surveys
were collated and corrected with respect to this base map, which made it possible to
weed out errors in various data sets as shown in Figure 5.3.
a) Precambrian granitic rocks forming the hill tracts and isolated hillocks and
b) Quaternary alluvium occupying the valleys, deposited over the uneven eroded
and faulted basement of granitic rocks.
The granite and granite gneissic rocks are well foliated and jointed, allowing deep
weathering along the joint and fault planes, and covered in most places by 1 to 3m
thick ferrugenous soil capping.
The Quaternary alluvium, perhaps, form the flood plain deposit of the Brahmaputra
River. This deposit could be classified into five aggradational units based on lithological
characters, state of weathering, order of superposition and unconformity between
them. They are, in order of increasing antiquity.
The Geotechnical properties of these units are of major concern in the work of seismic
microzonation.
The rocks of active flood plains consist of alternate layers of silt/fine sand and clay/
silty clay. The rocks of Digaru Surface are represented mostly by silt and fine sand
88 Seismic Microzonation Atlas
and are underlain by the Bordang Surface below an unconformity. The Bordang Surface
consists of white silty clay at the top and medium to coarse grained sand at the
bottom. The Sonapur Surface constitute the oldest unit of all the fluvial deposits of
Brahmaputra valley, mostly exposed near the foot hills showing contact with the
pediment forming colluviums or with the granite rocks. It comprises of bedded sand,
silt, and clay in varying proportion with maximum amount of clay. The Pediment
Surface is formed of weathered and eroded alluvial materials deposited mostly along
the foothills.
It is difficult to distinguish these litho-units at depth from the study of sludge collected
during digging of boreholes up to 120m depths for ground water exploration. However,
it has been seen that sand, silt, clay, and gravel alternate in irregular proportion with
extensive lateral variations. It may be mention here that seismic resistivity sounding
surveys by GSI could identify three layers of rocks at depth with resistivity
characteristics of 200, 100 and 25 Ohm-meter in the central and western part of the
valley.
Figure 5.3 Geology and Geomorphology map of Guwahati Region (Base Map)
The overall topography of the area is rugged with high relief due to presence of steep
sided hillocks carved out of the Meghalaya Plateau that occupies the southern and
the eastern fringe of the area under consideration. The denuded and continuous hill
tracts from Rani to Khanapara RF in the south and the Amchang Hills RF in the east
rise in altitude from 200 to 400m above MSL. Isolated hillocks within the valley
occasionally rise up to 300m above MSL. The highest hillock is located in the north-
western part of the area (Silapahar-381m) and the lowest elevated hillock, Odalbakra
– 145m, lies in the center of the area. The relative relief is high varying between 80
and 300m. The general elevation of the valley area varies from 25 to 50m above MSL.
a) Denuded hills,
It is important in the present context that many of the hillocks have been cut into
small terraces and are occupied for habitation, and the swampy lands/water bodies
have been filled up in many places for construction of houses due to tremendous
pressure of population. A few of the hillocks have scarp or steeply inclined face.
They are Engineering College hill, Gitanagar hill and Japorigong hill. The isolated
hillocks have NE-SW trending dendritic pattern of drainage of moderate density. Due
to deep weathering along joints and faults the hill slopes have become unstable
making them vulnerable to rock and debris slides, especially during heavy rains
accompanied by earthquake shaking as shown in Figure 5.3.
Basement contour map was prepared from the data provided by GSI obtained through
resistivity survey and that from the DGM obtained through drilling for ground water
Figures 5.4, 5.5 and 5.6. In general, depth of basement increases from east to west
in south bank and reaches up to –300m. In north bank basement depth reaches up
to –600m in Singimari – Dadara area and attains the shape of a circular basin in the
northwest corner of the area.
The basement forms N100W-S100E trending steep ‘V’ shaped valley with maximum
depth of –150m in between the Fatasil and Kalapahar hillocks that is traversed by a
fault. Another NE-SW trending valley occurs in between the engineering college hill
and the Fatasil hill that extends up to the junction of Neelachal hill and the Fatasil hill
with maximum depth of –250m. This valley is also underlain by a fault. A NE-SW
trending valley passing along the Tepar beel by the western side of the Japorigong
hills reaches the depth of –100m and it is relatively wide. An E-W deep valley passes
from Panjabri to Kahabari across the Deepor beel. Towards east the basins has
highs and lows but to the west the basin gradually get deeper and attains the depth
of –300m. A fault also runs along this valley. Thus it is seen that the valley area has
variable depth of basement having steep gradient along some zones like the western,
eastern and the southern margins of the Fatasil hill; western margin of the Kalapahar
hill and the northern periphery of the Rani-Khanapara hill tract. These zones will have
pronounced basin edge effect during earthquake shaking Figures 5.7 and 5.8 present
Basement Contour and zonation map of Guwahati Region.
Chapter 5 91
5.5 SEISMOTECTONICS
The city of Guwahati is located in an area surrounded on all sides by highly active
tectonic blocks as shown in Figure 5.9. Generally speaking the area is buttressed in
between the Himalayan collision zone to the north and the northeast, Indo-Myanmar
subduction interface of Indian plate to the east and the Meghalaya Plateau – Mikir
hills tectonic block to its south. Strictly speaking the Guwahati area falls in the
domain of Meghalaya Plateau and Mikir hill block. Juxtaposition of ongoing collision-
subduction tectonic processes has made the area one of the most intense seismic
zones of the world.
Analysis of contemporary tectonics in the region reveals that active Himalayan frontal
thrusts and cross faults cutting across these thrusts have been generating many
relatively shallow, small and moderate earthquakes. Strike-slip movements along a
NW-SE trending fault (Po Chu) in the Mishmi block had produced the 1950 Great
94 Seismic Microzonation Atlas
Assam earthquake (Mw 8.7) inflicting catastrophic damage in the Upper Assam
area. All the other thrusts/faults in this block viz. Mishmi thrust, Lohit thrust and
Tidding suture can be classified as capable faults.
The active subduction process along the Indo-Myanmar mobile belt and the conjugate
faults lying across this belt has been producing many large and major earthquakes
that shook the Guwahati region.
The tectonic block of Maghalaya Plateau-Mikir hill represents the northeastern most
exposed element of the Indian shield, occupying a crucial position in between the
northern collision and the eastern subduction zones of the Indian plate. This block is
under tremendous stress and is seismically active. The Guwahati area being located
in the northern margin of this block is vulnerable to severe earthquake damage. The
N-S faults cutting across the Plateau such as the Jamuna or Dhubri fault, Dhudnoi/
Chedrang fault, and Kulsi fault in the Meghalaya Plateau are very active; NW-SE
Kopili fault passing in between the Plateau and the Mikir hill and the NW-SE Bomdila
fault passing along the northern margin of Mikir hill, both traversing across the
Himalayan thrust and fold belt as well as the Indo-Myanmar mobile belt are very
important.
The Jamuna or Dhubri fault has been the source for 1931 (Mw 7.1) Dhubri earthquake.
Movement along the Dhudnoi or Chedrang fault generated the 1897 Great Assam
earthquake. The Kopili fault has the record of producing the 1869 Cachar earthquake
(Mw>7) and 1943 earthquake (Mw>7). Moreover, recent recording of earthquake events
clearly demonstrate that Kopili fault is highly active at present. Another important
active fault is the NE-Sylhet and its associated faults falling in the tectonic domain of
Bengal basin. This fault had generated the 1918 Srimangal earthquake (Mw 7.6).
Similar earthquake may affect the Guwahati city. The E-W Dauki fault system, a
regional structure of great importance, though seems to be dormant at present may
produce large earthquakes that may affect the Guwahati city. Focal mechanism
solutions of the past earthquakes reveal that most of the events were due to strike
slip motion in these terrains (Nandy and Dasgupta, 1991; Nandy, 2001).
Thus, it is seen that the Guwahati city area is vulnerable to catastrophic near source
great and large earthquakes. It may be mentioned here that recent release of stress
Chapter 5 95
along a 1200km long subduction interface in the southern part of the Indo-Maynmar-
Andaman-Sunda subduction zone by the 26th December, 2004, Mw 9.0 earthquake
has made it highly probable that next rupture may take place along the northern
sector of the subduction zone in Indo-Myanmar region.
The thickness of weathered zone/soil is fairly high in the granite hillocks in and
around the Guwahati city. Many of the landslide incidences (baring a few in the
Kalapahar area) are due to anthropogenic activities. A total of six landslides of
different categories have been recorded. Most of them are slump type followed by
debris slide and rock fall. Except one located SW of Deepor Beel, all are fresh
slides. Slides located at Kalapahar, Kahangkar basti, Rupnagar and 10th mile on G-
S road are of planer failure type whereas others are circular slip type.
Twenty two landslides are concentrated in and around the Guwahati City, especially
in its central part. Kalapahar area experienced 10 slides followed by Dhirenpara
having 4 slides. Area along the G-S road has the record of 9 landslides.
A component subjected to a solitary load will be considered in the first instance. This
load is interpreted in the context of the component’s nature and duty. Therefore, load
usually implies a transverse force in the case of a beam component, or a longitudinal
compressive force in a column, or a torque in the case of a shaft, or a pressure in a
fluid containment vessel, and so on.
There are two completely different manifestations of the load which have important
consequences for the component:
· The extrinsic actual load is the load exerted on the component by its surrounds.
· The intrinsic maximum load is the largest load that the component can withstand
without failure; the maximum load is a property of the component, a function of
its dimensions and material.
Chapter 5 101
Clearly, a component is safe only if the actual load applied to the component does
not exceed the component’s inherent maximum sustainable load. The degree of
safety is usually expressed by the safety factor ‘n’ given as,
Accordingly the study region is classified in two zones namely safe and unsafe. The
Factor of Safety contour and distribution map is shown in Figures 5.17 and 5.18
respectively.
6.1 INTRODUCTION
5. Geomorphological characterization.
The matrix developed by pair-wise comparisons between the factors can be used to
derive the individual normalized weights of each factor. It is performed by calculating
the principal eigenvector of the matrix. This results in a matrix of values that are in
the range of 0 to 1 and sum to ‘1’ in each column. The weights for each attribute can
be calculated by averaging the values in each row of the matrix. These weights will
also sum to ‘1’ and can be used in deriving the weighted sums of rating or scores for
each region of cells or polygon of the mapped layers.
Since the values within each thematic map/layer vary significantly, they are classified
into various ranges or types, which are known as the features of a layer. These
features are then assigned ratings or scores within each layer, normalized to ensure
that no layer exerts an influence beyond its determined weight. Therefore, a raw
rating for each feature of every layer is allocated initially on a standard scale such as
1 to 10 and then normalized using the relation,
R j - Rmin
xj = ................................................................. (6.1)
Rmax - Rmin
where, Rj is the raw score, Rmin and Rmax are the minimum and maximum scores of a
particular layer.
Chapter 6 109
For Seismic Microzonation and Hazard delineation of the above themes, both
Geomophological and Seismological themes are reclassified into a 1st phase
geohazard map and 2nd phase Seismic Microzonation map with PGA distribution for
a SEM of Mw 8.7. A typical two phase Microzonation procedure from hazard zonation
to regional hazard zonation mapping on GIS platform and finally to Seismic
Microzonation is shown in Figure 6.2.
The Geological vector layers that have been used for microzonation include Geology
and Geomorphology (GG), Basement (BS), Landslide hazard (LS) and Land Use
(LU), while the seismological themes are Shear wave velocity (Vs30), Site Response
(SR), Peak Ground Acceleration(PGA), Predominant Frequency (PF), and Factor of
Safety (FS).
The Geological and Seismological themes are weighted in scale of 9:1 depending on
their contribution to the seismic hazard, the highest being attached to Geological &
Geomorphological layer with a normalized weight of 0.2000. Basement has got the
next weightage with a normalized value of 0.1778; Landslide hazard and Landuse
have got the next weightage with normalized value of 0.1556 and 0.1333 whereas the
Vs30, PGA, SR, PF and FS are assigned the values 0.1111, 0.0889, 0.0667, 0.0444
and 0.0222 respectively.
A matrix of pair-wise comparisons (ratio) between the factors is built, which is used
to derive the individual normalized weights of each factor. The pair-wise comparison
is performed by calculating the principal eigen vector of the matrix and the elements
of the matrix are in the range of 0 to 1 summing to ‘1’ in each column. The weights for
each theme can be calculated by averaging the values in each row of the matrix.
These weights will also sum to ‘1’ and can be used in deriving the weighted sum of
rating or scores of each region of cells or polygons of the mapped layers. Since the
values within each thematic map/layer vary significantly, those are classified into
various ranges or types known as the features of a layer. These features are then
assigned ratings (r) or scores within each layer, normalized to 0-1 as shown in Table
6.1 below.
Chapter 6 111
Factor of 0.0222 £1 1 1
Safety (FS) >1 0 0
The seismic hazard zonation map is obtained through the integration of all the above
themes using the following relation
Five zones are mapped as shown in Figure 6.3 where average Hazard index is 0.55,
0.44, 0.35, 0.26 and 0.16. We termed these zones as very high, high, moderate, low
and very low hazard regions. Figure 6.4 shows microzonation map of Guwahati region
overlaid by PGA computed by Green’s Function approximation.
114 Seismic Microzonation Atlas
Figure 6.3 Seismic Microzonation Map of Guwahati Region (using PGA computed by
F-K integration)
Figure 6.4 Seismic Microzonation Map of Guwahati Region with PGA (Green’s Function
Approach) contours overlaid
CHAPTER 7
7.1 INTRODUCTION
Buildings and other structures in Guwahati city are facing highest seismic risk.
Therefore, collection of strong motion data and analysis of site response of Guwahati
city is essential for better estimation of hazard level in the Greater Guwahati region
and this will provide valuable input for improvement of earthquake resistant design of
structures. The geology and seismotectonics of the region have been well documented
by Nandy (2001). Figure 7.1 shows the seismicity map of the NE India and adjoining
region.
In this work, twelve numbers of Strong Motion Accelerograph stations are to be installed
at twelve sites in different part of the Greater Guwahati City area for recording free-
field ground acceleration during earthquakes. These twelve sites have been selected
on the basis of sub-soil and bedrock characteristics. RCC sheds as shown in Figure
7.2 have been constructed in nine sites and in all such sites accelerographs (Figure
7.3) have already been installed. At one of the sites, IIT Guwahati, first event has
been recorded on 18th September 2005. The source of the event has been located at
Indo-Myanamar border (24.64°N, 94.81°E). The magnitude and focal depth of the
event were Mw = 5.7 and 82km respectively. The peak ground acceleration of the
recorded motion is 0.013g.
116 Seismic Microzonation Atlas
Figure 7.1 Seismicity map of the NE India and adjoining region (after Nandy, 2001)
Chapter 7 117
Figure 7.2 SMA station at IIT Guwahati Figure 7.3 SM Recorder (Etna)
Subsequently, six events have been recorded in different strong motion stations of
the array established under the ongoing project. Details of the sources of these
events along with the corresponding magnitude have been furnished in the Table 7.1.
The sample ground motions recorded during these events are shown in the Figures
7.4 - 7.11.
Date Magnitude
Epicenter Focal Depth (km)
(Mw)
Figure 7.6 Strong motion accelerograms recorded at station – Cotton College, Guwahati
(26.187580°N, 91.743896°E) on 16-11-2006
Chapter 7 121
Figure 7.7 Strong motion accelerograms recorded at station – IIT Guwahati Hill Top
(26.187580o N, 91.743896 o E) on 06-11-2006
122 Seismic Microzonation Atlas
Figure 7.8 Strong motion accelerograms recorded at station – IIT Guwahati (MED)
(26.187481°N, 91.690567°E) on 06-11-2006
Chapter 7 123
Figure 7.9 Strong motion accelerograms recorded at station – CEO: Irrigation Guwahati
(26.184713°N, 91.772697°E) on 06-11-2006
124 Seismic Microzonation Atlas
In engineering practice, acceleration time history is the most sought after information
during an earthquake. It is well known that high frequency accelerations resulting
from complex faulting and medium irregularities are impossible to
Chapter 7 127
Figure 7.12 Average normalized pseudo-acceleration response spectrum for AEC Guwahati
site (26.141911°N, 91.661072°E)
Chapter 7 129
Figure 7.14 Average normalized pseudo-acceleration response spectrum for Cotton College
Guwahati site (26.187580°N, 91.743896°E)
Chapter 7 131
Figure 7.15 Average normalized pseudo-acceleration response spectrum for IIT Guwahati
hill top site (26.193830°N, 91.692047°E)
132 Seismic Microzonation Atlas
Figure 7.17 Average normalized pseudo-acceleration response spectrum for CEO- Irrigation
Guwahati site (26.184713°N, 91.772697°E)
134 Seismic Microzonation Atlas
Figure 7.18 Average normalized pseudo-acceleration response spectrum for RRL Guwahati
site (26.158409°N, 91.735542°E)
Chapter 7 135
Figure 7.19 Average normalized pseudo-acceleration response spectrum for S.D. Kalakhetra
Guwahati site (26.132675°N, 91.821747°E)
136 Seismic Microzonation Atlas
8.1 INTRODUCTION
The primary emphasis in disaster management has been placed on the performance
of structures, the failure of which would pose direct risk to life and property. Survey
and assessment of existing building stocks for earthquake vulnerability risk is
necessary to formulate the seismic hazard map of a city.
The seismic vulnerability study comprises mainly review of the existing buildings of
Guwahati in the light of guidelines for earthquake resistant construction in India,
construction practices being adopted in Guwahati urban area, building typologies,
designing of questionnaire for detail survey of buildings of the Guwahati municipal
area, selection of representative building samples for detailed analysis and NDT,
and creation of database. Subsequently, seismic vulnerability of existing building
stock estimated quantitatively and qualitatively. The quantitative approach covers
demand-capacity computation, while qualitative procedure estimates structural
scores using Rapid Screening Procedures (RSP). The results are mapped using
ArcInfo and GIS, which are later synergized with seismic hazard microzonation to
deliver seismic risk.
The specifications, which are predominantly seen in and around Guwahati are
underlined and put in bold letters.
Type-B: Ordinary brick building – Brick Masonry Wall – 6”X6” corner columns
with lintel bend and tie – timber trussed CI sheet roof, buildings of the large
block and prefabricated type, half-timbered structures, building in natural hewn stone
The majority of houses in the villages around Guwahati urban area is of ‘Assam
Type’ categorized in TYPE-B - (a) brick masonry wall, 6”X6” R.C. Column, lintel
band, tie and with Timber truss. (b) bamboo reinforced Biomass Wall cladding with
roofs made of thatch of bamboo supported on bamboo purloins. This building has
very low vulnerability.
been used for carrying out risk analysis given in the district-wise tables. The
Vulnerability Atlas, thus, provides ready macro-level information for use by the
authorities for natural disaster mitigation and preventive actions. The types of housing
as existing in each district has been taken from the 1991 Census of India and
categorized from vulnerability consideration. The vulnerability of these types to
various intensities of the hazards was estimated by and the damage risk in each
district has been presented in a separate table. The area of the district prone to
various hazard intensities has also been shown. In the present study we achieved
a preliminary seismic risk map using the PGA computed by F-K integration method
as shown in Figure 8.5. The Microzonation map was integrated with Demographic
distribution of the region. Table 8.1 is showing the normalized feature rating for risk
estimation of the region. The preliminary seismic risk map presented here is
recommended for decision making purpose.
Indian buildings built over past two decades are seismically deficient because of
lack of awareness regarding seismic behavior of structure, constant upgradation of
knowledge as regards earthquake resistant design and construction. Also seismic
design is not practiced in most of the buildings being built. It calls for seismic
evaluation of existing building stocks in an area.
Evaluation is a complex process, which has to consider not only the design of
building but also the deterioration of the material and damage caused to the building,
if any. The difficulties faced in the seismic evaluation of a building are manifold.
There is no reliable information/database available for existing building stock,
construction practices, in-situ strength of material and components of the building.
The seismic evaluation mainly relies on set of general evaluation statements. The
unavailability of a reliable estimate of earthquake parameters, to which the building
is expected to be subjected during its residual life poses another challenge.
Probabilistic approach to evolve needful parameters, would call for elaborate studies.
Hence, for preliminary appraisal, the ground motion parameters available in the
present code (BIS, 2002) have been estimated at the macro level. As regards the
effect of local soil conditions, which are known to greatly modify the earthquake
ground motion, experiences of ground accentuation and data generated through
collateral studies on site response have been considered. Also, in view of above
constraints, the present study is limited to seismic evaluation of representative
buildings of different typology viz. Type-A (Mud/RR Masonry, Adobe), Type-B (Brick
Chapter 8 143
Masonry Buildings), and Type-C (RCC Buildings), and projects a generalized pattern
of building response to future seismic ground motion in different wards/zones of
Guwahati urban area.
The methodology adopted in carrying out the survey and indexing damage potential
to various buildings is in line with the proposed methodology by Ministry of Home
Affairs (MHA), India and is based on the classification of buildings as per MSK
Intensity scale. The Rapid Visual Screening Form as shown in Table 8.1 was used
for collecting building data. A total of 4500 nos. of sample buildings spread over 60
ward were surveyed and the building data has been electronically warehoused for
future analytical work.
The study found that most of the residential buildings in Guwahati are constructed
based on socio-economic consideration rather than engineering approach. The
majority of building stock (about 80%) composed of Type-C buildings. The TYPE-A
and TYPE-B buildings in residential area (about 20 %) are though non-engineered,
they have earthquake resistant construction features inbuilt in them.
144 Seismic Microzonation Atlas
TABLE 8.1
The region delineated under Guwahati Metropolitan Area measures 262 sq. km.
includes Guwahati Municipal Corporation (GMC) area, North Guwahati Town
Committee, Amingaon Census Town and 21 revenue villages which are under :
Abhaypur, Rudreswar, Namati Jalah, Gouripur, Silamohekaiti, Tilingaon, Shila,
Ghorajan, Mikirpara, Kahikuchi, Kahikuchi, Mirjapur, Jugipara, Borjhar, Garal Gaon,
Dharapur, Janisimalu & Jansimalu (NC), Kalitakuchi(NC), Kharghuli, Bonda,
Bondagaon and Bonda Grant, (I&II) and Birkuchi. In addition to the above area,
there is Narangi Cantonment area, Guwahati Refinery (IOCL) area, and NF railway
Colony.
Table 8.2: The weights and normalized ratings for integrations of risk assessment
Figure 8.5 Preliminary Seismic Population Risk map of Guwahati urban center
(Using the PGA computed by wave number integration method)
Apart from filling-up of questionnaire for the selected buildings, surveyor has to
inspect the health of structure critically to assess its seismic resistance. In the
process, surveyor has to face several difficulties. The foremost problem is of
uncovering the structure. In many buildings the structure is concealed by architectural
finishes, and the surveyor will have to get into attics, and crawl spaces.
Non-availability of plans, and design calculations are yet another problem, and is
particularly frustrating with respect to reinforced concrete work. Assessing material
quality and associated allowable stresses are also difficult preposition, and one has
to rely on local available reports/information or otherwise one has to go for destructive
testing, which is seldom possible. Destructive and non-destructive testing of
reinforced concrete and masonry elements are necessary to determine strength
and quality of construction. The rebound hammer was used to assess the
150 Seismic Microzonation Atlas
The sample survey was carried out for about 4500 buildings spread over 60 wards
of Guwahati municipal area, out of which about 5% are of Type A, 15 % are of Type
B and 80% are of Type-C.
For quantitative approach, DCR computation has been used for mainly Type C
buildings and later it is related with the possible failure modes. Few Type B buildings
were taken up for DCR analysis to understand the vulnerability of typical masonry
structure existing in the area.
Chapter 8 151
Since earthquake is a random process, all the load bearing walls in a structure are
to be evaluated for their shear resistance. The demand placed by an earthquake
i.e. lateral forces at various levels, as per IS:1893-2002, along with gravity load
calculations were carried out for sample buildings, and later check in terms of
Demand Capacity ratio (DCR) for shear resistance, combined stress, overturning,
and stability of non-structural failures for long and short walls. The capacity of wall
is defined as its allowable stress depending upon mortar type in accordance with
the relevant codal provisions. The DCR greater than unity, indicates that the building
is seismically vulnerable in respective criterion, whereas DCR less than one implies
the building to be safe under earthquake loads. As indicated earlier, earthquake
demands for better shear resistance and hence the DCR in shear should be less
than one, otherwise the building will have diagonal (X) cracking. The DCR greater
than one for combined stresses means that the building is not even designed for
gravity loads and would lead to collapse on seismic shaking. The failures in
overturning corroborate falling of walls. The check for non-structural element implies
the falling hazard of parapet wall. The above analogy has been used to estimate
seismic vulnerability in terms of various failure modes i.e. collapse, excessive
cracking, falling of walls including parapet walls. The city of Guwahati fortunately
does not have many load bearing structures, and the only area where such structures
are located are the NF Railway colony, where all the official residences are on load
bearing wall structure ranging from one to even three stories.
been considered for analysis. The analysis directly computes member end forces
and then each member is designed for worst load combination. The design module
of analysis engine gives the longitudinal and transverse reinforcement for each
member. This reinforcement corresponds to the demand of a member due to
earthquake forces, whereas the actual reinforcement provided in a particular member
would correspond to capacity. In order to calculate the DCRs, the calculated
reinforcement of structural members has been compared with provided
reinforcement. The DCRs for longitudinal and transverse reinforcement reflect
DCRs for flexure and shear of member. The DCRs calculated for flexure and
shear give the idea about inherent ductility and strength of member to ensure safety
& serviceability during severe shocks.
The DCR greater than one for flexure indicates that the longitudinal reinforcement
in columns & beams are inadequate leading to failure. The possibility of failure of
such buildings is excessive cracking leading to collapse. Whereas DCR greater
than one in shear indicates that the lateral ties provided are not sufficient leading to
brittle failures i.e. catastrophic failure. In this case, there is possibility of diagonal
cracking in structural elements. The check for non-structural element implies the
falling hazard of parapet wall. Based upon above analogy, DCRs for flexure, shear
and non-structural members leading to estimate seismic vulnerability in terms of
failure modes i.e. excessive cracking, diagonal cracking and falling hazard
respectively for few representative RC buildings have been computed.
Chapter 8 153
Figure 8.6 DCR ANALYSIS - Few TYPE-C Buildings (Carried out for 25 buildings)
154 Seismic Microzonation Atlas
The prognostic damage scenario of a ward reflects the structural and non-structural
damages induced in the existing building stocks. The damage scenario of a ward
given here is based on representative building surveyed for different building
typologies and not final and conclusive. It would require the inputs from all the
working subgroups to finally project prognostic damage scenario.
Based on the survey and analysis of data, the seismic vulnerability of Guwahati
urban area obtained through qualitative approach. It is found that Type-A, Type-B
& Type-C buildings are 10%, 40%, & 30% vulnerable respectively from engineering
(design and construction) perpectives.
In order to present the prognostic damage scenario for Guwahati urban area using
quantitative approach, the failure modes of different type of buildings are being
collated. All the Type-A houses are 10% vulnerable since they are built from materials,
which has very low vulnerability to earthquake damage. Unlike in other parts of
India – rural housing in general in northeastern are very light and devoid of use of
load bearing walls and stones. The roofing is very light rendering them very safe
against earthquake. For Type-B buildings, the postulated failure modes have been
categorized as Excessive Cracking (EC); Falling of Walls (FW); Falling Hazard of
non-structural members (FH); and combination thereof - Excessive Cracking +
Falling of Wall (EC+FW); Excessive Cracking + Falling Hazard (EC+FH); Falling of
Wall + Falling Hazard (FW + FH); Excessive Cracking + Falling of Wall + Falling
Hazard (EC + FW + FH); and safe buildings (which do not have any failure). At the
first instance, ward wise seismic vulnerability is being derived, and later the ensemble
is projected to present prognostic damage scenario for Guwahati urban area.
Similarly, the various failure modes for assessing seismic vulnerability of Type-C
buildings are identified as excessive cracking (EC), diagonal cracking (DC); falling
hazard of non-structural members (FH); and combination thereof and safe buildings.
The prognostic damage scenario for Type-C buildings in Guwahati urban area also
being worked out and mapped to present Prognostic Seismic Vulnerability Map of
Buildings in Guwahati. The vulnerability map can be effectively used to project the
risk associated with existing building stock in Guwahati Urban area. Further, these
Chapter 8 155
maps may act as guidance for future planning, risk reduction and disaster mitigation
and management. The Prognostic Seismic Vulnerability Map needs to be
concurrently and continuously updated with the availability of more geo-scientific
and engineering data.
CHAPTER 9
Concluding Remarks
The seismic ground motion hazard mapped in the Guwahati Region is microzoned
with local and regional site conditions incorporated through GIS. The process of
overlaying, union and finally integration of various geomorphological and
seismological thematic maps are complicated spatial operations that are optimally
performed on GIS environment. The union of geological and geomorphological units,
landuse, landslide hazard and basement zonation provide the background site
condition of the Guwahati Region on which the seismological (site response,
predominant frequency, peak ground acceleration) and geotechnical
(Vs30 and Factor of Safety) attributes are overlaid.
The seismicity analysis of northeast India has been performed and the entire northeast
is classified into four source zones namely zone 1 (Himalayan source), Zone 2 (Mishmi
block), Zone 3 (Eastern Boundary Thrust Zone) and Zone 4 (Shillong plateau). Each
has a potential of generating a great earthquake of magnitude greater than eight and
that can devastate the entire region. An examination of the spatial distributions reveals
relatively higher b values on the upper Tibet-Himalayan area. Low values are seen in
the thrust zones especially in the Central Eastern Boundary Zone of Arakan Yoma
Range indicating high stress buildup in the zone. The historical Shillong earthquake
of Mw 8.7 is employed as the scenario earthquake. The presented Seismic
Microzonation of Guwahati region is, therefore, a seismic scenario for an earthquake
magnitude Mw 8.7 that represents a near source which may nucleate from the focus
of the 1897 Great Shillong Earthquake historically known for its catastrophic effect.
The Microzonation map given in Figure 6.3 obtained by integrating geological and
seismological themes presents five hazard zones namely very low with the average
hazard level 0.16, low with average hazard level 0.26, moderate with average hazard
158 Seismic Microzonation Atlas
level of 0.35, high with hazard level 0.44 and very high having hazard level greater
than 0.55.
Very high hazard zone comes under the site class IIIA where average shear wave
velocity (Vs30) is 228m/s with site amplification is larger than 5.5 and average
predominant frequency is 1.15Hz. Surface geology of this zone is predominantly
active flood plain with pediment, Sonapur surface and some patches of Bordang
surface. In this zone PGA level is very high and greater than 0.81g. Western part of
Guwahati region comes under this zone that covers the area namely Bhattapara,
Garal, Bulla, Singibari, Azpa, Baharbari, Khanapara etc.
High hazard zone comes under site class IIIB where average shear wave velocity
(Vs30) is 260m/s and average predominant frequency is 2.93Hz with average site
amplification of 4.81. Here the surface geology of this zone is predominantly Sonapur
surface, water bodies and pediments. Average PGA level in the zone is 0.52g. This
zone resemble the very high hazard zone and it covers area like Lachitpur, Teteliya,
Mirzapur, Dekapara etc. Guwahati airport incidentally falls within the high hazard
zone.
Moderate hazard zone comes under site class IIIC with average shear wave velocity
(Vs30) is 296m/s. Surface geological feature of the zone is predominantly Bordang
surface with average site amplification of 3.1. Average predominant frequency and
PGA are 5.41Hz and 0.38g respectively. Moderate hazard zone covers the entire
Bordang surface and its habitants like Dispur, Gopinath Nagar, Rehabari, Shantipur,
Chatribari with important installations namely AMTRON, BSNL, Guwahati University
and Irrigation department.
Low and very low hazard zone are under site class IIID which is almost covered by
Hillocks. Average shear wave velocity in the zone is 340m/s while average predominant
frequency is 7.1 Hz. This zone has very low PGA level of 0.11g with low site
amplification of the value less than 1.5. This hazard zone encompasses the region
near eastern flank of Brahmaputra River, in the periphery of hills like Fatasil hill,
Kalapahar hill and some area in the eastern part of the region.
Chapter 9 159
The hazard and risk map presented here may be useful for land use planning and
making hazard mitigation decisions. This map is generally better spatial representation
of seismic hazard including site specific analysis, and may be used for recognizing
hazardous areas at regional scale. The geologic site condition map is an initial model
to describe areas that may exhibit different seismic shaking characteristics during
future earthquakes.
Bibliography
Abrahamson, N.A., and W.J. Silva (1997), Empirical response spectral attenuation
relations for shallow crustal earthquakes, Seism. Res. Lett., 69, 94-127.
Acharyya, S.K., and K.K. Ray (1977), Geology of Darjeeling Sikkim Himalaya: Guide
to Excursion No. 4, 4th International Gondwana Symposium, New Delhi, 1-25.
Aki, K. (1965), Maximum likelihood estimate of b in the formula log N = a - bM and its
confidence limits, Bull. Earthq. Res. Inst. Tokyo Univ., 43, 237-239.
Aki, K. (1969), Analysis of the seismic coda of local earthquakes as scattered waves,
J. Geophys. Res., 74, 615-631.
Aki, K. (1988), Local site effects on strong ground motion, Earthquake Engineering
and Soil Dynamics II – Recent Advances in Ground Motion Evaluation, June 27-30,
Park City, Utah.
Aki, K. (1993), Local site effects on weak and strong ground motion, Tectonophysics,
218, 93-111.
Aki, K., and B. Chouet (1975), Origin of coda waves: Source, attenuation and
scattering effects, J. Geophys. Res., 80, 3322-3342.
Aki, K., and K. Irikura (1991), Characterization and mapping of earthquake shaking
for seismic zonation, Proc. 4th Intern. Conf. Seismic Zonation, EERI, Oakland, 1,
61-110.
Alvarez, L., F. Vaccari, and G.F. Panza (1999), Deterministic seismic zoning of eastern
Cuba, Pure Appl. Geophys., 156, 469–486.
Ambraseys, N.N. (2001), The earthquake of 10 July 1894 in the Gulf of Izmit (Turkey)
and its relation to the earthquake of 17 August 1999, J. Seismol., 5, 117-128.
Ambraseys, N.N., and C.F. Finkel (1995), The Seismicity of Turkey and Adjacent
Areas, EREN, Istanbul, p. 240.
Ambraseys, N., and R. Bilham (2000), A note on the Kangra Ms=7.8 earthquake of 4
April 1905., Curr. Sci., 79, 101-106.
Ambraseys, N.N., K.A. Simpson, and J.J. Bommer (1996), The Prediction of Horizontal
Response Spectra in Europe, Earthq. Eng. Struct. Dyn., 25, 371-400.
Anderson, J.G., and J. E. Luco (1983), Parametric study of near-field ground motions
for oblique-slip and dip-slip dislocation models, Bull. Seism. Soc. Am., 73, 45-57.
Anderson, J.G., P. Bodin, J. N. Brune, J. Prince, S.K. Singh, R. Quaas, and M. Onate
(1986), Strong ground motion from the Michoacan, Mexico, earthquake, Science,
233, 1043-1049.
Atakan, K., and J. Havskov (1996), Local site effects in the northern North Sea based
on single-station spectral ratios of OBS recordings, Terra Nova, 8, 22-33.
ATC-21 (1988), Rapid Visual Screening of Buildings for Potential Seismic Hazards:
A Handbook, Applied Technology Council, Redwood city, CA, USA.
ATC-40 (1996), Seismic Evaluation and Retrofit of Concrete Buildings, Vol. 1, Applied
Technology Council, Redwood city, CA, USA.
Athanasopoulos, G.A., P.C. Pelekis, and E.A. Leonidou (1999), Effects of surface
topography on seismic ground response in the Egion (Greece) 15 June 1995
earthquake, Soil Dyn. Earthq. Eng., 18, 135-149.
Atkinson, G.M., and D.M. Boore (1995), New ground motion relations for eastern
North America, Bull. Seis. Soc. Am., 85, 17-30.
Atkinson, G.M., and D.M. Boore (1997), Some comparisons between recent ground
motion relations, Seis. Res. Lett., 68, 24-40.
Auden, J.B. (1959), Earthquakes in relation to the Damodar Valley Project, Proc. 1st
Symp. Earthq. Eng. University of Roorkee, Roorkee.
Aviles, C.A., C.H. Scholz, and J. Boatwright (1987), Fractal Analysis applied to
Characteristic segments of the San Andreas fault, J. Geophys. Res., 92, 331-344.
Aviles, J., and L.E. Perez-Rocha (1998), Site effects and soil-structure interaction in
the Valley of Mexico, Soil Dyn. Earthq. Eng., 17, 29-39.
Avouac, J., and P. Tapponnier (1993), Kinetic model of active deformation in central
Asia, Geophys. Res. Lett., 20, 895-898.
Bakir, S, H. Sucuoglu, and T. Yilmaz (2002), An overview of local site effects and the
associated building damage in Adapazari during the 17 August 1999 Izmit earthquake,
Bull. Seis. Soc. Am., 91, 509-526.
164 Seismic Microzonation Atlas
Balassanian, S.Y., A.H. Martirosyan, S.N. Nazaretian, A.R. Arakelian, A.S. Avanessian
V.A. Igumnov, and E. Ruttener (1998), Seismic Hazard Assessment in Armenia, Nat.
Haz., 18, 227-236.
Banerjee S.K. (1957), Earthquakes in the Himalayan Region, Ind. Assoc. for the
Cultivation of Sci., Calcutta.
Bapat, A., R.C. Kulkarni, and S.K. Guha (1983), Catalogue of earthquakes in India
and neighbourhood from historical period up to 1979, Indian Soc. Earthq. Tech.,
Roorkee, India.
Bard, P.Y. (1994), Effects of surface geology on ground motion: recent results and
remaining issues, Proc. 10th Europ. Conf. on Earthq. Eng., G. Duma (Ed.), Vienna,
305–323.
Bard, P.Y. (1997), Local effects on strong ground motion: basic physical phenomena
and estimation methods for microzoning studies, SERINA-Seismic Risk: An Integrated
Seismological, Geotechnical and Structural Approach, ITSAK (Ed.), Thessaloniki.
Bard, P.Y. (2000), Local effects on strong ground motion: Basic physical phenomena
and estimation methods for microzonation studies, Lecture notes.
Bard, P.Y., and F.J. Chávez-García (1993), On the decoupling of surficial sediments
from surrounding geology at Mexico City, Bull. Seism. Soc. Am., 83, 1979-1991.
Bath, M., and S.J. Duda (1964), Earthquake volume, Fault plane area, seismic energy,
strain, deformation and related quantities, Ann. Geofis., 17, 353-368.
Ben-Menahem, A., E Aboodi, and R. Schild (1974), Source of the Great Assam
Earthquake and Interplate Wedge Motion, Phys. Earth and Planet. Int., 9, 265-289.
Beresnev, I.A., K.L. Wen, and Y.T. Yeh (1995), Seismological evidence for nonlinear
elastic ground behavior during large earthquakes, Soil Dyn. Earthq. Eng., 14, 103-
114.
Berge-Thierry, C., P. Bernard, and A. Herrero (2001), Simulating strong ground motion
with the ‘k’ kinematic source model: An application to the seismic hazard in the
Erzincan basin, Turkey, J. Seismol., 5, 85-101.
Bhatia, S.C., M.R. Kumar, and H.K. Gupta (1999), A probabilistic seismic hazard
map of India and adjoining regions, Ann. Geofis., 42, 1153–1164.
Bibliography 165
Bhattacharya, P.M., and Kayal, J.R. (2003), Mapping the b-value and its correlation
with the fractal dimension in the Northeast region of India, J. Geol. Soc. Ind., 62,
680–695.
Bilham, R. (1995), Location and magnitude of the 1833 Nepal earthquake and its
relation to the rupture zones of contiguous Great Himalayan Earthquakes, Curr.
Sci., 69, 155-187.
Bilham, R. (1998), Slip parameters for the Rann of Kachchh, India, 16 June 1819,
earthquake, quantified from contemporary accounts, Coastal Tectonics, I. S. Stewart
and C. Vita-Finzi (Eds.), 146, 295-318
Bilham, R., and P. England (2001), Plateau pop-up during the Great 1897 Assam
earthquake, Nature, 410, 806-809.
Bilham, R., and V.K. Gaur (2000), The geodetic contribution to the study of
seismotectonics in India, Curr. Sci., 79, 1259-1269.
Bilham, R., F. Blume, R. Bendick, and V.K. Gaur (1998), The geodetic constraints
on the translation and deformation of India: implications for future great Himalayan
earthquakes, Curr. Sci., 74, 213-229.
Bilham, R., R. Bendick, and K. Wallace, (2003), Flexure of the Indian Plate and
intraplate earthquakes, Proc. Indian Acad. Sci. (Earth Planet Sci.),112,1-14 .
Bilham, R., V.K. Gaur, and P. Molnar (2001), Himalayan Seismic Hazard, Science,
293, 1442-1444.
Boatwright, J. (1978), Detailed spectral analysis of two small New York State
earthquakes, Bull. Seism. Soc. Am., 68, 1117-1131.
166 Seismic Microzonation Atlas
Boatwright, J., L.C. Seekins, T.E. Fumal, Liu Hsi-Ping, and C.S. Mueller (1991a),
Ground motion amplification in the Marina district, Bull. Seism. Soc. Am., 81, 1980-
1997.
Boatwright, J., J.B. Fletcher, and T.E. Fumal (1991b), A general inversion scheme for
source, site and propagation characteristics using multiply recorded sets of moderate-
sized earthquakes, Bull. Seism. Soc. Am., 81, 1754-1782.
Bolt, B.A. (1993), Earthquake and Geological discovery, Scientific American Library,
New York.
Bommer, J.J. (2003), Uncertainty about the Uncertainty in Seismic Hazard Analysis,
Eng. Geol., 70, 165-168.
Boore, D.M., W.B. Joyner, and T.E. Fumal (1997), Equations for estimating horizontal
response spectra and peak acceleration from western North American earthquakes:
A summary of recent work, Seism. Res. Lett., 68, 128-153.
Boore, D.M. (1973), The Pacoima Dam accelerogram of the February 9, 1971, San
Fernando earthquake implications of local site topography, Bull. Seism. Soc. Am.,
63, 1603-1609.
Boore, D.M., W.B. Joyner, and T.E. Fumal (1993), Estimation of response spectra
and peak ground accelerations from western North America earthquakes: an interim
report, U. S. Geol. Surv. Open-File report, 93-509.
Boore, D.M., W.B. Joyner, and T.E. Fumal (1994), Estimation of response spectra
and peak ground accelerations from western North America earthquakes: an interim
report, part 2, U. S. Geol. Surv. Open-File report, 94-127.
Borcherdt, R.D. (1970), Effects of local geology on ground motion near San Francisco
Bay, Bull. Seism. Soc. Am., 60, 29-61.
Bordet, P. (1973), On the position of the Himalayan Main Central Thrust within Nepal
Himalayas, paper presented in Seminar on Geodynamics of the Himalayan Region,
Hyderabad, 140-155.
Bouchon, M., and J. S. Barker (1996), Seismic response of a hill: The Example of
Tarzana, California, Bull. Seism. Soc. Am., 86, 66-72.
Bouckovalas, G.D. and G.P. Kouretzis (2001), Stiff soil amplification effects in the 7
September 1999 Athens (Greece) earthquake, Soil Dyn. Earthq. Eng., 21, 671-687.
Bouhadad, Y., and N. Laouami (2002), Earthquake Hazard Assessment in the Oran
Region (Northwest Algeria), Nat. Haz., 26, 227-243.
Brune, J.N. (1970), Tectonic Stress and the spectra of seismic shear waves from
earthquakes, J. Geophys. Res., 75, 4997-5009.
Bullore M.D. (1904), The Seismic Phenomena in British India and Their Correction
with Its Geology, Mem. Geol. Surv. Ind., p. 3.
Bune, V.I., N.A. Vvedenskaya, I.V. Gorbunova, N.V. Kondorskaya, N.S. Landyreva,
and I.V. Fedorova (1970), Correlation of MLH and mpv by Data of the Network of
Seismic Stations of the USSR, Geophy. J. Int., 19, 533–542.
Campbell, K.W. (1990), Empirical prediction of near-source soil and soft-rock ground
motion for the Diablo Canyon power plant site, San Luis Obispo County, California,
Report Prepared for Lawrence Livermore National Laboratory, Dames and Moore,
Evergreen, Colorado.
Campbell, K.W. (1992), Recommended models for predicting strong ground motion
and MMI in the San Francisco bay area, Report Prepared for EQE International,
Dames and Moore, Evergreen, Colorado.
Campbell, K.W. (1993), Empirical prediction of near-source ground motion from large
earthquakes, Proceedings, International Workshop on Earthquake Hazard and Large
Dames in the Himalaya, Jan, 15-16, 1993, New Delhi, V.K. Gaur (Ed.), The Indian
National Trust for Art and Cultural Heritage (INTACH), New Delhi, 93-103.
Campbell, K.W. (2003), Prediction of Strong Ground Motion using the Hybrid Empirical
Method and its use in the development of Ground-Motion (Attenuation) relations in
eastern North America, Bull. Seism. Soc. Am., 93, 1012-1033.
Castro, R.R., F. Pacor, A. Sala, and C. Petrungaro (1996), S-wave attenuation and
site effects in the region of Friuli, Italy, J. Geophys. Res., 101, 22355-22369.
Castro, R.R., J.G. Anderson, and S.K. Singh (1990), Site response, attenuation and
source spectra of S-waves along the Guerrero, Mexico, subduction zone, Bull. Seism.
Soc. Am., 80, 1481-1503.
Chandler, A.M., L.S. Chan, and N.T.K. Lam (2001), Deterministic seismic hazard
parameters and engineering risk implications for the Hong Kong region, J. Asian
Earth Sci., 20, 59-72.
Chauhury H.M., and H.N. Srivastava (1976), Seismicity and Focal Mechanism of
Some Recent Earthquakes in Northeast India, Ann. Geoph., 29-41.
Chavez-García, F.J., L.R. Sanchez, and D. Hatzfeld (1996), Topographic site effects
and HSVR: A comparison between observations and theory, Bull. Seism. Soc. Am.,
86, 1559–1573.
Chin, B., and K. Aki (1991), Simultaneous study of the Source, Path and Site effects
on strong ground motion during the 1989 Loma Prieta earthquake: A preliminary
result on pervasive nonlinear site effects, Bull. Seism. Soc. Am., 81, 1859-1884.
Coppersmith, K.J., P.C. Thenhaus, J.E. Ebel, W.K. Wedge, J.H. Williams, and N.C.
Hester (1993), Regional Seismotectonic Settings, Seismic Hazard Assessment in
the central and eastern United states, Monograph 1: Hazard Assessment, S.T.
Algermissen and G.A Bollinger (Eds.), Central United States Earthquake Consortium,
Memphis, TN, 35-80.
Cornell, C.A. (1968), Engineering Seismic Risk Analysis, Bull. Seism. Soc. Am.,
58, 1583-1606.
Costa, G., G.F. Panza, P. Suhadolc, and F. Vaccari (1993), Zoning of the Italian
territory in terms of expected peak ground acceleration derived from complete synthetic
seismograms, Geophysical Exploration in Areas of Complex Geology, R. Cassinis,
K. Helbig, and G.F. Panza (Eds.), J. Appl. Geophys., 30, 149–160.
Curry, J.R., D.G. Moore, L.A. Lawver, F.J. Emmel, R.W. Raitt, M. Henry, and R.
Kieckhefer (1979), Tectonics of the Andaman Sea and Burma, Am. Assoc. Petr.
Geol. Mem., 29, 189-198.
Dahle, A., H. Bungum, and L.G. Dvamme (1990), Attenuation models inferred from
intraplate earthquake recordings, Earthq. Eng. Struct. Dyn., 19, 1125-1141.
De, R., and J.R. Kayal (2003), Seismotectonic Model of the Sikkim Himalaya:
Constraint from microearthquake surveys, Bull. Seism. Soc. Am., 91, 1395-1400.
DeMets, C., R.G. Gordon, D.F. Argus, and S. Stein (1994), Effect of recent version to
the geomagnetism reversal time scale on estimates of current plate motions, Geophys.
Res. Lett., 21, 2191-2194.
Downes, G.L. (1995), Atlas of Isoseismal Maps of New Zealand Earthquakes, Institute
of Geological and Nuclear Sciences, Lower Hutt.
Bibliography 171
Duan, B., D.D. Oglesby and S.K. Park (2003), A new optimization scheme for site
response estimation, Geophys. Res. Abstracts, 5, p. 11296.
Dutta, T.K. (1964), Seismicity of Assam - Zones of Tectonic Activity, Bull. N.G.R.I., 2,
p.152.
Dutta, T.K., and M.M. Saikia (1976) The Eastern Limit of the Himalayan Orogenic
Belt - the Indo-Burman Orogen and its Geodynamic Development. Himalayan Geology,
6, p. 303.
Duval, A.M., J.P, Meneroud, S. Vidal, and P.Y. Bard (1995), Usefulness of microtremor
measurements for site effect studies, Proc. 10th Euro. Conf. on Earthq. Engn., Duma
(Ed.), Rotterdam, 521-527.
EERI (1995), The Hyogo-ken Nanbu Earthquake January 17, 1995 Preliminary
Reconnaissance Report, Earthquake Engineering Research Institute, Oakland, Calif.
El-Sayed, A., F. Vaccari, and G.F. Panza (2001), Deterministic seismic hazard in
Egypt, Geophys. J. Int., 144, 555-567.
Epstein B., and C. Lomnitz (1966), A Model for the Occurrence of Large Earthquakes,
Nature, 211 954.
Erdik, M., Y.A. Biro, T. Onur, K. Sesetyan, and G. Birgoren (1999), Assessment of
earthquake hazard in Turkey and neighboring regions, Annali di Geofisica, 42, 1125-
1138.
172 Seismic Microzonation Atlas
Esteva, L., C. Lomnitz, and E. Rosenblueth (1958), Seismic Risk and Engineering
decision, Elsevier, New York.
Everdeen J.F. (1970), Study of Regional Seismicity and Associated Problem, Bull.
Siesm. Soc. Am., 60, 393.
Field, E.H., and K.H. Jacob (1993), The theoretical response of sedimentary layers
to ambient seismic noise, Geophys. Res. Lett., 20, 2925-2928.
Field, E.H., and K.H. Jacob (1995), A comparison and test of various site response
estimation techniques, including three that are not reference site dependent, Bull.
Seism. Soc. Am., 85, 1127 - 1143.
Field, E.H., A.C. Clement, V. Aharonian, P.A. Friberg, L. Carroll, T.O. Babaian, S.S.
Karapetian, S.M. Hovanessian, and H.A. Abramian (1995a), Earthquake site-response
study in Giumri (formerly Leninakan), Armenia, using ambient noise observations,
Bull. Seism. Soc. Am., 85, 349-353.
Field, E.H., J.G. Anderson, T.L. Henyey, D.D. Jackson, W.B. Joyner, Y. Lee, H.
Magistrale, B. Minster, K.B. Olsen, J.H. Steidl, L.A. Wald, and C.J. Wills (2000),
Accounting for site effects in probabilistic seismic hazard analyses of southern
California: Overview of the SCEC Phase III Report, Bull. Seism. Soc. Am., 90, S1 -
S31.
Field, E.H., K.H. Jacob, and S.E. Hough (1992), Earthquake Site Response
Estimation: A weak motion case study, Bull. Seism. Soc. Am., 82, 2283-2307.
Bibliography 173
Field, E.H., P.A. Johnson, I.A. Beresnev, and Y. Zeng (1997), Nonlinear ground-
motion amplification by sediments during the 1994 Northridge earthquake, Nature,
390, 599-602.
Field, E.H., S.E. Hough, K.H. Jacob, and P.A. Friberg (1995), Site response in
Oakland, California, near the failed section of the Nimitz freeway
Intern. Jour. of Rock Mechanics and Mining Sciences and Geomechanics Abstracts,
32, 361A-362A.
Finn, W.D., and A.M. Nichols (1988), Seismic response of long-period sites: lessons
from the September 19, 1985, Mexican earthquake, Can. Geotech. J., 25, 28-137.
Fitch, T.J. (1970), Earthquake Mechanism in the Himalayan, Burmese and Anadaman
Regions and Continental Tectonics of Central Asia., J. Geophys. Res., 75, 2699.
Fitch, T.J. (1972), Plate convergance, transcurrent faults and internal deformation
adjacent to southeast Asia and the western pacific, Geophys. Res., 77, 4432-4460.
Frankel, A. (1995a), Mapping Seismic Hazard in the central and eastern United
States, Bull. Seism. Soc. Am., 66, 8-21.
Frankel, A. (1996), National Seismic Hazard Maps: Documentation June 1996, U.S.
Geol. Surv. Open-File Rep., 96-532.
Fukushima, Y., and T. Tanaka (1990), A new attenuation relation for peak horizontal
ground acceleration of strong ground motion in Japan, Bull. Seism. Soc. Am., 80,
757-783.
174 Seismic Microzonation Atlas
Fukushima, Y., and T. Tanaka (1991), A new attenuation relation for peak horizontal
acceleration of strong earthquake ground motion in Japan, Shimizu Tech. Res. Bull.,
10, 1-11.
Fumal, T.E., and J.C. Tinsley (1985), Mapping shear-wave velocities of near-surface
geologic materials. Evaluating Earthquake Hazards in the Los Angeles Region —
An Earth-Science Perspective, J.I. Ziony (Ed.), US Geological Survey Prof. Paper,
1360, 127-149
Gahalaut, V.K., and R. Chander (1997), On interseismic elevation changes and strain
accumulation for great thrust earthquakes in the Nepal Himalaya, Geophys. Res.
Lett., 24, 1011-1014.
Gansser, A. (1964), Geology of the Himalayas, Interscience Publ., New York, p. 289.
Gaur, V.K., and R.K.S. Chouhan (1968), Quantitative measures of seismicity applied
to Indian regions, Bull. Ind. Soc. Earthq. Tech., 5, 63–78.
Gayskiy, Y.N., and A.P. Katok (1965), Primeneniye te0rii ekstremalnykh znazcheniy
dlya ocenki povtoryaemosti silnykh zemletryaseniy. Sb. Dinamika zemnoy kory AN
SSSR, Nauka, M. 9.
Geli, L., P.Y. Bard, and B. Jullien (1988), The effect of topography on earthquake
ground motion: a review and new results, Bull. Seism. Soc. Am., 78, 42-63.
Gibowicz, S.J., and S. Lasocki (2001), Seismicity induced by mining: ten years
later, Adv. Geophy., R. Dmoska, and B. Saltzmann (Eds.), 44, 39-181.
GSI (1939), The Bihar-Nepal earthquake of 1934, Geol. Surv. India Mem., 7, p. 280.
GSI (1993), Bihar-Nepal earthquake, August 20, 1988, in Geol. Surv. India, Sp. Pub.,
D.R. Nandy, A.K. Choudhury, C. Chakraborty, and P.L. Narula (Eds.), 31, p. 104.
GSI (2000), Seismotectonic Atlas of India and its Environs, P.L. Narula, S.K. Acharyya,
and J. Banerjee (Eds.), Geo. Surv. Ind. Spl. Pub., Kolkata.
Bibliography 175
Gueguen, P., J.L. Chatelain, B. Guillier, H. Yepes, and J. Egred (1998), Site effect
and damage distribution in Pujili (Ecuador) after the 28 March 1996 earthquake, Soil
Dyn. Earthq. Eng., 17, 329-334.
Guha, S.K. (1962), Seismic regionalization of India, Proc. 2nd Symp. Earthq. Eng.,
Roorkee, 191–207.
Gupta, I.D., and M.D. Trifunac (1988), Attenuation of intensity with epicentral distance
in India, Soil Dyn. Earthq. Eng., 7, 162-169.
Gupta, S.C., and A. Kumar (2002), Seismic wave attenuation characteristics of three
Indian regions: A comparative study, Curr. Sci., 82, 407-413.
Gutenberg, B., and C.F. Richter (1944), Frequency of earthquakes in California, Bull.
Seism. Soc. Am., 34, 185–188.
Gutenberg, B., and C.F. Richter (1954), Seismicity of Earth and Associated
Phenomena, Princeton University Press, Princeton, New Jersey.
Gutierrez, C., and S.K. Singh (1992), A site effect study in Acapulco, Guerrero,
Mexico: comparison of results from strong-motion and microtremor data, Bull. Seism.
Soc. Am., 82, 642-659.
Gyi, L.T. (1973), Seismic Zoning Map of the Union of Burma, Bull. Int. Inst. Seis. and
Earth. Eng., 9, p. 218.
Hanna, K.C., and R.B. Culpepper (1998), GIS and Site Design: New Tools for Design
Professionals, John Wiley & Sons Ltd., New York, USA, p. 240.
Harkrider, D.G. (1976), Potentials and Displacements for Two Theoretical Seismic
Sources, Geophy. J. R. Astr. Soc., 47, 97–133.
Hartzell, S.A., A. Leeds, A. Frankel, and J. Michael (1996), Site response for urban
Los Angeles using aftershocks of the Northridge earthquake, Bull. Seism. Soc. Am.,
86, 168-192.
176 Seismic Microzonation Atlas
Hartzell, S.A., D. Carver, and R.A. Williams (2001), Site response, shallow shear-
wave velocity and damage in Los Gatos, California, from the 1989 Loma Prieta
earthquake, Bull. Seism. Soc. Am., 91, 468-478.
Hartzell, S.H. (1992), Site response estimation from earthquake data, Bull. Seism.
Soc. Am., 82, 2308-2327.
Haskell, N.A. (1960), Crustal reflection of plane SH waves, J. Geophys. Res., 65,
4147-4150.
Haskell, N.A. (1963), Radiation pattern of Rayleigh waves from a fault of arbitrary dip
and direction of motion in a homogeneous medium, Bull. Seis. Soc. Am., 53, 619 -
642.
Haskell, N.A. (1964), Total energy and energy spectral density of elastic wave radiation
from propagating faults, Bull. Seis. Soc. Am., 54, 1811-1841
Havskov, J., and L. Ottemoeller (2000), SEISAN earthquake analysis software, Seis.
Res. lett., 70, 532-534.
Heaton, T.H., F. Tajimaand, and A.W. Mori (1986), Estimating ground motions using
recorded accelerograms, Surv. Geophys., 8, 25-83.
Helmberger, D.V. (1968), The crust-mantle transition in the Bering Sea, Bull. Seism.
Soc. Am., 58, 179-214.
Hong, P.N. (2001), Probabilistic Seismic Hazard Assessment along the Southeastern
Coast of Vietnam, Nat. Haz., 24, 53-74.
Hough, S.E., R.D, Borcherdt, P.A. Friberg, R. Busby, E. Filed, and K.H. Jacob (1990),
The role of sediment-induced amplification in the collapse of the Nimitz freeway
during the October 17, 1989 Loma Prieta earthquake, Nature, 344, 853-855.
Hwang, H., and J.R. Huo (1997), Attenuation relations of ground motion for rock and
soil sites in eastern United States, Soil Dyn. Earthq. Eng., 16, 363-372.
IMD (2000), Chamoli earthquake of March 29, 1999 and its aftershocks, Meteorological
Monograph, India Meteorological Department, Seismology No. 2/2000.
Iyengar, R.N. (2000), Seismic Status of Delhi Megacity, Curr. Sci., 78, 568–574.
Iyengar, R.N., and S.T.G. Raghukanth (2005), Strong Ground Motion Estimation during
the Kutch, India Earthquake, PAGEOPH (In Press).
Johnston, A.C., K.J. Coppersmith, L.R. Kanter, and C.A. Cornell (1994), The
earthquakes of stable continent interiors, Vol. 1, Assessment of large Earthquake
Potential, Electric Power Research Institute, Palo Alto, CA.
Joshi, A. (2003), Predicting strong motion parameters for the Chamoli earthquake of
28th March 1999, Garhwal Himalaya, India, from simplified finite fault model, J. Seis.,
7, 1-17.
Joyner, W.B., and D.M. Boore (1981), Peak horizontal acceleration and velocity from
strong motion records including records from the 1979 Imperial Valley, California,
earthquake, Bull. Seism. Soc. Am., 71, 2011-2038.
Joyner, W.B., and D.M. Boore (1988), Measurement, characteristics and prediction
of strong ground motion, Proc. Earthq. Eng. Soil Dyn. II, Geotechnical Div/ASCE,
Park City, Utah, 43-102.
Kagan, Y., and L. Knopoff (1976), Statistical search for non-random features of the
seismicity of strong earthquakes, Phys. Earth Planet. Interiors, 12, 291-318.
Kagan, Y.Y. (1991), Seismic moment distribution, Geophys. J. Int., 106, 123-134.
Kagan, Y.Y. (1993), Statistics of characteristics earthquakes, Bull. Seism. Soc. Am,
83, 7-24
Kagan, Y.Y. (1997), Statistical aspects of Parkfield earthquake sequence and Parkfield
prediction experiment, Tectonophysics, 270, 207-219.
Kaila, K.L., and D. Sarkar (1977), Earthquake intensity attenuation pattern in India,
Proc. Symp. Analysis of Seismicity and Seismic Risk, 173-191, Liblice,
Czechoslovakia.
Kaila, K.L., and M. Rao (1979), Seismic zoning maps of Indian subcontinent, Geophys.
Res. Bull.,17, 293–301.
Kanamori, H., and D.L. Anderson (1975), Theoretical basis of some empirical relations
in seismology, Bull. Seism. Soc. Am., 65, 1073-1095.
Karnik V. (1972), A Note on the Morphology and Activity of Seismic Zones in the
Aegean Region, Nature Phys. Sc., 239, p. 72.
Bibliography 179
Karnik V., and Z. Schenkova (1974), Application of the Theory of Largest Values to
Earthquake Occurrence in the Balkan Region, Studia geoph, et geod., 18, p. 134.
Karnik, V., and C. Radu (1977), Time Variation of Earthquake Activity in the Balkan
Region, Studia geoph, et geod, 21, p. 285.
Kato, K., K. Aki, and M. Takemura (1995), Site amplification from coda waves:
Validation and Application to S-wave Site response, Bull. Seism. Soc. Am., 85, 467-
477.
Kayabali K. (1996), Soil liquefaction evaluation using shear wave velocity, Eng. Geology,
44, 121-127.
Kayal, J.R., S.S. Arefiev, S. Barua, D. Hazarika, N. Gogoi, A. Kumar, S.N. Chowdhury,
and S. Kalita (2006), Shillong plateau earthquakes in northeast India region: complex
tectonic model, Curr. Sci., 91, 109-114.
Kayal, J.R. (1996), Precursor seismicity, foreshocks and aftershocks of the Uttarkashi
earthquake of October 20, 1991 at the Garhwal Himalaya, Tectonophys., 263, 339-
345.
Kayal, J.R. (2001), Microearthquake activity in some parts of the Himalaya and the
tectonic model, Tectonophys., 339, 331-351.
Kayal, J.R., O.P. Singh, P.K. Chakraborty, and G. Karunakar (2002), Aftershock
sequence of the Chamoli earthquake of March 1999 in the Garhwal Himalaya, Bull.
Seism. Soc. Am., 93, 109-117.
Kebede, F., and T. van Eck (1997), Probabilistic seismic hazard assessment for the
Horn of Africa based on seismotectonic regionalisation, Tectonophys., 270, 221-237.
Keffer (1984) Landslides caused by earthquakes, Geol. Soc. Am. Bull., 95,
406 – 421.
180 Seismic Microzonation Atlas
Khattri, K., M. Wyss, V.K. Gaur, S.N. Saha, and V.K. Bansal (1985), Seismic Gap
in Assam, Northeast India, Bull. Seis. Soc. Am., 75, 1469-1470.
Khattri, N., and M. Wyss (1978), Precursory Variation of Seismicity Rate in the
Assam Area, Indi. Geol., 6, 615.
Khattri, K.N. (1987), Great earthquakes, seismicity gaps and potential for earthquake
disaster along the Himalaya plate boundary. Earthquake Prediction, K. Mogi and
K.N. Khattri (Eds.), Tectonophys., 138, 79-92.
Khattri, K.N. (1999), An evaluation of earthquakes hazard and risk in northern India,
Him. Geo., 20, 1-46.
Khattri, K.N., A.M. Rogers, D.M. Perkins, and S.T. Algermissen (1984). A seismic
hazard map of India and adjacent areas, Tectonophys., 108, 93-134.
Khattri, K.N., and M. Weiss (1978), Precursory variation of seismic rate in Assam
Area, India, Geology, 6, 685-688.
Kijko, A., and B. J. Mitchell (1983), Multimode Rayleigh wave attenuation and Q in
the crust of the Barents shelf, J. Geophys. Res., 88, 3315-3328.
King, J.L., and B.E. Tucker (1984), Observed variations of earthquake motion across
a sediment-filled valley, Bull. Seism. Soc. Am., 74, 137-151.
Knopoff, L., and Kagan Y. (1977), Analysis of the theory of extremes as applied to
earthquake problems, Geophys. Res., 82, 5647-5657
Kobayashi, K. (1980), A method for presuming deep ground soil structures by means
of longer period microtremors, Proc. of the 7th WCEE, Sept. 8-13, Istanbul, Turkey,
1, 237-240.
Konno, T., K. Kato, and R.L. Nigbor (1999), Soil site response, Nuclear Eng. Design,
192, 285-294.
Bibliography 181
Koravos, G.Ch., I.G. Main, T.M. Tsapanos, and R.M.W. Musson (2003), Maximum
earthquake magnitudes in the Aegean area constrained by tectonic moment release
rates, Geophys. J. Int., 152, 94-112.
Korte, G.B. (1997), The GIS Book: Understanding the Value and Implementation of
Geographic Information Systems, OnWord Press, Albany, New York, p. 440.
Krinitzsky, E.L. (2003), How to combine deterministic and probabilistic methods for
assessing earthquake hazards, Eng. Geology, 70, 157-163.
Lachet, C., and P.Y. Bard (1994), Numerical and theoretical investigations on the
possibilities and limitations of Nakamura’s technique, J. Phy. Earth, 42, 377-397.
Lachet, C., M. Bouchon, N. Theodulidis, and P.Y. Bard (1995), Horizontal to vertical
ratio and geological conditions, Proc. 10th Eur. Conf. Earthq. Eng., Duma (Ed.),
Rotterdam, 285-289.
Langston, C.A. (1979), Structure under Mount Rainier, Washington, inferred from
teleseismic body waves, J. Geophys. Res., 84, 4749 - 4762.
Lapajne, J.K., B.S. Motnikar, and P. Zupancic (1996), Preliminary Seismic Hazard
Maps of Slovenia, Nat. Haz., 14, 155-164.
182 Seismic Microzonation Atlas
Legrand, D., A. Cisternas, and L. Dorbath (1996), Multifractal analysis of the 1992
Erzincan aftershock sequence, Geophy. Res. Lett., 23, 933-936.
Lermo, J., and F.J. Chávez-García (1992), Site effect evaluation using microtremors:
a review (abstract), EOS, 73, 352.
Lermo, J., and F.J. Chávez-García (1993), Site effect evaluation using spectral ratios
with only one station, Bull. Seism. Soc. Am., 83, 5, 1574-1594.
Lermo, J., and F.J. Chávez-García (1994), Are microtremors useful in site response
evaluation?, Bull. Seism. Soc. Am., 84, 1350-1364.
Leydecker, G., and J.R. Kopera (1999), Seismological hazard assessment for a site
in Northern Germany, an area of low seismicity, Eng. Geology, 52, 293-304.
Lindley, G.T., and R.J. Archuleta (1995), Variation of seismic site effects in the Santa
Cruz Mountains, California, Int. Jour. Rock Mech. Mining Sce. Geomech.
Abstracts, 32, 362A-362A.
Longley, P., and M. Batty (1997), Spatial Analysis: Modelling in a GIS Environment,
John Wiley & Sons Ltd., New York, USA, p. 392.
Lyon-Caen, H., and P. Molner (1985), Gravity anomalies, flexure of the Indian plate
and structure, support and evolution of the Himalaya and Ganga basin, Tectonics, 4,
513-538.
Lyubushin, A.A., T.M. Tsapanos, V.F. Pisarenko, and G.C. Koravos (2002), Seismic
Hazard for Selected Sites in Greece: A Bayesian Estimate of Seismic Peak Ground
Acceleration, Nat. Haz., 25, 83-98.
Main, I.G. (1996), Statistical physics, Seismogenesis and seismic hazard, Rev.
Geophys., 34, 433-462
Main, I.G., D. Irving, R. Musson, and A. Reading (1999), Constraints on the frequency-
magnitude relation and maximum magnitudes in the UK from observed seismicity
and glacio-isostatic recovery rates, Geophys. J. Int., 137, 535-550.
Bibliography 183
Makropoulos, K.C., and P.W. Burton (1985a), Seismic Hazard in Greece I. Magnitude
Recurrence, Tectonophys., 117, 105-257.
Makropoulos, K.C., and P.W. Burton (1985b), Seismic Hazard in Greece. II. Ground
Acceleration, Tectonophys., 117, 259-294.
Mäntyniemi, P., V. Mârza, A. Kijko, and P. Retief (2003), A New Probabilistic Seismic
Hazard Analysis for the Vrancea (Romania) Seismogenic Zone, Nat. Haz., 29, 371-
385.
Marble, D.F., and D.J. Pequet (1983), Geographic Information System and Remote
Sensing, Manual. Rem. Sens. Am. Soc., Photogrammetry, Falls Church, VA.
Marrara, F., and P. Suhadolc (1998), Site amplifications in the city of Benevento
(Italy): comparison of observed and estimated ground motion from explosive sources,
J. Seism., 2, 125-143.
Matsuoka, M., and S. Midorikawa (1994), The digital national land information and
seismic microzoning, 22nd Symp. Earthq. Ground Motion, 23-34 (in Japanese with
English abstract).
Mc Harg, I. (1969), Design with Nature, National History Press, New York.
McKenzie, D.P., and J.G. Sclatcr (1971), Evolution of the Indian Ocean since the
Late Cretaceous, Geoph. J. R. Astr. Soc., 24, p. 437.
Milne, W. A., and A.G. Davenport (1969), Distribution of Earthquake Risk in Canada,
Bull. Seism. Soc. Am., 59, 729.
Mithal, R.S., and L.S. Srivastava (1959), Geotectonic position and earthquakes of
Ganga-Brahmputra region, Proc. 1st Symp. Earthquake Eng., Roorkee.
184 Seismic Microzonation Atlas
Mitra, S., K. Priestley, V.K. Gaur, S.S. Rai, and J. Haines (2006), Variation of Rayleigh
wave group velocity dispersion and seismic heterogeneity of the Indian crust and
uppermost mantle, Geophys. J. Int., 164, 88–98
Molnar, P., and P. Taponnier (1975), Conozoic tectonics of Asia, Effects of continental
collision, Science, 189, 419-426.
Molnar, P., and P. Tapponnier (1977), Relation of tectonics of eastern Chiena to the
India-Eurasia collision, application of slipline field theory to large scale continental
tectonics, Geology, 5, 212-216.
Molnar, P., and W.P. Chen (1982), Seismicity and mountain building, Mountain Building
Processes, K. Hsu (Ed.), Academic Press, New York, 41-57.
Molnar, P., T.J. Fitch, and F.T. Wu (1973), Fault Plane Solutions of Shallow
Earthquakes and Contemporary Tectonics, Asia. Earth Planet. Sci. Lett., 19, 101.
Muço, B., F. Vaccari, G. Panza, and N. Kuka (2002), Seismic zonation in Albania
using a deterministic approach, Tectonophys., 344, 277-288.
Musson, R.M.W., and P.W. Winter (1996), Seismic Hazard Maps for the U.K., Nat.
Haz., 14, 141-154.
Bibliography 185
Nandy, D.R. (1973), Geology and Structural Lineaments of the Luhit Himalaya
(Arunachal Pradesh) and Adjoining Area; A Tectonic Interpretation, Seminar on
Geodynamics of the Himalayan Region, Hyderabad.
Nath, S.K. (2004), Seismic Hazard Mapping and Microzonation in the Sikkim Himalaya
through GIS Integration of Site Effects and Strong Ground Motion Attributes, Nat.
Haz., 31, 319-342.
Nath, S.K., M. Vyas, I. Pal, and P. Sengupta (2004), A Hazard Scenario in the
Sikkim Himalaya from Seismotectonics, Spectral Amplification, Source
Parameterization and Spectral Attenuation Laws using Strong Motion Seismometry,
J. Geophys. Res., 110, B01301, doi: 10.1029/2004/2004JB003199.
Nath, S.K., M. Vyas, I. Pal., and P. Sengupta (2005), A Seismic Hazard Scenario In
The Sikkim Himalaya from Seismotectonics, Spectral Amplification, Source
Parameterization, and Spectral Attenuation Laws Using Strong Motion Seismometry,
J. Geophys. Res., 110, B01301, doi: 10.1029/2004jb003199.
Nath, S.K., P. Sengupta, and J.R. Kayal (2002), Determination of Site Response at
Garhwal Himalaya from the aftershock sequence of 1999 Chamoli Earthquake, Bull.
Seism. Soc. Am., 92, 1071-1081.
186 Seismic Microzonation Atlas
National Bureau of Soil Survey (1994), Master Plan for Irrigation Development, Sikkim,
Govt. of Sikkim, Agricultural Finance Corporation Ltd, Report, 2, 35-47.
Nuannin, P., O. Kulhanek, and L. Persson (2005), Spatial and temporal b value
anomalies preceding the devastating off coast of NW Sumatra earthquake of December
26, 2004, Geophys. Res. Lett., 32, L11307, doi: 10. 1029/2005GL022679.
Oldham, R.D. (1882), A catalogue of Indian earthquakes from the earliest time to the
end of A.D. 1869, by the late Thomas Oldham edited by R. D. Oldham, Mem. Geol.
Surv. India, 19, 163-215.
Oldham, R.D. (1899), Report on the great earthquake of 12th June 1897 (incl. the
reports by P. Bose, G.Grimes, H.Hayden, T. LaTouche and E. Vredenburg), Mem.
Geol. Survey of India, 29, p.1379, Calcutta.
Oldham, R.D. (1899), Report on the Great Earthquake of 12th June, 1897, Mem.
Geol. Surv. Ind., 29, p. 1379, Calcutta.
Ordaz, M., and S.K. Singh (1992), Source spectra attenuation of seismic waves
from Mexican earthquakes and evidence of amplification in the hill zone of Mexico
City, Bull. Seism. Soc. Am., 82, 24-43.
Bibliography 187
Panza, G.F., F. Vaccari, G. Costa, P. Suhadolc, and D. Faeh (1996), Seismic input
modelling for zoning and microzoning, Earthq. Spectra, 12, 529–566.
Parolai, S., D. Bindi, and P. Augliera (2000), Application of the generalized inversion
technique (GIT) to a microzonation study: Numerical simulations and comparison
with different site-estimation techniques, Bull. Seism. Soc. Am., 90, 286-297.
Parvez, I.A., A.A. Gusev, G.F. Panza, and A.G. Petukhin (2001), Preliminary
determination of the interdependence among strong motion amplitude, earthquake
magnitude and hypocentral distance for the Himalayan region, Geophys. J. Int., 144,
577–596.
Parvez, I.A., F. Vaccari, and G.F. Panza (2003), A deterministic seismic hazard map
of India and adjacent areas, Geophys. J. Int., 155, 489-508.
Parvez, I.A., F. Vaccari, and G.F. Panza (2004), Site-specific microzonation study in
Delhi metropolitan city by 2-D modeling of SH and P-SV waves, Pure Appl. Geophys.,
161, 1165-1184.
Parvez, I.A., G.F. Panza, A.A. Gusev, and F. Vaccari (2002), Strong-motion Amplitude
in Himalayas and a Pilot Study for the Deterministic First-order Microzonation in a
Part of Delhi City, Curr. Sci., 82, 158–166.
188 Seismic Microzonation Atlas
Pedersen, H., B. Le. Brun, D. Hatzfeld, M. Campillo, and P.Y. Bard (1994), Ground-
Motion amplitude across ridges, Bull. Seism. Soc. Am., 84, 1786-1800.
Pendse, C.G. (1948), Earthquakes in India and Neighbourhood, India Met. Dept. Sci.
Notes, 10, p.177.
Petersen, M.D., W.A. Bryant, C.H. Cramer, M.S. Reichle, and C.R. Real (1997),
Seismic ground-motion hazard mapping incorporating site effects for Los Angeles,
Orange and Ventura counties, California: a geographical information system
application. Bull. Seism. Soc. Am., 87, 249-255.
Phillips, W.S., and K. Aki (1986), Site amplification of coda waves from local
earthquakes in central California, Bull. Seism. Soc. Am., 76, 627-648.
Ponce, V.M. (1989), Engineering Hydrology, Principles and Practices, Prentice Hall,
Englewood Cliffs, New Jersey.
Porter, L.D., and D.J. Leeds (2000), Correlation of strong ground motion for the 1994
Northridge earthquake, Geophys. J. Int.,143, 376-388.
Radulian, M., F. Vaccari, N. Mândrescu, G.F. Panza, and C.L. Moloveanu (2000),
Seismic hazard of Romania: deterministic approach, PAGEOPH,157, 221–247.
Raina, V.K., and B.S. Srivastava (1980), A reappraisal of the geology of the Sikkim
lesser Himalaya, Stratigraphy and Correlation of lesser Himalayan formation, K.S.
Valdiya and S. B. Bhatia (Eds.), 201-220, Hindustan Publishing Corporation, Delhi.
Rajendran, K., C.P. Rajendran, S.K. Jain, C.V.R. Murty, and J.N. Arlekar (2000), The
Chamoli earthquake, Garhwal Himalaya: Field observations and implications for
seismic hazard, Curr. Sci., 78, 45-51.
Raju, A.T.R. (1968), Geological Evaluation of Assam and Cambay Basins of India,
Bull. Am. Assoc. Petrol., 52, 2422.
Rastogi, B.K. (2000), Chamoli earthquake of magnitude 6.6 on 29th March 1999, J.
Geol. Soc. India, 55, 505-514.
Reasenberg, P.A. and L.M. Jones (1989), Earthquake hazard after a mainshock in
California, Science, 243, 1173-1176.
Richter, C.F. (1958), Elementary Seismology, San Francisco, W.H. Freeman and
Co., p. 345.
Riepl, J., P.Y. Bard, D. Hatzfeld, C. Papaioannou, and S. Nechtschein (1998), Detailed
evaluation of site response estimation methods across and the sedimentary valley of
Volvi (EURO-SEISTEST), Bull. Seism. Soc. Am., 88, 488-502.
Romanelli, F., and F. Vaccari (1999), Site response estimation and ground motion
spectrum scenario in the Catania area, J. Seismol., 3, 311-326.
Romeo, R., and A. Pugliese (2000), Seismicity, seismotectonics and seismic hazard
of Italy, Eng. Geology, 55, 241-266.
Saaty, T.L. (1980), The Analytic Hierarchy Process, McGraw-Hill International, New
York, NY, USA.
Saaty, T.L. (1990), An Exposition of the AHP in Reply to the Paper ‘Remarks on the
AnalyticHierarchy Process, Management Science, 36, 259-268.
Sadigh, K., C.Y. Chang, N.A. Abrahamson, S.J. Chiou, and M.S. Power (1993),
Specification of Long-Period Ground Motions: Updated Attenuation Relationships for
Rock Site Conditions and Adjustment Factors for Near-Fault Effects, In Proc. ATC-
17-1 Seminar on Seismic Isolation, Passive Energy Dissipation and Active Control,
Vol. I, March 11-12, San Francisco, Applied technology Council, Redwood City, CA,
11-23.
Sadigh, K., C.Y. Chang, Z.A. Egan, F. Makdisi, and R.R. Youngs (1997), Attenuation
Relationships for Shallow Crustal Earthquakes Based on California Strong Motion
Data, Seism. Res. Lett., 68, 180-189.
Safak, E. (1995), A new model to simulate site effects, Proc. 10th European Conference
on Earthquake Engineering, Duma (Ed.), Rotterdam, 297-303.
Safak, E. (2001), Local site effects and dynamic soil behavior, Soil Dyn. Earthq.
Eng., 21, 453-458.
Sander, P. (1998), Hard rock aquifers in arid and semi arid zones, Water Resources
of Hard Rock Aquifers in Arid and (Semi-Arid) Zones-Studies and Reports in Hydrology,
J.W. Lloyd (Ed.), UNESCO Publishing, Paris, France, 58, p. 224
Schenkova, Z., and V. Karnik (1978), The Third Asymptotic Distribution of Largest
Magnitudes in the Balkan Earthquake Provinces, Pure Appl. Geoph., 116, 1314-
1325
Schenkova, Z., and V. Schenk (1975), Return Periods of Earthquakes and Trends of
Seismic Activity, Pure Appl. Geoph., 113, p 683.
Schnabel, B., J. Lysmer, and H.B. Seed (1972), SHAKE – A Computer Program for
Earthquake Response Analysis of Horizontally Layered Sites, Report EERC, p. 7212.
Schnabel, P.B., (1973), Effects of Local Geology and Distance from Source on
Earthquake Ground Motion, Ph.D. Thesis, University of California, Berkeley, California.
Seale, S.H., and R.J. Archuleta (1989), Site amplification and attenuation of strong
ground motion, Bull. Seis. Soc. Am., 79, 1673-1696.
Seed, H.B., and I.M. Idriss (1970), Soil Moduli and Damping Factors for Dynamic
Response Analyses, Rep. No. EERC-70/10, Earthquake Engineering Research Center,
University of California, Berkeley, California
Seekins, L.C., L. Wennerberg, L. Margheriti, and H.P. Liu (1996). Site amplification
at five locations in San Francisco, California; a comparison of S-waves, codas and
microtremors, Bull. Seism. Soc. Am., 86, 627-635.
Sen, Z., and K. Al-Suba’I (2001), Seismic Hazard Assessment in the Tihamat Asir
Region, Southwestern Saudi Arabia, Math. Geology, 33, 967-991.
Shankar, R., and P. Narula (1999), Chamoli earthquake of 29th March 1999- A
preliminary appraisal, Ind. Minerals, 52, 141-150.
Sharma, M.L (1998). Attenuation relationship for estimation of peak ground horizontal
acceleration using data from strong motion arrays in India, Bull. Seism. Soc. Am.,
88, 1063-1069.
Shi, Y., and B.A. Bolt (1982), The standard error of the magnitude-frequency b value,
Bull. Seism. Soc. Am., 72, 1677-1687.
Bibliography 191
Shiono, K., Y. Ohta, and K. Kudo (1979), Observation of 1 to 5 sec microtremors and
their applications to earthquake engineering, Part VI: existence of Rayleigh wave
components, Jour. Seism. Soc. Japan, 32, 115-124 (in Japanese with English
abstract).
Shoji, Y., and M. Kamiyama (2002), Estimation of local site effects by generalized
inversion scheme using observed records of ‘Small-Titan’, Soil Dyn. Earthquake Eng.,
22, 855-864.
Singh, S.K., E. Mena, and R. Castro (1988), Some aspects of source characteristics
of the 19 September 1985 Michoacan earthquake and ground motion amplification in
and near Mexico City from strong motion data, Bull. Seism. Soc. Am., 78, 451-477.
Singh, S.K., W.K. Mohanty, B.K. Bansal, and G.S. Roonwal (2002), Ground motion
in Delhi from future large/great earthquakes in the central seismic gap of the Himalayan
Arc, Bull. Seism. Soc. Am., 92, 555-569.
Sousa, M.L., and C.S. Oliveira (1996), Hazard Mapping Based on Macroseismic
Data Considering the Influence of Geological Conditions, Nat.Haz., 14, 207-226.
Speidel, D.H., and P.H. Mattson (1997), Problems for Probabilistic Seismic Hazard
Analysis, Nat. Haz., 16, 165-179.
Spudich, P., W.B. Joyner, A.G. Lindh, D.M. Boore, B.M. Margaris, and J.B. Fletcher
(1999), SEA99: A revised ground motion prediction relation for use in extensional
tectonic regimes, Bull. Seism. Soc. Am., 89, 1156-1170.
192 Seismic Microzonation Atlas
Star, S., K.C. McGwire, and J.E. Estes (1997), Integration of Geographic Information
Systems and Remote Sensing. Cambridge University Press, Cambridge, UK, p.
248.
Sykes, L.R., E.S. Bruce, and C.H. Scholz (1999), Rethinking Earthquake Prediction,
Pure Appl. Geophys., 155, 207–232.
Talebian, M. (1996), Seismic risk analysis of Japan: peak ground acceleration and
uniform risk response spectra, Internl. Jour. of Rock Mechanics and Mining Sciences
and Geomechanics Abstracts, 33, 218A-218A.
Tandon, A.N. (1954), Study of Great Assam Earthquake of August 15, 1950 and Its
Aftershocks, Ind. J. Meteorol. Geophys., 5, 95-137.
Tandon, A.N. (1956), Zones of India liable to earthquake damage, Ind. J. Meteorol.
Geophys.,10, 137–146.
Tandon, A.N., and H.N. Srivastava (1975), Focal Mechanism of Some Recent
Himalayan Earthquakes and Regional Plate Tectonics, Bull. Seis. Soc. Am., 65,
963-969.
Teves-Costa, P., L. Matias, and P.Y. Bard (1996), Seismic behavior estimation of thin
alluvium layers using microtremor recordings, Soil Dyn. Earthquake Eng.,15, 201-
209.
Than, K. (1975), Earthquake Mechanisms Pressure and Tension Axes in Indian Ocean
and Related Regions, Bull. Int. Inst. Seis. and Earth. Eng., 11, p. 1.
Bibliography 193
Theodulidis, N.P., and B.C. Papazachos (1992), Dependence of strong ground motion
on magnitude-distance, site geology and macroseismic intensity for shallow
earthquakes in Greece: I, Peak horizontal acceleration, velocity and displacement,
Soil Dyn. Earthq. Eng., 11, 387-402.
Thingbaijam, K.K.S., S.K. Nath, A. Yadav, A. Raj, M.Y. Walling, and W.K. Mohanty
(2007), Recent Seismicity in Northeast India and its Adjoining Region, J. Seism.,
Communicated.
Tonouchi, K., T. Sakayama, and T. Imai (1983), S wave velocity in the ground and the
damping factor, Bull. Int. Assoc. Eng. Geologists, 26-27, 327-333.
Toro, G.R, N.A. Abrahamson, and J.F. Schneider (1997), Model of strong ground
motions from earthquakes in central and eastern North America: Best estimates and
uncertainties, Seism. Res. Lett., 68, 41-57.
Toro, G.R., W.J. Silva, R.K. McGuire, and R.B. Herrmann (1992), Probabilistic seismic
hazard mapping of the Mississippi embayment, Seism. Res. Lett., 63, 449-475.
Trifunac, M.D. (1976), Preliminary empirical model for scaling Fourier amplitude
spectra of strong ground acceleration in terms of earthquake magnitude, source-to
station distance and recording site condition, Bull. Seism. Soc. Am., 66, 1343-1373.
Trifunac, M.D., and M.I. Todorovska (2000), Can aftershock studies predict site
amplification factors? Northridge, CA, Earthquake of 17 January 1994, Soil Dyn.
Earthq. Eng.,19, 233–251.
Tsapanos, T.M. (2003), A seismic hazard scenario for the main cities of Crete island,
Greece, Geophys. J. Int., 153, 403-408.
194 Seismic Microzonation Atlas
Tsapanos, T.M., and P.W. Burton (1991), Seismic hazard evaluation for specific seismic
regions of the world, Tectonophys., 194, 153-169.
Tumarkin, A.G., and R.J. Archuleta (1997), Recent Advances in Prediction and
Processing of Strong Ground Motions, Nat. Haz., 15, 199-215.
Usami, T. (1981), Nihon Higai Jishin Soran (List of Damaging Japanese Earthquakes),
University of Tokyo Press (in Japanese).
Utsu, T. (1965), A method for determining the value of b in the formula log N=a-bM
showing the magnitude–frequency relation for the earthquakes. Geophy. Bull.
Hokkaido Univ., 13, 99-103 (in Japanese with English abstract).
Vorobieva, I.A. (1999), Prediction of a subsequent large earthquake, Phy. Earth Planet.
Interiors, 111, 197-206.
Vorobieva, I.A., and G.F. Panza (1993), Prediction of the occurrence of related strong
earthquakes in Italy, Pure Appl. Geophys., 141, 25-41.
Vulnerability Atlas of India (1999), Building Materials & Technology Promotion Council,
New Delhi.
Vyas, M., S.K. Nath, I. Pal, P. Sengupta and W.K. Mohanty (2004), GSHAP Revisited
for the Prediction of Maximum Credible Earthquake in the Sikkim Region, ACTA
Geophys. Polo., 53, 143-152.
Wang, C.Y., and R.B. Herrmann (1980), A numerical study of P-, SV- and SH wave
generation in a plane layered medium, Bull. Seism. Soc. Am., 70, 1015-1036.
Bibliography 195
Warren, N.W., and G.V. Latham (1970), An Experimental Study of Thermally Induced
Microfracturing and its Relation to Volcanic Seismicity, J. Geophys Res., 75, 4455-
4464.
Wells, D.L., and K.L. Coppersmith (1994), New empirical relationships among
magnitude, rupture length, rupture width, rupture area and surface displacement,
Bull. Seism. Soc. Am., 84, 974-1002
Wood, H.O. (1908), Distribution of apparent intensity in San Francisco in The California
Earthquake of April 18, 1906, Report of the State Earthquake Investigation
Commission, Carnegie Inst. of Wash. Publ., 87, Washington, D.C., 220-245.
Yegulalp, T. M., and J.T. Kuo (1974), Statistical prediction of the occurrence of
maximum magnitude earthquakes, Bull. Seism. Soc. Am., 64, 393-414.
Yilmaztürk, A., and P.W. Burton (1999), An evaluation of seismic hazard parameters
in southern Turkey, J. Seismol., 3, 61-81.
Youngs, R. R., S.J. Chiou, W.J. Silva, and J.R. Humphrey (1997), Strong ground
motion attenuation relationships for subduction zone earthquakes, Seism. Res. Lett.,
68, 58-73.
Zivciæ, P., P. Suhadolc, and F. Vaccari (2000), Seismic zonation of Slovenia based
on deterministic hazard computation, Pure Appl. Geophys., 157, 171–184.
ANNEXURE I
Sample SPT
Description Depth Strata Collected N- Ncr Remarks
(m) D/S U/S 15cm 15cm 15cm Value
(m) (m)
16.00
Coarse sandy clay with 16.50 22 28 * refusal 28
mix silt from 16m to 17.5m 17.00
19.00
19.50 26 * * refusal 37
20.00
Blue sandy clay with
mix silt from 17.5 m to 21.00 21.00 24 * * refusal 36
29.0 m
22.00
22.50 34 * * refusal 35
23.00
25.00
25.50 32 * * refusal 32
26.00
28.00
28.50 30 * * refusal 30
29.00
Fine sand with mix silt
30.00
294
Soil test borehole chart
Sample SPT
Description Depth Strata Collected N- Ncr Remarks
(m) D/S U/S 15cm 15cm 15cm Value
(m) (m)
0.00
Filled up soil upto 2.0m
1.00
1.50 0 1 1 2 3
Bluish Clay from 2m to 2.00
2.8m
3.00 3.00 1 2 3 5 7
Greyish brown clay from
2.8m to 4.6m 4.00
4.50 3 5 6 11 14
Redish grey clay from 5.00 5.00
4.6m to 5.0m
6.00 6.00 3 5 7 12 14
7.00
Greyish brown clay from 7.50 2 3 6 9 9
5.0m to 10.0m 8.00
9.00 9.00 3 3 5 8 8
10.00 10.00
greyish brown clay with 10.50 3 5 7 12 12
mix silt from 11.00
10.0m to 11.8m
12.00 12.00 3 6 10 16 15
Brownish silty clay from
11.8 m to 15.3 m 13.00
13.50 4 7 9 16 14
14.00
15.00 15.00 4 9 23 32 21
295
Soil test borehole chart
Sample SPT
Description Depth Strata Collected N- Ncr Remarks
(m) D/S U/S 15cm 15cm 15cm Value
(m) (m)
16.00
Brownish fine sand with 16.50 8 12 20 32 20
mix silt from 17.00
15.3m to 25.0m
18.00 18.00 8 12 23 35 21
19.00
19.50 9 11 22 33 20
20.00
22.00
22.50 20 28 * refusal 25
23.00
25.00
25.50 29 * * refusal 24
Brownish coarse sand with 26.00
mix silt from 25.0m to 27.0m
27.00 27.00 32 * * refusal 23
28.00
Fine sand with mix silt 28.50 30 * * refusal 23
from 27.0m. To 30.0m 29.00
30.00
296
Soil test borehole chart
Sample SPT
Description Depth Strata Collected N- Ncr Remarks
(m) D/S U/S 15cm 15cm 15cm Value
(m) (m)
0.00
3.00 3.00 2 4 6 10 17
Bluish / redish sandy clay
from 3.3 to 5.6 m 4.00
4.50 6 8 10 18 26
5.00
6.00 6.00 4 6 8 14 19
Brownish/ bluish silty clay
from 5.6m to 12.0m 7.00
7.50 4 7 8 15 19
8.00
9.00 9.00 7 12 14 26 31
10.00
10.50 8 12 13 25 28
11.00
12.00 12.00 4 8 11 19 21
15.00 15.00 6 9 14 23 23
297
Soil test borehole chart
Sample SPT
Description Depth Strata Collected N- Ncr Remarks
(m) D/S U/S 15cm 15cm 15cm Value
(m) (m)
16.00
Bluish silty clay from 15.5m 16.50 7 11 17 28 27
to 20.0m 17.00
18.00 18.00 12 21 27 49 45
19.00
19.50 18 27 * refusal 44
20.00
22.00
Sand from 20.0m to 26.8 m 22.50 22 30 * refusal 28
23.00
25.00
25.50 12 38 * refusal 26
26.00
27.00 27.00 10 15 22 37 27
30.00
298
Soil test borehole chart
Sample SPT
Description Depth Strata Collected 15cm 15cm 15cm N- Ncr Remarks
(m) D/S U/S Value
(m) (m)
0.00
Filled up soil upto 1.0m
1.00
Greyish blue clay from 1.50 1 2 2 4 7
from 1.0m to 3.15 m 2.00
3.00 3.00 2 8 8 16 24
Greyish blue sandy clay
from 3.15m to 4.65 m 4.00
4.50 3 1 2 3 4
5.00
Bluish sandy clay with mix 5.50
silt from 4.65m to 10.0m 6.00 6.00 0 0 0 0 0
7.00
7.50 0 1 1 2 2
8.00
9.00 9.00 2 3 4 7 7
10.00
Bluish silty clay from 10.0m 10.50 3 3 5 8 8
to 12.15 m 11.00
12.00 12.00 3 8 16 24 19
15.00 15.00 9 7 10 17 15
299
Soil test borehole chart
Sample SPT
Description Depth Strata Collected 15cm 15cm 15cm N- Ncr Remarks
(m) D/S U/S Value
(m) (m)
19.00
Fine sand with mix silt from 19.50 13 18 33+ refusal 27
19.2m to 20.3m 20.00
22.00
Brown coarse sand from 22.50 15 25 25+ refusal 26
20.3m to 26.0m 23.00
25.00
25.50 20 31 19+ refusal 24
26.00
gravel found
couldnot be
28.50 50 blows of
29.00 DCPT at
40cm.
30.00
300
Soil test borehole chart
Sample SPT
Description Depth Strata Collected N- Ncr Remarks
(m) D/S U/S 15cm 15cm 15cm Value
(m) (m)
0.00
Bluish / Redish clay upto
5.5m 1.00
1.50 1 2 5 7 12
2.00
3.00 3.00 4 9 9 18 26
3.50
4.00
4.50 5 9 10 19 25
5.00
6.00 6.00 8 9 12 21 26
9.00 9.00 4 7 9 16 17
10.00
10.50 6 8 11 19 18
11.00
Brownish clayey sand from
9.0m to 22.0m 12.00 12.00 8 10 11 21 19
13.00
13.50 8 12 14 26 21
14.00
301
Soil test borehole chart
Sample SPT
Description Depth Strata Collected N- Ncr Remarks
(m) D/S U/S 15cm 15cm 15cm Value
(m) (m)
16.00
16.50 15 25 25+ refusal 30
17.00
19.00
19.50 20 31 19+ refusal 28
20.00
22.00
22.50 25 35 15+ refusal 27
23.00
DCPT
29.00 performed.
30.00 30.00
302
ANNEXURE IV
Physical and Shear Parameters of Sediment as Obtained from 200 Boreholes
Bh_No. Depth From To Bulk Wet Saturated Dry SpGr Void Porosity% LL% PL% PI C Kg/cm2 φ
(m) (m) (m) Den Den Den Den R
1 2.00 1.97 1.99 1.57 2.63 0.67 40.00 42.80 23.62 19.18 0.420 7.0
3.50 1.98 1.98 1.59 2.62 0.65 39.00 44.00 22.87 21.13 0.430 6.5
5.00 1.98 1.99 1.61 2.63 0.63 39.00 43.60 21.56 22.04 0.450 7.0
8.00 2.00 1.99 1.62 2.61 0.61 38.00 47.20 22.48 24.72 0.440 6.0
12.00 1.99 2.03 1.66 2.62 0.58 37.00 40.70 23.51 19.19 0.360 10.0
15.00 2.00 2.00 1.62 2.62 0.62 38.00 45.90 22.63 23.27 0.450 7
Bh_No. Depth From To Bulk Wet Saturated Dry SpGr Void Porosity% LL% PL% PI C Kg/cm2 φ
(m) (m) Den Den Den Den R
Bh_No. Depth From To Bulk Wet Saturated Dry SpGr Void Porosity% LL% PL% PI C Kg/cm2 φ
(m) (m) (m) Den Den Den Den R
4 NO DATA
Bh_No. Depth From To Bulk Wet Saturated Dry SpGr Void Porosity% LL% PL% PI C Kg/cm2 φ
(m) (m) (m) Den Den Den Den R
5 0.00 1.00 2.12 1.60 2.60 0.62 35.08 22.51 12.57 0.65 0
1.00 3.50 2.00 1.55 2.60 0.67 34.70 21.95 12.75 0.62 0
3.50 4.50 1.96 1.50 2.60 0.73 34.10 21.60 12.50 0.58 0
7.00 9.00 2.07 1.65 2.60 0.58 32.70 21.33 11.37 0.98 0
Bh_No. Depth From To Bulk Wet Saturated Dry SpGr Void Porosity% LL% PL% PI C Kg/cm2 φ
(m) (m) (m) Den Den Den Den R
6 3.5 1.87 1.45 2.60 0.793 44.00 46.30 23.75 22.55 0.20 7.5
11.0 1.96 1.59 2.60 0.635 39.00 48.80 23.67 24.63 0.48 7.0
14.0 1.97 1.60 2.61 0.631 39.00 48.30 23.67 24.63 0.50 6.5
Bh_No. Depth From To Bulk Wet Saturated Dry SpGr Void Porosity% LL% PL% PI C Kg/cm2 φ
(m) (m) (m) Den Den Den Den R
7 1.5 1.83 1.86 1.39 2.61 0.877 47.00 43.80 24.68 19.12 0.14 7.0
3.0 1.85 1.88 1.43 2.62 0.832 45.00 41.50 0.13 8.5
5.0 1.84 1.87 1.41 2.60 0.844 46.00 45.30 23.87 21.43 0.16 5.0
8.0 1.83 1.87 1.42 2.61 0.838 45.00 48.70 22.95 25.75 0.15 6.0
10.5 2.01 2.14 1.83 2.67 0.459 31.00 0.08 28.0
13.5 2.03 2.05 1.71 2.60 0.520 34.00 47.60 23.22 24.3 0.48 7.0
Annexure IV
Bh_No. Depth From To Bulk Wet Saturated Dry SpGr Void Porosity% LL% PL% PI C Kg/cm2 φ
(m) (m) (m) Den Den Den Den R
1.50 1.95 42 20 22
3.00 3.45 42 20 22
4.50 4.95 42 20 22
6.00 6.45 42 20 22
7.50 7.95 46 22 24
9.00 9.45 46 22 24
10.50 10.95 46 22 24
12.00 12.45 46 22 24
305
ANNEXURE V
1 312.6 (2) 308.3 (2) 315.2 (3) 372.8 (2) 355.6 (2)
2 261.0 (2) 192.3 (1) 248.9 (3) 274.4 (3) 269.3 (3)
3 312.6 (2) 266.5 (1) 293 (1) 315.2 (3) 294.8 (3)
4 269.6 (2) 265.3 (2) 248.9 (3) 284.6 (3) 274.4 (3)
5 316.9 (2) 240.0(1) 294.8 (3) 289.7 (3) 279.5 (3)
6 208.2 (1) 248.1 (2) 239.5 (2) 239.5 (2) 279.5 (3)
7 282.4 (1) 286.8 (2) 274.4 (3) 289.7 (3) 274.4 (3)
8 282.5 (2) 213.2 (3) 284.6 (3) 289.7 (3) 274.4 (3)
9 224.1 (1) 312.6 (2) 299.9 (3) 279.5 (3) 269.3 (3)
10 286.8 (2) 305.0 (3) 299.9 (3) 284.6 (3) 269.3 (3)
11 213.5 (1) 259.1 (3) 305 (3) 289.7 (3) 279.5 (3)
12 282.5 (2) 208.1 (3) 238.7 (3) 289.7 (3) *
13 213.5 (1) 315.2 (3) 299.9 (3) 284.6 (3) 274.4 (3)
14 197.6 (1) 289.7 (3) 299.9 (3) 284.6 (3) 269.3 (3)
15 291.1 (2) 312.6 (2) 294.8 (3) 274.4 (3) 269.3 (3)
16 187.0 (1) 228.5 (3) 259.1 (3) 284.6 (3) 279.5 (3)
17 176.4 (1) 176.4 (1) 171.1 (1) 295.4 (2) 279.5 (3)
18 165.8 (1) 176.4 (1) 197.6 (1) 295.4 (2) 279.5 (3)
19 160.5 (1) 239.5 (2) 252.4 (2) 299.9 (3) 274.4 (3)
308 Seismic Microzonation Atlas
20 181.7 (1) 176.4 (1) 187 (1) 295.4 (2) 259.1 (3)
21 155.2 (1) 181.7 (1) 187 (1) 254 (3) 278.2 (2)
22 208.2 (1) 197.6 (1) 304 (2) 294.8 (3) 279.5 (3)
23 261.0 (2) 245.3 (1) 254 (3) 274.4 (3) 274.4 (3)
24 286.8 (2) 240 (1) 277.1 (1) 294.8 (3) 284.6 (3)
25 304.0 (2) 286.8 (2) 202.9 (1) 181.7 (1) 243.8 (3)
26 229.4 (1) 229.4 (1) 289.7 (3) 284.6 (3) 269.3 (3)
27 229.4 (1) 250.6 (1) 289.7 (3) 284.6 (3) 269.3 ( 3)
28 181.7 (1) 208.2 (1) 299.9 (3) 284.6 (3) 274.4 (3)
29 248.1 (2) 197.6 (1) 213.5 (1) 289.7 (3) 279.5 (3)
30 252.4 (2) 213.5 (1) 351.3 (2) 289.7 (3) 274.4 (3)
31 181.7 (1) 187 (1) 213.5 (1) 218.8 (1) 274.4 (3)
32 171.1 (1) 224.1 (1) 176.4 (1) 282.5 (1) 377.1 (3)
33 293 (1) 234.7 (1) 265.3 (1) 213.5 (1) 224.1 (1)
34 218.8 (1) 181.7 (1) 218.8 (1) 208.2 (1) 271.8 (1)
35 208.2 (1) 266.5 (1) 294.8 (3) 279.5 (3) 269.3 (3)
36 250.6 (1) 266.5 (1) 321.2 (2) 208.2 (1) 310.1 (3)
37 171.1 (2) 233.6 (3) 305 (3) 294.8 (3) 279.5 (3)
38 273.9 (2) 266.5 (1) 202.9 (1) 273.9 (2) *
39 139.3 (1) 224.1 (2) 264.2 (2) 291.1 (3) 278.2 (3)
40 149.9 (1) 187 (1) 310.1 (3) 289.7 (3) 274.4 (3)
41 171.1 (1) 165.8 (1) 245.3 (1) 334.1 (2) 284.6 (3)
42 261 (1) 250.6 (2) 265.3 (2) 279.5 (3) 254 (1)
43 202.9 (1) 265.3 (2) 284.6 (3) 284.6 (3) 274.4 (3)
44 213.5 (1) 269.3 (3) 289.7 (3) 279.5 (3) 264.2 (3)
45 171.1 (1) 176.4 (1) 265.3 (2) 208.2 (1) 240 (1)
46 202.9 (1) 294.8 (3) 284.6 (3) 274.4 (3) 264.2 (3)
Annexure V 309
47 176.4 (1) 187 (1) 269.3 (3) 284.6 (3) 274.4 (3)
48 282.5 (2) 304 (2) 224.1 (1) 240 (1) 261.2 (1)
49 197.6 (1) 294.8 (3) 294.3 (3) 289.7 (3) 274.4 (3)
50 266.5 (1) 208.2 (1) 299.9 (3) 289.7 (3) 274.4 (3)
51 245.3 (1) 293 (1) 240 (1) 294.8 (3) 248.9 (3)
52 235.2 (2) 299.7 (2) 279.5 (3) 289.7 (3) 274.4 (3)
53 213.5 (1) 192.3 (1) 347 (2) 284.6 (3) 274.4 (3)
54 218.8 (1) 255.9 (1) 359.9 (2) 250.6 (1) 229.4 (1)
55 192.3 (1) 181.7 (1) 208.2 (1) 305 (3) 274.4 (3)
56 208.2 (1) 271.8 (1) 282.4 (1) 289.7 (3) 274.4 (3)
57 155.2 (1) 240 (1) 274.4 (3) * *
58 256.7 (2) 192.3 (1) 279.5 (3) 284.6 (3) 269.3 (3)
59 176.4 (1) 197.6 (1) 291.1 (2) 269.3 (3) *
60 192.3 (1) 312.6 (2) 240 (1) 255.9 (1) 325.5 (2)
61 187 (1) 192.3 (1) 252.4 (2) 334.1 (2) 274.4 (3)
62 208.2 (1) 266.5 (1) 229.4 (1) 261.2 (1) 286.8 (2)
63 218.8 (1) 261.2 (1) 351.3 (2) 294.8 (3) 279.5 (3)
64 248.1 (2) 273.9 (2) 291.1 (2) 274.4 (3) 243.8 (3)
65 181.7 (1) 273.9 (2) 208.1 (3) 284.6 (3) *
66 176.4 (1) 282.4 (1) 377.1 (2) 359.9 (2) 282.4 (1)
67 252.4 (2) 265.3 (1) 192.3 (1) 218.3 274.4 (3)
68 208.2 (1) 273.9 (2) 289.7 (3) 274.4 (3) 274.4 (3)
69 234.7 (1) 277.1 (1) 240 (1) 338.4 (2) 325.5 (2)
70 278.2 (2) 316.9 (3) 187 (1) 208.2 250.6
71 187 (1) 278.2 (2) 286.8 240 (2) 325.5
72 234.7 (1) 381.4 (2) 304 (2) 250.6 (1) 224.1 (1)
73 197.6(1) 186.8 269.6 (2) 269.6 (2) 284.6 (3)
310 Seismic Microzonation Atlas
74 176.4 (1) 197.6 (1) 228.5 (3) 294.8 (3) 264.2 (3)
75 250.6 (1) 234.7 (1) 271.8 (1) 338.4 (2) 274.4 (3)
76 192.3 (1) 304 (2) 269.6 (2) 295.4 274.4 (3)
77 273.9 (2) 269.6 (2) 335.6 (3) 338.4 (2) 321.2 (2)
78 368.5 (2) 368.5 (2) 396.8 (3) 310.1 (3) *
79 282.5 (2) 213.5 (1) 278.2 (2) 304 (1) *
80 269.6 (2) 325.5 (2) 233.6 (3) 304 (2) 359.9 (2)
81 321.2 (1) 314.2 (1) 299.7 (1) 321.2 (2) 274.4 (3)
82 228.5 (3) 279.5 (3) 299.9 (3) 248.1 (2) 274.4 (3)
83 202.9 (1) 282.4 (1) 299.9 (3) 250.6 (1) 325.5 (2)
84 368.5 (2) 359.9 (2) 320.3 (3) 284.6 (3) 334.1 (2)
85 192.8 (3) 282.5 (2) 202.9 (1) 192.3 (2) 234.7 (1)
86 187(1) 282.5(2) 224.1(1) 274.4(3) 325.5 (2)
87 277.1 (1) 240 (1) 334.1 248.9 (2) 279.5 (3)
88 273.9 (2) 273.9 (2) 224.1 (1) 192.3 (1) 295.4 (2)
89 261 (2) 308.3 (2) 338.4 (2) 277.1 (1) 266.5 (1)
90 248.9 329.8 269.6 (2) 192.8 (3) *
91 187 (1) 202.9 (1) 181.7 (1) 314.2 (1) 279.5 (3)
92 282.5 (2) 282.5 (2) 269.6 (2) 279.5 (3) 294.8 (3)
93 234.7 (1) 240 (1) 295.4 (2) 312.6 (2) 284.6 (3)
94 160.5 (1) 250.6 (1) 243.8 (2) 171.1 (1) 240 (1)
95 197.6 (1) 291.1 (2) 282.5 (2) 266.5 (1) 277.1 (1)
96 187 (1) 202.9 (1) 187 (1) 299.7 (2) 305 (1)
97 202.9 (1) 265.3 (2) 234.7 (1) 330.1 (1) 277.1 (1)
98 213.5 (1) 338.4 (2) 291.1 (2) 330.1 (1) 271.8 (1)
99 224.1 (1) 255.9 (1) 213.5 (1) 271.8 (1) 279.5 (3)
100 250.6 (1) 197.6 (1) 197.6 (1) 289.7 (3) 319.5 (1)
Annexure V 311
128 218.8 (1) 218.8 (1) 321.2 (2) 329.8 (2) 321.2 (2)
129 218.8 (1) 381.4 (2) 299.9 (3) 284.6 (3) 274.4 (3)
130 273.9 (2) 229.4 (1) 274.4 (3) 266.5 (1) 279.5 (3)
131 176.4 (1) 279.5 (3) 304 (2) 284.6 (3) *
132 218.8 (1) 261.2 (1) 299.9 (3) 284.6 (3) 294.8 (3)
133 165.8 (1) 314.2 (1) 305 (3) 289.7 (3) *
134 299.7 (1) 340.7 (3) 305 (3) 284.6 (3) 274.4 (3)
135 240 (1) 325.4 (3) 192.3 (1) 261.2 (1) 279.5 (3)
136 187 (1) 224.1 (1) 293 (1) 289.7 (3) 187 (1)
137 181.7 (1) 303.6 (1) 245.3 (3) * *
138 181.7 (1) 208.2 (1) 310.1 (3) * *
139 197.6 (1) 181.7 (1) 255.9 (1) 197.6 (1) 350.9 (3)
140 278.2 (2) 279.5 (3) 208.2 (1) 278.2 (2) 338.4 (2)
141 282.4 (1) 325.4 (3) 305 (3) 321.2 (2) 274.4 (3)
142 192.3 (1) 197.6 (1) 284.6 (3) 284.6 (3) 427.4 (3)
143 208.2 (1) 312.6 (2) 342.7 (2) 284.6 (3) 279.5 (3)
144 192.3 (1) 261.2 (1) 305 (3) 284.6 (3) 274.4 (3)
145 213.5 (1) 181.7 (1) 259.1 (3) 289.7 (3) 274.4 (3)
146 218.8 (1) 224.1 (1) 289.7 (3) 284.6 (3) *
147 286.8 (2) 299.9 (3) 295.4 (2) 255.9 (1) 250.6 (1)
148 181.7 (1) 171.1 (1) 192.3 (1) 289.7 (3) 279.5 (3)
149 202.9 (1) 218.8 (1) 229.4 (1) 266.5 (1) *
150 295.4 (2) 329.8 (2) 310.1 (3) 289.7 (3) 274.4 (3)
151 208.2 (1) 261.2 (1) 299.9 (3) 289.7 (3) 274.4 (3)
152 202.9 (1) 286.8 (2) 305 (3) 284.6 (3) 269.3 (3)
153 261.2 (1) 312.6 (2) 299.9 (3) 289.7 (3) *
154 165.8 (1) 176.4 (1) 254 (1) 294.8 (3) 279.5 (3)
Annexure V 313
155 208.2 (1) 202.9 (1) 266.5 (1) 284.6 (3) 223.4 (3)
156 197.6 (1) 208.2 (1) 330.5 (3) 335.6 (3) 315.2 (3)
157 187 (1) 213.5 (1) 224.1 (1) 202.9 (1) 245.3 (1)
158 208.2 (1) 278.2 (2) 259.1 (3) 299.9 (3) *
159 213.5 (1) 218.8 (1) 287.7 (1) 294.8 (3) 279.5 (3)
160 176.4 (1) 218.8 (1) 229.4 (1) 218.8 (1) 208.2 (1)
161 229.4 (1) 273.9 (2) 181.7 (1) 213.5 (1) 299.7 (2)
162 187 (1) 229.4 (1) 213.5 (1) 289.7 (3) 279.5 (3)
163 338.4 (2) 424.4 (2) 250.6 (1) 308.3 (1) 295.4 (1)
164 245.3 (1) 286.8 (2) 304 (2) 294.8 (3) 284.6 (3)
165 229.4 (1) 197.6 (1) 255.9 (1) 284.6 (3) 274.4 (3)
166 282.5 (2) 299.7 (2) 310.1 (3) 289.7 (3) 279.5 (3)
167 208.2 (1) 261.2 (1) 315.2 (3) 289.7 (3) 279.5 (3)
168 181.7 (1) 229.4 (1) 316.9 (2) 289.7 (3) 274.4 (3)
169 176.4 (1) 197.6 (1) 238.7 (3) * *
170 181.7 (1) 245.3 (1) 228.5 (3) 224.1 (1) 213.5 (1)
171 240 (1) 286.8 (2) 269.6 (2) 229.4 (1) 291.1 (1)
172 213.5 (1) 187 (1) 187 (1) 277.1 (1) *
173 202.9 (1) 213.5 (1) 261.2 (1) 213.5 (1) 325.5 (2)
174 269.6 (2) 269.6 (2) 286.8 (2) 316.9 (2) 316.9 (2)
175 197.6 (2) 278.2 (2) 299.7 (2) 224.1 (1) 284.6 (3)
176 187 (1) 250.6 (1) 234.7 (1) 213.5 (1) 243.8 (3)
177 269.6 (2) 234.7 (1) 338.4 (2) 289.7 (3) 279.5 (3)
178 187 (1) 208.2 (1) 291.1 (2) 286.8 (2) *
179 192.3 (1) 269.3 (3) 305 (3) 284.6 (3) 274.4 (3)
180 223.4 (3) 245.3 (1) 255.9 (1) 224.1 (1) *
181 224.1 (1) 213.5 (1) 294.8 (3) * *
314 Seismic Microzonation Atlas
182 208.2 (1) 265.3 (2) 295.4 (2) 299.7 (2) 312.6 (2)
183 218.8 (1) 287.7 312.6 (2) 321.2 (2) 368.5 (2)
184 282.4 (1) 372.8 (2) 234.7 (1) 334.4 (1) 315.2 (3)
185 218.8 (1) 286.8 (2) 299.9 (3) 284.6 (3) 279.5 (3)
186 155.2 (1) 192.3 (1) 208.2 (1) 312.3 (2) 329.8 (2)
187 187 (1) 238.7 (3) 310.1 (3) 345.8 (3) 345.8 (3)
188 240 (1) 282.4 (1) 261.2 (1) 305 (3) 254 (3)
189 277.1 (1) 287.7 (1) 277.1 (1) 289.7 (3) 289.7 (3)
190 187 (1) 381.4 (2) 271.8 (1) * *
191 224.1 (1) 334.1 (2) 305 (3) 289.7 (3) *
192 181.7 234.7 304 (2) 271.8 261.2
193 229.4 (1) 202.9 (1) 265.3 (2) 304 (2) 329.8 (2)
194 197.6 (1) 264.2 (3) 310.1 (3) * *
195 165.8 (1) 208.2 (1) 305 (3) 284.6 (3) 274.4 (3)
196 245.3 (1) 320.3 (3) 305 (3) 284.6 (3) 274.4 (3)
197 165.8 (1) 213.5 (1) 312.6 (2) 338.4 (2) 314.2 (1)
198 261.2 (1) 321.2 (2) 250.6 (1) 359.9 (2) 303.6 (1)
199 230.9 269.6 (2) 310.1 (3) * *
200 213.5 334.1 308.3 312.6 (2) 274.4 (3)
Shear wave
Sl. No. Latitude (°N) Longitude (°E) velocity
(Vs30)(m/sec)
Shear wave
Sl. No. Latitude (°N) Longitude (°E) velocity
(Vs30)(m/sec)
Shear wave
Sl. No. Latitude (°N) Longitude (°E) velocity
(Vs30)(m/sec)
Shear wave
Sl. No. Latitude (°N) Longitude (°E) velocity
(Vs30)(m/sec)
Shear wave
Sl. No. Latitude (°N) Longitude (°E) velocity
(Vs30)(m/sec)
Shear wave
Sl. No. Latitude (°N) Longitude (°E) velocity
(Vs30)(m/sec)
Shear wave
Sl. No. Latitude (°N) Longitude (°E) velocity
(Vs30)(m/sec)
1 N 26.14167° E 91.66125° 2.04 2.04 2.04 2.04 1.85 1.85 1.87 1.91 2.20 2.20 2.20 2.20 2.20 2.20 2.25 2.25 2.25 2.30 2.30 2.35
2 N 26.15736° E 91.67229° 1.15 1.30 1.97 2.25 2.29 2.19 2.25 2.25 2.25 2.31 1.99 1.99 1.99 2.25 2.25 2.25 2.25 2.26 2.26 2.30
3 N 26.15273° E 91.65732° 1.50 1.40 1.58 1.58 1.58 1.76 1.76 1.65 1.70 1.71 1.80 1.94 2.20 2.20 2.20 2.25 2.25 1.89 1.91 2.20
4 N 26.14808° E 91.64718° 1.49 1.76 1.72 2.24 2.24 2.24 2.24 1.95 2.04 2.65 2.00 2.00 2.09 2.25 2.25 2.25 2.25 2.25 2.25 2.25
5 N 26.14073° E 91.63915° 1.71 1.71 1.71 1.75 1.75 1.75 1.76 1.60 1.69 2.20 2.21 2.21 2.22 2.22 2.23 2.23 2.23 2.24 2.24 2.24
6 N 26.13753° E 91.62938° 1.14 1.45 1.78 1.58 1.58 1.82 1.82 1.82 1.82 1.82 1.82 1.78 1.78 1.78 1.79 1.79 1.98 2.10 2.20 2.21
7 N 26.12998° E 91.62280° 1.86 1.86 1.81 1.80 1.81 1.82 1.82 1.82 1.60 1.64 1.98 2.20 2.20 2.21 2.22 2.22 2.24 2.24 2.25 2.26
8 N 26.12088° E 91.61617° 1.94 1.94 1.94 1.90 1.42 1.42 1.50 1.64 1.65 1.98 2.20 2.21 2.21 2.22 2.22 2.23 2.23 2.24 2.25 2.25
9 N 26.10162° E 91.61213° 1.80 1.85 1.89 1.85 1.90 1.92 1.95 1.99 2.20 2.21 2.22 2.23 2.23 2.25 2.25 2.25 2.26 2.28 2.30 2.30
10 N 26.09853° E 91.60620° 1.99 1.95 1.95 1.93 1.93 1.99 2.00 2.20 2.21 2.21 2.22 2.23 2.24 2.25 2.25 2.25 2.25 2.26 2.28 2.30
11 N 26.09957° E 91.59685° 1.71 1.71 1.70 1.44 1.51 1.50 1.76 1.82 1.98 2.20 2.20 2.23 2.23 2.23 2.24 2.24 2.25 2.26 * *
12 N 26.10955° E 91.60723° 1.81 1.81 1.84 1.84 1.32 1.32 1.44 1.46 1.51 1.60 1.77 2.20 2.20 2.23 2.23 * * * * *
13 N 26.10744° E 91.59333° 1.35 1.60 1.71 1.81 1.85 1.90 2.20 2.20 2.24 2.24 2.24 2.25 2.25 2.25 2.26 2.26 2.27 * * *
14 N 26.11728° E 91.59867° 1.80 1.82 1.85 1.60 1.74 1.95 2.10 2.20 2.20 2.23 2.23 2.25 2.26 2.26 2.26 2.27 2.28 2.28 2.30 2.30
15 N 26.13875° E 91.60713° 1.82 1.82 1.82 1.82 1.82 1.82 1.72 2.20 2.21 2.22 2.22 2.23 2.23 2.24 2.24 2.25 2.26 2.27 2.28 2.30
16 N 26.12150° E 91.60505° 1.50 1.60 1.60 1.61 1.72 1.70 1.70 1.75 1.80 1.81 1.90 1.90 1.91 2.20 2.20 2.23 2.25 2.26 2.26 2.22
17 N 26.13850° E 91.62132° 1.32 1.86 1.86 1.86 1.76 1.80 1.80 1.80 1.81 1.80 1.80 2.01 2.01 2.01 2.21 2.21 2.23 2.24 2.25 2.25
18 N 26.14743° E 91.63057° 1.81 1.81 1.79 1.79 1.79 1.79 1.79 1.81 1.81 1.81 1.82 1.83 1.85 1.76 2.21 2.21 2.22 2.23 2.23 2.25
19 N 26.15033° E 91.63953° 1.93 1.93 1.90 1.90 1.90 1.87 1.87 1.88 1.88 1.88 1.89 2.20 2.23 2.25 2.26 2.26 2.27 2.28 2.30 2.30
20 N 26.16930° E 91.66760° 2.01 2.01 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.82 1.82 1.82 1.83 1.83 1.76 2.10 2.19 1.92 1.92 2.21
324
Bore DENSITY (gm/c.c)
hole ID Latitude Longitude
21 N 26.16420° E 91.64974° 1.62 1.62 1.66 1.66 1.66 1.66 1.66 1.66 1.67 1.67 1.67 1.68 1.69 1.76 2.21 1.68 1.67 1.68 1.69 1.75
22 N 26.16909° E 91.68080° 1.10 1.81 1.81 1.81 1.81 1.81 1.82 1.50 1.20 1.81 1.81 2.20 2.20 2.23 2.25 2.25 2.26 2.26 2.26 2.27
23 N 26.14820° E 91.66824° 1.81 1.90 1.91 1.91 1.93 1.93 1.93 1.93 1.93 1.93 2.21 2.22 2.22 2.23 2.23 2.24 2.24 2.25 2.25 2.26
24 N 26.15969° E 91.68063° 1.44 1.76 1.76 1.76 1.71 1.79 1.82 1.82 1.78 1.88 1.96 2.20 2.19 2.19 2.20 2.22 2.22 * * *
25 N 26.16035° E 91.68913° 1.13 2.16 1.81 1.86 1.89 1.82 1.83 1.79 1.89 1.89 1.81 1.81 1.79 1.83 1.84 1.83 1.84 1.60 1.96 2.20
26 N 26.14802° E 91.67838° 1.45 1.82 1.90 1.90 1.90 1.90 1.81 1.84 1.90 2.20 2.20 2.23 2.23 2.24 2.26 2.28 2.29 2.30 2.30 2.30
27 N 26.13695° E 91.67352° 1.85 1.85 1.85 1.85 1.85 1.85 1.89 2.20 2.20 2.23 2.24 2.25 2.26 2.27 2.27 2.30 2.30 2.30 2.31 2.31
29 N 26.15994° E 91.71165° 1.82 1.82 1.84 1.90 1.92 1.92 1.83 1.89 1.83 1.83 1.83 2.20 2.22 2.23 2.25 2.26 2.27 * * *
30 N 26.15052° E 91.68468° 1.10 1.67 1.68 1.68 1.68 1.79 1.79 1.80 2.20 2.21 2.22 2.23 2.24 2.24 2.25 2.26 2.26 2.26 2.27 2.27
31 N 26.17373° E 91.72914° 1.13 1.65 1.76 1.76 1.78 1.75 1.81 1.81 1.81 1.75 2.00 2.00 2.00 2.00 2.00 2.20 2.20 2.28 2.25 2.30
32 N 26.17079° E 91.72370° 1.12 1.70 1.70 1.63 1.63 1.63 1.78 1.80 1.80 1.80 1.80 1.80 1.70 1.71 1.71 1.71 2.20 2.20 * *
33 N 26.17092° E 91.76926° 1.97 1.97 1.97 1.97 1.97 1.95 1.95 1.95 1.95 1.50 1.55 1.55 1.57 1.70 1.71 1.71 1.70 1.75 1.76 1.76
34 N 26.16844° E 91.72948° 1.45 1.82 1.85 1.71 1.80 1.71 1.71 1.71 1.71 1.71 1.71 1.71 1.50 1.40 1.51 1.65 2.20 2.20 2.20 2.00
35 N 26.15103° E 91.73013° 1.93 1.93 1.86 1.86 1.86 1.86 2.20 2.20 2.20 2.25 2.26 2.27 2.27 2.28 2.30 2.30 2.30 2.30 * *
36 N 26.18119° E 91.77602° 1.10 1.60 1.60 1.84 1.92 1.96 2.19 1.92 1.96 2.20 1.92 1.95 1.92 1.86 1.92 2.10 2.14 * * *
37 N 26.15991° E 91.73655° 1.20 1.45 1.50 1.50 1.50 1.50 1.79 1.85 1.80 2.20 2.20 2.20 2.25 2.25 2.25 2.25 2.26 2.28 2.30 2.30
38 N 26.11085° E 91.78908° 1.12 1.44 1.46 2.02 2.02 1.85 1.92 1.91 1.75 1.75 1.75 1.75 1.76 1.78 1.46 1.48 * * * *
39 N 26.17793° E 91.77400° 1.98 1.98 1.98 1.98 1.98 1.98 2.10 2.10 2.10 2.20 2.20 1.80 1.80 1.70 1.91 1.72 1.71 1.72 1.73 1.73
40 N 26.16747° E 91.74498° 1.10 1.10 1.76 1.76 1.78 1.78 1.78 1.79 2.20 2.20 2.20 2.20 2.23 2.25 2.26 2.27 2.27 2.28 2.28 2.30
41 N 26.17194° E 91.73618° 1.62 1.62 1.73 1.73 1.65 1.65 1.65 1.65 1.65 1.65 1.65 1.65 1.65 2.20 2.20 2.23 2.23 2.25 2.26 2.26
42 N 26.17378° E 91.77418° 1.82 1.88 1.88 1.88 1.88 1.89 1.87 1.97 1.87 1.85 1.86 1.85 1.88 2.20 2.20 2.24 2.25 1.81 1.82 1.82
43 N 26.17139° E 91.73836° 1.88 1.90 1.90 1.90 1.90 1.90 1.88 1.88 1.88 2.20 2.20 2.23 2.23 2.23 2.25 2.25 2.26 2.26 2.27 2.27
44 N 26.14195° E 91.72871° 1.98 1.98 1.98 1.98 1.72 1.78 2.21 2.21 2.22 2.23 2.24 2.24 2.25 2.26 2.26 2.27 2.27 2.27 2.28 2.28
Annexure VII
Bore DENSITY (gm/c.c)
hole ID Latitude Longitude
45 N 26.17572° E 91.73512° 1.70 1.70 1.70 1.70 1.45 1.45 1.45 1.45 1.51 1.51 1.45 1.53 1.55 1.53 1.52 1.53 1.56 1.61 1.76 1.78
46 N 26.12717° E 91.72135° 1.93 1.93 1.95 1.95 1.77 2.20 2.20 2.20 2.20 2.24 2.24 2.25 2.25 2.25 2.26 2.27 2.28 2.28 2.30 2.30
47 N 26.18357° E 91.73866° 1.88 1.86 1.88 1.89 1.96 1.96 1.93 1.93 1.95 2.20 2.21 2..21 2.21 2.21 2.21 2.23 2.23 2.24 2.25 2.25
48 N 26.16684° E 91.77177° 1.82 1.82 1.83 1.83 1.84 1.83 1.92 1.92 1.92 1.92 1.92 1.92 1.92 1.92 1.95 1.98 2.20 2.20 2.22 2.22
49 N 26.18992° E 91.75242° 1.15 1.84 1.84 1.84 1.84 1.86 1.86 1.60 1.61 1.65 2.20 2.20 2.22 2.23 2.23 2.23 2.24 2.24 2.25 2.25
50 N 26.18618° E 91.74863° 1.77 1.77 1.77 1.78 1.78 1.78 1.78 1.78 1.60 1.95 2.20 2.20 2.20 2.22 2.22 2.23 2.24 2.24 2.24 2.25
51 N 26.16974° E 91.75236° 1.44 1.55 1.60 1.75 1.75 1.73 1.76 1.73 1.73 1.81 1.81 1.81 1.73 2.20 2.20 2.23 2.25 1.80 1.81 1.89
52 N 26.17605° E 91.75524° 1.12 1.90 1.90 1.90 1.90 1.90 1.86 1.86 1.86 1.81 1.81 1.91 2.20 2.20 2.20 2.20 2.25 2.25 2.25 2.30
53 N 26.18457° E 91.74117° 1.81 1.81 1.81 1.81 1.77 1.77 1.77 1.85 2.20 2.20 2.20 2.20 2.20 2.25 2.25 2.25 2.25 2.20 * *
54 N 26.18155° E 91.77837° 1.97 1.97 1.97 1.97 1.97 1.97 1.95 1.95 1.99. 2.20 2.20 1.85 1.80 1.80 1.80 1.85 1.58 1.80 1.85 1.88
55 N 26.18472° E 91.75443° 1.44 1.78 1.78 1.71 1.71 1.71 1.71 1.74 1.79 1.79 1.80 1.81 2.20 2.20 2.20 2.25 2.25 2.25 2.25 2.32
56 N 26.18798° E 91.73660° 1.40 1.88 1.88 1.88 1.88 1.79 1.79 1.75 0.50 1.91 2.20 2.20 2.20 2.20 2.20 2.25 2.25 2.25 2.30 2.30
57 N 26.18906° E 91.79023° 1.12 1.12 1.59 1.59 1.59 1.65 1.65 2.20 2.20 2.20 2.25 * * * * * * * * *
58 N 26.18773° E 91.74272° 1.41 1.70 1.70 1.70 1.70 1.70 2.27 2.27 2.27 2.27 2.27 2.27 2.27 2.25 2.30 2.30 2.30 2.30 2.35 2.35
59 N 26.18668° E 91.76568° 1.12 1.75 1.75 1.70 1.98 2.27 2.27 1.71 1.71 1.91 1.91 1.75 2.20 2.20 2.20 2.25 * * * *
60 N 26.18640° E 91.78783° 1.13 1.92 1.92 1.94 1.94 1.94 1.94 1.94 2.20 1.85 1.85 1.71 1.71 1.71 1.71 1.85 2.20 2.20 2.20 *
61 N 26.16159° E 91.77460° 1.44 2.08 2.08 2.08 2.08 1.85 1.85 1.83 1.85 1.81 1.79 1.79 1.78 2.20 2.20 2.20 2.25 2.25 2.25 2.30
62 N 26.18573° E 91.76258° 1.40 1.71 1.71 1.94 1.93 1.93 1.93 2.00 2.03 2.03 2.03 2.03 2.03 2.09 2.09 2.12 2.00 1.81 2.03 2.20
63 N 26.18482° E 91.77630° 1.89 1.89 1.91 1.90 1.90 2.00 2.00 2.05 2.20 2.20 2.21 2.21 2.21 2.21 2.25 2.25 2.25 2.30 2.30 2.35
64 N 26.15980° E 91.77876° 1.88 1.88 1.89 1.87 1.87 1.87 1.88 1.70 1.77 1.85 1.70 1.73 1.84 2.20 2.20 1.81 1.81 1.82 1.81 1.81
65 N 26.18535° E 91.77317° 1.50 1.70 1.68 1.65 1.69 1.76 1.76 1.76 1.76 1.75 1.76 1.77 1.91 1.92 2.20 2.20 * * * *
66 N 26.18437° E 91.76755° 1.85 2.05 1.60 2.07 2.07 2.07 2.07 2.20 2.20 2.25 2.25 1.68 1.70 2.20 2.20 2.20 2.20 2.20 2.20 2.25
67 N 26.16165° E 91.78170° 1.75 1.77 1.76 1.89 1.89 1.85 1.89 1.89 1.89 1.85 1.85 1.85 1.92 1.89 2.20 2.20 2.20 2.20 2.25 2.25
325
68 N 26.18387° E 91.74662° 1.88 1.88 1.88 1.88 1.84 1.84 1.85 1.89 2.20 2.20 2.20 2.20 2.20 2.25 2.25 2.25 2.30 2.32 2.33 2.35
326
Bore DENSITY (gm/c.c)
hole ID Latitude Longitude
69 N 26.17929° E 91.76470° 1.12 1.19 1.60 2.07 2.08 2.08 1.91 1.74 1.79 1.79 1.72 2.20 2.12 2.11 2.20 2.20 2.20 2.21 2.22 2.22
70 N 26.17230° E 91.77992° 1.40 1.92 1.92 1.92 1.87 1.87 1.87 1.89 1.75 1.87 1.87 1.87 1.83 1.83 1.81 1.83 1.83 2.20 1.87 1.91
71 N 26.16425° E 91.77650° 1.15 1.15 1.90 1.90 1.90 1.90 1.84 1.88 1.88 1.88 1.87 1.87 1.86 1.60 1.97 2.20 2.27 2.27 2.28 2.28
72 N 26.17422° E 91.75918° 1.91 1.91 1.91 1.91 1.91 2.20 2.21 2.20 2.21 1.92 2.12 1.92 1.95 1.96 1.95 2.20 2.21 1.94 1.92 1.94
73 N 26.15813° E 91.78677° 1.80 1.80 1.80 1.79 1.65 1.66 1.61 1.59 1.67 1.65 1.65 1.57 1.65 1.65 1.65 2.20 2.20 2.24 2.26 2.26
74 N 26.14147° E 91.79957° 1.13 1.75 1.81 1.81 1.81 1.93 1.93 1.55 2.17 2.17 1.45 1.44 2.20 2.22 2.23 2.23 2.24 1.81 1.98 2.21
75 N 26.16657° E 91.76322° 1.89 1.95 1.94 1.93 1.92 1.92 1.92 1.91 1.91 1.92 1.91 1.91 1.92 2.20 2.20 2.20 2.21 2.21 2.23 2.25
77 N 26.13491° E 91.80510° 1.68 1.98 1.98 1.98 1.97 1.99 2.09 1.93 1.84 2.20 2.20 2.00 2.07 2.07 2.20 2.20 2.20 2.25 1.92 2.20
78 N 26.17299° E 91.76301° 1.12 1.12 1.93 1.44 1.92 1.56 1.57 1.57 2.20 2.20 2.20 2.20 2.25 2.25 2.25 2.30 * * * *
79 N 26.14380° E 91.78204° 1.70 1.90 1.90 1.85 1.84 1.84 1.84 1.84 1.93 1.93 1.93 1.98 1.98 2.00 * * * * * *
80 N 26.12820° E 91.80667° 1.44 1.90 1.90 1.89 1.92 1.97 1.90 1.97 1.94 1.97 1.79 1.80 1.77 1.79 2.20 2.20 1.96 1.95 1.90 1.91
81 N 26.15613° E 91.78868° 1.80 2.03 2.03 2.03 2.09 2.09 2.20 2.07 2.03 1.97 1.91 1.95 1.77 2.10 2.19 2.20 2.20 2.21 2.22 2.23
82 N 26.13284° E 91.78539° 1.58 1.80 1.99 1.88 1.89 1.89 2.21 2.21 2.21 2.21 2.23 2.23 1.82 1.80 2.07 1.82 1.82 2.20 2.22 2.23
83 N 26.12945° E 91.80169° 1.69 1.71 1.71 1.92 1.92 1.92 1.92 1.92 1.92 2.20 2.21 1.81 1.83 1.83 1.84 2.20 2.21 2.21 2.22 2.23
84 N 26.15262° E 91.78667° 1.12 1.82 1.92 1.95 2.15 1.92 1.81 2.20 2.21 2.20 2.18 2.19 2.20 2.20 2.21 2.21 2.17 1.85 1.82 1.80
85 N 26.12091° E 91.78821° 1.80 1.81 1.80 1.12 1.51 1.51 1.57 1.65 1.65 1.65 1.92 1.92 1.92 1.80 1.81 1.81 1.81 1.81 2.20 2.21
86 N 26.14098° E 91.78667° 1.10 1.68 1.86 1.91 1.91 1.71 1.71 1.71 1.80 1.68 2.20 1.71 1.71 2.20 2.20 2.21 2.21 1.89 * *
87 N 26.15343° E 91.79107° 1.68 1.98 1.98 2.07 1.92 1.97 2.10 2.09 1.91 2.04 2.07 1.91 1.88 1.92 1.81 1.82 2.20 2.20 2.21 2.25
88 N 26.11379° E 91.78621° 1.45 1.68 1.70 1.72 1.72 1.70 1.70 1.72 1.98 1.76 1.71 1.72 1.72 1.72 1.72 1.72 1.71 1.72 2.20 2.21
89 N 26.13844° E 91.77660° 1.44 1.80 1.86 1.86 1.86 1.86 1.89 1.98 2.20 2.12 1.89 1.88 2.20 2.23 2.24 2.25 1.99 * * *
90 N 26.17570° E 91.77977° 1.76 2.05 1.91 2.05 2.10 2.05 2.05 2.11 2.12 2.04 2.04 2.04 2.04 2.04 * * * * * *
91 N 26.12110° E 91.76988° 1.44 1.71 1.82 1.69 1.70 1.75 1.75 1.75 1.81 1.88 1.88 1.88 2.20 2.22 2.23 2.25 2.25 2.26 2.28 2.30
92 N 26.13567° E 91.79333° 2.33 2.33 2.33 2.20 2.02 2.02 2.02 1.94 1.94 1.94 1.94 1.96 1.95 2.20 2.20 2.20 2.25 2.25 2.25 2.35
Annexure VII
Bore DENSITY (gm/c.c)
hole ID Latitude Longitude
93 N 26.13999° E 91.76104° 1.65 1.70 1.78 1.78 1.78 1.78 1.94 1.94 1.69 1.69 1.69 1.69 1.79 1.70 1.80 1.85 2.20 2.20 2.20 2.20
94 N 26.11780° E 91.77639° 1.85 1.96 1.96 1.15 2.20 2.20 1.76 1.75 1.75 1.76 1.71 1.75 1.70 1.75 1.76 1.85 1.86 1.89 2.20 2.20
95 N 26.12552° E 91.78876° 1.54 1.90 1.80 1.80 1.78 1.91 1.92 1.85 1.85 1.80 1.75 1.87 2.20 2.20 2.20 2.20 2.20 2.25 2.25 2.29
96 N 26.14356° E 91.77450° 1.44 1.85 1.84 1.76 1.90 2.08 2.08 1.85 1.84 1.85 1.80 1.79 1.79 1.78 1.78 1.78 1.96 2.25 * *
97 N 26.12259° E 91.78823° 1.75 1.97 1.97 1.96 1.96 1.96 1.96 1.96 1.96 1.98 1.98 2.20 2.20 1.78 1.78 1.76 1.78 1.76 2.20 2.25
98 N 26.12951° E 91.79430° 1.20 1.20 1.86 1.86 1.86 1.87 1.90 1.87 1.87 1.85 1.76 1.78 1.81 2.17 2.17 2.18 1.94 2.20 2.20 2.21
99 N 26.14289° E 91.77021° 1.45 1.77 1.90 1.80 1.81 2.10 2.10 2.10 1.82 1.68 1.82 1.82 1.82 1.89 2.20 2.21 2.20 2.23 2.25 2.25
100 N 26.16052° E 91.75964° 1.83 1.83 1.83 1.83 1.78 1.78 1.78 1.79 1.79 1.79 1.71 1.73 1.70 1.70 2.20 2.20 1.79 1.79 1.79 1.79
101 N 26.13330° E 91.82640° 1.76 1.92 2.20 2.20 1.67 1.92 1.90 1.92 1.95 1.92 1.77 1.76 1.76 1.76 2.20 2.21 2.23 * * *
102 N 26.11461° E 91.79327° 1.41 1.73 1.73 2.10 1.92 2.21 2.21 1.92 1.90 1.92 1.92 2.20 2.22 1.92 2.03 1.92 2.21 2.22 2.24 2.24
103 N 26.15819° E 91.81835° 1.68 1.68 1.67 1.67 1.67 1.62 1.62 1.44 1.44 1.44 1.71 1.71 2.20 2.21 2.22 2.22 2.24 2.26 2.26 2.30
104 N 26.13457° E 91.80951° 1.12 1.65 1.65 1.77 1.89 1.80 1.92 1.77 1.77 1.89 1.66 1.67 1.77 1.95 2.21 2.24 2.25 2.26 2.28 2.30
105 N 26.12234° E 91.80431° 1.65 1.70 1.77 1.76 2.21 2.22 1.95 1.89 1.89 1.04 1.52 1.92 1.78 1.78 1.93 1.92 2.19 2.21 2.22 2.23
106 N 26.15590° E 91.82724° 1.65 1.65 1.76 1.77 1.97 1.97 1.97 1.97 1.97 1.97 1.71 1.71 1.90 1.96 2.21 2.21 2.23 2.23 2.24 2.26
107 N 26.16614° E 91.82574° 1.61 1.61 1.78 1.69 1.69 1.77 1.92 1.93 1.76 1.76 1.71 1.71 1.97 2.19 2.19 2.26 2.26 2.10 2.24 2.26
108 N 26.13953° E 91.82994° 1.62 1.76 1.80 1.92 1.92 1.95 1.76 2.20 2.20 1.80 1.76 1.76 1.96 2.20 2.21 2.22 2.19 2.20 2.20 2.20
109 N 26.15677° E 91.83789° 1.92 2.02 2.02 1.92 1.92 1.88 2.19 1.44 1.44 1.75 1.76 1.47 1.47 1.47 1.92 1.95 1.95 1.95 2.20 2.22
110 N 26.17557° E 91.83054° 1.10 1.10 1.52 1.69 1.81 1.62 1.62 1.69 1.69 1.91 1.93 1.82 1.79 2.21 2.22 * * * * *
111 N 26.13730° E 91.82093° 1.44 1.76 1.79 1.80 1.75 1.89 1.89 1.89 2.20 2.20 1.74 1.74 1.92 1.92 1.92 1.92 1.92 1.92 1.92 1.93
112 N 26.13135° E 91.67260° 1.79 1.79 1.79 1.79 1.52 1.60 2.09 2.20 2.20 2.20 2.25 2.25 2.27 2.28 2.30 2.31 2.32 2.32 2.34 2.35
113 N 26.18258° E 91.80232° 1.60 1.79 1.73 1.73 1.73 1.78 1.80 2.20 2.20 1.86 1.85 1.87 1.90 2.20 2.20 2.20 * * * *
114 N 26.13245° E 91.81193° 1.65 1.92 1.92 1.68 1.85 1.85 1.85 1.77 1.78 1.58 1.87 1.87 1.88 1.88 1.88 1.89 2.20 2.20 2.20 2.25
115 N 26.12685° E 91.68861° 1.83 1.83 1.83 1.83 1.76 1.76 1.76 1.80 1.83 1.83 1.83 1.90 1.90 1.76 1.76 1.78 1.78 2.20 2.20 2.20
327
116 N 26.12910° E 91.67760° 1.71 1.17 1.74 1.80 1.80 1.93 2.20 2.20 2.20 2.20 2.20 2.25 2.25 2.25 2.25 2.25 2.25 1.91 1.93 *
328
Bore DENSITY (gm/c.c)
hole ID Latitude Longitude
117 N 26.14862° E 91.81223° 1.86 1.19 1.19 2.20 2.20 1.97 1.80 1.84 1.80 1.68 1.68 1.87 1.87 2.10 1.98 2.20 2.20 2.25 2.17 2.17
118 N 26.11473° E 91.74764° 1.51 1.45 1.45 2.02 2.02 2.01 2.01 2.01 2.01 1.76 1.80 1.80 1.80 2.20 2.20 2.20 2.25 2.25 2.20 2.30
119 N 26.11428° E 91.74064° 1.45 1.19 1.88 1.30 1.30 1.88 1.90 2.20 2.20 1.86 1.86 1.86 1.86 1.92 2.20 2.20 2.20 2.23 2.23 2.24
120 N 26.18015° E 91.81043° 1.95 1.60 1.75 1.85 1.80 1.95 1.85 1.85 1.85 1.85 1.85 1.82 1.82 1.82 1.82 1.82 1.43 2.20 2.20 2.25
121 N 26.10395° E 91.75164° 1.62 1.76 1.65 1.94 1.76 1.77 1.80 1.80 1.95 1.96 2.19 2.19 1.81 1.81 1.81 1.92 1.93 1.93 2.18 2.19
122 N 26.11338° E 91.72913° 1.21 1.76 1.68 1.68 2.21 2.22 2.22 2.24 2.25 2.26 2.28 2.30 2.30 2.21 2.23 2.25 2.25 2.25 2.22 2.22
123 N 26.17800° E 91.80292° 1.67 1.01 1.67 1.76 1.80 1.80 1.92 1.92 1.93 1.93 1.84 1.84 1.92 1.92 1.95 1.94 1.92 2.22 2.24 2.25
124 N 26.11935° E 91.81679° 1.60 1.60 1.95 1.80 1.80 1.92 1.93 1.92 1.93 1.92 1.75 1.75 1.75 1.75 2.20 1.44 1.80 1.80 2.21 2.22
126 N 26.18392° E 91.79392° 1.60 1.76 1.76 1.76 1.68 1.92 1.94 1.93 2.21 1.99 2.21 2.21 2.18 2.19 2.19 2.20 2.20 2.21 2.22 2.24
127 N 26.12281° E 91.82168° 1.69 1.75 1.80 1.84 1.80 1.84 1.94 1.95 1.92 1.84 1.84 1.84 1.84 1.84 1.94 1.92 2.00 * * *
128 N 26.12236° E 91.74914° 1.76 1.81 1.81 2.21 1.81 1.84 1.76 1.97 1.82 2.00 2.20 2.21 2.21 2.23 2.23 2.25 2.22 2.22 2.23 2.26
129 N 26.11922° E 91.68511° 1.83 1.83 1.61 1.80 1.99 1.99 2.21 2.22 2.22 2.22 2.23 2.23 2.23 2.24 2.24 2.25 2.25 2.26 2.28 *
130 N 26.10215° E 91.79817° 1.80 1.82 1.76 1.76 1.98 1.93 1.93 1.78 1.78 2.19 2.20 2.21 2.21 1.82 2.15 2.17 2.19 2.20 2.20 2.20
131 N 26.13000° E 91.82118° 1.95 1.95 1.80 1.76 1.80 1.97 2.19 2.21 1.97 1.93 1.66 2.20 2.22 2.25 2.27 2.30 2.30 * * *
132 N 26.13719° E 91.66910° 1.84 1.92 1.44 1.44 1.39 1.80 2.20 2.20 2.20 2.23 2.23 2.25 2.26 2.28 2.30 2.23 2.24 2.25 2.27 2.30
133 N 26.13719° E 91.79117° 1.31 1.68 1.68 1.68 1.60 1.60 2.17 2.17 2.17 2.17 2.20 2.20 2.25 2.25 2.25 2.25 * * * *
134 N 26.10799° E 91.71913° 1.44 1.44 1.89 1.93 2.07 2.20 2.21 2.23 2.25 2.25 2.25 2.25 2.27 2.27 2.28 2.30 2.31 2.31 2.33 2.35
135 N 26.12057° E 91.69611° 1.61 1.85 1.65 1.75 2.17 2.20 2.20 2.20 1.61 1.61 1.61 1.61 1.61 1.61 2.20 2.20 2.25 2.25 2.25 2.30
136 N 26.14931° E 91.67260° 1.35 1.75 1.61 1.61 1.80 1.80 1.92 1.99 2.01 1.97 2.20 2.23 2.25 2.25 2.27 2.30 1.47 1.52 1.52 1.52
137 N 26.11415° E91.71882° 1.15 1.85 1.85 1.85 1.85 1.78 2.20 2.21 1.80 1.99 1.97 * * * * * * * * *
138 N 26.11652° E 91.70812° 1.70 1.70 1.77 1.75 1.76 1.79 1.80 1.85 2.20 2.20 2.20 * * * * * * * * *
139 N 26.20272° E 91.68222° 1.13 1.74 1.74 1.68 1.76 1.80 1.90 1.80 1.82 1.92 1.82 1.95 1.92 1.80 1.92 1.82 2.20 * * *
140 N 26.12770° E 91.81221° 1.14 1.85 1.85 1.85 1.85 2.20 2.20 1.82 1.82 1.82 1.85 1.85 1.87 1.87 1.85 1.89 1.89 1.91 2.20 2.20
Annexure VII
Bore DENSITY (gm/c.c)
hole ID Latitude Longitude
141 N 26.11113° E 91.79567° 1.93 1.93 1.93 2.18 2.18 2.20 2.20 2.20 2.25 2.25 2.25 1.85 1.85 1.95 2.20 2.20 2.20 2.25 2.25 2.30
142 N 26.13567° E 91.79333° 1.83 1.83 1.83 1.84 1.84 1.84 1.84 1.84 1.89 1.89 2.20 2.20 2.20 2.25 2.25 2.30 2.30 2.30 2.35 2.35
143 N 26.14258° E 91.66810° 1.85 1.72 1.72 1.65 1.89 1.91 1.91 2.20 2.20 2.25 2.25 2.25 2.25 2.27 2.27 2.27 2.30 * * *
144 N 26.12057° E 91.71813° 1.44 1.99 1.99 1.99 1.99 1.85 1.81 1.93 2.20 2.20 2.21 2.21 2.23 2.23 2.25 2.30 2.30 2.30 * *
145 N 26.19374° E 91.71749° 1.87 1.87 1.87 1.87 1.87 1.87 1.78 1.78 1.78 1.78 2.20 2.20 2.20 2.25 2.25 2.25 2.30 2.30 2.30 2.35
146 N 26.20335° E 91.67313° 1.75 1.80 1.95 1.84 1.80 1.98 1.80 1.95 1.97 2.19 2.20 2.23 2.23 2.25 2.27 2.30 * * * *
147 N 26.12813° E 91.81169° 1.41 1.94 1.94 1.93 1.93 2.20 2.20 1.90 1.90 1.91 1.90 1.97 1.99 1.96 2.10 1.89 1.91 1.91 1.91 *
148 N 26.21956° E 91.73471° 2.25 2.25 2.25 2.26 2.26 1.49 1.68 1.68 1.80 1.91 1.91 2.20 2.22 2.24 2.25 2.25 2.27 2.28 2.30 *
149 N 26.22154° E 91.71588° 1.77 1.77 1.77 1.80 1.79 1.79 1.92 1.92 1.90 1.92 1.93 2.04 2.04 2.04 2.07 2.20 * * * *
150 N 26.13763° E 91.65009° 1.80 1.80 1.80 1.20 1.20 1.65 1.75 2.20 2.20 2.20 2.25 2.28 2.30 2.30 2.30 2.32 2.35 2.35 2.36 *
151 N 26.21292° E 91.72890° 1.79 1.79 1.79 1.79 1.80 1.80 1.80 1.81 2.21 2.22 2.24 2.24 2.26 2.26 2.27 2.28 2.30 2.30 * *
152 N 26.20357° E 91.71428° 1.51 1.84 1.85 1.83 1.84 1.82 1.89 1.91 2.20 2.21 2.22 2.23 2.25 2.27 2.27 2.28 2.30 2.31 2.31 *
153 N 26.22010° E 91.69906° 1.80 1.86 1.86 1.99 1.95 1.90 1.91 2.20 2.20 2.20 2.23 2.25 2.25 2.25 * * * * * *
154 N 26.19980° E 91.71488° 1.68 1.68 1.11 1.11 1.53 1.53 1.53 1.53 1.80 2.20 2.21 2.22 2.24 2.25 2.26 2.28 2.28 2.30 2.30 2.31
155 N 26.18705° E 91.70547° 1.79 1.79 1.92 1.84 1.83 1.84 1.86 1.80 1.86 2.20 2.21 2.22 2.23 2.25 2.27 2.28 2.28 2.30 2.32 2.35
156 N 26.21166° E 91.70928° 1.79 1.78 1.87 1.87 1.87 1.87 1.87 1.88 1.90 2.20 2.20 2.23 2.23 2.25 2.26 2.27 2.28 2.30 2.18 2.20
157 N 26.18680° E 91.67834° 1.40 1.40 1.49 1.49 1.49 1.49 1.89 1.89 1.61 1.61 1.61 1.95 1.82 1.84 1.96 1.96 1.99 2.01 2.01 1.99
158 N 26.19978° E 91.72659° 1.68 1.58 1.58 1.76 1.80 1.82 1.92 1.44 1.52 1.72 2.20 2.21 2.25 2.27 * * * * * *
159 N 26.19697° E 91.72006° 1.76 1.82 1.77 1.77 1.70 1.70 1.70 1.72 1.75 2.20 2.22 2.24 2.24 2.25 2.26 2.27 2.30 * * *
160 N 26.11761° E 91.82617° 1.65 1.68 1.99 1.99 1.99 1.99 1.97 1.99 1.97 1.99 1.92 1.97 1.99 1.97 1.95 1.95 1.97 1.92 1.97 1.99
161 N 26.22248° E 91.70677° 1.78 1.78 1.78 1.80 1.82 1.80 1.92 1.91 1.76 1.76 1.80 1.92 1.93 1.92 1.95 1.99 2.01 2.03 2.05 2.05
162 N 26.21036° E 91.73950° 1.68 1.78 1.78 1.44 1.52 1.60 1.95 1.30 1.68 1.76 1.78 1.80 1.60 1.76 2.21 2.22 2.22 2.23 2.24 2.25
163 N 26.21117° E 91.68965° 1.73 1.73 1.73 1.97 2.12 2.15 2.13 2.19 2.17 1.95 2.02 2.00 1.99 2.20 1.97 1.65 1.65 2.00 2.02 *
329
164 N 26.21898° E 91.72659° 1.73 1.73 1.92 1.95 1.58 1.52 1.52 1.58 1.68 1.68 1.75 1.76 2.20 2.21 2.22 2.24 2.27 * * *
330
Bore DENSITY (gm/c.c)
hole ID Latitude Longitude
165 N 26.20183° E 91.72322° 1.82 1.92 1.86 1.81 1.81 1.81 1.81 1.83 1.83 1.85 2.21 2.21 2.22 2.23 2.24 2.25 2.26 2.26 2.28 2.30
166 N 26.17940° E 91.74526° 1.08 1.11 1.81 1.81 1.81 1.52 1.67 2.19 2.20 2.21 2.20 2.21 2.21 2.22 2.23 2.24 2.25 2.27 * *
167 N 26.20672° E 91.70496° 1.39 1.44 1.34 1.34 1.58 1.52 1.58 1.87 2.21 2.21 2.22 2.22 2.23 2.24 2.25 2.26 2.26 2.27 2.28 2.30
168 N 26.21031° E 91.71638° 1.60 1.60 1.60 1.60 1.60 1.58 1.58 1.58 1.58 1.65 2.21 2.21 2.22 2.23 2.24 2.25 2.26 2.27 2.28 2.30
169 N 26.18579° E 91.69515° 1.12 1.12 1.43 1.43 1.72 1.72 1.80 1.86 1.90 * * * * * * * * * * *
170 N 26.18240° E 91.68619° 1.68 1.68 1.68 1.83 1.41 1.41 1.97 1.95 1.99 1.80 1.92 1.92 1.92 1.95 1.97 1.99 1.82 1.97 1.99 2.01
171 N 26.18240° E 91.69445° 1.68 1.66 2.09 2.09 2.09 1.60 1.60 1.52 1.52 1.46 1.44 1.54 1.99 1.95 1.97 1.99 1.70 1.76 1.70 1.84
172 N 26.15342° E 91.69371° 1.65 1.85 1.85 1.81 1.85 1.85 1.85 1.85 1.31 1.31 1.82 1.80 1.80 1.89 2.20 * * * * *
174 N 26.18406° E 91.66617° 1.10 1.19 1.50 1.50 1.44 1.44 1.46 1.44 1.64 1.52 1.50 1.50 1.68 2.10 1.87 1.76 1.76 1.97 1.98 2.03
175 N 26.15161° E 91.69443° 1.10 1.68 1.67 1.70 1.70 1.70 1.70 1.50 1.90 1.92 1.68 1.68 1.65 2.18 2.19 2.20 2.20 2.20 2.21 2.22
176 N 26.18217° E 91.67849° 1.20 1.53 1.53 1.53 1.53 1.55 1.55 1.55 1.56 1.76 1.86 1.92 1.80 1.65 1.55 1.55 1.68 * * *
177 N 26.18918° E 91.70176° 1.32 1.68 1.82 1.52 1.12 1.12 1.69 1.60 1.68 2.21 2.21 2.22 2.23 2.24 2.25 2.25 2.27 2.28 2.30 *
178 N 26.15486° E 91.70041° 1.68 1.70 1.80 1.80 1.80 1.48 1.78 1.99 1.93 1.93 1.94 1.55 1.59 1.60 1.77 * * * * *
179 N 26.18989° E 91.74735° 1.20 1.20 1.48 1.48 1.48 1.61 1.67 2.21 2.21 2.22 2.23 2.23 2.25 2.26 2.26 2.27 2.28 2.28 2.30 2.31
180 N 26.19969° E 91.69440° 1.54 1.54 1.54 1.92 1.92 1.54 1.92 1.92 1.97 1.99 2.04 1.95 1.92 1.76 1.92 2.21 * * * *
181 N 26.16086° E 91.74903° 1.85 2.05 2.05 2.05 2.05 1.95 1.80 1.95 2.20 2.21 2.22 * * * * * * * * *
182 N 26.19750° E 91.75987° 1.56 1.56 1.56 1.56 1.44 1.44 1.50 1.32 1.52 1.60 1.60 1.58 1.58 1.60 1.76 1.76 1.82 1.85 1.98 1.87
183 N 26.18907° E 91.75780° 1.72 1.72 1.92 1.80 1.76 2.05 2.01 2.00 2.01 1.92 2.05 2.04 1.92 2.01 2.21 2.21 * * * *
184 N 26.16324° E 91.76584° 1.70 1.85 1.85 1.97 1.97 2.12 2.17 1.82 1.82 1.82 1.82 1.60 2.19 2.20 2.20 2.21 2.22 2.24 2.12 *
185 N 26.18040° E 91.75342° 1.65 1.82 1.80 1.80 1.80 1.80 1.80 1.97 2.12 2.19 2.19 2.20 2.20 2.21 2.22 2.24 2.24 2.26 * *
186 N 26.18488° E 91.75740 ° 1.37 1.37 1.68 1.55 1.68 1.68 1.80 1.92 1.92 1.92 1.85 1.88 1.90 1.90 2.20 2.21 2.21 2.23 2.24 *
187 N 26.16015° E 91.74101° 1.10 1.12 1.92 1.92 1.92 1.80 1.97 2.12 2.17 2.20 2.21 2.22 2.24 * * * * * * *
188 N 26.16733° E 91.75239° 1.60 1.77 1.73 1.73 1.73 1.73 1.73 2.19 2.20 1.92 1.95 2.20 2.20 1.99 2.19 2.19 2.20 1.67 1.77 1.70
Annexure VII
Bore DENSITY (gm/c.c)
hole ID Latitude Longitude
189 N 26.17849° E 91.75786° 1.10 1.30 1.87 1.87 1.87 1.87 1.93 1.93 1.93 1.93 2.19 2.20 2.20 2.21 2.21 2.23 * * * *
190 N 26.14839° E 91.74357° 1.75 1.77 1.77 1.77 1.85 1.99 1.99 2.02 2.02 2.03 * * * * * * * * * *
191 N 26.17135° E 91.74653° 1.10 1.85 1.80 1.80 1.92 1.68 2.19 2.20 2.21 2.21 2.21 2.22 2.22 2.24 2.25 * * * * *
192 N 26.16959° E 91.77882° 1.10 1.10 1.12 1.77 1.77 1.85 1.90 1.60 1.60 1.60 1.73 1.73 1.95 1.95 2.02 2.20 2.20 1.98 1.98 1.98
193 N 26.14772° E 91.74207° 1.68 1.85 2.11 2.11 2.11 1.92 1.80 1.85 1.85 1.50 1.44 1.49 1.85 1.60 2.19 2.20 2.20 2.20 * *
194 N 26.15284° E 91.73533° 1.10 1.68 1.83 1.76 1.92 2.18 2.19 2.19 2.20 2.20 * * * * * * * * * *
195 N 26.18082° E 91.74796° 1.68 1.81 1.81 1.68 1.68 1.80 1.80 2.09 2.19 2.19 2.20 2.20 2.21 2.21 2.23 2.23 2.24 2.25 2.25 2.26
196 N 26.13603° E 91.73539° 1.60 1.45 1.45 2.18 1.92 2.19 2.20 2.21 2.21 2.22 2.23 2.25 2.25 2.26 2.27 2.27 2.28 2.28 2.30 2.30
197 N 26.13700° E 91.74526° 1.01 1.65 1.75 1.75 1.60 1.80 1.97 1.85 1.87 1.95 2.02 2.19 2.20 2.20 2.20 2.24 * * * *
198 N 26.12811° E 91.73819° 1.12 1.92 1.92 1.97 1.97 1.92 1.99 1.95 1.89 1.89 2.20 2.20 2.24 2.25 2.26 2.27 2.30 2.17 2.15 2.17
199 N 26.17859° E 91.82757° 1.08 1.40 1.40 1.40 1.65 1.65 1.80 1.68 2.20 2.20 2.20 2.21 * * * * * * * *
200 N 26.18666° E 91.84341° 1.80 1.40 1.40 2.02 2.05 1.99 1.98 2.00 2.01 2.17 1.94 1.60 1.60 2.19 2.19 2.20 2.20 2.21 2.22 2.25
331
ANNEXURE VIII