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.     

Variation of intrinsic and scattering attenuation of seismic waves with depth in the Bam region,East-Central Iran

The attenuation properties of the lithosphere in the Bam region, East-Central Iran, have been investigated. For this purpose, 42 local earthquakes having focal depths less than 25 km have been used. The quality factor of coda waves (Q_c) has been estimated using the single back-scattering model. The quality factors Q_p, Q_d (P and direct S-waves) have been estimated using the extended-coda normalization method. Q_i and Q_s (the intrinsic and scattering attenuation parameters) have been estimated for the region. The values of Q_p, Q_d, Q_c, Q_i and Q_s show a dependence on frequency in the range of 1.5–24 Hz for the Bam region. The average frequency-dependent relationships estimated for the region are Q_p=(36±6)f^(1.03±0.06), Q_d=(59±8)f^(1.00±0.03), Q_c=(79±5)f^(1.01±0.04), Q_s=(131±4)f^(1.01±0.04) and Q_i=(104±6)f^(1.01±0.05). A comparison between Q_i and Q_s shows that intrinsic absorption is predominant over scattering.The variation of Q has also been estimated at different lapse times to observe heterogeneities variation with depth. The variation of Q with frequency and lapse time shows that the lithosphere becomes more homogeneous with depth.The estimated Q_o values at different stations suggest a low value of Q indicating a heterogeneous and attenuative crust beneath the entire region.     

Estimation of Q p and Q s of Kinnaur Himalaya

The attenuation characteristics of the Kinnaur area of the North West Himalayas were studied using local earthquakes that occurred during 2008–2009. Most of the analyzed events are from the vicinity of the Panjal Thrust (PT) and South Tibetan Detachment Thrust, which are well-defined tectonic discontinuities in the Himalayas. The frequency-dependent attenuation of P and S waves was estimated using the extended coda normalization method. Data from 64 local earthquakes recorded at 10 broadband stations were used. The coda normalization of the spectral amplitudes of P and S waves was done at central frequencies of 1.5, 3, 6, 9, and 12 Hz. Q _p increases from about 58 at 1.5 Hz to 706 at 12 Hz, and Q _s increases from 105 at 1.5 Hz to 1,207 at 12 Hz. The results show that the quality factors for both P and S waves (Q _p and Q _s) increase as a function of frequency according to the relation Q?=?Q _o f ⁿ , where Q _o is the corresponding Q value at 1 Hz frequency and “n” is the frequency relation parameter. We obtained Q _p?=?(47?±?2)f ^(1.04±0.04) and Q _s?=?(86?±?4)f ^(0.96±0.03) by fitting power law dependency model for the estimated values of the entire study region. The Q ₀ and n values show that the region is seismically very active and the crust is highly heterogeneous. There was no systematic variation of values of Q _p and Q _s at different frequencies from one tectonic unit to another. As a consequence, average values of these parameters were obtained for each frequency for the entire region, and these were used for interpretation and for comparison with worldwide data. Q _p values lie within the range of values observed for some tectonically active regions of the world, whereas Q _s values were the lowest among the values compared for different parts of the world. Q _s/Q _p values were >1 for the entire range of frequencies studied. All these factors indicate that the crust is highly heterogeneous in the study region. The high Q _s/Q _p values also indicate that the region is partially saturated with fluids.     

Coda Q Estimates in the Andaman Islands Using Local Earthquakes

The attenuation property of Andaman Island has been investigated analyzing coda waves from 57 local earthquakes in the magnitude range of 2.0–4.9, using the single backscattering model. These earthquakes waveforms, recorded on five broadband seismographs sited over the island from north to south during Nov. 2003 to March 2004, have been used to calculate the frequency dependent Coda Q (Q _c ) applying the time domain coda-decay method. The Coda Q, computed at central frequencies from (0.5–12) Hz and five-lapse time windows from 40 to 80 s, progressively increases from 105 f ^0.88 in the north Andaman to 135 f ^0.79 in the south Andaman with an average of 119 f ^0.80. The average Q _c values vary from 75 ± 42 at 0.5 Hz to 697 ± 54 at 12 Hz central frequency for 40 s lapse time window, while for 80 s lapse time window its variation is from 117 ± 38 at 0.5 Hz to 1256 ± 115 at 12 Hz. The Q _c estimated at different lapse times manifests a significant variation from 122f ^0.75 to 174f ^0.73, corresponding to lapse time window lengths of 40 and 80 s, respectively. The variation of Q _c with frequency, lapse time and also with the location of seismograph reflects the marked structural and compositional inhomogeneity with depth along the Andaman Islands. These observations are well correlated with the seismicity pattern and distinct high angle subduction beneath the island.     

