Temporal variation of codaQ during Hualien earthquake of 1986 in eastern Taiwan

The May 20, 1986, Hualien earthquake sequence occurred on the northeastern coast of Taiwan. TheM=6.1 (GS mb) mainshock was followed by a large number of closely-clustered aftershocks with the largest being anM=5.5 event. One seismic station, TWD of the Taiwan Telemetered Seismographic Network, is located in the surface projection of the source region and provides excellent recordings of the entire earthquake sequence. These recordings, plus events occurring in the same source area preceding the mainshock, offer a unique opportunity to study the spatial and temporal variations of codaQ in a region of active subduction. A simple technique is devised that uses the envelope of the coda waveform to enable a quick determination of the codaQ from drum records. For recordings with a peak power at about 8 Hz, the following findings have been obtained: 1. The ambient codaQ near an active subduction region was as low as 145; 2. There was no significant decrease in codaQ within the period beginning one year and four months prior to the mainshock; 3. There was a significant drop of codaQ immediately after the mainshock; this drop lasted approximately two days before returning to the ambient level; 4. CodaQ values varied with focal depth.     

Temporal and spatial variation on codaQ −1 associated with the North Palm Springs earthquake of July 8, 1986

A temporal and spatial change of codaQ ^–1 associated with the occurrence of the North Palm Springs earthquake of July 8, 1986 was studied by using 242 small local earthquakes in the vicinity of the mainshock. We found that the codaQ ^–1 of earthquakes which occurred before the mainshock was significantly higher than that of the aftershocks in the mainshock area while the codaQ ^–1 for the surrounding area remained almost constant throughout 1986. CodaQ ^–1 was determined separately for the lapse time windows of 10 to 20 sec. and 15 to 40 sec. for the period from 1981 to 1987. The result for the time window 10 to 20 sec. showed a peak in codaQ ^–1 before the time of mainshock at all frequencies. The peak appeared earlier at lower frequencies. There was no significant change in codaQ ^–1 for the time window 15 to 40 sec., probably because the change was restricted to a small area.     

Lapse time and frequency-dependent attenuation characteristics of coda waves in the Northwestern Himalayas

We analyzed the local earthquakes waveform recorded on a broadband seismic network in the northwestern Himalayan Region to compute lapse time and frequency dependence of coda Q (Q _c). The observed Q _c values increase with increasing lapse time at all frequency bands. The increase in Q _c values with lapse time is attributed to an increase in Q _c with depth. This implies that attenuation decreases with increasing depth. The approximate radius of medium contributing to coda generation varies from 55 to 130 km. By comparing the Q _c values with those from other regions of the world, we find that they are similar to those obtained from tectonically active regions. The estimated Q _c values show a frequency-dependent relationship, Q _c = Q ₀ f ⁿ , where Q ₀ is Q _c at 1 Hz and n represents degree of frequency dependence. They represent the level of heterogeneity and tectonic activity in an area. Our results show that northwest Himalayas are highly heterogeneous and tectonically very active. Q ₀ increases from 113 ± 7 to 243 ± 10 and n decreases from 1.01 ± 0.05 to 0.85 ± 0.03 when lapse time increases from 30 to 70 s. As larger time window sees the effect of deeper part of the Earth, it is concluded that Q ₀ increases and n decreases with increasing depth; i.e., heterogeneity decreases with depth in the study area.     

Coda-Q before and after the 1986 Andreanof Islands earthquake

Coda-Q values before and after the May 7, 1986 Andreanof Islands earthquake in the central Aleutian subduction zone were calculated from microearthquake seismograms on the basis of Aki's single backscattering model. Digital waveforms from 329 microeathquakes in 5.4 years before the mainshock and 40 aftershocks recorded by 11 stations of the Central Aleutians Seismic Network were utilized. Before the mainshock, coda-Q in the rupture zone west of the mainshock epicenter was 15% higher than theQ outside the rupture zone to the west. A lowQ in the region east of the mainshock epicenter is inferred, but lacks solid evidence. The highQ region accords with the part of the rupture zone where most seismic moment was released. During approximately two months after the mainshock, coda-Q in the rupture zone decreased about 10%. PrecursoryQ changes were not found in the 5.4 years before the mainshock. It is deduced that a coda-Q precursor, if it existed, was of small magnitude (<10%), or its duration was either shorter than one year or longer than four years.     

