Bulletin of the Seismological Society of America
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Bulletin of the Seismological Society of America; February 2006; v. 96; no. 1; p. 348-354; DOI: 10.1785/0120040086
© 2006 Seismological Society of America
This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (3)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Aleqabi, G. I.
Right arrow Articles by Wysession, M. E.
Right arrow Search for Related Content
GeoRef
Right arrow GeoRef Citation

Short Note

QLg Distribution in the Basin and Range Province of the Western United States

Ghassan I. Aleqabi1 and Michael E. Wysession1

1 Department of Earth and Planetary Sciences
Washington University
St. Louis, Missouri 63130-4899


    Abstract
 Top
 Abstract
 Introduction
 Data
 QLg Tomography
 Results
 Discussion and Conclusions
 Acknowledgments
 References
 
Q values for the direct Lg phase (Formula ) were estimated for event-station paths that lie entirely within the western United States, particularly in the southern Great Basin and surrounding areas. The Qo (Formula at 1 Hz) values were estimated by fitting synthetic spectra to observed Lg spectra using a genetic algorithm technique. We have created a tomographic image of the variations of Qo for a part of the southwestern United States. The image shows that Qo varies between about 234 and 312 at 1 Hz, with an average of 267. The lowest Qo occurs in the northwest part of the Basin and Range Province, where extensional deformation has occurred since the Mesozoic. Qo values start to increase toward the Colorado Plateau to the east and continue to gradually increase northward and decrease southward. We also find relatively high frequency dependence to the Formula values, with a 1D mean of {eta} = 0.57.


    Introduction
 Top
 Abstract
 Introduction
 Data
 QLg Tomography
 Results
 Discussion and Conclusions
 Acknowledgments
 References
 
We present the results of a study of crustal attenuation in the southwestern United States using the Lg waves from regional earthquakes and nuclear tests (Fig. 1). The region of study encompasses part of the Basin and Range Province, as well as the southeastern part of the Sierra-Nevadas. The diversity of tectonic regimes in the western United States causes variations in the propagation and attenuation of seismic waves. Mitchell (1995) showed that the time elapsed since the most recent significant tectonic activity in a region is important in determining the attenuation properties there. Areas that have undergone more recent tectonism tend to be more attenuating. In addition, factors that increase attenuation of regional seismic phases include the presence of a thick cover of low Q sediments and velocity gradients (rather than a sharp interface) at the crust–mantle boundary (Mitchell and Cong, 1998). Attenuation is very sensitive to increases in temperature, and so attenuation tomography is well suited for identifying thermal anomalies.


Figure 001
View larger version (38K):
[in this window]
[in a new window]
 
Figure 1. Ray paths used in the tomographic inversion. Earthquakes are represented by diamonds, NTS explosions by circles, and Lawrence Livermore National Laboratory stations by triangles. Numbers indicate each event identification number (table 3 in Al-Eqabi et al., 2001).

 

Studies of Lg propagation from earthquake and explosion sources in the Basin and Range Province have reported regional variations in Formula values. In a prior study we calculated the moment and corner frequency for 40 seismic sources (Al-Eqabi et al., 2001). In this article we also calculate Qo (Formula at 1 Hz) and {eta} (the frequency dependence) and relate the crustal anelastic properties to the tectonic evolution and geology of the Basin and Range Province. We further apply the back-projection inversion technique to obtain a tomographic map of Qo variations in the Basin and Range Province and its surroundings.


    Data
 Top
 Abstract
 Introduction
 Data
 QLg Tomography
 Results
 Discussion and Conclusions
 Acknowledgments
 References
 
Our dataset consists of 111 vertical-component seismograms from 40 regional seismic events between 1988 and 1994 (table 3 in Al-Eqabi et al., 2001). Fifteen of these events were Nevada Test Site (NTS) nuclear explosions, and 25 were western U.S. earthquakes. To process the data, Lg spectra were isolated in a manner similar to that used by Chael (1987) and Atkinson and Mereu (1992). The individual source-station Qo values were obtained by applying the Lg spectral fit by a genetic algorithm (LGSFGA) technique (Al-Eqabi et al., 2001). The LGSFGA method was developed to analyze the Lg amplitude spectrum in terms of four independent variables: the seismic moment, corner frequency of the seismic source, average path attenuation at 1 Hz (Qo), and attenuation frequency dependence ({eta}). A genetic algorithm was used to efficiently search the parameter space and find optimum combinations of the four parameters to produce a calculated spectrum that best fit the observed Lg spectrum. The resulting Qo and {eta} values along with event information are given in Table 1. While the genetic algorithm is efficient in finding solutions that provide the least misfit to the observed spectra, it must be remembered that noise in the signals and effects from structural heterogeneities can result in a trade-off between the attenuation values and the corner frequency. There are other factors, such as geometric spreading, that can also provide a source of contamination for the individual path attenuation values. However, a consistent technique has been applied across all of the data, and so contaminating effects have been addressed as evenly as possible, but they are always a concern.


