Research on codaQ distribution in Beijing and its surrounding regions

Based on the single scattering model of seismic coda waves, we have calculated theQ-factor in Beijing and its surrounding regions by means of calculating the power density spectrum in frequency domain with a fixed time window. The digital seismic data of 69 earthquakes from Beijing Telemetered Seismographic Network are used. These earthquakes were recorded from January 1, 1989 to December 31, 1990 at 20 stations. This paper shows the variations of the codaQ-factors in the studied region with different sites, frequency and lapse time, and the temporal change of the codaQ-factors in these two years. The results indicate that codaQ-factor depends strongly on the lapse time and frequency. It is assumed that whenQ _C=Q ₀f^η, for the three time windows of 15–30s, 30–60s and 60–90s, the average values ofQ ₀ are 48, 115 and 217; and the average values ofη are 0. 89, 0.91 and 0.74, respectively.     

Research on the regionalQ values of coda waves for the aftershocks in Luquan,Yunnan of 1985

Based on the scattering coda model by which local and regional earthquakes are interpreted (K. Aki, 1969), and using observational coda data of 68 aftershocks of the 1985 Luquan, Yunnan earthquake registered by the VGK seismographs installed at 12 stations in the Yunnan regional short-period network, theQ-values of coda waves are calculated respectively for 6 time intervals. It is observed that within the frequency range of 0.40–1.65 Hz of the observed data, theQ-values are closely related with the frequencies and the calculated codaQ ranges between 80–240 with the coefficient of frequency dependence η=0.45. The calculated source factorsB(f> _p) of the coda waves which indicate the scattering strength are mostly within the order 10^?23–10^?24. Areas with lowQ-values present high scattering. It should be noted that by comparing data obtained before and after the Luquan earthquake, clear changes can be detected in theQ-values measured at stations close to the epicentral region, and that theQ-values of the aftershock coda are less than about one half of the pre-shock values. It may be mentioned that the time-dependent regional variations of theQ-values might possibly bring about practical significance in earthquake prediction. Moreover, aftershock focal parameters are determined. Through discussions on the quantitative relations between the focal parameters, we get: 1gE=1.59M _L+ 11.335;E=(2.10 × 10^?5)M ₀; length of focal rupturea=0.40?0.80 km for 3.0≤M _L<5.0 events; stress drop Δσ=(6.0–130) ×10⁵ Pa. Through interpretation of the data, we have also learned the important characteristics that there is no linear relation between the stress drops and the earthquake magnitudes.     

Comparison of amplitude decay rates in reflection,refraction, and local earthquake records

Three types of seismic data recorded near Coalinga, California were analyzed to study the behavior of scattered waves: 1) aftershocks of the May 2, 1983 earthquake, recorded on verticalcomponent seismometers deployed by the USGS; 2) regional refraction profiles using large explosive sources recorded on essentially the same arrays above; 3) three common-midpoint (CMP) reflection surveys recorded with vibrator sources over the same area. Records from each data set were bandpassed filtered into 5 Hz wide passbands (over the range of 1–25 Hz), corrected for geometric spreading, and fit with an exponential model of amplitude decay. Decay rates were expressed in terms of inverse codaQ (Q _c ^–1 ).Q _c ^–1 values for earthquake and refraction data are generally comparable and show a slight decrease with increasing frequency. Decay rates for different source types recorded on proximate receivers show similar results, with one notable exception. One set of aftershocks shows an increase ofQ _c ^–1 with frequency.Where the amplitude decay rates of surface and buried sources are similar, the coda decay results are consistent with other studies suggesting the importance of upper crustal scattering in the formation of coda. Differences in the variation ofQ _c ^–1 with frequency can be correlated with differences in geologic structure near the source region, as revealed by CMP-stacked reflection data. A more detailed assessment of effects such as the depth dependence of scattered contributions to the coda and the role of intrinsic attenuation requires precise control of source-receiver field geometry and the study of synthetic seismic data calculated for velocity models developed from CMP reflection data.     

Lateral Variations of Lg Coda Q in Southern Mexico

