Triggering of earthquakes during the 2000 Papua New Guinea earthquake sequence

Size: px
Start display at page:

Download "Triggering of earthquakes during the 2000 Papua New Guinea earthquake sequence"

Transcription

1 JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 112,, doi: /2006jb004480, 2007 Triggering of earthquakes during the 2000 Papua New Guinea earthquake sequence Sun-Cheon Park 1 and Jim Mori 1 Received 3 May 2006; revised 20 September 2006; accepted 17 October 2006; published 8 March [1] A sequence of large earthquakes occurred in the New Britain/New Ireland region of Papua New Guinea in late The sequence started with a M w 6.8 earthquake along the New Britain Trench on 29 October. About 20 days later a M w 8.2 earthquake occurred in the New Ireland region on 16 November and produced large strike-slip surface displacements and tsunamis. Following the M w 8.2 event, two large earthquakes (M w 7.5) occurred along the nearby New Britain Trench on 16 and 17 November. Furthermore, small triggered events were observed over a wide area outside of the rupture zones of the large earthquakes, with different focal mechanisms from those of the major events. There is likely some mechanism(s) that triggered this remarkable sequence, and we document the details of the spatial and temporal patterns of the events. We investigated if static stress changes can explain the initiation of the large earthquakes and also some groups of triggered smaller events. There are mixed results. The occurrences of the two M 7.5 major events may be explained by the static stress changes; however, there are also some earthquakes that are not consistent with triggering by static stress changes. There may be multiple mechanisms, including static and dynamic stress changes, that are needed to explain the complicated sequence of earthquakes. Citation: Park, S.-C., and J. Mori (2007), Triggering of earthquakes during the 2000 Papua New Guinea earthquake sequence, J. Geophys. Res., 112,, doi: /2006jb Introduction [2] Papua New Guinea is one of the most seismically active regions in the world. Earthquakes, including many with magnitudes greater than 7, occur along the New Britain Trench, where the Solomon Sea Plate is subducting beneath the South Bismarck and Pacific plates (Figure 1). There is a triple junction south of New Ireland with the Bismarck Sea Seismic Lineation, which is a transform boundary between the South Bismarck Plate and the Pacific Plate (or North Bismarck Plate). In the studied area bounded by 150 E to 156 E and 2 S to 7 S, 17 earthquakes with magnitude greater than 7 occurred since 1970, including the remarkable sequence in late 2000 that included M w 8.2, M w 7.5, and M w 7.4 events. The latter two events are not aftershocks in the standard sense since they occurred well outside of the rupture area of the M w 8.2 earthquake in a different tectonic setting with different focal mechanisms. Stars in Figure 2 indicate the locations of the large earthquakes from the International Seismological Centre (ISC) and National Earthquake Information Center (NEIC) of the U.S. Geological Survey (USGS) catalogues. [3] This sequence is one of the clearest examples of the apparent triggering of very large earthquakes and we present details of the timing and spatial extents of the ruptures for 1 Disaster Prevention Research Institute, Kyoto University, Kyoto, Japan. Copyright 2007 by the American Geophysical Union /07/2006JB004480$09.00 the major events so that the triggering mechanisms can be evaluated. Over the past decade, many studies have addressed the occurrence of aftershocks and triggered earthquakes and there has been significant progress in clarifying some of the explanations for the induced seismicity [Steacy et al., 2005]. Mechanisms have included static stress triggering [e.g., King et al., 1994; Toda et al., 1998; Díez et al., 2005], dynamic triggering [e.g., Gomberg, 2001; Gomberg et al., 2003], rate-and-state friction [e.g., Dieterich, 1994; Parsons, 2005], and the effects of fluid flow [e.g., Gavrilenko, 2005]. Among these proposed mechanisms and processes, static stress triggering is one common way to evaluate the influence of large earthquakes on the surrounding area. Therefore we investigate the effects of static stress changes from the large earthquakes and evaluate their effectiveness in explaining the triggering of following earthquakes in the sequence. Geist and Parsons [2005] also studied the static stress triggering from the M w November event on the following M w 7.5 event. Using the fault models derived in this study, we extend this analysis to the 17 November event and other smaller earthquakes that are triggered over a wide region. Our waveform modeling of the larger events produces good constraints on the depths and faulting geometries. [4] A strong sequence of earthquakes occurred in Papua New Guinea in late 2000 (Figure 3). A M w 6.8 earthquake started the sequence along the New Britain Trench at 0837:08 UTC on 29 October Normal aftershocks followed this event and there were also earthquakes during the following 20 days in the New Ireland and New Britain 1of14

2 Figure 1. Tectonic setting of Papua New Guinea after Tregoning et al. [1998]. In the studied area, indicated by the dotted box, the Solomon Sea plate is subducting beneath the South Bismarck and Pacific plates. There is a triple junction with a transform fault, Bismarck Sea Seismic Lineation, south of New Ireland. region, which were few hundred kilometers from the 29 October event and near the locations of subsequent large earthquakes (Figure 3a). [5] At 0454:56 UTC on 16 November 2000, a great earthquake (M w 8.2) occurred in the New Ireland region. The epicenter was located in the St. George Channel between New Ireland and New Britain Islands. The aftershock distribution crosses southern New Ireland with locations that are coincident with the Weitin-Kamdaru fault (Figure 3b). There were observed strike-slip surface displacements of up to 5 m along the fault on the east side of New Ireland (I. Itikarai, personal communication, 2001). GPS observation also showed left-lateral movement up to 8 m along the Weitin-Kamdaru fault [Tregoning et al., 2005]. Tsunamis with heights up to 3 m were reported along the coasts of New Ireland, New Britain, and Bougainville (ISC catalogue). We name this earthquake the 16 November A event. [6] Another large earthquake (M w 7.5) occurred about 3 hours after the New Ireland event at 0742:16 UTC on 16 November along the New Britain Trench. We name this earthquake the 16 November B event. The following seismicity (Figure 3c) stretches about 100 km south of the rupture areas of these two major earthquakes. [7] A third large earthquake (M w 7.4) occurred along the New Britain Trench about 40 hours after the New Ireland event at 2101:56 UTC on 17 November. Figure 3d shows that following these major earthquakes, there are events in the Bismarck Sea west of New Ireland and continued seismicity south of the 16 November B event, which are outside of the immediate aftershock zones. These large earthquakes appear to strongly affect the stress field in this region and trigger many events over a wide area. Increased seismicity is noticeable here compared to the background seismicity (Figure 2) for 6 months before the sequence. Furthermore, the occurrences of the major earthquakes are probably related to the previously occurring seismicity. 2of14

3 Figure 2. Large earthquakes and background seismicity in studied area. Stars indicate large earthquakes (M > 7) since 1970: 1, 14 July 1971 shallow M7.9; 2, 26 July 1971 shallow M7.9; 3, 20 July 1975 shallow M7.9; 4, 18 March 1983 intermediate M7.9; 5, 10 May 1985 shallow M7.3; 6, 3 July 1985 shallow M7.4; 7, 30 December 1990 intermediate M7.5; 8, 28 September 1991 shallow M7.1; 9, 16 August 1995 shallow M7.8; 10, 16 August 1995 shallow M7.2; 11, 29 April 1996 shallow M7.5; 12, 10 May 1999 intermediate M7.1; 13, 16 May 1999 intermediate M7.1; 14, 16 November 2000 shallow M8.2; 15, 16 November 2000 shallow M7.8; 16, 17 November 2000 shallow M8.2; 17, 9 September 2005 intermediate M7.7. Circles and triangles indicate shallow seismicity from 1 May 2000 before the 29 October event for about 6 months. Earthquake data are from ISC and USGS catalogues. [8] Table 1 lists the source parameters for the studied events taken from the USGS catalogue. Hereafter event identifications shown in Table 1 will be used in the text. 2. Slip Distributions 2.1. Data and Methods [9] To determine the slip distributions for the three major earthquakes, we carried out waveform inversions using teleseismic P waveforms. We obtained digital P wave data for stations at distance ranges of about 30 and 90 from the Incorporated Research Institutions for Seismology (IRIS) data center. Original data were transformed to displacement using the station instrument files provided by IRIS, and high-passed filtered at 0.01 Hz. We used as many stations as possible (12 to 15) for the analyses; however, the azimuthal coverage is not optimum because there are far fewer station toward the south and east. Figure 4 shows the station distributions used for teleseismic waveform inversions for the three earthquakes. [10] To estimate the slip distributions, we used a multiple time window inversion [Hartzell and Heaton, 1983]. In this method, the fault is divided into smaller subfaults and Green s functions are calculated for the propagation response from each subfault to each station. We calculated Green s functions using a program by Langston and Helmberger [1975]. A least squares inversion was then used to determine the slip on each subfault by minimizing the fit between calculated waveforms and the observed data. Green s functions for slip both parallel and perpendicular to the strike of the fault are used, so that the slip direction of each subfault is determined. A smoothing matrix was added to the inversion to increase the stability of the solutions spatially and temporally. After testing various values, we found values of the smoothing that give a reasonable trade-off between the smoothness of the solution and a good spatial resolution The 16 November A Event [11] For setting up the inversion, the fault plane was divided into 20 4 subfaults. It has a length of 200 km and a downdip width of 42.6 km that reaches a depth of 40 km from the surface. The hypocenter is located at the northwestern end of the fault at a depth of 15 km. By trial and error we tried numerous models with a range of strike and dip angles for the fault plane, including multiple seg- 3of14

