Scattering attenuation and the fractal geometry of fracture systems

Scattering of seismic waves can be shown to have a frequency dependenceQ ^{–1 3–v} if scattering is produced by arrays of inhomogeneities with a 3D power spectrumW _3D(k) k ^–v. In the earth's crust and upper mantle the total attenuation is often dominated by scattering rather than intrinsic absorption, and is found to be frequency dependent according toQ ^–1 , where –1<–0.5. IfD ₁ is the fractal dimension of the surface of the 3D inhomogeneities measured on a 2D section, then this corresponds respectively to 1.5<D ₁1.75, since it can be shown that =2(D ₁–2). Laboratory results show that such a distribution of inhomogeneities, if due to microcracking, can be produced only at low stress intensities and slow crack velocities controlled by stress corrosion reactions. Thus it is likely that the earth's brittle crust is pervaded by tensile microcracks, at least partially filled by a chemically active fluid, and preferentially aligned parallel to the maximum principal compressive stress. The possibility of stress corrosion implies that microcracks may grow under conditions which are very sensitive to pre-existing heterogeneities in material constants, and hence it may be difficult in practice to separate the relative contribution of crack-induced heterogeneity from more permanent geological heterogeneities.By constrast, shear faults formed by dynamic rupture at critical stress intensities produceD ₁=1, consistent with a dynamic rupture criterion for a power law distribution of fault lengths with negative exponentD. The results presented here suggest empirically thatD ₁-1/2(D+1), thereby providing the basis for a possible framework to unify the interpretation of temporal variations in seismicb-value (b-D/2) and the frequency dependence of scattering attenuation ().This is PRIS contribution 046.     

Is the single scattering model invalid for the coda excitation at long lapse time?

Supposing that the distribution of scatterers in a three-dimensional medium is not uniform but fractally homogeneous with fractal dimensionD, we have made the dimensional analysis for the temporal decay of the multiple scattering energy density at the hypocenter.The number of scatters in a sphere of radiusR is assumed to be proportional toR ^D . Then, the energy density of thekth order scattering decays according to the [(D–2)k–3]th power of lapse time. A fractal dimension ofD=3 corresponds to the uniform distribution. If 2<D3, multiple scattering terms of orderk2 dominate over the single scattering term (k=1) at long lapse time. IfD=2, energy density of every order decays according to the — 3rd power of lapse time. The single scattering model survives on conditionD<2; the single scattering term dominates over the higher order multiple scattering terms even at long lapse time, since the negative power of lapse time fork=1 is the smallest of all.     

Unified approach to amplitude attenuation and coda excitation in the randomly inhomogeneous lithosphere

A unified model is proposed for explaining the frequency dependent amplitude attenuation and the coda wave excitation on the basis of the single scattering process in the randomly inhomogeneous lithosphere. Adopting Birch's law and a direct proportion between density and wave velocity, we statistically describe the inhomogeneous medium by one random function characterized by the von Karman autocorrelation function. We calculate the amplitude attenuation from the solid angle integral of scattered wave energy on the basis of the Born approxiimation after subtracting the travel-time fluctuation effect caused by slowly varying velocity inhomogeneities. This subtraction is equivalent to neglect energy loss by scattering within a cone around the forward direction. The random inhomogeneity of the von Karman autocorrelation function of order 0.35 with the mean square fractional fluctuation of 7.2×10^–3 1.3×10^–2 and the correlation distance of 2.15.1 km well explains observed backward scattering coefficientg and the ratioQ _P ^–1 /Q _S ^–1 , and observed and partially conjecturedQ _S ^–1 for frequencies between 0.5 Hz and 30 Hz.     

Fractal approach of the temporal earthquake distribution in the Hellenic arc-trench system

The time clustering of earthquakes occurring in the Hellenic arc-trench system is quantitatively analyzed by means of the fractal dimension,D, of their time distribution in the time intervals of 1950–1985 (M _s >-4.5) and 1964–1985 (M _s 4.0). The results obtained imply that scale-invariant clustering holds over very large scale lengths of time,T, with 2²–2⁸T (in min) 2²⁰–2²², depending on the seismotectonic segment considered. In all segments a common feature is the relation between theD ₁,D ₂ andD ₃-values found for shallow, intermediate-depth and all-depth shocks, respectively:D ₃>D₁>D₂. TheD-values found for shallow shocks range between 0.137 and 0.191 with the exception of the Ionian Islands and Cretan segments where anomalously high values (D=0.221–0.251) have been determined. We discuss possible seismotectonic interpretations of the results.     

