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.     

AVO investigations of shallow marine sediments

Amplitude‐variation‐with‐offset (AVO) analysis is based on the Zoeppritz equations, which enable the computation of reflection and transmission coefficients as a function of offset or angle of incidence. High‐frequency (up to 700 Hz) AVO studies, presented here, have been used to determine the physical properties of sediments in a shallow marine environment (20 m water depth). The properties that can be constrained are P‐ and S‐wave velocities, bulk density and acoustic attenuation. The use of higher frequencies requires special analysis including careful geometry and source and receiver directivity corrections. In the past, marine sediments have been modelled as elastic materials. However, viscoelastic models which include absorption are more realistic. At angles of incidence greater than 40°, AVO functions derived from viscoelastic models differ from those with purely elastic properties in the absence of a critical angle of incidence. The influence of S‐wave velocity on the reflection coefficient is small (especially for low S‐wave velocities encountered at the sea‐floor). Thus, it is difficult to extract the S‐wave parameter from AVO trends. On the other hand, P‐wave velocity and density show a considerably stronger effect. Attenuation (described by the quality factor Q) influences the reflection coefficient but could not be determined uniquely from the AVO functions. In order to measure the reflection coefficient in a seismogram, the amplitudes of the direct wave and the sea‐floor reflection in a common‐midpoint (CMP) gather are determined and corrected for spherical divergence as well as source and streamer directivity. At CMP locations showing the different AVO characteristics of a mud and a boulder clay, the sediment physical properties are determined by using a sequential‐quadratic‐programming (SQP) inversion technique. The inverted sediment physical properties for the mud are: P‐wave velocity α=1450±25 m/s, S‐wave velocity β=90±35 m/s, density ρ=1220±45 kg/m³, quality factor for P‐wave Q_P=15±200, quality factor for S‐wave Q_S=10±30. The inverted sediment physical properties for the boulder clay are: α=1620±45 m/s,β=360±200 m/s,ρ=1380±85 kg/m³,Q_P=790±660,Q_S=25±10.     

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.     

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.     

Estimation of frequency dependent coda wave attenuation structure at the vicinity of Cairo Metropolitan Area

Estimation of seismic wave attenuation in the shallow crust in terms of coda wave Q structure previously investigated in the vicinity of Cairo Metropolitan Area was improved using seismograms of local earthquakes recorded by the Egyptian National Seismic Network. The seismic wave attenuation was measured from the time decay of coda wave amplitudes on narrow bandpass filtered seismograms based on the single scattering theory. The frequency bands of interest are from 1.5 to 18 Hz. In general, the values obtained for various events recorded at El-Fayoum and Wadi Hagul stations are very similar for all frequency bands. A regional attenuation law Q _c = 85.66 f ^0.79 was obtained.     

Estimating permeability from field measurements of seismic attenuation in fractured chalk

Broadband (100–4000 Hz) cross‐hole seismic data have been acquired at a borehole test site where extensive hydrological investigations have previously been performed, including in situ estimates of permeability. The rock type is homogeneous chalk and fractures and bedding planes have been identified from well logs. High values of seismic attenuation, Q= 22 ≤ 27 ≤ 33, were observed over a 10 m depth interval where fracture permeability values of 20–50 darcy had been recorded. An attempt has been made to separate the attenuation due to scattering and intrinsic mechanisms. The estimated values of intrinsic attenuation, Q= 31 ≤ 43 ≤ 71, have been reproduced using a number of current theories of seismic‐wave propagation and fluid‐flow‐induced seismic attenuation in cracked and fractured media. A model that considers wavelength‐scale pressure gradients is the preferred attenuation mechanism. Model parameters were obtained from the hydro‐geological and seismic data. However, we conclude that it is not possible to use seismic Q to measure rock permeability remotely, principally because of the inherent uncertainties arising from model parameterisations.     

Anelasticity of the crust and upper mantle beneath north and central India from the inversion of observed Love and Rayleigh wave attenuation data