4 Figure 3. The 2000 earthquake sequence in Papua New Guinea. (a) Time period from the M w 6.8 event on 29 October 2000 up to the time immediately before the 16 November A event. (b) Time from the 16 November A event (M w 8.2) at 0454:56 UTC up to the time prior to the 16 November B event. (c) Time from the 16 November B event (M w 7.5) at 0742:16 UTC up to the time prior to the 17 November event. (d) Time from the 17 November event (M w 7.4) at 2101:56 UTC until 30 November. Open stars and solid rectangles indicate the epicenters and fault planes of the large earthquakes, respectively, occurred in the time interval. Gray solid stars and dotted rectangles indicate the previously occurring epicenters and fault planes. Circles and triangles are used to distinguish the hypocentral depth with gray for earthquakes occurred in the previous time interval. Sizes of the symbols are dependent on magnitude. Earthquake data are from the USGS catalogue. ments with different strikes, and found that a model with strike of 320 and dip of 70 provided the best fit to the data. Figure 5 shows the slip distribution (Figure 5a) and waveform fits (Figure 5b) for this event. The vectors show the direction and amount of displacement of the hanging wall at each subfault in Figure 5a. In Figure 5b, solid lines are observed waveforms and dotted lines are the synthetic seismograms. Synthetics in northeastern stations do not fit well the observed waveforms; however, they reproduce the overall amplitudes. The average value of rakes for all subfaults was 58.3 (Figure 5c). The source time function (Figure 5d) shows a long moment release with a total duration of 88 s. We tested various values of a fixed rupture velocity, which ranged from 2.0 to 3.0 km/s, and found the best solution to be 2.4 km/s. The largest slips are about 10 to 11.5 m. The total moment obtained is dyn cm, which is equivalent to M w 8.2. Table 1. Source Parameters of Earthquakes Studied a Event Origin Time, UTC Latitude Longitude Depth, km Magnitude M w 29 Oct 0837: Nov A 0454: Nov B 0742: Nov 2101: a Data from USGS. 4of14

5 Figure 4. Station distribution used for the waveform inversions for the (a) 16 November A, (b) 16 November B, and (c) 17 November events. [12] This event shows a complicated slip distribution having both strike-slip and dip-slip components. Strike-slip motions are seen in the upper portion of the fault about 40 to 110 km southeast of the hypocenter, consistent with the mostly strike-slip surface fault displacements seen on New Ireland. Especially those of the uppermost subfaults around 100 km from the hypocenter correspond to the large (about 5 m) surface displacement that was observed along the Weitin-Kamdaru fault on the southeastern coast of New Ireland (I. Itikarai, personal communication, 2001). Farther to the southeast, the slips gradually rotate to dip slip. The large shallow dip-slip movements on the southeastern 5of14

6 Figure 5. Slip distribution of the 16 November A event. (a) Slip distribution. Star indicates the hypocenter, and vectors show the direction and amount of displacement. Slip motions gradually change from strike slip to dip slip for subfaults from the northwest to southeast. (b) Fit of synthetic (dotted) and observed (solid) waveforms. Scales are the same for all stations. (c) Focal mechanism; strike 320, dip 70, average rake 58. (d) Source time function (horizontal scale in s, vertical scale in dyn cm/s). portion may provide an explanation for the observed tsunamis on the southeast coast of New Ireland The 16 November B Event [13] For this event and the 17 November event, we assumed a fault plane that is consistent with shallow thrust events on the New Britain Trench. Almost all moderate and large earthquakes located near the subduction interface with this typical mechanism are observed to occur on the shallowly dipping plane. Also relocated aftershocks by Tregoning et al. [2005] showed that both the 16 November B and 17 November events were plate boundary thrust events on the New Britain subduction zone defined by 20 dips on the Wadati-Benioff zones. [14] The fault plane of this event was divided into 10 7 subfaults. Figure 6 shows the slip distribution (Figure 6a) and waveforms (Figure 6b). The length and width of the fault plane that we inverted for are 100 km and 95.8 km, respectively, and the fault reaches 36 km from 8 km depth. The hypocenter is located close to the center of the fault at a depth of 26 km. These fault dimensions are determined by trial and errors so that the area is large enough that there are 6of14

7 Figure 6. Slip distribution of the 16 November B event. (a) Slip distribution. Star indicates the hypocenter. (b) Fit of synthetic (dotted) and observed (solid) waveforms. Scales are the same for all stations. (c) Focal mechanism; strike 240, dip 17, average rake 70. (d) Source time function (horizontal scale in s, vertical scale in dyn cm/s). no large values of slip on the edge, yet small enough so that there are not large areas of zero slip which can give unstable solutions. The strike and dip angles were determined as 240 and 17 and the average value of rakes for all subfaults was 69.8 (Figure 6c). According to the source time function (Figure 6d), seismic moment was released mainly around the hypocenter for totally 35 s. We found the best solution to be for a rupture velocity of 2.7 km/s. [15] The slip distribution is smaller and simpler than for the 16 November A event. There is an area of large slip very close to the hypocenter and the largest slip is about 8.5 m. 7of14

8 Figure 7. Slip distribution of the 17 November event. (a) Slip distribution. Star indicates the hypocenter. (b) Fit of synthetic (dotted) and observed (solid) waveforms. Scales are the same for all stations. (c) Focal mechanism; strike 230, dip 25, average rake 57. (d) Source time function (horizontal scale in s, vertical scale in dyn cm/s). The total moment obtained is dyn cm, which is equivalent to M w The 17 November Event [16] We divided the fault plane into 10 7 subfaults for the inversion. Figure 7 shows the slip distribution (Figure 7a) and waveforms (Figure 7b). The inverted fault plane has a length of 100 km and width of 82.8 km. It is located deeper than the 16 November B event and extends from a depth of 20 km to 55 km. The hypocenter is located southeast from the center of the fault and at a depth of 36 km. The strike and dip angles were determined as 230 8of14

9 Table 2. Focal Mechanisms and Magnitudes of the Three Major Earthquakes Obtained in This Study Comparing to Those From USGS MT and Harvard CMT (HRV) Solutions Source (Strike, Dip, Rake) a Depth, km M w 16 November A Event This study (320, 70, 58) USGS (180, 14, 143), (306, 82, 79) HRV (328, 43, 3), (236, 88, 133) November B Event This study (240, 17, 70) USGS (288, 24, 122), (73, 70, 76) HRV (253, 15, 93), (70, 75, 89) November Event This study (230, 25, 57) USGS (281, 32, 94), (97, 58, 88) HRV (230, 24, 64), (78, 68, 101) a Rake angles of this study are average values. and 25, respectively, and the average value of rakes was 56.5 (Figure 7c). The source time function (Figure 7d) has mainly one pulse and its total duration is 49 s. The best solution was found to be for a rupture velocity of 2.5 km/s. [17] The slip distribution is fairly simple showing an area of large slip including the hypocenter. The largest slips are about 4 m. The total moment obtained is dyn cm, which is equivalent to M w Comparison With Other Studies [18] We compare our results from the teleseismic waveform inversions with those from USGS moment tensor (MT) and Harvard centroid moment tensor (CMT) solutions [Dziewonski et al., 1981] in Table 2. Focal mechanisms obtained in this study are more similar to those of the Harvard CMT than the USGS MT solutions, although our moment magnitudes (M w ) are smaller than those of the Harvard solutions except for the 16 November A event. [19] The Earthquake Research Institute, University of Tokyo, also provided body waveform inversions of the 16 November A event (EIC note 94, 2000). Their results showed that the earthquake is left-lateral strike-slip event with focal mechanism of (strike, dip, rake) = (143, 85, 2) and moment magnitude (M w ) is 8.0. Their slip distribution has the largest slip at about 120 km to the southeast from the hypocenter. This slip distribution is somewhat different from our result which has the largest slip at about 80 km far, probably due to the difference of preferred rupture velocity. They also cited the possibility of a dip-slip subevent about 2 min after the initial rupture, which may be consistent with the dip-slip motion at southeast portion of our slip distribution. [20] Rham and Das [2003] relocated aftershocks of the 16 November A earthquake, including the 16 November B and 17 November events. They noted that the 16 November A event is a shallow left-lateral strike-slip event with a fault plane having strike 320 and vertical dip, and the rupture propagated unilaterally from northwest to southeast for about 120 km. This is generally consistent with our results if we consider only the strike-slip portions of our model. They also reported that most of the larger aftershocks (M w > 6) including the 16 November B and 17 November events have thrust mechanisms, which are consistent with our values. 3. Static Stress Changes [21] To evaluate the effects of static stress triggering, we calculated the Coulomb failure stress changes (Ds f ) that resulted from the larger earthquakes (the 29 October and 16 November A and B) in the region of the following earthquakes. This method for calculating a static failure potential has been used recently for many earthquakes [e.g., King et al., 1994; Stein, 1999]: Ds f ¼ Dt þ m 0 Ds n where Dt is the change in shear stress, Ds n is the change in normal stress, and m 0 is the value of the apparent coefficient of friction. Using the slip distributions obtained in section 2, we can calculate the shear and normal stress changes for all locations in the surrounding area. To calculate these stresses, we used the program written by Okada [1992] and calculated the Coulomb failure stress changes using equation (1). Positive values of the Coulomb failure stress changes (Ds f ) mean that static stress increases and brings the specified fault closer to failure. On the other hand, the negative values mean that static stress decreases and tends to inhibit failure on the fault. [22] We calculated the Coulomb failure stress changes at the hypocenter of the triggered event because once a rupture begins, the dynamic stress changes associated with the rupture front are much larger than the static stress changes. In this point of view, triggering of both large and small earthquakes may be treated similarly, that is, we may not need to distinguish different effects of static stress triggering depending on earthquake size. The mechanism of rupture growth to a large earthquake after initiation is another problem not addressed in this paper. [23] We first calculated the effects of the static stress changes from the large earthquakes on the following large earthquakes, starting with the M w October event and continuing for the 16 November A, 16 November B, and 17 November events. We also estimated the Coulomb failure stress changes from the 16 November A and B events in the region of the smaller triggered earthquakes. Triggered earthquakes are identified as seismicity outside of the immediate rupture areas of the major events having different focal mechanisms. Figure 8 shows the focal mechanisms of the major earthquakes, their aftershocks and three groups of triggered events that occurred between 16 and 30 November. Solid rectangles are the horizontal projections of the fault planes of the three major earthquakes obtained in section 3. For aftershocks and triggered events, we used the hypocenters determined by USGS and focal mechanisms both of the USGS and Harvard solutions Triggering of Major Earthquakes [24] For the fault geometry of the 29 October event (M w 6.8), we used (strike, dip, rake) = (322, 63, 90) obtained from the Harvard CMT solutions and focal depth of 50 km from the USGS catalogue. Then we calculated the regional static stress changes caused by this earthquake. The static stress perturbation patterns showed stress ð1þ 9of14

