Self-organized Fractal Seismicity and b Value of Aftershocks of the 2001 Bhuj Earthquake in Kutch (India)

— The devastating intraplate earthquake of M_w 7.7 of 26 January, 2001 took place along the south-dipping reverse fault in the lower crust ( 23 km) of Kutch, Gujarat, India, obliterating some 14,000 people. The aftershock activity has ensued for three years. We analyzed 997 aftershocks of M 3.0 to study the b value and fractal correlation dimensions in time and space. The b value is found to be 0.8 ± 0.03 from the Gutenberg-Richter relation and 0.76 ± 0.02 from the maximum-likelihood, suggesting a typical value for the intraplate region. The spatial correlation is 1.71 ± 0.02, indicating that events are approaching a two-dimensional region. Further, the temporal correlation dimension is estimated to be 0.78 ± 0.02, confirming the structure is mono-fractal in time domain. The depth section of b value shows a peak at 15–35 km depth range coinciding with the maximum occurrence of aftershocks ( 47%), which is inferred as a fluid-filled highly fractured rock matrix with fractures of high density. It will be important to note that tomographic results also suggest a low Vp, low Vs and a large Poissons ratio for the same depth range, further confirming this inference. Additionally, we have studied the variation of D₂^s and b value with time. During the first two months of aftershock activity the results show a marked negative correlation between spatial correlation dimension D₂ (large) and b value (low), indicating the predominance of large events associated with weak clustering. The negative correlation means the stress release along faults of a larger surface area. After two months the fractal dimension (D₂^s) and b value suggests a positive correlation implying more numerous smaller shocks with stress release along faults of a smaller surface area. This would indicate a reduced probability of large magnitude earthquakes due to fragmentation of the fault zone.Acknowledgement. The authors thank Dr. V.P. Dimri, Director, NGRI for his encouragement and kind permission to publish this work. The Department of Science and Technology, New Delhi supported this study.     

Estimation of Source Parameters for the Aftershocks of the 2001 Mw 7.7 Bhuj Earthquake, India

The source parameters for 213 Bhuj aftershocks of moment magnitude varying from 2.16 to 5.74 have been estimated using the spectral analysis of the SH- waveform on the transverse component of the three-componnet digital seismograms as well as accelerograms. The estimated stress drop values for Bhuj aftershocks show more scatter (Mo^{0.5 to 1} ∞ Δσ) toward the larger seismic moment values (log Mo ≥ 10^14.5 N-m, larger aftershocks), whereas, they show a more systematic nature (Mo³ ∞ Δσ) for smaller seismic moment (log Mo < 10^14.5 N-m, smaller aftershocks) values. This size dependency of stress drop has also been seen from the relation between our estimated seismic moment and source radius, however, this size-dependent stress drop is not observed for the source parameter estimates for the other stable continental region earthquakes in India and around the world. The estimated seismic moment (Mo), source radius (r) and stress drop (Δσ) for aftershocks of moment magnitude 2.16 to 5.74 range from 1.95 × 10¹² to 4.5 × 10¹⁷ N-m, 239 to 2835 m and 0.63 to 20.7 MPa, respectively. The near-surface attenuation factor (k) is found to be large of the order of 0.03 for the Kachchh region, suggesting thick low velocity sediments beneath the region. The estimated stress drop values show an increasing trend with the depth indicating the base of seismogenic layer (as characterized by larger stress drop values (>15 MPa)) lying in 22–26km depth range beneath the region. We suggest that the concentration of large stress drop values at 10–36km depth may be related to the large stress/strain associvated with a brittle, competent intrusive body of mafic nature.     

Crustal shear-wave splitting in the epicentral zone of the 2001 Mw 7.7 Bhuj earthquake, Gujarat, India

