Aftershocks of 26th january 2001 Bhuj earthquake and seismotectonics of the Kutch region

The 26th January 2001 Bhuj earthquake was followed by intense aftershock activity. Aftershock data from United States Geological Survey (USGS) utilized in this study encompasses three months period from 26th January to 26th April 2001. Epicenters of the aftershock are plotted on a map depicting active faults. All the aftershocks of magnitude > 5 and 70% of those ranging between magnitude 3 and 5 are confined to an area resembling a horseshoe pattern with a pointed end towards NE. The other 20% of magnitude 3 to 5 are enclosed within an almost parallel boundary. Only 10% are found to be beyond this limiting boundary. 50% of the recorded after-shocks took place within the first week of the main event and this study reveals that the basic characteristic pattern of aftershock activity can be determined on the basis of the data of only one week. Four major NW-SE trending active faults are mapped in the Kutch region. They define the western limit of Cambay structure and also mark the western limit of Dharangadhra and Wadhwan basins along the SE continuation in Saurashtra. These faults separate the Kutch region into two geologically different blocks. On the SW side the mapped horseshoe pattern gets characteristically truncated along the western most fault, which is characterized by a strike-slip movement in the south and vertical movement in the north. The present study has revealed that the epicenter of the 26th January earthquake is located in the vicinity of the Bhachau township, close to the intersection with the Kutch mainland fault. Furthermore, it has been noticed that most of the epicenters of the aftershock are confined in the intersectional area of the Kutch mainland fault and the NW-SE faults.     

Statistical parameters of Bhuj earthquake sequence of January 26th, 2001

An intraplate earthquake of magnitude (M _c) 6.9 (Anon 2001a) struck Bhuj and the adjoining region of Kachchh in Gujarat on January 26th, 2001 at about 0316 hrs (GMT) and was followed by a number of aftershocks. The epicentre of this earthquake was located at 23.4‡N and 70.28‡E close to the Kachchh mainland fault. The intensity observed around the epicenter was X on the MSK scale. A study of 531 aftershocks, in the magnitude range of 3.0–5.7, recorded at Vadodara Seismological Observatory till March 31st, 2001 has been carried out and various statistical parameters calculated. The total energy released during the study period is calculated to be 8.2 × 10¹⁴ joule. Sudden occurrence of the main shock without any foreshock in the same tectonic system is a unique feature of this sequence. Theb- value (0.86), value of M₀-M₁ (1.2), high M₁/M₀ (0.89) and high value of the decay constanth (0.91), all support the tectonic origin of the present study.     

Fault plane solutions of the january 26th, 2001 bhuj earthquake sequence

A 12-station temporary microearthquake network was established by the Geological Survey of India for aftershock monitoring of the January 26th, 2001 Bhuj earthquake (M _w 7.6) in the Kutch district of Gujarat state, western India. The epicentres of the aftershocks show two major trends: one in the NE direction and the other in the NW direction. Fault-plane solutions of the best-located and selected cluster of events that occurred along the NE trend, at a depth of 15–38 km, show reverse faulting with a large left-lateral strike-slip motion, which are comparable with the main-shock solution. The NW trending upper crustal aftershocks at depth <10 km, on the other hand, show reverse faulting with right-lateral strike-slip motion, and the mid crustal and lower crustal aftershocks, at a depth of 15–38 km, show pure reverse faulting as well as reverse faulting with right-lateral and left-lateral strike-slip motions; these solutions are not comparable with the main-shock solution. It is inferred that the intersection of two faults has been the source area for stress concentration to generate the main shock and the aftershocks.     

Study of the epicentral trends and depth sections for aftershocks of the 26th january 2001, Bhuj earthquake in Western India