10 Figure 8. Focal mechanisms of the large studied earthquakes, aftershocks, and triggered events occurred until 30 November. Stars indicate the epicenters of the large earthquakes, and rectangles are horizontal projection of fault planes of the three major earthquakes. Focal mechanisms of the three major earthquakes are obtained in this study, and those of the 29 October event, aftershocks and triggered earthquakes are from Harvard CMT solutions. Groups 1 to 3, which are assumed to be triggered earthquakes, are identified by focal mechanisms, which are different from those of the three major earthquakes. Earthquake parameters of events 1 to 6 are given in Table 4. increases around the hypocenters of the 16 November A and 17 November events and decreases around the 16 November B event, but all the changes were very small (order of 10 4 MPa). For a plate velocity of 15 cm/yr [Tregoning et al., 1998], the accumulating stress is about 0.3 MPa/yr, so that the static stress changes are equivalent to the tectonic stress that accumulates in just a few hours. Therefore static stress changes from the 29 October event can have little effect in the regions of the following large earthquakes. [25] We next estimated the Coulomb failure stress changes from the larger events in the sequence. Figure 9 shows the effects from the M w November A event for the region of the 16 November B event. The Coulomb failure stress changes are calculated for a depth of 26 km and m 0 = 0.6. We tested various values for the apparent coefficient of friction from 0.1 to 0.7. The stress changes patterns did not show significant differences from that for m 0 = 0.6, and the values of the stress increases around the hypocenter were nearly the same. Note that the hypocenter and almost all of the rupture area of the 16 November B event are located within a positive area of Figure 9. This result indicates that the event occurred in an area of static stress increases, promoting to faulting. The value of the stress increases around the hypocenter is about 0.29 MPa. We also calculated the static stress changes at different depths from 10 to 34 km for average rake angles and found that most of the slip area falls within a region with fairly large stress changes (more than about 0.1 MPa). [26] Figure 10 shows the combined Coulomb failure stress changes from the 16 November A and B events, estimated for the fault plane orientation of the 17 November event. These calculations are done for a depth of 36 km and m 0 = 0.6. The summed values of the Coulomb failure stress changes around the hypocenter of the 17 November event are about 0.03 MPa. There is an increase of the static stress around the hypocenter produced mainly by the 16 November A event and the effect of the 16 November B event, which is negative, is relatively small. Table 3 summarizes the results of the static stress changes for the major earthquakes. 10 of 14

11 Figure 9. Static stress changes from the 16 November A event for the fault plane of the 16 November B event. Star indicates the epicenter of the 16 November B event, and the rectangle is the fault plane of the 16 November A event. Static stress changes show strong increase around the hypocenter of the 16 November B event Triggering of Smaller Earthquakes [27] If static stress changes are effective in triggering the large earthquakes, they should also trigger smaller earthquakes. Using hypocentral data from the USGS, we identified three groups of triggered events that occurred outside of the regions of the fault planes of the three major earthquakes from 16 November through 30 November, as shown in Figure 8. Group 1 had 9 events (compared to 7 events that occurred in the same region over the past year) with oblique strike-slip mechanisms that were located along the Bismarck Sea Seismic Lineation (Figure 1). Group 2 had 70 events (4 events over the past year) with normal fault mechanisms that are outer rise events along the New Britain Trench. Group 3 had 37 events (2 events over the past year) Figure 10. Combined static stress changes from the 16 November A and B events for the fault plane of the 17 November event. Star indicates the epicenter of the 17 November event, and the rectangles are the fault planes of the 16 November A and B events. Static stress changes show increases around the hypocenter of the 17 November event. 11 of 14

12 Table 3. Summary of Static Stress Change Estimates for the Major Earthquakes a Receiver Plane Source 29 Oct 16 Nov A 16 Nov B 16 Nov A and B Dominant Event 16 Nov A + 16 Nov B + 17 Nov Nov A a A plus indicates that the Ds f was conducive to faulting while a minus indicates that Ds f was inconsistent with the faulting. See the text for the details. that occur within the New Britain subduction zone with strike-slip mechanisms. We calculated the Coulomb failure stress changes from the 16 November A and B events for the appropriate focal mechanism in the region of each group, and the results with source parameters are listed in Table 4. The static stress changes from the 16 November B event are relatively small and do not affect the results for the groups 1 and 2. We tested for both nodal planes of the USGS and Harvard focal mechanisms. [28] As shown in Table 4, earthquakes in the Bismarck Sea in group 1 have positive values of the Coulomb failure stress changes from both of the 16 November A and B events. Outer rise earthquakes in group 2, which are normal fault events, experienced static stress decreases from the 16 November A and B events, though static stress changes from the 16 November B event can be either increases or decreases depending on the details of the assumed mechanism. For the strike-slip events in group 3, a static stress decrease was observed. However, the hypocenter is located at almost the boundary between static stress increase Table 4. Results of Static Stress Change Estimates for Smaller Earthquakes a Date Origin Time, UTC Location (Latitude, Longitude) Depth, Focal km M w Mechanism Group 1, Event 1 21 Nov : , (26, 77, 177) A (117, 87, 13) A (4, 58, 163) A (103,76, 33) A (26, 77, 177) B (117, 87, 13) B (4, 58, 163) B (103, 76, 33) B Group 1, Event 2 23 Nov : , (2, 42, 166) A (102, 81, 49) A (2, 42, 166) B (102, 81, 49) B Group 2, Event 3 17 Nov : , (139, 38, 71) A (296,54, 104) A (130,44, 106) A (333,48, 75) A (139,38, 71) B (296,54, 104) B (130,44, 106) B (333,48, 75) B Group 2, Event 4 22 Nov : , (183, 53, 73) A (337, 41, 111) A (137, 48, 89) A (315, 42, 91) A (183, 53, 73) B (337, 41, 111) B (137, 48, 89) B (315, 42, 91) B Group 2, Event 5 26 Nov : , (134, 45, 85) A (307, 45, 95) A (134, 45, 85) B (307, 45, 95) B Group 3, Event 6 26 Nov : , (171, 66, 2) A (262, 88, 156) A (171, 66, 2) B (262, 88, 156) B a Origin time, location, depth, and M w are from USGS catalogue, and those from Harvard solutions are used when they were not determined by USGS. Focal mechanisms tested are two failure planes from USGS and Harvard solutions, but only Harvard ones are given when they were unavailable in USGS solutions. CFF denotes calculated Coulomb failure stress changes. A and B in source event column denotes the 16 November A and B events, respectively. CFF, MPa Source Event 12 of 14