Particle size distribution of cataclastic fault materials from Southern California: A 3-D study

The particle size distributions of fault gouge from the San Andreas, the San Gabriel, and the Lopez Canyon faults in Southern California were measured using sieving and Coulter-Counter techniques over a range of particle sizes from 2 m to 16 mm. The distributions were found to be power law (fractal) for the smaller fragments and log-normal by mass for sizes near and above the peak size. The apparent fractal dimensionD of the smaller particles in gouge samples from the San Andreas fault, the San Gabriel fault and the Lopez Canyon gouge were 2.4–3.6, 2.6–2.9 and 2.4–3.0, respectively. The averageD for the Lopez Canyon gouge was 2.7±0.2, which is in agreement with earlier studies of this gouge using planar 2-D sections. The fractal dimension of the finer fragments from all three faults is observed to be correlated with the peak fragment size, with finer gouges tending to have a largerD. A computer automaton is used to show that this observation may be explained as resulting from a fragmentation process which has a grinding limit at which particle reduction stops.     

Microtopographic evolution of mineral surfaces as a tool to identify and date young fault scarps in bedrock

Faulting that results in surface ruptures through bedrock can be particularly difficult to date. For example, stratigraphic control on the age of faulting, based on the age of the bedrock, often leaves unacceptably large uncertainty on the age of the faulting. From a paleoseismological perspective, there is a clear need to determine if a bedrock fault scarp is actually a young feature. For young fault ruptures that create fresh mineral surfaces, analysis of microtopography developed by weathering of the mineral surface may provide a quantifiable method for determining the fault age. The direct quantitative measurement of mineral surface microtopography using Atomic Force Microscopy affords a novel method to study the rupture ages of active faults. The method for using microtopographic evolution of mineral surfaces depends on three conditions. The first condition is that freshly exposed mineral cleavage surfaces, which can be described geometrically as planes, are formed during a rupture event. The formation of these fresh surfaces is analogous to the initiation of a weathering ‘clock’ that defines time t=0. Following cleavage formation dissolution of the planar mineral surface occurs. The rate of dissolution for a mineral species under given climatic conditions, governs the rate of mineral surface alteration. Thus as dissolution proceeds, the roughness of the mineral surface increases. We suggest that the progression of microtopographic roughness over time, which can be estimated by computing quantitative statistics derived from digital mineral surface topography, will systematically vary until a steady state surface topography is reached. The fractal dimension, D_f, is one such measure of surface roughness where, D_f at time t=0 is 2. The dissolution of the mineral surface increases the fractal dimension as the removal of material proceeds. We posit that somewhere between D_f=2 and D_f=3, the microtopography reaches a steady state. Therefore, in the pre-steady state stage of surface roughness, the quantitative measure of roughness of the mineral may serve as a measure of time elapsed since faulting. The period of time this initial stage of surface roughening represents is dependent on the mineral and as a consequence, its dissolution rate, in a specific set of environmental conditions. The time elapsed since fault rupture and grain cleavage can also be estimated from the measurement of the volume of material removed through dissolution. If part of the original cleavage surface remains and can be identified then AFM measurements of the surface microtopography can be used to calculate the dissolved volume per unit area.     

Clustering properties of Etna seismicity during 1981–1991

Space and time clustering properties ofseismic activity, affecting Etna Volcano (Italy)during 1981–1991, are investigated by fractaldimension analysis. Very interesting volcanic andseismic activity occurred within this time interval.Temporal evolution of the time fractal dimension D _t calculated on a moving window, revealscorrelation with the eruptive processes at differenttime scales confirming results obtained for a differenttime span (De Rubeis et al., 1997). Spatial fractaldimension D _s shows to be negativelycorrelated with the time fractal dimension D _t, suggesting a peculiar dynamic patternassociated with volcanic processes.     

Fractal research of fault gouge

ＦｒａｃｔａｌｒｅｓｅａｒｃｈｏｆｆａｕｌｔｇｏｕｇｅＳＨＵＮ－ＭＥＩＳＨＡＯ（邵顺妹）ａｎｄＪＩＮ－ＣＨＡＮＧＺＯＵ（邹瑾敞）ＥａｒｔｈｑｕａｋｅＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅｏｆＬａｎｚｈｏｕ，ＳｔａｔｅＳｅｉｓｍｏｌｏｇｉｃａｌＢｕｒｅａ...     

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.     

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.     

Scattering and anelastic attenuation of seismic energy in northern Greece