The Geological Survey of India (GSI) established a twelve-station temporary microearthquake (MEQ) network to monitor the aftershocks in the epicenter area of the Bhuj earthquake (M _w7.5) of 26th January 2001. The main shock occurred in the Kutch rift basin with the epicenter to the north of Bhachao village, at an estimated depth of 25 km (IMD). About 3000 aftershocks (M _d ≥ 1.0), were recorded by the GSI network over a monitoring period of about two and half months from 29th January 2001 to 15th April 2001. About 800 aftershocks (M _d ≥ 2.0) are located in this study. The epicenters are clustered in an area 60 km × 30 km, between 23.3‡N and 23.6‡N and 70‡E and 70.6‡E. The main shock epicenter is also located within this zone. Two major aftershock trends are observed; one in the NE direction and other in the NW direction. Out of these two trends, the NE trend was more pronounced with depth. The major NE-SW trend is parallel to the Anjar-Rapar lineament. The other trend along NW-SE is parallel to the Bhachao lineament. The aftershocks at a shallower depth (<10km) are aligned only along the NW-SE direction. The depth slice at 10 km to 20 km shows both the NE-SW trend and the NW-SE trend. At greater depth (20 km–38 km) the NE-SW trend becomes more predominant. This observation suggests that the major rupture of the main shock took place at a depth level more than 20 km; it propagated along the NE-SW direction, and a conjugate rupture followed the NW-SE direction. A N-S depth section of the aftershocks shows that some aftershocks are clustered at shallower depth ≤ 10 km, but intense activity is observed at 15–38 km depth. There is almost an aseismic layer at 10–15 km depth. The activity is sparse below 38 km. The estimated depth of the main shock at 25 km is consistent with the cluster of maximum number of the aftershocks at 20–38 km. A NW-SE depth section of the aftershocks, perpendicular to the major NE-SW trend, indicates a SE dipping plane and a NE-SW depth section across the NW-SE trend shows a SW dipping plane. The epicentral map of the stronger aftershocksM ≥ 4.0 shows a prominent NE trend. Stronger aftershocks have followed the major rupture trend of the main shock. The depth section of these stronger aftershocks reveals that it occurred in the depth range of 20 to 38 km, and corroborates with a south dipping seismogenic plane.     

Behaviour of masonry structures during the Bhuj earthquake of january 2001

A variety of masonry structures suffered damage during the recent Bhuj earthquake. Some of the traditional masonry structures had no earthquake resistant features and suffered considerable damage. This paper attempts to evaluate the behaviour of masonry structures based on the type of masonry used in places like Bhuj, Anjar, Bhachau, Morbi, Samakhyali and several other places. Quite a few masonry buildings had used earthquake resistant features like lintel bands and corner reinforcements. The cracking and failure patterns of such buildings have also been examined. The paper concludes with a discussion on the relevance of the current codal provisions for earthquake resistance of masonry structures and the direction of further research in the area.     

Geoseismological investigation of the January 26, 2001 Bhuj earthquake in western peninsular India

Of the intraplate seismic events, the January 26, 2001 Bhuj earthquake (Mw 7.7) would be remembered as one of the deadliest, in which 13,805 human lives were lost, 0.177 million injured and a total of 1,205,198 houses were fully or partly damaged in 16 districts of Gujarat state with an estimated overall loss of Rs. 284, 23 million. The brunt of the calamity was borne by five districts, namely Kachchh, Ahmadabad, Rajkot, Jamnagar and Surendranagar, where 99?% of the total casualties and damage occurred. In the neighbouring parts of Sindhh Province of Pakistan, 40 human casualties were reported, and some buildings cracked in the Karachi city as well. In the Kachchh district of Gujarat state, the telecommunication links and power supply were totally disrupted, road and rail links partially impaired and water supply snapped at many places. The Bhuj airbase had to be closed for some time due to damage to the infrastructure. The macroseismic survey carried out by the Geological Survey of India in an area as large as 1.2 million?sq?km indicated an epicentral intensity as high as X on the MSK scale in an area of 780?sq?km in the central part of Kachchh rift basin. Apart from damages to civil structures, the January 26 earthquake induced conspicuous terrain deformation in the form of liquefaction features, structural ground deformation and low-order slope failures that were mainly prevalent within the higher intensity isoseists. Liquefaction occurred in an area of about 50,000?sq?km. The extensive plains of Rann of Kachchh, the marshy tracts of the Little Rann and the shallow groundwater table zones of Banni Land provided the most conducive geotechnical environments for the development of seismites. The liquefaction activity was profuse in seismic intensity zones X and IX, widespread in intensity VIII, subdued in intensity VII and stray in intensity VI. The common forms of liquefaction were sand blows/boils, ground fissures, craters, lateral spreading and slumping. Ground deformation of tectonic origin was witnessed in the epicentral tract. Such features, though much less subdued in comparison with the 1819 large earthquake (Mw 7.8) in region, occurred along the Kachchh Mainland fault (KMF) and along a transverse lineament, referred to as Manfara?CKharoi fault. The manifestations were in the form of fractures, displacement of strata, linear subsidence, upheaval, formation of micro-basins/micro-ridges, ripping off of rock surface, and at places violent forms of liquefaction. The localities where coseismic deformations were observed include Bodhormora, Sikra, Vondh, Chobari, Manfara and Kharoi. The 2001 event has demonstrated the role of local geology in influencing the ground motion characteristics and, therefore, the hazard estimation.     