13 and decrease. Therefore the result will strongly depend on the details of the hypocentral determinations and focal mechanisms Summary of Static Stress Changes [29] We summarize the results of our analyses of static stress changes and triggered earthquakes as follows: (1) Effect of the M w October event in the region of the following major earthquakes was too small to be significant, (2) effect of the M w November A event had a positive stress change of about 0.29 MPa in the region of the M w November B event (similar to conclusion of Geist and Parsons [2005]), (3) combined effects of the 16 November A and B events had a positive stress change of about 0.03 MPa in the region of the M w November event, (4) combined effects of the 16 November A and B events had a positive stress change of about to 0.05 MPa for the strike-slip events in the Bismarck Sea group, (5) combined effects of the 16 November A and B events have a negative stress change of about MPa for the outer rise normal-faulting events, and (6) combined effects of the 16 November A and B events have a mixed stress change, depending on the faulting geometries for the strike-slip events in subduction zone. 4. Discussion [30] For the five cases of triggered events summarized above, where the static stress changes are large enough to potentially be the triggering mechanism, three examples have positive results. With this mixed result and small sample, it is difficult to make a definitive statement about the validity of static stress triggering. Taking into account the fact that for this sequence the areas of the positive triggering are larger than the areas of negative triggering (Figures 9 and 10), we can say that the hypothesis of static stress triggering does not make a large improvement over random chance for explaining the locations of triggered earthquakes. [31] It seems unlikely that this sequence is a coincidence of large events, and we think there is some triggering mechanism(s) that relates the M w 8.2, M w 7.5, and M w 7.4 earthquakes, along with the smaller clusters that occurred in the wide surrounding area. If static stress triggering does not provide a good general explanation for all of the triggered events, what is the mechanism? We do not have any strong evidence for other types of mechanisms, but the possibility of dynamic triggering [e.g., Gomberg et al., 2003] from the seismic waves has been proposed as a mechanism that might be able to explain all the triggered events. Furthermore, the study of Pollitz and Johnston [2006] suggests that even a large percentage of ordinary aftershocks may be due to dynamic triggering. [32] It is intriguing that there was also a M6 event on 7 November, near the epicenter (within 25 km) of the subsequent M w November event (Figure 3). One might speculate that the region of the 17 November event was close to failure so that it could have been triggered by the M6 event. There have been some similar examples in regions of triggered seismicity where smaller foreshocks led to a larger triggered earthquake within hours to days. These include the M6.5 Big Bear event following the M7.4 Landers earthquake [King et al., 1994], the M5.5 event following the M western Tottori, Japan, earthquake [Ohmi et al., 2002], the M7.4 event following the M Kii Peninsula, Japan, earthquake [Park and Mori, 2005]. [33] The triggered normal faulting earthquakes in group 2 are tensional earthquakes near the outer rise. As described by Christensen and Ruff [1988], the tensional stresses perpendicular to the arc in this region could be increased by large slip downdip on the subduction interface. Another example of such triggering of tensional earthquakes occurred after the 1990 M w 7.0 event in Costa Rica [Protti et al., 1995]. [34] Lay and Kanamori [1980] point out that the occurrences of large earthquakes in this region are often characterized by doublets, so that triggering of large earthquakes has often observed in the past. The example involving the largest earthquakes is for the two M8 events on 14 July and 26 July 1971 (Figure 2) on adjacent segments of the New Britain Trench. Using slip distributions calculated by Park and Mori [2007], we estimated the static stress changes from the first event in the region of the hypocenter of the second event. We found that the 26 July event occurred in a positive area with static stress changes of about 0.03 to 0.06 MPa, with the location very close to the boundary between stress increase and decrease. So the evaluations of the effectiveness of static stress changes depend on the details of the fault and slip geometry. [35] Our evaluation of static stress changes as a mechanism for the triggering of earthquakes in the New Britain/ New Ireland region has mixed results. The complicated sequences of interrelated earthquakes cannot all be explained by static stress changes, although the occurrences of some large and small earthquakes are consistent with changes in the static stress. At this stage it is difficult to judge if static stress changes are part of the explanation for the triggered earthquakes, or if the dominant mechanism is completely different. 5. Conclusions [36] We used teleseismic P wave data to obtain slip distributions for the M w 8.2 New Ireland earthquake in 2000 and the two following large earthquakes (M w 7.5 on 16 November and M w 7.4 on 17 November). The slip distribution for the New Ireland event showed large amounts of strike-slip and dip-slip displacements along a rupture length of about 200 km. The slip motions changed gradually from strike slip in the northwest to dip slip in the southeast. The two large following events to the south along the New Britain Trench had smaller and simpler slip distributions. For both events, the slip mainly occurred in a region close to the hypocenter. [37] We calculated the Coulomb failure stress changes to evaluate the effects of the static stress triggering on the initiation of the large earthquakes and groups of smaller triggered events. We found that the occurrences of the two M 7.5 earthquakes were consistent with the static changes; however, also examining other smaller earthquakes, static stress triggering does not provide a general explanation for all of the triggered events in the sequence. It is likely that dynamic stress triggering or dynamic and static stress triggering play an important role in explaining the 13 of 14

14 interrelated occurrences of the large and small earthquakes of this sequence. [38] Acknowledgments. We thank Marino Protti, the Associate Editor, and an anonymous reviewer for their critical comments. We used digital waveform data from the Incorporated Research Institutions for Seismology data center. Some of figures in this paper were made using the General Mapping Tool (GMT) software [Wessel and Smith, 1998]. References Christensen, D. H., and L. J. Ruff (1988), Seismic coupling and outer rise earthquakes, J. Geophys. Res., 93, 13,421 13,444. Dieterich, J. (1994), A constitutive law for rate of earthquake production and its application to earthquake clustering, J. Geophys. Res., 99, Díez, M., P. C. La Femina, C. B. Connor, W. Strauch, and V. Tenorio (2005), Evidence for static stress changes triggering the 1999 eruption of Cerro Negro Volcano, Nicaragua and regional aftershock sequences, Geophys. Res. Lett., 32, L04309, doi: /2004gl Dziewonski, A. M., T.-A. Chou, and J. H. Woodhouse (1981), Determination of earthquake source parameters from waveform data for studies of global and regional seismicity, J. Geophys. Res., 86, Gavrilenko, P. (2005), Hydromechanical coupling in response to earthquakes: On the possible consequences for aftershocks, Geophys. J. Int., 161, Geist, E. L., and T. Parsons (2005), Triggering of tsunamigenic aftershocks from large strike-slip earthquakes: Analysis of the November 2000 New Ireland earthquake sequence, Geochem. Geophys. Geosyst., 6, Q10005, doi: /2005gc Gomberg, J. (2001), The failure of earthquake failure models, J. Geophys. Res., 106, 16,253 16,263. Gomberg, J., P. Bodin, and P. A. Reasenberg (2003), Observing earthquakes triggering in the near field by dynamic deformations, Bull. Seismol. Soc. Am., 93, Hartzell, S. H., and T. H. Heaton (1983), Inversion of strong ground motion and teleseismic waveform data for the fault rupture history of the 1979 Imperial Valley, California earthquake, Bull. Seismol. Soc. Am., 73, King, G. C. P., R. S. Stein, and J. Lin (1994), Static stress changes and the triggering of earthquakes, Bull. Seismol. Soc. Am., 84, Langston, C. A., and D. V. Helmberger (1975), A procedure for modeling shallow dislocation sources, Geophys. J. R. Astron. Soc., 42, Lay, T., and H. Kanamori (1980), Earthquake doublets in the Solomon Islands, Phys. Earth Planet. Inter., 21, Ohmi, S., K. Watanabe, T. Shibutani, N. Hirano, and S. Nakao (2002), The 2000 western Tottori earthquake Seismic activity revealed by the regional seismic networks, Earth Planets Space, 54, Okada, Y. (1992), Internal deformation due to shear and tensile faults in a half-space, Bull. Seismol. Soc. Am., 82, Park, S.-C., and J. Mori (2005), The 2004 sequence of triggered earthquakes off the Kii peninsula, Japan, Earth Planets Space, 57, Park, S.-C., and J. Mori (2007), Are asperity patterns persistent? Implication from large earthquakes in Papua New Guinea, J. Geophys. Res., 112, B03303, doi: /2006jb Parsons, T. (2005), A hypothesis for delayed dynamic earthquake triggering, Geophys. Res. Lett., 32, L04302, doi: /2004gl Pollitz, F. F., and M. J. S. Johnston (2006), Direct test of static stress versus dynamic stress triggering of aftershocks, Geophys. Res. Lett., 33, L15318, doi: /2006gl Protti, M., et al. (1995), The March 25, 1990 (M w = 7.0, M L = 6.8), earthquake at the entrance of the Nicoya Gulf, Costa Rica: Its prior activity, foreshocks, aftershocks, and triggered seismicity, J. Geophys. Res., 100, 20,345 20,358. Rham, D. J., and S. Das (2003), Triggering of thrust aftershocks on the New Britain Trench by the November 16, 2000 M w 8.0 New Ireland strike-slip earthquake, paper presented at the General Assembly, Int. Union of Geol. and Geophys., Sapporo, Japan. Steacy, S., J. Gomberg, and M. Cocco (2005), Introduction to special section: Stress transfer, earthquake triggering, and time-dependent seismic hazard, J. Geophys. Res., 110, B05S01, doi: /2005jb Stein, R. S. (1999), The role of stress transfer in earthquake occurrence, Nature, 402, Toda, S., R. S. Stein, P. A. Reasenberg, J. H. Dieterich, and A. Yoshida (1998), Stress transferred by the 1995 M w = 6.9 Kobe, Japan, shock: Effect on aftershocks and future earthquake probabilities, J. Geophys. Res., 103, 24,543 24,565. Tregoning, P., et al. (1998), Estimation of current plate motions in Papua New Guinea from Global Positioning System observations, J. Geophys. Res., 103, 12,181 12,203. Tregoning, P., M. Sambridge, H. McQueen, S. Toulmin, and T. Nicholson (2005), Tectonic interpretation of aftershock relocations in eastern Papua New Guinea using teleseismic data and the arrival pattern method, Geophys. J. Int., 160, Wessel, P., and W. H. F. Smith (1998), New, improved version of the Generic Mapping Tools released, Eos Trans. AGU, 79(47), 579. J. Mori and S.-C. Park, Disaster Prevention Research Institute, Kyoto University, Gokasho, Uji, Kyoto, , Japan. (suncheon@eqh.dpri. kyoto-u.ac.jp) 14 of 14

Centroid moment-tensor analysis of the 2011 Tohoku earthquake. and its larger foreshocks and aftershocks

Centroid moment-tensor analysis of the 2011 Tohoku earthquake. and its larger foreshocks and aftershocks Earth Planets Space, 99, 1 8, 2011 Centroid moment-tensor analysis of the 2011 Tohoku earthquake and its larger foreshocks and aftershocks Meredith Nettles, Göran Ekström, and Howard C. Koss Lamont-Doherty