We report here crustal shear-wave anisotropy, ranging from 1% to 10.76% with an average of 2.4% in the aftershock zone of the 2001 Bhuj earthquake, Gujarat, India, from a study of leading shear-wave polarization directions (LPSDs), which vary on average from NNW–SSE to E–W with a delay of 0.07–0.14 s. The delays in the NNW–SSE to NE–SW directions observed at seven stations, near the seismogenic fault, suggest cracks parallel to the direction of the maximum horizontal regional compressional stress prevailing in the region, suggesting a dilatancy-induced anisotropy resulting from approximately stress-aligned parallel vertical micro-cracks. In contrast, the LPSDs at Ramvav, Rapar and Vondh stations, away from the seismogenic fault, are fault parallel, approximately E–W and almost orthogonal to the stress-aligned polarizations inferred elsewhere. The maximum average time delay of 0.14 s is observed at Lodai, where the fast polarization direction is found to be N338°W. This has been observed from anisotropic poro-elastic (APE) modelling and observations that these are 90° flips in shear-wave polarization, resulting from propagation through micro-cracks containing fluids at critically high pore-fluid pressure surrounding the hypocenter of the 2001 mainshock. The presence of high pore-fluid pressure in the seismogenic fault zone could also explain the observed scatter in shear-wave time delays. Further, the coincidence of the N–S trending intrusive bodies (as inferred from tomographic studies in the area) with the N–S direction of regional maximum horizontal compressional stress supports the interpretation of stress-aligned vertical extensive-dilatant anisotropic (EDA) cracks. The depth distribution of the estimated anisotropy (1–10.76%), b-values and stress drop values suggests an increase at 18–30 km depths, which could be attributed to high pore-fluid pressures resulting from a fluid-filled fractured rock matrix or open micro-cracks (characterized by high crack density and high porosity) coinciding with a low velocity zone (at 18–30 km depths) as delineated from tomographic studies in the area.     

Coda Q in the Kachchh Basin, Western India Using Aftershocks of the Bhuj Earthquake of January 26, 2001

Q_C-estimates of Kachchh Basin in western India have been obtained in a high frequency range from 1.5 to 24.0 Hz using the aftershock data of Bhuj earthquake of January 26, 2001 recorded within an epicentral distance of 80 km. The decay of coda waves of 30 sec window from 186 seismograms has been analysed in four lapse time windows, adopting the single backscattering model. The study shows that Q_c is a function of frequency and increases as frequency increases. The frequency dependent Q_c relations obtained for four lapse-time windows are: Q_c=82 f^1.17 (20–50 sec), Q_c=106 f^1.11 (30–60 sec), Q_c=126f^1.03 (40–70 sec) and Q_c=122f^1.02 (50–80 sec). These empirical relations represent the average attenuation properties of a zone covering the surface area of about 11,000, 20,000, 28,000 and 38,000 square km and a depth extent of about 60, 80, 95, 110 km, respectively. With increasing window length, the degree of frequency dependence, n, decreases marginally from 1.17 to 1.02, whereas Q₀ increases significantly from 82 to 122. At lower frequencies up to 6 Hz, Q_c⁻¹ of Kachchh Basin is in agreement with other regions of the world, whereas at higher frequencies from 12 to 24 Hz it is found to be low.     

Fractal dimension of faults network in the upper silesian coal basin (Poland): Preliminary studies

Fractal analysis of faults network, tremor foci spatial distribution as well as the Gutenberg-Richter relationship could further explain whether the biggest seismic events are connected with recent tectonic activity. Fractality of fault systems geometry, as a first step of the analysis, was tested fro a part of the USCB embodying the main structural units. The cluster analysis and the box counting methods were employed.The calculated fractal dimension of fault network was 1.98 for the whole area yet for considered structural units it was close to 1.6. The results point to similarity of studied fault pattern to river network. Faults within selected tectonic units make separate sets which have a distinct geometry and origin. The value of 1.6 is an upper limit to the fracture geometry of rocks that can be explained on the basis of Griffith energy balance concept.     

An Appraisal of the 2001 Bhuj Earthquake (Mw 7.7, India) Source Zone: Fractal Dimension and b Value Mapping of the Aftershock Sequence

