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.     

Simultaneous inversion of the aftershock data of the 1993 Killari earthquake in Peninsular India and its seismotectonic implications

The aftershock sequence of the September 30th, 1993 Killari earthquake in the Latur district of Maharashtra state, India, recorded by 41 temporary seismograph stations are used for estimating 3-D velocity structure in the epicentral area. The local earthquake tomography (LET) method of Thurber (1983) is used. About 1500P and 1200S wave travel-times are inverted. TheP andS wave velocities as well asV _P/V_Sratio vary more rapidly in the vertical as well as in the horizontal directions in the source region compared to the adjacent areas. The main shock hypocentre is located at the junction of a high velocity and a low velocity zone, representing a fault zone at 6–7 km depth. The estimated average errors ofP velocity andV _P/V_Sratio are ±0.07 km/s and ±0.016, respectively. The best resolution ofP and S-wave velocities is obtained in the aftershock zone. The 3-D velocity structure and precise locations of the aftershocks suggest a ‘stationary concept’ of the Killari earthquake sequence.     

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.     

Focal mechanism dependence of static stress triggering of earthquakes

We perform a global statistical investigation into the problem of Coulomb stress triggering of earthquakes by using the Harvard CMT catalogue of shallow earthquakes from 1976 to 2001. We consider ‘earthquake pairs’, i.e., pairs of successive earthquakes occurring near each other with similar focal mechanisms, and address the problem of whether or not the change of the Coulomb failure stress (CFS) produced by the first earthquake of the pair ‘encourages’ the occurrence of the second one. An interesting feature is that such a Coulomb stress-triggering phenomenon has an apparent focal mechanism dependence: for thrust earthquakes, a more evident triggering effect can be observed.     

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.     

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.     

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.     

A study of the 2 December 2002 (M5.5) Vartholomio (western Peloponnese, Greece) earthquake and of its largest aftershocks

Broadband data from the Greek National Seismological Network are used to study the moderate size (M5.5) earthquake, which occurred on 2 December 2002 near the town of Vartholomio, in western Peloponnese (Greece). Time domain moment tensor inversion applied to retrieve the focal mechanism of the mainshock and of three of the larger aftershocks of the sequence, revealed almost pure strike-slip faulting along NW–SE or NE–SW trending nodal planes. The relative source time functions for the mainshock, obtained from an empirical Green's function analysis, do not reveal any clear directivity to any of the stations. A careful observer might suggest directivity towards NW, if any. Optimum values are 0.4 s for the rise time and 2.7 km/s for the rupture velocity. The spatial and temporal distribution of fault slip showed that the major part of the resolved slip occurred beneath the mainshock's epicenter, 20 km underneath the western coast of Peloponnese. This probably accounts for the considerable damage observed to the nearby towns. The resolution between the two nodal planes does not permit an identification of the fault plane; however the statistics on the slip distribution model, the preliminary analysis of aftershock locations and macroseismic data favour the NW–SE trending plane as the fault plane, which is connected with sinistral strike-slip motions. These are the first implications for sinistral strike-slip motions in this area and more data are needed in the future to get more reliable resolution of the motions.     

First observation of microspherule from the infratrappean Gondwana sediments below Killari region of Deccan LIP,Maharashtra (India) and possible implications

A rare occurrence of a microspherule has been found in the infratrappean sediments, encountered below 338 m thick Deccan volcanic cover in KLR-1 scientific borehole, drilled in the epicentral zone of the 1993 Killari earthquake (Maharashtra, India). Palynological studies of the sediments indicate their age as Early Permian (Asselian, 298–295 Ma) for deposition. Transmission electron microscope studies reveal that the spherule from the infratrappeans, is having a similar composition to that of the Neoarchean amphibolite to granulite facies mid crustal basement. The spherule is non-spherical in nature, containing mostly FeO (10.70 ± 0.20 wt.%), CaO (13.8 ± 0.5 wt.%), Al₂O₃ (7.78 ± 0.30 wt.%), MgO (6.47 ± 0.3 wt.%), SiO₂ (47.46 ± 0.50 wt.%), TiO₂ (2.47 ± 0.3 wt.%), K₂O (1.89 ± 0.20 wt.%), and Cl (0.33 ± 0.05 wt.%). Since the Fe composition of the spherule is almost same as the basement rock (10.5 wt.%), and the chlorine content is also in the same range as the basement (0.04–0.24 wt.%), it would suggest possibility of an extraterrestrial impact over the Indian terrain during the erstwhile Gondwana sedimentation period that may be associated with the Permian–Triassic mass extinction, the most severe one in the Earth's history.     