The fundamental mode Love and Rayleigh waves generated by earthquakes occurring in Kashmir, Nepal Himalaya, northeast India and Burma and recorded at Hyderabad, New Delhi and Kodaikanal seismic stations are analysed. Love and Rayleigh wave attenuation coefficients are obtained at time periods of 15–100 seconds, using the spectral amplitude of these waves for 23 different paths along northern (across Burma to New Delhi) and central (across Kashmir, Nepal Himalaya and northeast India to Hyderabad and Kodaikanal) India. Love wave attenuation coefficients are found to vary from 0.0003 to 0.0022 km^–1 for northern India and 0.00003 km^–1 to 0.00016 km^–1 for central India. Similarly, Rayleigh wave attenuation coefficients vary from 0.0002 km^–1 to 0.0016 km^–1 for northern India and 0.00001 km^–1 to 0.0009 km^–1 for central India. Backus and Gilbert inversion theory is applied to these surface wave attenuation data to obtainQ ^–1 models for the crust and uppermost mantle beneath northern and central India. Inversion of Love and Rayleigh wave attenuation data shows a highly attenuating zone centred at a depth of 20–80 km with lowQ for northern India. Similarly, inversion of Love and Rayleigh wave attenuation data shows a high attenuation zone below a depth of 100 km. The inferred lowQ value at mid-crustal depth (high attenuating zone) in the model for northern India can be by underthrusting of the Indian plate beneath the Eurasian plate which has caused a low velocity zone at this shallow depth. The gradual increase ofQ ^–1 from shallow to deeper depth shows that the lithosphere-asthenosphere boundary is not sharply defined beneath central India, but rather it represents a gradual transformation, which starts beneath the uppermost mantle. The lithospheric thickness is 100 km beneath central India and below that the asthenosphere shows higher attenuation, a factor of about two greater than that in the lithosphere. The very lowQ can be explained by changes in the chemical constitution taking place in the uppermost mantle.     

Study of theQ-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.     

Empirical model for estimating Fourier amplitude spectra of ground acceleration in Taiwan region

A collection of ground‐motion recordings (1070 acceleration records) of moderate (5.1⩽M_L⩽6.5) earthquakes obtained during the execution of the Taiwan Strong Motion Instrumentation Program (TSMIP) since 1991 was used to study source scaling model and attenuation relations for a wide range of earthquake magnitudes and distances and to verify the models developed recently for the Taiwan region. The results of the analysis reveal that the acceleration spectra of the most significant part of the records, starting from S‐wave arrival, can be modelled accurately using the Brune's ω‐squared source model with magnitude‐dependent stress parameter Δσ, that should be determined using the recently proposed regional relationships between magnitude (M_L) and seismic moment (M₀) and between M₀ and Δσ. The anelastic attenuation Q of spectral amplitudes with distance may be described as Q=225 ƒ^1.1 both for deep (depth more than 35 km) and shallow earthquakes. The source scaling and attenuation models allow a satisfactory prediction of the peak ground acceleration for magnitudes 5.1⩽M⩽6.5 and distances up to about 200 km in the Taiwan region, and may be useful for seismic hazard assessment. Copyright © 2000 John Wiley & Sons, Ltd.     

Borehole study of compressional and shear attenuation of basalt flows penetrated by the Brugdan and William wells on the Faroes shelf

We investigated the seismic attenuation of compressional (P‐) and converted shear (S‐) waves through stacked basalt flows using short‐offset vertical seismic profile (VSP) recordings from the Brugdan (6104/21–1) and William (6005/13–1A) wells in the Faroe‐Shetland Trough. The seismic quality factors (Q) were evaluated with the classical spectral ratio method and a root‐mean‐square time‐domain amplitude technique. We found the latter method showed more robust results when analysing signals within the basalt sequence. For the Brugdan well we calculated effective Q estimates of 22–26 and 13–17 for P‐ and S‐waves, respectively, and 25–33 for P‐waves in the William well. An effective Q_S/Q_P ratio of 0.50–0.77 was found from a depth interval in the basalt flow sequence where we expect fully saturated rocks. P‐wave quality factor estimates are consistent with results from other VSP experiments in the North Atlantic Margin, while the S‐wave quality factor is one of the first estimates from a stacked basalt formation using VSP data. Synthetic modelling demonstrates that seismic attenuation for P‐ and S‐waves in the stacked basalt flow sequence is mainly caused by one‐dimensional scattering, while intrinsic absorption is small.     

Determination of moment magnitude from acceleration spectrum

The theoretical acceleration spectrum of observation site has been obtained from source acceleration spectrum derived from scaling law, using attenuation modelQ=Q _ν f ^η . A comparison of a set of theoretical acceleration spectra with observation spectra has been made, and we have obtained the attenuation model for observation site and seismic moment magnitude. We obtain thatQ _o=300,η=0.25 for Wuqia area, Xinjing Zizhizhou, and seismic moment magnitudes of 18 greater aftershocks of Wuqia earthquake occurred in 1985. In order to obtain seismic moment magnitued conveniently, three functional tables of acceleration spectra at 1Hz as the distances for variousQ value have been made. The seismic moment magnitude can be quickly measured from acceleration spectrum at 1Hz according to these tables (epicenter has to be known). The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, 435–445, 1992.     

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

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

