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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Rupture mechanism of the 1993 Killari earthquake, India: constraints from aftershocks and static stress change

The Killari earthquake of September 29, 1993 (Mw=6.2) in peninsular India triggered several aftershocks that were recorded by a network of 21 stations. We computed the change in regional static stress caused by coseismic slip on the earthquake rupture and correlated it with the aftershocks with a view to constrain some of the rupture parameters of this earthquake. We evaluated the six available estimates of fault plane solutions for this earthquake and concluded that reverse slip on a 42° dipping, N112° trending fault, which extends up to the surface from a depth of 7 km, produces maximum correlation between the increased static stress and aftershock distribution. Our analysis suggests that the majority of coseismic slip occurred on the part of the rupture that lies in the depth range of 3–6.5 km.     

An insight into crack density,saturation rate,and porosity model of the 2001 Bhuj earthquake in the stable continental region of western India

The 2001 Bhuj earthquake (Mw 7.6) source zone is examined in the light of crack density (ε), saturation rate (ξ) and porosity parameter (ψ) using new data set derived from a large aftershock sequence recorded by the Gujarat seismic network (GSNet) during November, 2006–December, 2009. Processes of rupture initiations of the mainshock and its aftershock sequence are better understood by synthesizing the dynamic snapshots of the source zone using the new dataset. Pattern of crustal heterogeneities associated with high-ε, high-ξ and high-ψ anomalies at depths varying from 20 km to 25 km is similar to those of earlier study by Mishra and Zhao (2003). The anomalous zone is found extended distinctly by 50–60 km in the lateral direction, indicating the reinforcement of cracks and fractured volume of rock matrix due to long aftershock sequence since 2001 Bhuj earthquake in the source area. It is inferred that the presence of a fluid-filled fractured rock matrix with super saturation may have affected the structural and seismogenic strengths of the source zone and is still contributing significantly to the geneses of earthquakes in and around the source zone. Anomalous pattern of high-ε with wider distribution of high-ξ indicates the existence of micro-cracks in the lower crust, while high-ψ suggests the cementation of cracks through permeation of residual magma/metamorphic fluids into the hypocenter zone. The results suggest that the existence of residual fluids in the fractured rock matrix in the mid to lower crust might have played a key role in triggering the 2001 mainshock and is still responsible for its continued long aftershock sequences.     

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.     

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.     

Local site effects in Latur, India: an empirical study based on the aftershocks of the Killari earthquake of September 29, 1993

The site amplification is estimated at five seismic stations of the Latur region using the horizontal to vertical spectral ratios of 33 aftershocks of the main Killari earthquake of September 29, 1993 (UTC). Spectral amplifications, ranging from a factor of 2–6 are found to vary with frequency at different places. Significant amplification is found at four sites within the Latur region, at Basavakalyan, Kasgi, Killari, and Mudgad Eakoji villages. Our results show a positive correlation between the site amplification and the damage pattern in area. The pattern and the nature of the site amplification estimated in the present study corroborates also with the analytical models and the borehole data indicating alternating layers of unconsolidated sediments and basaltic rocks.     

Quaternary fault segmentation and interaction in the epicentral area of the 1561 earthquake (Mw = 6.4), Vallo di Diano, southern Apennines, Italy

The main structural characteristics of the Caggiano and Polla faults, exposed in the epicentral area of the 1561 earthquake (Mw = 6.4), southern Italy, have been investigated in detail to assess their spatial and temporal properties, and to evaluate their seismogenic potential. These right stepping normal faults show an overlap of about 7 km and an across strike separation of about 4 km. The geometric relationships between the Caggiano and Polla faults, but also the displacement distribution along each fault, demonstrate that they have been strongly interacting throughout the Pleistocene. Nevertheless, geological evidence of Holocene tectonic activity was mainly recognized along the Caggiano Fault (faulted late glacial deposits) and in the southernmost part of the Polla Fault (faulted deposits of probably Late Pleistocene age). This suggests that the Caggiano Fault can be considered as the most tectonically active fault in the Vallo di Diano Fault System. By calculating Coulomb stress changes, we have constrained modes of mechanical interactions between the two faults in a scenario compatible with the 1561 earthquake. This approach allows us to argue that both the Caggiano and the Polla Faults are probably linked at depth, and part of the same seismogenic structure which may be potentially responsible for composite ruptures with magnitude ≥ 6.5.     

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.     

3-D seismic structure of the source area of the 1993 Latur, India, earthquake and its implications for rupture nucleations