We examined seismic characteristics, b value and fractal dimension of the aftershock sequence of the January 26, 2001 Bhuj earthquake (Mw 7.7) that occurred in the Kutch failed rift basin, western margin of the Stable Continental Region (SCR) of India. A total of about 2,000 events (M?≥?2.0) were recorded within two and a half months, immediately after the main shock. Some 795 events were precisely relocated by simultaneous inversion. These relocated events are used for mapping the frequency-magnitude relation (b value) and fractal correlation dimension (Dc) to understand the seismic characteristics of the aftershocks and the source zone of the main shock. The surface maps of the b value and Dc reveal two distinct tectonic arms or zones of the V-shaped aftershock area, western zone and eastern zone. The b value is relatively higher (~1.6) in the western zone compared to a lower value (~1.4) in the eastern zone. The Dc map also shows a higher value (1.2–1.35) in the western zone compared to a lower Dc (0.80–1.15) in the eastern zone; this implies a positive correlation between Dc and b value. Two cross sections, E–W and N–S, are examined. The E–W sections show similar characteristics, higher b value and higher Dc in the western zone and lower in the eastern zone with depth. The N–S sections across the fault zones, however, show unique features; it imaged both the b and Dc characteristics convincingly to identify two known faults, the Kutch Mainland fault and the South Wagad fault (SWF), one stepping over the other with a seismogenic source zone at depth (20–35?km). The source zone at depth is imaged with a relatively lower b and higher Dc at the ‘fault end’ of the SWF showing a negative correlation. These observations, corroborated with the seismic tomography as well as with the proposed geological/tectonic model, shed a new light to our understanding on seismogenesis of the largest SCR earthquake in India in the recent years.     

Relocation of aftershocks of the 2001 Bhuj earthquake: A new insight into seismotectonics of the Kachchh seismic zone,Gujarat, India

In view of an anomalous crust–mantle structure beneath the 2001 Bhuj earthquake region, double-difference relocations of 1402 aftershocks of the 2001 Bhuj earthquake were determined, using an improved 1D velocity model constructed from 3D velocity tomograms based on data from 10 to 58 three-component seismograph stations. This clearly delineated four major tectonic features: (i) south-dipping north Wagad fault (NWF), (ii and iii) south-dipping south Wagad faults 1 and 2 (SWF₁, SWF₂), and (iv) a northeast dipping transverse fault (ITF), which is a new find. The relocated aftershocks correlate satisfactorily with the geologically mapped and inferred faults in the epicentral region. The relocated focal depths delineate a marked variation to the tune of 12 km in the brittle–ductile transition depths beneath the central aftershock zone that could be attributed to a lateral variation in crustal composition (more or less mafic) or in the level of fracturing across the fault zone. A fault intersection between the NWF and ITF has been clearly mapped in the 10–20 km depth range beneath the central aftershock zone. It is inferred that large intraplate stresses associated with the fault intersection, deepening of the brittle–ductile transition to a depth of 34 km due to the presence of mafic/ultramafic material in the crust–mantle transition zone, and the presence of aqueous fluids (released during the metamorphic process of eclogitisation of lower crustal olivine-rich rocks) and volatile CO₂ at the hypocentral depths, might have resulted in generating the 2001 Bhuj earthquake sequence covering the entire lower crust.     

分形几何在岩土力学研究中的过去、现在与未来

首先简略地介绍了分形几何理论及其方法。从岩土材料结构的定量描述、流体在岩土体中的渗流问题、岩土材料强度的分形模型和分形空间的力学特征四个方面叙述了分形理论在岩土力学研究中取得的成果及应用。指出只有将岩土力学系统在欧氏空间的各种基本变量、原理和方法向分形空间推广和拓展,才能产生对岩土材料非线性力学行为的新认识,这也是分形理论在岩土力学研究中进一步深化的基础。     

Modelling of Strong Ground Motions from 1991 Uttarkashi,India, Earthquake Using a Hybrid Technique

We present a simple and efficient hybrid technique for simulating earthquake strong ground motion. This procedure is the combination of the techniques of envelope function (Midorikawa et al. Tectonophysics 218:287–295, 1993) and composite source model (Zeng et al. Geophys Res Lett 21:725–728, 1994). The first step of the technique is based on the construction of the envelope function of the large earthquake by superposition of envelope functions for smaller earthquakes. The smaller earthquakes (sub-events) of varying sizes are distributed randomly, instead of uniform distribution of same size sub-events, on the fault plane. The accelerogram of large event is then obtained by combining the envelope function with a band-limited white noise. The low-cut frequency of the band-limited white noise is chosen to correspond to the corner frequency for the target earthquake magnitude and the high-cut to the Boore’s f _max or a desired frequency for the simulation. Below the low-cut frequency, the fall-off slope is 2 in accordance with the ω² earthquake source model. The technique requires the parameters such as fault area, orientation of the fault, hypocenter, size of the sub-events, stress drop, rupture velocity, duration, source–site distance and attenuation parameter. The fidelity of the technique has been demonstrated by successful modeling of the 1991 Uttarkashi, Himalaya earthquake (Ms 7). The acceptable locations of the sub-events on the fault plane have been determined using a genetic algorithm. The main characteristics of the simulated accelerograms, comprised of the duration of strong ground shaking, peak ground acceleration and Fourier and response spectra, are, in general, in good agreement with those observed at most of the sites. At some of the sites the simulated accelerograms differ from observed ones by a factor of 2–3. The local site geology and topography may cause such a difference, as these effects have not been considered in the present technique. The advantage of the technique lies in the fact that detailed parameters such as velocity-Q structures and empirical Green’s functions are not required or the records of the actual time history from the past earthquakes are not available. This method may find its application in preparing a wide range of scenarios based on simulation. This provides information that is complementary to the information available in probabilistic hazard maps.     

