RUSSIAN JOURNAL OF EARTH SCIENCES VOL. 10, ES5002, doi:10.2205/2008ES000303, 2008
Analysis of seismicity in North IndiaV. Gitis, and E. Yurkov Institute for Information Transmission Problems, RAS, Moscow, Russia B. Arora, S. Chabak, and N. Kumar Wadia Institute of Himalayan Geology, Dehradun, India P. Baidya India Meteorological Department, New Delhi, IndiaContents
Abstract[1] The geotectonic information on North India region has been processed using GIS (Geographic Information Systems), GeoProcessor for extracting the spatial seismotectonic properties. According to spatial modeling of the principle seismotectonic parameters the region under study can be divided into four zones as Kangra-Chamba, Garhwal-Kumaon, Punjab reentrant and Delhi-Haridwar. The b-value is low (0.5) for the Kangra-Chamba and Garhwal-Kumaon region, which shows that the regions have less capacity to withstand the developing stress. The elliptical region between Main Central thrust (MCT) and Main Boundary Thrust (MBT) are sensitive to the tidal force. 1. Introduction[2] Since last 20 years, the Wadia Institute of Himalayan Geology (WIHG), Dehradun, India has been supporting the multidisciplinary research in exploring the structure and history of the Himalaya based on seismotectonic investigations. The institute has collected a complementary dataset of geophysical time series at different locations in Northern India (26o-34o N and 74o-82o E) containing the information of seismology, tectonics deformation. The Himalaya region has many thrust faults capable of producing the earthquakes of magnitude 8.0 or greater. Some of these faults are very prominent and visible at the surface but some faults are hidden. These faults are mostly concentrated in the mountain region producing hundreds of earthquakes every year though most are too small to detect. Recently an earthquake of M=7.6 occurred at Muzzfrabad on 8 October 2005 at northwest boundary of present study region [Rao et al., 2006].
[3] Himalaya, the arcuate mountain belt of complex geotectonic setup stretching about 2400 km
long in east-west direction with variable width of 230 to 320 km is formed due to the convergent
movement of two plates of the earth's lithosphere. The Indian and Asian continental plates
collided some 50 m.y. ago [Le Fort, 1975] resulting lithosphere deformation and modification of
the seismotectonic model of the region with the span of time. The seismotectonic investigations
have been done by many authors [Seeber and Armbruster, 1981; Ni and Barazangi, 1984;
Thakur and Kumar, 2002; Kayal, 2007;
Bollinger et al., 2007; Kumar et al., 2008, (in press)] and for the region of
NW Himalaya; we have well documented information of great ( M [4] This work addresses the results on development of geoinformation technology oriented to seismotectonic problem domain [Gitis and Ermakov, 2004] and application of the technology to spatial analysis of seismicity in North India region.
[5] Web-GIS GeoProcessor 1.5 [Gitis et al., 1998; Gitis, 2004]
and GIS SeismoTide [Yurkov and Gitis, 2005] were used. Further we will consider the
initial data and geoinformation tools, after
then we present three main results of data exploration: spatial patterns of main seismic
parameters, some relationships between seismicity and components of tidal force, and spatial
analysis of relationship between strong earthquakes with
M 2. Geotectonic Model of North India Region
[7] The HFT or MFT (Main Frontal Thrust) is the southern most, younger and neotectonically active thrust of the Himalayan region giving a topographical break of the Himalaya (the Outer Himalaya), against the alluvial Indo-Gangetic Plains. The Outer Himalaya comprises the low altitude Siwalik Hills with flat-floored structural valleys consisting of about 9500 m thick Cenozoic sedimentary pile characterized by folds and faults [Thakur, 1992]. The Indo-Gangetic plains of thick alluvial lie to the south of this discontinuity that is formed due to thick deposited sediments by the Indus, the Ganga and other Rivers. The Delhi-Hardwar ridge is the important structural feature in the Ganga basin which is tectonically and seismically active as compared to the rest of the basin. [8] The evolutionary model of Himalaya [Le Fort, 1975] states that MBT is the younger tectonic discontinuity as compare to MCT, which is more active currently. However, both the MCT and MBT has been treated as the contemporaneous features in the steady-state model of Seeber and Armbruster [1981] and these merge with each other at depth with a common detachment surface. The MBT is the tectonic boundary between the Lesser Himalaya nappes lying to its north and the Tertiaries of the Outer Himalaya foreland basin to the south. The majority of the earthquakes in the NW Himalaya are concentrated in the zone between MBT and MCT of shallow focal depth and great Himalayan earthquakes are originated at the surface of the detachment, which represents the upper surface of the underthrusting Indian plate with apparent northward dip of about 15o [Ni and Barazangi, 1984]. [9] The MCT is the major tectonic discontinuity, which divide the two contrasting structures as Lesser and Higher Himalaya having contrasting stratigraphic and tectonic features. The relative movement of the blocks across this tectonic discontinuity has caused the development of crustal buckles in which the Palaeogene and Neogene sediments have deposited. At depth the MCT represents a ductile shear zone and tectonically it is a duplex shear zone having three distinct thrust planes as MCT I, MCT II and MCT III. According to the degree of metamorphism, lithostratigraphy and tectonic setting, Valdiya (1980) has named these thrusts as Chail (MCT I, lower thrust), Jutogh (MCT II, middle thrust) and Vaikrita (MCT III, upper thrust). [10] The ITSZ is the discontinuity along which the initial contact between the two plates occurred in the Early Eocene (60-70 my ago). The Tethys Sea, which was existed between two continents, vanished to the end of Eocene and after that the ITSZ converted to active plate boundary. The region between ITSZ and the MCT is composed of the Tethyan sedimentary zone in the north and the central crystalline zone in the south. The suture is made up of imbricated mélanges of flysch sediments, pillow lavas, volcanics and the basic and ultrabasic rocks, which are cut by steep faults [Allegre et al., 1984]. According to Burg et al., [1984], the Tethyan slab consists of several thrust sheets with distinct ages and deformations. Mainly normal faulting is existed in the high Himalaya and the normal fault zone separates the highly metamorphosed crystalline rock and a transitional zone from the overlying sedimentary sequence. 3. Initial Data Set and Geoinformation Modelling Tools
[11] The following initial data were used for geoinformation analysis: historical and instrumental
earthquake catalogue of Wadia Institute of Himalayan Geology from 1552 till 30.01.2005, digital
elevation model in grid
30'' [12] Geoinformation modeling tools were presented by Web-GIS GeoProcessor 1.5 and GIS SeismoTide. [13] Web-GIS GeoProcessor 1.5 (http://www.iitp.ru/projects/geo) is implemented as Java 1.5 applet. The applet can be used for decision support in environmental zonation, natural hazard and risk assessment, and exploration of natural resources. The system provides the following facilities of complex geodata analysis:
[14] GIS SeismoTide operates in the Matlab environment. It was designed specially for calculating tidal force characteristics, assessing the measure of correlation between seismicity and tidal force, and composing the maps of prevailing seismic activity. Each element of such a map shows phases (positive or negative) of tidal force components corresponding to higher recurrence rates of earthquakes. These maps facilitate the identification of seismically sensitive areas to the tidal force. Detection of the area is based both on visual analysis of the prevailing seismicity spatial pattern and on a statistical measure relationship between the seismic activity predominance and the different tidal force phases. 4. Spatial Modelling of Seismological Parameters[15] Earthquake catalogue of Wadia Institute of Himalayan Geology includes 2628 seismic events occurred from 1552 till 30.01.2005. Seismological network was significantly improved in 1988. Therefore 1856 earthquakes from 01.01.1999 till 30.01.2005 were used for spatial modeling. [16] These data were used for compiling the following three grid based models:
[17] To compile the models the seismic data are scanned with moving spatio-temporal window.
Seismic parameters are estimated by the earthquakes with the epicenters inside the window.
Evidently, the accuracy of the estimates depends on a number of the epicenters. It is supposed
that a set of epicenters inside the window is statistically uniform. For that the window must cover
homogeneous area in seismotectonic sense. Increasing the window size leads to doubtful
estimates because of mixing evidences from different probabilistic distributions. Decreasing the
window leads to uncertainties because of a small number of a sample set. A compromise consists
in adaptive fitting of the window size: the more density of seismic events the less the window
size.
[18] Moving window in our case is a cylinder with fixed height equal temporal interval from
01.01.1999 till 30.01.2005 and with radius which adaptively varies from
R1 to
R2. Minimal
number of seismic events in window
N min and nominal number of events
N are specified.
[19] The adaptive window algorithm is in following. Axis of the cylinder coincide with a
grid point
(l, j, t). If a number of events in cylinder of radius
R1 is
n [20] Minimal representative magnitude of earthquakes m0 is the left boundary of the interval, in which Gutenberg-Richter low is hold. Value m0 depends on number of seismic stations, their spatial distribution, and sensitivity. The method of estimation of m0 is based on testing of statistical hypothesis. The method was developed by Pisarenko, [1989], algorithm was developed by Smirnov, [1995], and together with Smirnov was modified for GIS GeoTime [Gitis et al., 1994].
![]() [22] Maximum likelihood estimate of b-value is following:
![]()
where
b(l, j, t) = (1Nm0
Nm0 = Nm0(l, j, t) is a number of earthquakes with
m
![]()
[24] Seismic activity
A(l, j, t) is defined as normalized according to spatio-temporal window, the
number of earthquakes adjusted to the magnitude
m
where
Nm0(l, j, t) is a number of earthquakes with the magnitudes
m
![]() 5. Relationship Between Seismicity and Components of Tidal Force[26] Relationship between tidal force and frequency of earthquakes was investigated. Tidal force was calculated from a model of the gravitational interaction between the Earth and Moon. We neglect the solar component of the tidal force because, on average, it is more than two times less of the lunar component. The disregard of this component does not hinder the establishment of the fact of existence of statistical relations between seismicity and earth tides.
[28] Two phases: positive ("+'') and negative ("-'') were considered for each characteristic of tidal force. The phase was considered positive if the characteristic exceeds its average value for the long-term period. Otherwise the phase was considered negative. The phase changes a sign one or two times per day for the first five characteristics and two times per month for the daily variation A r. Duration of each phase depends on a sort of the characteristic and from latitude of a place of supervision and does not depends on a longitude. [29] For revealing correlations between seismicity and tidal force the special statistical measure was used
where f+ = N+/T+ and f- = N-/T- are the frequencies of earthquakes for "+'' and "-'' phases of tidal force correspondingly, N+ and N- are the numbers of earthquakes for "+'' and "-'' phases, T+ and T- are the total duration of intervals for "+'' and "-'' phase. The sign of s variable specifies in what phase ("+'' or "-'') earthquakes occur more often. Measure of s possesses statistical properties: if N+ and N- are sufficiently great, then value s has standard normal distribution and so its big values (for example, 2 or 3) are "significant'', i.e. argue in favor of correlation between seismicity and tidal force.
[30] Earthquakes with the magnitudes
m [31] Under close examination of the results some "positive'' or "negative'' cells form clusters of the cells. Such clusters are approximated by ellipses. Frequencies f+, f- and measure s are calculated for the seismic events which fit in the ellipses.
[33] The interpretation of this result is based on the assumption that a unidirectional action of the tectonic and tidal pressures increases probability the earthquake. We shall notice that vertical tidal force in a negative phase causes compression of the earth crust. Setting the intensity of earthquakes during a negative phase of characteristic A r is raised, it is necessary to conclude, that area AI is in a condition of primary vertical tectonic compression. 6. Spatial Analysis of Relationship Between Strong Earthquakes With M≥6 and Digital Elevation Model
[35] This is the zone where Himalayan mid-crustal ramp under the southern Higher Himalaya has been proposed on the basis of various geophysical studies [Gahalaut and Kalpna, 2001; Yeats and Thakur, 1998]. Throughout Himalaya and mainly in NW Himalaya, a narrow belt between Higher and Lesser Himalaya shows intense seismic activity, the region is close to ramp and better known as Himalayan Seismic Belt (HSB). Yeats and Thakur, [1998] has given a fault-bend fold model where the megathrust drives southward and upward over the ramp and the axial surface of the fold is active at crustal sacle only. The central Himalaya moves southward as a fault bend fold due to slip by earthquakes that nucleate the detachment under Himalaya of this region. That study also shows maximum relief over the structural ramp indicating a high uplift rates over the ramp. Therefore a high uplift rates in the central Himalaya above the frontal ramp is due to thrust movement during megathrust earthquakes identified by recent compiled seismic catalogue. 7. Discussion
[36] Presently, most of the analysis is done for the seismic data of January 1999 to January 2005
with the improved quality of the catalogue for this period. The grid based modeling of spatial
analysis of seismicity (Figure 2) shows that the region around Delhi has lowest value (2.0) of the
minimal representative magnitudes, while its value is highest (2.6) in the region of Punjab
reentrant and SW of Garhwal-Kumaon region. The seismic activity of earthquakes with
m=3 is
highest to the SW of Garhwal-Kumaon region and lowest in the Delhi-Haridwar region (Figure 3),
While the b-value pattern is opposite at these two regions (Figure 4). The b-value is around 0.5 for
Kangra-Chamba region and less than 0.5 for the Garhwal-Kumaon region. The low b-value of
these two regions can be stated that the regions have less capability to withstand the developing
stress. The Lunar component of tidal force is compared with the seismic activity of the region in
Figure 6 that shows the area AI seismically active to the daily variation of tidal force
A r. The level
of statistics
s in this area is high enough to testify essential distinction between frequencies of
earthquakes for two compared phases. This result manifested that area AI is in a condition of
primary vertical tectonic compression. The standard deviation of the Earth surface elevations
(consequences of geotectonic deformation) with RMS value more than 500 m is positively
correlated with the strong earthquakes
M Acknowledgments[37] This work was carried out under ILTP cooperation program. It was supported in part by ILTP, by the Russian Foundation for Basic Research, project no. 06-07-89139 and by Basic research program of Presidium of RAS No 15, section "Electronic Earth". We are grateful to colleagues in the projects for helpful discussions. ReferencesAllegre, C., J., and, et al. (1984), Structure and evolution of the Himalaya-Tibet orogenic belt, Nature, 307, 17, doi:10.1038/307017a0. [CrossRef] Bollinger, L., F. Perrier, J.-P. Avouac, S. Sapkota, U. Gautam, and D. R. Tiwari (2007), Seasonal modulation of seismicity in the Himalaya of Nepal, Geophys. Res. Lett., 34, L08304. Burg, J. P., M. Guiraurd, G. M. Chenn, and G. C. Li (1984), Himalayan metamorphism and deformations in the north Himalaya Belt (Southern Tibet, China), Earth Planet. Sci. Lett., 69, 391, doi:10.1016/0012-821X(84)90197-3. [CrossRef] Gahalaut, V. K., and Kalpna (2001), Himalayan mid-crustal ramp, Curr. Sci., 81, (12), 1641. Gitis, V. (2004), Experience of spatio-temporal seismotectonic data mining in multidisciplinary measurements, in: Proceedings of European Seismological Commission XXIX General Assembly, p. 146, ESC, Potsdam. Gitis, V., A. Dovgyallo, B. Osher, and T. Gergely (1998), GeoNet: an information technology for WWW on-line intelligent Geodata analysis, in: Proceedings of the 4th EC-GIS Workshop, p. 124, Joint Research Centre of European Commission, Hungary. Gitis, V., and B. Ermakov (2004), Fundamentals of spatio-temporal forecasting in geoinformatics, 256 pp., Fizmatlit, Moscow. Gitis, V. G., B. V. Osher, S. A. Pirogov, A. V. Ponomarev, G. A. Sobolev, and E. F. Jurkov (1994), A System for Analysis of Geological Catastrophe Precursors, Journal of Earthquake Prediction Research, 3, 540. Kayal, J. R. (2007), Recent large earthquakes in India: Seismotectonic Perspective, IAGR Memoir, 10, 189. Le Fort, P. (1975), Himalaya: the collided range, present knowledge of the continent arc, Am. J. Sci., 275, 7. Ni, J., and M. Barazangi (1984), Seismotectonics of the Himalayan collision zone: Geometry of the Underthrsuting Indian Plate beneath the Himalaya, J. Geophys. Res., 89, 1147. Pisarenko, V. (1989), About frequency-magnitude relationship of earthquakes, in: Discret properties of geological environment, p. 47, Nauka, Moscow. Rao, N. P., P. Kumar, Kalpna, T. Tsukuda, and D. Ramesh (2006), The devastating Muzaffarabad earthquake of 8 October 2005: New insights into Himalayan seismicity and tectonics, Gondwana Research, 9, (4), 365, doi:10.1016/j.gr.2006.01.004. [CrossRef] Seeber, L., and J. G. Arnbruster (1981), Great detachment earthquakes along the Himalayan arc and long forecasting, in: Earthquake Prediction - an international review, 4, p. 259, Am. Geophys. Union., Maurice Ewing Ser., Washington. Smirnow, V. (1995), Earthquake recurrence and seismic regime parameters, Volcanology and seismology, (3), 59. Thakur, V. C. (1992), Geological Map of Western Himalaya, 22 pp., Wadia Inst. Him. Geol., Dehradun. Thakur, V. C., and S. Kumar (2002), Seismotectonics of Chamoli Earthquake of March 29, 1999 and Earthquake Hazard Assessment of Garhwal-Kumaon Region, NW Himalaya, Himalayan Geology, 23, (1-2), 113. Yeats, R. S., and V. C. Thakur (1998), Reassessment of earthquake hazard based on a fault-bend fold model of the Himalayan plate-boundary fault, Curr. Sci., 74, (3), 220. Yurkov, E., and V. Gitis (2005), Relation between Seismicity and Phases of Tidal Waves, Izvestiya, Physics of the Solid Earth, 41, (4), 255. Received 4 June 2008; accepted 10 June 2008; published 18 June 2008. Keywords: seismicity, tectonics, geographic information systems, Himalaya. Index Terms: 0500 Computational Geophysics; 0530 Computational Geophysics: Data presentation and visualization; 7230 Seismology: Seismicity and tectonics; 8010 Structural Geology: Fractures and faults. ![]() Citation: 2008), Analysis of seismicity in North India, Russ. J. Earth Sci., 10, ES5002, doi:10.2205/2008ES000303. (Copyright 2008 by the Russian Journal of Earth SciencesPowered by TeXWeb (Win32, v.2.0). |