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. 176-187; DOI: 10.1785/0120050001
© 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 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 (1)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Singh, S. K.
Right arrow Articles by Iglesias, A.
Right arrow Search for Related Content
GeoRef
Right arrow GeoRef Citation

Article

An Estimate of Shear-Wave Q of the Mantle Wedge in Mexico

S. K. Singh1, J. F. Pacheco1, D. García2 and A. Iglesias3

1 Instituto de Geofisica
Universidad Nacional Autónoma de México
Del. Coyoacán, C.P. 04510
Mexico, D.F.
 (S.K.S., J.F.P.)

2 Departmento Geofisica y Meteorología
Facultad de Ciencias Físicas
Universidad Complutense de Madrid
Av. Ciudad Universitaria
28040 Madrid, Spain
 (D.G.)

3 Instituto de Ingeniería
Universidad Nacional Autonóma de México
Del. Coyoacán, C.P. 04510
Mexico, D.F.
 (A.I.)


    Abstract
 Top
 Abstract
 Introduction
 Data
 Analysis
 Results
 Discussion and Conclusions
 Appendix
 
We utilize earthquake recordings at two broadband seismographs located on the coast of the Gulf of Mexico to estimate Q of the mantle wedge. From the data at a station situated in the Laguna Verde nuclear power plant (LVIG), about 100 km northwest of Veracruz, and at the eastern edge of the Mexican Volcanic Belt, we estimate an upper bound of shear-wave Q of the mantle wedge, Q(f) ~ 120f0.75 (0.1 ≤ f ≤ 10 Hz), as compared with Q(f) = 251 f0.58 for the average path through subducted slab and continental lithosphere. Curiously, the estimated Q of the mantle wedge at station SCIG, which is located on the Yucatan Block, is about the same as the average Q in southern Mexico. There are several possible explanations for the difference: (1) very few recordings at SCIG, (2) a site effect at SCIG masking the effect of low Q of the mantle wedge, and (3) relatively high-mantle Q beneath the Yucatan Block. The third possibility is supported by surface-wave tomography that reveals a thick, cold, mantle lithosphere below the Yucatan Block and near absence of mantle lithosphere in the backarc region of central Mexico. A higher density of seismographs along the gulf coast is needed to resolve these issues. Our study predicts diminished ground motions at the Laguna Verde nuclear power plant if the seismic waves pass through the mantle wedge or traverse below the Popocatepetl or Orizaba volcanoes.


    Introduction
 Top
 Abstract
 Introduction
 Data
 Analysis
 Results
 Discussion and Conclusions
 Appendix
 
It is well known that inslab seismicity in the Michoacan and Guerrero segments of the Mexican subduction zone ceases before reaching the volcanic front. In these segments, the subducted Cocos plate bends with a small initial dip, then unbends and becomes subhorizontal (Suarez et al., 1990; Singh and Pardo, 1992; Pardo and Suarez, 1995). The maximum depth of the inslab events in the Michoacan and Guerrero segments is less than about 80 km. There is no seismic definition of the slab below this depth. As a consequence, the location and geometry of the mantle wedge in this region is poorly known. A shallow-dipping (~20°) slab without unbending has been proposed in central Oaxaca to explain the observed seismicity (Singh et al., 1985; Pardo and Suá rez, 1995). The inslab seismicity in this segment extends to about 400 km from the trench and reaches a depth of about 140 km. Thus, it is possible to delineate the mantle wedge here.

Low Q in the mantle wedge was postulated and confirmed in the early days of plate tectonics (e.g., Oliver and Isacks, 1967; Baranzagi et al., 1973). In recent years, the Q of the mantle wedge below Japan has been mapped in great detail (see, e.g., Takanami et al., 2000; Tsumura et al., 2000), thanks to the favorable location of seismic stations with respect to inslab seismicity and a high density of these stations. Low Q of the mantle wedge implies high attenuation of seismic waves reaching the backarc. An impressive demonstration of this comes from the extensive recordings of recent Japanese inslab earthquakes. A rather low Q for S waves (~60) in the mantle wedge is required to explain recorded seismic waves on the western seaboard of Japan, across the volcanic front, as compared with eastern Japan where the high-frequency seismic waves propagate through the subducted slab and the continental lithosphere with much less attenuation (T. Furumura, personal comm., 2004).