—Broad band digital three-component data recorded at UNM, a GEOSCOPE station, were used to estimate Lg coda Q for 34 medium size (3.9 ≤m _b ≤ 6.3) earthquakes with travel paths laying in different geological provinces of southern Mexico in an effort to establish the possible existence of geological structures acting as wave guides and/or travel paths of low attenuation between the Pacific coast and the Valley of Mexico. The stacked spectral ratio method proposed by XIE and NUTTLI (1988) was chosen for computing the coda Q. The variation range of Q ₀ (Q at 1?Hz) and the frequency dependence parameter η estimates averaged on the frequency interval of 0.5 to 2?Hz for the regions and the three components considered are: i) Guerrero region 173 ≤Q¯ ₀≤ 182 and 0.6 ≤Q¯ ₀ ≤ 0.7, ii) Oaxaca region 183 ≤Q¯ ₀ ≤ 198 and 0.6 ≤Q¯ ₀ ≤ 0.8, iii) Michoacan-Jalisco region 187 ≤Q¯ ₀ ≤ 204 and 0.7 ≤Q¯ ₀ ≤ 0.8 and iv) eastern portion of the Transmexican Volcanic Belt (TMVB) 313 ≤Q ₀≤ 335 and η = 0.9. ¶The results show a very high coda Q for the TMVB as compared to other regions of southern Mexico. This unexpected result is difficult to reconcile with the geophysical characteristics of the TMVB, e.g., low seismicity, high volcanic activity and high heat flow typical of a highly attenuating (low Q) region. Visual inspection of seismograms indicates that for earthquakes with seismic waves traveling along the TMVB, the amplitude decay of Lg coda is anomalously slow as compared to other earthquakes in southern Mexico. Thus, it seems that the high Q value found does not entirely reflect the attenuation characteristics of the TMVB but it is probably contaminated by a wave-guide effect. This phenomenon produces an enhancement in the time duration of the Lg wave trains travelling along this geological structure. This result is important to establish the role played by the transmission medium in the extremely long duration of ground motion observed during the September 19, 1985 Michoacan earthquake. ¶The overall spatial distribution of coda Q values indicates that events with focus in the Michoacan-Jalisco and Oaxaca regions yield slightly higher values than those from Guerrero. This feature is more pronounced for the horizontal component of coda Q. A slight dependence of average coda Q ^?1 on earthquake focal depth is observed in the frequency range of 0.2 to 1.0?Hz approximately on the horizontal component. Deeper (h > 50?km) events yield lower values of Q ^?1 than shallower events. For frequencies higher than 1.0?Hz no clear dependence of Q ^?1 on focal depth is observed. However, due to the estimates uncertainties this result is not clearly established.     

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.     

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.     

Spatial variation of crustal codaQ in California

Coda wave data from California microearthquakes were studied in order to delineate regional fluctuations of apparent crustal attenuation in the band 1.5 to 24 Hz. Apparent attenuation was estimated using a single back scattering model of coda waves. The coda wave data were restricted to 30 s following the origin time; this insures that crustal effects dominate the results as the backscattered shear waves thought to form the coda would not have had time to penetrate much deeper. Results indicate a strong variation in apparent crustal attenuation at high frequencies between the Franciscan and Salinian regions of central California and the Long Valley area of the Sierra Nevada. Although the codaQ measurements coincide at 1.5 Hz (Q _c =100), at 24 Hz there is a factor of four difference between the measurements made in Franciscan (Q _c =525) and Long Valley (Q _c =2100) with the Salinian midway between (Q _c =900). These are extremely large variations compared to measures of seismic velocities of comparable resolution, demonstrating the exceptional sensitivity of the high frequency codaQ measurement to regional geology. In addition, the frequency trend of the results is opposite to that seen in a compilation of codaQ measurements made worldwide by other authors which tend to converge at high and diverge at low frequencies, however, the worldwide results generally were obtained without limiting the coda lengths and probably reflect upper mantle rather than crustal properties. Our results match those expected due to scattering in random media represented by Von Karman autocorrelation functions of orders 1/2 to 1/3. The Von Karman medium of order 1/3 corresponding to the Franciscan codaQ measurement contains greater amounts of high wavenumber fluctuations. This indicates relatively large medium fluctuations with wavelengths on the order of 100 m in the highly deformed crust associated with the Franciscan, however, the influence of scattering on the codaQ measurement is currently a matter of controversy.     

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.     

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.     

Research on coda Q distribution in Beijing and its surrounding regions

Based on the single scattering model of seismic coda waves, we have calculated theQ-factor in Beijing and its surrounding regions by means of calculating the power density spectrum in frequency domain with a fixed time window. The digital seismic data of 69 earthquakes from Beijing Telemetered Seismographic Network are used. These earthquakes were recorded from January 1, 1989 to December 31, 1990 at 20 stations. This paper shows the variations of the codaQ-factors in the studied region with different sites, frequency and lapse time, and the temporal change of the codaQ-factors in these two years. The results indicate that codaQ-factor depends strongly on the lapse time and frequency. It is assumed that whenQ _C=Q ₀f^η, for the three time windows of 15–30s, 30–60s and 60–90s, the average values ofQ ₀ are 48, 115 and 217; and the average values ofη are 0. 89, 0.91 and 0.74, respectively. Contribution No. 95A0009, Institute of Geophysics, SSB, China. This work is a contract subject 85-04-01-02 of the State Seismological Bureau, China.     

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.     

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.     

<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.     

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.     

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.     

CodaQ as a combination of scattering and intrinsic attenuation: Numerical simulations with the boundary integral method