More information

Centroid-moment-tensor analysis of the 2011 off the Pacific coast of Tohoku Earthquake and its larger foreshocks and aftershocks

Centroid-moment-tensor analysis of the 2011 off the Pacific coast of Tohoku Earthquake and its larger foreshocks and aftershocks LETTER Earth Planets Space, 63, 519 523, 2011 Centroid-moment-tensor analysis of the 2011 off the Pacific coast of Tohoku Earthquake and its larger foreshocks and aftershocks Meredith Nettles, Göran Ekström,

More information

COULOMB STRESS CHANGES DUE TO RECENT ACEH EARTHQUAKES

COULOMB STRESS CHANGES DUE TO RECENT ACEH EARTHQUAKES COULOMB STRESS CHANGES DUE TO RECENT ACEH EARTHQUAKES Madlazim Physics Department, Faculty Mathematics and Sciences of Surabaya State University (UNESA) Jl. Ketintang, Surabaya 60231, Indonesia. e-mail:

More information

SOURCE 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 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 information

Slip 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 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 information

RELOCATION 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 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 information

Analysis of Seismological and Tsunami Data from the 1993 Guam Earthquake

Analysis of Seismological and Tsunami Data from the 1993 Guam Earthquake PAGEOPH, Vol. 144, Nos. 3/4 (1995) 0033-4553/95/040823-1551.50 + 0.20/0 9 1995 Birkh/iuser Verlag, Basel Analysis of Seismological and Tsunami Data from the 1993 Guam Earthquake YUICHIRO TANIOKA, 1 KENJI

More information

Earthquake Doublet Sequences: Evidence of Static Triggering in the Strong Convergent Zones of Taiwan

Earthquake Doublet Sequences: Evidence of Static Triggering in the Strong Convergent Zones of Taiwan Terr. Atmos. Ocean. Sci., Vol. 19, No. 6, 589-594, December 2008 doi: 10.3319/TAO.2008.19.6.589(PT) Earthquake Doublet Sequences: Evidence of Static Triggering in the Strong Convergent Zones of Taiwan

More information

Source of the July 2006 West Java tsunami estimated from tide gauge records

Source of the July 2006 West Java tsunami estimated from tide gauge records GEOPHYSICAL RESEARCH LETTERS, VOL. 33, L24317, doi:10.1029/2006gl028049, 2006 Source of the July 2006 West Java tsunami estimated from tide gauge records Yushiro Fujii 1 and Kenji Satake 2 Received 13

More information

Source rupture process of the 2003 Tokachi-oki earthquake determined by joint inversion of teleseismic body wave and strong ground motion data

Source rupture process of the 2003 Tokachi-oki earthquake determined by joint inversion of teleseismic body wave and strong ground motion data LETTER Earth Planets Space, 56, 311 316, 2004 Source rupture process of the 2003 Tokachi-oki earthquake determined by joint inversion of teleseismic body wave and strong ground motion data Yuji Yagi International

More information

Coulomb stress change for the normal-fault aftershocks triggered near the Japan Trench by the 2011 M w 9.0 Tohoku-Oki earthquake

Coulomb stress change for the normal-fault aftershocks triggered near the Japan Trench by the 2011 M w 9.0 Tohoku-Oki earthquake Earth Planets Space, 64, 1239 1243, 2012 Coulomb stress change for the normal-fault aftershocks triggered near the Japan Trench by the 2011 M w 9.0 Tohoku-Oki earthquake Tamao Sato 1, Shinya Hiratsuka

More information

Widespread Ground Motion Distribution Caused by Rupture Directivity during the 2015 Gorkha, Nepal Earthquake

Widespread Ground Motion Distribution Caused by Rupture Directivity during the 2015 Gorkha, Nepal Earthquake Widespread Ground Motion Distribution Caused by Rupture Directivity during the 2015 Gorkha, Nepal Earthquake Kazuki Koketsu 1, Hiroe Miyake 2, Srinagesh Davuluri 3 and Soma Nath Sapkota 4 1. Corresponding

More information

Synthetic Seismicity Models of Multiple Interacting Faults

Synthetic Seismicity Models of Multiple Interacting Faults Synthetic Seismicity Models of Multiple Interacting Faults Russell Robinson and Rafael Benites Institute of Geological & Nuclear Sciences, Box 30368, Lower Hutt, New Zealand (email: r.robinson@gns.cri.nz).

More information

Empirical Green s Function Analysis of the Wells, Nevada, Earthquake Source

Empirical Green s Function Analysis of the Wells, Nevada, Earthquake Source Nevada Bureau of Mines and Geology Special Publication 36 Empirical Green s Function Analysis of the Wells, Nevada, Earthquake Source by Mendoza, C. 1 and Hartzell S. 2 1 Centro de Geociencias, Universidad

More information

Source Characteristics of Large Outer Rise Earthquakes in the Pacific Plate

Source Characteristics of Large Outer Rise Earthquakes in the Pacific Plate Source Characteristics of Large Outer Rise Earthquakes in the Pacific Plate T. Sasatani, N. Takai, M. Shigefuji, and Y. Miyahara Hokkaido University, Sapporo, Japan W. Kawabata Electric Power Development

More information

Seismic Activity near the Sunda and Andaman Trenches in the Sumatra Subduction Zone

Seismic Activity near the Sunda and Andaman Trenches in the Sumatra Subduction Zone IJMS 2017 vol. 4 (2): 49-54 International Journal of Multidisciplinary Studies (IJMS) Volume 4, Issue 2, 2017 DOI: http://doi.org/10.4038/ijms.v4i2.22 Seismic Activity near the Sunda and Andaman Trenches

More information

Coulomb stress changes due to Queensland earthquakes and the implications for seismic risk assessment

Coulomb stress changes due to Queensland earthquakes and the implications for seismic risk assessment Coulomb stress changes due to Queensland earthquakes and the implications for seismic risk assessment Abstract D. Weatherley University of Queensland Coulomb stress change analysis has been applied in

More information

Rupture Characteristics of Major and Great (M w 7.0) Megathrust Earthquakes from : 1. Source Parameter Scaling Relationships

Rupture Characteristics of Major and Great (M w 7.0) Megathrust Earthquakes from : 1. Source Parameter Scaling Relationships Journal of Geophysical Research Solid Earth Supporting Information for Rupture Characteristics of Major and Great (M w 7.0) Megathrust Earthquakes from 1990-2015: 1. Source Parameter Scaling Relationships

More information

Spatial and Temporal Distribution of Slip for the 1999 Chi-Chi, Taiwan, Earthquake

Spatial and Temporal Distribution of Slip for the 1999 Chi-Chi, Taiwan, Earthquake Bulletin of the Seismological Society of America, 91, 5, pp. 1069 1087, October 2001 Spatial and Temporal Distribution of Slip for the 1999 Chi-Chi, Taiwan, Earthquake by Kuo-Fong Ma, Jim Mori, Shiann-Jong

More information

Scaling relations of seismic moment, rupture area, average slip, and asperity size for M~9 subduction-zone earthquakes

Scaling 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 information

Coseismic slip distribution of the 1946 Nankai earthquake and aseismic slips caused by the earthquake

Coseismic 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 information

Earth Tides Can Trigger Shallow Thrust Fault Earthquakes

Earth Tides Can Trigger Shallow Thrust Fault Earthquakes Earth Tides Can Trigger Shallow Thrust Fault Earthquakes Elizabeth S. Cochran, 1 * John E. Vidale, 1 Sachiko Tanaka 2 1 Department of Earth and Space Sciences and Institute of Geophysics and Planetary

More information

Slip Partition of the 26 December 2006 Pingtung, Taiwan (M 6.9, M 6.8) Earthquake Doublet Determined from Teleseismic Waveforms

Slip Partition of the 26 December 2006 Pingtung, Taiwan (M 6.9, M 6.8) Earthquake Doublet Determined from Teleseismic Waveforms Terr. Atmos. Ocean. Sci., Vol. 19, No. 6, 567-578, December 2008 doi: 10.3319/TAO.2008.19.6.567(PT) Slip Partition of the 26 December 2006 Pingtung, Taiwan (M 6.9, M 6.8) Earthquake Doublet Determined

More information

EARTHQUAKE SOURCE PARAMETERS FOR SUBDUCTION ZONE EVENTS CAUSING TSUNAMIS IN AND AROUND THE PHILIPPINES

EARTHQUAKE SOURCE PARAMETERS FOR SUBDUCTION ZONE EVENTS CAUSING TSUNAMIS IN AND AROUND THE PHILIPPINES EARTHQUAKE SOURCE PARAMETERS FOR SUBDUCTION ZONE EVENTS CAUSING TSUNAMIS IN AND AROUND THE PHILIPPINES Joan Cruz SALCEDO Supervisor: Tatsuhiko HARA MEE09186 ABSTRACT We have made a set of earthquake source

More information

Multi-planar structures in the aftershock distribution of the Mid Niigata prefecture Earthquake in 2004

Multi-planar structures in the aftershock distribution of the Mid Niigata prefecture Earthquake in 2004 LETTER Earth Planets Space, 57, 411 416, 2005 Multi-planar structures in the aftershock distribution of the Mid Niigata prefecture Earthquake in 2004 Shigeki Aoki 1, Masaki Nishi 2, Koji Nakamura 2, Tetsuo

More information

revised October 30, 2001 Carlos Mendoza