A media-based assessment of damage and ground motions from the january 26th, 2001<Emphasis Type="Italic">M</Emphasis> 7.6 Bhuj,India earthquake

We compiled available news and internet accounts of damage and other effects from the 26th January, 2001, Bhuj earthquake, and interpreted them to obtain modified Mercalli intensities at over 200 locations throughout the Indian subcontinent. These values are used to map the intensity distribution using a simple mathematical interpolation method. The maps reveal several interesting features. Within the Kachchh region, the most heavily damaged villages are concentrated towards the western edge of the inferred fault, consistent with western directivity. Significant sedimentinduced amplification is also suggested at a number of locations around the Gulf of Kachchh to the south of the epicenter. Away from the Kachchh region intensities were clearly amplified significantly in areas that are along rivers, within deltas, or on coastal alluvium such as mud flats and salt pans. In addition we use fault rupture parameters inferred from teleseismic data to predict shaking intensity at distances of 0–1000 km. We then convert the predicted hard rock ground motion parameters to MMI using a relationship (derived from internet-based intensity surveys) that assigns MMI based on the average effects in a region. The predicted MMIs are typically lower by 1–2 units than those estimated from news accounts. This discrepancy is generally consistent with the expected effect of sediment response, but it could also reflect other factors such as a tendency for media accounts to focus on the most dramatic damage, rather than the average effects. Our modeling results also suggest, however, that the Bhuj earthquake generated more high-frequency shaking than is expected for earthquakes of similar magnitude in California, and may therefore have been especially damaging.     

Postseismic relaxation due to Bhuj earthquake on January 26, 2001: possible mechanisms and processes

Earthquakes cause static stress perturbations in the nearby crust and mantle. Obeying rheological laws, this stress relaxes in a time frame of months to years with the spatial extent of few km to hundreds of km. While postseismic relaxation associated with major inter-plate earthquakes is well established, there have been few opportunities to explore its occurrence following intraplate earthquakes. The M _w 7.6 Bhuj earthquake on January 26, 2001 in western India is considered to be an intraplate event and provided a unique opportunity to examine post-earthquake relaxation processes sufficiently away from plate boundaries. To study the characteristics of transient postseismic deformation, six Global Positioning System campaigns were made at 14 sites. The postseismic transients were delineated after removing plate motions from the position time series. Postseismic deformation has been observed at all the sites in the study area. During 2001?C2007, the site closest to the epicenter exhibited postseismic deformation of about 30 and 25?mm in the north and east components, respectively. Time series of the NS and EW components of the postseismic transients can be fitted to both logarithmic and exponential functions. Close to the epicenter, the logarithmic function fits well to the initial transient, and an exponential function fits well to the later phases. The remaining sites (located east and west of the epicentral region) exhibited significantly diminished north?Csouth relaxation. Rapidly decaying afterslip and poroelastic mechanisms seem to be responsible for postseismic relaxation in the vicinity of epicenter during the initial period subsequent to the Bhuj earthquake. Postseismic relaxation by viscoelastic flow below the seismogenic zone seems to affect displacements across the entire Bhuj region. This paper presents the characteristics of postseismic transients and deformation processes in the scenario of the highly heterogeneous crust in the Bhuj region.     

Post-seismic deformation associated with the 2001 Bhuj earthquake

