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.     

Influence of attenuation and site on microearthquakes’ spectra in Shillong region, of Northeast, India: A case study

We examine the influence of attenuation and site on the spectra of microearthquakes having origin within the Shillong region. The ratios of spectral amplitudes at lower and higher frequencies are measured for three different stations at varying epicentral distances to estimate Q value for both P- and S-wave in near and sub-surface layer. The average estimates of Q _P and Q _S are found to be 178 and 195. The ratio of Q _S to Q _P emerges to be greater than unity in major parts of the Shillong area, suggesting dominance of dry crust prevailing in Shillong region. The variation in corner frequencies for these spectra is inferred to be characteristics of the site. Besides, the disparity in spectral content with reference to hard rock site yields the inference that the incoming seismic signals get amplified considerably while traversing from southern part to northeastern part of Shillong, best outlined at 2 to 5 Hz, which is well corroborated by the existing lithology.     

Anisotropic P-wave attenuation measured from a multi-azimuth surface seismic reflection survey

A system of aligned vertical fractures produces azimuthal variations in stacking velocity and amplitude variation with offset, characteristics often reported in seismic reflection data for hydrocarbon exploration. Studies of associated attenuation anisotropy have been mostly theoretical, laboratory or vertical seismic profiling based. We used an 11 common‐midpoint‐long portion of each of four marine surface‐seismic reflection profiles, intersecting each other at 45° within circa 100 m of a common location, to measure the azimuthal variation of effective attenuation, Q⁻¹_eff and stacking velocity, in a shallow interval, about 100 m thick, in which consistently orientated vertical fracturing was expected due to an underlying salt diapirism. We found qualitative and quantitative consistency between the azimuthal variation in the attenuation and stacking velocity, and published amplitude variation with offset results. The 135° azimuth line showed the least apparent attenuation (1000 Q⁻¹_eff= 16 ± 7) and the fastest stacking velocity, hence we infer it to be closest to the fracture trend: the orthogonal 45° line showed the most apparent attenuation (1000Q⁻¹_eff= 52 ± 15) and slowest stacking velocity. The variation of Q⁻¹_eff with azimuth φ is well fitted by 1000Q⁻¹_eff = 34 − 18cos[2(φ+40°)] giving a fracture direction of 140 ± 23° (±1SD, derived from ‘bootstrapping’ fits to all 11⁴ combinations of individual common‐midpoint/azimuth measurements), compared to 134 ± 47° from published amplitude variation with offset data. The effects of short‐window spectral estimation and choices of spectral ratio bandwidth and offset ranges used in attenuation analysis, individually give uncertainties of up to ±13° in fracture direction. This magnitude of azimuthal variation can be produced by credible crack geometries (e.g., dry cracks, radius 6.5 m, aspect ratio 3 × 10⁻⁵, crack density 0.2) but we do not claim these to be the actual properties of the interval studied, because of the lack of well control (and its consequences for the choice of theoretical model and host rock physical properties) and the small number of azimuths available here.     

Body-wave dispersion in the Friuli focal area

The differential dispersion of P- and S- body waves is studied in the Friuli area. We estimate theoretically the respective contributions of the source, from dislocation models, and of the propagation, from the Futterman model. We demonstrate the possible existence of differential dispersion, which is interpreted as being due to attenuation. The quality factor Q_P deduced from this hypothesis shows a regional variation from 15 to 210. Results also indicate a possible dependence of Q_P on frequency in the range 5–40 Hz.     

Coda Q c Attenuation and Source Parameter Analysis in Friuli (NE Italy) and its Vicinity

—?The digital data acquired by 16 short-period seismic stations of the Friuli-Venezia-Giulia seismic network for 56 earthquakes of magnitude 2.3–4.7 which occurred in and near NE Italy have been used to estimate the coda attenuation Q _c and seismic source parameters. The entire area under study has been divided into five smaller regions, following a criterion of homogeneity in the geological characteristics and the constrains imposed by the distribution of available events. Standard IASPEI routines for coda Q _c determination have been used for the analysis of attenuation in the different regions showing a marked anomaly in the values measured across the NE border between Friuli and Austria for Q ₀ value. A large variation exists in the coda attenuation Q _c for different regions, indicating the presence of great heterogeneities in the crust and upper mantle of the region. The mean value of Q _c (f) increases from 154–203 at 1.5?Hz to 1947–2907 at 48?Hz frequency band with large standard deviation estimates.¶Using the same earthquake data, the seismic-moment, M ₀, source radius, r and stress-drop, Δσ for 54 earthquakes have been estimated from P- and S-wave spectra using the Brune's seismic source model. The earthquakes with higher stress-drop (greater than 1?Kbar) occur at depths ranging from 8 to 14?km.     