View this table:
[in this window]
[in a new window]
 
Table 1 Earthquakes and Explosions in Western United States with Results of LGSFGA Technique
 

Figure 1 shows the ray-path distribution within the area of study. Data coverage across the Basin and Range Province and its surroundings is not uniform because the earthquakes are mostly located in the southern Basin and Range and California coastal region, and the explosions are clustered inside the NTS. The coverage is best in the southern Basin and Range, while the northeastern part of the Basin and Range is poorly sampled.


    QLg Tomography
 Top
 Abstract
 Introduction
 Data
 QLg Tomography
 Results
 Discussion and Conclusions
 Acknowledgments
 References
 
Utilizing the assembled Qo values for each event-station pair (Table 1), we used the back-projection tomographic method (Xie and Mitchell, 1990a; Cong and Mitchell, 1998; L. Cong, personal comm., 2000) to create a tomographic image of the distribution of Qo over the region of study. This technique requires that the area of study be divided into a grid of rectangular cells, so we divided the area spanning latitudes 34–41° N and longitudes 113–120° W into 49 1° x cells (Nc). Then we calculated Qn, the Qo value calculated for the nth time series n = 1,2,3,... , Nd, for each cell m according to


Formula 001

(1)
where lmn is the path length in each cell, Qm is the Qo for each cell, {epsilon}n is a term signifying the error in measurement and modeling of Formula , Nd is the total number of seismic records, and Ln is the event-station path length as defined by


Formula 002

(2)
At the beginning of the inversion (i = 0), Formula was set to 270, which is the median value between the lowest and the highest Qo for all the paths. For the nth record, a residual term Formula is calculated using


Formula 003

(3)
A new estimate of Qm for each grid, Formula , is found by back- projecting Formula into the inverse of Qm. The residual term is redistributed back to the cells using


Formula 004

(4)
After experimentation, we found that 20 iterations (i = 20) were sufficient to provide stable Qm values.


    Results
 Top
 Abstract
 Introduction
 Data
 QLg Tomography
 Results
 Discussion and Conclusions
 Acknowledgments
 References
 
The inversion yielded a tomographic map of Qo for the Basin and Range Province (Fig. 2). The Qo map exhibits moderate variations in Qo within the area of study, with an average of 267 and a range of block values between 234 and 312. The result suggests a southwest to northeast trend in the attenuation. The highest Qo values (286–312) occur in the northeastern and eastern areas of the NTS, and they decrease smoothly to the southwest. The poor path coverage in the northeastern part of the region might have biased the Formula values we obtained for this area. All values higher than the average lie to the east and north of the NTS, while the lowest values are concentrated to the south and southwest of the NTS. The higher Formula values seem to be associated with the east-central Basin and Range and western part of the Colorado Plateau. This pattern of Formula is consistent with the distribution of Lg coda Q variations obtained by Baqer and Mitchell (1998).


Figure 002
View larger version (122K):
[in this window]
[in a new window]
 
Figure 2. Figure 002 from tomographic inversion of Lg amplitudes. Darker shades represent areas of relatively high Figure 002, and lighter shades show areas of low Figure 002, which probably indicates high crustal temperatures. The shades and Figure 002 contour lines are superimposed on a topography gradient map; dashed boundaries distinguish various tectonic provinces.

 

Mitchell (1995) and Mitchell and Cong (1998) found that crustal Qo values of any region are directly proportional to the length of time elapsed since the most recent crustal event that produced substantial metamorphism and associated fluid release. They considered regions (including the Basin and Range Province) with Q for Lg coda (Formula ) values between 250 and 333 as still tectonically active. Thus the relatively low Qo values found in this study would characterize a region that is still tectonically or orogenically active.