Tectonic stress field in the epicentral zone of the Latur earthquake of 1993

In situ stress measurements by hydraulic fracturing were carried out in the 617 m deep borehole specially drilled in the epicentral zone of the 1993 Latur earthquake for the purpose of research. The stress measurements carried out at 592 m depth in this borehole are the deepest of all such measurements made so far in the Indian shield. The maximum and minimum principal horizontal stresses (S _H max andS _h min) have been derived from the hydrofracture data using the classical method. TheS _{H max} andS _h min are found to be 16.5 and 9.6 MPa at 373 m depth, and 25.0 and 14.1 MPa at 592 m depth, indicating that the vertical gradients ofS _hmax andS _hmin in the epicentral zone are 39 MPa/km and 21 MPa/km respectively. The principal horizontal stresses in the epicentral zone are comparable with those at Hyderabad and 30% higher than in most other comparable intra-continental regions. Analysis of the results indicate that the stresses in the focal region of the 1993 Latur earthquake have not undergone any significant change following its occurrence and this is in agreement with a similar inference drawn from the seismic data analysis. It appears that the Latur earthquake was caused due to rupturing of the overpressured fault segment at the base of the seismogenic zone.     

静钻根植竹节桩承载力及荷载传递机制研究

静钻根植竹节桩是由预应力竹节桩和桩周水泥土构成的一种新型组合桩基,该新型桩不仅承载性能较好,而且可以大量减少桩基施工过程中所产生的泥浆污染。通过静钻根植竹节桩和钻孔灌注桩的静荷载对比试验及埋设在竹节桩桩身上的应变计对桩身轴力进行测量,分析了静钻根植竹节桩桩身的轴力分布情况以及侧摩阻力的分布,用有限元软件ABAQUS对静钻根植竹节桩进行三维建模计算,详细地分析这种新型组合桩的荷载传递机制。结合现场试验与模拟计算可以得到:在软土地区,静钻根植竹节桩这种新型组合桩的承载力比普通钻孔灌注桩要高;静钻根植竹节桩桩身变形由预制桩所控制,竹节桩与桩外围水泥土近似变形协调;竹节桩竹节的存在对组合桩承载性能有着极其重要的作用;在软土中静钻根植桩侧摩阻力是灌注桩侧摩阻力的1.05¹.10倍。     

由多个地震震源机制解求福建东南沿海地区平均应力场

基于震源断层面的空间取向和滑动方向,写出力轴张量在地理坐标系中的表达式,进而给出计算平均力轴张量及主值的方法,即通过求解相应的本征方程得到。对使用多个震源机制解中T,B,P轴参数计算平均应力场的正确性,以许忠淮用滑动方向拟合法反演富蕴、唐山地区平均应力场数据进行计算和验证。选择具有地震构造意义地块或地震带内大量地震的震源机制解研究区域平均应力场。根据福建东南沿海地区47次主要地震的震源机制解,给出了该区应力张量的定量分析结果。该方法物理意义明确,使由大量地震震源机制解资料分析构造应力场的方法走向定量化,并为定量探索区域应力场的整体动态变化成为可能。     

苏鲁-大别造山带岩石圈应力场、构造运动特征以及超高压变质带折返机制的研究