The relative contribution of scattering (Q _s ^–1 ) and intrinsic (Q _i ^–1 ) attenuation to the totalS-wave attenuation for the frequencies of 1.5, 3.0, 6.0 and 12.0 Hz has been studied by applying the radiative energy transfer theory, Data of local earthquakes which occurred in northern Greece and were recorded by the permanent telementered network of the Geophysical Laboratory of the University of Thessaloniki have been used. The results show that in this area the scattering attenuation is dominant over all frequencies while intrinsic attenuation is significantly lower. The estimatedQ _s ^–1 andQ _i ^–1 values have frequency dependences off ^–0.72 andf ^–0.45, respectively. The frequency dependence ofQ _s ^–1 is the same as that of the codaQ _c ^–1 , obtained by applying the single scattering model, which probably implies that the frequency dependence of the coda wave attenuation is attributed to the frequency dependence of the scattering attenuation.Q _c ^–1 values are very close to scattering attenuation for short lapse times, (10–20 sec), and intermediate between scattering and intrinsic attenuation for the longer lapse times, (50–100 sec). This difference is explained as the result of the depth-dependent attenuation properties and the multiple scattering effects.     

Fractal dimension of fault systems in Japan: Fractal structure in rock fracture geometry at various scales

Based on fault maps, whether or not the fracture geometry of rocks is self-similar, was examined by using a box-counting algorithm. The statistical self-similarity (fractal structure) of the fault fracture systems holds well at the scale of about 2 to 20 km. The fractal dimension in Japan varied from 1.05 to 1.60. The fractal dimension is about 1.5–1.6 at the central part of the Japan Arc, and decreases with distance from the center. At a smaller scale, the fractal structure also holds well in the rock fracture geometry. The fractal dimension of the North Izu Peninsula fault system (branching faults) is 1.49 at the scale of 0.625 to 10 km, the fractal dimension of rock fracture geometry at the scale order of 10^–1 to 10^–2 meters is about 1.49–1.61. The upper limit of the fractal dimension of rock fracture geometry is about 1.6, judging from the estimation of fractal dimension on actual fracture geometry of rocks. This value may impose a restraint on modeling of faulting and the fracture process of rocks.     

Changes in the early part of the seismic coda due to localized scatterers: The estimation ofQ in a stope environment

A single scattering model was used to analyse the temporary changes in the mean density of scattered waves in a discrete random medium. The model of the mean energy density, originally proposed bySato (1977) for spherical radiation and isotropic scattering, has been modified and applied to a medium in which the scatterers are confined to a specified volume. The time variation of the early part of the mean energy density function for the different source durations was investigated. The dominant effect on the theoretical mean energy density is caused by the specified volume containing scatterers. The duration of the source pulse influences the early part of the coda fort/t ₀<1.2, wheret is the lapse time measured from the source origin time, andt ₀is arrival time of the body wave.The analysis of the coda signal of micro-events occurring immediately in front of the face enables us to estimate the size of the fracture zone induced by the stope. The model of the mean energy density of coda for a medium containing scatterers close to the seismic source was used to analyse a large number of events recorded close to an advancing mine face in a deep level gold mine in South Africa. The coda decay rate has two trends: the first, with a steep decay of coda, is produced by a larger deviation of rock parameters and/or larger size of the scatterers; the second trend, which decays more slowly, has the corresponding mean-free path ranging from 20 m to 200 m. The analysis indicates that the rock mass about 15–20 m from the stope contains a large proportion of fractured and blocked rock, which is the source of scattering. The scattering of theS-wave was much stronger and more stable, with the mean-free path varying from 11 m to 45 m. This is due to the shorter wavelength of theS wave in comparison with theP wave. The quality factor for theP coda wave varies from 30 to 100 in the fracture zone of stope and outside this zone it has a value of 300. The quality factor of theS wave varies from 20 to 78 in the equivalent volume. For rock surrounding the stope the ratioQ _sp ^–1 /Q _ss ^–1 varied from 0.31 to 0.69. This suggests that the radii of scatterers are smaller than 3.5 m.     

Observations of teleseismicP wave coda for underground explosions