The Latur earthquake (M_w 6.1) of 29 September 1993 is a rare stable continental region (SCR) earthquake that occurred on a previously unknown blind fault. In this study, we determined detailed three-dimensional (3-D) P- and S-wave velocity (V_p, V_s) and Poisson's ratio (σ) structures by inverting the first P- and S-wave high-quality arrival time data from 142 aftershocks that were recorded by a network of temporary seismic stations. The source zone of the Latur earthquake shows strong lateral heterogeneities in V_p, V_s and σ structures, extending in a volume of about 90 × 90 × 15 km³. The mainshock occurred within, but near the boundary, of a low-V_p, high-V_s and low-σ zone. This suggests that the structural asperities at the mainshock hypocenter are associated with a partially fluid-saturated fractured rock in a previously unknown source zone with intersecting fault surfaces. This might have triggered the 1993 Latur mainshock and its aftershock sequence. Our results are in good agreement with other geophysical studies that suggest high conductivity and high concentration of radiogenic helium gas beneath the source zone of the Latur earthquake. Our study provides an additional evidence for the presence of fluid related anomaly at the hidden source zone of the Latur earthquake in the SCR and helps us understand the genesis of damaging earthquakes in the SCR of the world.     

Structure and peculiarities of the seismic regime in the source zone of the Takoe earthquake on September 1, 2001 (M W 5.2)

The Takoe earthquake (M _W 5.2) occurred between two en-echelon segments of the active Aprelovskii fault on September 1, 2001, and was accompanied by an earthquake swarm, which was successfully recorded by a local network of digital seismic stations located on the southern part of Sakhalin Island. Modern methods were applied to relocate the parameters of the sources for the earthquake swarm event and significantly specify their spatial distribution and relations to the structural-geological features of the complex system of interacting faults. New data on the correlation between the source mechanism and the modern geodynamic setting in the southern part of Sakhalin were obtained.     

Sr、He同位素指示河西走廊中西部地震活动性

河西走廊地处青藏活动地块东北部边缘,为祁连山地震构造带和南北地震构造带的交汇部位,历史上曾发生过多次中强以上的破坏性地震。2008年5月12日汶川80级特大地震发生后,河西地震带是否会发生中强震是一个较受关注的问题。文中分析了河西走廊张掖、武威、酒泉盆地地下水中的Sr、He同位素特征。结果显示,这3个地区的地下水经历了3种不同的循环路径,不是源自黑河、石羊河、疏勒河的地表径流,而是祁连山降水在出山前下渗深大断裂后的补给。这表明河西走廊盆地内的断裂带现今活动性仍然很强,玉门、高台、民乐、永昌等位于断裂带交汇处或附近的区域未来发生中、强震的可能性较大,应加强对这些地区地下水的实时监测。     

应用综合物探方法探测新疆萨尔苏克外围及深部金铜矿床

新疆萨尔苏克外围金铜矿床位于新疆阿舍勒大型铜矿床西部,处于多拉纳萨依弧沟间隙构造带与阿舍勒岛弧构造带之间,两个构造单元以F6-1玛尔卡库里深大断裂为界。依据本区成矿规律及地层特点,在萨尔苏克外围采用地面高精度磁测、激电中梯、激电测深、CSAMT、TEM组合勘探,在钻孔定位依据上利用激电测深圈定中浅部矿化体,通过研究CSAMT方法成果在浅部与激电测深成果对应关系,进一步研究CSAMT成果对深部矿化体的反应规律。依据CSAMT方法,大致了解了深部可能存在的隐伏岩体、控矿构造、硅化带等的规模和产状变化情况,从而提供深部找矿信息。在萨尔苏克外围Ⅱ号矿化蚀变带上布设的两个钻探验证孔均见铜矿体,预示萨尔苏克外围及深部找矿潜力巨大。     

Finite element modelling of elastic intraplate stresses due to heterogeneities in crustal density and mechanical properties for the Jabalpur earthquake region,central India

Deep lower crustal intraplate earthquakes are infrequent and the mechanism of their occurrence is not well understood. The Narmada-Son-lineament region in central India has experienced two such events, the 1938 Satpura earthquake and the 1997 Jabalpur earthquake, having a focal depth of more than 35 km. We have estimated elastic stresses due to the crustal density and mechanical properties heterogeneities along the Hirapur-Mandla profile passing through the Jabalpur earthquake region to analyse conditions suitable for the concentration of shear stresses in the hypocentral region of this earthquake. Elastic stresses have been computed by a finite element method for a range of material parameters. The results indicate that the shear stresses generated by the density heterogeneities alone are not able to locally enhance the stress concentration in the hypocentral region. The role of mechanical properties of various crustal layers is important in achieving this localization of stresses. Among a range of material parameters analysed, the model with a mechanically strong lower crust overlying a relatively weak sub-Moho layer is able to enhance the stress concentration in the hypocentral region, implying a weaker mantle in comparison to the lower crust for this region of central India.