Relocation of Early and Late Aftershocks of the 2001 Bhuj Earthquake Using Joint Hypocentral Determination (JHD) Technique: Implication toward the Continued Aftershock Activity for more than Four Years

We employed layered model joint hypocentral determination (JHD) with station corrections to improve location identification for the 26 January, 2001 M_w 7.7 Bhuj early and late aftershock sequence. We relocated 999 early aftershocks using the data from a close combined network (National Geophysical Research Institute, India and Center for Earthquake Research Institute, USA) of 8–18 digital seismographs during 12–28 February, 2001. Additionally, 350 late aftershocks were also relocated using the data from 4–10 digital seismographs/accelerographs during August 2002 to December 2004. These precisely relocated aftershocks (error in the epicentral location<30 meter, error in the focal depth estimation < 50 meter) delineate an east-west trending blind thrust (North Wagad Fault, NWF) dipping (~ 45°) southward, about 25 km north of Kachchh main land fault (KMF), as the causative fault for the 2001 Bhuj earthquake. The aftershock zone is confined to a 60-km long and 40-km wide region lying between the KMF to the south and NWF to the north, extending from 2 to 45 km depth. Estimated focal depths suggest that the aftershock zone became deeper with the passage of time. The P- and S-wave station corrections determined from the JHD technique indicate that the larger values (both +ve and -ve) characterize the central aftershock zone, which is surrounded by the zones of smaller values. The station corrections vary from −0.9 to +1.1 sec for the P waves and from −0.7 to +1.4 sec for the S waves. The b-value and p-value of the whole aftershock (2001–2004) sequences of M_w ≥ 3 are estimated to be 0.77 ± 0.02 and 0.99 ± 0.02, respectively. The p-value indicates a smaller value than the global median of 1.1, suggesting a relatively slow decay of aftershocks, whereas, the relatively lower b-value (less than the average b-value of 1.0 for stable continental region earthquakes of India) suggests a relatively higher probability for larger earthquakes in Kachchh in comparison to other stable continental regions of the Indian Peninsula. Further, based on the b-value, mainshock magnitude and maximum aftershock magnitude, the Bhuj aftershock sequence is categorized as the Mogi's type II sequence, indicating the region to be of intermediate level of stresses and heterogeneous rocks. It is inferred that the decrease in p-value and increase in aftershock zone, both spatially as well as depth over the passage of time, suggests that the decay of aftershocks perhaps could be controlled by visco-elastic creep in the lower crust.     

A quantitative study of the energy release in the aftershocks of the Bhuj earthquake, 2001, India, using Lg phase

The devastating earthquake on 26 January 2001 at Bhuj, India, resulted in large-scale death and destruction of properties of several million US dollars. The moment magnitude of the earthquake was 7.7 and its maximum focal intensity exceeded X in MM scale. The rate of aftershocks of this earthquake, recorded at Gauribidanur seismic array station (GBA), shows a monotonic decay with time superposed with oscillations. For the Indian continent the Lg phase is a prominent arrival at regional distances. The estimate of Lg amplitude is obtained by optimally fitting the Lg wave train to a exponential decay curve. The logarithm of these amplitudes and logarithm of root mean square (rms) value of actual amplitudes of the Lg are calibrated with USGS m_b to create a local m_bLg magnitude scale. The energy released from these aftershocks is calculated from the rms value of Lg phase. The plot of cumulative energy release with time follows the power law of the form t^p, superposed with oscillations. The exponent of the power law, p, is estimated both by a time-window scanning method and by an interpolation method. The value of p is 0.434 for time-window scanning method and 0.432 for the interpolation method. The predominant periods found in the oscillatory part of the cumulative energy, obtained by differencing the observed from the power law fit, are 10.6, 7.9, 5.4, 4.6 and 3.5 h for time-window scanning method. The corresponding periods for interpolation method are 13.4, 11.5, 7.4, 4.2, 3.5, 2.6 and 2.4 h.     

