ASPERITIES RECOGNIZED FROM SEISMICITY PATTERN CHANGE PRECEDING INTERPLATE MEGATHRUST

Shozo Matsumura National Research Institute for Earth Science and Disaster Prevention

Abstract

Recently, rupture processes of large earthquakes can be analyzed in detail to discriminate largely slipping area, that is, asperity by using various kinds of observation data. However such a manner is valid only to the that has already happened. In the view of seismic disaster prevention, we want to know where the asperity of a future earthquake should be located before its breakage. Replying to this request, we treat a case of the Tokai assumed earthquake. Analyzing the seismicity rate change, we could discriminate a patch-like pattern covering the assumed seismogenic zone. We interpret this pattern to reflect inhomogeneity in strength such that a stronger part can stand as it is, while a weaker part should slip to release the corresponding load. Then, the stress concentration might focus on the former one, which can be distinguished as an activated area in seismicity rate. And so, we recognize it as an asperity. If this is the case, we can extract a potential asperity for a future earthquake.

Introduction

The Tokai district, located around the central part of the Japan islands is well known as one of the most dangerous places in Japan since 1976, when Ishibashi(1976) gave a warning to the public that a mega-thrust earthquake had been impending around the Suruga bay, that is the next Tokai earthquake. This earthquake is regarded to belong to the historical Tokai earthquake series, the latest one of which was the 1944 Tonankai earthquake. The problem is that the fault area of the latest one did not reach inside the Suruga Bay. Based on this evidence, Ishibashi concluded that the Tokai area around the Suruga Bay was left as a potential seismogenic area, having been accumulating stress since the 1854 Ansei Tokai earthquake, and as a result, had been grown enough to arouse M8 earthquake. According to Ishibashi’s warning, the government organization, the Central Disaster Management Council (CDMC) started to establish a national program to promote a countermeasure act, one of the Figure 1. Seismogenic zone assumed for main subjects in which is to make an the next Tokai earthquake assessment of the influence brought by the earthquake.

Earthquake model

Ishibashi proposed his original prospective view for the seismogenic zone of the next Tokai earthquake by investigating the historical evidences that appeared with the occurrence of the 1854 Ansei Tokai earthquake. This zone has a rectangle shape enclosing the west Suruga Bay and the coastal area of the prefecture (see Fig.1). This was adopted as an official model by CDMC in 1978. Based on this model, various kinds of countermeasure acts had been set up by the government, the local autonomous bodies, and the public offices. However, it has already passed more than 20 years without any occurrence of the focused one. During the while, various kinds of observation networks have been developed covering the Tokai area in order to monitor and catch some anomalous signals even if it might be very small. One of representatives is a high sensitive seismic network, which was constructed and has been operated by our institute, the National Research Institute for Earth Science and Disaster Prevention (NIED). Another is the GPS network constructed and managed by the Geographical Survey Institute (GSI). Both observation networks could provide remarkable results, which led us to a new point of view for the impending Tokai earthquake. This resulted in a new version of the potential seismogenic zone as shown in the figure, an egg-plant like zone enclosed by the broken line.

Figure 2. Asperity model The next step is to evaluate the intensity of the placed inside the strong ground motion generated due to the earthquake. In assumed seismogenic order to estimate it through a numerical procedure, it is zone necessary to know the information respecting the asperities, that is, locations and sizes of the asperities. The committee of CDMC has proposed a procedure to set it up, that is called “the Irikura’s recipe”. In case of the Tokai earthquake, the whole seismogenic zone is divided into three segments, for each of which a couple of asperities are placed both in the inland zone, and in the oceanic zone, respectively. The total area of asperities is adjusted to be 30% of the whole. The result is shown in Fig.2. This was utilized for evaluation of strong motions and tsunami intensities. Then, some problems will be brought on; How much practical is such a model ? Is it impossible to recognize the asperities before occurrence of the main earthquake ?

2003 Tokachi-Oki Earthquake

We have a good example to answer to the proposed problem. The 2003 Tokachi-Oki earthquake (M8.0) happened around the off-coastal region of south-eastern Hokkaido on Sep. 26 in 2003, as a succeeding one of the 1952 Tokachi-Oki earthquake (M8.2), both of which belong to the sequence of interplate megathrust earthquakes between the Pacific and the North-American plates. Figure 3 shows the epicenter and the major asperities of the earthquake, analyzed and shown by Yagi(2004), and Tohoku University(2004).

Although any significant crustal movement such as slow-slip or pre-slip was never detected, a clear change of the seismicity could be recognized preceding the main earthquake. The seismicity rate inside the fault zone of the main earthquake indicated a slight but clear quiescence progressing since the middle of 1998, about 5 years preceding the occurrence (Katsumata and Kasahara, 2004). Figure 4(a) shows hypocentral distribution of microearthquakes occurred before the main earthquake, and (b) Figure 3. Slip distribution of the 2003 magnitude distribution for those inside the Tokachi-Oki earthquake. Yagi’s result and forward slip analysis based on GPS measurement. The star indicates the epicenter.