Estimation of Coda Wave Attenuation for the National Capital Region, Delhi, India Using Local Earthquakes

Attenuation of seismic waves is very essential for the study of earthquake source parameters and also for ground-motion simulations, and this is important for the seismic hazard estimation of a region. The digital data acquired by 16 short-period seismic stations of the Delhi Telemetric Network for 55 earthquakes of magnitude 1.5 to 4.2, which occurred within an epicentral distance of 100 km in an area around Delhi, have been used to estimate the coda attenuation Q_c. Using the Single Backscattering Model, the seismograms have been analyzed at 10 central frequencies. The frequency dependence average attenuation relationship Q_c = 142f ^1.04 has been attained. Four Lapse-Time windows from 20 to 50 seconds duration with a difference of 10 seconds have been analyzed to study the lapse time dependence of Q_c. The Q_c values show that frequency dependence (exponent n) remains similar at all the lapse time window lengths. While the change in Q₀ values is significant, change in Q₀ with larger lapsetime reflects the rate of homogeneity at the depth. The variation of Q_c indicates a definitive trend from west to east in accordance with the geology of the region.     

Coda-Q before and after the 1983 Misasa earthquake of M 6.2, Tottori Prefecture,Japan

To detect any temporal change in coda-Q associated with an earthquake of M 6.2, digitized data from analog magnetic records of microearthquakes for a small confined region just around its hypocenter have been processed, with special attention to measuring errors. The M 6.2 earthquake took place at the Misasa town, Tottori Prefecture, in southwest Japan with preceding anomalous seismicity changes. TheQ values were elevated 20% around three years before the M 6.2 event for frequencies of 5 to 20 Hz, and tended to decrease around two years before. Data with high quality show undulated temporal variation with a period of 5–6 years before the earthquake, which is correlated with the regional seismic activity. There is a possibility that the observedQ change stands for precursory effect of the medium-scale earthquake. The aftershock sequence also shows an undulated temporal variation ofQ with a period of 150 days for around 10Hz, correlating also with the seismic activity. The fracturing processes by numerous microearthquakes may be responsible for the attenuation property of randomly scattered seismic waves.     

Estimating seismic moments and Lg Q using Lg spectra

A multi-event and multi-station inverse method is presented in the paper to simultaneously estimate the seismic moments (M ₀) and source corner frequencies (f _c) of several Jiashi (Xinjiang, China) earthquakes, as well as the apparent Lg Q models for the paths from Jiashi to eight seismic stations (WMQ, AAK, TLG, MAKZ, KUR, VOS, ZRN and CHK) in Central Asia. The resultant seismic moments correlate well with the M ₀ values obtained by Harvard University using the centroid moment tensor (CMT) inversion and the surface-wave magnitudes as well. After the correction by a typical value of average radiation coefficient for regional SV waves, the M ₀ values from Lg spectral inversion are still close to the corresponding values obtained from CMT inversion. The obtained apparent Q _0Lg values (Lg Q at 1 Hz) are consistent with the tectonic features of corresponding propagation paths. The Q _0Lg values are 351±87, 349±86 and 300±27 for the paths from Jiashi to AAK, TLG and MAKZ, respectively. They are smaller than Q _0Lg values for the paths to KUR, VOS, ZRN and CHK, which are 553±72, 569±58, 550±57 and 603±65, respectively. These results agree with the condition that the paths to AAK, TLG and MAKZ mainly propagate through the mountainous Tianshan area where relatively strong seismic activities and large variations of topography are exhibited, while the paths to KUR, VOS, ZRN and CHK mainly propagate through the stable area of Kazak platform. The Q _0Lg value for the path to WMQ is 462±56. This is also in agreement with the condition that the path to WMQ is basically along the border area between Tianshan Mountain and Tarim Basin, and along this path the variations of topography and crustal thickness are moderate in comparison with that along the path to MAKZ. Foundation item: Foundation of Verification Researches for Army Control Technology (413290102).     