revised October 30, 2001 Carlos Mendoza Earthquake Sources in the circum-caribbean Region Puerto Rico Tsunami Mitigation and Warning Program Federal Emergency Management Agency Preliminary Report: Task 3 revised October 30, 2001 Carlos Mendoza

More information

RELATION BETWEEN RAYLEIGH WAVES AND UPLIFT OF THE SEABED DUE TO SEISMIC FAULTING

RELATION BETWEEN RAYLEIGH WAVES AND UPLIFT OF THE SEABED DUE TO SEISMIC FAULTING 13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 24 Paper No. 1359 RELATION BETWEEN RAYLEIGH WAVES AND UPLIFT OF THE SEABED DUE TO SEISMIC FAULTING Shusaku INOUE 1,

More information

Moment tensor inversion of near source seismograms

Moment 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 information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY 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 information

Onto what planes should Coulomb stress perturbations be resolved?

Onto what planes should Coulomb stress perturbations be resolved? JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 110,, doi:10.1029/2004jb003356, 2005 Onto what planes should Coulomb stress perturbations be resolved? Sandy Steacy, Süleyman S. Nalbant, and John McCloskey School

More information

Magnitude 7.5 NEW BRITAIN REGION, PAPUA NEW GUINEA

Magnitude 7.5 NEW BRITAIN REGION, PAPUA NEW GUINEA A magnitude 7.5 earthquake struck off the eastern coast of Papua New Guinea on Sunday, approximately 54 km (33 miles) southeast of Kokopo. Residents reported strong ground shaking for about five minutes.

More information

GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L19604, doi: /2004gl020366, 2004

GEOPHYSICAL 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 information

RELOCATION OF LARGE EARTHQUAKES ALONG THE PHILIPPINE FAULT ZONE AND THEIR FAULT PLANES

RELOCATION OF LARGE EARTHQUAKES ALONG THE PHILIPPINE FAULT ZONE AND THEIR FAULT PLANES RELOCATION OF LARGE EARTHQUAKES ALONG THE PHILIPPINE FAULT ZONE AND THEIR FAULT PLANES Rey M. Lumbang MEE12608 Supervisor: Nobuo Hurukawa ABSTRACT We relocated large magnitude (Mw 7.0) earthquakes that

More information

ON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL

ON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL 1 Best Practices in Physics-based Fault Rupture Models for Seismic Hazard Assessment of Nuclear ON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL Hideo AOCHI

More information

Rapid source characterization of the 2011 M w 9.0 off the Pacific coast of Tohoku Earthquake

Rapid source characterization of the 2011 M w 9.0 off the Pacific coast of Tohoku Earthquake LETTER Earth Planets Space, 63, 529 534, 2011 Rapid source characterization of the 2011 M w 9.0 off the Pacific coast of Tohoku Earthquake Gavin P. Hayes 1,2 1 U.S. Geological Survey, National Earthquake

More information

Seismic Source Mechanism

Seismic Source Mechanism Seismic Source Mechanism Yuji Yagi (University of Tsukuba) Earthquake Earthquake is a term used to describe both failure process along a fault zone, and the resulting ground shaking and radiated seismic

More information

Effect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics

Effect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics Effect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics Naoyuki Kato (1) and Tomowo Hirasawa (2) (1) Geological

More information

Estimation of S-wave scattering coefficient in the mantle from envelope characteristics before and after the ScS arrival

Estimation 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 information

Scaling of apparent stress from broadband radiated energy catalogue and seismic moment catalogue and its focal mechanism dependence

Scaling of apparent stress from broadband radiated energy catalogue and seismic moment catalogue and its focal mechanism dependence Earth Planets Space, 53, 943 948, 2001 Scaling of apparent stress from broadband radiated energy catalogue and seismic moment catalogue and its focal mechanism dependence Z. L. Wu Institute of Geophysics,

More information

Chapter 2. Earthquake and Damage

Chapter 2. Earthquake and Damage EDM Report on the Chi-Chi, Taiwan Earthquake of September 21, 1999 2.1 Earthquake Fault 2.1.1 Tectonic Background The island of Taiwan is located in the complex junction where the Eurasian and Philippine

More information

The Size and Duration of the Sumatra-Andaman Earthquake from Far-Field Static Offsets

The Size and Duration of the Sumatra-Andaman Earthquake from Far-Field Static Offsets The Size and Duration of the Sumatra-Andaman Earthquake from Far-Field Static Offsets P. Banerjee, 1 F. F. Pollitz, 2 R. Bürgmann 3 * 1 Wadia Institute of Himalayan Geology, Dehra Dun, 248001, India. 2

More information

Short Note Source Mechanism and Rupture Directivity of the 18 May 2009 M W 4.6 Inglewood, California, Earthquake

Short Note Source Mechanism and Rupture Directivity of the 18 May 2009 M W 4.6 Inglewood, California, Earthquake Bulletin of the Seismological Society of America, Vol. 100, No. 6, pp. 3269 3277, December 2010, doi: 10.1785/0120100087 Short Note Source Mechanism and Rupture Directivity of the 18 May 2009 M W 4.6 Inglewood,

More information

Source processes of the 2009 Irian Jaya, Indonesia, earthquake doublet

Source processes of the 2009 Irian Jaya, Indonesia, earthquake doublet Earth Planets Space, 62, 475 481, 2010 Source processes of the 2009 Irian Jaya, Indonesia, earthquake doublet Natalia Poiata, Kazuki Koketsu, and Hiroe Miyake Earthquake Research Institute, University

More information

A GLOBAL MODEL FOR AFTERSHOCK BEHAVIOUR

A GLOBAL MODEL FOR AFTERSHOCK BEHAVIOUR A GLOBAL MODEL FOR AFTERSHOCK BEHAVIOUR Annemarie CHRISTOPHERSEN 1 And Euan G C SMITH 2 SUMMARY This paper considers the distribution of aftershocks in space, abundance, magnitude and time. Investigations

More information

LETTER Earth Planets Space, 56, , 2004

LETTER Earth Planets Space, 56, , 2004 LETTER Earth Planets Space, 56, 353 357, 2004 Deep seismic activities preceding the three large shallow earthquakes off south-east Hokkaido, Japan the 2003 Tokachi-oki earthquake, the 1993 Kushiro-oki

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION SUPPLEMENTARY INFORMATION DOI: 1.138/NGEO177 The Long Precursory Phase of Most Large Interplate Earthquakes Supplementary Information Supplementary Methods 1. Data and Classification We select the earthquakes

More information

Source characteristics of large deep earthquakes: Constraint on the faulting mechanism at great depths

Source characteristics of large deep earthquakes: Constraint on the faulting mechanism at great depths JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 108, NO. B2, 2091, doi:10.1029/2002jb001948, 2003 Source characteristics of large deep earthquakes: Constraint on the faulting mechanism at great depths Rigobert Tibi

More information

TSUNAMI CHARACTERISTICS OF OUTER-RISE EARTHQUAKES ALONG THE PACIFIC COAST OF NICARAGUA - A CASE STUDY FOR THE 2016 NICARAGUA EVENT-

TSUNAMI CHARACTERISTICS OF OUTER-RISE EARTHQUAKES ALONG THE PACIFIC COAST OF NICARAGUA - A CASE STUDY FOR THE 2016 NICARAGUA EVENT- TSUNAMI CHARACTERISTICS OF OUTER-RISE EARTHQUAKES ALONG THE PACIFIC COAST OF NICARAGUA - A CASE STUDY FOR THE 2016 NICARAGUA EVENT- Amilcar Cabrera Supervisor: Yuichiro TANIOKA MEE16718 ABSTRACT Nicaragua

More information

Tsunami waveform analyses of the 2006 underthrust and 2007 outer-rise Kurile earthquakes

Tsunami waveform analyses of the 2006 underthrust and 2007 outer-rise Kurile earthquakes Author(s) 2008. This work is licensed under a Creative Commons License. Advances in Geosciences Tsunami waveform analyses of the 2006 underthrust and 2007 outer-rise Kurile earthquakes Y. Tanioka 1, Y.

More information