Numerical modelling ofSH wave seismograms in media whose material properties are prescribed by a random distribution of many perfectly elastic cavities and by intrinsic absorption of seismic energy (anelasticity) demonstrates that the main characteristics of the coda waves, namely amplitude decay and duration, are well described by singly scattered waves in anelastic media rather than by multiply scattered waves in either elastic or anelastic media. We use the Boundary Integral scheme developed byBenites et al. (1992) to compute the complete wave field and measure the values of the direct waveQ and coda wavesQ in a wide range of frequencies, determining the spatial decay of the direct wave log-amplitude relation and the temporal decay of the coda envelope, respectively. The effects of both intrinsic absorption and pure scattering on the overall attenuation can be quantified separately by computing theQ values for corresponding models with (anelastic) and without (elastic) absorption. For the models considered in this study, the values of codaQ ^–1 in anelastic media are in good agreement with the sum of the corresponding scatteringQ ^–1 and intrinsicQ ^–1 values, as established by the single-scattering model ofAki andChouet (1975). Also, for the same random model with intrinsic absorption it appears that the singly scattered waves propagate without significant loss of energy as compared with the multiply scattered waves, which are strongly affected by absorption, suggesting its dominant role in the attenuation of coda waves.     

Variation of Seismic Coda Wave Attenuation in the Garhwal Region, Northwestern Himalaya

Seismic coda wave attenuation ( $ Q_{\text{c}}^{ - 1} $ ) characteristics in the Garhwal region, northwestern Himalaya is studied using 113 short-period, vertical component seismic observations from local events with hypocentral distance less than 250?km and magnitude range between 1.0 to 4.0. They are located mainly in the vicinity of the Main Boundary Thrust (MBT) and the Main Central Thrust (MCT), which are well-defined tectonic discontinuities in the Himalayas. Coda wave attenuation ( $ Q_{\text{c}}^{ - 1} $ ) is estimated using the single isotropic scattering method at central frequencies 1.5, 3, 5, 7, 9, 12, 16, 20, 24 and 28?Hz using several starting lapse times and coda window lengths for the analysis. Results show that the ( $ Q_{\text{c}}^{ - 1} $ ) values are frequency dependent in the considered frequency range, and they fit the frequency power law ( $ Q_{\text{c}}^{ - 1} \left( f \right) = Q_{0}^{ - 1} f^{ - n} $ ). The Q ₀ (Q _c at 1?Hz) estimates vary from about 50 for a 10?s lapse time and 10?s window length, to about 350 for a 60?s lapse time and 60?s window length combination. The exponent of the frequency dependence law, n ranges from 1.2 to 0.7; however, it is greater than 0.8, in general, which correlates well with the values obtained in other seismically and tectonically active and highly heterogeneous regions. The attenuation in the Garhwal region is found to be lower than the Q _c ^?1 values obtained for other seismically active regions of the world; however, it is comparable to other regions of India. The spatial variation of coda attenuation indicates that the level of heterogeneity decreases with increasing depth. The variation of coda attenuation has been estimated for different lapse time and window length combinations to observe the effect with depth and it indicates that the upper lithosphere is more active seismically as compared to the lower lithosphere and the heterogeneity decreases with increasing depth.     

Attenuation of coda waves in northern Greece

The single scattering model has been applied for the estimation of codaQ values for local earthquakes that occurred in northern Greece during the period 1983–1989 and recorded by the telemetered network of the Geophysical Laboratory of the University of Thessaloniki. CodaQ estimations were made for four frequency bands centered at 1.5 Hz, 3.0 Hz, 6.0 Hz and 12.0 Hz and for the lapse time windows 10–20 sec, 15–30 sec, 20–45 sec, 30–60 sec and 50–100 sec. The codaQ values obtained show a clear frequency dependence of the formQ _c =Q ₀ f ⁿ , whileQ ₀ andn depend on the lapse time window.Q ₀ was found equal to 33 andn equal to 1.01 for the time window of 10 to 20 sec, while for the other windowsQ ₀ increased from 60 to 129, withn being stable, close to 0.75. This lapse time dependence is interpreted as due to a depth dependent attenuation. The high attenuation and the strong frequency dependence found are characteristic of an area with high seismicity, in agreement with studies in other seismic regions.     

Attenuation of coda waves in the contact zone between the Dinarides and the Adriatic Microplate

The attenuation of coda waves is analysed for nine seismic stations in the area of convergent motion of the Adriatic microplate and the Dinarides. The frequency dependent coda quality factor of the form Q_c = Q₀ fⁿ is estimated for up to seven central frequencies (1.5, 3, 6, 9, 12, 18 and 24 Hz) and for 21 successive 30 s long time windows. Q₀ was found to increase from 68–353 for short lapse times of 20–50 s, to 158–373 for lapse times of 90–100 s. Parameter n is observed to vary between 0.46 and 0.89, with a pronounced tendency to decrease with increasing Q₀, and vice versa. Both Q₀ and n seem to stabilize for lapse times larger than 50–80 s, indicating that a change in rock properties controlling coda attenuation occurs at depths of about 100–160 km. The spatial distribution of observed Q₀ is well correlated with observed seismicity and inferred tectonic activity. In particular, we observe significant negative correlation of Q₀ with the peak ground acceleration (PGA) estimate for the return period of 475 years. The degree of frequency dependence n, is the smallest for stations on the islands, where the crust is the thinnest. The largest n is observed over the thickest crust in the region, where the Moho lies at depths of over 55 km.