versão impressa ISSN 0016-7169
Geofís. Intl v.48 n.3 México jul./set. 2009
Bodywave attenuation near the rupture of the 1887 Sonora, México, earthquake (Mw 7.5)
R. R. Castro1*, C. I. Huerta1, O. Romero2, C. Jacques3, A. Hurtado3 and A. I. Fernández1
1 Centro de Investigación Científica y de Educación Superior de Ensenada, División Ciencias de la Tierra, Departamento de Sismología, km 107 Carretera TijuanaEnsenada, 22860 Ensenada, Baja California, México *Corresponding author: firstname.lastname@example.org
2 Department of Geological Sciences, University of Texas, El Paso, TX799680555.
3 Instituto de Geología, Universidad Nacional Autónoma de México, Estación Regional del Noroeste, A. P. 1039, 83000 Hermosillo, Sonora, México.
Received: August 5, 2008
Accepted: March 18, 2009
Analizamos ondas de cuerpo de sismos locales y regionales registrados por la Red Sísmica del Noreste de Sonora, México (RESNES) usadas para determinar funciones empíricas de atenuación para las ondas P y S que describen el decaimiento de las amplitudes espectrales a frecuencias discretas entre 0.4 y 63 Hz. Los sismos seleccionados tienen su fuente cerca de la ruptura superficial del sismo de Sonora del 3 de mayo de 1887 (Mw 7.5), ubicado en la porción sur de la provincia tectónicofisiográfica de cuencas y cordilleras. Nuestros datos provienen de las componentes verticales y horizontales de sismogramas registrados a distancias hipocentrales de 25 a 150 km. Usamos funciones de atenuación obtenidas por Castro et al. (2008) con un modelo noparamétrico invirtiendo amplitudes espectrales observadas. Las curvas de atenuación reportadas indican que las amplitudes espectrales decaen más rápido durante los primeros 40 km de las trayectorias fuenteestación y muestran otro cambio en el porcentaje de decaimiento a aproximadamente 80 km. En este artículo separamos el efecto de Q de la corteza superior e inferior de la atenuación total: usamos los primeros 40 km de las curvas de atenuación para estimar Q de la corteza superior y el segmento de 4080 km para estimar Q de la corteza inferior. Encontramos que en la corteza superior la dependencia exponencial de Q con la frecuencia es la misma para ambas ondas P y S, específicamente Q ~ f 1.1, pero QS es un factor de 3 mayor, indicando que las ondas P experimentan mayor atenuación que las ondas S. En la corteza inferior Q es más compleja, ya que muestra una dependencia con la frecuencia más fuerte a bajas frecuencias (f < 3 Hz) y una débil dependencia con la frecuencia para las altas frecuencias (f > 6 Hz).
Palabras clave: Atenuación sísmica, Sonora, México.
We analized bodywaves of local and regional earthquakes recorded on the Seismic Network of Northeastern Sonora, México (RESNES) are used to determine empirical attenuation functions for P and S waves in the interval of 0.463 Hz. The epicenters are near to the surface rupture of the May 3, 1887 (MW 7.5) Bavispe, Sonora earthquake. Seismograms of vertical and horizontal components at of ground motion hypocentral distances of 25 to 150 km were used in a nonparametric model. The reported attenuation curves indicate that the spectral amplitudes decay faster in the first 40 km of the path and show different decay rate at approximately 80 km. We found that in the upper crust the exponential dependence of Q with frequency is the same for P and S waves, namely Q ~ f 1.1, but QS is larger by a factor of 3. In the lower crust Q is more complex, showing a strong frequency dependence below 3 Hz and a weaker frequency dependence above 6 Hz.
Key words: Seismic attenuation, Sonora, México.
The Sonora earthquake of May 3, 1887 is the largest historical seismic event occurred in the southern Basin and Range province. This normal fault earthquake destroyed the town of Bavispe, and had a maximum vertical slip at the surface of 5.1 m (Suter, 2006). Based on an endtoend rupture trace of 101.8 km, Suter (2006) estimated a momentmagnitude of Mw 7.5, in agreement with MS =7.5 estimated by Slemmons (1977) from maximum ground displacement. Collapse of adobe structures killed about fifty persons and intensities of III on the MMI scale were observed as far south as México city (Aguilera, 1920; Dubois and Smith, 1980). DuBois and Sbar (1979) estimated an average radius of felt reports of 700 km.
Seismic attenuation in the southern Basin and Range province for regional earthquakes is important for hazard evaluations. Condori (2006) and Castro et al. (2008) determined local and regional attenuation curves and described the decay of spectral amplitude with distance. Wich is faster near the rupture area of the 1887 earthquake compared with the surrounding region. They estimated the quality factor Q for both P and S waves using a homogeneous attenuation model.
The present paper discusses the role of seismic attenuation in the observed regional distribution of large intensities generated by the 1887 Sonora earthquake. We use local attenuation curves to study the variability of Q with depth for a twolayer model. Our model is based on the velocity model by Harder and Keller (2000) obtained from a wideangle seismic profile near the border between Sonora and Arizona. The uppermost layer of their model consists of sedimentary and volcanic rocks. The second layer, with a thickness of 22 km, corresponds to the crystalline upper crust. The base of the lower crust (third layer) is at 38 km depth, identified as the Moho.
We used a subset of the data in Condori (2006) and Castro et al. (2008) (Fig. 1). They used records from RESNES (Red Sísmica del Noreste de Sonora) to determine nonparametric attenuation functions (NAF) for P and S waves. Castro et al. (2008) had separated the data into events located northwest of the array and another with local events located in the rupture zone of the 1887 earthquake. In this paper we only use the attenuation functions obtained by Castro et al. (2008) from local earthquakes.
The RESNES network (Fig. 1)monitors the three major normal faults that ruptured during the 1887 earthquake: Pitáycachi, Teras, and Otates at north, central, and southern segments of the rupture. The RESNES network consists of nine digital fourcomponent stations (Castro et al., 2002; Romero et al., 2004; Castro et al., 2008). For epicentral location some stations from the NARSBaja (Trampert et al., 2003) and RESBAN arrays located in Sonora are also used. Fig. 2 shows a 3component acceleration record, and the vertical component of ground velocity, from a local earthquake recorded at Ojos Calientes, located at the northern end of the array, next to the MéxicoUS border.
The data set used by Castro et al. (2008) included earthquakes with hypocentral distances between 25 and 150 km and codamagnitudes between 1.8 and 3.2. Regional attenuation curves were obtained from events with hypocentral distances greater than 100 km as measured from the center of the network. Local NAF were found for earthquakes located near the fault zone that ruptured in 1887, within a radius of 100 km from the center of the network. These empirical attenuation functions describe spectral amplitude decay at 23 different frequencies between 0.4 and 63.1 Hz.
In this paper we use the NAF reported by Castro et al. (2008) in the region of the 1887 earthquake. Fig. 3 shows a sample of these attenuation functions for 8 different frequencies between 1 and 50 Hz for P (solid lines) and S waves (dashed lines). Note that the rate of amplitude decay with distance is not constant. However, at about 80 km the decay functions become approximately constant with distance, probably due to less severe geometrical spreading. At short distances (r < 40 km), the raypaths sample the shallow part of the crust and the spectral amplitudes decay faster than at longer distances, suggesting higher attenuation at shallow depths.
We model the first 40 km of the NAF as a product of three functions representing different attenuation mechanisms
where A(f,r) represents the attenuation curve at hypocentral distance r and at a fixed frequency f. The first righthand term accounts for geometrical spreading, the second term for amplitude decrease due to total Q (intrinsic plus scattering Q), and the third term for nearsurface attenuation. The parameter v is the average wave velocity in the upper crust, 3.23 km/s for S waves and 5.75 km/s for P waves. These average velocities were taken from Harder and Keller (2000). Based on their velocity model, we estimate that sourcestation paths in the distance range r< 40 km sample the upper 2325 km of the crust. In equation (1) R=(r10) and Ko is the average decay parameter.
The geometrical spreading function is approximated
where r ' is the distance at which the geometrical spreading starts having less effect upon the amplitude decay. We find r'=80 km based on the change of rate of amplitude decay in the attenuation functions (Fig. 3). This distance also marks the transition between bodywave and surfacewave geometrical spreading. In many regions this transition starts near 100 km, but Fig. 3 suggests that in northeastern Sonora the geometrical spreading transition may begin at a shorter distance range. Equations (1) and (2) are normalized at 10 km, the minimum sourcestation distance recorded by RESNES (see Fig. 2 in Castro et al., 2008). To determine the NAF, Castro et al., (2008) used acceleration spectra from selected events at hypocentral distances between 25 and 150 km. They assumed that at r = 0 there is not attenuation; consequently the attenuation functions were constrained to cross the origin. Thus, the NAF describe the spectral amplitude decay within the 0150 km range (Fig. 3).
Parameters Q and Ko may be estimated at each frequency by taking the logarithms of both sides of equation (1):
where a( R ) = log A(f,r)logG(r), m=π f log e /Q v, and b= π f log eKo. The slope m is estimated by a least square fit of equation (3). For a given frequency f, Q can be estimated as
Attenuation in the lower crust was obtained by correcting the original NAF for upper crust attenuation using the values of Q(f) and Ko(f) estimated from equations (4) and (5). We used the 4080 km segment of the corrected NAF to estimate total Q in the lower crust from the parametric model
where Al(f,r) are the amplitudes of the NAF corrected for upper crust attenuation, including near surface attenuation. We used an average velocity (v2) of 6.6 and 3.7 km/s for P and S waves, respectively, after Harder and Keller (2000). Q2 is the sum of intrinsic and scattering Q in the lower crust. We linearize equation (6) by taking logarithms and we search for values of b and Q2 that yield the best leastsquares solution for each frequency considered (e.g. Castro et al., 2004). We estimate that the 4080 km segment of the NAF yields the attenuation in the lower crust; thus we are sampling down to approximately 50 km depth.
Results and Discussion
A comparison between P and Swave NAF (Fig. 3) suggests that P waves may be more strongly attenuated than S waves. This observation agrees with the estimates of total Q from the 2layer model described above. From equations (1) and (6) and the tradeoff between geometrical spreading and Q, our Q estimates depend on the definition of geometrical spreading adopted. Also, our estimates of Q are based on the observed rate of amplitude decay represented by the NAF (Eq. 4), wich may be influenced by factors unrelated to attenuation. Thus in the 60100 km range some postcritical reflections may influence the spectral amplitude decay, although we selected time windows containing mainly direct arrivals to calculate the acceleration spectra.
Fig. 4 shows the estimates of Q for the upper and lower crust (Q > 0 ), and the Qfrequency relations from a leastsquares fit of the Q estimates. We were able to make estimates of QP in the frequency band 0.663 Hz and for QS between 1.0 Hz and 63 Hz. Fig. 4 shows estimates of Q for P waves (triangles) and for S waves (squares). In the upper crust (left frame) the exponential dependence of Q with frequency is the same for both, waves namely Q~f1.1,but Q is higher by a factor of 3 for P compared with S waves. In the lower crust (right frame) the frequency dependence of Q is more complicated and differs for P and S waves. In general, the frequency dependence of Q is stronger at low frequencies. At high frequencies, f > 7 Hz, QP is approximately constant, while QS remains frequency dependent. It is interesting that, as in the upper crust, P waves attenuate more strongly than S waves in the lower crust as a result of scattering and intrinsic Q. In the lower crust QP < 0 between 2.5 and 5.0 Hz for values of the geometrical spreading exponent b close to 1.0 (Fig. 5). For this frequency band Pwave geometrical spreading must be less severe (b smaller than 1.0), but the tradeoff between Q and b yields a better solution for Q < 0. Fig. 5 shows the exponent b in equation (6) as a function of frequency for P (triangles) and S (squares). In general the exponent b < 1.0. between 2 and 20 Hz b for S waves takes values close to 0.8.
The possible cause of the dependence of Q with frequency has been debated by different authors (Knopoff, 1964; Tsai and Aki, 1969; Aki and Chouet, 1975; Rautian and Khalturin, 1978; Aki, 1980; among others). Aki and Chouet (1975), suggested that the apparent frequency dependence of coda Q may be due to strong changes of Q with depth. The frequency dependence of Q has been also related to tectonic activity. For instance, Aki (1980) observed that in Hawaii, Q increases only slightly with frequency and is lower than in other areas. Under the Kilauea, Ellsworth (1977) observed a weaker heterogeneity in the mantle and this may cause less efficient bodywave scattering and also may explain the weak dependence of Q with frequency. In this volcanic region the ductility associated with the presence of magma and the high temperature could reduce smallscale heterogeneities that generate seismic wave scattering (Aki, 1980).
The southern Basin and Range in Sonora is a tectonic region under extension. High temperatures in the lower crust may play a role by reducing highfrequency scattering and caussing a weaker frequency dependence of Q, as in Fig. 4. Tectonic extension and magmatism in the Basin and Range province has generated an increase in surface heat flow relative to other more stable tectonic provinces such as the Colorado Plateau (Parsons, 1995).
The increased heat flow may be caused by lithospheric thinning (e.g. Crough and Thompson, 1976; Lachenbruch and Sass, 1978; Artyushkov and Batsanin, 1984). In contrast, in the upper crust, near the rupture zone of the 1887 Sonora earthquake, the presence of faults and the highly fractured rocks can generate strong bodywave scattering and lower Q. Comparing the values of total Q between upper and lower crust (Fig. 4), it is clear that attenuation is stronger in the upper crust. The implication is that seismic waves arriving at epicentral distances larger than 40 km will have less attenuated amplitudes because the rays will tend to travel along the lower crust, where Q is higher. This may help explain why, the felt intensity generated by the 1887 earthquake extended as far as it did. DuBois and Sbar (1979) estimated that the average radius of felt reports was about 700 km. These results also illustrate the important role played by the values of Q in the crust.
The NAF curves and the estimates of total Q obtained for the crust indicate that P waves attenuate more than S waves near the rupture zone of the 1887 Sonora earthquake. Sourcestation rays traveling in the lower crust attenuate less, due to a less severe scattering and internal friction, than rays that propagate in the upper crust. Consequently, at short epicentral distances the intensities attenuate faster than at intermediate (r > 40 km) and regional (r > 100 km) distances. Finally, the weaker highfrequency dependence of Q observed in the lower crust may be due to the ductility of the deeper rocks and high temperature that reduces smallscale heterogeneities.
We are thankful for the comments and insights of referee A which help us to improve the manuscript. This research was funded by the Mexican National Council for Science and Technology (CONACYT), grant 59216. Funding to maintain in operation the seismic stations of the RESB AN and NARSBaja networks comes from CONACYT, grants 24507 and 62116.LuisInzunza,AntonioMendoza,. Luis Inzunza, Antonio Mendoza, and Arturo PérezVertti provided technical support. The last version of this paper was prepared while one of the authors (RRC) was a UC MEXUSCONACYT visiting scholar fellow in the University of California, San Diego, SIOIGPP.
Aguilera, J. G., 1920. The Sonora earthquake of 1887, Bull. Seism. Soc. Am. 10, 556. [ Links ]
Aki, K. and B. Chouet, 1975. Origin of coda waves: source, attenuation and scattering effects, J. Geophys. Res. 80, 33223342. [ Links ]
Aki, K., 1980. Scattering and attenuation of shear waves in the lithosphere, J. Geophys. Res. 85, 64966504. [ Links ]
Artyushkov, Y. V. and S. F. Batsanin, 1984. Change in the thermal regime of the Eartht's crust associated with the approach of anomalous mantle to its lower boundary, Izvestiya, Earth Physics, 20, 887891. [ Links ]
Castro, R. R., O. M. Romero and M. Suter, 2002. Red sísmica para el monitoreo de la sismicidad del sistema de fallas normales del noreste de Sonora, GEOS 22, 379. [ Links ]
Castro, R. R., M. R. Gallipoli and M. Mucciarelli, 2004. An attenuation study in southern Italy using local and regional earthquakes recorded by the seismic network of Basilicata, Ann. Geofísica 47, 15971608. [ Links ]
Castro, R. R., C. Condori, O. Romero, C. Jacques and M. Suter, 2008. Seismic attenuation in northeastern Sonora, México, Bull. Seism. Soc. Am. 98, 722732. [ Links ]
Condori, C., 2006. Estudio de atenuación sísmica de la región noreste de Sonora, Tesis de Maestría, División Ciencias de la Tierra, CICESE, 88 pp. [ Links ]
Crough, S. T. and G. A. Thompson, 1976. Thermal model of continental lithosphere, J. Geophys. Res., 81, 48574862. [ Links ]
DuBois, S. and M. L. Sbar, 1979. The northern Sonora earthquake of 1887: hazard implications for Arizona, Geol. Soc. Am. Abstracts with Programs 11, 416. [ Links ]
DuBois, S. and A. W. Smith, 1980. The 1887 earthquake in San Bernardino Valley, Sonora, State of Arizona, Bureau of Geology and Mineral Technology, Special Paper 3, 112pp. [ Links ]
Ellsworth, W. L. 1977. Threedimensional structure of the crust and mantle beneath the island of Hawaii, Ph.D. thesis, 327 pp., Mass. Inst. of Technol., Cambridge. [ Links ]
Harder, S. and G. R. Keller, 2000. Crustal structure determined from a new wideangle seismic profile in southwestern New México, New México Geol. Soc. Guidebook 57,7578. [ Links ]
Knopoff, L., 1964. Q Rev. Geophys. 2, 625660. [ Links ]
Lachenbruch, A. H. and J. H. Sass, 1978. Models of an extending lithosphere and heat flow in the Basin and Range province, in Cenozoic tectonics and regional geophysics of the western Cordillera, Smith, R. B. and G. P. Eaton (Editors), Geol. Soc. Am. Mem. 152, 209250. [ Links ]
Parsons, T., 1995. The Basin and Range province, in Continental Rifts: Evolution, Structure, Tectonics, K.H. Olsen (Editor), Elsevier, Amsterdam, 277324 [ Links ]
Rautian, T. G. and V. I. Khalturin, 1978. The use of coda for determination of the earthquake source spectrum, Bull. Seism. Soc. Am. 68, 923948. [ Links ]
Romero, O. M., C. Jacques and R. R. Castro, 2004. Análisis de la sismicidad detectada por la red sismológica del noreste de Sonora, GEOS 24, 230. [ Links ]
Slemmons, D., 1977. State of the art for assessing earthquake hazards in the United Status: faults and earthquake magnitude, Rep. 6, U.S. Army Waterways Expt. Station, Vicksburg, Mississippi, 129 pp. [ Links ]
Suter, M. and J. Contreras, 2002. Active tectonics of northeastern Sonora, México (southern Basin and Range province) and the 3 May 1887 Mw 7.4 earthquake, Bull. Seism. Soc. Am. 92,581589. [ Links ]
Suter, M., 2006. Contemporary studies of the 3 May 1887 Mw 7.5 Sonora, México (Basin and Range Province) earthquake, Seism. Res. Lett. 77, 134147. [ Links ]
Trampert, J., H. Paulsen, A. Van Wettum, J. Ritsema, R. Clayton, R. Castro, C. Rebollar and A. PérezVertti, 2003. New array monitors seismic activity near the Gulf of California in México, EOS, Trans. Am. Geoph. Union 84, 2932. [ Links ]
Tsai, Y. B. and K. Aki, 1969. Simultaneous determination of seismic moment and attenuation of seismic surface waves, Bull. Seismol. Soc. Am. 59, 275287. [ Links ]