Attenuation of P and pP waves in Vrancea area – Romania

The seismic attenuation in the Vrancea region (Romania) is investigated from teleseismic recordings of P and pP waves during the four major, intermediate-depth Romanian events that occurred since the onset of digital instrumentation. Most stations are located in Canada and in the United States, being equipped with a variety of sensors, especially short-period ones. The amplitude spectral ratio method is used, assuming no frequency dependence of the Q _P factor in the range 0.2–2 Hz. No apparent correlation between the derived attenuation value and the type of recording sensor is observed. Lateral variations of the attenuation are obtained, with a very low Q _P area (values down to 33) located in the northwestern part of the Vrancea seismogenic volume. For the stations with different azimuth angles in relation to the epicentral area, Q _P values routinely exceed 200. Most likely, the low attenuation values are related to an upwelling mantle material located immediately beneath the crust, but limited in depth to at least 100 km.     

2-D Image of Seismic Attenuation beneath the Deep Seismic Sounding Profile "Quartz," Russia

—We present a 2-D image of the upper mantle attenuation using nuclear explosion data from the ultra-long refraction/reflection profile "Quartz." Our analysis is based on a modified common spectrum technique followed by least-squares inversion for Q and iterative ray tracing in the velocity structure obtained earlier. The resulting attenuation structure corroborates the earlier model for northern Eurasia, as well as our recent estimate based on the analysis of the long-range P _n phase, and provides significantly more detail than the existing models. The resulting upper mantle attenuation structure is characterised by Q values ranging from 400 to 1800. Down to the depths of 150–190, and probably 400 km, the attenuation increases horizontally in SE direction, away from the Baltic Shield. Our model exhibits strong 2-D, vertical and horizontal attenuation contrasts. A high-attenuation layer in the depth range of 120–150 to 160–180 km can apparently be associated with the presence of a partial melts within the base of the lithosphere.     

Time-domain forward and inverse modeling of lossy soils with frequency-independent Q for near-surface applications

In this paper, we are concerned with a full-waveform-based methodology that allows the simultaneous imaging of the soil's stiffness and attenuating properties, using solely the soil's surficial response to probing waves.To date, field observations of small-strain wave attenuation in geomaterials at moderate spatial scales suggest that a commonly used metric of intrinsic and apparent attenuation, the seismic quality factor Q, is frequency-independent for a wide part of the frequency spectrum, including the frequency range of interest to seismic applications. We discuss first the forward simulation of waves in near-surface soil deposits directly in the time-domain using simplified models that adequately approximate nearly frequency-independent Q. To this end, we first review various attenuation models that aim at reproducing the frequency-independent Q behavior, and conclude, supported by site analyses, that, even though a generalized Maxwell body with eight Maxwell elements in parallel (GMB8) provides the best fit to frequency-independent Q, we favor a version of it with fewer parameters (GMB2), in order to reduce modeling complexity, while still retaining good agreement with the GMB8 model.We report on forward site analyses that lend credence to the choice of the GMB2 simplified model. We, then, use the GMB2 constitutive relation in the context of full-waveform inversion, and report on numerical experiments that lead to the imaging of the soil's properties in heterogeneous semi-infinite domains.     

Seismic Attenuation Computed from GLIMPCE Reflection Data and Comparison with Refraction Results

—Instantaneous frequency matching has been used to compute differential t* values for seismic reflection data from the Great Lakes International Multidisciplinary Program on Crustal Evolution (GLIMPCE) experiment. The differential attenuation values were converted to apparent Q ^?1 models by a fitting procedure that simultaneously solves for the interval Q ^?1 values using non-negative least squares. The bootstrap method was then used to estimate the variance in the interval Q ^?1 models. The shallow Q ^?1 structure obtained from the seismic reflection data corresponds closely with an attenuation model derived using instantaneous frequency matching on seismic refraction data along the same transect. This suggests that the effects of wave propagation and scattering on the apparent attenuation are similar for the two data sets. The Q ^?1 model from the reflection data was then compared with the structural interpretation of the reflectivity data. The highest interval Q ^?1 values (>0.01) were found near the surface, corresponding to the sedimentary rock sequence of the upper Keweenawan. Low Q ^?1 values (<0.0006) are found beneath the Midcontinent rift’s central basin. In addition to structural interpretation, seismic attenuation models derived in this way can be used to correct reflection data for dispersion, frequency and amplitude effects, and allow for improved imaging of the subsurface.     

A Summary of Attenuation Measurements from Borehole Recordings of Earthquakes: The 10 Hz Transition Problem

—Earthquake seismograms recorded by instruments in deep boreholes have low levels of background noise and wide signal bandwidth. They have been used to extend our knowledge of crustal attenuation both in the near-surface and at seismogenic depths. Site effects are of major importance to seismic hazard estimation, and the comparison of surface, shallow and deep recordings allows direct determination of the attenuation in the near-surface. All studies to date have found that Q is very low in the near-surface (～ 10 in the upper 100 m), and increases rapidly with depth. Unlike site amplification, attenuation at shallow depths exhibits little dependence on rock-type. These observations are consistent with the opening of fractures under decreasing lithostatic pressure being the principal cause of the severe near-surface attenuation. Seismograms recorded in deep boreholes are relatively unaffected by near-surface effects, and thus can be used to measure crustal attenuation to higher frequencies (≥ 100 Hz) than surface recordings. Studies using both direct and coda waves recorded at over 2 km depth find Q to be high (～ 1000) at seismogenic depths in California, increasing only weakly with frequency between 10 and 100 Hz. Intrinsic attenuation appears to be the dominant mechanism. These observations contrast with those of the rapidly increasing Q with frequency determined from surface studies in the frequency range 1 to 10 Hz. Further work is necessary to constrain the factors responsible for this apparent change in the frequency dependence of Q, but it is clearly unwise to extrapolate Q estimates made below about 10 Hz to higher frequencies.     

Worldwide variations in the attenuative properties of the upper mantle as determined from spectral studies of short-period body waves

A worldwide study of short-period teleseismic body wave spectra shows that the high frequency falloff rates of spectra are correlated with the tectonic type of the source and receiver regions and with source depth. The data indicate, in a consistent manner, that the main cause for such variations is the lateral variation of Q in the upper mantle as well as change of Q with depth. Using the internal consistency checks provided by redundancies in the data set other effects such as crustal, site dependent distortion of the spectra, source effects and instrument non-linearity can be ruled out as significant factors influencing the $\text{t1}$ estimates obtained. The results indicate high attenuation in the upper mantle under tectonic regions and new oceans. Long-period regional attenuation studies indicate similar variations in mantle Q among the types of regions mentioned but yield significantly lower Q estimates in all areas. The short- and long-period attenuation results can be reconciled only by assuming a frequency dependent Q that increases with frequency along all types of paths, such that the relative differences in Q along various types of paths retain the same sign over the short- and long-period bands.     

Investigation on regional attenuation of Vrancea (Romania) intermediate-depth earthquakes

This paper aims at investigating possible regional attenuation patterns in the case of Vrancea(Romania) intermediate-depth earthquakes.Almost 500 pairs of horizontal components recorded during 13 intermediate-depth Vrancea earthquakes are employed in order to evaluate the regional attenuation patterns.The recordings are grouped according to the azimuth with regard to the Vrancea seismic source and subsequently,Q models are computed for each azimuthal zone assuming similar geometrical spreading.Moreover,the local soil amplification which was disregarded in a previous analysis performed for Vrancea intermediate-depth earthquakes is now clearly evaluated.The results show minor differences between the four regions situated in front of the Carpathian Mountains and considerable differences in attenuation of seismic waves between the forearc and backarc regions(with regard to the Carpathian Mountains).Consequently,an average Q model of the type Q(f) = 115×f~(1.25) is obtained for the four forearc regions,while a separate Q model of the type Q(f) = 70×f~(0.90) is computed for the backarc region.These results highlight the need to evaluate the seismic hazard of Romania by using ground motion models which take into account the different attenuation between the forearc/backarc regions.     

Q P andQ S in Beijing and Jianchuan,Yunnan areas

Events from the December 1982 Huairou County, Beijing, and the July 1982 Jianchuan, Yunnan earthquake series were recorded at one station in Beijing City and at four stations in Eryuan area, Yunnan, respectively. Dividing the spectra (for P and S waves) from the smaller events in the series by spectra (for corresponding P and S waves) from the larger events in cach region, we have determined the high frequency source spectral decay rate to be ω^?1in both cases. Through trial and error method, we are able to determine the appropriate constantQ _P andQ _S that corrects the individual spectra to the proper high frequency decay rate. It is found that aQ _P of 800 and aQ _S of 550 can adequately compensate for the attenuation of the waves in Beijing area. For Jianchuan, Yunnan area the corresponding values are 900 and 400.     

Estimation of the Intrinsic Absorption and Scattering Attenuation in Northeastern Venezuela (Southeastern Caribbean) Using Coda Waves