Mapping of Q and other physical properties of the mantle wedge in Mexico has been difficult because of the lack of suitably distributed seismic stations in the backarc region and because the ideal location for some of these seismographs would be in the Gulf of Mexico. Only three seismic stations are currently in operation along the gulf coast: LVIG, TUIG, and SCIG (Fig. 1). Each of these stations is equipped with a broadband seismograph and an accelerograph. TUIG is a highly noisy site located on thick gulf-coast sediments. For this reason, the seismograms from this site are of limited use. SCIG is located on limestone. LVIG, situated in the Laguna Verde nuclear power plant facility, is located on basaltic flows underlain by vulcanites that overlie a granitic inclusion. Logistical difficulties and instrumental failures, however, have resulted in few recordings at these stations, especially at SCIG.


Figure 001
View larger version (65K):
[in this window]
[in a new window]

 
Figure 1. Map showing locations and focal mechanisms of the earthquakes analyzed in this study. Gray zone delineates Mexican Volcanic Belt. Isodepth contours of the Benioff zone (modified from Singh and Mortera, 1991; Pardo and Suarez, 1995) are marked (dashed where inferred). LVIG and PNIG are seismic stations at Laguna Verde nuclear power plant and Pinotepa Nacional, respectively. Seismograms from TUIG, located on a noisy site, are not used in the analysis. SCIG recorded only three of the events listed in Table 1. A to D refer to the four groups into which the events have been divided. Straight dashed lines, denoted A to D, from the trench to LVIG indicate the locations of the vertical sections on which the events are projected in Figure 4.

 


View this table:
[in this window]
[in a new window]

 
Table 1 Earthquakes Analyzed in This Study*
 
In this study we analyze available seismograms of 14 inslab and 2 interplate earthquakes recorded at LVIG and other stations of the Mexican seismological and accelerometric networks. Three of these 16 events were also recorded at SCIG. Our goal is to obtain a preliminary estimation of Q of the mantle wedge. The seismic waves from these earthquakes to LVIG and SCIG travel, at least in part, through the mantle wedge (Fig. 1) where Q is expected to be lower than in the subducted slab and continental lithosphere.


    Data
 Top
 Abstract
 Introduction
 Data
 Analysis
 Results
 Discussion and Conclusions
 Appendix
 
The 16 earthquakes analyzed in this study are listed in Table 1 and shown in Figure 1. Our analysis includes all interplate and inslab earthquakes of Mexico that were recorded at LVIG. Only three of these events were also recorded at SCIG. We use these data in our study. Some recordings at SCIG are from interplate and inslab events of Chiapas (with epicenters east of 94° W; Fig. 1). Because of a very sparse network of stations in this region, it is not possible to determine whether the Q is anomalous for a wave path to SCIG. For this reason, these events have been excluded from this study.

Two of the 16 events analyzed here are shallow, thrust- faulting, interplate earthquakes (1 January 2004 and 14 June 2004); all others are inslab earthquakes (10 normal-faulting and 4 steeply dipping thrust-faulting inslab events). Our analysis is based on recordings from "hard" rock sites.

The events have been divided in four groups, A to D, indicated in Figure 1. Figure 2 shows the distances at which each event was recorded at all stations in the Mexican seismological and accelerometric networks. Figure 3 illustrates the paths from the epicenters to the stations. The events of each group were projected on a vertical plane that extends from the trench to LVIG. These sections are shown in Figure 4.


Figure 002
View larger version (28K):
[in this window]
[in a new window]

 
Figure 2. The distance coverage of the data available for events of groups A to D. For each event the source spectra retrieved from LVIG and, if available, from PNIG (marked by filled circles and triangles, respectively) are compared with the median spectrum obtained from all other stations of Mexican seismological and accelerometric stations (shown by open circles). Group A and C events provide an estimate of Q of the mantle wedge. Station SCIG is marked by filled rectangles in the figure.

 

Figure 003
View larger version (96K):
[in this window]
[in a new window]

 
Figure 3. (Top) Paths of inslab events used by García et al. (2004) to estimate Q and source spectra. Only 5 of the 268 records used in the study were from LVIG. Dots, earthquakes; triangles, stations. (Bottom) Paths of events included in the present study but not in García et al. (2004). The median spectra of these events were computed excluding data from LVIG and SCIG and are not appreciably affected by paths through the wedge.

 