Frequency-Dependent Attenuation of Coda Waves in the Crust in Southwest Anatolia (Turkey)

The attenuation of coda waves in the earth’s crust in southwest (SW) Anatolia is estimated by using the coda wave method, which is based on the decrease of coda wave amplitude in time and distance. A total of 159 earthquakes were recorded between 1997 and 2010 by 11 stations belonging to the KOERI array. The coda quality factor Q _c is determined from the properties of scattered coda waves in a heterogeneous medium. Firstly, the quality factor Q ₀ (the value of Q _c at 1 Hz.) and its frequency dependency η are determined from this method depending on the attenuation properties of scattered coda waves for frequencies of 1.5, 3.0, 6.0, 8.0, 12 and 20 Hz. Secondly, the attenuation coefficients (δ) are estimated. The shape of the curve is controlled by the scattering and attenuation in the crustal volume sampled by the coda waves. The average Q _c values vary from 110 ± 15 to 1,436 ± 202 for the frequencies above. The Q ₀ and η values vary from 63 ± 7 to 95 ± 10 and from 0.87 ± 0.03 to 1.04 ± 0.09, respectively, for SW Anatolia. In this region, the average coda Q–f relation is described by Q _c = (78 ± 9)f ^0.98±0.07 and δ = 0.012 km^?1. The low Q ₀ and high η are consistent with a region characterized by high tectonic activity. The Q _c values were correlated with the tectonic pattern in SW Anatolia.     

On the Causes of Frequency-Dependent Apparent Seismological Q

Variability of the Earth’s structure makes a first-order impact on attenuation measurements which often does not receive adequate attention. Geometrical spreading (GS) can be used as a simple measure of the effects of such structure. The traditional simplified GS compensation is insufficiently accurate for attenuation measurements, and the residual GS appears as biases in both Q ₀ and η parameters in the frequency-dependent attenuation law Q(f) = Q ₀ f ^η . A new interpretation approach bypassing Q(f) and using the attenuation coefficient χ(f) = γ + πf/Q _e(f) resolves this problem by directly measuring the residual GS, denoted γ, and effective attenuation, Q _e. The approach is illustrated by re-interpreting several published datasets, including nuclear-explosion and local-earthquake codas, Pn, and synthetic 50–300-s surface waves. Some of these examples were key to establishing the Q(f) concept. In all examples considered, χ(f) shows a linear dependence on the frequency, γ ≠ 0, and Q _e can be considered frequency-independent. Short-period crustal body waves are characterized by positive γ _SP values of (0.6–2.0) × 10^?2 s^?1 interpreted as related to the downward upper-crustal reflectivity. Long-period surface waves show negative γ _LP ≈ ?1.9 × 10^?5 s^?1, which could be caused by insufficient modeling accuracy at long periods. The above γ values also provide a simple explanation for the absorption band observed within the Earth. The band is interpreted as apparent and formed by levels of Q _e ≈ 1,100 within the crust decreasing to Q _e ≈ 120 within the uppermost mantle, with frequencies of its flanks corresponding to γ _LP and γ _SP. Therefore, the observed absorption band could be purely geometrical in nature, and relaxation or scattering models may not be necessary for explaining the observed apparent Q(f). Linearity of the attenuation coefficient suggests that at all periods, the attenuation of both Rayleigh and Love waves should be principally accumulated at the sub-crustal depths (~38–100 km).     

Body wave attenuation characteristics in the crust of Alborz region and North Central Iran