苏鲁-大别造山带是华北、华南地块之间的大地构造交界带.本文分析研究了我国东部地区的3000余年的地震活动性,并根据1918～2006年问苏鲁-大别及其周围地区发生的1000余个地震的震源机制解,系统研究了地壳应力场和构造运动的区域特征,探讨了超高压变质带的折返机制.结果表明,苏鲁带以及华北地区受到太平洋板块向欧亚板块俯冲挤压的同时,受到从贝加尔湖经过大华北直到琉球海沟的广阔地域里存在着的统一的、方位为170°的引张应力场的控制.右旋走滑运动是苏鲁地区和郯庐断裂带的现代构造运动的主要特征.地震发生类型多为右旋走滑型或右旋走滑正断层型地震.华南地区构造应力场主要表现为,受北西西向运动的菲律宾板块向欧亚板块碰撞挤压运动所产生的方位为110～120°挤压应力的控制.华北、华南地块之间现代地壳应力场的分界线,西部与秦岭带大致相符,在大别东开始则逐渐偏离大地构造,到黄山附近向东南偏转,在温州附近转为向东延伸,最终穿过东海直至琉球海沟.研究结果还表明,苏鲁-大别超高压变质带的折返运动机制,即致使大量超高压变质岩折返到上地壳或地表的岩石圈应力场背景原因为,中生代以来大华北地区存在着较强扩张应力场的主控作用.此外,岩石圈地幔的蘑菇云构造增强了华北地区扩张应力场及扩张构造运动,导致形成大量深裂谷、裂隙,深部物质上涌.苏鲁及东大别地带处于或接近蘑菇云构造运动发育地区的中心,深部地幔物质的上涌导致并加剧了超高压变质岩折返到上地壳或地表,形成了世界著名的苏鲁-别超高压变质带.     

New constraints on the current stress field and seismic velocity structure of the eastern Yilgarn Craton from mechanisms of local earthquakes

The Yilgarn Craton has hosted some of the largest earthquakes within the Australian continent in the last 100 years. Earthquakes have mainly been studied in the western part of the craton, and are thought to result from the reactivation of Precambrian structures in an E–W compressive regional stress field imposed by plate-scale processes. Here we present moment tensor solutions for three recent moderate-sized earthquakes around the town of Kalgoorlie that are inconsistent with E–W compression, but instead suggest E–W extension in the eastern Yilgarn Craton. Waveforms of earthquakes at Boulder (M_W = 4.0, 20 April 2010), Kalgoorlie (M_W = 4.3, 26 February 2014) and Coolgardie (M_W = 3.9, 31 October 2014) were inverted for moment tensors. All three earthquakes were shallow (centroid depth ≤4 km) normal-faulting events that occurred along roughly N–S-striking planes, either with a steep westward or a relatively shallow eastward dip. The robustness of the retrieved mechanisms has been thoroughly tested, employing different earth models, assuming different locations for the earthquakes and using different period bands for the inversion. The fit of synthetic long-period waveforms to the observations was in all cases substantially improved by assuming a two-layered crust with high S wavespeeds (about 3.9–4 km/s) overlying substantially slower material. Since there is independent evidence from active source profiles for a P velocity increase between the upper and lower crust, a large difference in v_p/v_s ratio between upper and lower crust is the only way to explain both lines of evidence. This vertical contrast could represent a dominance of felsic material in the upper crust, and substantially more mafic material in the lower crust. Taken together, our results also appear to imply that the regional stress field is E–W extensive in the Kalgoorlie area, and possibly for the entire Kalgoorlie Terrane. This is contrary to current assumptions from continent-scale stress modelling. That the orientations of rupture planes roughly align with the regional structural grain could indicate that Archean structures are reactivated in response to the current stress field.     

The 20 May 2016 Petermann Ranges earthquake: centroid location,magnitude and focal mechanism from full waveform modelling

