Migration of low frequency tremors revealed from multiple array analyses in western
|
|
- Ambrose Nicholson
- 5 years ago
- Views:
Transcription
1 Migration of low frequency tremors revealed from multiple array analyses in western Shikoku, Japan Tomotake Ueno, Takuto Maeda*, Kazushige Obara, Youichi Asano, and Tetsuya Takeda National Research Institute for Earth Science and Disaster Prevention *Now at the Center for Integrated Disaster Information Research, Interfaculty Initiative in Information Studies, The University of Tokyo Short title: MIGRATION OF TREMORS BY MULTIPLE ARRAYS Tomotake Ueno National Research Institute for Earth Science and Disaster Prevention 3-1 Tennodai, Tsukuba, Ibaraki , Japan 1
2 Abstract Multiple array observation above a belt-like tremor zone was conducted to investigate the detailed location and migration of tremor activity in western Shikoku, Japan. In March 2007, an episodic tremor and slip event occurred, and highly coherent waveforms were recorded at three arrays. Multiple signal classification analysis for the data from each array enabled measuring precise arrival directions. The majority of tremor signals suggested relatively low slowness. The arrival directions of tremor signals were used to locate tremor sources by the grid search method. Tracking the tremor activity showed that the tremor migrated within several hours in the northeast-southwest direction over a distance of km, and its migration velocity was km/h. This migration velocity is more rapid than the mean velocity of 0.5 km/h over the whole tremor episode lasting several days. Such a short timescale migration may represent fluctuation of slip acceleration during the slow slip event. Whenever a tremor migrates southwestward, very-low-frequency earthquakes occur in the vicinity of the tremor migration terminus. This indicates that the tremor migration is related to the occurrence of very-low-frequency earthquakes and slow slip events. 2
3 1. Introduction In the Nankai subduction zone in southwest Japan, low-frequency tremors have been detected [Obara, 2002] with high sensitivity seismograph networks (Hi-net) operated by the National Research Institute for Earth Science and Disaster Prevention (NIED) [Okada et al., 2004; Obara et al., 2005]. The tremors are located approximately km deep along the strike of the subducting Philippine Sea plate (PHS). In particular, in the western Shikoku district, significant tremor episodes recur at six-month intervals, and they are accompanied by slow slip events (SSEs) [Obara et al., 2004; Hirose and Obara, 2005]. Such episodic tremor and slip are also observed on the Cascadia margin in the northwestern part of the North American plate [Rogers and Dragert, 2003]. Moreover, similar phenomena have been detected in some subduction zones, including Costa Rica, Mexico, Alaska, and New Zealand [e.g., Schwartz and Rokosky, 2007]. In southwest Japan, very-low-frequency (VLF) earthquakes accompanied by tremors and SSEs have also been detected [Ito et al., 2007], and they reflect intermediate frictional properties at the deeper part of the megathrust seismogenic asperity on the subducting plate interface. During tremor and other coincident phenomena that have a characteristic period of several days, the tremor sources usually migrate along the strike of the subducting plate 3
4 with a mean velocity of approximately 10 km/day, which was inferred from the source locations by the envelope correlation method (ECM) [Obara, 2002]. A similar tremor migration velocity was estimated at the Cascadia subduction zone [Kao et al., 2006]. On the other hand, on a very short timescale, low-frequency earthquakes, which are distinct tremors composed of impulsive signals, migrated along the dip direction of the slab with a velocity of ~45 km/h during a tremor burst, as measured by waveform correlation analysis [Shelly et al., 2007]. Such a variation of the estimated velocity suggests that there may exist a scale dependency of the migration velocity associated with the duration of slow events [e.g., Ide et al., 2007; Matsuzawa et al., 2009]. The detailed migration of tremor activity is not well understood because the ECM is based on the smoothed envelope pattern obtained at Hi-net stations separated by a mean distance of approximately 25 km, and it is spatially limited for the hypocentral location of tremors. To locate tremor sources with high accuracy, dense seismic array observation is effective because wavefield analysis can be performed. In Cascadia, La Rocca et al. [2005] observed well-correlated tremor waveforms with three array observations deployed in the northern Washington and British Columbia region. La Rocca et al. [2008] pointed out that small-aperture array analysis is useful to analyze signals from a single tremor source, and 4
5 more complex and detailed analyses are required to analyze multiple-source tremors. To analyze the multiple signals in tremor wavetrains, the Multiple Signal Classification (MUSIC) method [Schmidt, 1986] is suitable, as it is for studying volcanic tremors [e.g., Saccorotti et al., 2004]. In the present study, we carried out a dense seismic array observation in western Shikoku, Japan. We located epicenters of tremor sources using arrival directions estimated from MUSIC during the episode of March From the distribution of the tremor sources, the migration of tremors and the relationship among tremors, VLF earthquakes, and SSEs are discussed. 2. Data and analysis 2-1. Dense array observation Western Shikoku is the most active tremor region and SSE area on the belt-like tremor zone in southwest Japan, and tremor and SSE episodes occur approximately every six months. To record signals from the tremor sources, dense seismic arrays were deployed at three sites in western Shikoku separated by approximately 20 km, as shown in Figure 1. The observation period was from February to May 2007 because a previous tremor episode 5
6 occurred in September 2006, and we anticipated the next tremor would occur during this observation period. As expected, a tremor and SSE episode occurred in March Each array was composed of 23 vertical-component seismometers with a natural frequency of 2 Hz and 3 three-component seismometers with a natural frequency of 1 Hz. The receiver configuration of each array is illustrated in Figure 2. The mean spacing between the receivers was approximately 30 m. The continuous waveform data were recorded at a sampling interval of 0.01 s Array processing To detect a coherent signal and arrival direction from dense seismic arrays, various methods have been developed to estimate the slowness vector [e.g., Niedel and Tarner, 1971; Rost and Thomas, 2002]. Here we adopted a two-step approach to detect the tremor signal. First, we seek a coherent wave propagating to the array as a candidate for the tremor signal. Then, the MUSIC method [Schmidt, 1986] is applied to estimate an accurate backazimuth of the propagation for each tremor signal, which is used to determine the epicentral location of the tremor. The method of tremor epicenter determination is described in the section
7 Preprocessing The tremor signal is dominant in the horizontal records compared to the vertical ones. However, as shown in Figure 3, the tremor signal is clearly coherent and the tremor amplitude is two or more times larger than the amplitude of background noise in the vertical records, perhaps because the shear-vertical component of the tremor signals might be recorded. Therefore, in our analysis, we used only vertical seismometers at each array. The continuous waveform data observed at the three arrays were processed by a bandpass filter with a frequency range of 2 8 Hz to eliminate background noise and to enhance the tremor signal with a dominant frequency of around 2 Hz [Obara and Hirose, 2006]. To detect a coherent wavefield containing the tremor signal, we evaluated the cross-correlation coefficients of the waveforms recorded at neighboring receivers for each 10-s time window every 2 s. By comparing the waveforms with and without tremor signals, the averaged cross-correlation coefficients for the tremor signals are much higher than those without tremors (Figure 3). For the processing described herein, we adopt array records whose averaged cross-correlation coefficient is higher than 0.75 as candidates for a tremor signal. 7
8 Estimation of arrival direction by the MUSIC method Tremor often occurs at some source locations simultaneously during the time period of SSE episode. In order to estimate an arrival direction of the coherent waveform from the strongest tremor source in a given time window, we used the MUSIC method [Schmidt, 1986] because this method can decompose arrival directions for multiple sources with high accuracy from a contaminated wavefield [e.g., Goldstein and Archuleta, 1987]. If we were using other traditional methods in the array analysis, the detection of the most dominant arrival direction might have been less accurate. For well-correlated wavetrains in the preprocessing step, we calculated a cross-covariance matrix for each array s data in a 10-s time window at a frequency range of 2 8 Hz. Then, eigenvalue decomposition was performed on the cross-covariance matrix. To decompose the covariance matrix into signal and noise subspaces, the MUSIC method requires a priori indication of the number of eigenvalues that correspond to signals. For the case of tremor with ambiguous arrivals, because the difference of eigenvalues between signal and background noise was small in comparison to a regular earthquake, it was difficult to determine the number of eigenvalues that correspond to tremor signals. Hence, 8
9 we have an expedient criterion as follows: the sum of selected eigenvalues from the largest signal must exceed 95% of the total sum of eigenvalues. Then, we search the MUSIC spectrum in the slowness space for the most dominant backazimuth and slowness. An example of the MUSIC spectra from each array at which the coherent signals were obtained (Figure 3b) are plotted in Figure 4. Regardless of the complexity of the tremor waveforms, clear peaks in the slowness space appear at each array. In order to evaluate the error of the MUSIC spectrum, we adapt the delete-1 Jackknife method [Tichelaar and Ruff, 1989] for the dataset, by removing a waveform recorded at a single station in each array. Examples of estimated peaks and errors of the MUSIC spectra are shown with the white circle and bars in the right panel of Figure 4. From this method, we measure the stability of backazimuth and slowness for incoming waves. In order to exclude distant signals and coherent noises originating near the surface, we eliminated the backazimuth and slowness for which the MUSIC spectrum peak is greater than 0.4 s/km. The uncertainty in the backazimuth of the MUSIC spectrum peak estimated by the delete-1 Jackknife method tends to increase with decreasing slowness. Therefore, we choose a well-determined slowness that satisfies the following threshold: the error in slowness and backazimuth are smaller than 0.2 s/km and 60, respectively, for peak 9
10 slowness of more than 0.25 s/km; and these values must be smaller than 0.4 s/km and 140, respectively, for peak slowness of less than 0.25 s/km Tremor source detection with multiple array data We determine the focus point of the three array backazimuths as the epicenter of the tremor source based on an assumption that tremors occur at the plate interface [Shelly et al., 2006; Ghosh et al., 2009]. For each time window, the epicenters of tremor sources are determined by the grid search method. The backazimuths for each grid point are evaluated using the following equations to estimate the epicentral location: Δ n obs cal obs cal φn φn ( φn φn 180) =, (1) obs cal obs cal 360 φn φn ( φn φn > 180) Δ r = Δ 2 n, (2) n Δ sd = 1 n 1 n ( Δ Δ) n 2 Δ = 1 n n Δ n, (3) 10
11 obs cal where and φ denote the observed and calculated backazimuths for the n th array, φ n n respectively. Δn denotes the difference between the observed and calculated values. Δ r, Δ, and Δ sd denote the sum of squares, the mean value, and the standard deviation of Δ n, respectively. A point where the minimum value of Δ r and Δ sd are less than 30 and 10 is treated as an epicenter of a tremor source. The reliability of the epicentral location decreases as the distance from the three arrays increases using only backazimuths. Furthermore, some sources are located at the edge of the search area. Because these sources were focused over 100 km in epicentral distance, we chose tremor sources less than 100 km from the nearest array determined by the above location process. Then, to obtain a stable estimate for tracing the tremor migration, we calculate the centroid of each tremor per hour by averaging individual epicentral locations with more than 15 events within an hour. The error of this mean location is evaluated by a bootstrap simulation [Efron and Tibshirani, 1986] using the following equations: B N B ˆ 1 = L b ( n) B n, (4) b= 1 n= 1 11
12 2 B N 1 ( ) ˆ Lb n B b= 1 N n= 1 σ = B, (5) B 1 where N is the number of data. L b (n) denotes the b th bootstrap sample consisting of a random sample selected from the detected latitudes or longitudes of tremor sources each hour. B denotes the bootstrap frequency, which is 2,000 trials in the present study. 3. Results The backazimuths of coherent signals estimated well by the MUSIC method at each array from 9 20 March 2007 are shown in Figure 5. As an indicator of epicentral locations of tremors and their migration, we adopt a statistical "mode" that is the most frequent value in the hourly backazimuth dataset within a 10 degree bin. In Figure 5, the mode is shown by a square, and frequency magnitude indicates tremor activity. Here, we detect coherent tremor signals continuing to arrive at the three arrays from March. In this period, several tremor signals frequently appeared per minute, as shown in Figure 3b, even if they were recorded by vertical component seismometers. Out of this tremor period, many cross-correlation coefficients of incoming waves were less than 0.75, as shown in 12
13 Figure 3a Characteristics of array estimation The tremors were clearly detected as coherent signals with wide ranges of backazimuths and low slowness, indicated by the dark gray region in Figure 5. Before the afternoon of 13 March, a small number of coherent signals were detected as a weak tremor. Then coherent signals were clearly detected from the afternoon of the 13th. This indicates that the active tremor started. In the active tremor period, there were two stages because of halting of tremor in the evening of the 15th; the first is from the 13th to the 15th and the second is from the 16th to the 18th. From the afternoon of the 13th to the morning of the 14th, the mode of backazimuths clearly migrated from the north (approximately 340 ) to the west (approximately 240 ) until the evening of the 14th in array C (Figure 5c). Then, the migration of backazimuths in Figure 5c repeated in this region until the evening of the 15th. In the same time period, although scattered backazimuths were obtained in array A, the mode of backazimuths changed around 200 from 240 on the morning of the 14th, and scattered between 0 and 60 from the 14th to the 15th (Figure 5a). In array B, tremor 13
14 signals are mainly distributed around a backazimuth of 0 with a very large mode size on the 13th, and the signals started to have various backazimuths between 0 and 180 from the beginning of the 14th (Figure 5b). The wide distribution of backazimuths means that the sources are located underneath the array, because backazimuth estimation is often unstable for a near-vertical incident wave. Therefore, the area of deep tremor activity migrated from the vicinity of array A to array B from the 13th to the 14th. This migration was consistent with the backazimuth change in array C. On the evening of the 15th, the tremor activity became weak for a few hours as shown in Figure 5. On the 16th, active tremor began again. After the evening of the 16th, many tremor signals were suddenly detected at the three arrays. The tremor activity was concentrated around a backazimuth of 200 in array B, whereas in arrays A and C the backazimuth also clustered around 240. This indicates that the tremor direction of each array jumped southwestward Epicentral location of tremors determined by multiple arrays The distributions of polar histograms of backazimuths are illustrated in Figure 6. The histograms indicate that many signals propagated from between arrays A and B. The 14
15 centroids of tremor activity were distributed along the strike of the PHS slab [Ueno et al., 2008], which is 60 km long and 20 km wide in the vicinity of 40 km deep in the depth contour map of PHS (Figure 7a). The tremor epicenters concentrated between arrays A and B were well located in our analysis; however other tremors, scattered far from the arrays, were less accurately located. A similar distribution was also obtained by ECM as shown in Figure 7b; however, tremor epicenters were relatively more scattered in the vicinity of the arrays. The time history of the spatial location of tremors projected on longitude and latitude axes is presented in Figure 8. The vertical bar shows the error distribution based on the bootstrap method of equation (5). The tremor activity gradually migrated southwestward from the northern part of western Shikoku during the episode. The tremor migrating southwest seems to include a repeated movement along the same route. Such a migration pattern is reliable because the error in the epicentral locations of tremor during the period from the 13th to the 16th is half as much as during the other period as shown in Figure 8. In particular, from the 13th to the 14th and the 14th to the 15th, two southwestward migrations through a similar route were clearly recognized because the migration distance exceeds the error estimation. The first migration has a 1 km/h velocity 15
16 over a distance of 12 km for 12 hours, from the evening of the 13th to the morning of the 14th. The second migration has a 2.5 km/h velocity for six hours over a distance of 15 km from the evening of the 14th to the morning of the 15th in the same region. And then, on the evening of the 16th, the tremor started to migrate steeply southwestward, far from the array stations, as shown in Figures 7 and 8. The steep migration appeared as a jump in the tremor direction (Figure 5). 4. Discussion 4-1. Scale variation of migrating tremor The tremor episode in March 2007 is characterized by a horizontal migration of the epicenters with a mean velocity of roughly 0.5 km/h estimated from the traditional ECM shown by broken lines in Figure 8b, as in the other tremor episodes in western Shikoku [Obara, 2009]. We detected a short timescale tremor migration pattern with a faster velocity of km/h over a distance of km for several hours. The sharp onset of tremor at the migration front and persisting tremor activity behind the front are interpreted as the rupture initiation of the SSE [Obara, 2009] and continuation of the slip, respectively [Obara, 2009; Maeda and Obara, 2009]. In the present study, we found that 16
17 the tremor activity migrated more than once along the same path within one episode, as shown in Figures 8. The repeating tremor migration in the short timescale might be a seismic response to the local stress change caused by the small-scale fluctuation in the spatio-temporal distribution of the slip velocity at the rupture front of the SSE and/or by the continued slip after the passage of the rupture front. Various timescales of tremor migration exist within the same region. An increased migration velocity for a short time was predicted by a scaling relationship [Ide et al., 2007]. The migration velocity in the short timescale in this study also supports the scaling relationship between characteristic duration and velocity [Matsuzawa et al., 2009]. The faster migration of the tremor burst in the same region [Shelly et al., 2007] might provide an endmember for this scaling relationship Relationship between the tremor and other coincident phenomena In March 2007, the tremor activity, VLF earthquakes, and the SSE occurred simultaneously at the western Shikoku region like in other episodes. The spatio-temporal locations of the VLF earthquakes determined by Ito et al. [2009] and Matsuzawa et al. [2009] during the SSE are indicated by diamonds in Figures 7 and 8. The VLF earthquake 17
18 locations appeared to be around the southwestern edge of the tremor activity in the early mornings of the 14th and the 15th. Comparing the spatio-temporal locations of the tremor and the VLF earthquakes, the tremor seems to migrate toward the VLF earthquake location, which is the western edge of the tremor migration shown in Figure 8a. The VLF earthquake is considered as a seismic slip with slow rupture velocity caused by the stress concentration at a VLF patch surrounded by a slow slip region [Ito et al., 2007]. Consequently, the stress was concentrated in a patch at the western edge because of the fluctuation of the SSE reflected by the repeating tremor migration on March, which triggered the VLF earthquakes at this western edge. A similar situation was observed on the 17th. On the evening of the 16th, after the tremor abruptly migrated southwestward, VLF earthquakes occurred at the end of the tremor migration. Our results suggest that the occurrence of VLF earthquakes is related to the stress concentration caused by small-scale fluctuations of SSE slips that can be detected by short timescale tremor migration. Here, we compare the tremor migration and the daily slip rate of the SSE obtained by time-dependent inversion [Hirose and Obara, 2010] as shown in Figure 9. The cumulative slip distribution, VLF earthquakes, and tremor activity that occurred in March 2007 are presented in Figure 9a. All of the events were concentrated in the northwest 18
19 portion of western Shikoku. This result was also consistent with the area exhibiting a high-energy radiation of tremors [Maeda and Obara, 2009]. Prior to the SSE occurrence and before the 13th, weak and scattered tremor activity appeared around the northern portion of array B (Figure 9b d). Once the slip began and its SSE area became large, the number of tremor sources also increased (Figure 9d f). The tremor activity appears to be located repeatedly along the slow slip contours, and it migrates toward the source area of the VLF earthquakes, as shown in Figure 9f i. On the evening of the 15th, the centroid of the tremor activity could not be obtained for Figure 8a because the tremor activity abruptly became very weak for a few hours, and then the tremor became active again as shown in Figure 5. On the evening of the 16th, the tremor activity location abruptly migrated southwestward as if by jumping, as shown in Figures 5 and 8, and the slip rate contour of the SSE also simultaneously moved to the southwestern area (Figure 9h and i). These phenomena of halting and jumping were also found in Cascadia [Kao et al., 2006, 2007]. Since the active area of the VLF earthquakes and the SSE were located far from our arrays, especially after the 17th as shown by ECM in Figure 8b, the centroid locations of tremor activity had large uncertainties, as shown in Figures 7a, 8a, and 9i. While the 19
20 locations are not well-resolved, the migration pattern is still discernible above the errors. Finally, the tremor activities scattered and gradually ceased after the conclusion of the SSE (Figures 9j l). Based on the above discussion, we conclude that the short timescale of tremor migration is a small-scale fluctuation within the overall slow slips, and it is related to the occurrence of VLF earthquakes. This type of rupture fluctuation and the relationships among detailed spatio-temporal distributions of each slow event may be keys to disclosing a series of slow events. 5. Conclusions We succeeded in detecting the migration of tremors to investigate the arrival directions of tremor signals by conducting a multiple array analysis on continuous wave data. This array experiment was the first such observation and analysis of tremor activity in the western Shikoku region. Combining the results of the array analyses, we obtained the following information regarding the March 2007 episode: 1. Highly coherent tremor signals were obtained, and tremor sources were detected with the multiple array analyses. 2. The majority of the apparent slowness values detected by the three arrays during 20
21 the active tremor period were relatively low. The incident angles of wavetrains from the tremor source to the arrays were nearly vertical. The distributions of backazimuths indicated the direction of the tremor sources. 3. A number of tremor sources were located in northwestern Shikoku. Furthermore, we obtained a highly reliable distribution of tremor epicenters for each hour. These results permitted repeating tremor migrations to be detected in a short time. The velocity of the short timescale tremor migration was approximately km/h over a distance of km for this episode. The migration velocity for the whole episode was roughly 0.5 km/h. 4. The spatio-temporal distribution of the tremor sources was correlated to the VLF earthquakes and the SSE that occurred in the same region. The short timescale migration of the tremor suggests a small-scale fluctuation within the overall slow slip that was related to the occurrences of the VLF earthquake. 21
22 Acknowledgments We thank Dr. H. Hirose for providing figures and advice regarding the short term SSE that occurred in We also thank Dr. T. Matsuzawa for his advice concerning a very-low-frequency event and the scaling relationship between characteristic duration and velocity of a migrating tremor. Dr. Y. Ito helped us prepare and deploy the three arrays and discussed the relationships among slow events with us. Dr. B. Enescu discussed our work and provided us with helpful comments. We drew all figures using Generic Mapping Tools [Wessel and Smith, 1998]. Our study was supported by all Hi-net members. Professor J. Vidale, an anonymous reviewer, and the associate editor provided useful comments to improve our manuscript. 22
23 References Efron, B., and R. Tibshirani (1986), Bootstrap Method for Standard Errors, Confidence Intervals, and Other Measures of Statistical Accuracy, Statistical Science, 1, Ghosh, A., J. E. Vidale, J. R. Sweet, K. C. Creager, and A. G. Wech (2009), Tremor patches in Cascadia revealed by seismic array analysis, Geophys. Res. Lett., 36, L17316, doi: /2009gl Goldstein, P., and R. Archuleta (1987), Array analysis of seismic signals, Geophys. Res. Lett., 14, Hirose, H., and K. Obara (2005), Repeating short- and long-term slow slip events with deep tremor activity around the Bungo channel region, southwest Japan, Earth Planets Space, 57, Hirose, H., and K. Obara (2010), Recurrence behavior of short-term slow slip and correlated non-volcanic tremor episodes in western Shikoku, southwest Japan, J. Geophys. Res., doi: /2008JB006050, in press. Ide, S., G. C. Beroza, D. R. Shelly, and T. Uchide (2007), A scaling law for slow earthquakes, Nature, 447, 76 79, doi: /nature Ito, Y., K. Obara, K. Shiomi, S. Sekine, and H. Hirose (2007), Slow earthquakes coincident 23
24 with episodic tremors and slow slip events, Science, 315, Ito, Y., K. Obara, T. Matsuzawa, and T. Maeda (2009), Very-low-frequency earthquake as small asperity on plate boundary in transition zone from locked to aseismic slip, J. Geophys. Res., 114, B00A13, doi: /2008jb Kao, H., S. Shan, G. Rogers, and H. Dragert (2007), Migration characteristics of seismic tremors in the northern Cascadia margin, Geophys. Res. Lett., 34, L03304, doi: /2006gl Kao, H., S. Shan, H. Dragert, G. Rogers, J. F. Cassidy, K. Wang, T. S. James, and K. Ramachandran (2006), Spatialtemporal patterns of seismic tremors in northern Cascadia, J. Geophys. Res., 111, B03309, doi: /2005jb La Rocca, M., D. Galluzzo, S. Malone, W. McCausland, G. Saccorotti, and E. Del Pezzo (2008), Testing Small-Aperture Array Analysis on Well-Located Earthquakes, and Application to the Location of Deep Tremor, Bull. Seism. Soc. Am., 98, , doi: / La Rocca, M., W. McCausland, D. Galluzzo, S. Malone, G. Saccorotti, and E. Del Pezzo (2005), Array measurements of deep tremor signals in the Cascadia subduction zone, Geophys. Res. Lett., 32, L21319, doi: /2005gl
25 Maeda, T., and K. Obara (2009), Spatio-temporal distribution of seismic energy radiation from low-frequency tremor in western Shikoku, Japan, J. Geophys. Res., 114, B00A09, doi: /2008jb Matsuzawa, T., K. Obara, and T. Maeda (2009), Source duration of deep very-low-frequency earthquakes in western Shikoku, Japan, J. Geophys. Res., 114, B00A11, doi: /2008jb Niedel and Tarner (1971), Semblance and other coherency measures for multichannel data, Geophysics, 36, Obara, K. (2002), Nonvolcanic deep tremor associated with subduction in southwest Japan, Science, 296, Obara, K. (2009), Phenomenology of deep slow earthquake family in southwest Japan --Spatio-temporal characteristics and segmentation--, J. Geophys. Res., submitted to this special section. Obara, K., and H. Hirose (2006), Non-volcanic deep low-frequency tremors accompanying slow slips in the southwest Japan subduction zone, Tectonophys., 417, Obara, K., H. Hirose, F. Yamamizu, and K. Kasahara (2004), Episodic slow slip events accompanied by non-volcanic tremors in southwest Japan subduction zone, Geophys. 25
26 Res. Lett., 31, L23602, doi: /2004gl Obara, K., K. Kasahara, S. Hori, and Y. Okada (2005), A densely distributed high-sensitivity seismograph network in Japan: Hi-net by National Research Institute for Earth Science and Disaster Prevention, Rev. Sci. Instrum., 76, , doi: / Okada, Y., K. Kasahara, S. Hori, K. Obara, S. Sekiguchi, H. Fujiwara, and A. Yamamoto (2004), Recent progress of seismic observation networks in Japan Hi-net, F-net, K-NET and KiK-net, Earth Planets Space, 56, xv xxviii. Rogers, G., and H. Dragert (2003), Episodic tremor and slip on the Cascadia subduction zone: The chatter of silent slip, Science, 300, Rost and Thomas (2002), Array seismology: Methods and applications, Rev. Geophys., 40 (3), 1008, doi: /2000rg Saccorotti, G., L. Zuccarello, E. Del Pezzo, J. Ibanez, and S. Gresta (2004), Quantitative analysis of the tremor wavefield at Etna Volcano, Italy, J. Volcanol. Geotherm. Res., 136, Schmidt, R. (1986), Multiple emitter location and signal parameter estimation, IEEE Trans. Antennas Propag., 34,
27 Shelly, D. R., G. C. Beroza, and S. Ide (2007), Non-volcanic tremor and low-frequency earthquake swarms, Nature, 446, , doi: /nature Shelly, D. R., G. C. Beroza, S. Ide, and S. Nakamula (2006), Low-frequency earthquakes in Shikoku, Japan and their relationship to episodic tremor and slip, Nature, 442, , doi: /nature Schwartz, Y. S., and J. M. Rokosky (2007), Slow slip events and seismic tremor at circum-pacific subduction zones, Rev. Geophys., 45, RG3004, doi: /2006rg Tichelaar, B. W., and L. J. Ruff (1989), How good are our best models? Jackknifing, Bootstrapping, and Earthquake Depth, EOS Trans. AGU, 70(20), 593. Ueno, T., T. Shibutani, and K. Ito (2008), Configuration of the continental Moho and Philippine Sea slab in Southwest Japan derived from receiver function analysis: Relation to subcrustal earthquakes, Bull. Seism. Soc. Am., 98, Wessel, P., and W. H. F. Smith (1998), New, improved version of Generic Mapping Tools released, EOS Trans. AGU, 79(47),
28 Figure captions Figure 1. Array map and distribution of tremors in southwest Japan. Triangles indicate the locations of each array, and inverse triangles denote Hi-net stations. Circles indicate the epicenters of tremors detected using ECM [Obara, 2002] from 9 20 March The depth contour map of the upper part of the PHS slab derived from receiver functions [Ueno et al., 2008] is indicated by dashed lines. Figure 2. Receiver configurations of each array. An open triangle indicates a vertical sensor. A filled triangle indicates a 3-component sensor. Figure 3. Examples of waveform records (a) without and (b) with tremor waveforms recorded by array A at 22:10 on the 13th and the 14th of March 2007 (JST), respectively, in upper panels. Lower panels illustrate the averaged correlation coefficient (A.C.), indicated by triangles. Horizontal axis represents time (s). Figure 4. Examples of waveforms recorded by arrays A (a), B (b), and C (c). Traces were recorded at the same time as the upper panel in Figure 3b. Horizontal axes represent time 28
29 (s). Right panels illustrate the MUSIC spectra for unshaded areas above the heavy bars. Vertical and horizontal axes represent slowness Px (s/km) and Py (s/km), respectively. White circles and bars show the mean peak of the normalized MUSIC spectrum and its error estimated by the delete-1 Jackknife method. Figure 5. Time series of arrival direction for arrays A (a), B (b), and C (c). Slowness values for each backazimuth are shown in gray scale. Vertical axes represent the backazimuth (degree). Horizontal axes represent time (day). Dots indicate peak values of the MUSIC spectrum. Open squares show the backazimuth mode per hour and their size depends on frequency. Figure 6. Distribution of tremor sources and polar histograms for each array. Triangles indicate locations of arrays. Gray dots denote tremor sources that were located using multiple array analysis. The sector width of the diagram is 10. The radiated axis is normalized to unity. Figure 7. Distribution of tremor activity per hour (circle colored to indicate day) of (a) the 29
30 array analysis in this study, and (b) ECM [Obara, 2002]. VLF earthquakes (diamonds denote day) and the depth contour map of the upper part of the PHS slab (dashed lines) [Ueno et al., 2008] are indicated. Black triangles denote locations of arrays. Figure 8. Changes in tremor sources for longitude (upper panel) and latitude (lower panel) of (a) the array analysis in this study and (b) ECM [Obara, 2002]. Dots indicate the centroid location of tremor activity for each hour. White diamonds denote VLF earthquakes. The error bar was estimated using the bootstrap method. The two velocities of 1 and 2.5 km/h shown in (a) indicate a migrating tremor velocity from the evening of the 13th to the morning of the 14th and from the evening of the 14th to the morning of the 15th, respectively. Broken lines in (b) show the migration velocity of 0.5 km/h. Figure 9. Mean tremor activity per hour and VLF earthquakes for the total study (a) and each day (b) (l). For the SSE, (a) illustrates the cumulative slip distribution with contours in mm and (d) (j) illustrate the daily slip rate distributions with contours in mm/y [Hirose and Obara, 2010]. Black triangles indicate array locations. Diamonds and circles denote VLF earthquakes and tremor activity, respectively. Note that the color scale is different for 30
31 (a) (days) and the others (hours). 31
32
33
34
35
36
37
38
39
40
Migration of low frequency tremors revealed from multiple array analyses in western Shikoku, Japan
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2008jb006051, 2010 Migration of low frequency tremors revealed from multiple array analyses in western Shikoku, Japan Tomotake Ueno, 1 Takuto Maeda,
More informationArray-analysis of Tremors in Shikoku Triggered by the 2012 Sumatra Earthquake
Array-analysis of Tremors in Shikoku Triggered by the 2012 Sumatra Earthquake Tianyi Li 1 Instructor: Prof. Kazushige Obara 2 1. Department of Geophysics, Peking University 2. Earthquake Research Institute,
More informationDepth-dependent slip regime on the plate interface revealed from slow earthquake activities in the Nankai subduction zone
2010/10/11-14 Earthscope Workshop Depth-dependent slip regime on the plate interface revealed from slow earthquake activities in the Nankai subduction zone Kazushige Obara, ERI, Univ. Tokyo Recurrence
More informationAlong strike variations in short term slow slip events in the southwest Japan subduction zone
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2008jb006059, 2010 Along strike variations in short term slow slip events in the southwest Japan subduction zone Shutaro Sekine, 1,2 Hitoshi Hirose,
More informationDetection and location of non-volcanic tremor beneath the Central Range in Taiwan
Detection and location of non-volcanic tremor beneath the Central Range in Taiwan Aileen Zhang Advisor: Prof. Kazushige Obara Introduction Volcanic tremor is a long-duration episode of weak seismic motions
More informationAn autocorrelation method to detect low frequency earthquakes within tremor
GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L16305, doi:10.1029/2008gl034560, 2008 An autocorrelation method to detect low frequency earthquakes within tremor Justin R. Brown, 1 Gregory C. Beroza, 1 and David
More informationAn earthquake like magnitude frequency distribution of slow slip in northern Cascadia
GEOPHYSICAL RESEARCH LETTERS, VOL. 37,, doi:10.1029/2010gl044881, 2010 An earthquake like magnitude frequency distribution of slow slip in northern Cascadia Aaron G. Wech, 1,2 Kenneth C. Creager, 1 Heidi
More informationThe Non-volcanic tremor observation in Northern Cascadia. Hsieh Hsin Sung 3/22
The Non-volcanic tremor observation in Northern Cascadia Hsieh Hsin Sung 3/22 Reference Kao, H., S. J. Shan, H. Dragert, and G. Rogers (2009), Northern Cascadia episodic tremor and slip: A decade of observations
More informationA Brownian walk model for slow earthquakes
GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L17301, doi:10.1029/2008gl034821, 2008 A Brownian walk model for slow earthquakes Satoshi Ide 1 Received 29 May 2008; revised 10 July 2008; accepted 16 July 2008;
More informationScaling relationship between the duration and the amplitude of non-volcanic deep low-frequency tremors
GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L07305, doi:10.1029/2007gl029391, 2007 Scaling relationship between the duration and the amplitude of non-volcanic deep low-frequency tremors Tomoko Watanabe, 1 Yoshihiro
More informationTremor asperities in the transition zone control evolution of slow earthquakes
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 117,, doi:10.1029/2012jb009249, 2012 Tremor asperities in the transition zone control evolution of slow earthquakes Abhijit Ghosh, 1 John E. Vidale, 2 and Kenneth
More informationTremor patches in Cascadia revealed by seismic array analysis
GEOPHYSICAL RESEARCH LETTERS, VOL. 36, L17316, doi:10.1029/2009gl039080, 2009 Tremor patches in Cascadia revealed by seismic array analysis Abhijit Ghosh, 1 John E. Vidale, 1 Justin R. Sweet, 1 Kenneth
More informationSeismic and geodetic constraints on Cascadia slow slip
1 2 3 4 5 6 Seismic and geodetic constraints on Cascadia slow slip Aaron G. Wech 1, Kenneth C. Creager 1, & Timothy I. Melbourne 2 1 University of Washington, Department of Earth and Space Science, Box
More informationvolcanic tremor and Low frequency earthquakes at mt. vesuvius M. La Rocca 1, D. Galluzzo 2 1
volcanic tremor and Low frequency earthquakes at mt. vesuvius M. La Rocca 1, D. Galluzzo 2 1 Università della Calabria, Cosenza, Italy 2 Istituto Nazionale di Geofisica e Vulcanologia Osservatorio Vesuviano,
More informationNon-volcanic tremor resulting from the combined effect of Earth tide and slow slip event
1 Non-volcanic tremor resulting from the combined effect of Earth tide and slow slip event Ryoko Nakata 1, Naoki Suda 1 & Hiroshi Tsuruoka 2 1 Department of Earth and Planetary Systems Science, Hiroshima
More informationSegmentation in episodic tremor and slip all along Cascadia
Segmentation in episodic tremor and slip all along Cascadia Michael R. Brudzinski and Richard M. Allen Geology 35 (10) 907-910, 2007, doi: 10.1130/G23740A.1 Data Repository: Methods for Automated Data
More informationAseismic slip and low-frequency earthquakes in the Bungo channel, southwestern Japan
GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L769, doi:1.19/3gl19381, Aseismic slip and low-frequency earthquakes in the Bungo channel, southwestern Japan Shinzaburo Ozawa, 1 Yuki Hatanaka, 1 Masaru Kaidzu,
More informationLETTERS. Non-volcanic tremor driven by large transient shear stresses
Vol 448 2 August 27 doi:1.138/nature617 Non-volcanic tremor driven by large transient shear stresses LETTERS Justin L. Rubinstein 1, John E. Vidale 1, Joan Gomberg 2, Paul Bodin 1, Kenneth C. Creager 1
More informationGeneration mechanism of slow earthquakes: Numerical analysis based on a dynamic model with brittle ductile mixed fault heterogeneity
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 116,, doi:10.1029/2010jb008188, 2011 Generation mechanism of slow earthquakes: Numerical analysis based on a dynamic model with brittle ductile mixed fault heterogeneity
More informationModel for triggering of non-volcanic tremor by earthquakes. Naum I. Gershenzon 1,2, Gust Bambakidis 1. Highway Dayton, OH 45435
Model for triggering of non-volcanic tremor by earthquakes Naum I. Gershenzon 1,, Gust Bambakidis 1 1 Physics Department, Wright State University, 364 Colonel Glenn Highway Dayton, OH 45435 Department
More informationDynamic Triggering Semi-Volcanic Tremor in Japanese Volcanic Region by The 2016 Mw 7.0 Kumamoto Earthquake
Dynamic Triggering Semi-Volcanic Tremor in Japanese Volcanic Region by The 016 Mw 7.0 Kumamoto Earthquake Heng-Yi Su 1 *, Aitaro Kato 1 Department of Earth Sciences, National Central University, Taoyuan
More informationMigration process of very low-frequency events based on a chain-reaction model
Migration process of very low-frequency events based on a chain-reaction model and its application to the detection of preseismic slip for megathrust earthquakes Keisuke Ariyoshi* 1, Toru Matsuzawa, Jean-Paul
More informationArray analysis and precise source location of deep tremor in Cascadia
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2008jb006041, 2010 Array analysis and precise source location of deep tremor in Cascadia Mario La Rocca, 1 Danilo Galluzzo,
More informationNonvolcanic deep tremor associated with subduction in Southwest Japan. Kazushige Obara (NIED)
Nonvolcanic deep tremor associated with subduction in Southwest Japan Kazushige Obara (NIED) One-hour record chart at the station IKTH in Shikoku Island, Southwest Japan (2001/8/17 4 a.m.) Time (min) Time
More informationLocating nonvolcanic tremors beneath the San Andreas Fault using a station pair double difference location method
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 37,, doi:10.1029/2010gl043577, 2010 Locating nonvolcanic tremors beneath the San Andreas Fault using a station pair double difference location
More informationQuantifying the time function of nonvolcanic tremor based on a stochastic model
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2009jb000829, 2010 Quantifying the time function of nonvolcanic tremor based on a stochastic model Satoshi Ide 1 Received 24 December 2009; revised
More informationG 3. AN ELECTRONIC JOURNAL OF THE EARTH SCIENCES Published by AGU and the Geochemical Society
Geosystems G 3 AN ELECTRONIC JOURNAL OF THE EARTH SCIENCES Published by AGU and the Geochemical Society Research Letter Volume 8, Number 10 30 October 2007 Q10014, doi:10.1029/2007gc001640 ISSN: 1525-2027
More informationCoseismic slip distribution of the 1946 Nankai earthquake and aseismic slips caused by the earthquake
Earth Planets Space, 53, 235 241, 2001 Coseismic slip distribution of the 1946 Nankai earthquake and aseismic slips caused by the earthquake Yuichiro Tanioka 1 and Kenji Satake 2 1 Meteorological Research
More informationModulation of Tremor Amplitudes by Tidal Stresses in Cascadia Amanda Klaus Manuscript for Master s project. Abstract
Modulation of Tremor Amplitudes by Tidal Stresses in Cascadia Amanda Klaus Manuscript for Master s project Abstract We recorded the 2010 and 2011 episodic tremor and slip events in Northern Cascadia on
More informationLETTERS. Low-frequency earthquakes in Shikoku, Japan, and their relationship to episodic tremor and slip
Vol 442 13 July 2006 doi:10.1038/nature04931 Low-frequency earthquakes in Shikoku, Japan, and their relationship to episodic tremor and slip David R. Shelly 1, Gregory C. Beroza 1, Satoshi Ide 2 & Sho
More informationJOURNAL OF GEOPHYSICAL RESEARCH, VOL. 114, XXXX, DOI: /,
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 114, XXXX, DOI:10.1029/, Short-term and long-term tremor migration patterns of the Cascadia 2004 tremor and slow slip episode using small aperture seismic arrays W
More informationTohoku University, Sendai , Japan 3 Seismological Laboratory, California Institute of Technology, Pasadena, California , USA
Earth Planets Space, 64, 693 702, 2012 Migration process of very low-frequency events based on a chain-reaction model and its application to the detection of preseismic slip for megathrust earthquakes
More informationVery low frequency earthquakes excited by the 2004 off the Kii peninsula earthquakes: A dynamic deformation process in the large accretionary prism
LETTER Earth Planets Space, 57, 321 326, 25 Very low frequency earthquakes excited by the 2 off the Kii peninsula earthquakes: A dynamic deformation process in the large accretionary prism Kazushige Obara
More informationEpisodic tremors and slip in Cascadia in the framework of the Frenkel Kontorova model
GEOPHYSICAL RESEARCH LETTERS, VOL. 38,, doi:10.1029/2010gl045225, 2011 Episodic tremors and slip in Cascadia in the framework of the Frenkel Kontorova model Naum I. Gershenzon, 1,2 Gust Bambakidis, 1 Ernest
More informationFluid related deep low-frequency earthquakes resonant with the Rayleigh waves from the 2004 Sumatra-Andaman earthquake
JOURNAL OF GEOPHYSICAL RESEARCH, VOL.???, XXXX, DOI:10.1029/, 1 2 3 Fluid related deep low-frequency earthquakes resonant with the Rayleigh waves from the 2004 Sumatra-Andaman earthquake Masatoshi Miyazawa
More informationSlip distributions of the 1944 Tonankai and 1946 Nankai earthquakes including the horizontal movement effect on tsunami generation
Slip distributions of the 1944 Tonankai and 1946 Nankai earthquakes including the horizontal movement effect on tsunami generation Toshitaka Baba Research Program for Plate Dynamics, Institute for Frontier
More informationDetection of seismic events triggered by P-waves from the 2011 Tohoku-Oki earthquake
Earth Planets Space, 64, 1223 1229, 2012 Detection of seismic events triggered by P-waves from the 2011 Tohoku-Oki earthquake Masatoshi Miyazawa 1,2 1 Disaster Prevention Research Institute, Kyoto University,
More informationEstimation of S-wave scattering coefficient in the mantle from envelope characteristics before and after the ScS arrival
GEOPHYSICAL RESEARCH LETTERS, VOL. 30, NO. 24, 2248, doi:10.1029/2003gl018413, 2003 Estimation of S-wave scattering coefficient in the mantle from envelope characteristics before and after the ScS arrival
More informationEvolution of seismic signals and slip patterns along subduction zones: insights from a friction lab scale experiment.
Evolution of seismic signals and slip patterns along subduction zones: insights from a friction lab scale experiment. Christophe Voisin, Jean-Robert Grasso, Eric Larose, François Renard * Laboratoire de
More informationRELOCATION OF THE MACHAZE AND LACERDA EARTHQUAKES IN MOZAMBIQUE AND THE RUPTURE PROCESS OF THE 2006 Mw7.0 MACHAZE EARTHQUAKE
RELOCATION OF THE MACHAZE AND LACERDA EARTHQUAKES IN MOZAMBIQUE AND THE RUPTURE PROCESS OF THE 2006 Mw7.0 MACHAZE EARTHQUAKE Paulino C. FEITIO* Supervisors: Nobuo HURUKAWA** MEE07165 Toshiaki YOKOI** ABSTRACT
More informationNon-volcanic deep low-frequency tremors accompanying slow slips in the southwest Japan subduction zone
Tectonophysics 417 (2006) 33 51 www.elsevier.com/locate/tecto Non-volcanic deep low-frequency tremors accompanying slow slips in the southwest Japan subduction zone Kazushige Obara *, Hitoshi Hirose National
More informationNon-volcanic Tremor: A Window into the Roots of Fault Zones
Non-volcanic Tremor: A Window into the Roots of Fault Zones Justin L. Rubinstein, David R. Shelly, and William L. Ellsworth Abstract The recent discovery of non-volcanic tremor in Japan and the coincidence
More informationSendai Earthquake NE Japan March 11, Some explanatory slides Bob Stern, Dave Scholl, others updated March
Sendai Earthquake NE Japan March 11, 2011 Some explanatory slides Bob Stern, Dave Scholl, others updated March 14 2011 Earth has 11 large plates and many more smaller ones. Plates are 100-200 km thick
More informationLong-period Ground Motion Characteristics of the Osaka Sedimentary Basin during the 2011 Great Tohoku Earthquake
Long-period Ground Motion Characteristics of the Osaka Sedimentary Basin during the 2011 Great Tohoku Earthquake K. Sato, K. Asano & T. Iwata Disaster Prevention Research Institute, Kyoto University, Japan
More informationUnlocking the Secrets of Slow Slip in Cascadia Using Low- Frequency Earthquakes. Justin R. Sweet. A dissertation
Unlocking the Secrets of Slow Slip in Cascadia Using Low- Frequency Earthquakes Justin R. Sweet A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy
More informationHigh resolution receiver function imaging of the seismic velocity discontinuities in the crust and the uppermost mantle beneath southwest Japan
LETTER Earth Planets Space, 55, 59 64, 2003 High resolution receiver function imaging of the seismic velocity discontinuities in the crust and the uppermost mantle beneath southwest Japan Makiko Yamauchi
More informationRemote triggering of tremor along the San Andreas Fault in central California
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 114,, doi:10.1029/2008jb006049, 2009 Remote triggering of tremor along the San Andreas Fault in central California Zhigang Peng, 1 John
More informationMoment tensor inversion of near source seismograms
Moment tensor inversion of near source seismograms Yuji Yagi and Naoki Nishimura ABSTRACT We construct a program set for estimating moment tensor solution using near source seismograms. We take the effect
More informationSlip acceleration generates seismic tremor like signals in friction experiments
Slip acceleration generates seismic tremor like signals in friction experiments Dimitri Zigone, Christophe Voisin, Eric Larose, Francois Renard, Michel Campillo To cite this version: Dimitri Zigone, Christophe
More informationModeling short and long term slow slip events in the seismic cycles of large subduction earthquakes
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2010jb007566, 2010 Modeling short and long term slow slip events in the seismic cycles of large subduction earthquakes Takanori Matsuzawa, 1 Hitoshi
More informationSpatio-temporal variation in slip rate on the plate boundary off Sanriku, northeastern Japan, estimated from small repeating earthquakes
Spatio-temporal variation in slip rate on the plate boundary off Sanriku, northeastern Japan, estimated from small repeating earthquakes T. Matsuzawa, N. Uchida, T. Igarashi *, N. Umino, and A. Hasegawa
More informationEstimation of deep fault geometry of the Nagamachi-Rifu fault from seismic array observations
Earth Planets Space,,, Estimation of deep fault geometry of the Nagamachi-Rifu fault from seismic array observations Ayako Nakamura, Youichi Asano, and Akira Hasegawa Research Center for Prediction of
More informationScaling relations of seismic moment, rupture area, average slip, and asperity size for M~9 subduction-zone earthquakes
GEOPHYSICAL RESEARCH LETTERS, VOL. 4, 7 74, doi:1.12/grl.976, 213 Scaling relations of seismic moment, rupture area, average slip, and asperity size for M~9 subduction-zone earthquakes Satoko Murotani,
More informationMoment release rate of Cascadia tremor constrained by GPS
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 114,, doi:10.1029/2008jb005909, 2009 Moment release rate of Cascadia tremor constrained by GPS Ana C. Aguiar, 1 Timothy I. Melbourne, 1
More informationCoseismic slip distribution of the 2005 off Miyagi earthquake (M7.2) estimated by inversion of teleseismic and regional seismograms
Coseismic slip distribution of the 2005 off Miyagi earthquake (M7.2) estimated by inversion of teleseismic and regional seismograms Tadashi Yaginuma 1, Tomomi Okada 1, Yuji Yagi 2, Toru Matsuzawa 1, Norihito
More informationA hierarchy of tremor migration patterns induced by the interaction of brittle asperities mediated by aseismic slip transients
A hierarchy of tremor migration patterns induced by the interaction of brittle asperities mediated by aseismic slip transients J.-P. Ampuero (Caltech Seismolab), H. Perfettini (IRD), H. Houston and B.
More informationOccurrence of quasi-periodic slow-slip off the east coast of the Boso peninsula, Central Japan
LETTER Earth Planets Space, 9, 11 1, Occurrence of quasi-periodic slow-slip off the east coast of the Boso peninsula, Central Japan Shinzaburo Ozawa, Hisashi Suito, and Mikio Tobita Geographical Survey
More informationJCR (2 ), JGR- (1 ) (4 ) 11, EPSL GRL BSSA
Dun Wang ( ) In collaboration with: Hitoshi Kawakatsu, Jim Mori, Kazuki Koketsu, Takuto Maeda, Hiroshi Tsuroka, Jiancang Zhunag, Lihua Fang, and Qiang Yao School of Geosciences, China University of Geosciences
More informationComplicated repeating earthquakes on the convergent plate boundary: Rupture processes of the 1978 and 2005 Miyagi-ken Oki earthquakes
Complicated repeating earthquakes on the convergent plate boundary: Rupture processes of the 1978 and 2005 Miyagi-ken Oki earthquakes Changjiang Wu 1 and Kazuki Koketsu Earthquake Research Institute, University
More informationDeep low-frequency tremor that correlates with passing surface waves
1 2 Deep low-frequency tremor that correlates with passing surface waves Masatoshi Miyazawa 3 Disaster Prevention Research Institute, Kyoto University, Uji, Kyoto, Japan. Emily E. Brodsky 4 Department
More informationImaging an asperity of the 2003 Tokachi-oki earthquake using a dense strong-motion seismograph network
Geophys. J. Int. (2008) 172, 1104 1116 doi: 10.1111/j.1365-246X.2007.03702.x Imaging an asperity of the 2003 Tokachi-oki earthquake using a dense strong-motion seismograph network R. Honda 1,S.Aoi 2, H.
More informationAVERAGE AND VARIATION OF FOCAL MECHANISM AROUND TOHOKU SUBDUCTION ZONE
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 24 Paper No. 414 AVERAGE AND VARIATION OF FOCAL MECHANISM AROUND TOHOKU SUBDUCTION ZONE Shunroku YAMAMOTO 1 Naohito
More informationA Prototype of Strong Ground Motion Prediction Procedure for Intraslab Earthquake based on the Characterized Source Model
A Prototype of Strong Ground Motion Prediction Procedure for Intraslab Earthquake based on the Characterized Source Model T. Iwata, K. Asano & H. Sekiguchi Disaster Prevention Research Institute, Kyoto
More informationEffects of subsurface structures of source regions on long period ground motions observed in the Tokyo Bay area, Japan
Uetake Earth, Planets and Space (27) 69:7 DOI.86/s4623-7-6-x LETTER Open Access Effects of subsurface structures of source regions on long period ground motions observed in the Tokyo Bay area, Japan Tomiichi
More informationSUPPLEMENTARY INFORMATION
doi:10.1038/nature11492 Figure S1 Short-period Seismic Energy Release Pattern Imaged by F-net. (a) Locations of broadband seismograph stations in Japanese F-net used for the 0.5-2.0 Hz P wave back-projection
More informationGEOPHYSICAL RESEARCH LETTERS, VOL. 31, L19604, doi: /2004gl020366, 2004
GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L19604, doi:10.1029/2004gl020366, 2004 Characteristic seismic activity in the subducting plate boundary zone off Kamaishi, northeastern Japan, revealed by precise
More informationFEASIBILITY STUDY ON EARTHQUAKE EARLY WARNING SYSTEM FOR THE CITY OF LIMA, PERU, USING A NEWLY DEPLOYED STRONG-MOTION NETWORK
FEASIBILITY STUDY ON EARTHQUAKE EARLY WARNING SYSTEM FOR THE CITY OF LIMA, PERU, USING A NEWLY DEPLOYED STRONG-MOTION NETWORK Cinthia CALDERON MEE1771 Supervisor: Takumi HAYASHIDA Toshiaki YOKOI ABSTRACT
More informationAn intermediate deep earthquake rupturing on a dip-bending fault: Waveform analysis of the 2003 Miyagi-ken Oki earthquake
GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L24619, doi:10.1029/2004gl021228, 2004 An intermediate deep earthquake rupturing on a dip-bending fault: Waveform analysis of the 2003 Miyagi-ken Oki earthquake Changjiang
More informationNumerical simulation of seismic cycles at a subduction zone with a laboratory-derived friction law
Numerical simulation of seismic cycles at a subduction zone with a laboratory-derived friction law Naoyuki Kato (1), Kazuro Hirahara (2) and Mikio Iizuka (3) (1) Earthquake Research Institute, University
More informationFigures S1 S4 show the measurements taken from the synthetic vespagrams where a)
Figures S1 S4 show the measurements taken from the synthetic vespagrams where a) is the differential travel time versus the Dʺ discontinuity height, b) is the differential travel time versus δv S, c) is
More informationRupture process of the 2005 West Off Fukuoka Prefecture, Japan, earthquake
LETTER Earth Planets Space, 8, 87 92, 26 Rupture process of the 2 West Off Fukuoka Prefecture, Japan, earthquake Haruo Horikawa Active Fault Research Center, National Institute of Advanced Industrial Science
More informationEffect of the Emperor seamounts on trans-oceanic propagation of the 2006 Kuril Island earthquake tsunami
GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L02611, doi:10.1029/2007gl032129, 2008 Effect of the Emperor seamounts on trans-oceanic propagation of the 2006 Kuril Island earthquake tsunami S. Koshimura, 1 Y.
More informationGeophysical Research Letters. Supporting Information for
Geophysical Research Letters Supporting Information for A P-wave based, on-site method for earthquake early warning S. Colombelli(1), A. Caruso(1), A. Zollo(1), G. Festa(1) and H. Kanamori(2) (1) Department
More informationSupporting Online Material for
www.sciencemag.org/cgi/content/full/1131692/dc1 Supporting Online Material for Localized Temporal Change of the Earth s Inner Core Boundary This PDF file includes: Materials and Methods Figs. S1 to S3
More informationEffects of Surface Geology on Seismic Motion
4 th IASPEI / IAEE International Symposium: Effects of Surface Geology on Seismic Motion August 23 26, 2011 University of California Santa Barbara ESTIMATION OF INTERSTATION GREEN S FUNCTIONS IN THE LONG-PERIOD
More informationSlow earthquakes and low frequency tremor along the Apennines, Italy
ANNALS OF GEOPHYSICS, VOL. 51, N. 2/3, April/June 2008 Slow earthquakes and low frequency tremor along the Apennines, Italy Roberto Scarpa ( 1 ) ( 2 ), Antonella Amoruso ( 3 ), Luca Crescentini ( 3 ),
More informationHigh frequency identification of non volcanic tremor triggered by regional earthquakes
GEOPHYSICAL RESEARCH LETTERS, VOL. 37,, doi:10.1029/2010gl044660, 2010 High frequency identification of non volcanic tremor triggered by regional earthquakes Aurélie Guilhem, 1 Zhigang Peng, 2 and Robert
More information29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies MODELING P WAVE MULTIPATHING IN SOUTHEAST ASIA
MODELING P WAVE MULTIPATHING IN SOUTHEAST ASIA Ali Fatehi and Keith D. Koper Saint Louis University Sponsored by the Air Force Research Laboratory ABSTRACT Contract No. FA8718-06-C-0003 We have used data
More informationHitoshi Hirose (1), and Kazuro Hirahara (2) Abstract. Introduction
Three dimensional simulation for the earthquake cycle at a subduction zone based on a rate- and state-dependent friction law: Insight into a finiteness and a variety of dip-slip earthquakes Hitoshi Hirose
More informationThe 2006 slow slip event and nonvolcanic tremor in the Mexican subduction zone
GEOPHYSICAL RESEARCH LETTERS, VOL. 37,, doi:10.1029/2010gl045424, 2010 The 2006 slow slip event and nonvolcanic tremor in the Mexican subduction zone Vladimir Kostoglodov, 1 Allen Husker, 1 Nikolai M.
More informationDouble-difference relocations of the 2004 off the Kii peninsula earthquakes
LETTER Earth Planets Space, 57, 357 362, 25 Double-difference relocations of the 24 off the Kii peninsula earthquakes Bogdan Enescu 1, James Mori 1, and Shiro Ohmi 1 1 Disaster Prevention Research Institute
More informationABSTRACT SPATIAL AND TEMPORAL PATTERNS OF NON-VOLCANIC TREMOR ALONG THE SOUTHERN CASCADIA SUBDUCTION ZONE. by Devin C. Boyarko
ABSTRACT SPATIAL AND TEMPORAL PATTERNS OF NON-VOLCANIC TREMOR ALONG THE SOUTHERN CASCADIA SUBDUCTION ZONE by Devin C. Boyarko Episodic tremor and slip (ETS), the spatial and temporal correlation of slow
More informationDeep low-frequency earthquakes near the downward extension of the seismogenic fault of the 2000 Western Tottori earthquake
Earth lanets Space, 56, 1185 1189, 2004 Deep low-frequency earthquakes near the downward extension of the seismogenic fault of the 2000 Western Tottori earthquake Shiro Ohmi, Issei Hirose, and James J.
More informationA slow slip event in the Tokai area detected by tilt and seismic observation and its possible recurrence
Earth Planets Space, 57, 917 923, 2005 A slow slip event in the Tokai area detected by tilt and seismic observation and its possible recurrence Eiji Yamamoto, Shozo Matsumura, and Tadashi Ohkubo National
More informationSource Process and Constitutive Relations of the 2011 Tohoku Earthquake Inferred from Near-Field Strong-Motion Data
Source Process and Constitutive Relations of the 2011 Tohoku Earthquake Inferred from Near-Field Strong-Motion Data Kunikazu Yoshida, Anatoly Petukhin & Ken Miyakoshi Geo-Research Institute, Japan Koji
More informationSOURCE MODELING OF RECENT LARGE INLAND CRUSTAL EARTHQUAKES IN JAPAN AND SOURCE CHARACTERIZATION FOR STRONG MOTION PREDICTION
SOURCE MODELING OF RECENT LARGE INLAND CRUSTAL EARTHQUAKES IN JAPAN AND SOURCE CHARACTERIZATION FOR STRONG MOTION PREDICTION Kimiyuki Asano 1 and Tomotaka Iwata 2 1 Assistant Professor, Disaster Prevention
More information2008 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
STRUCTURE OF THE KOREAN PENINSULA FROM WAVEFORM TRAVEL-TIME ANALYSIS Roland Gritto 1, Jacob E. Siegel 1, and Winston W. Chan 2 Array Information Technology 1 and Harris Corporation 2 Sponsored by Air Force
More informationFault Specific, Dynamic Rupture Scenarios for Strong Ground Motion Prediction
Fault Specific, Dynamic Rupture Scenarios for Strong Ground Motion Prediction H. Sekiguchi Disaster Prevention Research Institute, Kyoto University, Japan Blank Line 9 pt Y. Kase Active Fault and Earthquake
More informationNonvolcanic tremor observed in the Mexican subduction zone
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L07305, doi:10.1029/2007gl032877, 2008 Nonvolcanic tremor observed in the Mexican subduction zone Juan S. Payero, 1,2 Vladimir Kostoglodov,
More informationSource characterization of induced earthquakes by the 2011 off Tohoku, Japan, earthquake based on the strong motion simulations
Source characterization of induced earthquakes by the 2011 off Tohoku, Japan, earthquake based on the strong motion simulations K. Somei & K. Miyakoshi Geo-Reserch Institute, Osaka, Japan SUMMARY: A great
More informationoverlie the seismogenic zone offshore Costa Rica, making the margin particularly well suited for combined land and ocean geophysical studies (Figure
Chapter 1 Introduction Historically, highly destructive large magnitude (M w >7.0) underthrusting earthquakes nucleate along the shallow segment of subduction zone megathrust fault, and this region of
More informationFracture induced shear wave splitting in a source area of triggered seismicity by the Tohoku-oki earthquake in northeastern Japan.
Fracture induced shear wave splitting in a source area of triggered seismicity by the Tohoku-oki earthquake in northeastern Japan Masahiro Kosuga 1 1. Corresponding Author. Professor, Graduate School of
More informationScattered waves from low-frequency earthquakes and plate boundary structure in northern Cascadia
GEOPHYSICAL RESEARCH LETTERS, VOL. 40, 4238 4243, doi:10.1002/grl.50826, 2013 Scattered waves from low-frequency earthquakes and plate boundary structure in northern Cascadia Robert L. Nowack 1 and Michael
More informationNegative repeating doublets in an aftershock sequence
LETTER Earth Planets Space, 65, 923 927, 2013 Negative repeating doublets in an aftershock sequence X. J. Ma and Z. L. Wu Institute of Geophysics, China Earthquake Administration, 100081 Beijing, China
More informationGROUND MOTION SPECTRAL INTENSITY PREDICTION WITH STOCHASTIC GREEN S FUNCTION METHOD FOR HYPOTHETICAL GREAT EARTHQUAKES ALONG THE NANKAI TROUGH, JAPAN
GROUND MOTION SPECTRAL INTENSITY PREDICTION WITH STOCHASTIC GREEN S FUNCTION METHOD FOR HYPOTHETICAL GREAT EARTHQUAKES ALONG THE NANKAI TROUGH, JAPAN Masayuki YOSHIMI 1, Yasuto KUWAHARA 2, Masayuki YAMADA
More informationEffects of Surface Geology on Seismic Motion
4 th IASPEI / IAEE International Symposium: Effects of Surface Geology on Seismic Motion August 23 26, 2011 University of California Santa Barbara TESTING THREE-DIMENSIONAL BASIN STRUCTURE MODEL OF THE
More informationGEOPHYSICAL RESEARCH LETTERS, VOL. 35, L14308, doi: /2008gl034461, 2008
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35,, doi:10.1029/2008gl034461, 2008 Tomographic evidence for hydrated oceanic crust of the Pacific slab beneath northeastern Japan: Implications
More informationStudy megathrust creep to understand megathrust earthquakes
1 Study megathrust creep to understand megathrust earthquakes Kelin Wang Pacific Geoscience Centre, Geological Survey of Canada, kelin.wang@canada.ca Introduction Once upon a time, there was a belief that
More informationRepeating short- and long-term slow slip events with deep tremor activity around the Bungo channel region, southwest Japan
Earth Planets Space, 57, 961 972, 5 Repeating short- and long-term slow slip events with deep tremor activity around the Bungo channel region, southwest Japan Hitoshi Hirose and Kazushige Obara National
More informationCorrelation of tremor activity with tidal stress in the northern. Cascadia subduction zone
Correlation of tremor activity with tidal stress in the northern Cascadia subduction zone Anthony Lambert 1, Honn Kao 1, Garry Rogers 1 and Nicholas Courtier 1 1 Geological Survey of Canada, Pacific Geoscience
More information