Questions? Call 800-624-6242

| Items in cart [0]

## (NAS Colloquium) Earthquake Prediction: The Scientific Challenge (1996) National Academy of Sciences (NAS)

### Citation Manager

. "Nonuniformity of the constitutive law parameters for shear rupture and quasistatic nucleation to dynamic rupture: A physical model of earthquake generation processes." (NAS Colloquium) Earthquake Prediction: The Scientific Challenge. Washington, DC: The National Academies Press, 1996.

Please select a format:

 Page 3802

The following HTML text is provided to enhance online readability. Many aspects of typography translate only awkwardly to HTML. Please use the page image as the authoritative form to ensure accuracy.

Proceedings of the National Academy of Sciences of the United States of America

Laboratory data (Fig. 3) show, however, that both ∆τb and Dc are an increasing function of σn, so that the assumption of ∆τb being constant is not justified for more rigorous discussion. From Eqs. 4, 5, and 20,

[21]

which shows that Xc is directly proportional to λc, and that Xc is independent of σn if M=1. If M≠1, Xc depends not only on λc, but also on σn. However, the dependence of {m(σn)}1−M on σn is modest, so that Xc is practically prescribed by λc alone. One can thus conclude from Eqs. 19 and 21 that Xc and Lc are virtually prescribed and scaled by λc alone. If it is tentatively assumed that the critical size of the nucleation zone is attained at a value in the V/Vs range 0.1–0.5, and that τi/∆τb=0.5−0.8, then Lc=(4.9−9.7)×103 λc from Eqs. 19 and 21 by considering that μ=2×104 MPa for the granite sample used, and assuming that M=1 and ∆τb0=3 MPa (see Fig. 3). This theoretical estimate agrees with the experimental result (Fig. 4).

During the nucleation, the energy is exclusively consumed in and around the zone of nucleation, since the nucleation is the process (or zone) in which slip failure deformation is concentrated and accelerated. Accordingly, dynamic instabilities of small to microscopic scales (microseismicity or acoustic emission) are induced and activated during the nucleation process, when the rupture growth resistance on the fault varies on small to microscopic scales. This has been demonstrated by recent laboratory experiments (unpublished data). It has also been found that seismic b values decrease as the nucleation proceeds (unpublished data). These observations suggest that microcrack interactions in the fault zone play a crucial role in the nucleation process.

### Earthquake Rupture Nucleation and Immediate Precursors

We have seen that geometrical and/or mechanical inhomogeneities play a crucial role in the shear rupture nucleation process. Such inhomogeneities prevail in the brittle seismogenic layer in the Earth’s crust. The seismogenic layer contains a large number of preexisting faults of microscopic to macroscopic scales, and therefore the seismogenic layer is inherently inhomogeneous. In addition, a preexisting fault itself in the Earth’s crust exhibits geometrical irregularities and mechanical inhomogeneities of various scales on the fault surfaces and in the fault zones. This strongly suggests that earthquake dynamic rupture is necessarily preceded by a quasistatic to quasidynamic nucleation process. However, whether or not a sizable zone of the nucleation appears prior to earthquake dynamic instability depends on how nonuniformly the constitutive law parameters prevail on the fault in the lithosphere (18).

The shear rupture nucleation model presented in this paper (see also ref. 18) can explain why and how short-term (or immediate) precursors are intrinsically related to the earthquake nucleation that proceeds quasistatically to quasidynamically prior to the mainshock dynamic rupture and how essential it is in carrying immediate foreshocks that the rupture growth resistance is distributed nonuniformly on a local to small scale in the fault zone in the brittle regime. The model shows that immediate foreshock activity is a part of the mainshock earthquake nucleation (18, 22). This provides a physical explanation for observations commonly made for decades that immediate foreshocks are concentrated in the vicinity of the epicenter of the pending mainshock. Whether or not immediate foreshocks occur during the mainshock nucleation depends on how the rupture growth resistance varies on a local to small scale in the nucleation zone (18, 22).

During the nucleation, local shear stresses decrease gradually in the breakdown zone, and at the same time the corresponding premonitory slip also proceeds in the zone, since the shear strength degrades with ongoing slip in the nucleation zone. Premonitory stress (or strain) changes are also observed outside (but adjacent to) the nucleation zone; that is, local shear stresses on the remaining unslipped parts adjacent to the nucleation zone increase with time because the unslipped segments must bear extra stress loads that have been sustained by the slipped parts. These gradual, accelerating changes in both local stress and slip are inevitable precursors that occur locally in or adjacent to the zone of the rupture nucleation. However, no such precursory slip and stress degradation necessarily occurs in a region distant from the nucleation zone, and the precursory deformation and stress changes can locally be confined in (or adjacent to) the zone of the nucleation. This suggests that the key to the short-term (or immediate) earthquake prediction is to identify where the nucleation occurs on the fault that has the potential to cause a major earthquake.

If the critical size of the nucleation zone for a real major earthquake is small enough, it may be meaningless to regard the nucleation process as an effective tool for the immediate prediction. In this sense, it is important to know how large is the critical size of the nucleation zone for a real major earthquake. The critical size of the nucleation zone has recently been estimated for a number of major earthquakes (B.Shibazaki and M.Matsu’ura, personal communication; see also refs. 22 and 23), showing that Lc for major earthquakes with M=7.0−7.7 is 5−10 km.