—Northeastern Venezuela has been studied in terms of coda wave attenuation using seismograms from local earthquakes recorded by a temporary short-period seismic network. The studied area has been separated into two subregions in order to investigate lateral variations in the attenuation parameters. Coda-Q ^?1 (Q _c ^?1) has been obtained using the single-scattering theory. The contribution of the intrinsic absorption (Q _i ^?1) and scattering (Q _s ^?1) to total attenuation (Q _t ^?1) has been estimated by means of a multiple lapse time window method, based on the hypothesis of multiple isotropic scattering with uniform distribution of scatterers. Results show significant spatial variations of attenuation the estimates for intermediate depth events and for shallow events present major differences. This fact may be related to different tectonic characteristics that may be due to the presence of the Lesser Antilles subduction zone, because the intermediate depth seismic zone may be coincident with the southern continuation of the subducting slab under the arc.     

Low CodaQcin the Epicentral Region of the 2001 Bhuj Earthquake ofMw7.7

On 26 January, 2001 (03:46:55,UT) a devastating intraplate earthquake of M_w 7.7 occurred in a region about 5 km NW of Bhachau, Gujarat (23.42°N, 70.23°E). The epicentral distribution of aftershocks defines a marked concentration along an E-W trending and southerly dipping (45°) zone covering an area of (60 × 40) km². The presence of high seismicity including two earthquakes of magnitudes exceeding 7.7 in the 200 years is presumed to have caused a higher level of shallow crustal heterogeneity in the Kutch area; a site lying in the seismic zone V (zone of the highest seismicity for potentially M8 earthquakes) on the seismic zoning map of India. Attenuation property of the medium around the epicentral area of the Bhuj earthquake covering a circular area of 61,500 km² with a radius of 140 km is studied by estimating the coda-Q_c from 200 local earthquakes of magnitudes varying from 3.0–4.6. The estimated Q₀ values at locations in the aftershock zone (high seismicity) are found to be low in comparison to areas at a distance from it. This can be attributed to the fact that seismic waves are highly scattered for paths through the seismically active and fractured zone but they are well behaved outside the aftershock zone. Distribution of Q₀ values suggests that the local variation in Q₀ values is probably controlled by local geology. The estimated Q₀ values at different stations suggest a low value of Q=(102 ± 0.80)^*f^{(0.98 ± 0.02)} indicating an attenuative crust beneath the entire region. The frequency-dependent relation indicates a relatively low Q_c at lower frequencies (1–3 Hz) that can be attributed to the loss of energy due to scattering attenuation associated with heterogeneities and/or intrinsic attenuation due to fluid movement in the fault zone and fluid-filled cracks. The large Q_c at higher frequencies may be related to the propagation of backscattered body waves through deeper parts of the lithosphere where less heterogeneity is expected. Based on the attenuation curve estimated for Q₀=102, the ground acceleration at 240 km distance is 13% of 1 g i.e., 0.13 g agreeing well with the ground acceleration recorded by an accelerograph at Ahmedabad (0.11 g). Hence, it is inferred that the Q₀ value obtained from this study seems to be apt for prediction of ground motion for 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.     

Q β structure of the crust and upper mantle in the eastern Sino-Korean paraplatform

Based on the long period surface wave data recorded by the China Digital Seismograph Network (CDSN), theQ _R of fundamental mode Rayleigh wave with periods from 10 s to 146 s is determined for the eastern Sino-Korean paraplatform in this paper. TheQ _β models of the crust and upper mantle are respectively obtained for the 4 paths, with the aid of stochastic inverse method. It shows that in the eastern Sino-Korean paraplatform, the average crustalQ _β is about 200, and that there exists a weak attenuation layer in the middle crust (about 10–20 km deep) which is possibly related to earthquake-prone layer. A strong attenuation layer (lowQ) of 70 km thick extensively exists in the uppermost mantle, with the buried depth about 80 km. The averageQ _R of fundamental mode Rayleigh wave is between the value of stable tectonic region and that of active tectonic region, and much close to the latter.     

利用ML振幅研究地壳横波Q值Ⅱ:Q横向变化特征

利用中国地震年报(1984～2004年)测定M_L所用的54000余条区域地震最大振幅数据,采用二维层析成像方法反演中国大陆地壳介质品质因子Q₀结构.反演结果表明,中国大陆地壳介质品质因子Q₀平均为370,Q₀横向变化量可达±200.在整个研究区解的分辨率能够达到3°×3°,大部分地区能够达到2°×2°.Q₀的变化显示出与地质构造分区相关的特征：低衰减 (高Q₀值) 异常区主要分布在塔里木盆地、四川盆地、鄂尔多斯高原及华南地区,鄂尔多斯高原北部Q₀值最高;高衰减异常区(低Q₀值)主要分布在青藏高原以东的川滇地区、华北平原、渤海湾,最低Q₀值在渤海湾地区,天山、祁连山和汾渭地堑Q₀也呈较低异常变化.总体看来,构造上稳定的克拉通地台的介质品质因子Q₀较高,而构造上活动的地区往往Q₀较低,显现出区域的衰减变化特征与地壳结构和区域构造之间的密切关系.     

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.     

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