fault area. It is recognized that those earthquakes of M2.5 and greater are almost uniformly detected. Figure 5 shows the cumulative earthquake frequency counted from 1990, where earthquakes of M2.5 and greater and depth less than 100km are sampled, and processed through declustering. It is found that a clear depression of the rate commenced in 1998, and is somewhat recovered toward the occurrence. We consider that this is an epoch when some change has started in the locking state, and has been extending toward the entire fault zone. In order to investigate the spatial pattern of the change, we depicted a contour map of rate change between both seismicities before and after the epoch. Figure 6 shows the result, where 120 activated zones (shaded), and 150% zones (black-painted) are extracted. Comparing this with Fig.3, we are impressed in that there is a good correspondence of both locations between the activated zones and the asperities. Then, is there any necessity in such correspondence ?

Figure 4. (a) Hypocenter distribution, and (b) magnitude distribution before the occurrence of 2003 Tokachi-Oki earthquake.

Figure 5. Cumulative frequency of earthquakes (>M2.5) inside the fault zone of the Tokachi-oki earthquake.

We interpret this as follows. It might be probable that there is intrinsic Figure 6. Activated zones derived from inhomegeneity in strength of locking state on rate change between both seismicities the fault zone. However, it is not distinctive in before and after the epoch in 1998. the early stage of stress accumulation. In approaching the critical stage, some stress release will progress partly in the portions where strength is relatively weak, due to the rise of stress level. Even in this stage, the remaining strong portions will retain this slip to prevent an overall slip. At that time when most of the stress is supported by particular portions of strong locking, stress concentration will progress, and is going to appear. Then, this can be reflected on the seismicity change, and the asperity corresponding to the strong locking will be recognized as activated zones caused by stress concentration. This means that the potential asperity will be discovered by analyzing the seismicity pattern change inside the fault zone preceding the occurrence of the main earthquake.

Tokai Region

Again, we will discuss on the locking state in the Tokai region. We have been performing microearthquake observation and continuing to monitor the seismicity change covering this area more than two decades. In recent years, some delicate but significant changes of the seismicity were detected as similar like the Tokachi-oki case. Figure 7 shows hypocentral distribution in the Tokai region, where several blocks are selected for monitoring the seismicity. The first two, (a) and (b) are the top and the bottom of the layered activities, respectively. A small Figure 7. Hypocenter distribution in the rectangular area (c), surrounding the northeast Tokai region. Blocks a, b, c, and d are coast of Lake Hamana is also sampled. referred in Fig.8. Adding these, we sampled a distorted rectangular area extending to the northwest of Lake Hamana (d). The last one was selected as a reference to monitor the motion of the plate convergence. For these four blocked areas, cumulative number of earthquakes were counted, and displayed in Fig.8, where counting started from June 1986. Here, earthquakes of magnitude M1.5 and greater were sampled, and declusterd in order to reject step-like changes due to earthquake swarms and aftershocks. In these four graphs, we can recognize quiescence in recent years commonly for blocks (a), (b), and (c), while cannot find it only for (d), which is belonging to the extended portion of the subducted slab. The focal mechanism of earthquakes in (d) shows predominantly the normal fault type sustained by the stress with approximate EW extension (Matsumura, 1997). According to Ukawa(1982), this stress is attributed to the lateral stretching force generated from geometrical deformation of the subducted slab. In this case, the constant seismicity rate of (d) suggests that the convergence rate of the plate must be constant during the entire period. Conversely, it is suggested that the quiescence found in (a), (b), and (c) is not due to the stagnation of the plate movement itself, but should be attributed to some changes in the locking state, because all of them are located inside the assumed locked zone, and their focal mechanisms are consistent to the seismogenic stress caused by the interplate locking. However, there are some time differences in commencement of each quiescence, for example, in (a) and (c), it started from the end of 1996 to the beginning of 1997, while in (b) it started with a slight delay from August 1999. On the other hand, almost simultaneously with these seismic quiescences, some other change was detected in surface ground movement by GPS survey since the latter half of 2000. As a result, this was interpreted to be caused by a slow-slip progressing on the plate interface just beneath Lake Hamana (Ozawa et al., 2002).

Compiling all the current anomalous phenomena on a single figure of Fig.9, we summarize the analyzed results with the following

interpretation. First, we performed Figure 8. Cumulative frequency of earthquakes the analysis of extracting the since 1986. activated zones. In the figure, those