We report the results of GPS measurements of post-seismic deformation due to the 2001 Bhuj earthquake in the Kachchh region, western India. The estimated horizontal velocity vectors in ITRF05 are in the range of 48?C49?mm/year in N46?C50°E. The observed velocity at the Gandhinagar permanent site, a far off site from the earthquake source region and probably unaffected by the post-seismic deformation, is 49?±?1?mm/year in N47°E, which is consistent with the predicted motion of Indian plate at Gandhinagar. At other sites in the source region, transient post-seismic deformation is found to be low; it attenuated rapidly within 3?C4?years of the earthquake and is much low now. Our results support the idea that mantle rheology is weak in the region.     

Poroelastic relaxation and aftershocks of the 2001 Bhuj earthquake, India

We analyse aftershocks of the 26 January 2001 Bhuj earthquake, India, that were recorded for 10 weeks following the mainshock. We calculate undrained or instantaneous pore pressure and change in Coulomb stress due to the earthquake and their poroelastic relaxation in the following 10 weeks period. Almost all aftershocks occurred in the region of coseismic dilatation. In the subsequent period, pore pressure increased through relaxation in the dilatation region which further modified coseismic Coulomb stress. Maximum increase in pore pressure is estimated to be about 0.7 MPa in 60 days time following the mainshock. Correlation between the zones of increased pore pressure and postseismic Coulomb stress with that of aftershocks, suggests a definite role of fluid diffusion in their delayed triggering.     

Isoseismals for the Kutch earthquake of 26th January 2001

Damage observed during the Kutch earthquake of 26th January 2001 is classified according to the type of structure. This damage pattern has been used to draw the isoseismal map on the basis of the intensity scale given in IS 1893-2002 which covers only traditional one to two storey houses. Having thus estimated intensities, and observed damage to other types of structures, a new criterion was adopted which is included in this paper. Areas that showed maximum damage to structures have been assigned maximum intensity X.     

Temporal change in the structural fabric during pre- and post-26 January 2001 Bhuj earthquake,using remote sensing data

The Bhuj region, Kutch, India, is included in the highly seismic zonation map of India. The Kutch is a rift basin and so far has experienced three major earthquakes that are due to reverse fault mechanism, which in turn have been ascribed to compressive stresses. Origin of these stresses is considered to be due to north–south convergence of the Indian Plate with the Tibetan plate, and this has placed the entire Indian Plate under the compressive stress regime. Analysis of the stress pattern in the Bhuj region, therefore, has been carried out by extracting lineaments with the help of remote sensing data for the pre- and post-earthquake periods of 26 January 2001 earthquake. For this purpose, the area has been segmented into four sectors. The lineament frequency and the percent frequency from each sector and also for the whole area have been worked out. Resolution of stress on the principle of triaxial ellipsoid has been worked out for each sector and also for the whole area. There results a temporal change in the stress pattern in each sector and also for the whole area. However, the direction of horizontal maximum compressive stress for the whole area appears to be in N 10°E in the pre-earthquake period that has changed to N 10°W in the post-earthquake period. Thus, temporal change in the horizontal maximum compressive stress direction as N 23°E, inferred by Gowd et al. (J Geophy Res 97:11879–11888, 1992) to N 10°E prior to and N 10°W in the post-earthquake period, as inferred from lineament analysis and near parallelism of the lineament maxima with that of the North Kathiawar Fault and the Chambal Jamnagar Lineament along with the longer axis of the isoseismals of the Bhuj 2001 earthquake indicates a modification in the structural fabric of the region as well as a possibility of development of a major plane of weakness.     

Source parameters of the 2001 Mw 7.7 Bhuj earthquake,Gujarat, India,aftershock sequence