Spatial and temporal rupture process of the January 26, 2001, Gujarat, India, M_S=7.8 earthquake

Introduction An earthquake of MS=7.8 occurred near the Gujarat of India on January 26, 2001, which was one of the most deadly earthquakes since there was the record in the Indian history (Bendick, et al, 2001; Gupta, et al, 2001). The USGS of USA determined the origin time of the earthquake to be 3h16min41s (UTC), and the epicenter location to be 70.32篍, 23.40篘. Shortly after the earthquake, the moment tensor solutions or focal mechanisms and other related parameters were offered by s…     

2001年1月26日印度古吉拉特邦7.8级地震综述

2001年1月26日印度古吉拉特邦发生7．8级地震，造成生命和财产的严重损失，这是该地区历史上发生的第二次大地震，引起国际社会的普遍关注，本文综述了这次地震的情况，内容包括：震源参数，发震背景，震灾情况，救援情况，地震造成的社会和经济影响以及值得思考的几个问题。     

Spatial and temporal rupture process of the January 26, 2001, Gujarat, India, M S=7.8 earthquake

The source parameters, such as moment tensor, focal mechanism, source time function (STF) and temporal-spatial rupture process, were obtained for the January 26, 2001, India, M _S=7.8 earthquake by inverting waveform data of 27 GDSN stations with epicentral distances less than 90°. Firstly, combining the moment tensor inversion, the spatial distribution of intensity, disaster and aftershocks and the orientation of the fault where the earthquake lies, the strike, dip and rake of the seismogenic fault were determined to be 92°, 58° and 62°, respectively. That is, this earthquake was a mainly thrust faulting with the strike of near west-east and the dipping direction to south. The seismic moment released was 3.5×10²⁰ Nm, accordingly, the moment magnitude M _W was calculated to be 7.6. And then, 27 P-STFs, 22 S-STFs and the averaged STFs of them were determined respectively using the technique of spectra division in frequency domain and the synthetic seismogram as Green’s functions. The analysis of the STFs suggested that the earthquake was a continuous event with the duration time of 19 s, starting rapidly and ending slowly. Finally, the temporal-spatial distribution of the slip on the fault plane was imaged from the obtained P-STFs and S-STFs using an time domain inversion technique. The maximum slip amplitude on the fault plane was about 7 m. The maximum stress drop was 30 MPa, and the average one over the whole rupture area was 7 MPa. The rupture area was about 85 km long in the strike direction and about 60 km wide in the down-dip direction, which, equally, was 51 km deep in the depth direction. The rupture propagated 50 km eastwards and 35 km westwards. The main portion of the rupture area, which has the slip amplitude greater than 0.5 m, was of the shape of an ellipse, its major axis oriented in the slip direction of the fault, which indicated that the rupture propagation direction was in accordance with the fault slip direction. This phenomenon is popular for strike-slip faulting, but rather rare for thrust faulting. The eastern portion of the rupture area above the initiation point was larger than the western portion below the initiation point, which was indicative of the asymmetrical rupture. In other words, the rupturing was kind of unilateral from west to east and from down to up. From the snapshots of the slip-rate variation with time and space, the slip rate reached the largest at the 4th second, that was 0.2 m/s, and the rupture in this period occurred only around the initiation point. At the 6th second, the rupture around the initiation point nearly stopped, and started moving outwards. The velocity of the westward rupture was smaller than that of the eastward rupture. Such rupture behavior like a circle mostly stopped near the 15th second. After the 16th second, only some patches of rupture distributed in the outer region. From the snapshots of the slip variation with time and space, the rupture started at the initiation point and propagated outwards. The main rupture on the area with the slip amplitude greater than 5 m extended unilaterally from west to east and from down to up between the 6th and the 10th seconds, and the western segment extended a bit westwards and downwards between the 11th and the 13th seconds. The whole process lasted about 19 s. The rupture velocity over the whole rupture process was estimated to be 3.3 km/s. Foundation item: 973 Project (G1998040705) from Ministry of Science and Technology, P. R. China, and the National Science Foundation of China under grant No.49904004. Contribution No. 02FE2026, Institute of Geophysics, China Seismological Bureau.     