The earlyP wave coda (5–15 sec after the first arrival) of underground explosions at the Nevada Test Site is studied in the time domain using 2082 teleseismic short-period recordings, with the intent of identifying near-source contributions to the signals in the frequency range 0.2–2.0 Hz. Smaller magnitude events tend to have relatively high coda levels in the 0.4–0.8 Hz frequency band for both Yucca Flat and Pahute Mesa explosions. Coda complexity in this low-frequency passband is negatively correlated with burial depth for Pahute Mesa events but is only weakly correlated with depth for Yucca Flat events. Enhanced excitation of relatively long-period scattered waves for smaller, less deeply buried events is required to explain this behavior. Coda complexity in the 0.8–1.1 Hz band is positively correlated with magnitude and depth for Pahute Mesa events, but has no such dependence for Yucca Flat events. This may result from systematic variations between the spectra of direct signals and coda arrivals caused bypP interference for the largest events, all of which were detonated at Pahute Mesa. Another possible explanation is a frequency-dependent propagation effect on the direct signals of the larger events, most of which were located in the center of the mesa overlying strong lateral velocity gradients in the crust and upper mantle. Event average complexity varies spatially for both test sites, particularly in the 0.8–1.1 Hz band, providing evidence for frequency-dependent focussing or scattering by near-source structure. Both the direct arrivals and the early coda have strong azimuthal amplitude patterns that are produced by defocussing by mantle heterogeneity. The direct arrivals have stronger coherent azimuthal patterns than the early coda for Pahute Mesa events, indicating more pronounced deep crustal and shallow mantle defocussing for the direct signals. However, for Yucca Flat events the direct arrivals have less coherent azimuthal patterns than the coda, suggesting that a highly variable component of near-source scattering preferentially affecting the downgoing energy is superimposed on a pattern produced by mantle heterogeneity that affects the entire signal. This complicated behavior of the direct arrivals may be the result of triplications and caustics produced by the complex basement structure known to underlie the Yucca Flat test site. The presence of strong azimuthal patterns in the early coda indicates that source strength estimates based on early coda are subject to biases similar to those affecting estimates based on direct arrivals.     

Attenuation characteristics of coda waves and estimation of Q c values in eastern China

For short-period near-earthquake records in eastern China, from the empirical attenuation formula of coda ground motion amplitudeA with timeτ: lgA=G−2. 235 lgτ, using the single scattering theory modified with epicentral distance, we obtain the curve family of corrected coda amplitudeA _c(r,t), andω/2Q _c values for each time interval of coda. From this,Q _c(f,h) values, which correspond to each observational average frequency and sampling depth, are calculated. The results substantially agree with those observationalQ _c values in Yunnan, Beijing and central Asia.     

A model for attenuation and scattering in the Earth's crust

The mechanisms contributing to the attenuation of earthquake ground motion in the distance range of 10 to 200 km are studied with the aid of laboratory data, coda wavesRg attenuation, strong motion attenuation measurements in the northeast United States and Canada, and theoretical models. The frequency range 1–10 Hz has been studied. The relative contributions to attenuation of anelasticity of crustal rocks (constantQ), fluid flow and scattering are evaluated. Scattering is found to be strong with an albedoB ₀=0.8–0.9 and a scattering extinction length of 17–32 km. The albedo is defined as the ratio of the total extinction length to the scattering extinction length. TheRg results indicate thatQ increases with depth in the upper kilometer or two of the crust, at least in New England. CodaQ appears to be equivalent to intrinsic (anelastic)Q and indicates that thisQ increases with frequency asQ=Q _o f ⁿ , wheren is in the range of 0.2–0.9. The intrinsic attenuation in the crust can be explained by a high constantQ (500Q _o2000) and a frequency dependent mechanism most likely due to fluid effects in rocks and cracks. A fluid-flow attenuation model gives a frequency dependence (QQ _o f ^0.5) similar to those determined from the analysis of coda waves of regional seismograms.Q is low near the surface and high in the body of the crust.     

Temporal change in scattering and attenuation associated with the earthquake occurrence—A review of recent studies on coda waves

The study of coda waves has recently attracted increasing attention from seismologists. This is due to the fact that it is viewed as a new means by which the stress accumulation stage preceding a large earthquake can be measured, since the scattering paths nearly uniformly cover a fairly large region around the focus and observation stations, compared with the direct ray paths. To date, we have had many reports on the temporal variation of the relation between coda duration and amplitude magnitude, and that of the coda attenuationQ _c ^–1 which is estimated from coda amplitude decay. Some of these have shown a precursor-like behavior; however, others seem to have shown a coseismic change. We have critically reviewed these reports, and discussed what these observational facts tell us about the change in the heterogeneous crust. We found significant temporal variations, not only in the mean but also in the scatter ofQ _c ^–1 , associated with the mainshock occurrence. The formation of new cracks, the reopening and growing of existing cracks, the interaction of these cracks, and the pore water movement through these cracks might correspond to such variations. In addition, we may expect an inhomogeneous distribution of crack clusters in a fairly large region, compared with the aftershock region. The gradual appearance of such crack clusters seems to be the most plausible mechanism by which coda decay gradients are caused to largely scatter in the stress accumulation stage.     