While the observed range of Qo values in the Basin and Range Province indicates a region of active deformation, these variations could also be due to other factors such as intrinsic crustal attenuation, increased relative stability and age of volcanism across the Basin and Range Province, and variations in the accumulated sediment thickness, permeability, and age. Two other competing mechanisms may also have caused enhanced attenuation of seismic energy in the Basin and Range Province. One possibility is the presense of interstitial fluids, due to either reactions at depth caused by high geothermal gradients or the circulation of fluids from shallow sources driven by the geothermal gradient (Wyllie, 1988). The second explanation for enhanced attenuation in the Basin and Range Province associates the region's large- scale extension, higher-than-normal heat-flow values, elevated upper-mantle temperatures, and abundant intrusive and extrusive igneous activity during the Cenozoic with the presence of fluid inclusions and partial melt.

Substantial frequency dependence of attenuation is also observed in the area as indicated by the {eta} values (Table 1). The {eta} values predominantly vary between 0.4 and 0.8 across the region, with an average of 0.57 and most falling in the range of 0.40 to 0.62. Studies have associated low {eta} values (<0.3) with shields and higher values with tectonically active regions (Mitchell, 1995; Baumont et al., 1999). Our higher {eta} values, suggesting a significant frequency-dependent attenuation of Formula , are therefore expected since the region is tectonically active.


    Discussion and Conclusions
 Top
 Abstract
 Introduction
 Data
 QLg Tomography
 Results
 Discussion and Conclusions
 Acknowledgments
 References
 
A preliminary study using Lg waves from 40 regional seismic events has investigated the crustal attenuation structure in the southwestern United States. Individual path attenuation values were determined using a new genetic algorithm that simultaneously fits the source spectra for path attenuation (Formula and the frequency dependence {eta}) as well as corner frequency and magnitude (Al-Eqabi et al., 2001). These Formula values have been inverted for the crustal attenuation structure of a portion of the Basin and Range region.

While uncertainties in the attenuation measurements prevent a direct tectonic interpretation for each anomaly on the map, the major trends of the inversion are likely robust, and agree well with previous studies. The Formula values for individual paths vary between 140 and 350 for frequencies from 0.3 to 10 Hz, and when inverted, the mean value for our Formula map (267) agrees very well with previous estimates, which are listed in Table 2. For instance, the study of Benz et al. (1997) obtained a Q-value of 235 ± 11 for the Basin and Range Province. Several studies have determined Q for Lg waves (Formula ) or coda (Formula ) in the Basin and Range Province. Reported Formula and Formula values vary between 138 and 774, but most studies put the Q-values between 200 and 280. In addition, the significant southwest–northeast trend in increasing Q, ranging between 234 and 312, agrees with the previous work of Baqer and Mitchell (1998).


View this table:
[in this window]
[in a new window]
 
Table 2 Summary of Attenuation Parameters in the Basin and Range Province
 

We tested the success of the tomographic attenuation model in fitting the amplitude data by examining the reduction of the variance of the individual data. We first computed the variance of the data relative to a uniform model with Q = 267 everywhere (summing the squares of the differences between the individual path Q–1 values and the averaged Q–1 = 1/267). We then determined Q–1 for each path through the tomographic model, by summing Q–1 for each block segment along a path, and computed the variance of the data path Q–1 values and the tomographic-model Q–1 values. Using the Q–1 values computed from the tomographic model reduced the variance of the data by 11% to a level of 89% of the variance that occurred just using the average Q-value of 267.

The majority of {eta} values fall between 0.4 and 0.8 (Table 1), suggesting a strong frequency-dependent attenuation in the Basin and Range Province. The low Q-values also indicate a region that is still tectonically active. The lowest Qo values occur in the south and southwest part of the Basin and Range Province where active tectonic deformation is occurring today. Higher values occur in the east-central Basin and Range and the western margin of the Colorado Plateau. The areas with high Qo values underwent volcanism between 17 and 55 m.y.a., while volcanism in low Qo regions occurred less than 17 m.y.a. (Blackwell, 1978).


    Acknowledgments
 Top
 Abstract
 Introduction
 Data
 QLg Tomography
 Results
 Discussion and Conclusions
 Acknowledgments
 References
 
We thank William Walter for providing the set of Lawrence Livermore National Laboratory seismograms, and for a very thorough review. We also thank Robert B. Herrmann for providing us with programs to isolate the Lg spectra, Keith Koper for the genetic algorithm code used in this study, Patrick Shore for computer support, and Brian Shiro for his review of this work. This research was supported by the National Science Foundation, Grant Number EAR-9629018, and the David and Lucile Packard Foundation.