Distribution of<Emphasis Type="Italic">L</Emphasis><Subscript>g</Subscript> coda<Emphasis Type="Italic">Q</Emphasis> in the Chinese continent and its adjacent regions

785 traces of vertical components from shallow earthquakes recorded by 10 CDSN (Chinese Digital Seismographic Network) stations and 5 GSN (Global Seismographic Network) stations were collected to study the attenuation characteristics ofL _g coda in the Chinese continent and its adjacent regions. The records were processed first using the stack spectral ratio method to obtain the average values ofQ ₀ (Q at 1Hz) and η, the frequency dependence, ofL _g coda in the ellipses corresponding to the paths. The back-projection technique was then employed to obtain the tomographic maps ofQ ₀ and η values, and the distribution of their errors. Results indicate that in the studied areaQ ₀ varies between 200 and 500. The lowest value ofQ ₀ exists in the Yun-nan-Tibetan region, while the highest value ofQ ₀ occurs in the southern edge of Siberian platform. η varies between 0.3 and 0.8. For most part of the studied area η varies inversely withQ ₀.     

Estimation of coda wave attenuation in East Central Iran

The attenuation of coda waves, Q _c , has been estimated in Zarand, Jiroft, and Bam regions of east central Iran using a single back-scattering model of S-coda envelopes. For this purpose, the recordings of 97 earthquakes by three seismic networks and a local strong ground motion network have been used. In this research, the frequency-dependent Q _c values are estimated at central frequencies of 1.5, 3, 6, 8, 12, 16, and 24 Hz using different lapse time windows from 20 to 60 s. The frequency-dependent relationships obtained are for Zarand, for Jiroft, and for Bam region. From the strong ground motion data, we obtain the relation . The Q _c frequency-dependent relationship for the entire region of east central Iran from all data (both seismograms and accelerograms) is . The average Q _c values estimated and their frequency dependent relationships correlate well with a highly heterogeneous and highly tectonically active region. Results also show that the attenuation is higher in Bam region compared to Zarand and Jiroft regions.     

Three-dimensional Q structure in Jiashi earthquake region of Xinjiang

3-D S-waveQ structure in Jiashi earthquake region is inverted based on the attenuation of seismic waves recorded from earthquakes in this region in 1998 by the Research Center of Exploration Geophysics (RCEG), CSB, and a rough configuration of deep crustal faults in the earthquake region is presented. First, amplitude spectra of S-waves are extracted from 450 carefully-chosen earthquake records, called observed amplitude spectra. Then, after instrumental and site effect correction, theoretical amplitude spectra are made to fit observed amplitude spectra with nonlinear damped least-squares method to get the observed travel time overQ, provided that earthquake sources conform to Brune’s disk dislocation model. Finally, by 3-D ray tracing method, theoretical travel time overQ is made to fit observed travel time overQ with nonlinear damped least-squares method. In the course of fitting, the velocity model, which is obtained by 3-D travel time tomography, remains unchanged, while onlyQ model is modified. When fitting came to the given accuracy, the ultimateQ model is obtained. The result shows that an NE-trending lowQ zone exists at the depths of 10–18 km, and an NW-trending lowQ zone exists at the depths of 12–18 km. These roughly coincide with the NE-trending and the NW-trending low velocity zones revealed by other scientists. The difference is that the lowQ zones have a wider range than the low velocity zones. Foundation item: Joint Seismological Science Foundation of China (957-07-414) and State Key Basic Research Development and Programming Project (95-13-02-02). Contribution No. RCEG200105, Research Center of Exploration Geophysics, China Seismological Bureau.     

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.     

Attenuation of Coda Waves at the Changbaishan Tianchi Volcanic Area in Northeast China

