Synthetic sensitivity analysis of high frequency radiation of 2011 Tohoku-Oki (M W 9.0) earthquake

Transcription

1 Earthq Sci (214) 27(4): DOI 1.17/s RESEARCH PAPER Synthetic sensitivity analysis of high frequency radiation of 211 Tohoku-Oki (M W 9.) earthquake Haoran Meng Yongshun John Chen Received: 3 December 213 / Accepted: 16 May 214 / Published online: 26 June 214 The Seismological Society of China, Institute of Geophysics, China Earthquake Administration and Springer-Verlag Berlin Heidelberg 214 Abstract Frequency-dependent rupture behavior of subduction zone interplate megathrust faults has been observed by back-projection method in different frequency bands (from.5 to 5 Hz). It has been suggested that the down-dip region of the Tohoku megathrust radiated strongly at high frequencies ([1 Hz) compared with that of the up-dip region. By assuming the same source time function of each fault patch, we perform a synthetic sensitivity analysis to compare the energy received from the shallower parts (and further way from the receiver sites) with that from the deeper parts (and closer to the receiver sites) of the rupture. Our results indicate that regional onshore recordings are probably not adequate to constrain the presence of far-off shore high frequency radiations because of the strong attenuation of this region. Keywords Tohoku-Oki earthquake High frequency radiation Frequency-dependent rupture behavior 1 Introduction Frequency-dependent rupture process of the 211 M W 9. Tohoku-Oki earthquake (e.g., Lay et al., 212) has been examined using multi-frequency (.5 5 Hz) back-projection imaging (e.g., Koper et al. 211a, b; Kiser and Ishii 211; Wang and Mori 211; Ishii 211; Huang et al., 212; Meng et al, 212) and compressive sensing studies (e.g., Yao et al., 211), inferred to correspond to depth-varying rupture properties. This interpretation was corroborated by H. Meng Y. J. Chen (&) Institute of Theoretical and Applied Geophysics, School of Earth and Space Sciences, Peking University, Beijing 1871, China johnyc@pku.edu.cn analyzing seismic source spectra of interplate megathrust events in frequency bands up to 1 Hz. (Ye et al., 213). Analyses of strong motion data also indicate that higher frequency ([1 Hz) energy was radiated strongly at the down-dip region rather than at the up-dip region of the Tohoku-Oki megathrust fault (e.g., Ide et al. 211; Iwakiri and Hoshiba 212). Investigations of this phenomenon are important both for understanding of variations in basic source physics across the megathrust and the distance decay of strong shaking. Even though the robust results of teleseismic back-projection imaging and source spectra analyses strongly indicate frequency-dependent rupture process using data with relative low frequency bands (\1 Hz), the results of analyzing near-field data with higher frequency band could provide supplement about the intrinsic differences in rupture behavior across the megathrust, so as to help us improving dynamic source models. Therefore, it is significant to carefully examine whether using near-field data could solidly support the conclusion of back-projection. However, due to poor station coverage of near field observation (Fig. 1a) and geometry of the low angle interplate megathrust fault (Fig. 1c), other important factors, such as geometrical and anelastic attenuation as well as diffraction of seismic waves should also be taken into account in these analyses. Because of the geometrical settings (Fig. 1c; Table 1) of the Japan Trench megathrust and the stations on shore, high frequency radiations from the shallower part of the rupture attenuate much more along the path compared with that from deeper parts (Fig. 1c; Table 1). Rupture inversion results, however, showed that slip of the Tohoku earthquake (Fig. 2), increases to maximum from the down-dip region to trench (e.g., Ide et al. 211; e.g., Yue and Lay 213). The tradeoff between attenuation and slip distribution makes it less certain about

2 356 Earthq Sci (214) 27(4): (a) IWTH IWTH N (b) j 1,8 1,2 1,1 2,1 MYGH N IWTH14 MYGH N IWTH27 MYGH FKSH FKSH MYGH8 FKSH17 FKSH km 8,8 9,8 8,1 9,1 Trench IBRH N IBRH14 IBRH IBRH CHBH14 CHBH t (s) 35 N Slip (m) 2 km 139 E 14 E 141 E 142 E 143 E 144 E 16,8 24 km 15,1 16,2 16,1 i (c) Coast Overriding plate low-angle megathrust fault up-dip Trench East Oceanic plate down-dip rupture zone Fig. 1 a (Left) Strong motion records high-pass filtered above 1 Hz at nine stations distributed along the rupture. The maximum amplitude in centimeters per second per second is shown at the top-right of each waveform. (Right) Slip distribution. The S-wave arrival times from several source locations, denoted with colored stars and with time after the rupture initiation also shown, are plotted on the waveforms with length inversely proportional to the source-station distance, modified from (Ide et al. 211). b Mesh of the rupture region, model blocks of the central positions of 3 km 9 3 km patches. i and j increase along strike and along dip separately. c Schematic cross-section of interplate megathrust fault Table 1 Geometry of the interplate megathrust fault Depth (km) Distance from trench (km) Dip angle () [27 [ distribution pattern of high frequency radiation based only on simple frequency analysis of the data observed. It requires carefully quantitative analysis of the data taking into account also the wave propagation geometry and attenuation effects. Using a finite fault model (e.g., Yue and Lay 213) and synthetic seismograms, we have performed a synthetic sensitivity analysis to investigate whether data recorded from stations along the east coast of Honshu island could be analyzed to resolve the possible high frequency ([1 Hz) radiation from the up-dip region of the megathrust. By placing testing receivers in a 2D region from the island to the trench, we finally present the result of detectability of high frequency radiations ([1 Hz) in the shallower parts of the rupture around the Tohoku region.

3 Earthq Sci (214) 27(4): (a) Slip (m) (b) Moment density (112N/m) Distance along strike (km) Distance along strike (km) Fig. 2 a Inverted rupture pattern on the fault plane for the Tohoku-Oki earthquake, modified from (Yue and Lay 213). b Moment density distribution calculated from (a) and PREM model. Value of distance along dip increases from the trench to down-dip (a) Vertical Q=inf. (b) Horizontal Q=inf. (c) Total Q=inf. Distance along strike (km) Distance along strike (km) Distance along strike (km) Fig. 3 Synthetic test of received high frequency energy distribution with uniform slip of the megathrust fault when Q is infinite, station is set just above the southwest corner of the megathrust fault. a Vertical component, b Horizontal component, c Total. Apparently, the rsum k of vertical component is the highest so as to be the most sensitive component. Value of distance along dip increases from the trench to down-dip

4 358 Earthq Sci (214) 27(4): (a) IWTH14 (b) IWTH27 (c) MYGH8 (d) IBRH14 Distacne along strike (km) Distacne along strike (km) Distacne along strike (km) Distacne along strike (km) (e) CHBH Distacne along strike (km) Q=2 Q=4 Q=6 Q=8 Q=

5 Earthq Sci (214) 27(4): b Fig. 4 ORHFED patterns with different Q structures of station IWTH14 (a), IWTH27 (b), MYGH8 (c), IBRH14 (d), CHBH14 (e), north to south, respectively. E ð k i;j Þ of a certain fault patch, which has been normalized from to 1, represents the high frequency energy received by station k rather than the energy radiated from that fault patch. Value of distance along dip increases from the trench to downdip consideration that (1) geodetic and seismic data are jointly inverted in solution, yielding better slip estimation and resolution than other inversions using one type of data only. (2) This estimated slip distribution pattern is similar to that of (Ide et al. 211), what we will have our results to compare with. To generate synthetic seismograms with sampling rate of 4 Hz and duration of about 2 s, we adopt the wavenumber integration method because of its higher efficiency and broader frequency bands than spectral element or finite difference method. For simplicity and feasibility, we only consider layered homogeneous media following the PREM model between the source region and receivers. Due to the lack of data constraints, we do not consider lateral variations of Q structure. Although Q S /Q P generally ranges from 1 to 3 for frequencies higher than 1 Hz (e.g., Sato et al., 213), after propagating from source to receiver, the dominant energy of seismic wave mostly depends on shear wave or surface wave, while the effect of Q P is negligible. As a result, we simply assume a homogeneous Q structure where Q = Q S = Q P, independent of frequency. By varying Q, we obtain different Green s function G k ði;jþðqþ of media. We then high pass filter the accelerograms above 1 Hz and integrate the k-th station record from the fault patch (i,j) into its L2 norm E ð k i;j Þ, i.e., Z Te k i;jþ ¼ A k ði;j Þ ðtþ 2 dt ð2þ Fig. 5 The relation between rsum k and Q value. IWTH14 to CHBH14 are results of different stations from north to south, respectively. The black dashed line denotes the cut ratio.1 2 Methods To test the sensitivity of high frequency radiation, we parameterize the source region into 16 by 8 rectangle fault patches (Fig. 1b). Accelerograms of k-th station A k ði;jþ ðtþ from source patch (i,j) can be expressed as the convolution of Green s function G k ði;jþ ðtþ, source time function M (i,j)s(t) and instrumental response I(t): A k ði;jþ ðtþ ¼Gk ði;jþ ðtþ M ði;jþ d 2 SðtÞ dt 2 IðtÞ ð1þ where M (i,j) denotes the moment of each fault patch, S(t) is the is the Heaviside function to ensure production of abundant high frequency energy. We simply assume the same S(t) of each fault patch in order to investigate only the propagation effects. In this study we use the finite-fault slip model of (Yue and Lay 213) derived from a joint inversion of high-rate geodesy and seismic data. This is based on the T b where T b and T e denotes the begin and end time of the record of an earthquake. According to our synthetic test, the k i;jþ of vertical component is more sensitive to high frequency energy radiated from the shallower parts of the megathrust fault (Fig. 3) than the horizontal components. Plotting k i;jþ onto the fault plane, we can examine at each station the observed high frequency energy distribution (ORHFED) pattern on the fault plane (Fig. 4). 3 Discussion We calculate the ORHFED patterns for stations along the east coast of Honshu island shown in Fig. 1a, with the Q value ranges from 2 to 1, (Fig. 4). Generally speaking, when Q \ 6, received energy from the down-dip region of the megathrust dominates, while the received energy form the up-dip region area begins to be comparable with the energy from the down-dip region area as Q exceeding 6 for northern stations along the coast. The general trend indicates that, for the tradeoff between attenuation and slip distribution, the latter begins to

6 36 Earthq Sci (214) 27(4): (a) Stations 138 E 14 E 142 E 144 E 146 E (b) Q = 1 (c) Q = E 14 E 142 E 144 E 146 E 138 E 14 E 142 E 144 E 146 E 4 N 4 N 4 N (d) Q = 3 (e) Q = 4 (f) Q = E 14 E 142 E 144 E 146 E 138 E 14 E 142 E 144 E 146 E 138 E 14 E 142 E 144 E 146 E 4 N 4 N 4 N (g) Q = 6 (h) Q = 7 (i) Q = E 14 E 142 E 144 E 146 E 138 E 14 E 142 E 144 E 146 E 138 E 14 E 142 E 144 E 146 E 4 N 4 N 4 N (j) Q = E 14 E 142 E 144 E 146 E (k) Q = E 14 E 142 E 144 E 146 E 4 N 4 N

7 Earthq Sci (214) 27(4): b Fig. 6 a Distribution of testing stations. b k The ^r sum k maps for different Q values, where ^r sum k denotes the logarithm of rk sum. If the region has a ^r sum k lower than -1, station in this region could not resolve the high frequency energy radiated form the shallower parts of the megathrust fault. The write contours denote the slip distribution of Tohoku-Oki earthquake inverted by (Yue and Lay 213) Event 1 P Event 4 dominate as Q increases. To better evaluate the tradeoff, we introduce a term r k sum, i.e., r k sum ¼ P P updip P downdip k i;jþ k i;jþ ¼ P 8 P 4 i¼1 j¼1 P 8 P 8 i¼1 j¼5 k i;jþ k i;jþ S S t (s) Fig. 7 Synthetic accelerograms from event pair Event 1 & Event 4 (Table 2). The ratio of maximum amplitude of Event 1 to that of Event 4 is Both accelerograms are calculated with Q = 5, epicentral distance = 3 km and azimuth = 252, respectively. P and S arrivals are marked by light gray arrows Table 2 Parameters of synthetic events and test stations Event No. Focal depth (km) Dip angle () Group 1 E E E Group 2 E E E Other parameters of synthetic events: strike = 22, rake = 9. Events of Group1 are shallower than the Group 2. We avoid setting events shallower than 15 km due to few moderate earthquakes happen that area. Source time function of each event is the Heaviside function in order to ensure abundant high frequency radiation ð3þ where rsum k denotes the ratio between the sum of energy k i;jþ from up-dip region to that of down-dip region. As Q increases, rsum k increases steadily (Fig. 5). When rk sum is less than the cutoff ratio of.1, the total energy received by k- station from the up-dip region is at least an order of magnitude less than that from the down-dip region. In other words, we propose that when rsum k is less than the cutoff ratio of.1, the high frequency energy from the up-dip region is unresolvable in a complicate megathrust source rupture process because seismic phases overlap in time domain. After examining of the rsum k curves of stations shown in Fig. 1a, we conclude that whether these stations along the coastline can resolve the high frequency energy radiated from up-dip region depends on the Q value in our model. Generally speaking, when Q is lower than 5, the high frequency energy radiated from the shallower parts is unresolvable using data from stations located along the east coastline of the Honshu island due to the dominance of attenuation effect, irrespective of the amount of high frequency energy radiated from the maximum slip in the updip region area. Technically, non-geometrical seismic wave attenuation is usually caused by two mechanisms, scattering due to lateral heterogeneity and intrinsic attenuation. Unlike intrinsic attenuation which is mainly associated with anelastic processes, scattering redistributes energy instead of causing energy loss per cycle (e.g., Sato et al., 213). In this study, we use wavenumber integration in a layered homogeneous media to generate Green s function without taking scattering into consideration. Therefore, for a certain value of Q, we may overestimate k i;jþ because the L2 norms of the seismic records decrease considerably after shifting energy from direct arrivals back to the coda, and so does the rsum k. Since only anelastic attenuation is considered in our model, it actually represents less attenuation effect in real media. Q -1 S generally takes the value on the order of.1 around 1 Hz and on the order of.1 at 2 Hz following the frequency-dependent relationship of Q -1 S f -n, where the exponent n ranges from. to 1. and varies regionally, with depth in the crust (e.g., Sato et al., 213; Romanowicz and Mitchell 27). Q S dominates in our model because of the amplitude of S wave or surface wave is much larger than P wave. Especially in calculating the L2 norms of the records, the k i;jþ affected by Q P is negligible. Hong (21) investigated the Lg attenuation in east Asia with both continental and oceanic environments. Their results suggest significantly low quality factor at 1 Hz of 1 or less in the Sea of Japan including our study area. As a result, attenuation of shear wave is strong in crust. The distributed cracks caused by the high seismicity of the

8 362 Earthq Sci (214) 27(4): (a) 1 1 Q=2 Q=5 E4/E1 E5/E2 E6/E Q=2 Q=5 Q=2 Q=5 spectral ratio (b) 1 1 Q=2 Q=2 Q=5 Q=5 1 1 Q=2 Q=2 Q=5 Q=5 1 1 Q=2 Q=2 Q=5 Q=5 spectral ratio (c) 1 1 Q=2 Q=2 Q=5 Q=5 1 1 Q=2 Q=2 Q=5 Q=5 1 1 Q=2 Q=2 Q=5 Q=5 spectral ratio Fig. 8 Spectral ratios of testing events from deeper parts and those from shallower parts of the megathrust fault. a Ratios of events pairs E4 & E1, E5 & E2, and E6 &E3 (Table 2.), direct distance from event to station is 2 km, spectra are calculated from 5 or 1-s-long signals. b and c are same as (a), but direct distance from event to station are 3 and 4 km. Black lines denote model with Q = 2, whereas the gray lines demote model with Q = 5 hanging wall result in highly intrinsic friction and scattering, which probably is the main reason for high attenuation. Together with the general attenuation behavior of the crust, we would expect significant low Q S (probably less than 5) in 1 to 2 Hz in this region even through the average frequency dependency parameter n is still unclear. As a result, the stations along the east coast of Honshu Island probably could not resolve the high frequency ([1 Hz) radiation from the up-dip region due to the dominance of attenuation. For further discussion, we calculate ^r sum k as the receiver is placed in a 2D region from the island to the trench, where ^r sum k ¼ lgðrk sum Þ. In the plot regions with ^rk sum lower than -1, the high frequency energy radiated from the updip region of the megathrust fault cannot be resolved (Fig. 6). As Q increases, the area with ^r sum k lower than -1 shrinks. Received high frequency energy from the shallower parts is not resolvable by analyzing data recorded by stations on Honshu island when Q \ 5. As Q increases, the high frequency energy radiated from the shallower

9 Earthq Sci (214) 27(4): parts becomes resolvable at greater regions for the receiver station to be located at, moving from the east coast toward the west coast on the island. Spectra of moderate earthquakes from various depths of the megathrust fault are usually analyzed to study the regional difference of the source radiation. (Iwakiri and Hoshiba 212) analyze the frequency-dependent trend of spectral ratio (up to 2 Hz) has been analyzed to indicate the possible different source time functions of interplate earthquakes located from the up-dip region to down-dip region. However, the main factors in affecting the wave spectra include both regional difference in the source radiation and the path effects. (e.g., Iwakiri and Hoshiba 212). (Ye et al. 213) use empirical Green s function (EGF) method to suppress propagation effects. The fairly reliable source spectra can only be extracted up to 5 Hz by carefully deconvoluting the EGF events because of the clearly increasing noise levels. In addition, the energy of surface waves could be much stronger rather than body waves for shallower events (Fig. 7). While calculating the spectral ratio of an event from the down-dip region to one from the up-dip region, surface wave could decrease the value in relative low frequency band compared to high frequency band, which could result in the biased estimation that deeper events radiate more strongly at high frequency. Here, we set six hypothetical testing events into two groups from up-dip region to down-dip region of the megathrust fault (Table 2). With a radius of 2, 3 or 4 km, we let station evenly distribute along a circle with azimuths range from 22 to 22 at a 1 interval. We then sum up the spectrums of signals at the stations for a window of 1-s-long and calculate the ratio of the spectra of event pair E4 & E1, E5 &E2, and E6 &E3. Figure 8 indicates the spectral ratios with different direct distances and Q structures in frequency band from.2 to 2 Hz. Despite the same source time function for each event pair, the spectral ratios increase considerably with frequency, indicating that lower spectral ratios in low frequency mainly result from the dramatic amplitude difference of surface wave vibrated with sources in different focal depths. Overall, while using spectral ratios to analyze the source properties, frequency-dependent trend of spectral ratio could be explained by the combination of difference in source time functions, propagation effects and the strong bias introduced by surface waves. Extracting reliable broadband source spectra (up to 2 Hz) of the interplate events requires high quality nearfault records or analysis of the high-resolution attenuation structure for the high-frequency range ([1 Hz). One of the most challenging missions of seismic observation is to overcome the difficulty of maintaining an ocean bottom seismic array from the island to the trench. According to our result (Fig. 6), hypothetical stations close to the trench are extremely sensitive to the high frequency radiated from the up-dip region. OBS (ocean bottom seismometers) data near the trench have significant advantages in analyzing the depth-varying rupture process. It could be verified by broadband OBS data in the future that whether the up-dip region area could generate as much high frequency ([1 Hz) radiation as down-dip region does. 4 Conclusions The results of our synthetic tests indicate that the possible high frequency ([1 Hz) energy radiated from the up-dip region could hardly be resolved by stations because of the high attenuation of this region, whereas the OBS stations located near the megathrust fault can detect the up-dip region generating high frequency radiation. In addition, when analyzing spectral ratios of moderate earthquakes from various depths of the megathrust fault, one should be careful about the bias introduced by the strong low frequency surface wave of earthquakes from shallower parts of the fault, which could result in erroneous interpretation that deeper events radiate more strongly at high frequency. Acknowledgments This work made use of GMT, SAC. and COMPUTER PROGRAMS IN SEISMOLOGY software. We thank Professor Z. Shen, Y. Zhou, J. Ning, Z. Ge, Haiming Zhang and two anonymous reviewers for providing thoughtful discussions and comments. This work is supported by NSFC grants (Nos and ). References Hong (21) Lg attenuation in a region with both continental and oceanic environments. Bull Seismol Soc Am. doi:1.1785/ Huang Y, Meng L, Ampuero J-P (212) A dynamic model of the frequency-dependent rupture process of the 211 Tohoku-Oki earthquake. Earth Planets Space. doi:1.547/eps Ide S, Baltay A, Beroza GC (211) Shallow dynamic overshoot and energetic deep rupture in the 211 M W 9. Tohoku-Oki Earthquake. science. doi:1.1126/science.1272 Ishii M (211) High-frequency rupture properties of the M W 9. off the Pacific coast of Tohoku Earthquake. Earth Planets Space. doi:1.547/eps Iwakiri K, Hoshiba M (212) High-frequency ([1 Hz) content of the initial fifty seconds of waveforms from the 211 Off the Pacific Coast of Tohoku Earthquake. Bull Seismol Soc Am. doi: / Kiser E, Ishii M (211) The 21 M W 8.8 Chile earthquake: triggering on multiple segments and frequency-dependent rupture behavior. Geophys Res Lett. doi:1.129/211gl4714 Koper KD, Hutko AR, Lay T (211a) Along-dip region variation of teleseismic short-period radiation from the 11 March 211

10 364 Earthq Sci (214) 27(4): Tohoku earthquake (M W 9.). Geophys Res Lett. doi:1.129/ 211GL49689 Koper KD, Hutko AR, Lay T, Ammon CJ, Kanamori H (211b) Frequency-dependent rupture process of the 211 M W 9. Tohoku Earthquake: comparison of short-period P wave backprojection images and broadband seismic rupture models. Earth Planets Space. doi:1.547/eps Lay T, Kanamori H, Ammon CJ, Koper KD, Hutko AR, Ye L, Yue H, Rushing TM (212) Depth-varying rupture properties of subduction zone megathrust faults. J Geophys Res. doi:1.129/ 211JB9133 Meng L, Ampuero JP, Luo Y, Wu W, Ni S (212) Mitigating artifacts in back-projection source imaging with implications on frequency-dependent properties of the Tohoku-Oki Earthquake. Earth Planets Space. doi:1.547/eps Romanowicz B, Mitchell BJ (27) Deep earth structure-q of the Earth from crust to core. Elsevier B.V, Amsterdam Sato H, Fehler M, Maeda T (213) Seismic wave propagation and scattering in the heterogeneous earth: second edition. Springer Press, New York Wang D, Mori J (211) Frequency-dependent energy radiation and fault coupling for the 21 M W 8.8 Maule, Chile, and 211 M W 9. Tohoku, Japan, earthquakes. Geophys Res Lett. doi:1. 129/211GL49652 Yao H, Gerstoft P, Shearer PM, Mecklenbräuker C (211) Compressive sensing of the Tohoku-Oki M W 9. earthquake: frequencydependent rupture modes. Geophys Res Lett. doi:1.129/ 211GL49223 Ye L, Lay Thorne, Kanamori Hiroo (213) Ground shaking and seismic source spectra for large earthquakes around the megathrust fault offshore of northeastern Honshu, Japan. Bull Seismol Soc Am. doi:1.1785/ Yue H, Lay T (213) Source rupture models for the M W Tohoku Earthquake from joint inversions of high-rate geodetic and seismic data. Bull Seismol Soc Am. doi:1.1785/