For a given Lc of 5–10 km, λc=0.5−2 m from Eqs. 19 and 21. This estimate suggests that for earthquakes of M=7.0−7.7, the breakdown process behind the tip of propagating rupture is virtually governed by a characteristic length of the order of 1 m, which is 104–105 times greater than λc for shear rupture in the laboratory.

I am deeply grateful to Professor L.Knopoff and the other organizers for inviting me to present this paper and for their courtesy during the Colloquium.

1. Ohnaka, M., Kuwahara, Y. & Yamamoto, K. (1987) Tectonophysics 144, 109–125.

2. Ohnaka, M. & Kuwahara, Y. (1990) Tectonophysics 175, 197–220.

3. Ohnaka, M. & Yamashita, T. (1989) J. Geophys. Res. 94, 4089–4104.

4. Rabinowicz, E. (1965) Friction and Wear of Materials (Wiley, New York).

5. Bowden, F.B. & Tabor, D. (1954) The Friction and Lubrication of Solids (Clarendon, Oxford).

6. Byerlee, J.D. (1967) J. Appl Phys. 38, 2928–2934.

7. Ohnaka, M. (1975) J. Phys. Earth 23, 87–112.

8. Okubo, P.G. & Dieterich, J.H. (1986) Geophys. Monogr. Am. Geophys. Union 37, 25–35.

9. Aviles, C.A., Scholz, C.H. & Boatwright, J. (1987) J. Geophys. Res. 92, 331–344.

10. Okubo, P.G. & Aki, K. (1987) J. Geophys. Res. 92, 345–355.

11. Ohnaka, M. (1995) in Theory of Earthquake Premonitory and Fracture Processes, ed. Teisseyre, R. (Polish Sci. Publ., Warsaw, Poland).

12. Ohnaka, M. (1990) Can. J. Phys. 68, 1071–1083.

13. Brown, S.R. & Scholz, C.H. (1985) J. Geophys. Res. 90, 12575–12582.

14. Kuwahara, Y., Ohnaka, M., Yamamoto, K. & Hirasawa, T. (1985) Annual Meet. Seis. Soc. Japan Progr. Abstr. 2, 110.

15. Ida, Y. (1973) Bull. Seis. Soc. Am. 63, 959–968.

16. Biegel, R.L., Wang, W., Scholz, C.H. & Boitnott, G.N. (1992) J. Geophys. Res. 97, 8951–8964.

17. Matsu’ura, M., Kataoka, H. & Shibazaki, B. (1992) Tectonophysics 211, 135–148.

18. Ohnaka, M. (1992) Tectonophysics 211, 149–178.

19. Shibazaki, B. & Matsu’ura, M. (1995) J. Geophys. Res., SUBMITTED.

20. Papageorgiou, A.S. & Aki, K. (1983) Bull. Seis. Soc. Am. 73, 953–978.

21. Kuwahara, Y., Ohnaka, M., Yamamoto, K. & Hirasawa, T. (1986) Annual Meet. Seis. Soc. Japan Progr. Abstr. 2, 233.

22. Yamashita, T. & Ohnaka, M. (1991) J. Geophys. Res. 96, 8351–8367.

23. Ohnaka, M., Yoshida, S. & Shen, L.-F. (1994) Workshop on Earthquake Source Modeling and Fault Mechanics, IASPEI General Assembly, Wellington, New Zealand.

24. Ohnaka, M. (1993) J. Phys. Earth 41, 45–56.

25. Ellsworth, W.L. & Beroza, G.C. (1995) Science, in press.

 Page 3802
 Front Matter (R1-R2) Earthquake prediction: The scientific challenge (3719-3720) Earthquake prediction: The interaction of public policy and science (3721-3725) Initiation process of earthquakes and its implications for seismic hazard reduction strategy (3726-3731) Intermediate- and long-term earthquake prediction (3732-3739) Scale dependence in earthquake phenomena and its relevance to earthquake prediction (3740-3747) Intermediate-term earthquake prediction (3748-3755) A selective phenomenology of the seismicity of Southern California (3756-3763) The repetition of large-earthquake ruptures (3764-3771) Hypothesis testing and earthquake prediction (3772-3775) What electrical measurements can say about changes in fault systems (3776-3780) Geochemical challenge to earthquake prediction (3781-3786) Implications of fault constitutive properties for earthquake prediction (3787-3794) Nonuniformity of the constitutive law parameters for shear rupture and quasistatic nucleation to dynamic rupture: A physical model of earthquake generation processes (3795-3802) Rock friction and its implications for earthquake prediction examined via models of Parkfield earthquakes (3803-3810) Slip complexity in earthquake fault models (3811-3818) Dynamic friction and the origin of the complexity of earthquake sources (3819-3824) Slip complexity in dynamic models of earthquake faults (3825-3829) The organization of seismicity on fault networks (3830-3837) Geometric incompatibility in a fault system (3838-3842)