21 earthquakes recorded by a temporary seismic network in the Changbaishan Tianchi volcanic area in Northeast China operated during the summer of 2002 and 2003 were analyzed to estimate the S coda attenuation. The attenuation quality factor Q_c was estimated using the single scattering attenuation model of Sato (1977) in the frequency band from 4 to 24 Hz. All the events studied in this paper occurred at depths from 2 to 6 km with ML of 1.4–2.8. The epicentral distances are less than 25 km. For all events which occurred near the Tianchi Lake (caldera), the Q_c patterns obtained at the stations near the lake are similar, and the Q_c values are relatively small. At the stations located about 15 km east of the Tianchi Lake, however, the average Q_c is significantly higher. For an event which occurred 25km from the lake to the west, Q_c patterns derived at the stations near the lake are quite similar to the above mentioned Q_c for stations located in the east. Further study shows that Q_c value in the north and central areas of the volcano is relatively lower than that in the surrounding area. Compared to other volcanic areas in the world, the average Q_c of the Changbaishan Tianchi volcanic area is obviously lower. The deep seismic sounding and teleseismic receiver function studies indicated more than one lower velocity layer in the crust. The MT studies suggested the presence of high conductive bodies beneath the area. We interpret the strong attenuation of coda waves near the Changbaishan Tianchi volcano as being possibly related to high temperature medium caused by shallow magma chambers.     

Study of the Q -value in and around Minle Basin

The digital seimograph network set up by China and France in Zhangye of China had been operated in 1988. The Zhangye network is situated in the middle segment of Hexi corridor and Qilina mountain, which was regarded as a monitoring earthquake area. Using the records of Zhangye digital network theQ-values in and around Minle basin have been measured. The results of this study showed that theQp-values range from 500 to 780, andQs-values range from 230 to 460. TheQ-values of inside of Minle basin are higher than that around the basin. The greater parts of moderate and strong earthquakes occurred along tectonic belts around the Minle basin. Moreover, TheQ-values increased with the depth of penctration of wave ray. The attenuation of S wave is stronger than P wave in shallow layer of crust. Some problems ofQ-value change versus time before and after Sunan eathquake (M _s=5.7) have been also disscussed. These results can be applied to study and to monitor seismic danger of the Minle monitoring area. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 296–302, 1993.     

Attenuation and source parameters at Deception Island (South Shetland Islands,Antarctica)

Deception Island is the most active volcano of the South Shetland Islands-Antarctic Peninsula group, experiencing eruptions in 1967, 1969 and 1970. Local attenuation derived from coda analysis and source parameters derived from Brune's model, for well located seismic events, have been studied in order to complement the available geophysical information. Results show abnormally lowQ ₀ values and an abnormally high frequency dependence, as well as large dispersion. These factors are strongly dependent on the path travelled by the seismic wave. Retrieved values of the source parameters (stress drop, seismic moment and source radius), are again abnormally low compared to world-wide average values, for example, those obtained for the Oroville, California aftershock series between June and September, 1976. These results are consistent with some aspects of the geology of Deception Island, such as the very high degree of fracturing and faulting, and the existence of a strong hydrothermal alteration affecting most of the subaerial volcanic rocks. Moreover, the pattern defined for the lateral variations ofQ ₀ shows minimum values in the inner bay of the island, close to the most recent eruption vents. A large reduction in spectral amplitudes over a particular frequency range occurs in several observations, corresponding to the path crossing the zone of highest attenuation. This observation suggests the existence of a hot magmatic intrusion produced during the most recent eruption, and coincides with the superficial low density mass distribution obtained from the gravimetric model and the long wave magnetic field component obtained from magnetic surveys. The width of this intrusion is estimated to be about 200 m, in agreement with the previous results obtained analyzing residuals of the location of seismic events.     

Distribution of Lg coda Q in Xinjiang and its adjacent regions

In this study, we collected 1 156 broadband vertical components records at 22 digital seismic stations in Xinjiang region, ürümqi station, and 7 stations in the adjacent regions during the period of 1999–2003. The records were firstly processed by the stacked spectral ratio method to obtain Q ₀ (Q at 1 Hz) and the frequency correlation factor η corresponding to each path. Based on the results, the distribution images of Q ₀ and η in 1°×1° grids for Xinjiang region were gained by the back-projection technique. The results indicate that Q ₀ is high (300–450) in the Tarim platform and marginal Siberian platform, while Q ₀ is low (150–250) in the southern regions as west Kunlun fold system and Songpan-Ganzi fold system. In the northern regions as Junggar fold system and Tianshan fold system, Q ₀ is also low (250–300) and η varies between 0.5 and 0.9. Foundation item: National Natural Science Foundation of China (49974012) and Joint Seismological Science Foundation of China (604004).     