Figure 004
View larger version (41K):
[in this window]
[in a new window]

 
Figure 4. Projection of the events of each group on a vertical plane extending from the trench to the station LVIG (see Fig. 1). All events are located either in the subducted slab or on the plate interface. For some events the plots do not show this. This is because the events are projected on sections that are not normal to the trench and the isodepth contours of the Benioff zone. Open circles, normal fault; filled circles, reverse fault.

 

Group A comprises five relatively deep inslab earthquakes that occurred near the northern Oaxaca-Veracuz border and in the Isthmus of Tehuantepec. The only three recordings available at SCIG are from events in this group. The waves from these sources travel through the mantle wedge and the continental lithosphere before reaching LVIG (Figs. 1 and 4) and SCIG. Group B consists of four inslab events that occurred near the coast of central Oaxaca and one shallow, thrust-faulting interplate event (depth H = 20 km). It includes the 30 September 1999 (Mw 7.4, H = 47 km) normal-faulting earthquake (Singh et al., 2000; Hernandez et al., 2001). The wave path from group B events to LVIG is probably mostly confined to the continental lithosphere. Group C consists of two inslab events: the 15 June 1999 (Mw 6.9, H = 61 km) Tehuacan earthquake (Singh et al., 1999) and the 21 July 2000 (Mw 5.9, H = 50 km) Copalillo earthquake (Iglesia et al., 2002). Group D includes three inslab, normal-faulting events located below the Balsas basin and one shallow, thrust-faulting interplate event (H = 17 km).

From Figures 1 and 4 it is at once clear that the waves from group A and C events to LVIG spend a larger fraction of their total travel time in the mantle wedge than group B and D events. Indeed, the events of groups B and D may barely sample the mantle wedge.


    Analysis
 Top
 Abstract
 Introduction
 Data
 Analysis
 Results
 Discussion and Conclusions
 Appendix
 
To examine whether the waves partly traveling through the mantle wedge are more attenuated than the average path, we reduced the observed S-wave spectrum at LVIG and SCIG to the source acceleration spectrum, S(f), and compared this spectrum with the median source spectrum estimated from the data of all the other stations. In computing S(f) of inslab earthquakes we used Q = 251f0.58, which is based on a recent study of inslab earthquakes of Mexico (García et al., 2004). If the Q of the mantle wedge is about the same as the Q for the average path (i.e., Q = 251f0.58), then the S(f) retrieved from LVIG and SCIG would be nearly equal to the median S(f); a lower or higher level of the S(f) at LVIG and SCIG would suggest a lower or higher Q along the path to LVIG and SCIG. For the two interplate earthquakes studied here (events B5 and D4, Table 1), we take Q = 273f0.66 (Ordaz and Sigh, 1992).

The Q from inslab events was estimated by García et al. (2004) for paths shown in Figure 3 (top). The figure suggests that the wave paths to stations, other than to LVIG, did not travel through the mantle wedge. Only 5 of the 268 records used in the study of García et al. (2004) were from LVIG. We recomputed Q by excluding data from LVIG. As expected, this resulted in an insignificant change in Q (from Q(f) = 251f0.58 obtained in García et al. to Q(f) = 256f0.59). Thus, we conclude that the Q reported in García et al. (2004) was not significantly affected by paths through the mantle wedge. Figure 3 (bottom) shows paths from events included in the present study but not in García et al. (2004). The median spectra of these events were computed, excluding data from LVIG and SCIG, and are not likely to be much affected by paths through the wedge. The median spectra of 6 of the 16 events are from García et al. (2004). The data of the 10 remaining events were processed in the same fashion as in García et al. (2004).

We follow a well-established procedure to compute S(f) (see, e.g., Singh et al., 1982; Ordaz and Singh, 1992; García et al., 2004) that we outline briefly here. S(f) is estimated from S-wave recordings on the two horizontal components. The horizontal components of the accelerations at each station were windowed. The window length was chosen such that it included the main S-wave arrival and ~85% of the total energy. The signals were then Fourier transformed, smoothed by a 1/6 octave-band filter, and 5% tapered. Let Ai(f, Ri) be the Fourier acceleration spectrum at station i located at a hypocentral distance of Ri. Then Ai(f, Ri) may be written as


Formula 001

(1)
where,


Formula 002

(2)
S(f), the source acceleration spectrum, may be written as


Formula 003