Estimation of Site Response in Kachchh,Gujarat, India,Region Using H/V Spectral Ratios of Aftershocks of the 2001 <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> 7.7 Bhuj Earthquake

Site response in the aftershock zone of 2001 Bhuj M_w 7.7 earthquake has been studied using the H/V spectral ratio method using 454 aftershocks (M_w 2.5–4.7) recorded at twelve three-component digital strong motion and eight three-component digital seismograph sites. The mean amplification factor obtained for soft sediment sites (Quaternary/Tertiary) varies from 0.75–6.03 times for 1–3 Hz and 0.49–3.27 times for 3–10 Hz. The mean amplification factors obtained for hard sediment sites (hard Jurassic/Mesozoic sediments) range from 0.32–3.24 times for 1–3 Hz and 0.37–2.18 times for 310 Hz. The upper bounds of the larger mean amplification factors for 1–3 Hz are found to be of the order of 3.13–6.03 at Chopadwa, Vadawa, Kavada, Vondh, Adhoi, Jahwarnagar and Gadhada, whereas, the upper bounds of the higher mean amplification factors at 3–10 Hz are estimated to be of the order of 2.00–3.27° at Tapar, Chopadwa, Adhoi, Jahwarnagar, Gandhidham and Khingarpur. The site response estimated at Bhuj suggests a typical hard-rock site behavior. Preliminary site response maps for 1–3 Hz and 310 Hz frequency ranges have been prepared for the area extending from 23–23.85 °N and 69.65–70.85°E. These frequency ranges are considered on the basis of the fact that the natural frequencies of multi-story buildings (3 to 10 floor) range between 1–3 Hz, while the natural frequencies for 1 to 3 story buildings vary from 3–10 Hz. The 1–3 Hz map delineates two distinct zones of maximum site amplification (>3 times): one lying in the NW quadrant of the study area covering Jahwarnagar, Kavada and Gadadha and the other in the SE quadrant of the study area with a peak of 6.03 at Chopadwa covering an area of 70 km × 50 km. While the 3–10 Hz map shows more than 2 times site amplification value over the entire study area except, NE quadrant, two patches in the southwest corner covering Bhuj and Anjar, and one patch at the center covering Vondh, Manfara and Sikara. The zones for large site amplification values (∼3 times) are found at Tapar, Chopadwa, Adhoi and Chobari. The estimated site response values show a good correlation with the distribution of geological formations as well as observed ground deformation in the epicentral zone.     

Automatic recognition of damaged town buildings caused by earthquake using remote sensing information: Taking the 2001 Bhuj, India, earthquake and the 1976 Tangshan, China, earthquake as examples

Introduction As we well know, the hazard of earthquake is very wide especially in cities. The conventionalmethods to investigate the damage are difficult to meet the requirements in applications. In recentyears, with the rapid development of remote sensing, especially the successful launch and applica-tion of high-resolution commercial remote sensing satellite, it has become possible to recognize andextract damage information by using remote sensing. The researchers at home and abroad hav…     

Simulation of the damage to a simple RC structure in the Bhuj,India, earthquake to estimate the level of input‐motion

Analysis of a simple reinforced concrete (RC) structure damaged by the Bhuj, India, earthquake was carried out to estimate the level of shaking in the epicentral region. For this, an attempt was made to estimate the level of input motion to cause inelastic behavior to the extent observed during the field visit. To consider the inelastic effects, both yielding of steel bars as well as crushing of the concrete cover has been investigated employing the hysteretic model known as the Fiber model. The only available record at Ahmedabad of the Bhuj earthquake and four additional earthquake records from Japan and California were used in the analysis. Considering simple scaling of input motion, the level of input motion to cause crushing and spall‐off of the concrete cover as observed in the field was estimated to be of the order of 6 times the original instrumental record obtained at 240 km away from the epicenter. The methodology proposed was promising in providing a useful quantitative indication of the level of shaking when instrumental records are not available. It was also noted that the design response spectrum specified in Indian seismic code IS1893: 1984 appears inadequate compared to the extent of shaking estimated. Copyright © 2005 John Wiley & Sons, Ltd.     