Multiwell study of seismic attenuation at the CO2CRC Otway project geosequestration site: Comparison of amplitude decay,centroid frequency shift and 1D waveform inversion methods

At the CO2CRC Otway geosequestration site, the abundance of borehole seismic and logging data provides a unique opportunity to compare techniques of Q (measure of attenuation) estimation and validate their reliability. Specifically, we test conventional time-domain amplitude decay and spectral-domain centroid frequency shift methods versus the 1D waveform inversion constrained by well logs on a set of zero-offset vertical seismic profiles. The amplitude decay and centroid frequency shift methods of Q estimation assume that a seismic pulse propagates in a homogeneous medium and ignore the interference of the propagating wave with short-period multiples. The waveform inversion explicitly models multiple scattering and interference on a stack of thin layers using high-resolution data from sonic and density logs. This allows for stable Q estimation in small depth windows (in this study, 150 m), and separation of the frequency-dependent layer-induced scattering from intrinsic absorption. Besides, the inversion takes into account band-limited nature of seismic data, and thus, it is less dependent on the operating frequency bandwidth than on the other methods. However, all considered methods of Q estimation are unreliable in the intervals where subsurface significantly deviates from 1D geometry. At the Otway site, the attenuation estimates are distorted by sub-vertical faults close to the boreholes. Analysis of repeated vertical seismic profiles reveals that 15 kt injection of the CO₂-rich fluid into a thin saline aquifer at 1.5 km depth does not induce detectable absorption of P-waves at generated frequencies 5–150 Hz, most likely because the CO₂ plume in the monitoring well is thin, <15 m. At the Otway research site, strong attenuation Q ≈ 30–50 is observed only in shaly formations (Skull Creek Mudstone, Belfast Mudstone). Layer-induced scattering attenuation is negligible except for a few intervals, namely 500–650 m from the surface, and near the injection interval, at around 1400–1550 m, where Q_scat ≈ 50–65.     

Universality of the Seismic Moment-frequency Relation

—We analyze the seismic moment-frequency relation in various depth ranges and for different seismic regions, using Flinn-Engdahl's regionalization of global seismicity. Three earthquake lists of centroid-moment tensor data have been used the Harvard catalog, the USGS catalog, and the Huang et al. (1997) catalog of deep earthquakes. The results confirm the universality of the β-values and the maximum moment for shallow earthquakes in continental regions, as well as at and near continental boundaries. Moreover, we show that although fluctuations in earthquake size distribution increase with depth, the β-values for earthquakes in the depth range of 0–500 km exhibit no statistically significant regional variations. The regional variations are significant only for deep events near the 660 km boundary. For declustered shallow earthquake catalogs and deeper events, we show that the worldwide β-values have the same value of 0.60 ± 0.02. This finding suggests that the β-value is a universal constant. We investigate the statistical correlations between the numbers of seismic events in different depth ranges and the correlation of the tectonic deformation rate and seismic activity (the number of earthquakes above a certain threshold level per year). The high level of these correlations suggests that seismic activity indicates tectonic deformation rate in subduction zones. Combined with the universality of the β-value, this finding implies little if any variation in maximum earthquake seismic moment among various subduction zones. If we assume that earthquakes of maximum size are similar in different depth ranges and the seismic efficiency coefficient, χ, is close to 100% for shallow seismicity, then we can estimate χ for deeper earthquakes for intermediate earthquakes χ≈ 5%, and χ≈ 1% for deep events. These results may lead to new theoretical understanding of the earthquake process and better estimates of seismic hazard.     

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.     

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.     

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.     

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.     

Viscoacoustic least-squares reverse-time migration of different-order free-surface multiples

Multiples have longer propagation paths and smaller reflection angles than primaries for the same source–receiver combination, so they cover a larger illumination area. Therefore, multiples can be used to image shadow zones of primaries. Least-squares reverse-time migration of multiples can produce high-quality images with fewer artefacts, high resolution and balanced amplitudes. However, viscoelasticity exists widely in the earth, especially in the deep-sea environment, and the influence of Q attenuation on multiples is much more serious than primaries due to multiples have longer paths. To compensate for Q attenuation of multiples, Q-compensated least-squares reverse-time migration of different-order multiples is proposed by deriving viscoacoustic Born modelling operators, adjoint operators and demigration operators for different-order multiples. Based on inversion theory, this method compensates for Q attenuation along all the propagation paths of multiples. Examples of a simple four-layer model, a modified attenuating Sigsbee2B model and a field data set suggest that the proposed method can produce better imaging results than Q-compensated least-squares reverse-time migration of primaries and regular least-squares reverse-time migration of multiples.