(3)
where M0(f) is the moment-rate spectrum. In the preceding equations, β = shear-wave velocity, {rho} = density, Q(f) = quality factor, R{theta}{phi} is the average radiation pattern (=0.55; Boore and Boatwright, 1984), F is the free-surface amplification (=2), P takes into account the partitioning of energy in the two horizontal components (=1/{surd}2), and M0 is the seismic moment of the earthquake. G(R) in equation (1) is the geometrical spreading term, which, for the inslab events, we take as G(R) = R (García et al., 2004). This form of G(R) implies dominance of body waves at all distances. For the interplate events, we take G(R) = R for R < 100 km and R1/2 for R ≥ 100 km (Ordaz and Singh, 1992). For inslab events we take β = 4.68 km/sec and {rho} = 3.2 g/cm3, whereas for the interplate events the corresponding values are taken as 3.5 km/sec and 2.8 g/cm3, respectively.


    Results
 Top
 Abstract
 Introduction
 Data
 Analysis
 Results
 Discussion and Conclusions
 Appendix
 
The source acceleration spectra of the events of groups A to D are shown in Figures 5, 6, 7, and 8. These figures illustrate the median ± one standard deviation source spectra (dashed curves). Henceforth, we denote the median source spectrum by SM(f) and the source spectrum obtained from the LVIG and SCIG recordings by SL(f) and SS(f), respectively.


Figure 005
Figure 005
View larger version (83K):
[in this window]
[in a new window]

 
Figure 5. The source spectrum obtained from the LVIG recording alone (continuous curve), and the median and ±one standard deviation source spectra (dashed curves) for the three events of group A. Dotted curves show source spectra retrieved from station PNIG. Left frames, all data corrected with Q(f) = 251f0.58; right frames, recordings at LVIG corrected with different Q(f) (given in the frames), for all other stations Q(f) = 251f0.58.

 

Figure 006
View larger version (35K):
[in this window]
[in a new window]

 
Figure 6. The source spectrum obtained from the LVIG recording alone (continuous curve), and the median and ±one standard deviation source spectra (dashed curves) for the events of group B. Dotted curves show source spectra retrieved from station PNIG. Note that a PNIG recording is not available for all events. Q(f) = 273f0.66 for the interplate B5 event; for all others events Q(f)=251f0.58.

 

Figure 007
View larger version (44K):
[in this window]
[in a new window]

 
Figure 7. Same as Figure 5 but for group C events.

 

Figure 008
View larger version (45K):
[in this window]
[in a new window]

 
Figure 8. Same as Figure 6 but for group D events. Here Q(f) = 273f0.66 for event D3, which is an interplate earthquake. For all other events Q(f) = 251f0.58.

 

For all five events of group A, SL(f) is consistently lower than SM(f) (see left frames of Fig. 5). A likely explanation for this difference is a smaller-than-average Q along that fraction of the path that traverses through the mantle wedge before reaching LVIG. If we assume that the lower Q in the mantle wedge is the correct explanation, then it is possible to establish a bound on its value. The source spectra from LVIG for events A2, A3, and A5 can be made to roughly coincide with the corresponding median spectra if Q along the path to LVIG is taken as Q(f) = 120f0.75 (right frames of Fig. 5). This provides an upper bound for the Q in the mantle wedge for the three events. The true Q(f) must be less because the waves to LVIG travel partly through the subducted slab and continental lithosphere. Let us assume that half of the total path to LVIG passes through the mantle wedge and the other half through the subducted slab and continental lithosphere. It is straightforward to infer that, in this case, the Q of the mantle wedge is roughly given by Q ~ 80f0.82.

Note that for SL(f) to match SM(f), the Q required for event A1 is much lower, Q(f) = 75f0.90, and for event A4 it is somewhat higher, Q(f) = 160f0.85. Local/regional and teleseismic recordings of event A1 show clear source directivity toward the southeast, away from LVIG. This may explain the lower spectral level at LVIG and, hence, the lower Q estimated for this as compared with the other three events. The cause for the relatively high Q for event A4 is not clear. We note that average S-wave radiation pattern of the group A events at LVIG is about 0.4, close to 0.55 taken in the calculations. Thus, a consistent bias in the radiation pattern cannot explain the observed lower spectral amplitudes at LVIG. Curiously, the source spectra of the three events retrieved from SCIG, SS(f), fall close to the median spectra SM(f). This suggests that the Q for paths to SCIG is close to the average Q or that there is a local site effect. We return to this issue in a later section.

