Numerical study on multi-scaling earthquake rupture
|
|
- Tracey Pitts
- 5 years ago
- Views:
Transcription
1 GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L266, doi:1.129/23gl1878, 24 Numerical study on multi-scaling earthquake rupture Hideo Aochi Institut de Radioprotection et de Sûreté Nucléaire, France Satoshi Ide Department of Earth and Planetary Science, University of Tokyo, Japan Received 2 September 23; revised 21 November 23; accepted 17 December 23; published 24 January 24. [1] A new numerical scheme using a renormalization and a 3D boundary integral equation method is proposed to simulate a multi-scaling dynamic rupture of earthquakes: How a small earthquake grows up to a large one in spatially heterogeneous field of critical slip-weakening distance D c (fracture energy G c )? We examine the case where D c grows according to a hypocentral distance L (D c / L b ). When b =1, we succeed to show numerically that a rupture propagates at a constant rupture speed in uniform initial stress field. This result still keeps the scaling relation of G c and D c inferred for earthquake size, however no scale-dependent initial process is required. The break of the proportional relation changes rupture speed as well as slip velocity to keep the energy balance. The rupture is accelerated up to a speed even faster than the shear wave velocity (b < 1) or naturally arrested (b > 1). INDEX TERMS: 721 Seismology: Earthquake parameters; 729 Seismology: Earthquake dynamics and mechanics; 726 Seismology: Theory and modeling; 81 Structural Geology: Fractures and faults. Citation: Aochi, H., and S. Ide (24), Numerical study on multi-scaling earthquake rupture, Geophys. Res. Lett., 31, L266, doi:1.129/ 23GL Introduction [2] When we discuss scaling problems in earthquake rupture, it is common to use fracture energy G c or critical slip displacement D c of the slip-weakening friction law, which is mathematically expressed in the form; tðu Þ ¼ t b ð1 u=d c ÞHð1 u=d c Þþt r : ð1þ Here H() is a Heaviside function and t b is the breakdown strength drop defined by the difference between peak and residual strengths, (t p t r ). The fracture energy follows the relation G c = 1 2 t b D c. Indeed, Ide and Takeo [1997] estimated G c of an order of 1 6 J/m 2 for earthquakes of magnitude around 7, while it is thought of an order of 1 4 J/m 2 in laboratory experiments [e.g., Scholz, 22]. It is also proposed that D c may reach in an order of 1 m for an earthquake of a magnitude 8, whereas it is actually observed in mm for rock experiments [e.g., Ohnaka, 23]. These D c -scaling and G c -scaling are consistent, as it is known in seismology that average stress change is about a few MPa independent of the scale. Shibazaki and Matsu ura [1992] modeled a concept of scale-dependent D c, leading to scale-dependent nucleation theory. Copyright 24 by the American Geophysical Union /4/23GL1878 [3] However, such parameters as D c and G c do not always need to keep a value of the same order in a fault system, because, in general, earthquakes repeatedly occur in a single fault system with a different magnitude of more than 6. This implies the existence of heterogeneity in D c or G c of an order of 1. This paper is concerned about this heterogeneity, which we shall think as an intrinsic feature on the fault to control important aspects of earthquake rupture. Although rupture behaviors based on energy balance at a crack tip have been studied analytically for mainly 2D models [e.g. Freund, 199], we need a flexible numerical scheme to solve more realistic problems. However, it has been very difficult to numerically do it. When we aim to model a large earthquake, it is customary to take a relatively larger model region, bigger grid size and longer time step than those for a small event. This is allowed because the scale-dependent character of D c or G c is explicitly or implicitly taken into account so that the relatively large grid and time step have resolution enough to follow the change in friction. A single grid for a large earthquake (a few hundred meters for a M7 event) is usually much larger than the rupture size of a small event (a hundred meter even for a M3 event). This is why no one can answer to the principal physical question how a large earthquake is different from small ones, or whether a small one grows up eventually to a large one [Ellsworth and Beroza, 199]. It has been impossible to figure out what really happens within a single grid of a large earthquake especially in its initial stage. This is due to the fact in how the grid size is chosen to study small and large events. Alternatively, we shall first propose a new numerical procedure to handle efficiently multi-scaling rupture events, and then consider the effect of heterogeneous, still scale-dependent distribution of D c or G c. 2. Numerical Procedure [4] The boundary integral equation method (BIEM) for a single planar fault in unbounded, homogeneous elastic medium [Fukuyama and Madariaga, 1998] is very useful for our purpose. The stress field t(~x, t) is generally written in the form of the spatio-temporal convolution (BIE) of the Green s function and the slip velocity history _u(~x, t)onthe fault plane. Namely the stress t ijk on the fault element can be rewritten in the discrete form of the summation of the discrete kernel P (i l)(j m)(k n) and discrete slip velocity V lmn, which is supposed to be constant on each spatial grid (s s) during each time step (t). Here the suffix (i, j) and (l, m) represent the spatial grid number, and k and n are for the time step, respectively. L266 1of
2 Figure 1. Schematic illustration of renormalization procedure in BIEM. Slip velocity on a single renormalized grid requires the 16 spatial grids over 4 time steps. [] Since the discrete slip velocity is a box function both in space and in time, it allows us to renormalize its grid size and time step, by conserving the rate of seismic moment (slip velocity grid surface duration). Letting a new grid size s, a new time step t and the slip velocity renormalized on the new grid V *, we will see V i j k * s2 t ¼ V ijk s 2 t ; ð2þ where the summation is taken over the grids covering the same fault area and over the same time duration. This procedure is also graphically shown in Figure 1. When we scale up by a factor of 4, for example, then s =4s, t = 4t, and V * i j k = (1/64)V ijk. [6] On the other hand, it is physically required that the total fracture energy on the fault is also conserved through this renormalization process. We can write it down in the same form as above; G i j c s 2 ¼ G ij c s2 ; ð3þ where G c * i j is the renormalized fracture energy on a grid (i, j ) based on those G ij c in the smaller scale. In the case where we renormalize by a factor of 4, G c * i j is thus composed by 16 grids in the smaller scale. As we assume that t b is constant for all scales, G c in the above equation can be replaced by D c. [7] By repeating this procedure, we are now able to follow the transit process from a very microscopic event to a macroscopic earthquake rupture. 3. Self-Similar Rupture Propagation [8] We shall see what happens in different scales during a single rupture process. In spite of the concept of scaledependent D c or G c explained above, some people consider that rupture always starts at a very small point corresponding to a very small D c or G c as observed in rock experiments. To avoid the contradiction between two opinions, D c or G c should grow up following the progress of rupture propagation itself. In fact this was implied by Richards [1976] in his review of analytic solution of 3D self similar rupture propagation. Such growth of G c,orr (resistance)-curve in the fracture mechanics, is widely observed in tensional fracture experiments [e.g., Lawn, 1993]. Here we generally suppose a simple relation with constant b; D c / L b ; where L is the hypocentral distance. There are a few explications for this general scaling law (Andrews, personal communication, 23). One is that D c is intrinsically distributed in the fractal manner on the fault. Rupture may starts almost everywhere at the smallest patch with the smallest G c (D c ), so that its averaged feature from the beginning point increases with distance to some extent, such a way as the above equation. Another is that D c is not a fixed property on fault, but arises as a result of dynamic process itself during the rupture. Size dependence of G c in material fracture arises because energy is lost to plastic yielding in a zone that becomes larger as the fracture length becomes larger. Andrews [1976] reported such the similar effect in dynamic shear rupture. [9] Table 1 shows model parameters for different scale in the case of b = 1, where D c is proportional to L. The other elastic parameters are common through the simulations; rigidity of 32.4 GPa and P/S wave velocity of 6./3.46 km/s. Strength t p, residual strength t r and initial stress t are MPa, MPa and 3 MPa, respectively, for all stage of simulations (scale independent). At the minimum scale stage, we suppose that a finite frictionless crack of a radius of 1. m suddenly appears at t =, so that rupture initiates outward spontaneously. We use the fixed grid size and time step in each stage. Table 1. Model parameters used in each stage of the simulation Simulation Stage 1st... Grid Size s m 1.6 m km Time Step t.33 ms.133 ms s Critical Slip Distance D c (L ).1 mm mm... m Reference Distance L 1. m 6. m km Magnitude at 24 t D c value is given at distance L as a reference following Equation (4). ð4þ 2of
3 slip velocity, m/s.4 slip, m shear stress, MPa Renormalization x 4 Figure 2. Snapshot of a single rupture propagation at different stages of simulation with different scales. As shown in detail at the seventh stage, slip velocity after the renormalization is very small at the beginning, however dynamic rupture propagation actually occurs in the smaller scales. Slip velocity in blue is the one renormalized from the earlier stage, whose model region is shown as a small frame. [1] We simulate over a model region of spatial grids up to 24 time steps, and then renormalize the result by a factor of 4, according to Equation (2), into the first 6 time steps of the next stage. Figures 2 shows snapshots of a single rupture propagation in different scales. After 6 times of renormalizations, we arrive in the seventh stage whose dimension is 496 (= 4 6 ) times larger than the smallest stage. We first notice that the rupture propagates in the same way at any stage with a constant rupture velocity. The spatial patterns and the value of slip velocity and shear stress do not change for any scale, while total slip actually increases. We need no spatial heterogeneity in initial stress field to generate self-similar rupture propagation. Although we have not yet constrained any condition to stop this rupture, our self-similar rupture propagation can be arrested at any size to lead a particular size of earthquake. In other words, all size of earthquake may start with the minimum size of crack. In this way, the condition D c / L produces the self-similar rupture propagation. This numerical result is also consistent with theoretical implication [e.g., Richards, 1976]. We suppose a significantly small frictionless patch to nucleate this spontaneous rupture propagation, because we do not introduce such heterogeneity in strength as observed in the experiments [e.g., Ohnaka et al., 1987]. The self- Shear S tress [MPa] t= 6 t= 6 (28,32) (3,32) (x 1, x 2) = (32,32) t= 6 t= 6 t= 6 (28,3) (3,3) (32,3) t= 6 t= 6 (28,28) t= 6 (3,28) t= 6 (32,28) Normalized Fault Slip Figure 3. Comparison of stress-slip constitutive curves before and after renormalization at the stage. Grid position (32, 32) is the closest to the center of rupture. Bold line represents the stress-slip curve calculated by BIE, and broken bold line is for the constitutive law (1). Each of the circles shows the position at time step and 6, respectively. Spontaneous process after time step 6 leads to that two lines are identical. Grey thin lines are stress-slip curve on corresponding smaller grids before renormalization. 3of
4 1. Normalized Seismic Moment Release st 3rd Time Step Figure 4. Evolution of seismic moment release function. The lines except for the first and second stages are identical. similar rupture propagation shown here implies that any particular size of earthquake needs no scale-dependent nucleation process and that all may begin from the smallest patch. [11] We have to check the stability of this renormalization process. In order to follow correctly rupture progress, it is generally required that the cohesive zone behind the rupture front is expressed by plural grids, and that slip-weakening process (u < D c ) is also followed by plural time steps. When we carry out the renormalization every 16 time steps (instead of 24), which corresponds the first 4 time steps (instead of 6) in the next stage, we get the same results, too. But if we take much shorter time step for this renormalization, the rupture cannot continue correctly. As shown in the snapshot of the seventh stage of Figure 2, we cannot distinguish a rupture front clearly at its early stage (until 3 time steps), nor even any clear evidence of slip at the beginning, although the dynamic rupture actually occurs in the smaller scale. [12] We examine the consistency between V and fracture energy G c, as we renormalized them. The stress can be independently recalculated by either of the constitutive law (1) or of the BIE. Figure 3 shows the comparison among two descriptions of stress evolution after the renormalization and the relations corresponding to the same grid area before the renormalization. As grid is closer to the hypocenter, stress accumulation is apparently lower than the actually required peak strength t p, so that slip weakening rate t b /D c becomes less steep and apparent D c becomes longer. This is because the renormalization procedure has a Time Step Dc, cm S P Hypocentral Distance Figure. Comparison of rupture propagation for different b. The curves represent the evolution of rupture front in the inplane direction. Both axes are normalized by the minimum sizes of grid size and time step. D c is also shown by side. Dotted lines correspond to the distance L b. 4of
5 smoothing effect. However this discrepancy is localized just inside of the area that already ruptured, especially at the beginning. Furthermore, the fracture energy actually consumed is in total conserved with respect to the required one. Therefore we can conclude that this smoothing effect does not influence the later stages after the renormalization. [13] Finally we show the evolution of seismic moment release function, normalized at each stage, in Figure 4. Since we give a finite frictionless crack at t =, we observe somehow its influence at the beginning. However, the relation curve for the second stage starts with an infinitesimal value, and the lines are identical after the third stage. This also infers that rupture starts with a point for any scale and propagates outward with a constant speed, keeping the self-similarity. 4. Break in Self-Similarity [14] Next we consider the case of b 6¼ 1. As shown in Figure, we keep a linear relation (b = 1) for hypocentral distance shorter than L b (48 times of the smallest grid size, corresponding to 192 m from Table 1), but suppose different b (b = 1.1,.8 and., for example) for a longer distance than L b. Figure shows the evolution of rupture front in the in-plane direction. We plot the simulation results until the 6th stage, so that the smallest scale is invisible and discrete steps appeared in larger scales correspond to one spatial grid size at each stage. We can also confirm the self-similar rupture propagation at a constant rupture velocity for the case b = 1, seen in Figure 2. [1] In the case of b = 1.1, rupture propagation is decelerated gradually according to the energy balance around the rupture front, and finally ends up a magnitude 4.2 event. This explains how the earthquake size is determined in the term of fracture energy balance. The break toward b > 1 naturally arrests the rupture propagation. On the other hand, the rupture is accelerated for b < 1. We note that slip velocity also increases according to the rupture progress so as to keep energy balance totally around the rupture front. The case of b =. has an extreme feature that rupture velocity exceeds even the shear wave velocity. This is the situation usually treated in most rupture simulations with a fixed grid size and a scale-invariable D c. Super shear rupture is rather rare in experiments and field observations. This implies that the spatial heterogeneity of fracture energy in the real field is close to the case b = 1, and tends to break toward b > 1 in some scale for the arrest of rupture.. Summary [16] We developed a numerical scheme combining a renormalization process with the boundary integral equation method, in order to accomplish a single multi-scaling earthquake rupture process in the same framework. This allows us to follow how a small earthquake grows up to a large one by a several order of its scale. We showed the self-similar rupture propagation in the case where the critical slip-weakening distance D c is proportional to the hypocentral distance L. For this case, we did not need any heterogeneity in the initial stress. We confirmed numerically that a small event becomes larger as a cascading rupture style without preparing any scale-dependent nucleation process. Then we showed that the break of the proportional relation toward b > 1 naturally arrests the rupture, so that we can generate any size of earthquake in this framework. For the case of b < 1, we also observed no-self-similar behavior in rupture process, where the rupture accelerates even at a rupture velocity higher than the shear wave velocity. [17] Acknowledgments. We thank D. J. Andrews and an anonymous reviewer for their valuable comments. This work is contribution to the DaiDaiToku Project, MEXT, Japan. References Andrews, D. J. (1976), Rupture propagation with finite stress in antiplane strain, J. Geophys. Res., 81, Ellsworth, W. L., and G. C. Beroza (199), Seismic evidence for an earthquake nucleation phase, Science, 268, Freund, L. B. (199), Dynamic fracture mechanics, Cambridge Univ. Press, Cambridge, UK. Fukuyama, E., and R. Madariaga (1998), Rupture dynamics of a planar fault in a 3D elastic medium: Rate- and Slip-weakening friction, Bull. Seismol. Soc. Am., 88, Ide, S., and M. Takeo (1997), Determination of constitutive relations of fault slip based on seismic wave analysis, J. Geophys. Res., 12(B12), Lawn, B. (1993), Fracture of brittle solids, edition, Cambridge Univ. Press, Cambridge, UK. Ohnaka, M. (23), A constitutive scaling law and a unified comprehension for frictional slip failure, shear fracture of intact rock, and earthquake rupture, J. Geophys. Res., 18(B2), 28, doi:1.129/2jb123. Ohnaka, M., Y. Kuwahara, and K. Yamamoto (1987), Constitutive relations between dynamic physical parameters near a tip of the propagating slip zone during stick-slip shear failure, Tectonophys., 144, Richards, P.G. (1976), Dynamic motions near an earthquake fault: A threedimensional solution, Bull. Seismol. Soc. Am., 66, Scholz, C. H. (22), The mechanics of earthquake and faulting, edn., Cambridge Univ. Press, Cambridge, UK. Shibazaki, B., and M. Matsu ura (1992), Spontaneous processes for nucleation, dynamic propagation, and stop of earthquake rupture, Geophys. Res. Lett., 19(12), H. Aochi, Institut de Radioprotection et de Sûreté Nucléaire, IRSN/DEI/ SARG/BERSSIN, BP17, Fontenay-aux-Roses, 92262, France. (hideo. aochi@irsn.fr) S. Ide, Department of Earth and Planetary Science, University of Tokyo, Hongo, Bunkyo-ku, Tokyo, , Japan. (ide@eps.s.u-tokyo. ac.jp) of
Effect of normal stress during rupture propagation along nonplanar faults
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 107, NO. B2, 2038, 10.1029/2001JB000500, 2002 Effect of normal stress during rupture propagation along nonplanar faults Hideo Aochi and Raúl Madariaga Laboratoire
More informationThe role of fault continuity at depth in numerical simulations of earthquake rupture
Bull. Earthq. Res. Inst. Univ. Tokyo Vol. 12,**- pp. 1/2, The role of fault continuity at depth in numerical simulations of earthquake rupture Hideo Aochi* Laboratoire de Géologie, E cole Normale Supérieure
More informationON 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 informationImprovement in the Fault Boundary Conditions for a Staggered Grid Finite-difference Method
Pure appl. geophys. 63 (6) 977 99 33 553/6/9977 DOI.7/s-6-8- Ó Birkhäuser Verlag, Basel, 6 Pure and Applied Geophysics Improvement in the Fault Boundary Conditions for a Staggered Grid Finite-difference
More informationEFFECTS OF NON-LINEAR WEAKENING ON EARTHQUAKE SOURCE SCALINGS
Extended abstract for the 11th International Conference on Fracture 2005 1 EFFECTS OF NON-LINEAR WEAKENING ON EARTHQUAKE SOURCE SCALINGS J.-P. Ampuero Geosciences Department, Princeton University, USA
More informationOn the nucleation of creep and the interaction between creep and seismic slip on rate- and state-dependent faults
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L15303, doi:10.1029/2007gl030337, 2007 On the nucleation of creep and the interaction between creep and seismic slip on rate- and state-dependent
More informationSimulation of earthquake rupture process and strong ground motion
Simulation of earthquake rupture process and strong ground motion Takashi Miyatake (1) and Tomohiro Inoue (2) (1) Earthquake Research Institute, University of Tokyo, Yayoi, Bunkyo, Tokyo, 113-0032, Japan
More informationEXAMINATION ON CONSECUTIVE RUPTURING OF TWO CLOSE FAULTS BY DYNAMIC SIMULATION
EXAMINATION ON CONSECUTIVE RUPTURING OF TWO CLOSE FAULTS BY DYNAMIC SIMULATION M. Muto 1, K. Dan 1, H. Torita 1, Y. Ohashi 1, and Y. Kase 2 1 Ohsaki Research Institute, Inc., Tokyo, Japan 2 National Institute
More information3D MODELING OF EARTHQUAKE CYCLES OF THE XIANSHUIHE FAULT, SOUTHWESTERN CHINA
3D MODELING OF EARTHQUAKE CYCLES OF THE XIANSHUIHE FAULT, SOUTHWESTERN CHINA Li Xiaofan MEE09177 Supervisor: Bunichiro Shibazaki ABSTRACT We perform 3D modeling of earthquake generation of the Xianshuihe
More informationSpectral Element simulation of rupture dynamics
Spectral Element simulation of rupture dynamics J.-P. Vilotte & G. Festa Department of Seismology, Institut de Physique du Globe de Paris, 75252 France ABSTRACT Numerical modeling is an important tool,
More informationA constitutive scaling law and a unified comprehension for frictional slip failure, shear fracture of intact rock, and earthquake rupture
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 108, NO. B2, 2080, doi:10.1029/2000jb000123, 2003 A constitutive scaling law and a unified comprehension for frictional slip failure, shear fracture of intact rock,
More informationFracture surface energy of natural earthquakes from the viewpoint of seismic observations
Bull.Earthq.Res.Inst. Univ.Tokyo Vol. 12,**- pp. /30/ Fracture surface energy of natural earthquakes from the viewpoint of seismic observations Satoshi Ide* Department of Earth and Planetary Science, University
More informationScale Dependence in the Dynamics of Earthquake Rupture Propagation: Evidence from Geological and Seismological Observations
Euroconference of Rock Physics and Geomechanics: Natural hazards: thermo-hydro-mechanical processes in rocks Erice, Sicily, 25-30 September, 2007 Scale Dependence in the Dynamics of Earthquake Rupture
More informationSimulation of tensile crack generation by three-dimensional dynamic shear rupture propagation during an earthquake
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 108, NO. B3, 2144, doi:10.1029/2001jb001738, 2003 Simulation of tensile crack generation by three-dimensional dynamic shear rupture propagation during an earthquake
More informationDistinguishing barriers and asperities in near-source ground motion
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 110,, doi:10.1029/2005jb003736, 2005 Distinguishing barriers and asperities in near-source ground motion Morgan T. Page, Eric M. Dunham, 1 and J. M. Carlson Department
More informationDi#erences in Earthquake Source and Ground Motion Characteristics between Surface and Buried Crustal Earthquakes
Bull. Earthq. Res. Inst. Univ. Tokyo Vol. 2+,**0 pp.,/3,00 Di#erences in Earthquake Source and Ground Motion Characteristics between Surface and Buried Crustal Earthquakes Paul Somerville* and Arben Pitarka
More informationRupture dynamics with energy loss outside the slip zone
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 110,, doi:10.1029/2004jb003191, 2005 Rupture dynamics with energy loss outside the slip zone D. J. Andrews U.S. Geological Survey, Menlo Park, California, USA Received
More informationfriction friction a-b slow fast increases during sliding
µ increases during sliding faster sliding --> stronger fault --> slows sliding leads to stable slip: no earthquakes can start velocity-strengthening friction slow fast µ velocity-strengthening friction
More informationHeterogeneous Coulomb stress perturbation during earthquake cycles in a 3D rate-and-state fault model
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L21306, doi:10.1029/2008gl035614, 2008 Heterogeneous Coulomb stress perturbation during earthquake cycles in a 3D rate-and-state fault
More information1.0. Shear Strength ( τ τ c )/ τ Fault Slip (w/d c ) Peak Strength (τp τ c)/ τ 0 1.2
Evolution of contacting rock surfaces and a slip- and time-dependent fault constitutive law Hideo Aochi and Mitsuhiro Matsu'ura Department of Earth and Planetary Physics, University of Tokyo, Tokyo, Japan
More informationSlip-weakening behavior during the propagation of dynamic ruptures obeying rate- and state-dependent friction laws
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 108, NO. B8, 2373, doi:10.1029/2002jb002198, 2003 Slip-weakening behavior during the propagation of dynamic ruptures obeying rate- and state-dependent friction laws
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lectures & 3, 9/31 Aug 017 www.geosc.psu.edu/courses/geosc508 Discussion of Handin, JGR, 1969 and Chapter 1 Scholz, 00. Stress analysis and Mohr Circles Coulomb Failure
More informationVerification of the asperity model using seismogenic fault materials Abstract
Verification of the asperity model using seismogenic fault materials Takehiro Hirose*, Wataru Tanikawa and Weiren Lin Kochi Institute for Core Sample Research/JAMSTEC, JAPAN * Corresponding author: hiroset@jamstec.go.jp
More informationHitoshi Hirose (1), and Kazuro Hirahara (2) Abstract. Introduction
Three dimensional simulation for the earthquake cycle at a subduction zone based on a rate- and state-dependent friction law: Insight into a finiteness and a variety of dip-slip earthquakes Hitoshi Hirose
More informationOverestimates of Earthquake Prediction Efficiency in a "Post-Prediction"
Letter J. Phys. Earth, 41, 41-43, 1993 Overestimates of Earthquake Prediction Efficiency in a "Post-Prediction" State Shozo Matsumura National Research Institute for Earth Science and Disaster Prevention,
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Surface and body forces Tensors, Mohr circles. Theoretical strength of materials Defects Stress concentrations Griffith failure
More informationHeterogeneous distribution of the dynamic source parameters of the 1999 Chi-Chi, Taiwan, earthquake
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 108, NO. B5, 2232, doi:10.1029/2002jb001889, 2003 Heterogeneous distribution of the dynamic source parameters of the 1999 Chi-Chi, Taiwan, earthquake Wenbo Zhang,
More informationEarthquake nucleation. Pablo Ampuero Caltech Seismolab
Earthquake nucleation Pablo Ampuero Caltech Seismolab How do earthquakes start? Do small and large earthquakes start differently? Predictive value of earthquake onset and foreshock sequences? Seismological
More information1/22/2015. High velocity shear experiments with possible implications to earthquake physics
High velocity shear experiments with possible implications to earthquake physics Thanks: Amir Sagy Andrew Madden, David Lockner, Einat Aharonov Harry Green, Jay Finberg, Jefferson Chang, Shalev Siman Tov
More informationCriticality of Rupture Dynamics in 3-D
Pure appl. geophys. 157 (2000) 1981 2001 0033 4553/00/121981 21 $ 1.50+0.20/0 Criticality of Rupture Dynamics in 3-D RAUL MADARIAGA 1 and KIM B. OLSEN 2 Abstract We study the propagation of seismic ruptures
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lecture 16, 9 Nov. 2017 www.geosc.psu.edu/courses/geosc508 Energy Balance of dynamic rupture Crack tip stress field Frictional Rupture Fronts Meet in the lab (522
More informationFRICTIONAL HEATING DURING AN EARTHQUAKE. Kyle Withers Qian Yao
FRICTIONAL HEATING DURING AN EARTHQUAKE Kyle Withers Qian Yao Temperature Change Along Fault Mode II (plain strain) crack rupturing bilaterally at a constant speed v r Idealize earthquake ruptures as shear
More informationSeismic and aseismic processes in elastodynamic simulations of spontaneous fault slip
Seismic and aseismic processes in elastodynamic simulations of spontaneous fault slip Most earthquake simulations study either one large seismic event with full inertial effects or long-term slip history
More informationPulse like, crack like, and supershear earthquake ruptures with shear strain localization
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2009jb006388, 2010 Pulse like, crack like, and supershear earthquake ruptures with shear strain localization Eric G.
More informationA possible mechanism of M 9 earthquake generation cycles in the area of repeating M 7 8 earthquakes surrounded by aseismic sliding
LETTER Earth Planets Space, 63, 773 777, 2011 A possible mechanism of M 9 earthquake generation cycles in the area of repeating M 7 8 earthquakes surrounded by aseismic sliding Takane Hori 1 and Shin ichi
More informationRotation 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 informationCharacterization of nucleation during laboratory earthquakes
GEOPHYSICAL RESEARCH LETTERS, VOL. 40, 5064 5069, doi:10.1002/grl.50974, 2013 Characterization of nucleation during laboratory earthquakes S. Latour, 1 A. Schubnel, 1 S. Nielsen, 2,3 R. Madariaga, 1 and
More informationFriction. Why friction? Because slip on faults is resisted by frictional forces.
Friction Why friction? Because slip on faults is resisted by frictional forces. We first describe the results of laboratory friction experiments, and then discuss the implications of the friction constitutive
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2015) 202, 2194 2206 GJI Seismology doi: 10.1093/gji/ggv284 Rupture phase diagrams for a planar fault in 3-D full-space and half-space Jiankuan Xu, 1,2,3
More informationAn intermediate deep earthquake rupturing on a dip-bending fault: Waveform analysis of the 2003 Miyagi-ken Oki earthquake
GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L24619, doi:10.1029/2004gl021228, 2004 An intermediate deep earthquake rupturing on a dip-bending fault: Waveform analysis of the 2003 Miyagi-ken Oki earthquake Changjiang
More informationCalculation of the local rupture speed of dynamically propagating earthquakes
ANNALS OF GEOPHYSICS,,,, S; doi:./ag-79 Calculation of the local rupture speed of dynamically propagating earthquakes Andrea Bizzarri Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Bologna,
More informationDevelopment of a Predictive Simulation System for Crustal Activities in and around Japan - II
Development of a Predictive Simulation System for Crustal Activities in and around Japan - II Project Representative Mitsuhiro Matsu'ura Graduate School of Science, The University of Tokyo Authors Mitsuhiro
More informationFault Representation Methods for Spontaneous Dynamic Rupture Simulation. Luis A. Dalguer
Fault Representation Methods for Spontaneous Dynamic Rupture Simulation Luis A. Dalguer Computational Seismology Group Swiss Seismological Service (SED) ETH-Zurich July 12-18, 2011 2 st QUEST Workshop,
More informationMulticycle dynamics of nonplanar strike-slip faults
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 110,, doi:10.1029/2004jb003298, 2005 Multicycle dynamics of nonplanar strike-slip faults Benchun Duan and David D. Oglesby Department of Earth Sciences, University
More informationNumerical simulation of seismic cycles at a subduction zone with a laboratory-derived friction law
Numerical simulation of seismic cycles at a subduction zone with a laboratory-derived friction law Naoyuki Kato (1), Kazuro Hirahara (2) and Mikio Iizuka (3) (1) Earthquake Research Institute, University
More informationMaterial is perfectly elastic until it undergoes brittle fracture when applied stress reaches σ f
Material is perfectly elastic until it undergoes brittle fracture when applied stress reaches σ f Material undergoes plastic deformation when stress exceeds yield stress σ 0 Permanent strain results from
More informationConditions governing the occurrence of supershear ruptures under slip-weakening friction
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 112,, doi:10.1029/2006jb004717, 2007 Conditions governing the occurrence of supershear ruptures under slip-weakening friction Eric M. Dunham
More informationAccelerating energy release prior to large events in simulated earthquake cycles: implications for earthquake forecasting
Accelerating energy release prior to large events in simulated earthquake cycles: implications for earthquake forecasting Peter Mora and David Place QUAKES, Department of Earth Sciences, The University
More information3D modeling of the cycle of a great Tohoku oki earthquake, considering frictional behavior at low to high slip velocities
GEOPHYSICAL RESEARCH LETTERS, VOL. 38,, doi:10.1029/2011gl049308, 2011 3D modeling of the cycle of a great Tohoku oki earthquake, considering frictional behavior at low to high slip velocities B. Shibazaki,
More informationChallenges in earthquake physics and source imaging
Challenges in earthquake physics and source imaging Jean-Paul Ampuero and Nadia Lapusta (Caltech Seismolab) Main goals and current issues in earthquake dynamics The source imaging inverse problem Parallels
More informationSTRESS DROP AS A RESULT OF SPLITTING, BRITTLE AND TRANSITIONAL FAULTING OF ROCK SAMPLES IN UNIAXIAL AND TRIAXIAL COMPRESSION TESTS
Studia Geotechnica et Mechanica, Vol. 37, No. 1, 2015 DOI: 10.1515/sgem-2015-0003 STRESS DROP AS A RESULT OF SPLITTING, BRITTLE AND TRANSITIONAL FAULTING OF ROCK SAMPLES IN UNIAXIAL AND TRIAXIAL COMPRESSION
More informationEffective friction law for small scale fault heterogeneity in 3D dynamic rupture
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 116,, doi:10.1029/2010jb008118, 2011 Effective friction law for small scale fault heterogeneity in 3D dynamic rupture S. Latour, 1 M. Campillo, 1 C. Voisin, 2 I. R.
More informationBrittle Deformation. Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm
Lecture 6 Brittle Deformation Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm WW Norton, unless noted otherwise Brittle deformation EarthStructure (2 nd
More informationInfluence of material contrast on fault branching behavior
GEOPHYSICAL RESEARCH LETTERS, VOL. 38,, doi:10.1029/2011gl047849, 2011 Influence of material contrast on fault branching behavior Nora DeDontney, 1 James R. Rice, 1,2 and Renata Dmowska 2 Received 20 April
More informationMODELING OF HIGH-FREQUENCY WAVE RADIATION PROCESS ON THE FAULT PLANE FROM THE ENVELOPE FITTING OF ACCELERATION RECORDS
MODELING OF HIGH-FREQUENCY WAVE RADIATION PROCESS ON THE FAULT PLANE FROM THE ENVELOPE FITTING OF ACCELERATION RECORDS Yasumaro KAKEHI 1 SUMMARY High-frequency (higher than 1 Hz) wave radiation processes
More informationDepth variation of coseismic stress drop explains bimodal earthquake magnitude-frequency distribution
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L24301, doi:10.1029/2008gl036249, 2008 Depth variation of coseismic stress drop explains bimodal earthquake magnitude-frequency distribution
More informationProperties of Dynamic Earthquake Ruptures With Heterogeneous Stress Drop
Properties of Dynamic Earthquake Ruptures With Heterogeneous Stress Drop J-P Ampuero, J Ripperger, P M Mai Institute of Geophysics, ETH Zürich, Switzerland Earthquake rupture is a notoriously complex process,
More informationSimulated and Observed Scaling in Earthquakes Kasey Schultz Physics 219B Final Project December 6, 2013
Simulated and Observed Scaling in Earthquakes Kasey Schultz Physics 219B Final Project December 6, 2013 Abstract Earthquakes do not fit into the class of models we discussed in Physics 219B. Earthquakes
More informationPUBLICATIONS. Journal of Geophysical Research: Solid Earth
PUBLICATIONS RESEARCH ARTICLE Key Points: Supershear rupture transitions studied numerically For weaker faults, rupture speeds increase continuously to supershear For stronger faults, the rupture jumps
More informationSource Process and Constitutive Relations of the 2011 Tohoku Earthquake Inferred from Near-Field Strong-Motion Data
Source Process and Constitutive Relations of the 2011 Tohoku Earthquake Inferred from Near-Field Strong-Motion Data Kunikazu Yoshida, Anatoly Petukhin & Ken Miyakoshi Geo-Research Institute, Japan Koji
More informationCoupling of the random properties of the source and the ground motion for the 1999 Chi Chi earthquake
GEOPHYSICAL RESEARCH LETTERS, VOL. 32, L08311, doi:10.1029/2004gl022202, 2005 Coupling of the random properties of the source and the ground motion for the 1999 Chi Chi earthquake Daniel Lavallée 1 and
More informationAndrews D. J. ( 1973 ) : A numerical study of tectonic stress release by underground explosions, Bull. Seism. Soc. Am., 63, No. 4, pp.
Aki K. ( 1979 ) : Characterization of barriers on an earthquake fault, J. Geophys. Res., 84, No. B11, pp. 6140 6148 Andrews D. J. ( 1973 ) : A numerical study of tectonic stress release by underground
More informationSOURCE MODELING OF RECENT LARGE INLAND CRUSTAL EARTHQUAKES IN JAPAN AND SOURCE CHARACTERIZATION FOR STRONG MOTION PREDICTION
SOURCE MODELING OF RECENT LARGE INLAND CRUSTAL EARTHQUAKES IN JAPAN AND SOURCE CHARACTERIZATION FOR STRONG MOTION PREDICTION Kimiyuki Asano 1 and Tomotaka Iwata 2 1 Assistant Professor, Disaster Prevention
More informationPotency-magnitude scaling relations for southern California earthquakes with 1.0 < M L < 7.0
Geophys. J. Int. (2002) 148, F1 F5 FAST TRACK PAPER Potency-magnitude scaling relations for southern California earthquakes with 1.0 < M L < 7.0 Yehuda Ben-Zion 1, and Lupei Zhu 2 1 Department of Earth
More informationSynthetic 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 informationWhen you are standing on a flat surface, what is the normal stress you exert on the ground? What is the shear stress?
When you are standing on a flat surface, what is the normal stress you exert on the ground? What is the shear stress? How could you exert a non-zero shear stress on the ground? Hydrostatic Pressure (fluids)
More informationJournalofGeophysicalResearch: SolidEarth
JournalofGeophysicalResearch: SolidEarth RESEARCH ARTICLE Key Points: Relations between source dimensions and corner frequencies are established Stress drop estimates are affected by differences in source
More informationThe Mechanics of Earthquakes and Faulting
The Mechanics of Earthquakes and Faulting Christopher H. Scholz Lamont-Doherty Geological Observatory and Department of Earth and Environmental Sciences, Columbia University 2nd edition CAMBRIDGE UNIVERSITY
More informationPulse-like, crack-like, and supershear earthquake ruptures with shear strain localization
JOURNAL OF GEOPHYSICAL RESEARCH, VOL.???, XXXX, DOI:10.1029/, Pulse-like, crack-like, and supershear earthquake ruptures with shear strain localization Eric G. Daub, 1 M. Lisa Manning, 1,2 and Jean M.
More informationTitle from Slip-Velocity Functions on Ear.
Title Stress-Breakdown Time and Slip-Weak from Slip-Velocity Functions on Ear Author(s) Mikumo, Takeshi; Olsen, Kim Yuji B.; Fuk Citation Bulletin of the Seismological Socie 9(): 6-8 Issue Date - URL http://hdl.handle.net//9
More informationScaling relationship of initiations for moderate to large earthquakes
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 111,, doi:10.1029/2005jb003613, 2006 Scaling relationship of initiations for moderate to large earthquakes K. Sato 1 and J. Mori 1 Received 5 January 2005; revised
More informationSeismic Efficiency, Overshoot and Enhanced Dynamic Weaking of Fractures Associated with Stimulation in Heavy Oil Reservoirs
Seismic Efficiency, Overshoot and Enhanced Dynamic Weaking of Fractures Associated with Stimulation in Heavy Oil Reservoirs Lindsey N. Meighan, Ted Urbancic and Adam Baig ESG, Kingston, Ontario Summary
More informationInstabilities and Dynamic Rupture in a Frictional Interface
Instabilities and Dynamic Rupture in a Frictional Interface Laurent BAILLET LGIT (Laboratoire de Géophysique Interne et Tectonophysique) Grenoble France laurent.baillet@ujf-grenoble.fr http://www-lgit.obs.ujf-grenoble.fr/users/lbaillet/
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Standard Solids and Fracture Fluids: Mechanical, Chemical Effects Effective Stress Dilatancy Hardening and Stability Mead, 1925
More informationEarthquakes and Seismotectonics Chapter 5
Earthquakes and Seismotectonics Chapter 5 What Creates Earthquakes? The term Earthquake is ambiguous: Applies to general shaking of the ground and to the source of the shaking We will talk about both,
More informationDepth dependence of constitutive law parameters for shear failure of rock at local strong areas on faults in the seismogenic crust
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 112,, doi:10.1029/2006jb004260, 2007 Depth dependence of constitutive law parameters for shear failure of rock at local strong areas on
More informationDoes Aftershock Duration Scale With Mainshock Size?
GEOPHYSICAL RESEARCH LETTERS, VOL.???, NO., PAGES 1 16, Does Aftershock Duration Scale With Mainshock Size? A. Ziv A. Ziv, Ben-Gurion University of the Negev, Beer-Sheva, 84105, Israel. (e-mail: zival@bgu.ac.il)
More informationGeneration mechanism of slow earthquakes: Numerical analysis based on a dynamic model with brittle ductile mixed fault heterogeneity
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 116,, doi:10.1029/2010jb008188, 2011 Generation mechanism of slow earthquakes: Numerical analysis based on a dynamic model with brittle ductile mixed fault heterogeneity
More informationDynamic source inversion for physical parameters controlling the 2017 Lesvos earthquake
Dynamic source inversion for physical parameters controlling the 2017 Lesvos earthquake F. Kostka 1, E. Sokos 2, J. Zahradník 1, F. Gallovič 1 1 Charles University in Prague, Faculty of Mathematics and
More informationVariability of earthquake nucleation in continuum models of rate-and-state faults and implications for aftershock rates
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 113,, doi:10.1029/2007jb005154, 2008 Variability of earthquake nucleation in continuum models of rate-and-state faults and implications
More informationPulse-like, crack-like, and supershear earthquake ruptures with shear strain localization
JOURNAL OF GEOPHYSICAL RESEARCH, VOL.???, XXXX, DOI:.9/, Pulse-like, crack-like, and supershear earthquake ruptures with shear strain localization Eric G. Daub, M. Lisa Manning, and Jean M. Carlson Abstract.
More informationTwo ways to think about the dynamics of earthquake ruptures
Two ways to think about the dynamics of earthquake ruptures (1) In terms of friction (2) In terms of fracture mechanics Scholz describes conditions for rupture propagation (i.e. instability) via energy
More informationScaling relations of seismic moment, rupture area, average slip, and asperity size for M~9 subduction-zone earthquakes
GEOPHYSICAL RESEARCH LETTERS, VOL. 4, 7 74, doi:1.12/grl.976, 213 Scaling relations of seismic moment, rupture area, average slip, and asperity size for M~9 subduction-zone earthquakes Satoko Murotani,
More informationIntroduction Faults blind attitude strike dip
Chapter 5 Faults by G.H. Girty, Department of Geological Sciences, San Diego State University Page 1 Introduction Faults are surfaces across which Earth material has lost cohesion and across which there
More informationEstimation of S-wave scattering coefficient in the mantle from envelope characteristics before and after the ScS arrival
GEOPHYSICAL RESEARCH LETTERS, VOL. 30, NO. 24, 2248, doi:10.1029/2003gl018413, 2003 Estimation of S-wave scattering coefficient in the mantle from envelope characteristics before and after the ScS arrival
More informationGeomechanical Analysis of Hydraulic Fracturing Induced Seismicity at Duvernay Field in Western Canadian Sedimentary Basin
Geomechanical Analysis of Hydraulic Fracturing Induced Seismicity at Duvernay Field in Western Canadian Sedimentary Basin Suvrat P Lele 1, Timothy Tyrrell 2, Ganeswara R Dasari 1, William A Symington 1
More informationGeophysical Research Letters
Geophysical Research Letters 16 JANUARY 2007 Volume 34 Number 1 American Geophysical Union Toward a more objective determination of seismic hazards Andaman Islands deformation associated with the 2004
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lecture 20, 30 Nov. 2017 www.geosc.psu.edu/courses/geosc508 Seismic Spectra & Earthquake Scaling laws. Seismic Spectra & Earthquake Scaling laws. Aki, Scaling law
More informationEffect of the Emperor seamounts on trans-oceanic propagation of the 2006 Kuril Island earthquake tsunami
GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L02611, doi:10.1029/2007gl032129, 2008 Effect of the Emperor seamounts on trans-oceanic propagation of the 2006 Kuril Island earthquake tsunami S. Koshimura, 1 Y.
More informationProbabilities for Jumping Fault Segment Stepovers
GEOPHYSICAL RESEARCH LETTERS, VOL.???, XXXX, DOI:1.129/, Probabilities for Jumping Fault Segment Stepovers Bruce E. Shaw Lamont-Doherty Earth Observatory, Columbia University, New York James H. Dieterich
More informationDynamic rupture simulation of non-planar faults with a finite-difference approach
Geophys. J. Int. (2004) 158, 939 954 doi: 10.1111/j.1365-246X.2004.02291.x Dynamic rupture simulation of non-planar faults with a finite-difference approach V. M. Cruz-Atienza and J. Virieux Géosciences
More informationEarthquake and Volcano Deformation
Earthquake and Volcano Deformation Paul Segall Stanford University Draft Copy September, 2005 Last Updated Sept, 2008 COPYRIGHT NOTICE: To be published by Princeton University Press and copyrighted, c
More informationA new hybrid numerical scheme for modelling elastodynamics in unbounded media with near-source heterogeneities
A new hybrid numerical scheme for modelling elastodynamics in unbounded media with near-source heterogeneities Setare Hajarolasvadi Ahmed E. Elbanna https://www.youtube.com/watch?v=bltx92tuwha MOTIVATION
More informationInferring fault strength from earthquake rupture properties and the tectonic implications of high pore pressure faulting
Earth Planets Space, 54, 1173 1179, 22 Inferring fault strength from earthquake rupture properties and the tectonic implications of high pore pressure faulting Stephen A. Miller Geophysics Institute, ETH-Zurich
More informationFault Specific, Dynamic Rupture Scenarios for Strong Ground Motion Prediction
Fault Specific, Dynamic Rupture Scenarios for Strong Ground Motion Prediction H. Sekiguchi Disaster Prevention Research Institute, Kyoto University, Japan Blank Line 9 pt Y. Kase Active Fault and Earthquake
More informationRate and State-Dependent Friction in Earthquake Simulation
Rate and State-Dependent Friction in Earthquake Simulation Zac Meadows UC Davis - Department of Physics Summer 2012 REU September 3, 2012 Abstract To better understand the spatial and temporal complexity
More informationVelocity contrast along the Calaveras fault from analysis of fault zone head waves generated by repeating earthquakes
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L01303, doi:10.1029/2007gl031810, 2008 Velocity contrast along the Calaveras fault from analysis of fault zone head waves generated by
More informationCharacterizing Earthquake Rupture Models for the Prediction of Strong Ground Motion
Characterizing Earthquake Rupture Models for the Prediction of Strong Ground Motion Paul Somerville URS Corporation, 566 El Dorado Street, Pasadena, CA, 91101, USA Summary The uncertainty in estimates
More informationA Prototype of Strong Ground Motion Prediction Procedure for Intraslab Earthquake based on the Characterized Source Model
A Prototype of Strong Ground Motion Prediction Procedure for Intraslab Earthquake based on the Characterized Source Model T. Iwata, K. Asano & H. Sekiguchi Disaster Prevention Research Institute, Kyoto
More informationarxiv:physics/ v2 [physics.geo-ph] 18 Aug 2003
Is Earthquake Triggering Driven by Small Earthquakes? arxiv:physics/0210056v2 [physics.geo-ph] 18 Aug 2003 Agnès Helmstetter Laboratoire de Géophysique Interne et Tectonophysique, Observatoire de Grenoble,
More information