We present the estimated source parameters from SH-wave spectral modeling of selected 463 aftershocks (2002–06) of the 26 January 2001 Bhuj earthquake, the well-recorded largest continental intraplate earthquake. The estimated seismic moment (M_o), corner frequency (f_c), source radius (r) and stress drop (Δσ) for aftershocks of moment magnitude 1.7 to 5.6 range from 3.55×10¹¹ to 2.84×10¹⁷ N-m, 1.3 to 11.83 Hz, 107 to 1515 m and 0.13 to 26.7 MPa, respectively, while the errors in f_c and Δσ are found to be 1.1 Hz and 1.1 MPa, respectively. We also notice that the near surface attenuation factor (k) values vary from 0.02 to 0.03. Our estimates reveal that the stress drop values show more scatter (Mo^{0.5 to 1} is proportional to Δσ) toward the larger M_o values (≥10^14.5 N-m), while they show a more systematic nature (Mo³ is proportional to Δσ) for smaller M_o values (<10^14.5 N-m), which can be explained as a consequence of a nearly constant rupture radius for smaller aftershocks in the region. The large stress drops (= 10 MPa) associated with events on the north Wagad fault (at 15–30 km depth) and Gedi fault (at 3–15 km depth) can be attributed to the large stress developed at hypocentral depths as a result of high fluid pressure and the presence of mafic intrusive bodies beneath these two fault zones.     

Evidence for non-self-similarity in the M w 7.7 2001 Bhuj earthquake sequence

Natural Hazards - Scaled energy, apparent stress, seismic moments, stress drops and corner frequencies are measured through the Levenberg–Marquardt nonlinear inversion modeling of S-wave...     

Triggered/migrated seismicity due to the 2001 M w7.7 Bhuj earthquake, Western India

Paper describes triggered seismicity to 200?km distance and for a decade due to the 2001 M _w7.7 Bhuj earthquake. The Kachchh region is seismically one of the most active intraplate regions of the World due to the occurrence of two large earthquakes 1819 (M _w7.8) and 2001 (M _w7.7). Though, it has high hazard but was known to have low seismicity in view of the occurrence of fewer smaller shocks. However, the status seems to have changed after 2001. Besides the strong aftershock activity for over a decade, seismicity has spread to nearby faults in Kachchh peninsula and at several places southward for 200?km distance in Saurashtra peninsula. Beyond the rupture zone of the 2001 Bhuj earthquake, more than 40 mainshocks of M _w?~?3?C5 have occurred at 20 different locations, which is unusual. The increased seismicity is inferred to be caused by stress perturbation due to the 2001 Bhuj earthquake by viscoelastic process. In Saurashtra, over and above the viscoelastic stress increase, the transient stress increase by water table rise in monsoons seems to be affecting the timing of mainshocks and associated sequences of earthquakes.     

Tectonics and crustal structures related to Bhuj earthquake of January 26, 2001: based on gravity and magnetic surveys constrained from seismic and seismological studies

Gravity and magnetic data of the Kachchh basin and surrounding regions have delineated major E–W and NW–SE oriented lineaments and faults, which are even extending up to plate boundaries in the north Arabian Sea and western boundary of the Indian plate, respectively. The epicentral zone of Bhuj earthquake and its aftershocks is located over the junction of Rann of Kachchh and median uplifts viz. Kachchh mainland and Wagad uplifts, which are separated by thrust faults. Gravity data with constraints from the results of the seismic studies along a profile suggest that the basement is uplifted towards the north along thrust faults dipping 40–60° south. Similarly gravity and magnetic modeling along a profile across Wagad uplift suggest south dipping (50–60°) basement contacts separating rocks of high susceptibility and density towards the north. One of these contacts coincides with the fault plane of the Bhuj earthquake as inferred from seismological studies and its projection on the surface coincides with the E–W oriented north Wagad thrust fault. A circular gravity high in contact with the fault in northern part of the Wagad uplift along with high amplitude magnetic anomaly suggests plug type mafic intrusive in this region. Several such gravity anomalies are observed over the island belt in the Rann of Kachchh indicating their association with mafic intrusions. The contact of these intrusives with the country rock demarcates shallow crustal inhomogeneities, which provides excellent sites for the accumulation of regional stress. A regional gravity anomaly map based on the concept of isostasy presents two centers of gravity lows of −11 to −13 mGal (10⁻⁵ m/s²) representing mass deficiency in the epicentral region. Their best-fit model constrained from the receiver function analysis and seismic refraction studies suggest crustal root of 7–8 km (deep crustal inhomogeneity) under them for a standard density contrast of −400 kg/m³. It is, therefore, suggested that significant amount of stress get concentrated in this region due to (a) buoyant crustal root, (b) regional stress due to plate tectonic forces, and (c) mafic intrusives as stress concentrators and the same might be responsible for the frequent and large magnitude earthquakes in this region including the Bhuj earthquake of January 26, 2001.     

Performance of elevated tanks in<Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> 7.7 Bhuj earthquake of january 26th, 2001

The current designs of supporting structures of elevated water tanks are extremely vulnerable under lateral forces due to an earthquake and the Bhuj earthquake provided another illustration when a great many water tank stagings suffered damage and a few collapsed. The more popular shaft type stagings suffer from poor ductility of thin shell sections besides low redundancy and toughness whereas framed stagings consist of weak members and poor brace-column joints. A strength analysis of a few damaged shaft type stagings clearly shows that all of them either met or exceeded the strength requirements of IS:1893-1984, however, they were all found deficient when compared with requirements of the International Building Code. IS:1893-1984 is unjustifiably low for these systems which do not have the advantage of ductility and redundancy and are currently being underestimated at least by a factor of 3 and need an upward revision of forces immediately.     

Estimation of static stress changes after the 2001 Bhuj earthquake: Implications towards the northward spatial migration of the seismic activity in Kachchh,Gujarat

Spatial-temporal patterns of aftershocks of the 2001 Mw7.7 Bhuj earthquake during 2001–2008 reveal a northward spatial migration of seismic activity in the Kachchh seismic zone, which could be related with the loading stresses caused by the continued occurrences of aftershocks on the north Wagad fault (NWF), the causative fault of the 2001-mainshock. Aiming at explaining the observed northward migration of activity, we modelled the Coulomb failure stress change (DCFS) produced by the 2001-mainshock, the 2006 Mw5.6 Gedi fault (GF) and the 2007 Mw4.5 Allah bund fault (ABF) events on optimally oriented plane. A strong correlation between occurrences of earthquakes and regions of increased DCFS is obtained on the associated three faults i.e. NWF, ABF and GF. Predicted DCFS on the GF increased by 0.9 MPa at 3 km depth, where the 7^th March 2006 Mw5.6 event occurred, whereas predicted DCFS on the ABF increased by 0.07 MPa at 30 km depth, where the 15^th December 2007 Mw4.5 event occurred. Focal mechanism solutions of three events on the ABF have been estimated using the iterative inversion of broadband data from 5–10 stations, which are also constrained by the first P-motion data from 8–12 stations. These focal mechanism solutions for the ABF events reveal a dominant reverse movement with a strike-slip component along a preferred northwest or northeast dipping plane (∼50–70°). Focal mechanisms of the events on all the three fault zones reveal an N-S oriented P- axis or maximum principal stress in the region, which agrees with the prevailing N-S compression over the Indian plate. It is apparent that the northward migration of the static stress changes from the NWF, resulting from the occurrence 2001 Bhuj mainshock, might have caused the occurrence of the events on the GF and ABF during 2006–08.     

A simulation of earthquake induced undrained pore pressure changes with bearing on some soil liquefaction observations following the 2001 Bhuj earthquake

The Bhuj earthquake of January 26th, 2001, induced wide spread liquefaction within the Kachch peninsula. It has been pointed out that inundation due to soil liquefaction was short lived in some parts than in others in the affected region. Several geological, seismological and hydrological factors would have cumulatively contributed to these observed changes. We simulate in this article, undrained or short-term change in pore pressure in a poroelastic half space, in response to a simplified model of the Bhuj earthquake source. We find that the regions of relatively shorter lived inundation due to soil liquefaction may fall in the region where pore pressure responsible for soil liquefaction attributable to strong ground shaking was counteracted by pore pressure changes due to undrained poroelastic effect and vice versa.     

An eyewitness account of the Bhuj earthquake

The occurrence of a severe earthquake is a rare event with its effect localized in a limited region. There are no prior indications of its occurrence too; hence experiencing such an event is just a matter of chance, which the author had by virtue of his posting at Bhuj. This paper presents a detailed account of observations made in the wake of the Bhuj earthquake of January 26th 2001, describing physical and mental reactions during the earthquake, the post-quake scenario, the nature and cause of damage to buildings, the trend of aftershocks, various deformities, including ruptures and fissures on the surface of the earth, etc. which may be useful for a detailed study of the seismological activity in the region.