Source spectra of events in group B are shown in Figure 6. SL(f) values of all events of this group fall within the curves defining ± one standard deviation from the median source spectra. Thus, the Q for paths from group B events to LVIG is roughly the same as for an average path, a result anticipated earlier based on the vertical section in Figure 4.

The wave paths to LVIG from the two events of group C are similar to those for events in group A (Figs. 1 and 4). Indeed, the SL(f) values as compared with the SM(f) values are similar to those for the group A (compare Figs. 5 and 7). Taking Q = 120f0.75 for the total paths to LVIG brings the SL(f) of the two events close to the SM(f) (see right frames in Fig. 7).

Group D events are farthest from LVIG (see Figs. 1, 2, 3, and 4). They range in depth from 20 to 64 km. Three of these events are inslab earthquakes (D1 to D3), and one of them is an interplate event (D4). In general, SL(f) is less than SM(f) for f >0.8 Hz and greater than SM(f) at lower frequencies (Fig. 8). At first glance this result is surprising because we anticipated group B and D events to show a similar trend. The wave paths from events D1, D2, and D3 to LVIG, however, pass below Popocatepetl, presently an active volcano (Fig. 1). The evidence suggests high attenuation of S waves crossing the volcano (Shapiro et al., 2000). Thus, SL(f) < SM(f) for f >0.8 Hz for these events may be a consequence of low Q below Popocatepetl. We note that event D4 shows the same trend as events D1 to D3 (Fig. 8), although its wave path does not cross the Popocatepetl volcano but the Orizaba volcano (Fig. 1). A likely explanation is that the Q below Orizaba is also low.

Amax at LVIG
These results predict lower-than-expected peak accelerations, Amax, at LVIG for the events of groups A, C, and D, and close to the expected values for group B events. Figure 9 shows plot of residuals (logarithm to the base 10 of the ratio of observed to expected Amax) for events of each group. Here the expected Amax has been computed for inslab and interplate events from regression relations of García et al. (2004) and Ordaz et al. (1989), respectively. As Figure 9 illustrates, however, the residuals at LVIG are negative for events of groups A and C, and are close to zero for events in group B. The tendency for group D events is less clear. Most data for group D have negative residuals. The residuals at LVIG, however, are not more negative, as we expected, but seem to follow the general tendency of the data.


Figure 009
View larger version (34K):
[in this window]
[in a new window]

 
Figure 9. Residual Amax as a function of the hypocentral distance for events of the four groups. Stations LVIG and PNIG are marked by filled circles and triangles, respectively.

 

Less-than-expected Amax values at LVIG for group A and C events is a welcome result for the seismic hazard for Laguna Verde nuclear power plant. The critical earthquake in the design of Laguna Verde nuclear power plant, however, was a local Mw 6.6 crustal event. Q in the wedge has no impact on the expected ground motions from such an earthquake.

Low Q in the Mantle Wedge or an Extremely Competent, "Hard" LVIG Site?
The left frames in Figure 5 also include source spectra obtained from PNIG, a station located in Pinotepa Nacional on the Pacific coast of Mexico (see Fig. 1 for its location). We note that the level of source spectra retrieved from LVIG (Fig. 5) are lower than those from PNIG, which, in turn, fall near minus one standard deviation curves of the median spectra. Could the low levels of source spectra at LVIG and PNIG be a consequence of both of these sites being extremely hard, LVIG being harder than PNIG? If so, then this could also explain negative residuals of Amax at both LVIG and PNIG for group A and C events (Fig. 9). In this scenario, the Q of the mantle wedge would be about the same as the average Q through the slab and the continental lithosphere. An examination of Figure 6, however, does not support this hypothesis. As seen in this figure, the source spectra from LVIG for three of the group B events are now higher than that those from PNIG. The simplest explanation, consistent with the observations, is to assume that (1) the LVIG site is similar to the typical hard site of the network, (2) PNIG is an extremely hard site, and (3) the Q in the mantle wedge is relatively low. This would explain the spectral curves of group B events in Figure 6, with the exception of B5. The wave paths from these events to LVIG involve little, if any, propagation through the mantle wedge. For this group of events, the spectral levels essentially reflect the site condition; higher and lower levels for LVIG and PNIG, respectively, are a result of larger and smaller amplification of seismic waves at these sites. For group A and C events, the levels at LVIG are lower than those at PNIG, which, in turn, is near minus one standard deviation curves of the median spectra. The levels at PNIG are low because it is a very hard site with relatively small amplification. The much lower level at LVIG is attributed to the low Q in the mantle. The spectral curves for group D events can be explained in much the same way as those for the other groups, except that the lower levels of LVIG are now attributed to propagation below the Popocatepetl or Orizaba volcanoes, regions of low Q.

Our assumption that PNIG is situated on an extremely hard site is supported by geologic and geophysical data. PNIG is situated on Xolapa Precambrian metamorphic complex. The P-wave velocity below the station is high, about 6.2 km/sec, and the Bouguer gravity anomaly in the area is anomalously positive, ~60 mGal (Valdes et al., 1986). The records of local earthquakes at PNIG are also very simple (Singh et al., 1997; Pacheco and Singh, 1998), suggesting a simple crustal structure and the absence of a significant weathered superficial layer. As mentioned earlier, LVIG is located on a 3-million-year-old basaltic flow, underlain by vulcanites that overly a granitic inclusion. It seems reasonable to expect that the site effect at LVIG be similar to an average site of the networks.

We estimated site effects at LVIG and PNIG from the spectral ratio of the horizontal-to-vertical component (H/V) of S waves. The results, illustrated in Figure 10, suggest little, if any, site effect at PNIG and a somewhat larger one at LVIG than at PNIG. Roughly, then, Figure 10 supports a low Q of the mantle wedge.


Figure 010
View larger version (17K):
[in this window]
[in a new window]

 
Figure 10. Horizontal-to-vertical spectral ratio of S-wave group at LVIG, PNIG, and SCIG as function of frequency. Median and ± one standard deviation curves are shown.

 

Q of the Mantle Wedge along Paths to Yucatan Block
As mentioned earlier, the spectra of the three events of group A recorded at SCIG (Fig. 5) do not reveal low Q along this path, which partly traverses through the mantle wedge (Fig. 1). SCIG is located on the Yucatan block. It is possible that the low mantle-wedge Q effect is masked by a site effect at SCIG. In fact, the H/V plot at SCIG (Fig. 10) does indicate a relatively large site effect. It is also possible that both the site effect at SCIG and the absence of a low Q mantle wedge along this path are responsible for the spectral level at SCIG.

Recent results of surface-wave tomography of Mexico and Central America reveal low and high shear-wave velocities in the backarc region northwest and southeast of the Tehuantepec ridge, respectively (N. Shapiro and M. Ritzwoller, personal comm., 2002). These authors propose near absence of mantle lithosphere below the Mexican backarc and the existence of substantial lithosphere below the Yucatan block. For group A events, this model predicts diminished spectral levels at LVIG and average or higher levels at SCIG. Thus, low Q for paths to LVIG is supported by this model.


    Discussion and Conclusions
 Top
 Abstract
 Introduction
 Data
 Analysis
 Results
 Discussion and Conclusions
 Appendix
 
Our results show that S-wave spectra from inslab earthquakes of group A are anomalously attenuated as they reach LVIG. This is also true for group C events. A significant fraction of wave paths from the events of these two groups to LVIG passes through the mantle wedge. Thus, the diminished spectra at LVIG may be attributed to the low Q of the mantle wedge. Group A events yield an upper bound on Q of the mantle wedge: Q(f) = 120f0.75 (0.1 ≤ f ≤ 10 Hz); a realistic estimation may be Q(f) = 80f0.82. These Q values are smaller than Q(f) = 251f0.58, the average Q for Mexican inslab earthquakes with wave paths mostly confined to the subducted slab and the continental lithosphere (García et al., 2004).

The three events of group A recorded at SCIG (Fig. 5) do not reveal low Q along this path, although it partly traverses through the mantle wedge (Fig. 1). SCIG is located on the Yucatan block. We consider the data from SCIG as inconclusive however, because of the small number of recordings and possible contamination from the local site effect. Only more extensive deployment of seismographs along the coast of Yucatan Block can resolve this important issue.