Ground velocity records of the 20 May 2016 Petermann Ranges earthquake are used to calculate its centroid-moment-tensor in the 3?D heterogeneous Earth model AuSREM. The global-centroid-moment-tensor reported a depth of 12?km, which is the shallowest allowed depth in the algorithm. Solutions from other global and local agencies indicate that the event occurred within the top 12?km of the crust, but the locations vary laterally by up to 100?km. We perform a centroid-moment-tensor inversion through a spatiotemporal grid search in 3?D allowing for time shifts around the origin time. Our 3?D grid encompasses the locations of all proposed global solutions. The inversion produces an ensemble of solutions that constrain the depth, lateral location of the centroid, and strike, dip and rake of the fault. The centroid location stands out with a clear peak in the correlation between real and synthetic data for a depth of 1?km at longitude 129.8° and latitude –25.6°. A collection of acceptable solutions at this centroid location, produced by different time shifts, constrain the fault strike to be 304?±?4° or 138?±?1°. The two nodal planes have dip angles of 64?±?5° and 26?±?4° and rake angles of 96?±?2° and 77?±?5°, respectively. The southwest-dipping nodal plane with the dip angle of 64° could be seen as part of a near vertical splay fault system at the end of the Woodroffe Thrust. The other nodal plane could be interpreted as a conjugate fault rupturing perpendicular to the splay structure. We speculate that the latter is more likely, since the hypocentres reported by several agencies, including the Geoscience Australia, as well as the majority of aftershocks are all located to the northeast of our preferred centroid location. Our best estimate for the moment magnitude of this event is 5.9. The optimum centroid is located on the 20?km surface rupture caused by the earthquake. Given the estimated magnitude, the long surface rupture requires only ～4?km of rupture down dip, which is in agreement with the shallow centroid depth we obtained.     

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.     

2008年汶川大地震孕震、同震及震后变形和应力演化全过程的数值模拟

2008年MS 8.0级汶川大地震发生在具有复杂的地质构造背景、 强烈的地表起伏、 不均匀的弹性和黏性结构的龙门山断裂带上.由于震前地震活动性不够强烈且地表构造变形较小,龙门山断裂带的地震危险性在汶川地震之前被低估.从数值模拟的角度,建立黏弹性有限元模型,考虑了初始地形、 重力、构造加载、 黏弹性松弛等因素对2008年...     

Arsenic contamination in groundwater from parts of Ambagarh-Chowki block, Chhattisgarh, India: source and release mechanism

Arsenic contamination in tube-well water in Ambagarh-Chowki block, central India, is restricted to local areas confined within the N-S trending Dongargarh rift zone. Affected areas are preferentially located in acid volcanics, close to shear zones and also in granites. Dug-wells even in severely contaminated areas generally have As concentration ≤10 μg/l. But in Kaurikasa area, several tube-wells and dug-wells are severely polluted. Weathered rocks and soils are also enriched in As from severely contaminated areas. As preferentially occurs in iron-enriched soil and similarly altered biotite, chlorite in granite. As sorbed in hydrated iron oxide (HFO) that preferably occurs in acid-leachable fraction and possibly as coatings on kaolinite, illite and goethite in soil or as coatings and along cleavage traces on weathered biotite and chlorite. Reductive dissolution of HFO released sorbed As to groundwater and enriched it in Fe. Pyrite in volcanic and shear zone rocks, although locally As-bearing is a minor source of As in groundwater.     

Electrical structures of the source area of the 1999 Chi-Chi, Taiwan, earthquake: Spatial correlation between crustal conductors and aftershocks

Wide-band magnetotelluric (MT) data were collected on an east–west profile, approximately perpendicular to the local strike of the Chelungpu thrust, through the hypocentral area of the Chi-Chi earthquake for imaging the seismogenic structure. MT data were then inverted for two-dimensional resistivity models plus best-fitting static shift parameters using a nonlinear conjugate gradient algorithm that minimizes the sum of the normalized data misfits and the smoothness of the model. As shown in the inverted 2D resistivity models presented in this paper, an electrical conductor beside the hypocenter of the Chi-Chi earthquake indicates that deep-crustal fluids may participate in the rupture process of the Chi-Chi earthquake. A striking spatial correlation between the crustal conductor and occurrence of aftershocks beneath the Chelungpu fault suggests a postseismic pore pressure adjustment ongoing after the mainshock. Additionally, the hypocenter exhibits an electrical resistive zone, consistent very well with a predicted compact zone from a crustal deformation and transient fluid flow modeling.     

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.