Instability structures, synsedimentary faults and turbidites, witnesses of a Liassic seismotectonic activity in the Dauphiné Zone (French Alps): A case example in the Lower Pliensbachian at Saint-Michel-en-Beaumont

Instability structures, synsedimentary faults and turbidites have been studied in the Lower Pliensbachian succession of Saint-Michel-en-Beaumont, belonging to the Taillefer block, an ancient half-graben emplaced during the Liassic Tethyan rifting. Geometrical and mechanical analyses demonstrate that the instability structures occurred thanks to movements along spineless synsedimentary normal faults, when the turbiditic and limestone layers were already case-hardened and partly fractured by tension gashes even when the mudstones were still unlithified. Both the tension gashes and the synsedimentary faults are homogeneous in strike with the major regional faults and are in good agreement with the regional direction of extension for this period. The characters of the turbiditic beds, with erosive base, graded bedding, and incomplete Bouma sequence, are in favour of a seismic origin. Instability structures, spineless synsedimentary faults and turbiditic inflows are thus considered as seismites and interpreted as the result of high seismicity periods including some events with M > 5 in the general extensive ambiance of the Liassic Tethyan rifting. The analysis of the geometrical relationships between all these sedimentary features allows to distinguish the successive stage of occurrence of an instability structure, from the sedimentation of alternating marls and limestones, and sudden turbiditic inflows, then early case-hardening of the turbidites, until the important seismotectonic event generating the spineless normal faults, themselves triggering the fall of indurated blocks and locally the forming of breccias. The Ornon Fault, which constitutes the border of the Taillefer block, 15 km eastward, played a major role during the Liassic sedimentation and may represent the major seismic fault related to the seismites occurrence in the Beaumont basin.     

Estimation of Coda Wave Attenuation for the National Capital Region, Delhi, India Using Local Earthquakes

Attenuation of seismic waves is very essential for the study of earthquake source parameters and also for ground-motion simulations, and this is important for the seismic hazard estimation of a region. The digital data acquired by 16 short-period seismic stations of the Delhi Telemetric Network for 55 earthquakes of magnitude 1.5 to 4.2, which occurred within an epicentral distance of 100 km in an area around Delhi, have been used to estimate the coda attenuation Q_c. Using the Single Backscattering Model, the seismograms have been analyzed at 10 central frequencies. The frequency dependence average attenuation relationship Q_c = 142f ^1.04 has been attained. Four Lapse-Time windows from 20 to 50 seconds duration with a difference of 10 seconds have been analyzed to study the lapse time dependence of Q_c. The Q_c values show that frequency dependence (exponent n) remains similar at all the lapse time window lengths. While the change in Q₀ values is significant, change in Q₀ with larger lapsetime reflects the rate of homogeneity at the depth. The variation of Q_c indicates a definitive trend from west to east in accordance with the geology of the region.     

Tectonic accretion of a subducted intraoceanic remnant arc in Cretaceous Hokkaido, Japan, and implications for evolution of the Pacific northwest

Abstract   An accretionary complex, which contains fragments of a remnant island arc, was newly recognized in the Cretaceous accretionary terranes in Hokkaido, Japan. It consists of volcanics, volcanic conglomerate, intermediate to ultramafic intrusive rocks with island-arc affinity including boninitic rocks, accompanied by chert and deformed terrigenous turbidites. Compared with the results of modern oceanic surveys, the preserved sequence from island-arc volcanics to chert, via reworked volcanics, is indicative of intraoceanic remnant arc, because the sequence suggests an inactive arc isolated within a pelagic environment before its accretion. The age of a subducting oceanic crust can be discontinuous before and after a remnant-arc subduction, resulting in abrupt changes in accretion style and metamorphism, as seen in Cretaceous Hokkaido. Subduction of such an intraoceanic remnant arc suggests that the subducted oceanic plate in the Cretaceous was not an extensive oceanic plate like the Izanagi and/or Kula Plates as previously believed by many authors, but a marginal basin plate having an arc–back-arc system like the present-day Philippine Sea Plate.