The waves from events of groups B and D to LVIG mostly travel through the continental lithosphere. Hence, we expect the Q for paths to LVIG to be similar to the Q of the average path. Indeed, the results for group B agree with this expectation. This is not true, however, for group D events. The waves from these events travel below the active Popocatepetl or Orizaba volcano before reaching LVIG and hence suffer high attenuation. This results in diminished amplitudes at LVIG. Low Q below Popocatepetl has been reported in a previous study (Shapiro et al., 2000). Our study suggests that the Q may be low below the Orizaba volcano also.

We have offered a qualitative explanation of the observations which are based on a small data set. A detailed tomography of the Q of the mantle wedge will require a much larger data set and generation of synthetic seismograms through appropriate 2D and 3D structures of the Mexican subduction zone (see, e.g., Furumura and Singh, 2002).

Our study suggests reduced seismic hazard at the Laguna Verde nuclear power plant from earthquakes with wave paths through the low-Q mantle wedge or the active volcanoes of Popocatepetl and Orizaba. The same must be true for other sites along the coast of Gulf of Mexico. The critical design earthquake for the nuclear plant, however, was a local, crustal earthquake. Seismic hazard from such an event remains the same.


    Appendix
 Top
 Abstract
 Introduction
 Data
 Analysis
 Results
 Discussion and Conclusions
 Appendix
 
We are grateful to T. Furumura, N. Shapiro, and M. Ritzwoller for making their results available to us prior to publication. As always, we are thankful to technicians who, under rather adverse conditions, maintain networks operated by Instituto de Geofísica (UNAM), Instituto de Ingeniería (UNAM), and Cenapred. F. Mooser provided information about the geology of Laguna Verde site. Critical comments by Anton Dainty and the two anonymous reviewers are appreciated. The research was supported, in part, by DGAPA, UNAM project JN114305 and CONACyT project 42671-F. D. García was supported, in part, in Mexico by Programa Predoctoral UCM fellowship.

Manuscript received January 3, 2005

Baranzagi, M., B. Isacks, J. Oliver, J. Dubois, and G. Pascal (1973). Descent of lithosphere below New Hebrides, Tonga-Fiji and New Zealand, Nature242 ,98 –101.[CrossRef]

Boore, D. M., and J. Boatwright (1984). Average body-wave radiation coefficient, Bull. Seism. Soc. Am.74 ,1615 –1621.[Abstract/Free Full Text]

Furumura, T., and S. K. Singh (2002). Regional wave propagation from Mexican subduction zone earthquakes: the attenuation functions for interplate and inslab earthquakes, Bull. Seism. Soc. Am. 92,2110 – 2125.[Abstract/Free Full Text]

García, D., S. K. Singh, M. Herráiz, J. F. Pacheco, and M. Ordaz (2004). Inslab earthquakes of central Mexico: Q, source spectra and stress drop, Bull. Seism. Soc. Am.94 ,789 –802.[Abstract/Free Full Text]

Hernandez, B., N. M. Shapiro, S. K. Singh, J. F. Pacheco, F. Cotton, M. Campillo, A. Iglesias, V. Cruz, J. M. Gómez, and L. Alcántara (2001). Rupture history of September 30, 1999 intraplate earthquake of Oaxaca, Mexico (Mw=7.5) from inversion of strongmotion data, Geophys. Res. Lett.28 ,363 –366.[CrossRef][ISI][GeoRef]

Iglesias, A., S. K. Singh, J. F. Pacheco, and M. Ordaz (2002). A source and wave propagation study of the Copalillo, Mexico earthquake of July 21, 2000 (Mw=5.9): Implications for seismic hazard in Mexico City from inslab earthquakes, Bull. Seism. Soc. Am. 92,1060 –1071.[Abstract/Free Full Text]

Oliver, J., and B. Isacks (1967). Deep earthquake zones, anomalous structure in the upper mantle, and in the lithosphere, J. Geophys. Res. 72,4259 –4275.[ISI]

Ordaz, M., and S. K. Singh (1992). Source spectra and spectral attenuation of seismic waves from Mexican earthquakes, and evidence of amplification in the hill zone of Mexico City, Bull. Seism. Soc. Am. 82,24 –43.[Abstract/Free Full Text]

Ordaz, M., J. M. Jara, and S. K. Singh (1989). Riesgo sísmico y espectros de diseño en el estado de Guerrero, Mem. VIII Congreso Nac. de Ingeniería Sísmica, Acapulco, Mexico,D40 –D56.

Pacheco, J. F., and S. K. Singh (1998). Source parameters of two moderate earthquakes estimated from a single-station, near-source recording, and from MT inversion of regional data: a comparison of results, Geofis. Int.37 ,95 –102.

Pardo, M., and G. Suárez (1995). Shape of the subducted Rivera and Cocos plates in southern Mexico: seismic and tectonic implications, J. Geophys. Res.100 ,12,357 –12,373.[CrossRef]

Shapiro, N., S. K. Singh, A. Iglesias-Mendoza, V. Cruz-Atienza, and J. F. Pacheco (2000). Evidence of low Q below Popocatépetl, and its implication to seismic hazard in Mexico City, Geophys. Res. Lett. 27,2753 –2756.[CrossRef][ISI][GeoRef]

Singh, S. K., and F. Mortera (1991). Source-time functions of large Mexican subduction earthquakes, morphology of the Benioff zone and the extent of the Guerrero gap, J. Geophys. Res. 96,21,487 –21,502.

Singh, S. K., and M. Pardo (1992). Geometry of the Benioff zone and the state of stress in the overriding plate in central Mexico, Geophys. Res. Lett.20 ,1483 –1486.[ISI]

Singh, S. K., G. Suárez, and T. Domínguez (1985). The great Oaxaca earthquake of 15 January 1931: lithosphere normal faulting in the subducted Cocos Plate, Nature317 ,56 –58.[CrossRef][GeoRef]

Singh, S. K., R. Apsel, J. Fried, and J. N. Brune (1982). Spectral attenuation of SH-wave along the Imperial fault, Bull. Seism. Soc. Am.72 , 2003– 2016.[Abstract/Free Full Text]

Singh, S. K., J. F. Pacheco, F. Courboulex, and D. A. Novelo (1997). Source parameters of the Pinotepa Nacional, Mexico, earthquake of 27 March, 1996 (Mw = 5.4) estimated from near-field recordings of a single station, J. Seism.1 ,39 –45.[CrossRef]

Singh, S. K., M. Ordaz, L. Alcántara, N. Shapiro, V. Kostoglodov, J. F. Pacheco, S. Alcocer, C. Gutierrez, R. Quaas, T. Mikumo, and E. Ovando (2000). The Oaxaca earthquake of September 30, 1999 (Mw = 7.5): a normal-faulting event in the subducted Cocos Plate, Seism. Res. Lett.71 ,67 –78.

Singh, S. K., M. Ordaz, J. F. Pacheco, R. Quaas, L. Alcántara, S. Alcocer, C. Gutiérrez, R. Meli, and E. Ovando (1999). A preliminary report on the Tehuacan, Mexico earthquake of June 15, 1999 (MW = 7.0), Seism. Res. Lett. 70,489 –504.

Suárez, G., T. Monfret, G. Wittlinger, and C. David (1990). Geometry of subduction and the depth of the seismogenic zone in the Guerrero gap, Nature345 ,336 –338.[CrossRef][GeoRef]

Takanami, T., I. Selwyn Sacks, and A. Hasegawa (2000). Attenuation structure beneath the volcanic front in northeastern Japan from broad-band seismograms, Phys. Earth Planet. Interiors121 ,339 –357.[CrossRef]

Tsumura, N., S. Matsumoto, S. Horiuchi, and A. Hasegawa (2000). Three- dimensional attenuation structure beneath the northeastern Japan arc estimated from spectra of small earthquakes, Tectonophysics 319,241 –260.[CrossRef][ISI][GeoRef]

Valdés, C., W. D. Mooney, S. K. Singh, R. P. Meyer, C. Lomnitz, J. H. Luetgert, C. E. Helsley, B. T. R. Lewis, and M. Mena (1986). Crustal structure of Oaxaca, Mexico, from seismic refraction measurements, Bull. Seism. Soc. Am.76 ,547 –563.[Abstract/Free Full Text]




This article has been cited by other articles:


Home page
Bulletin of the Seismological Society of AmericaHome page
S. K. Singh, A. Iglesias, D. Garcia, J. F. Pacheco, and M. Ordaz
Q of Lg Waves in the Central Mexican Volcanic Belt
Bulletin of the Seismological Society of America, August 1, 2007; 97(4): 1259 - 1266.
[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 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 (1)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Singh, S. K.
Right arrow Articles by Iglesias, A.
Right arrow Search for Related Content
GeoRef
Right arrow GeoRef Citation


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS