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.     

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.     

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.     

Aftershock activity of Bhuj earthquake of january 26th, 2001

Following a large-sized Bhuj earthquake (M _s = 7.6) of January 26th, 2001, a small aperture 4-station temporary local network was deployed, in the epicentral area, for a period of about three weeks and resulted in the recording of more than 1800 aftershocks (-0.07 ≤M _L <5.0). Preliminary locations of epicenters of 297 aftershocks (2.0 ≤M _L <5.0) have brought out a dense cluster of aftershock activity, the center of which falls 20 km NW of Bhachau. Epicentral locations of after-shocks encompass a surface area of about 50 × 40 km² that seems to indicate the surface projection of the rupture area associated with the earthquake. The distribution of aftershock activity above magnitude 3, shows that aftershocks are nonuniformly distributed and are aligned in the north, northwest and northeast directions. The epicenter of the mainshock falls on the southern edge of the delineated zone of aftershock activity and the maximum clustering of activity occurs in close proximity of the mainshock. Well-constrained focal depths of 122 aftershocks show that 89% of the aftershocks occurred at depths ranging between 6 and 25 km and only 7% and 4% aftershocks occur at depths less than 5 and more than 25 km respectively. The Gutenberg-Richter (GR) relationship, logN = 4.52 - 0.89M_L, is fitted to the aftershock data (1.0<-M _L<5.0) and theb-value of 0.89 has been estimated for the aftershock activity.     

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.     

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.     

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.     

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.     

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.     

Satellite detection of carbon monoxide emission prior to the Gujarat earthquake of 26 January 2001

NOAA AVHRR images have clearly shown anomalous changes in land surface temperature associated with earthquakes in the past two decades. Soon after the Gujarat earthquake of January 26, 2001, an anomalous increase in land surface temperature was inferred from MODIS satellite data a few days prior to the main earthquake event. The cause of such an anomalous change in surface temperature prior to the earthquake is attributed to many probable phenomena, but no definite cause has been identified. In the present study, changes of a complementary nature were found of land surface temperature associated with the emission of CO from the epicentral region. The observed changes on land and atmosphere associated with the Gujarat earthquake of 26 January, 2001, show the existence of strong coupling between land, atmosphere and ionosphere.     

Intrinsic and scattering attenuation in western India from aftershocks of the 26 January, 2001 Kachchh earthquake

176 vertical-component, short period observations from aftershocks of the M_w 7.7, 26 January, 2001 Kachchh earthquake are used to estimate seismic wave attenuation in western India using uniform and two layer models. The magnitudes (M_w) of the earthquakes are less than 4.5, with depths less than 46 km and hypocentral distances up to 110 km. The studied frequencies are between 1 and 30 Hz. Two seismic wave attenuation factors, intrinsic absorption (Q_i^− 1) and scattering attenuation (Q_s^− 1) are estimated using the Multiple Lapse Time Window method which compares time integrated seismic wave energies with synthetic coda wave envelopes for a multiple isotropic scattering model. We first assume spatial uniformity of Q_i^− 1, Q_s^− 1 and S wave velocity (β). A second approach extends the multiple scattering hypothesis to media consisting of several layers characterized by vertically varying scattering coefficient (g), intrinsic absorption strength (h), density of the media (ρ) and shear wave velocity structure. The predicted coda envelopes are computed using Monte Carlo simulation. Results show that, under the assumption of spatial uniformity, scattering attenuation is greater than intrinsic absorption only for the lowest frequency band (1 to 2 Hz), whereas intrinsic absorption is predominant in the attenuation process at higher frequencies (2 to 30 Hz). The values of Q obtained range from Q_t = 118, Q_i = 246 and Q_s = 227 at 1.5 Hz to Q_t ≈ 4000, Q_i ≈ 4600 and Q_s ≈ 33,300 at 28 Hz center frequencies, being Q_t^− 1 a measure of total attenuation. Results also show that Q_i^− 1, Q_s^− 1 and Q_t^− 1 decrease proportional to f^−ν. Two rates of decay are clearly observed for the low (1 to 6 Hz) and high (6 to 30 Hz) frequency ranges. Values of ν are estimated as 2.07 ± 0.05 and 0.44 ± 0.09 for total attenuation, 1.52 ± 0.21 and 0.48 ± 0.09 for intrinsic absorption and 3.63 ± 0.07 and 0.06 ± 0.08 for scattering attenuation for the low and high frequency ranges, respectively. Despite the lower resolution in deriving the attenuation parameters for a two layered crust, we find that scattering attenuation is comparable to or smaller than the intrinsic absorption in the crust whereas intrinsic absorption dominates in the mantle. Also, for a crustal layer of thickness 42 km, intrinsic absorption and scattering estimates in the crust are lower and greater than those of the mantle, respectively.     

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.     

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.     

The Bhuj, India, earthquake: lessons learned for earthquake safety of dams on alluvium

The Bhuj, India, earthquake of 26 January 2001, M_s 7.9, caused dams built on alluvium to sustain damage ranging from cosmetic to severe. Major damage was caused almost entirely by soil liquefaction in the alluvium. The critical factor was the level of earthquake ground motion.

The Bhuj earthquake showed that peak horizontal accelerations (PHAs)≤0.2 g were generally safe. PHAs>0.2 g were hazardous, when unconsolidated granular foundation soils were water saturated. N values of <20 are indicative of susceptibility to soil liquefaction. The Bhuj experience showed that alluvial foundation soils, subject to a PHA>0.2 g, must be evaluated over the full area beneath a new dam and all soils deemed susceptible to liquefaction must be either removed or treated. For remediating an old dam, reliable options are removal and replacement of liquefiable alluvium beneath upstream and downstream portions of the dam, combined with building berms designed to provide stability for the dam should there be a strength loss in soils beneath the dam.     

Artificial neural network and liquefaction susceptibility assessment: a case study using the 2001 Bhuj earthquake data,Gujarat, India

This study pertains to prediction of liquefaction susceptibility of unconsolidated sediments using artificial neural network (ANN) as a prediction model. The backpropagation neural network was trained, tested, and validated with 23 datasets comprising parameters such as cyclic resistance ratio (CRR), cyclic stress ratio (CSR), liquefaction severity index (LSI), and liquefaction sensitivity index (LSeI). The network was also trained to predict the CRR values from LSI, LSeI, and CSR values. The predicted results were comparable with the field data on CRR and liquefaction severity. Thus, this study indicates the potentiality of the ANN technique in mapping the liquefaction susceptibility of the area.     

Internal geometry of reactivated and non-reactivated sandblow craters related to 2001 Bhuj earthquake,India: A modern analogue for interpreting paleosandblow craters

The liquefaction attributes and crater geometry related to 2001 Bhuj earthquake has been reconstructed by trenching along large known craters formed near Umedpar in Kachchh. The study characterises the liquefied sediments in a large reactivated crater and distinguishes it from a non-reactivated crater located nearby. These characteristics can help in the interpretation of large paleocraters formed as a result of earthquake induced liquefaction.     

Relocation of aftershocks, focal mechanisms and stress inversion: Implications toward the seismo-tectonics of the causative fault zone of Mw7.6 2001 Bhuj earthquake (India)

The HYPODD relocation of 1172 aftershocks, recorded on 8–17 three-component digital seismographs, delineate a distinct south dipping E–W trending aftershock zone extending up to 35 km depth, which involves a crustal volume of 40 km × 60 km × 35 km. The relocated focal depths delineate the presence of three fault segments and variation in the brittle–ductile transition depths amongst the individual faults as the earthquake foci in the both western and eastern ends are confined up to 28 km depth whilst in the central aftershock zone they are limited up to 35 km depth. The FPFIT focal mechanism solutions of 444 aftershocks (using 8–12 first motions) suggest that the focal mechanisms ranged between pure reverse and pure strike slip except some pure dip slip solutions. Stress inversion performed using the P and T axes of the selected focal mechanisms reveals an N181°E oriented maximum principal stress with a very shallow dip (= 14°). The stress inversions of different depth bins of the P and T axes of selected aftershocks suggest a heterogeneous stress regime at 0–30 km depth range with a dominant consistent N–S orientation of the P-axes over the aftershock zone, which could be attributed to the existence of varied nature and orientation of fractures and faults as revealed by the relocated aftershocks.     

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.     

Pre-seismic, co-seismic and post-seismic displacements associated with the Bhuj 2001 earthquake derived from recent and historic geodetic data

The 26th January 2001 Bhuj earthquake occurred in the Kachchh Rift Basin which has a long history of major earthquakes. Great Triangulation Survey points (GTS) were first installed in the area in 1856–60 and some of these were measured using Global Positioning System (GPS) in the months of February and July 2001. Despite uncertainties associated with repairs and possible reconstruction of points in the past century, the re-measurements reveal pre-seismic, co-seismic and post-seismic deformation related to Bhuj earthquake. More than 25 Μ-strain contraction north of the epicenter appears to have occurred in the past 140 years corresponding to a linear convergence rate of approximately 10 mm/yr across the Rann of Kachchh. Motion of a single point at Jamnagar 150 km south of the epicenter in the 4 years prior to the earthquake, and GTS-GPS displacements in Kathiawar suggests that pre-seismic strain south of the epicenter was small and differs insignificantly from that measured elsewhere in India. Of the 20 points measured within 150 km of the epicenter, 12 were made at existing GTS points which revealed epicentral displacements of up to 1 m, and strain changes exceeding 30 Μ-strain. Observed displacements are consistent with reverse co-seismic slip. Re-measurements in July 2001 of one GTS point (Hathria) and eight new points established in February reveal post-seismic deformation consistent with continued slip on the Bhuj rupture zone.