marked with A, B, and C (black-painted) are activated zones extracted from (b), which is the seismicity within the lower layer, that is, within the slab. Those of a, b, and c (shaded) are from (a) within the crust. The anomalous period selected for analyses is three years from October 2000 till September 2003, which is exactly corresponding to the period of slow-slip (They are still going on now). The slow-slip analyzed by the Geographical Survey Institute are also drawn in the figure. The longest arrow indicates the slow-slip during the three-year period located near Lake Hamana, which has exceeded 15 cm. It is found that the slow-slip area is localized to the west end of the inferred locked zone Figure 9. Activated zones in the Tokai region. and has not yet intruded into the A, B, and C are within the slab, and a, b, and c are interior. However, the simultaneous within the crust. Arrows are slow-slip progressing change in seisimicity is in progress on the plate interface. The star indicates epicenter inside the interior of the locked zone of the 1996 central Shizuoka-Pref. earthquake. by far from the area of the slow-slip. This implies that the place of the current progressing anomalies cannot be limited in the area of the present slipping but would cover even the interior of the locked zone. In Fig.9, we notice that there are three combinations of activated zones, A-a. B-b, and C-c, which are located inside the interior of the locked zone. Each set looks to align with a slight separating distance. Figure 10 provides an interpretation of such situation. As shown in the figure, when the stress is held on a particular portion with strong strength (that is an asperity), stress concentration should appear both in the top and bottom layers putting the asperity between them (Here, only compressional stress is taken into consideration). Looking under it from the sky, we will find a combination of stress concentration in the top and bottom layers with a separation putting an asperity between them. In Fig. 9, the arrangements of Figure 10. Image explaining (A-a) and (C-c) look to well match the above diagram, alignement of activated zones showing the alignment along the converging direction derived from stress concentration of the , suggesting that there is an induced in the slab and in the crust. asperity between the combinations. The arrangement of (B-b) conflicts with this image, although we don’t know the practical stress situation around this.

We can present another evidence showing the stress concentration on the activated zones. Figure 11 shows comparisons of magnitude distribution in the activated zones between both periods before and after the epochs of the seismicity change. The left figure is for those of 2003 Tokachi-oki

Figure 11. Comparisons of magnitude distribution in the activated zones between before (open triangle) and after (solid triangle) the epoch of the seismicity change. earthquake, the middle for the upper layer activity in the Tokai region (a) in Fig.8, and the right for the lower one (b). The open triangles correspond to the former period before the epoch, and the solid triangles to the later period. We can recognize that in all of these cases, the distribution of the latter one always indicates lower b-value than the former period. Such a change of b-value decrease is considered generally to represent existence of stress concentration.

Concluding Remarks

We have not yet experienced any megathrust earthquakes along the Nankai trough including the Tokai earthquakes after establishment of observations based on a modern technology. This means that it is essentially difficult to verify those hypotheses derived from various kinds of observations. However, it is required to make use of the real data in order not only to observe the phenomena, but to apply the result with some interpretations to another case in order to construct possible predictions. In case of the Tokai, we could have grasped the phenomena progressing on the assumed seismogenic zone, and have constructed a hypothesis respecting the existence of asperities. Based on an analogy to the case of 2003 Tokachi-oki earthquake, we could get valuable prospect respecting the coming process in Tokai as depicted in Fig.9, where three sets of activated zone are considered to reflect the places of potential asperities. It is interesting to compare it with Fig.2. At least three asperities of Fig.2 look well to match those of Fig.9. Such a result brings a great expectation that we could get an important prospect in evaluation of future earthquakes.

Referrences

Central Disaster Management Council, 2001, Report of the Technical Committee respecting the Tokai Earthquake, No.11 (in Japanese). Ishibashi, K., 1976, Reexamination of a great earthquake expected to occur in the Tokai district, central Japan – the great Suruga bay earthquake, Proc. Fall Meet. Seismol. Soc. Jpn., pp.30-34 (in Japanese). Katsumata, K. and M. Kasahara, 2004, Seismic quiescence and activation prior to the 2003 Tokachi-Oki earthquake, Proc. 2004 Japan Earth and Planetary Science Joint Meeting, S045-004. Matsumura, S., 1997, Focal zone of a future Tokai earthquake inferred from the seismicity pattern around the plate interface, Tectonophysics, 273, 271-291. Ozawa, S., M. Murakami, M. Kaizu, T.Tada, T. Sagiya, Y. Hatanaka, H. Yarai and T. Nishimura, 2002, Detection and monitoring of ongoing aseismic slip in the Tokai region, central Japan, Science, 298, 1009-1012. Tohoku University, 2004, On the M8.0 Tokachi-oki earthquake on September 26, 2003, Report Coordinating Committee Earthq. Predict., 71, 108-110 (in Japanese). Ukawa, M., 1982, Lateral stretching of the Philippine Sea plate subducting along the Nankai-Suruga trough, Tectonics, 1, 543-571. Yagi, Y., 2004, Source rupture process of the 2003 Tokachi-oki earthquake determined by joint inversion of teleseismic body wave and strong ground motion data, Earth Planets Space, 56, 311-316.