Attenuation of Coda Waves in the Saurashtra Region, Gujarat (India)

The attenuation characteristics based on coda waves of two areas—Jamnagar and Junagarh of Saurashtra, Gujarat (India)—have been investigated in the present study. The frequency dependent relationships have been developed for both the areas using single back scattering model. The broadband waveforms of the vertical components of 33 earthquakes (Mw 1.5–3.5) recorded at six stations of the Jamnagar area, and broadband waveforms of 68 earthquakes (Mw 1.6–5) recorded at five stations of the Junagarh area have been used for the analysis. The estimated relations for the Junagarh area are: Q _c?=?(158?±?5)f^(0.99±0.04) (lapse time : 20?s), Q _c?=?(170?±?4.4)f^(0.97±0.02) (lapse time : 30?s) and Q _c?=?(229?±?6.6)f^(0.94±0.03) (lapse time : 40?s) and for the Jamnagar area are: Q _c?=?(178?±?3)f^(0.95±0.05) (lapse time : 20?s), Q _c?=?(224?±?6)f^(0.98±0.06) (lapse time : 30?s) and Q _c?=?(282?±?7)f^(0.91±0.03) (lapse time : 40?s). These are the first estimates for the areas under consideration. The Junagarh area appears to be more attenuative as compared to the Jamnagar area. The increase in Q _c values with lapse time found here for both the areas show the depth dependence of Q _c as longer lapse time windows will sample larger area. The rate of decay of attenuation (Q ^?1) with frequency for the relations obtained here is found to be comparable with those of other regions of the world though the absolute values differ. A comparison of the coda-Q estimated for the Saurashtra region with those of the nearby Kachchh region shows that the Saurashtra region is less heterogeneous. The obtained relations are expected to be useful for the estimation of source parameters of the earthquakes in the Saurashtra region of Gujarat where no such relations were available earlier. These relations are also important for the simulation of earthquake strong ground motions in the region.     

Attenuation of Broadband P and S Waves in Tonga: Observations of Frequency Dependent Q

—Teleseismic broadband recordings of intermediate and deep focus earthquakes are used to quantify both compression (Q and shear(Q ) wave attenuation within the Lau backarc basin. A spectral-ratio method is employed to measure differential attenuation (t*) between the depth phases sS, pP, and sP and the direct S and P phases over the frequency band 0.05 and 0.5 Hz. We use a stacking algorithm to combine the spectra of several phase pairs from a single event, having similar azimuth and range, to obtain more robust t* measurements; these estimates are then used to compute the average Q above the focal depth. Q and Q are measured directly from the sS-S and pP-P phase pairs respectively, however, the interpretation of t* measured from sP-P requires assumptions about the ratio Q ?/Q . We find an empirical ratio of Q ?/Q = 1.93 for this region and use it to compute Q and Q from the Q _sP observations. We observe lateral and depth variations in both Q _β and Q _>alpha< beneath the tectonically active Lau Basin and the geologically older, inactive Lau Ridge and Fiji Plateau. The upper 200 km beneath the Central and Northern Lau Basin show a Q of 45–57 and a Q of 102–121, and Q appears to increase rapidly with depth. The upper 600 km beneath the Lau backarc basin has a Q of 118–138, while over the same depth interval we observe a higher Q of 139–161 beneath the Lau Ridge and Fiji Plateau. We also find Q of 235–303 beneath the northern Lau Basin and a higher Q of 292–316 beneath the Fiji Plateau and the Lau Ridge measured directly from pP-P phase pairs. These geographic trends in the broadband Q measurements correlate with our previous long-period estimates of Q in this region, however, the broadband measurements themselves are higher by about a factor of two. These observations suggest substantial frequency dependence of Q in the upper mantle, beginning at frequencies less than 1.0 Hz and consistent with the power-law form Q∝ω^α with between -0.1 and -0.3.     

赣南地区P波、S波随频率衰减关系研究

利用发生在赣南地区的183次ML≥1.0地震资料,采用尾波归一化方法对寻乌、安远和龙南3个台站的P波、S波品质因子QP和QS值进行了计算,分析了该地区地壳介质的衰减特征。研究结果表明,寻乌(流逝时间35 s)和龙南台(流逝时间40 s)QP值较可靠且值相近,拟合平均值为QP(f)=(20.76+1.67)*f(1.03±0.04),寻乌台和安远台(流逝时间均为21 s)拟合平均值为QS(f)=(18.86+1.46)*f(0.84±0.04)。Q值总体呈现低Q0和高η值的特点,这是构造运动活跃地区的特征体现,也和该地区中强地震比较活跃的特征相符。     

Attenuation of P,S, and coda waves in Koyna region,India

The attenuation properties of the crust in the Koyna region of the Indian shield have been investigated using 164 seismograms from 37 local earthquakes that occurred in the region. The extended coda normalization method has been used to estimate the quality factors for P waves and S waves , and the single back-scattering model has been used to determine the quality factor for coda waves (Q _c). The earthquakes used in the present study have the focal depth in the range of 1–9 km, and the epicentral distance vary from 11 to 55 km. The values of and Q _c show a dependence on frequency in the Koyna region. The average frequency dependent relationships (Q = Q ₀ f ⁿ) estimated for the region are , and . The ratio is found to be greater than one for the frequency range considered here (1.5–18 Hz). This ratio, along with the frequency dependence of quality factors, indicates that scattering is an important factor contributing to the attenuation of body waves in the region. A comparison of Q _c and in the present study shows that for frequencies below 4 Hz and for the frequencies greater than 4 Hz. This may be due to the multiple scattering effect of the medium. The outcome of this study is expected to be useful for the estimation of source parameters and near-source simulation of earthquake ground motion, which in turn are required in the seismic hazard assessment of a region.     

Frequency-dependent Attenuation of High-frequency P and S Waves in the Upper Crust in Western Nagano, Japan

—Borehole seismograms from local earthquakes in the aftershock region of the 1984 western Nagano Prefecture, Japan earthquake were analyzed to measure the frequency-dependent characteristics of P- and S-wave attenuation in the upper crust. The records from a three-component velocity seismometer at the depth of 145m exhibit high S/N-ratio in a wide frequency range up to 100 Hz. Extended coda normalization methods were applied to bandpass-filtered seismograms of frequencies from 25 to 102 Hz. For the attenuation of high-frequency P and S waves, our measurements show Q _P ^-1? 0.052?^-0.66 and Q _S ^-1? 0.0034?^-0.12 respectively. The frequency dependence of the quality factor of S waves is very weak as compared with that of P waves. The ratio of Q _P ^-1/Q _S ^-1 is larger than unity in the entire analyzed frequency range.     

Attenuation of P- and S-waves in the Chamoli Region, Himalaya, India

The attenuation properties of the crust in the Chamoli region of Himalaya have been examined by estimating the frequency-dependent relationships of quality factors for P waves (Q_α) and for S waves (Q_β) in the frequency range 1.5–24 Hz. The extended coda normalization method has been applied on the waveforms of 25 aftershocks of the 1999 Chamoli earthquake (M 6.4) recorded at five stations. The average value of Q_α is found to be varied from 68 at 1.5 Hz to 588 at 24 Hz while it varies from 126 at 1.5 Hz to 868 at 24 Hz for Q_β. The estimated frequency-dependent relations for quality factors are Q_α = (44 ± 1)f^(0.82±.04) and Q_β = (87 ± 3)f^(0.71±.03). The rate of increase of Q(f) for P and S waves in the Chamoli region is comparable with the other regions of the world. The ratio Q_β/Q_α is greater than one in the region which along with the frequency dependence of quality factors indicates that scattering is an important factor contributing to the attenuation of body waves in the region. A comparison of attenuation relation for S wave estimated here (Q_β = 87f^0.71) with that of coda waves (Q_c = 30f^1.21) obtained by Mandal et al. (2001) for the same region shows that Q_c > Q_β for higher frequencies (>8 Hz) in the region. This indicates a possible high frequency coda enrichment which suggests that the scattering attenuation significantly influences the attenuation of S waves at frequencies >8 Hz. This observation may be further investigated using multiple scattering models. The attenuation relations for quality factors obtained here may be used for the estimation of source parameters and near-source simulation of earthquake ground motion of the earthquakes, which in turn are required for the assessment of seismic hazard in the region.     

Probabilistic Assessment of Earthquake Hazard in Gujarat and Adjoining Region of India

The Gujarat and adjoining region falls under all four seismic zones V, IV, III and II of the seismic zoning map of India, and is one of the most seismically prone intracontinental regions of the world. It has experienced two large earthquakes of magnitude M _w 7.8 and 7.7 in 1819 and 2001, respectively and several moderate earthquakes during the past two centuries. In the present study, the probability of occurrence of earthquakes of M ≥ 5.0 has been estimated during a specified time interval for different elapsed times on the basis of observed time intervals between earthquakes using three stochastic models namely, Weibull, Gamma and Lognormal. A complete earthquake catalogue has been used covering the time interval of 1819 to 2006. The whole region has been divided into three major seismic regions (Saurashtra, Mainland Gujarat and Kachchh) on the basis of seismotectonics and geomorphology of the region. The earthquake hazard parameters have been estimated using the method of maximum likelihood. The logarithmic of likelihood function (ln L) is estimated and used to test the suitability of models in three different regions. It was found that the Weibull model fits well with the actual data in Saurashtra and Kachchh regions, whereas Lognormal model fits well in Mainland Gujarat. The mean intervals of occurrence of earthquakes are estimated as 40.455, 20.249 and 13.338 years in the Saurashtra, Mainland Gujarat and Kachchh region, respectively. The estimated cumulative probability (probability that the next earthquake will occur at a time later than some specific time from the last earthquake) for the earthquakes of M ≥ 5.0 reaches 0.9 after about 64 years from the last earthquake (1993) in Saurashtra, about 49 years from the last earthquake (1969) in Mainland Gujarat and about 29 years from the last earthquake (2006) in the Kachchh region. The conditional probability (probability that the next earthquake will occur during some specific time interval after a certain elapsed time from last earthquake) is also estimated and it reaches about 0.8 to 0.9 during the time interval of about 57 to 66 years from the last earthquake (1993) in Saurashtra region, 31 to 51 years from the last earthquake (1969) in Mainland Gujarat and about 21 to 28 years from the last earthquake (2006) in Kachchh region.     

Attenuation of coda waves in the Northeastern Region of India

Coda wave attenuation quality factor Qc is estimated in the northeastern region of India using 45 local earthquakes recorded by regional seismic network. The quality factor Qc was estimated using the single backscattering model modified by Sato (J Phys Earth 25:27–41, 1977), in the frequency range 1–18 Hz. The attenuation and frequency dependence for different paths and the correlation of the results with geotectonics of the region are described in this paper. A total of 3,890 Qc measurements covering 187 varying paths are made for different lapse time window of 20, 30, 40, 50, 60, 70, 80, and 90 s in coda wave. The magnitudes of the analyzed events range from 1.2 to 3.9 and focal depths range between 7 and 38 km. The source–receiver distances of the selected events range between 16 and 270 km. For 30-s duration, the mean values of the estimated Qc vary from 50 ± 12 (at 1 Hz) to 2,078 ± 211(at 18 Hz) for the Arunachal Himalaya, 49 ± 14 (at 1 Hz) to 2,466 ± 197 (at 18 Hz) for the Indo-Burman, and 45 ± 13 (at 1 Hz) to 2,069 ± 198 (at 18 Hz) for Shillong group of earthquakes. It is observed that Qc increases with frequency portraying an average attenuation relation for the region. Moreover, the pattern of Qc^− 1 with frequency is analogous to the estimates obtained in other tectonic areas in the world, except with the observation that the Qc^− 1 is much higher at 1 Hz for the northeastern region. The Qc^− 1 is about 10^− 1.8 at 1 Hz and decreases to about 10^− 3.6 at 18 Hz indicating clear frequency dependence. Pertaining to the spatial distribution of Qc values, Mikir Hills and western part of Shillong Plateau are characterized by lower attenuation.     

QLg tomography in Gujarat,Western India

We propose a novel Lg attenuation tomography model (Q_Lg tomography) for the state of Gujarat, Western India, using earthquake data recorded by the Gujarat Seismic Network, operated by the Institute of Seismological Research in Gandhinagar. The waveform dataset consist of 400 3-component recordings, produced by 60 earthquakes with magnitude (M_L) spanning from 3.6 to 5.1, recorded at 60 seismic stations having epicentral distances spanning between 200 and 500 km. Spectral amplitude decays for Lg wave displacement were obtained by generalized inversion at 17 frequencies spanning between 0.9 and 9 Hz. Lg wave propagation efficiency was measured by Lg/Pn spectral ratio categorizing as efficient ratio ≥6 for 86%, intermediate ratio of 3–6 for 10% and inefficient ratio <3 for 4% paths of total 400 ray paths. The earthquake size and quality of waveform recorded at dense network found sufficient to resolve lateral variation of Q_Lg in Gujarat.Average power-law attenuation relationship obtained for Gujarat as Q_Lg(f) = 234f^0.64, which corresponds to high attenuation in comparison to peninsular India shield region and other several regions around the world. Q_Lg tomography resolves the highly attenuating crust of extremely fractured Saurashtra region and tectonically active Kachchh region. The Gujarat average attenuation is also lying in between them. The low attenuation in Cambay and Narmada rift basins and extremely low attenuation in patch of Surendranagar area is identified. This study is the first attempt and can be utilized as pivotal criteria for scenario hazard assessment, as maximum hazard has been reported in highly attenuating tectonically active Kachchh region and in low attenuating Cambay, Narmada and Surendranagar regions. The site and source terms are also obtained along with the Q_Lg inversion. The estimated site responses are comparable with observed local geological condition and agree with the previously reported site amplifications at the same sites. The source terms are comparable with local magnitude estimated from Network. The M_w (Lg) is nearly equivalent to M_L (GSN) and the slight differences are noted for larger magnitude events.     

龙滩水库地区P波、S波和尾波衰减

利用发生在龙滩水库地区的144次ML≥1.0地震资料,研究了该地区地壳介质的衰减特征,并采用尾波归一化方法分别估算了P波、S波的品质因子QP和QS,用单次反散射模型计算了尾波品质因子QC．结果表明,该地区上述3个值与频率的依赖关系分别为QP（f）=（7.43±0.78）f(1.05±0.06),QS（f）=（13.08±3.20）f(1.76±0.15),QC（f）=（50.52±2.71）f(1.13±0.03),总体呈现低QO和高η值的特点．在所研究的频带范围1—16Hz内,QS／QP比值均大于1,表明该地区体波衰减中散射可能起着较大作用．QC与QS的比较结果表明,在小于8Hz时,QS＜QC;而大于8Hz时,则QS＞QC．这可能是由于介质的多次散射作用的结果．     

Stress Rotation in the Kachchh Rift Zone, Gujarat, India

Double difference relocations of the 1402 Kachchh events (2001–2006) clearly delineate two fault zones viz. south-dipping North Wagad fault (NWF) and almost vertical Gedi fault (GF). The relocated focal depths delineate a marked variation of 4 and 7 km in the brittle-ductile transition depths beneath GF and NWF, respectively. The focal mechanism solutions of 464 aftershocks (using 8–12 first motions) show that the focal mechanisms ranged between pure reverse and pure strike-slip except for a few pure dip-slip solutions. The stress inversions performed for five rectangular zones across the Kachchh rift reveal both clockwise and anticlockwise rotation (7–32°) in the σ₁ orientation within the rupture zone, favoring a heterogeneous stress regime with an average N-S fault normal compression. This rotation may be attributed to the presence of crustal mafic intrusives (5–35 km depth) in the rupture zone of the 2001 Bhuj main shock. Results suggest a relatively homogeneous stress regime in the GF zone favoring strike-slip motion, with a fault normal N-S compression.     

Over-exploitation of an alluvial aquifer in Gujarat,India

Abstract

Over-exploitation of an alluvial aquifer in India has occurred since tubewells were drilled into deeper aquifers. Field information concerning-the fall in both the water table and the piezometric heads in the deeper aquifers has been examined and a numerical model has been developed to represent a heavily exploited area. Predictions using the model suggest that the annual declines in the water table and the piezometric heads are accelerating as the water table moves into a zone having a lower specific yield. Careful monitoring of the movement of the water table is proposed.     

Magmatic differentiation of the volcanics at Pavagarh,Gujarat, India

Magmatic differentiation in Pavagarh adopted the courses corresponding to the olivine-basalt trend including olivine-basalt, mugearite, latite and pitch-stone and the normal basalt trend comprising tholeiitic basalts, andesite, dellenite, and rhyolite. The dellenite-rhyolite assemblage is not consanguineous. Their average composition corresponds with the petrochemical field of the basalt-andesite-dacite-rhyolite association or the calc-alkaline series of Nockolds. It lies in between plateau and tholeiitic basalts of different provinces and shows alkali enrichment in the initial phases of differentiation. The main factors governing the nature and course of differentiation are the changing concentration of the volatiles during different phases of magmatic activity, the different composition of the liquid during the successive phases of its eruption and the partial assimilation of the siliceous rocks forming the platform over which the present volcanics were laid down. Other factors are crystal settling (ankaramites), reaction of the early-formed crystals with the remaining magma (olivine-basalt) and mixing of the magmatic liquids belonging to different phases (mugearites). The chronological history of the magmatic activity in relation to the geological setting of the area shows that the rocks belong to a single parent magma closely akin to olivine-basalt or plateu basalt magma type which through differentiation into two lines of descent, has given rise to the present volcanic assemblage.     

甘肃地区S波非弹性衰减Q值研究

文中利用新建"十五"数字化地震观测网络项目记录到的甘肃地区地震波形资料,研究了甘肃地区S波的非弹性衰减Q值,得到甘肃东南地区Q值与频率的关系为Q(f)=374.0f0.3084;祁连山中东段地区Q值与频率的关系为Q(f)=693.6f0.2865。从分区Q值研究结果看,各分区之间存在较大差异,显示出了Q值空间分布的非均匀性特征,这种区域差异性与甘肃地区构造活动、地震活动、地壳速度结构及热活动状态的区域差异性有较好的一致性,活动强烈、高热活动以及低速区为相对低Q值区,地震波衰减较快。     

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.     

Body-wave Attenuation in the Region of Garda, Italy

We analyzed the spectral amplitude decay with hypocentral distance of P and S waves generated by 76 small magnitude earthquakes (M_L 0.9–3.8) located in the Garda region, Central-Eastern Alps, Italy. These events were recorded by 18 stations with velocity sensors, in a distance range between 8 and 120 km. We calculated nonparametric attenuation functions (NAF) and estimated the quality factor Q of both body waves at 17 different frequencies between 2 and 25 Hz. Assuming a homogeneous model we found that the Q frequency dependence of P and S can be approximated with the functions Q _P = 65 f ^0.9 and Q _S = 160 f ^0.6 , respectively. At 2 Hz the Q _S /Q _P ratio reaches the highest value of 2.8. At higher frequencies Q _S /Q _P varies between 0.7 and 1.7, suggesting that for this frequency band scattering may be an important attenuation mechanism in the region of Garda. To explore the variation of Q in depth, we estimated Q at short (r ≤ 30 km) and intermediate (35–90 km) distance paths. We found that in the shallow crust P waves attenuate more than S (1.3 < Q _S /Q _P < 2.5). Moreover, P waves traveling along paths in the lower crust (depths approximately greater than 30 km) attenuate more than S waves. To quantify the observed variability of Q in depth we considered a three-layer model and inverted the NAF to estimate Q in each layer. We found that in the crust Q increases with depth. However, in the upper mantle (~40–50 km depth) Q decreases and in particular the high frequency Q _S (f > 9 Hz) has values similar to those estimated for the shallow layer of the crust.     

石油地震纵横波速度层析成像

本文以分层均匀“多层介质模型”为依据,正演模拟多炮人工地面地震记录。其观测方式为:三炮位置分别在1500m,2500m,3500m处;道检距25m,140道(或道检距5m,700道),偏移距1000m。为提高地震反射层析成像的精度,运用“走时最小纵横波速度层析成像”与“波形基本解纵横波速度层析成像”相结合的方法,反演各层界面形状及各层纵波、横波速度,这是本文的特点。反演结果图与“多层介质模型”基本一致,相对误差小于5%,达到实际部门的要求。我们在数学上发展了“层析成像”的内容,如探索出求二维场强公式用一维基本解表示,因此大大提高了数值计算速度,具有一定的实用与理论价值。用“ProMAX软件”分析实际地震数据,获得地下11层垂直剖面图与纵横波速度层析成像结果。“项目程序系统”可供相关部门实际参考使用。     

The Effects of Attenuation and Site on the Spectra of Microearthquakes in the Shillong Region of Northeast India

Microearthquake spectra from the Shillong region are analyzed to observe the effect of attenuation and site on these spectra. The spectral ratio method is utilized to estimate the Q values for both P- and S-waves in the subsurface layer, wherein the ratio of spectral amplitudes at lower and higher frequencies are taken into consideration for three stations at varying epicentral distances. Average estimates of Q _P and Q _S are 178 and 195. The ratio of Q _S to Q _P is estimated to be greater than 1 in major parts of the Shillong area, which can be related to the dry crust prevailing in the Shillong region. Typically, the variation in corner frequencies for these spectra is inferred to be characteristic of the site. Simultaneously, observations from spectral content of local earthquakes recorded at two different stations with respect to the reference site yield greater amplification of incoming seismic signals in the frequency range of 2–5 Hz, which is found to be well supported by the existing local lithology pertinent to that region.     

Seismic Wave Attenuation in the Greater Cairo Region, Egypt

In the present study, a digital waveform dataset of 216 local earthquakes recorded by the Egyptian National Seismic Network (ENSN) was used to estimate the attenuation of seismic wave energy in the greater Cairo region. The quality factor and the frequency dependence for Coda waves and S-waves were estimated and clarified. The Coda waves (Q _c) and S-waves (Q _d) quality factor were estimated by applying the single scattering model and Coda Normalization method, respectively, to bandpass-filtered seismograms of frequency bands centering at 1.5, 3, 6, 12, 18 and 24?Hz. Lapse time dependence was also studied for the area, with the Coda waves analyzed through four lapse time windows (10, 20, 30 and 40?s). The average quality factor as function of frequency is found to be Q _c?=?35?±?9f ^0.9±0.02 and Q _d?=?10?±?2f ^0.9±0.02 for Coda and S-waves, respectively. This behavior is usually correlated with the degree of tectonic complexity and the presence of heterogeneities at several scales. The variation of Q _c with frequency and lapse time shows that the lithosphere becomes more homogeneous with depth. In fact, by using the Coda Normalization method we obtained low Q _d values as expected for a heterogeneous and active zone. The intrinsic quality factor (Q _i ^?1 ) was separated from the scattering quality factor (Q _s ^?1 ) by applying the Multiple Lapse Time Domain Window Analysis (MLTWA) method under the assumption of multiple isotropic scattering with uniform distribution of scatters. The obtained results suggest that the contribution of the intrinsic attenuation (Q _i ^?1 ) prevails on the scattering attenuation (Q _s ^?1 ) at frequencies higher than 3?Hz.