overlie the seismogenic zone offshore Costa Rica, making the margin particularly well suited for combined land and ocean geophysical studies (Figure

overlie 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 information

Tohoku-oki event: Tectonic setting

Tohoku-oki event: Tectonic setting Tohoku-oki event: Tectonic setting This earthquake was the result of thrust faulting along or near the convergent plate boundary where the Pacific Plate subducts beneath Japan. This map also shows the

More information

Magnitude 7.1 NEAR THE EAST COAST OF HONSHU, JAPAN

Magnitude 7.1 NEAR THE EAST COAST OF HONSHU, JAPAN Japan was rattled by a strong aftershock and tsunami warning Thursday night nearly a month after a devastating earthquake and tsunami flattened the northeastern coast. This earthquake can be considered

More information

Aftershocks are well aligned with the background stress field, contradicting the hypothesis of highly heterogeneous crustal stress

Aftershocks are well aligned with the background stress field, contradicting the hypothesis of highly heterogeneous crustal stress JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2010jb007586, 2010 Aftershocks are well aligned with the background stress field, contradicting the hypothesis of highly heterogeneous crustal stress

More information

A search for seismic radiation from late slip for the December 26, 2004 Sumatra-Andaman (M w = 9.15) earthquake

A search for seismic radiation from late slip for the December 26, 2004 Sumatra-Andaman (M w = 9.15) earthquake Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 33, L18305, doi:10.1029/2006gl027286, 2006 A search for seismic radiation from late slip for the December 26, 2004 Sumatra-Andaman (M w =

More information

Magnitude 7.9 SE of KODIAK, ALASKA

Magnitude 7.9 SE of KODIAK, ALASKA A magnitude 7.9 earthquake occurred at 12:31 am local time 181 miles southeast of Kodiak at a depth of 25 km (15.5 miles). There are no immediate reports of damage or fatalities. Light shaking from this

More information

An intermediate deep earthquake rupturing on a dip-bending fault: Waveform analysis of the 2003 Miyagi-ken Oki earthquake

An 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 information

Magnitude 7.6 & 7.4 SOLOMON ISLANDS

Magnitude 7.6 & 7.4 SOLOMON ISLANDS A magnitude 7.6 earthquake struck near the Solomon Islands on Sunday morning local time; there were no immediate reports of damage. The earthquake was centered 100 km (60 miles) south of Kira Kira, a town

More information

Hitoshi Hirose (1), and Kazuro Hirahara (2) Abstract. Introduction

Hitoshi 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 information

Rupture complexity of the M w 8.3 sea of okhotsk earthquake: Rapid triggering of complementary earthquakes?

Rupture complexity of the M w 8.3 sea of okhotsk earthquake: Rapid triggering of complementary earthquakes? GEOPHYSICAL RESEARCH LETTERS, VOL. 40, 1 6, doi:10.1002/grl.50977, 2013 Rupture complexity of the M w 8.3 sea of okhotsk earthquake: Rapid triggering of complementary earthquakes? Shengji Wei, 1 Don Helmberger,

More information

X-2 HIKIMA AND KOKETSU: THE 2004 CHUETSU, JAPAN, EARTHQUAKE We relocated the hypocenters of the 2004 Chuetsu earthquake sequence, Niigata, Japan, usin

X-2 HIKIMA AND KOKETSU: THE 2004 CHUETSU, JAPAN, EARTHQUAKE We relocated the hypocenters of the 2004 Chuetsu earthquake sequence, Niigata, Japan, usin GEOPHYSICAL RESEARCH LETTERS, VOL. 32, XXXX, DOI:1029/2005GL023588, Rupture processes of the 2004 Chuetsu (mid-niigata prefecture) earthquake, Japan: A series of events in a complex fault system Kazuhito

More information

RELOCATION OF LARGE EARTHQUAKES ALONG THE SUMATRAN FAULT AND THEIR FAULT PLANES

RELOCATION OF LARGE EARTHQUAKES ALONG THE SUMATRAN FAULT AND THEIR FAULT PLANES Synopses of Master Papers Bulletin of IISEE, 47, 25-30, 2013 RELOCATION OF LARGE EARTHQUAKES ALONG THE SUMATRAN FAULT AND THEIR FAULT PLANES Biana Rahayu Wulandari MEE11605 Supervisor: Nobuo HURUKAWA ABSTRACT

More information

Teleseismic waveform modelling of the 2008 Leonidio event

Teleseismic waveform modelling of the 2008 Leonidio event The 6 January 2008 (Mw6.2) Leonidio (southern Greece) intermediate depth earthquake: teleseismic body wave modelling Anastasia Kiratzi and Christoforos Benetatos Department of Geophysics, Aristotle University

More information

Case Study 2: 2014 Iquique Sequence

Case Study 2: 2014 Iquique Sequence Case Study 2: 2014 Iquique Sequence Overview Mw 8.2 earthquake on 1 April 2014 at 11:08:43 UTC Pictures of damage Seismicity maps Foreshock seismicity -> Main shock -> Aftershock 1 April 2014 11:08:43

More information

Materials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation.

Materials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation. Materials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation. In the process zone, stress amplitudes are poorly determined and much

More information

SOURCE PROCESS OF THE 2003 PUERTO PLATA EARTHQUAKE USING TELESEISMIC DATA AND STRONG GROUND MOTION SIMULATION

SOURCE PROCESS OF THE 2003 PUERTO PLATA EARTHQUAKE USING TELESEISMIC DATA AND STRONG GROUND MOTION SIMULATION Synopses of Master Papers Bulletin of IISEE, 47, 19-24, 2013 SOURCE PROCESS OF THE 2003 PUERTO PLATA EARTHQUAKE USING TELESEISMIC DATA AND STRONG GROUND MOTION SIMULATION Fabricio Moquete Everth* Supervisor:

More information

Rapid magnitude determination from peak amplitudes at local stations

Rapid magnitude determination from peak amplitudes at local stations Earth Planets Space, 65, 843 853, 2013 Rapid magnitude determination from peak amplitudes at local stations Akio Katsumata 1, Hiroshi Ueno 1, Shigeki Aoki 1, Yasuhiro Yoshida 2, and Sergio Barrientos 3

More information

Supplementary Materials for

Supplementary Materials for advances.sciencemag.org/cgi/content/full/4/3/eaao4915/dc1 Supplementary Materials for Global variations of large megathrust earthquake rupture characteristics This PDF file includes: Lingling Ye, Hiroo

More information

Determination of fault planes in a complex aftershock sequence using two-dimensional slip inversion

Determination of fault planes in a complex aftershock sequence using two-dimensional slip inversion Geophys. J. Int. (2001) 146, 134 142 Determination of fault planes in a complex aftershock sequence using two-dimensional slip inversion Rachel E. Abercrombie, 1,2, * Stephen Bannister, 1 Aasha Pancha,

More information

FOCAL MECHANISM DETERMINATION OF LOCAL EARTHQUAKES IN MALAY PENINSULA

FOCAL MECHANISM DETERMINATION OF LOCAL EARTHQUAKES IN MALAY PENINSULA FOCAL MECHANISM DETERMINATION OF LOCAL EARTHQUAKES IN MALAY PENINSULA Siti Norbaizura MAT SAID Supervisor: Tatsuhiko HARA MEE10505 ABSTRACT Since November 30, 2007, small local earthquakes have been observed

More information

1.3 Short Review: Preliminary results and observations of the December 2004 Great Sumatra Earthquake Kenji Hirata

1.3 Short Review: Preliminary results and observations of the December 2004 Great Sumatra Earthquake Kenji Hirata 1.3 Short Review: Preliminary results and observations of the December 2004 Great Sumatra Earthquake Kenji Hirata We give a brief review about observations and preliminary results regarding the 2004 great

More information

Effect of the Emperor seamounts on trans-oceanic propagation of the 2006 Kuril Island earthquake tsunami

Effect 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 information

Comment on A new estimate for present-day Cocos-Caribbean plate motion: Implications

Comment on A new estimate for present-day Cocos-Caribbean plate motion: Implications Comment on A new estimate for present-day Cocos-Caribbean plate motion: Implications for slip along the Central American volcanic arc by Charles DeMets Marco Guzmán-Speziale Juan Martín Gómez Unidad de

More information

Spatio-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 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 information

AVERAGE AND VARIATION OF FOCAL MECHANISM AROUND TOHOKU SUBDUCTION ZONE

AVERAGE 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 information

Crustal deformation by the Southeast-off Kii Peninsula Earthquake

Crustal deformation by the Southeast-off Kii Peninsula Earthquake Crustal deformation by the Southeast-off Kii Peninsula Earthquake 51 Crustal deformation by the Southeast-off Kii Peninsula Earthquake Tetsuro IMAKIIRE, Shinzaburo OZAWA, Hiroshi YARAI, Takuya NISHIMURA

More information

DETERMINATION OF SLIP DISTRIBUTION OF THE 28 MARCH 2005 NIAS EARTHQUAKE USING JOINT INVERSION OF TSUNAMI WAVEFORM AND GPS DATA

DETERMINATION OF SLIP DISTRIBUTION OF THE 28 MARCH 2005 NIAS EARTHQUAKE USING JOINT INVERSION OF TSUNAMI WAVEFORM AND GPS DATA Synopses of Master Papers Bulletin of IISEE, 47, 115-10, 013 DETERMINATION OF SLIP DISTRIBUTION OF THE 8 MARCH 005 NIAS EARTHQUAKE USING JOINT INVERSION OF TSUNAMI WAVEFORM AND GPS DATA Tatok Yatimantoro

More information

LETTER Earth Planets Space, 57, , 2005

LETTER Earth Planets Space, 57, , 2005 LETTER Earth Planets Space, 57, 345 35, 25 Estimation of the source model for the foreshock of the 24 off the Kii peninsula earthquakes and strong ground motion simulation of the hypothetical Tonankai

More information

FOCAL MECHANISMS OF SUBDUCTION ZONE EARTHQUAKES ALONG THE JAVA TRENCH: PRELIMINARY STUDY FOR THE PSHA FOR YOGYAKARTA REGION, INDONESIA

FOCAL MECHANISMS OF SUBDUCTION ZONE EARTHQUAKES ALONG THE JAVA TRENCH: PRELIMINARY STUDY FOR THE PSHA FOR YOGYAKARTA REGION, INDONESIA FOCAL MECHANISMS OF SUBDUCTION ZONE EARTHQUAKES ALONG THE JAVA TRENCH: PRELIMINARY STUDY FOR THE PSHA FOR YOGYAKARTA REGION, INDONESIA Myo Thant 1, Hiroshi Kawase 2, Subagyo Pramumijoyo 3, Heru Hendrayana

More information

SOURCE INVERSION AND INUNDATION MODELING TECHNOLOGIES FOR TSUNAMI HAZARD ASSESSMENT, CASE STUDY: 2001 PERU TSUNAMI

SOURCE INVERSION AND INUNDATION MODELING TECHNOLOGIES FOR TSUNAMI HAZARD ASSESSMENT, CASE STUDY: 2001 PERU TSUNAMI Paper No. TS-4-1 SOURCE INVERSION AND INUNDATION MODELING TECHNOLOGIES FOR TSUNAMI HAZARD ASSESSMENT, CASE STUDY: 2001 PERU TSUNAMI Bruno Adriano 1, Shunichi Koshimura 2 and Yushiro Fujii 3 ABSTRACT The

More information

Three Dimensional Simulations of Tsunami Generation and Propagation

Three Dimensional Simulations of Tsunami Generation and Propagation Chapter 1 Earth Science Three Dimensional Simulations of Tsunami Generation and Propagation Project Representative Takashi Furumura Authors Tatsuhiko Saito Takashi Furumura Earthquake Research Institute,

More information

Magnitude 7.5 NEW BRITAIN REGION, PAPUA NEW GUINEA

Magnitude 7.5 NEW BRITAIN REGION, PAPUA NEW GUINEA A magnitude 7.5 earthquake struck off the eastern coast of Papua New Guinea on Tuesday, approximately 130 km (81 mi) southsouthwest of Kokopo at a depth of 42 km (26.1 mi). There are reports of some structural

More information

Crustal deformation in Taiwan: Results from finite source inversions of six M w > 5.8 Chi-Chi aftershocks

Crustal deformation in Taiwan: Results from finite source inversions of six M w > 5.8 Chi-Chi aftershocks JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 109,, doi:10.1029/2003jb002606, 2004 Crustal deformation in Taiwan: Results from finite source inversions of six M w > 5.8 Chi-Chi aftershocks Wu-Cheng Chi 1 and Doug

More information

Lecture 20: Slow Slip Events and Stress Transfer. GEOS 655 Tectonic Geodesy Jeff Freymueller

Lecture 20: Slow Slip Events and Stress Transfer. GEOS 655 Tectonic Geodesy Jeff Freymueller Lecture 20: Slow Slip Events and Stress Transfer GEOS 655 Tectonic Geodesy Jeff Freymueller Slow Slip Events From Kristine Larson What is a Slow Slip Event? Slip on a fault, like in an earthquake, BUT

More information

RECIPE FOR PREDICTING STRONG GROUND MOTIONS FROM FUTURE LARGE INTRASLAB EARTHQUAKES

RECIPE FOR PREDICTING STRONG GROUND MOTIONS FROM FUTURE LARGE INTRASLAB EARTHQUAKES RECIPE FOR PREDICTING STRONG GROUND MOTIONS FROM FUTURE LARGE INTRASLAB EARTHQUAKES T. Sasatani 1, S. Noguchi, T. Maeda 3, and N. Morikawa 4 1 Professor, Graduate School of Engineering, Hokkaido University,

More information

Static stress transfer from the May 20, 2012, M 6.1 Emilia-Romagna (northern Italy) earthquake using a co-seismic slip distribution model

Static stress transfer from the May 20, 2012, M 6.1 Emilia-Romagna (northern Italy) earthquake using a co-seismic slip distribution model 2012 EMILIA EARTHQUAKES ANNALS OF GEOPHYSICS, 55, 4, 2012; doi: 10.4401/ag-6176 Static stress transfer from the May 20, 2012, M 6.1 Emilia-Romagna (northern Italy) earthquake using a co-seismic slip distribution

More information

FOCAL MECHANISM DETERMINATION USING WAVEFORM DATA FROM A BROADBAND STATION IN THE PHILIPPINES

FOCAL MECHANISM DETERMINATION USING WAVEFORM DATA FROM A BROADBAND STATION IN THE PHILIPPINES FOCAL MECHANISM DETERMINATION USING WAVEFORM DATA FROM A BROADBAND STATION IN THE PHILIPPINES Vilma Castillejos Hernandez Supervisor: Tatsuhiko Hara MEE10508 ABSTRACT We performed time domain moment tensor

More information

JCR (2 ), JGR- (1 ) (4 ) 11, EPSL GRL BSSA

JCR (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 information

27th Seismic Research Review: Ground-Based Nuclear Explosion Monitoring Technologies

27th Seismic Research Review: Ground-Based Nuclear Explosion Monitoring Technologies GROUND TRUTH LOCATIONS USING SYNERGY BETWEEN REMOTE SENSING AND SEISMIC METHODS-APPLICATION TO CHINESE AND NORTH AFRICAN EARTHQUAKES C. K. Saikia 1, H. K. Thio 2, D. V. Helmberger 2, G. Ichinose 1, and

More information

Earthquake Focal Mechanisms and Waveform Modeling

Earthquake Focal Mechanisms and Waveform Modeling Earthquake Focal Mechanisms and Waveform Modeling Rengin Gök Lawrence Livermore National Laboratory USA RELEMR Workshop İstanbul 2008 Gudmundar E. Sigvaldason The Dynamic Earth, USGS The size of the event

More information

Lesvos June 12, 2017, Mw 6.3 event, a quick study of the source

Lesvos June 12, 2017, Mw 6.3 event, a quick study of the source Lesvos June 12, 2017, Mw 6.3 event, a quick study of the source E. Sokos 1) and J. Zahradník 2) 1) University of Patras, Greece 2) Charles University, Prague, Czech Republic Report sent to EMSC on 21/06/2017

More information

Partial Breaking of a Mature Seismic Gap: The 1987 Earthquakes in New Britain

Partial Breaking of a Mature Seismic Gap: The 1987 Earthquakes in New Britain PAGEOPH, Vol. 136, No. 4 (1991) 0033-4553/91/040459-1951.50 + 0.20/0 9 I991 Birkh/iuser Verlag, Basel Partial Breaking of a Mature Seismic Gap: The 1987 Earthquakes in New Britain R. DMOWSKA, 1'2 L. C.

More information

Slab pull, slab weakening, and their relation to deep intra-slab seismicity

Slab pull, slab weakening, and their relation to deep intra-slab seismicity GEOPHYSICAL RESEARCH LETTERS, VOL. 32, L14305, doi:10.1029/2005gl022922, 2005 Slab pull, slab weakening, and their relation to deep intra-slab seismicity Susan L. Bilek Earth and Environmental Science

More information

Routine Estimation of Earthquake Source Complexity: the 18 October 1992 Colombian Earthquake

Routine Estimation of Earthquake Source Complexity: the 18 October 1992 Colombian Earthquake Bulletin of the Seismological Society of America, Vol. 84, No. 4, pp. 1266-1271, August 1994 Routine Estimation of Earthquake Source Complexity: the 18 October 1992 Colombian Earthquake by Charles J. Ammon,*

More information

The Mw 6.2 Leonidio, southern Greece earthquake of January 6, 2008: Preliminary identification of the fault plane.

The Mw 6.2 Leonidio, southern Greece earthquake of January 6, 2008: Preliminary identification of the fault plane. The Mw 6.2 Leonidio, southern Greece earthquake of January 6, 28: Preliminary identification of the fault plane. J. Zahradnik 1, E. Sokos 2, A.Serpetsidaki 2, and G A.Tselentis 2 1 Charles University in

More information

Tsunami Simulation of 2009 Dusky Sound Earthquake in New Zealand

Tsunami Simulation of 2009 Dusky Sound Earthquake in New Zealand Tsunami Simulation of 2009 Dusky Sound Earthquake in New Zealand Polina Berezina 1 Institute of Geology, Taras Shevchenko National University of Kyiv, Kyiv, Ukraine Supervisor: Prof. Kenji Satake Earthquake

More information

Rotation of the Principal Stress Directions Due to Earthquake Faulting and Its Seismological Implications

Rotation of the Principal Stress Directions Due to Earthquake Faulting and Its Seismological Implications Bulletin of the Seismological Society of America, Vol. 85, No. 5, pp. 1513-1517, October 1995 Rotation of the Principal Stress Directions Due to Earthquake Faulting and Its Seismological Implications by

More information

Earthquake Source. Kazuki Koketsu. Special Session: Great East Japan (Tohoku) Earthquake. Earthquake Research Institute, University of Tokyo

Earthquake Source. Kazuki Koketsu. Special Session: Great East Japan (Tohoku) Earthquake. Earthquake Research Institute, University of Tokyo 2012/9/24 17:20-17:35 WCEE SS24.4 Special Session: Great East Japan (Tohoku) Earthquake Earthquake Source Kazuki Koketsu Earthquake Research Institute, University of Tokyo 1 Names and features of the earthquake

More information

Preliminary slip model of M9 Tohoku earthquake from strongmotion stations in Japan - an extreme application of ISOLA code.

Preliminary slip model of M9 Tohoku earthquake from strongmotion stations in Japan - an extreme application of ISOLA code. Preliminary slip model of M9 Tohoku earthquake from strongmotion stations in Japan - an extreme application of ISOLA code. J. Zahradnik 1), F. Gallovic 1), E. Sokos 2) G-A. Tselentis 2) 1) Charles University

More information

Rupture Process of the Great 2004 Sumatra-Andaman Earthquake

Rupture Process of the Great 2004 Sumatra-Andaman Earthquake Rupture Process of the Great 2004 Sumatra-Andaman Earthquake Supporting Online Materials Submitted to Science, March 12, 2005 Charles J. Ammon 1, Ji Chen 2, Hong-Kie Thio 3, David Robinson 5, Sidao Ni

More information

Rupture process of the 2005 West Off Fukuoka Prefecture, Japan, earthquake

Rupture 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 information