Time domain solution for multiple scattering and the coda envelopes

This article summarizes work on multiple scattering based on models of media with randomly distributed scatterers. The scatterers are isotropic and statistically uniform. Measuring distance in terms of mean-free pathL _s and time in terms of the mean-free timesL _s/V, whereV is the velocity of scattered waves, we have more convenient dimensionless distance and time. It can be shown that after the dimensionless time equals 0.65 energy contributed from multiple scattering becomes predominant. Thus the later coda reflects the effect of multiple scattering rather than single scattering. Treating the seismic record, including starting and tail parts, as a whole, the diffusion theory predicts that at a dense distribution of scatterers and a small distance between source and receiver, codas reflect mainly intrinsicQ _i. Of course, this conclusion is coincident with the presumption of the diffusion theory,Q _s>Q _i. However, from a new integral equation of multiple scattering, which deals with the scattered waves and primary waves separately, the conclusion is similar but clearer. This article quotes the new expression for coda energy in two-dimensional space. It shows that if the receiver is close to the source, the coda decay reflects only intrinsicQ _i, then as the distance increases, effects of scatteringQ _s, are involved in the decay feature. The theoretical plots of coda decay show that it seems in most cases in the earthQ _i should not be smaller than one tenth ofQ _s.Project Sponsored by the Joint Earthquake Science Foundation of China.     

Attenuation of coda waves and Q c value beneath the Chengdu telemetered seismic network

Based on the single scattering model of coda power spectrum analysis, digital waveform data of 50 events recorded by the real-time processing system of the Chengdu telemetry network are analyzed to estimate the Q _c values of earth medium beneath the Chengdu telemetry network for several specified frequencies. It is found that the Q _c shows the frequency dependency in the form of Q _c = Q ₀ f ⁿ in the range of 1.0 to 20.0Hz. Estimated Q ₀ ranges from 60.83 to 178.05, and n is found to be 0.713 to 1.159. The average value of Q ₀ and n are 117 and 0.978 respectively. This result indicates the strong frequency dependency of the attenuation of coda waves beneath the Chengdu telemetry network. Comparing with the results obtained in other regions of the world, it is found that Q ₀ ⁻¹ value and its change with frequency are similar to those in regions with strong tectonic activity. This subject is supported by the Ministry of Personnel, China for partly sponsoring.     

Inferred attenuation from site effect-correctedT phases recorded on the island of Hawaii

Coda site amplification factors are used to eliminate the site effect from records of three circum-PacificT phases recorded by the Hawaiian Volcano Observatory (HVO) network on the island of Hawaii. ObservedT-phase amplitudes across the island generally decrease with increasing distance from the conversion point where acoustic waves in the SOFAR channel become seismic waves propagating through the crust. However, the decay of the observedT-phase signal across the island is not regular in regions of dense station coverage, in particular, the Kilauea caldera region. We divide the maximum observedT-phase amplitudes at a given station by the coda site amplification factor obtained for the same station and frequency band (3.0Hz); the distribution of these amplitudes reveals a smooth pattern over the entire island. The distance over which the site effect-correctedT-phase amplitude decreases by one-half, combined with the apparent velocity of propagation ofT phases across the island, allows for an approximate determination of near-surfaceQ over much of the island of Hawaii. We found a region of lowQ in the Kilauea summit area (Q≈30) and east rift zone (Q≈60) with considerably higherQ in the Kaoiki and northern portions of the island (Q≈150 to 200). The lowQ values obtained in the Kilauea summit region and east rift zone are significantly lower than estimates of codaQ in the same region, suggestingT phases may be sampling the earth's near-surface properties.     

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.     

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.