Self-organized Fractal Seismicity of Reservoir Triggered Earthquakes in the Koyna-Warna Seismic Zone,Western India

— Analysis of the Koyna-Warna earthquake catalog (1968–1996) shows that on an average there is a positive correlation between the b value (decrease) and fractal dimensions (decrease in both D₂^s and D₂^t) of earthquake epicenters 0.5 and 2.5 years prior to 1973 (M5.2) and 1980 (M5.5) events, respectively, except a negative correlation for about five years (1988–1993) prior to the 1993/1994 sequence (M5.4). This positive correlation indicates a weaker clustering, or that the epicenters tend to fill the two-dimensional plane. While the origin of the negative correlation seems to be that during periods of large events (low b value), there is strong clustering around the main shock epicenter (high fractal dimension). Interestingly, during the last year (1995–1996) of the studied period both the b value and correlation dimensions rose significantly, suggesting that stress release occurs through increased levels of low magnitude and increasingly scattered seismicity, suggesting an increased risk of larger magnitude events. Incidentally, during 2000 three earthquakes of magnitude M 5.0, one earthquake of M 4.0, 45 earthquakes of magnitude M 3.0–3.9, and several thousand earthquakes of M < 3 have occurred in the region. Thus it can be inferred that at local scales the relationship yields both positive and negative correlation that appears to be controlled by different modes of failure within the active fault complex.Acknowledgement. The authors are grateful to Dr. B.K. Rastogi of NGRI for providing the catalog of Koyna earthquakes and for useful scientific discussions. The comments of Dr. I. G. Main have improved the quality of paper for which we extend to him our sincere thanks. One of the authors (AOM) thanks the Third World Academy of Science and the Council of Scientific and Industrial Research, India for the Postdoctoral Fellowship award under which this work was carried out.     

A study on the characteristics of the information dimension D 1 of the temporal and spatial distributions of earthquakes in an active fault zone

The research of the information dimension (D ₁) in an active fault zone considers the contribution of each seismic event to information and reflects the characteristics of the temporal and spatial distributions of earthquakes from a new point of view, avoiding some short-comings of the research about the capacity dimension (D ₀). The results of calculation show that the information dimension of the temporal distribution in Xianshuihe active fault zone before Luhuo large earthquake isD ₁=0.1051. It is a consult creterion of large earthquakes in future in the fault zone. The information dimensions of the temporal distribution of the earthquakes in Anninghe active fault zone are respectivelyD ₁(t _N)=0.1363 (for the north section) andD ₁(t _S)=0.06710 (for the south section). The information dimensions of the spatial distribution are respectivelyD ₁(K _N)=1.053 (for the north section) andD ₁(K _S)=0.7758 (for the south section). The north section and the south section belong to two self-similar systems with different information dimensions respectively. The extent of the self-organization of seismic activity in the south section is higher than that in the north section. This is helpful for us to judge the major dangerous section in the key region of the seismic monitoring. The research about the information dimension of the temporal and the spatial distributions of earthquakes is significant for the exploration of active fault zones and seismic prediction. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 372–379, 1991. This paper is sponsored by the Chinese Joint Seismological Science Foundation. The English version is improved by Zhenwen An.     

Lateral variations of coda Q and attenuation of seismic waves in Southwest Anatolia

The seismic quality factor (Q _c) and the attenuation coefficient (δ) in the earth’s crust in southwest (SW) Anatolia are estimated by using the coda wave method based on the decrease of coda wave amplitude by time on the seismogram. The quality factor Q _o, 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. δ is determined from the observations of amplitude variations of seismic waves. In applying the coda wave method, firstly, a type curve representing the average pattern of the individual coda decay curves for 0.75, 1.5, 3.0, 6.0, 12.0, and 24.0 Hz values was estimated. Secondly, lateral variation of coda Q and the attenuation coefficients for three main tectonic patterns are estimated. The shape of the type curve is controlled by the scattering and attenuation in the crustal volume sampled by the coda waves. The Q _o and η values vary from 30 to 180 and from 0.55 to 1.25, respectively for SW Anatolia. In SW Anatolia, coda Q–f relation is described by and δ = 0.008 km⁻¹. These results are expected to help in understanding the degree of tectonic complexity of the crust in SW Anatolia.