Attenuation of P and S waves has been investigated in Alborz and north central part of Iran using the data recorded by two permanent and one temporary networks during October 20, 2009, to December 22, 2010. The dataset consists of 14,000 waveforms from 380 local earthquakes (2 < M _L < 5.6). The extended coda normalization method (CNM) was used to estimate quality factor of P (Q _P) and S waves (Q _S) at seven frequency bands (0.375, 0.75, 1.5, 3, 6, 12, 24 Hz). The Q _P and Q _S values have been estimated at lapse times from 40 to 100 s. It has been observed that the estimated values of Q _P and Q _S are time independent; therefore, the mean values of Q _P and Q _S at different lapse times have been considered. The frequency dependence of quality factor was determined by using a power-law relationship. The frequency-dependent relationship for Q _P was estimated in the form of (62 ± 7)f ^{(1.03 ± 0.07)} and (48 ± 5)f ^{(0.95 ± 0.07)} in Alborz region and North Central Iran, respectively. These relations for Q _S for Alborz region and North Central Iran have estimated as (83 ± 8)f ^{(0.99 ± 0.07)} and (68 ± 5)f ^{(0.96 ± 0.05)}, respectively. The observed low Q values could be the results of thermoelastic effects and/or existing fracture. The estimated frequency-dependent relationships are comparable with tectonically active regions.     

Near-source attenuation of high-frequency body waves beneath the New Madrid Seismic Zone

Attenuation characteristics in the New Madrid Seismic Zone (NMSZ) are estimated from 157 local seismograph recordings out of 46 earthquakes of 2.6?≤?M?≤?4.1 with hypocentral distances up to 60 km and focal depths down to 25 km. Digital waveform seismograms were obtained from local earthquakes in the NMSZ recorded by the Center for Earthquake Research and Information (CERI) at the University of Memphis. Using the coda normalization method, we tried to determine Q values and geometrical spreading exponents at 13 center frequencies. The scatter of the data and trade-off between the geometrical spreading and the quality factor did not allow us to simultaneously derive both these parameters from inversion. Assuming 1/R ^1.0 as the geometrical spreading function in the NMSZ, the Q _P and Q _S estimates increase with increasing frequency from 354 and 426 at 4 Hz to 729 and 1091 at 24 Hz, respectively. Fitting a power law equation to the Q estimates, we found the attenuation models for the P waves and S waves in the frequency range of 4 to 24 Hz as Q _P?=?(115.80?±?1.36) f ^{(0.495?±?0.129)} and Q _S?=?(161.34?±?1.73) f ^{(0.613?±?0.067)}, respectively. We did not consider Q estimates from the coda normalization method for frequencies less than 4 Hz in the regression analysis since the decay of coda amplitude was not observed at most bandpass filtered seismograms for these frequencies. Q _S/Q _P?>?1, for 4?≤?f?≤?24 Hz as well as strong intrinsic attenuation, suggest that the crust beneath the NMSZ is partially fluid-saturated. Further, high scattering attenuation indicates the presence of a high level of small-scale heterogeneities inside the crust in this region.     

Coda wave attenuation in the Parecis Basin,Amazon Craton,Brazil: sensitivity to basement depth