Manuscript received May 6, 2004


    References
 Top
 Abstract
 Introduction
 Data
 QLg Tomography
 Results
 Discussion and Conclusions
 Acknowledgments
 References
 

Al-Eqabi, G. I., K. Koper, and M. E. Wysession (2001). Source characterization of Nevada Test Site explosions and western U.S. earthquakes using Lg waves: implications for regional source discrimination, Bull. Seism. Soc. Am.91 ,140 –153.[Abstract/Free Full Text]

Atkinson, G. M., and R. F. Mereu (1992). The shape of ground motion attenuation curves in southeastern Canada, Bull. Seism. Soc. Am. 82,2014 –2031.[Abstract/Free Full Text]

Baqer, S., and B. J. Mitchell (1998). Regional variation of Lg coda Q in the continental United States and its relation to crustal structure and evolution, Pure Appl. Geophys. 153,613 –638.[CrossRef]

Baumont, D., A. Paul, S. Beck, and G. Zandt (1999). Strong crustal heterogeneity in the Bolivian Altiplano as suggested by attenuation of Lg waves, J. Geophys. Res. 104,20287 –20305.[CrossRef]

Benz, H. M., A. Frankel, and D. M. Boore (1997). Regional Lg attenuation for the continental United States, Bull. Seism. Soc. Am.87 ,606 –619.[Abstract/Free Full Text]

Blackwell, H. M. (1978) Heat flow and energy loss in the western United States, in Cenozoic Tectonics and Geophysics of the Western Cordillera, R. B. Smith and G. P. Eaton (Editors), Geol. Soc. Am. Memoir, Vol. 152,175 –208.

Chael, E. (1987). Spectral scaling of earthquakes in the Miramichi region of New Brunswick, Bull. Seism. Soc. Am.77 ,347 –365.[Abstract/Free Full Text]

Chavez, D. E., and K. F. Priestley (1986). Measurement of frequency dependent Lg attenuation in the Great Basin, Geophys. Res. Lett.13 ,551 –554.[ISI][GeoRef]

Cong, L., and B. J. Mitchell (1998). Lg coda Q and its relation to the geology and tectonics of the Middle East, Pure Appl. Geophys.153 ,563 –585.[CrossRef]

Mitchell, B. J. (1991). Frequency dependence QLg and its relation to crustal anelasticity in the Basin and Range Province, Geophys. Res. Lett.18 ,621 –624.[ISI][GeoRef]

Mitchell, B. J. (1995). Anelastic structure and evolution of the continental crust and upper mantle from seismic surface wave attenuation, Rev. Geophys.33 ,441 –462.[CrossRef]

Mitchell, B. J., and L. Cong (1998). Lg coda Q and its relation to the structure and evolution of the continents: a global perspective, Pure Appl. Geophys. 153,655 –663.[CrossRef]

Singh, S. K., and R. B. Herrmann (1983). Regionalization of crustal coda Q in the continental United States, J. Geophys. Res. 88,527 –538.[CrossRef][ISI]

Wyllie, P. J. (1988). Magma genesis, plate tectonics, and chemical differentiation of the Earth, Rev. Geophys.26 ,370 –404.

Xie, J., and B. J. Mitchell (1990a). A back-projection method for imaging large-scale lateral variations of Lg coda Q with application to continental Africa, Geophys. J. Int. 100,161 –181.[CrossRef]

Xie, J., and B. J. Mitchell (1990b). Attenuation of multiphase surface waves in the Basin and Range Province, part I, Lg and Lg coda, Geophys. J. Int. 102,121 –137.[CrossRef]




This article has been cited by other articles:


Home page
Bulletin of the Seismological Society of AmericaHome page
S. R. Ford, D. S. Dreger, K. Mayeda, W. R. Walter, L. Malagnini, and W. S. Phillips
Regional Attenuation in Northern California: A Comparison of Five 1D Q Methods
Bulletin of the Seismological Society of America, August 1, 2008; 98(4): 2033 - 2046.
[Abstract] [Full Text] [PDF]


Home page
Bulletin of the Seismological Society of AmericaHome page
R. R. Castro, C. Condori, O. Romero, C. Jacques, and M. Suter
Seismic Attenuation in Northeastern Sonora, Mexico
Bulletin of the Seismological Society of America, April 1, 2008; 98(2): 722 - 732.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (3)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Aleqabi, G. I.
Right arrow Articles by Wysession, M. E.
Right arrow Search for Related Content
GeoRef
Right arrow GeoRef Citation


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS