GEOPHYSICAL RESEARCH LETTERS, VOL. 36, L10308, doi:10.1029/2009gl037421, 2009 Apparent stress and corner frequency variations in the 1999 Taiwan (Chi-Chi) sequence: Evidence for a step-wise increase at M w 5.5 Kevin Mayeda 1 and Luca Malagnini 2 Received 22 January 2009; revised 13 April 2009; accepted 21 April 2009; published 28 May 2009. [1] Apparent stress and corner frequencies are measured for the Chi-Chi, Taiwan sequence beginning with the M w 7.6 mainshock on 20 September 1999. Using the coda source ratio method, we obtained stable source ratio estimates using broadband stations on Taiwan. We find the following: (1) For the mainshock and 3 of the larger aftershocks, apparent stress is clustered around 0.8 MPa (±0.1 MPa); (2) Events below M w 5.5 exhibit lower apparent stress with larger scatter, ranging between 0.08 and 0.8 MPa and are spatially variable; 3) The Brune [1970a, 1970b] omega-square source model fits the spectral shape for events 4.7 < M w < 7.6, however a step-wise break in self-similarity exists at M w 5.5. We hypothesize that larger events are subject to the regional state-of-stress, whereas smaller aftershocks are sensitive to the local state-of-stress from stressfield redistribution following the mainshock and/or fault zone lubrication that affects only larger events. Citation: Mayeda, K., and L. Malagnini (2009), Apparent stress and corner frequency variations in the 1999 Taiwan (Chi-Chi) sequence: Evidence for a step-wise increase at M w 5.5, Geophys. Res. Lett., 36, L10308, doi:10.1029/2009gl037421. 1. Introduction 1 Weston Geophysical Corporation, Lexington, Massachusetts, USA. 2 Istituto Nazionale di Geofisica e Vulcanologia, Rome, Italy. Copyright 2009 by the American Geophysical Union. 0094-8276/09/2009GL037421 [2] For many years and particularly over the past decade, earthquake source scaling studies have been undertaken in a variety of tectonic settings and in some cases, on the same data sets using different methodologies. More often than not, there is usually no clear consensus on scaling behavior and hence the underlying dynamics of faulting for that particular region. In addition to the academic pursuit of earthquake physics, there are the societal applications of differentiating explosions from natural earthquakes, predicting strong ground motion, and hazard mitigation. Though seismic network quality and station coverage has improved, and significant advancements have been gained in fault rupture models and high-resolution 3-D structure, the accurate prediction of strong ground motion requires a priori knowledge of the regional source scaling [e.g., Bodin et al., 2001; Akinci et al., 2001]. A recent study by Malagnini et al. [2008] demonstrates the necessity of inclusion of nonself similar scaling in order to accurately predict strong ground motion for the Colfiorito earthquake on September 26, 1997 (M w 6) in the central Apennines, Italy. [3] Nonetheless, remaining questions still exist such as: What is the current state-of-the-art?; Are systematic variations region-dependent?; Is the observed scatter real and, if so, what are the causes?; Can spatial variations in apparent stress be reliably mapped?; Are slow earthquakes, supershear earthquakes, and typical sub-shear tectonic earthquakes three different phenomena or part of a continuous spectrum of earthquake rupture behavior? [4] Aside from long-period regional and teleseismic waveform modeling estimates for M w, broadband studies that extrapolate back to the source have always been hampered by inadequate knowledge of scale-dependent path and site effects. Radiated energy and corner frequency estimation over a broad range of event sizes requires significant frequency-dependent corrections and assumptions, resulting in large error that makes interpretation of the results highly questionable. There are however a number of local and regional methods that circumvent the problem of path and site corrections, namely the empirical Green s function deconvolution method [e.g., Hough, 1997] and amplitude ratio techniques [e.g., Izutani, 2005]. These approaches have gained popularity because adjacent or co-located events recorded at common stations have shared path and site effects, which therefore cancel. In this paper we use a variant of the direct wave amplitude ratio, the coda ratio methodology of Mayeda et al. [2007] which has been shown to be roughly 3 times more stable than direct wave ratios and can be used with event pairs that are separated by several tens of kilometers, with little ill effects on amplitude ratio scatter. More recently, the method has been extended to datasets in Italy such as the San Giuliano di Puglia sequence in October November 2002 [Malagnini and Mayeda, 2008] and Colfiorito sequence of 1997 and 1998 [Malagnini et al., 2008]. [5] Over the past decade a number of studies have either suggested that earthquake scaling is constant with high scatter [e.g., Ide and Beroza, 2001; Prieto et al., 2006] or that scaling increases as a function of magnitude [e.g., Kanamori et al., 1993; Mayeda and Walter, 1996; Izutani, 2005]. However, due to large errors and regional variations, it has been difficult to definitively tell which of these two ideas is correct. Knowledge of scaling, whether self-similar or not, is fundamental to earthquake rupture simulations, seismic hazard prediction, and seismic discrimination, especially if there exists region-dependent variations in apparent stress. 2. Coda Amplitude Ratios [6] For this study we selected 18 events ranging between M w 4.7 and 7.6 recorded by 9 stations of the Broadband Array in Taiwan (BATS) (Figure 1). We applied the coda source methodology outlined by Mayeda et al. [2003] to L10308 1of5
the moment-rate functions for two events (1 and 2) is given by, 2 p=2 M 01 1 þ _M w = 1 ðwþ wc2 _M 2 ðwþ ¼ 2 p=2 ð1þ M 02 1 þ w = wc1 where M 0 is the seismic moment and w c is the angular corner frequency (2pf c ) and p is the high frequency decay rate. At the low frequency limit the source ratio shown in equation (1) is proportional to the ratio of the seismic moments [M 01 /M 02 ], whereas at the high frequency limit, equation (1) is asymptotic to [M 01 /M 02 ] ð 1 p 3Þ under selfsimilarity. If we follow the usual Brune [1970a, 1970b] omega-square model and set p = 2, the exponent of the high-frequency ratio becomes 1/3. However, it has been proposed by Kanamori and Rivera [2004] that the scaling between moment and corner frequency could take on the form, M o w ð 3þe c Þ ð2þ Figure 1. Map of Chi-Chi aftershocks that were used in this study. Symbol size is proportional to moment magnitude, M w, and shading is proportional to the log 10 [Apparent stress (Pa)]. compute stable M w s and validated these results against independent estimates from the Harvard CMT catalog for the largest events and selected M w s for the moderate-sized events from regional waveform modeling listed by Kao and Angelier [2001]. The coda ratio methodology is outlined by Mayeda et al. [2007] so we only give a brief processing description here. First, narrowband time-domain envelopes ranging between 0.03 and 8.0-Hz were made using the two horizontal components and log-averaged for additional stability and smoothed. Coda synthetic envelopes were then fit to the data for each station so that relative amplitudes could be measured using an L-1 fitting routine for each narrowband envelope [see Mayeda et al., 2003], then ratios were formed for all possible event pairs by subtracting the log 10 amplitudes for each station that recorded the event pair. The end of the coda window was based upon signal-tonoise ratio as well goodness of fit between the synthetic and the observed envelope. For each frequency band we also required that we had a minimum of 3 stations recording the event pairs where the maximum event spacing was limited to 60 km. We note that the vast majority of the events were within 30 km distance and we normally had a minimum of 6 stations recording the event pairs. We visually checked each frequency-dependent envelope and those that were prematurely short due to aftershocks in the coda or had dropouts were rejected. In total, this left us with 18 suitable ratios that had a minimum of roughly 6 stations to average over. [7] Assuming a simple single corner frequency source model [Aki, 1967; Brune, 1970a, 1970b], the ratio of where e represents the deviation from self-similarity and is usually thought to be a small positive number. For example, Walter et al. [2006] and Mayeda et al. [2005] found e to be close to 0.5 for the Hector Mine mainshock and its aftershocks using independent spectral methods. Mayeda et al. [2007] used the source spectrum portion of the Magnitude Distance Amplitude Correction (MDAC) methodology of Walter and Taylor [2001], which allows for the variation of the corner frequency that does not have to be self-similar. For example, w c ¼ ks 1=3 a ; s a ¼ s 0 M y 0 a M 0 M0 0 ; y ¼ e e þ 3 and k ¼ 16p ð3þ R 2 qfp z2 b 2 S a 5 S þ R2 qfs b 5 S where s a is the apparent stress [Wyss, 1970], s 0 a and M 0 0 are the apparent stress and seismic moment of the reference event, y is a scaling parameter, k is a constant related to the P and S- wave velocities at the source (a S, b S ), radiation pattern coefficient for P and S waves (R qfp and R qfs ), and scale factor (z) which is the ratio of the wave velocities. For the case of constant apparent stress, y = 0 and e = 0 in the previous equations, however, Mayeda and Walter [1996] found y = 0.25 for moderate-to-large earthquakes in the western United States and Mayeda et al. [2007] found y = 0.25 for the Hector Mine sequence using the coda ratio methodology. By using the corner frequency defined in equation (3) into equation (1), we can apply a grid search to find the parameters that best fit the individual spectral ratio data. We note however, the decay parameter p must be greater than 1.5 to keep the energy finite [e.g., Walter and Brune, 1993]. If the high-frequency decay value p were close to 1.5 it would be possible to nearly match the spectral ratios observed for Hector Mine, however, given that such shallow falloff is not consistent with most earthquake 2of5
Table 1. Source Parameter Information for Events Used in This Study Event M w f c Error f c MPa log 10 Ds app log 10 Ds app Error 992631747 7.60 0.037964 0.005857 5.85817 0.182852 992632312 4.85 0.711207 0.143167 5.54359 0.226294 992641812 5.05 0.384863 0.044540 5.05724 0.127113 992650013 6.30 0.169909 0.026288 5.86057 0.184166 992650048 5.85 0.302091 0.050970 5.93313 0.198450 992651217 5.05 0.761561 0.168219 5.92802 0.252971 992661244 5.05 0.598417 0.121341 5.61804 0.230245 992680842 5.15 0.531952 0.089953 5.62035 0.198095 992700355 4.70 0.690211 0.117706 5.28430 0.200888 992700720 4.95 0.505562 0.103510 5.24742 0.236861 992710546 4.85 0.693966 0.123247 5.51527 0.207271 992711537 4.85 0.425885 0.058248 4.88649 0.151852 992741247 5.15 0.467906 0.094100 5.44677 0.237266 992950217 5.90 0.264308 0.043720 5.83467 0.194255 993030826 5.20 0.453224 0.072976 5.48806 0.188683 993210734 5.15 0.479438 0.079193 5.48573 0.192368 001712155 5.00 0.655923 0.113220 5.66864 0.192658 002450924 4.95 0.431670 0.070715 5.04982 0.184263 observations [e.g., Hough, 2001] and independent methods [e.g., Mayeda et al., 2005; Walter et al., 2006], the preferred interpretation of Mayeda et al. [2007] was that the apparent stresses were systematically lower for the aftershocks than the mainshock. [8] Using as many as 9 stations, we formed the average spectral ratio for all Chi-Chi event pairs so long as the M w difference was at least 1.0. Next, we grid-searched using equations (1) and (3) assuming that the reference moment corresponded to an M w 5.0 event and the reference apparent stress was varied between 0.1 and 10 MPa. We note that our results are not dependent upon the choice of the reference moment. For every source pair, we obtain an estimate of the corner frequency for both events, then form averages and compute standard deviations. Next, we obtain estimates of the apparent stress from the spectral fits and the Brune stress drop from the corner frequency (see Table 1). As observed for other regions where this method has been applied, the coda spectral ratios for Chi-Chi are very stable, with average standard deviations around 0.1 0.2 for all frequencies. Though the Chi-Chi mainshock is complex with multiple sub-events [e.g., Kao and Chen, 2000; Chi et al., 2001], the broadband coda represents a time-domain convolution over the entire source duration. Due to crustal heterogeneity, a homogeneous scattered wavefield is quickly formed that is virtually free of source radiation and directivity effects. Figure 2 shows example spectral ratios between the mainshock and selected aftershocks. The codaderived source ratios exhibit little scatter and thus source parameters, such as corner frequency, will be better constrained when we fit the observed data with theoretical source models, such as the commonly used omega-square model [Aki, 1967; Brune, 1970a, 1970b]. As found in the Hector Mine sequence [e.g., Mayeda et al., 2007], there is an increase in the corner-frequency versus moment relation (see Figure 3). However, unlike Hector Mine, the Chi-Chi sequence was very rich in large magnitude aftershocks, in the M w 6+ range. From Figure 3 we see that the large events above M w 5.5, behave self-similarly whereas below this threshold, the event scatter is larger and the mean of the population decreases, representing a clear break in similarity. In an attempt to understand the observed variations, we also show in Figure 1 the apparent stress values for the aftershocks to see if there is any spatial pattern. Though there appears to be a hint that lower apparent stresses exist for events in the north, the region defined by high slip by finite-fault inversions by Chi et al. Figure 2. Example of coda-derived source ratios between the M w 7.6 mainshock and selected aftershocks. In each plot, we show the low and high frequency asymptotes from equation (1) as solid black horizontal lines and dashed lines represent the case when p = 1.5. Red lines represent ±1s about the average over a minimum of 6 stations and blue lines represent the best-fit theoretical ratio using equation (1). 3of5
Figure 3. Seismic moment (M o ) versus corner frequency (f c ) for the Chi-Chi events used in this study along with error estimates based on a bootstrap method. Lines of constant Brune stress drop are shown ranging between 0.1 and 100.0 MPa. We note that the mainshock and the three largest aftershocks appear to scale self-similarly, however a number of the aftershocks around M w 5 are quite variable, as evidenced in the spectral ratios in Figure 2. [2001], we cannot say anything definitively due to the complex deformation on Taiwan [Kao and Chen, 2000]. 3. Discussion and Conclusion [9] The averaging nature of coda waves has been shown to provide significantly lower amplitude variance than any traditional direct phase method. We have obtained stable source ratio estimates using broadband local and regional stations on Taiwan. We find the following: (1) For the M w 7.6 mainshock and 3 of the larger aftershocks (>M w 5.5), apparent stress is tightly clustered around 0.8 MPa (±0.1 MPa); (2) In contrast, events below moment magnitude M w 5.5 exhibit lower average apparent stress as well as larger scatter, ranging between 0.08 and 0.8 MPa and are spatially variable; (3) For this dataset, the Brune [1970a, 1970b] omega-square source model fits the spectral shape for events 4.7 < M w < 7.6, however a clear, step-wise break in self-similarity exists at around M w 5.5; (4) y =0is appropriate only for the mainshock and the largest aftershocks; (5) Using only the smaller events, extrapolation to larger events would yield unrealistically high apparent stresses and corner frequencies for the largest events; (6) We hypothesize that the larger events are subject to the average state-of-stress over a broader region, whereas the smaller aftershocks are more sensitive to the local stateof-stress resulting from stress-field redistribution following the mainshock. Recently, dynamic phenomena such as fault zone lubrication [e.g., Brodsky and Kanamori, 2001] as well as pore fluid pressurization (L. Malagnini, personal communication, 2008) have been proposed that may help to explain the observed step-wise change, though much work is needed. We note that other sequences where this method has been applied also suggests a step-wise break, albeit not as obvious due to insufficient numbers of large magnitude aftershocks. [10] In conclusion, the scaling results from this dataset are consistent with those from other sequences where this same method has been applied such as Hector Mine (CA), Wells (NV) (unpublished), San Giuliano (Italy), and Colfiorito (Italy). Rather than invoking a linear increase in Brune stress or apparent stress over the entire magnitude range, there appears to be a linear increase for events below M w 5.5 (non-self-similar), then constant above this magnitude (self-similar). However, that being said, there is strong evidence that the absolute levels of the stress are region dependent, as Hector Mine and Chi-Chi are roughly 5 and 10 times smaller, respectively, in Brune stress drop than the other 2 Italian sequences and Wells, NV. The compilation of more sequences is needed in conjunction with models that can account for the observed dynamic phenomena. [11] Acknowledgments. K. Mayeda was supported under Weston Geophysical subcontract GC19762NGD and AFRL contract FA8718-06- C-0024. Work by L. Malagnini was performed under the auspices of the Dipartimento della Protezione Civile, under contract S3 INGV-DPC (2007 2009), project: Valutazione rapida dei parametri e degli effetti dei forti terremoti in Italia e nel Mediterraneo. References Aki, K. (1967), Scaling law of seismic spectrum, J. Geophys. Res., 72, 1217 1231. Akinci, A., L. Malagnini, R. B. Herrmann, N. A. Pino, L. Scognamiglio, and H. Eyidogan (2001), High-frequency ground motion in the Erzincan region, Turkey: Inferences from small earthquakes, Bull. Seismol. Soc. Am., 91, 1446 1455. Bodin, P., L. Malagnini, and A. Akinci (2001), Ground-motion scaling in the Kachchh basin, India, deduced from aftershocks of the 2001 M w 7.6 Bhuj earthquake, Bull. Seismol. Soc. Am., 91, 1446 1455. Brodsky, E. E., and H. Kanamori (2001), Elastohydrodynamic lubrication of faults, J. Geophys. Res., 106, 16,357 16,374. Brune, J. N. (1970a), Tectonic stress and the spectra of seismic shear waves from earthquakes, J. Geophys. Res., 75, 4997 5009. Brune, J. N. (1970b), Tectonic stress and the spectra of seismic shear waves from earthquakes, J. Geophys. Res., 75, 4997 5009. (Correction, J. Geophys. Res., 76, 5002, 1971.) Chi, W.-C., D. Dreger, and A. Kaverina (2001), Finite-source modeling of the 1999 Taiwan (Chi-Chi) earthquake derived from a dense strongmotion network, Bull. Seismol. Soc. Am., 91, 1144 1157. Hough, S. E. (1997), Empirical Green s function analysis: Taking the next step, J. Geophys. Res., 102, 5369 5384. Hough, S. E. (2001), Empirical Green s function analysis of recent moderate events in California, Bull. Seismol. Soc. Am., 91, 456 467. Ide, S., and G. C. Beroza (2001), Does apparent stress vary with earthquake size?, Geophys. Res. Lett., 28, 3349 3352. Izutani, Y. (2005), Radiated energy from the mid Niigata, Japan, earthquake of October 23, 2004, and its aftershocks, Geophys. Res. Lett., 32, L21313, doi:10.1029/2005gl024116. Kanamori, H., and L. Rivera (2004), Static and dynamic scaling relations for earthquake and their implications for rupture speed and stress drop, Bull. Seismol. Soc. Am., 94, 314 319. Kanamori, H., E. Hauksson, L. K. Hutton, and L. M. Jones (1993), Determination of earthquake energy release and M L using TERRAscope, Bull. Seismol. Soc. Am., 83, 330 346. Kao, H., and J. Angelier (2001), Stress tensor inversion for the Chi-Chi earthquake sequence and its implications on regional collision, Bull. Seismol. Soc. Am., 91, 1028 1040. Kao, H., and W.-P. Chen (2000), The Chi-Chi earthquake sequence: Active, out-of-sequence thrust faulting in Taiwan, Science, 288, 2346 2349, doi:10.1126/science.288.5475.2346. Malagnini, L., and K. Mayeda (2008), High-stress strike-slip faults in the Apennines: An example from the 2002 San Giuliano earthquakes (southern Italy), Geophys. Res. Lett., 35, L12302, doi:10.1029/ 2008GL034024. Malagnini, L., L. Scognamiglio, A. Mercuri, A. Akinci, and K. Mayeda (2008), Strong evidence for non-similar earthquake source scaling in 4of5
central Italy, Geophys. Res. Lett., 35, L17303, doi:10.1029/ 2008GL034310. Mayeda, K., and W. R. Walter (1996), Moment, energy, stress drop, and source spectra of western United States earthquakes from regional coda envelopes, J. Geophys. Res., 101, 11,195 11,208. Mayeda, K., A. Hofstetter, J. L. O Boyle, and W. R. Walter (2003), Stable and transportable regional magnitudes based on coda-derived momentrate spectra, Bull. Seismol. Soc. Am., 93, 224 239. Mayeda, K., R. Gök, W. R. Walter, and A. Hofstetter (2005), Evidence for non-constant energy/moment scaling from coda-derived source spectra, Geophys. Res. Lett., 32, L10306, doi:10.1029/2005gl022405. Mayeda, K., L. Malagnini, and W. R. Walter (2007), A new spectral ratio method using narrow band coda envelopes: Evidence for non-selfsimilarity in the Hector Mine sequence, Geophys. Res. Lett., 34, L11303, doi:10.1029/2007gl030041. Prieto, G. A., D. J. Thomson, F. L. Vernon, P. M. Shearer, and R. L. Parker (2007), Confidence intervals of earthquake source parameters, Geophys. J. Int., 168, 1227 1234, doi:10.1111/j.1365-246x.2006.03257.x. Walter, W. R, and J. N. Brune (1993), Spectra of seismic radiation from a tensile crack, J. Geophys. Res., 98, 4449 4459. Walter, W. R., and S. R. Taylor (2001), A revised magnitude and distance amplitude correction (MDAC2) procedure for regional seismic discriminants: Theory and testing at NTS, Lawrence Livermore Natl. Lab. Rep. UCRL-ID-146882, Lawrence Livermore Natl. Lab., Livermore, Calif. (Available at http://www.llnl.gov/tid/lof/documents/pdf/240563.pdf) Walter, W. R., K. Mayeda, R. Gök, and A. Hofstetter (2006), The scaling of seismic energy with moment: Simple models compared with observations, in Earthquakes: Radiated Energy and the Physics of Faulting, Geophys. Monogr. Ser., vol. 170, edited by R. Abercrombie et al., pp. 25 41, AGU, Washington, D. C. Wyss, M. (1970), Stress estimates for South American shallow and deep earthquakes, J. Geophys. Res., 75, 1529 1544. L. Malagnini, Istituto Nazionale di Geofisica e Vulcanologia, via di Vigna Murata, 605, I-00143 Roma, Italy. K. Mayeda, Weston Geophysical Corporation, 181 Bedford Street, Suite1, Lexington, MA 02420, USA. (kmayeda@yahoo.com) 5of5