Small local earthquakes from two aftershock sequences in Porto dos Gaúchos, Amazon craton—Brazil, were used to estimate the coda wave attenuation in the frequency band of 1 to 24 Hz. The time-domain coda-decay method of a single backscattering model is employed to estimate frequency dependence of the quality factor (Q _c) of coda waves modeled using Q_c = Q₀ f^hQ_{\rm c} =Q_{\rm 0} f^\eta , where Q ₀ is the coda quality factor at frequency of 1 Hz and η is the frequency parameter. We also used the independent frequency model approach (Morozov, Geophys J Int, 175:239–252, 2008), based in the temporal attenuation coefficient, χ(f) instead of Q(f), given by the equation c(f)=g+\fracpfQ_e \chi (f)\!=\!\gamma \!+\!\frac{\pi f}{Q_{\rm e} }, for the calculation of the geometrical attenuation (γ) and effective attenuation (Q_e^-1 )(Q_{\rm e}^{-1} ). Q _c values have been computed at central frequencies (and band) of 1.5 (1–2), 3.0 (2–4), 6.0 (4–8), 9.0 (6–12), 12 (8–16), and 18 (12–24) Hz for five different datasets selected according to the geotectonic environment as well as the ability to sample shallow or deeper structures, particularly the sediments of the Parecis basin and the crystalline basement of the Amazon craton. For the Parecis basin Q_c = (98±12)f^(1.14±0.08)Q_{\rm c} =(98\pm 12)f^{(1.14\pm 0.08)}, for the surrounding shield Q_c = (167±46)f^(1.03±0.04)Q_{\rm c} =(167\pm 46)f^{(1.03\pm 0.04)}, and for the whole region of Porto dos Gaúchos Q_c = (99±19)f^(1.17±0.02)Q_{\rm c} =(99\pm 19)f^{(1.17\pm 0.02)}. Using the independent frequency model, we found: for the cratonic zone, γ = 0.014 s^− 1, Q_e^-1 = 0.0001Q_{\rm e}^{-1} =0.0001, ν ≈ 1.12; for the basin zone with sediments of ~500 m, γ = 0.031 s^− 1, Q_e^-1 = 0.0003Q_{\rm e}^{-1} =0.0003, ν ≈ 1.27; and for the Parecis basin with sediments of ~1,000 m, γ = 0.047 s^− 1, Q_e^-1 = 0.0005Q_{\rm e}^{-1} =0.0005, ν ≈ 1.42. Analysis of the attenuation factor (Q _c) for different values of the geometrical spreading parameter (ν) indicated that an increase of ν generally causes an increase in Q _c, both in the basin as well as in the craton. But the differences in the attenuation between different geological environments are maintained for different models of geometrical spreading. It was shown that the energy of coda waves is attenuated more strongly in the sediments, Q_c = (78±23)f^(1.17±0.14)Q_{\rm c} =(78\pm 23)f^{(1.17\pm 0.14)} (in the deepest part of the basin), than in the basement, Q_c = (167±46)f^(1.03±0.04)Q_{\rm c} =(167\pm 46)f^{(1.03\pm 0.04)} (in the craton). Thus, the coda wave analysis can contribute to studies of geological structures in the upper crust, as the average coda quality factor is dependent on the thickness of sedimentary layer.     

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.     

Region-specific <Emphasis Type="Italic">S</Emphasis>-wave attenuation for earthquakes in northwestern Iran

In this study, continuous wavelet transform is applied to estimate the frequency-dependent quality factor of shear waves, Q _S , in northwestern Iran. The dataset used in this study includes velocigrams of more than 50 events with magnitudes between 4.0 and 6.5, which have occurred in the study area. The CWT-based method shows a high-resolution technique for the estimation of S-wave frequency-dependent attenuation. The quality factor values are determined in the form of a power law as Q _S (f)?=?(147?±?16)f ^0.71?±?0.02 and (126?±?12)f ^0.73?±?0.02 for vertical and horizontal components, respectively, where f is between 0.9 and 12 Hz. Furthermore, in order to verify the reliability of the suggested Q _S estimator method, an additional test is performed by using accelerograms of Ahar-Varzaghan dual earthquakes on August 11, 2012, of moment magnitudes 6.4 and 6.3 and their aftershocks. Results indicate that the estimated Q _S values from CWT-based method are not very sensitive to the numbers and types of waveforms used (velocity or acceleration).     

Coda Q Estimates in the Koyna Region, India

—The coda Q, Q _c ?, have been estimated for the Koyna region of India. The coda waves of 76 seismograms from thirteen local earthquakes, recorded digitally in the region during July–August, 1996, have been analyzed for this purpose at nine central frequencies viz., 1.5, 2.0, 3.0, 4.0, 6.0, 8.0, 12.0, 16.0 and 24.0 Hz using a single backscattering model. All events with magnitude less than 3 fall in the epicentral distances less than 60 km and have focal depths which range from 0.86 to 9.43 km. For the 30 sec coda window length the estimated Q _c values vary from 81 to 261 at 1.5 Hz and 2088 to 3234 at 24 Hz, whereas the mean values of Q _c with the standard error vary from 148 ± 13.5 at 1.5 Hz to 2703 ± 38.8 at 24 Hz. Both the estimated Q _c values and their mean values exhibit the clear dependence on frequency in the region and a frequency dependence average attenuation relationship, Q _c = 96f ^1.09, has been obtained for the region, covering an approximate area of 11500 km² with the surfacial extent of about 120 km and depth of 60 km.¶Lapse time dependence of Q _c has also been studied for the region, with the coda waves analyzed at five lapse time windows from 20 to 60 sec duration with the difference of 10 sec. The frequency dependence average Q _c relationships obtained at these window lengths Q _c = 66f ^1.16 (20 sec), Q _c = 96f ^1.09 (30 sec), Q _c =131f ^1.04 (40 sec), Q _c = 148f ^1.04 (50 sec), Q _c = 182f ^1.02 (60 sec) show that the frequency dependence (exponentn) remains mostly stationary at all the lapse time window lengths, while the change in Q ₀ value is significant. Lapse time dependence of Q _c in the region is also interpreted as the function of depth.     

<Emphasis Type="Italic">Q</Emphasis><Subscript>c</Subscript> and <Emphasis Type="Italic">Q</Emphasis><Subscript>S</Subscript> wave attenuation of South African earthquakes

Quality factor Q, which describes the attenuation of seismic waves with distance, was determined for South Africa using data recorded by the South African National Seismograph Network. Because of an objective paucity of seismicity in South Africa and modernisation of the seismograph network only in 2007, I carried out a coda wave decay analysis on only 13 tectonic earthquakes and 7 mine-related events for the magnitude range 3.6?≤?M _L ?≤?4.4. Up to five seismograph stations were utilised to determine Q _c for frequencies at 2, 4, 8 and 16 Hz resulting in 84 individual measurements. The constants Q ₀ and α were determined for the attenuation relation Q _c(f)?=?Q ₀ f ^α . The result was Q ₀?=?396?±?29 and α?=?0.72?±?0.04 for a lapse time of 1.9*(t _s???t ₀) (time from origin time t ₀ to the start of coda analysis window is 1.9 times the S-travel time, t _s) and a coda window length of 80 s. This lapse time and coda window length were found to fit the most individual frequencies for a signal-to-noise ratio of at least 3 and a minimum absolute correlation coefficient for the envelope of 0.5. For a positive correlation coefficient, the envelope amplitude increases with time and Q _c was not calculated. The derived Q _c was verified using the spectral ratio method on a smaller data set consisting of nine earthquakes and one mine-related event recorded by up to four seismograph stations. Since the spectral ratio method requires absolute amplitudes in its calculations, site response tests were performed to select four appropriate stations without soil amplification and/or signal distortion. The result obtained for Q _S was Q ₀?=?391?±?130 and α?=?0.60?±?0.16, which agrees well with the coda Q _c result.     

The coda attenuation of the Yao’an area in Yunnan Province

An earthquake withM=6.5 happened on January 15, 2000 in Yao’an of Yunnan Province. After the earthquake, a temporary digital network with 6 detectors around the epicenter area was set up. 402 aftershocks were located more precisely. According to coda short recording observed, the coda averaging quality factor has been acquired via Sato’s single scattering model analyses,Q _c(f)=49f ^0.95,f=1.5～20.0 Hz, which has the attenuation characteristics of high structural active 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.     

Attenuation of High Frequency P and S Waves in the Gujarat Region,India

The local earthquake waveforms recorded on broadband seismograph network of Institute of Seismological Research in Gujarat, India have been analyzed to understand the attenuation of high frequency (2–25 Hz) P and S waves in the region. The frequency dependent relationships for quality factors for P (Q _P) and S (Q _S) waves have been obtained using the spectral ratio method for three regions namely, Kachchh, Saurashtra and Mainland Gujarat. The earthquakes recorded at nine stations of Kachchh, five stations of Saurashtra and one station in mainland Gujarat have been used for this analysis. The estimated relations for average Q _P and Q _S are: Q _P = (105 ± 2) f ^{0.82 ± 0.01}, Q _S = (74 ± 2) f ^{1.06 ± 0.01} for Kachchh region; Q _P = (148 ± 2) f ^{0.92 ± 0.01}, Q _S = (149 ± 14) f ^{1.43 ± 0.05} for Saurashtra region and Q _P = (163 ± 7) f ^{0.77 ± 0.03}, Q _S = (118 ± 34) f ^{0.65 ± 0.14} for mainland Gujarat region. The low Q (<200) and high exponent of f (>0.5) as obtained from present analysis indicate the predominant seismic activities in the region. The lowest Q values obtained for the Kachchh region implies that the area is relatively more attenuative and heterogeneous than other two regions. A comparison between Q _S estimated in this study and coda Q (Qc) previously reported by others for Kachchh region shows that Q _C > Q _S for the frequency range of interest showing the enrichment of coda waves and the importance of scattering attenuation to the attenuation of S waves in the Kachchh region infested with faults and fractures. The Q _S/Q _P ratio is found to be less than 1 for Kachchh and Mainland Gujarat regions and close to unity for Saurashtra region. This reflects the difference in the geological composition of rocks in the regions. The frequency dependent relations developed in this study could be used for the estimation of earthquake source parameters as well as for simulating the strong earthquake ground motions in the region.     

Estimation of site effects using strong motion data of BYTNet array in Turkey

Simultaneous estimation of effects of source, propagation path, and local site amplification was carried out using observed strong motion records in a frequency range from 0.8 to 20 Hz for the purpose of empirical evaluation of the local site effects in different geological conditions in the northwestern part of Turkey. The analyzed data are S-wave portions of 162 accelerograms from 39 shallow events observed at 14 sites of BYTNet array. A spectral separation method was applied to the observed S-wave spectra. The solutions for source spectra, inelasticity factor of propagation path for S-waves (Q _s-value), and factor of site amplification at each site were obtained in a least squares sense. In the analysis, we assumed that the factor of the site amplification at a reference site is the same as that of theoretical amplification of S-waves to the soil model whose bottom layer has an S-wave velocity around 2.15 km/s. The estimated Q _s-value of the propagation path is modeled as Q _s(f)?=?87.4f^0.78. The estimated site amplifications are characterized into three groups. The sites in the first group belong to rock site with no dominant peaks at a frequency range of 2 to 10 Hz. The second group of hard soil sites is characterized with moderately dominant peaks at a frequency of 5 Hz. The last group for soft soil sites has common peaks at a frequency of 4 Hz with larger amplitudes than those in the hard soil group. We, then, compare the amplifications with average S-wave velocity in top 30 m of the shallow S-wave profiles and proposed linear empirical formula between them at each frequency. We, furthermore, inverted the observed amplification factors into S-wave velocity and Q _s-value profiles of the deep soil over the basement.     

A study of Guptkashi,Uttarakhand earthquake of 6 February 2017 (<Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript>5.3) in the Himalayan arc and implications for ground motion estimation

The 2017 Guptkashi earthquake occurred in a segment of the Himalayan arc with high potential for a strong earthquake in the near future. In this context, a careful analysis of the earthquake is important as it may shed light on source and ground motion characteristics during future earthquakes. Using the earthquake recording on a single broadband strong-motion seismograph installed at the epicenter, we estimate the earthquake’s location (30.546° N, 79.063° E), depth (H?=?19 km), the seismic moment (M₀?=?1.12×10¹⁷ Nm, M _w 5.3), the focal mechanism (φ?=?280°, δ?=?14°, λ?=?84°), the source radius (a?=?1.3 km), and the static stress drop (Δσ _s ~22 MPa). The event occurred just above the Main Himalayan Thrust. S-wave spectra of the earthquake at hard sites in the arc are well approximated (assuming ω^?2 source model) by attenuation parameters Q(f)?=?500f^0.9, κ?=?0.04 s, and f_max?=?infinite, and a stress drop of Δσ?=?70 MPa. Observed and computed peak ground motions, using stochastic method along with parameters inferred from spectral analysis, agree well with each other. These attenuation parameters are also reasonable for the observed spectra and/or peak ground motion parameters in the arc at distances ≤?200 km during five other earthquakes in the region (4.6?≤?M _w ?≤?6.9). The estimated stress drop of the six events ranges from 20 to 120 MPa. Our analysis suggests that attenuation parameters given above may be used for ground motion estimation at hard sites in the Himalayan arc via the stochastic method.     

Attenuation coefficients of Rayleigh and <Emphasis Type="Italic">Lg</Emphasis> waves

Analysis of the frequency dependence of the attenuation coefficient leads to significant changes in interpretation of seismic attenuation data. Here, several published surface-wave attenuation studies are revisited from a uniform viewpoint of the temporal attenuation coefficient, denoted by χ. Theoretically, χ( f) is expected to be linear in frequency, with a generally non-zero intercept γ?=?χ(0) related to the variations of geometrical spreading, and slope dχ/df = π/Q _e caused by the effective attenuation of the medium. This phenomenological model allows a simple classification of χ( f) dependences as combinations of linear segments within several frequency bands. Such linear patterns are indeed observed for Rayleigh waves at 500–100-s and 100–10-s periods, and also for Lg from ~2 s to ~1.5 Hz. The Lg χ( f) branch overlaps with similar linear branches of body, Pn, and coda waves, which were described earlier and extend to ~100 Hz. For surface waves shorter than ~100 s, γ values recorded in areas of stable and active tectonics are separated by the levels of \(\gamma _{D} \approx 0.2 \times 10^{-3}\) s^???1 (for Rayleigh waves) and 8 ×10^???3 s^???1 (for Lg). The recently recognized discrepancy between the values of Q measured from long-period surface waves and normal-mode oscillations could also be explained by a slight positive bias in the geometrical spreading of surface waves. Similarly to the apparent χ, the corresponding linear variation with frequency is inferred for the intrinsic attenuation coefficient, χ _i , which combines the effects of geometrical spreading and dissipation within the medium. Frequency-dependent rheological or scattering Q is not required for explaining any of the attenuation observations considered in this study. The often-interpreted increase of Q with frequency may be apparent and caused by using the Q-based model of attenuation and following preferred Q( f) dependences while ignoring the true χ( f) trends within the individual frequency bands.     

Regional characteristics of coda attenuation after the lancang-Gengma earthquake

In the light of the single scattering model of coda originating from local earthquakes, and based on the aftershock coda registered respectively at the 4 short period stations installed near the foci shortly after theM7.6 Lancang andM7.2 Gengma earthquakes, this paper has tentatively calculated the rate of amplitude attenuation and theQ _c-value of the coda in the Lancang and Gengma areas using a newly-founded synthetic determination method. Result of the study shows the rate of coda amplitude attenuation demonstrates remarkable regional differences respectively in the southern and northern areas. The southern area presents a faster attenuation (Q _c=114), whereas the northern area shows a slower attenuation (Q _c=231). The paper also discusses the reasons causing such differences. Result of the study also suggests a fairly good linear relation between the coda source factorA _o(f) and the seismic moment and the magnitude. Using the earthquake scaling law, the following formulas can be derived: lgM ₀=lgA ₀(f)+17.6,M _D=0.67lgA ₀(f)+1.21 and logM ₀=1.5M _D+15.79. In addition, the rates of amplitude attenuationβ _s andβ _m are respectively calculated using the single scattering and multiple scattering models, and the ratioβ _s/β_m=1.20−1.50 is found for the results respectively from the two models. Finally, the mean free pathL of the S-wave scattering in the southern and northern areas are determined to be 54 km and 122 km respectively by the relations which can distinguish between the inherentQ _i and scatteringQ _s, testify to this areas having lowQ-values correspond to stronger scatterings. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, 71–82, 1992. This study is partly supported by the Seismological Science Foundation of the State Seismological Bureau of China, and the present English version of the paper is translated from its Chinese original by Wenyi Xia, Seismological Bureau of Yunnan Province.