Simultaneous inversion of regional wave spectra for attenuation and seismic moment in Scandinavia

Frequency-dependent regional wave attenuation is estimated for continental paths to the NORESS array in Norway. Regional Lg and Pn spectra from 186 events at ranges between 200 and 1400 km and local magnitudes between 1.1 and 4.8 are inverted for both seismic moment and apparent attenuation. The Lg spectra were inverted between 1 and 7 Hz, and the Pn spectra were inverted between 1 and 15 Hz. The method uses both the spectral and spatial decay of observed signal amplitudes to separate source and path contributions. The assumptions include the geometric spreading rate and the source spectrum to be uniquely defined by its long-period level. Most events considered have local magnitudes less than 3.0, so the source corner frequencies are near or beyond the upper limit of the inverted bandwidth. The Q results, particularly for Lg, are therefore not very sensitive to the details of our source parameterization. The inversion parameters are source moment (for each event), a constant relating corner frequency and moment for the entire data set, and two parameters describing a power law frequency dependence of Q in the region. For fixed source and spreading assumptions the inversion defines clear trade-offs among model parameters. These trade-offs are resolved by adding the constraint that the separately derived source parameters for Lg and Pn are consistent. The “preferred” estimates for the apparent attenuation are QLg(f) = 560f0.26 and QPn(f) = 325f0.48. These Q values correspond to assumed geometric spreading rates of r−0.5 for Lg and r−1.3 for Pn. For fixed Lg spreading, the Pn spreading rate is constrained by requiring that the relative Lg amplitude for earthquakes and explosions of the same moment be consistent with well-supported results from previous empirical studies. The relationship between the inverted seismic moment values and local magnitude is generally consistent with values from near-field studies. Since magnitude does not enter the inversion, this result lends considerable support to the derived Q models. Whatever the physical interpretation of the results, they certainly provide an accurate parameterization of observed amplitude spectra in this region. This is valuable for representing wave propagation in the region, and it provides important data for assessing the event monitoring capabilities of small regional networks.

[1]  K. Aki.  Scaling law of seismic spectrum , 1967 .

[2]  C. Langston.  Aspects of Pn and Pg propagation at regional distances , 1982 .

[3]  H. Hasegawa.  Attenuation of Lg waves in the Canadian Shield , 1985 .

[4]  T. C. Bache,et al.  Regional Wave Attenuation and Seismic Moment from the Inversion of NORESS Spectra. , 1987 .

[5]  E. S. Husebye,et al.  Seismic array configuration optimization , 1983 .

[6]  R. Herrmann,et al.  Spatial attenuation of the Lg wave in the Central United States , 1983 .

[7]  O. Nuttli,et al.  Yield estimates of nevada test site explosions obtained from seismic Lg waves , 1986 .

[8]  L. Knopoff,et al.  Interpretation of Lg , 1973 .

[9]  J. Brune.  Tectonic stress and the spectra of seismic shear waves from earthquakes , 1970 .

[10]  Robert B. Herrmann,et al.  Lg attenuation and source studies using 1982 Miramichi data , 1987 .

[11]  Robert B. Herrmann,et al.  Regionalization of crustal coda Q in the continental United States , 1983 .

[12]  Paul W. Pomeroy,et al.  Test ban treaty verification with regional data—A review , 1982 .

[13]  Thomas C. Hanks,et al.  b values and ω−γ seismic source models: Implications for tectonic stress variations along active crustal fault zones and the estimation of high‐frequency strong ground motion , 1979 .

[14]  J. R. Clements.  Intrinsic Q and its frequency dependence , 1982 .

[15]  H. Bungum,et al.  The Meløy earthquake sequence, northern Norway: Source parameters and their scaling relations , 1982 .

[16]  Jack F. Evernden,et al.  An Evaluation of Seismic Decoupling and Underground Nuclear Test Monitoring Using High-Frequency Seismic Data (Paper 5R0913) , 1986 .

[17]  B. Kennett,et al.  Guided wave propagation in laterally varying media — II. Lg-waves in north-western Europe , 1984 .

[18]  J. E. Glynn,et al.  Numerical Recipes: The Art of Scientific Computing , 1989 .

[19]  K. Priestley,et al.  Measurement of frequency dependent Lg attenuation in the Great Basin , 1986 .

[20]  D. E. Willis.  Comparison of seismic waves generated by different types of source , 1963 .

[21]  L. Frazer,et al.  High‐frequency seismic attenuation of oceanic P and S waves in the western Pacific , 1987 .

[22]  S. R. Taylor,et al.  Attenuation and scattering of broadband P and S waves across North America , 1986 .

[23]  O. Nuttli,et al.  Seismic wave attenuation and magnitude relations for eastern North America , 1973 .

[24]  J. L. Stevens,et al.  The physical basis of mb:Ms and variable frequency magnitude methods for earthquake/explosion discrimination. Topical report , 1982 .

[25]  O. Nuttli,et al.  On the attenuation of Lg waves in western and central Asia and their use as a discriminant between earthquakes and explosions , 1981 .

[26]  H. Hasegawa.  Lg spectra of local earthquakes recorded by the Eastern Canada Telemetered Network and spectral scaling , 1983 .

[27]  G. Panza,et al.  Lg, Li and Rg from Rayleigh Modes , 1975 .

[28]  John R. Murphy,et al.  Seismic characteristics of underground nuclear detonations Part I. Seismic spectrum scaling , 1971 .

[29]  Michel Campillo,et al.  Frequency-dependent attenuation in the crust beneath Central France from Lg waves: Data analysis and numerical modeling , 1985 .

[30]  Gordon F. West,et al.  A novel technique for measuring Lg attenuation—results from Eastern Canada between 1 to 10 hz , 1987 .

[31]  Robin K. McGuire,et al.  The character of high-frequency strong ground motion , 1981 .

[32]  J. Orcutt,et al.  Synthetic Pn and Sn phases and the frequency dependence of Q of oceanic lithosphere , 1987 .

[33]  T. J. Bennett,et al.  A discrimination analysis of short-period regional seismic data recorded at Tonto Forest Observatory , 1982 .

[34]  T. C. Bache,et al.  Spectral characteristics of regional phases recorded at Noress , 1988 .

[35]  B. Mitchell,et al.  Crustal Q structure in the United States from multi-mode surface waves , 1981 .

[36]  S. Day,et al.  Simulation of short period Lg, expansion of three-dimensional source simulation capabilities and simulation of near-field ground motion from the 1971 San Fernando, California, earthquake. Final report 1 Oct 79-30 Nov 80 , 1981 .

[37]  Robert B. Herrmann,et al.  Modeling some empirical vertical component Lg relations , 1983 .

[38]  E. Chael.  Spectral scaling of earthquakes in the Miramichi region of New Brunswick , 1987 .

[39]  Robert B. Herrmann,et al.  Spectral Characteristics of the Lg wave Generated by Central United States Earthquakes , 1975 .

[40]  T. C. Bache,et al.  Q for teleseismic P waves from central Asia , 1985 .

[41]  Charles S. Mueller,et al.  Source parameters from locally recorded aftershocks of the 9 January 1982 Miramichi, New Brunswick, earthquake , 1985 .

[42]  I. Gupta,et al.  Attenuation of ground motion in the Eastern United States , 1987 .

[43]  H. Bungum,et al.  Processing of regional seismic events using data from small-aperture arrays , 1984 .