0
Research Papers

Complex Dynamics of Spring-Block Earthquake Model Under Periodic Parameter Perturbations

[+] Author and Article Information
Srđan Kostić

Department of Geology,
University of Belgrade Faculty
of Mining and Geology,
Đušina 7,
Belgrade 11000 Serbia
e-mail: srdjan.kostic@rgf.bg.ac.rs

Nebojša Vasović

Department of Applied Mathematics,
University of Belgrade Faculty
of Mining and Geology,
Đušina 7,
Belgrade 11000 Serbia
e-mail: nebojsa.vasovic@rgf.bg.ac.rs

Igor Franović

Department of Theoretical Mechanics,
Statistical Physics, and Electrodynamics,
University of Belgrade Faculty of Physics,
Studentski Trg 12,
Belgrade 11000 Serbia
e-mal: igor.franovic@gmail.com

Kristina Todorović

Department of Physics and Mathematics,
University of Belgrade Faculty of Pharmacy,
Vojvode Stepe 450,
Belgrade 11000 Serbia
e-mail: kisi@pharmacy.bg.ac.rs

1Corresponding author.

Contributed by the Design Engineering Division of ASME for publication in the JOURNAL OF COMPUTATIONAL AND NONLINEAR DYNAMICS. Manuscript received June 10, 2013; final manuscript received December 11, 2013; published online March 13, 2014. Assoc. Editor: Ahmet S. Yigit.

J. Comput. Nonlinear Dynam 9(3), 031019 (Mar 13, 2014) (10 pages) Paper No: CND-13-1134; doi: 10.1115/1.4026259 History: Received June 10, 2013; Revised December 11, 2013

A simple model of earthquake nucleation that may account for the onset of chaotic dynamics is proposed and analyzed. It represents a generalization of the Burridge–Knopoff single-block model with Dieterich–Ruina's rate- and state-dependent friction law. It is demonstrated that deterministic chaos may emerge when some of the parameters are assumed to undergo small oscillations about their equilibrium values. Implementing the standard numerical methods from the theory of dynamical systems, the analysis is carried out for the cases having one or two periodically variable parameters, such that the appropriate bifurcation diagrams, phase portraits, power spectra, and the Lyapunov exponents are obtained. The results of analysis indicate two different scenarios to chaos. On one side, the Ruelle–Takens–Newhouse route to chaos is observed for the cases of limit amplitude perturbations. On the other side, when the angular frequency is assumed constant for the value near the periodic motion of the block in an unperturbed case, variation of oscillation amplitudes probably gives rise to global bifurcations, with immediate occurrence of chaotic behavior. Further analysis shows that chaotic behavior emerges only for small oscillation frequencies and higher perturbation amplitudes when two perturbed parameters are brought into play. If higher oscillation frequencies are assumed, no bifurcation occurs, and the system under study exhibits only the periodic motion. In contrast to the previous research, the onset of chaos is observed for much smaller values of the stress ratio parameter. In other words, even the relatively small perturbations of the control parameters could lead to deterministic chaos and, thus, to instabilities and earthquakes.

Copyright © 2014 by ASME
Your Session has timed out. Please sign back in to continue.

References

Ohnaka, M., 1993, “Critical Size of the Nucleation Zone of Earthquake Rupture Inferred From Immediate Foreshock Activity,” J. Phys. Earth, 41, pp. 45– 56. [CrossRef]
Ellsworth, W. L., and Beroza, G. C., 1995, “Seismic Evidence for an Earthquake Nucleation Phase,” Science, 268, pp. 851–855. [CrossRef] [PubMed]
Dodge, D. A., Beroza, G. C., and Ellsworth, W. L., 1996, “Detailed Observations of California Foreshock Sequences: Implications for the Earthquake Initiation Process,” J. Geophys. Res., 101, pp. 22,371–22,392. [CrossRef]
Abercrombie, R., and Mori, J., 1994, “Local Observations of the Onset of a Large Earthquake: 28 June 1992 Landers, California,” B. Seismol. Soc. Am.84, pp. 725–734.
Mori, J., and Kanamori, H., 1996, “Initial Rupture of Earthquakes in the 1995 Ridgecrest, California Sequence,” Geophys. Res. Lett., 23, pp. 2340–2437. [CrossRef]
Geller, R. J., 1976, “Body Force Equivalents for Stress-Drop Seismic Sources,” B. Seismol. Soc. Am., 66, pp. 1801–1804.
Ben-Menahem, A., and Singh, S. J., 2000, Seismic Waves and Sources, Springer-Verlag, New York.
Aki, K., and Richards, P. G., 2002, Quantitative Seismology, Theory and Methods, Freeman, San Francisco.
Gibowicz, S. J., and Kijko, A., 1994, An Introduction to Mining Seismology, Academic Press, Salt Lake City, UT.
Brace, W., and Byerlee, J., 1966, “Stick-Slip as a Mechanism for Earthquakes,” Science, 153(3739), pp. 990–992. [CrossRef] [PubMed]
Burridge, R., and Knopoff, L., 1967, “Model and Theoretical Seismicity,” B. Seismol. Soc. Am.57, pp. 341–371.
Carlson, J. M., and Langer, J. S., 1989, “Mechanical Model of an Earthquake Fault,” Phys. Rev. A, 40, pp. 6470–6484. [CrossRef] [PubMed]
Elmer, F.-J., 1997, “Self-Organized Criticality With Complex Scaling Exponents in the Train Model,” Phys. Rev. E, 56, pp. R6225–R6228. [CrossRef]
Galvanetto, U., 2002, “Some Remarks on the Two-Block Symmetric Burridge–Knopoff Model,” Phys. Lett. A, 293, pp. 251–259. [CrossRef]
Popp, K., and Stelter, P., 1990, “Stick-Slip Vibrations and Chaos,” Phil. Trans. R. Soc. Lond. A, 332, pp. 89–105. [CrossRef]
Galvanetto, U., 2001, “Some Discontinuous Bifurcations in a Two-Block Stick-Slip System,” J. Sound Vib., 248, pp. 653–669. [CrossRef]
Ranjith, K., and Rice, J. R., 1999, “Stability of Quasi-Static Slip in a Single Degree of Freedom Elastic System With Rate and State Dependent Friction,” J. Mech. Phys. Solids, 47, pp. 1207–1218. [CrossRef]
Scholz, C. H., 2002, The Mechanics of Earthquakes and Faulting, Cambridge University Press, Cambridge.
Montagne, R., and Vasconcelos, G. L., 2004, “Complex Dynamics in a One-Block Model for Earthquakes,” Physica A, 342, pp. 178–185. [CrossRef]
Erickson, B., Birnir, B., and Lavallee, D., 2008, “A Model for Aperiodicity in Earthquakes,” Nonlinear Proc. Geoph., 15, pp. 1–12. [CrossRef]
Pomeau, Y., and Le Berre, M., 2011, “Critical Speed-Up vs Critical Slow-Down: A New Kind of Relaxation Oscillation With Application to Stick-Slip Phenomena,” arXiv: 1107.3331v1.
Ruina, A. L., 1983, “Slip Instability and State Variable Friction Laws,” J. Geophys. Res., 88, pp. 10359–10370. [CrossRef]
Tullis, T. E., and Weeks, J. D., 1986, “Constitutive Behavior and Stability of Frictional Sliding in Granite,” Pure Appl. Geophys., 124, pp. 383–414. [CrossRef]
Dieterich, J. H., 1979, “Modeling of Rock Friction, 1. Experimental Results and Constitutive Equations,” J. Geophys. Res., 84, pp. 2161–2168. [CrossRef]
Madariaga, R., and Olsen, K. B., 2002, “Earthquake Dynamics,” International Handbook of Earthquake and Engineering Seismology, Part A., W. H. K.Lee, H.Kanamori, P. C.Jennings, and C.Kisslinger, eds., Academic Press, London, pp. 176–194.
Kawamura, H., Hatano, T., Kato, N., Biswas, S., and Chakrabarti, B. K., 2012, Statistical Physics of Fracture, Friction and Earthquake, Rev. Mod. Phys., 84, pp. 839–884. [CrossRef]
Kostić, S., Franović, I., Todorović, K., and Vasović, N., 2013, “Friction Memory Effect in Complex Dynamics of Earthquake Model,” Nonlinear Dynam., 73(3), pp. 1933–1943. [CrossRef]
Gomberg, J., Blanpied, M., and Beeler, N., 1997, “Transient Triggering of Near and Distant Earthquakes,” B. Seismol. Soc. Am., 87, pp. 294–309.
Perfettini, H., Schmittbuhl, J., and Cochard, A., 2003, “Shear and Normal Load Perturbations on a Two-Dimensional Continuous Fault: 2. Dynamic Triggering,” J. Geophys. Res.108(2409), pp. 1–16.
Burić, N., and Vasović, N., 1999, “A Simple Model of the Chaotic Immune Response,” Chaos Soliton. Fract.10, pp. 1185–1192. [CrossRef]
Dieterich, J. H., and Kilgore, B. D., 1994, “Direct Observation of Frictional Contacts: New Insights for State Dependent properties,” Pure Appl. Geophys., 143, pp. 283–302. [CrossRef]
Rice, J. R., Lapusta, N., and Ranjith, K., 2001, “Rate and State Dependent Friction and the Stability of Sliding Between Elastically Deformable Solids,” J. Mech. Phys. Solids, 49, pp. 1865–1898. [CrossRef]
Gomberg, J., Bodin, P., Savage, W., and Jackson, M. E., 1995, “Landslide Faults and Tectonic Faults, Analogs?—The Slumgullion Earthflow, Colorado,” Geology, 23, pp. 41–44. [CrossRef]
Perfettini, H., Schmittbuhl, J., Rice, J. R., and Cocco, M. J., 2001, “Frictional Response Induced by Time-Dependent Fluctuations of the Normal Loading,” Geophys. Res., 106, pp. 13455–13472. [CrossRef]
Gomberg, J., Beeler, N. M., Blanpied, M. L., and Bodin, P., 1998, “Earthquake Triggering by Transient and Static Deformations,” J. Geophys. Res., 103, pp. 24411–24426. [CrossRef]
Parsons, T., 2005, “A Hypothesis for Delayed Dynamic Earthquake Triggering,” Geophys. Res. Lett., 32(4), pp. 1–4. [CrossRef]
Wolf, A., Swift, J., Swinney, H., and Vastano, J., 1985, “Determining Lyapunov Exponents From a Time Series,” Physica D, 16, pp. 285–317. [CrossRef]
Rosenstein, M. T., Collins, J. J., and De Luca, C. J., 1993, “A Practical Method for Calculating Largest Lyapunov Exponents From Small Data Sets,” Physica D, 65, pp. 117–134. [CrossRef]
Erickson, B. A., Birnir, B., and Lavallée, D., 2011, “Periodicity, Chaos and Localization in a Burridge–Knopoff Model of an Earthquake With Rate-and-State Friction,” Geophys. J. Int., 187, pp. 178–198. [CrossRef]

Figures

Grahic Jump Location
Fig. 2

Bifurcations of the system (2) (or (3) if δi = 0 is set) under the variation of one of the parameters ε (1) and ξ (2). Orbital diagram is constructed for the section with plane θ = 1, and calculation step 0.01, showing the asymptotic dynamics after 8 × 106 time units. At each instance, the parameters held constant are awarded values that admit the equilibrium point, ε = 0.2, ξ = 0.6, and γ = 0.8.

Grahic Jump Location
Fig. 1

The Burridge–Knopoff block and spring model, represented by a slider coupled through a spring to a loader plate, which moves with velocity V

Grahic Jump Location
Fig. 3

Attractors of the system (3) (or (4) with δi = 0) in parameter plane ε-ξ. The remaining parameters are held fixed at values admitting the equilibrium point, as in Fig. 2 Corresponding time series and phase portraits for points 1 and 2 are shown in Fig. 4 EQ and PM are abbreviations for equilibrium state and periodic motion, respectively.

Grahic Jump Location
Fig. 4

Temporal evolution of variable V and the appropriate phase portraits for: (1) ε = 0.1, ξ = 1.0, and γ = 0.8 (equilibrium state); (2) ε = 0.5, ξ = 1, and γ = 0.8 (periodic motion)

Grahic Jump Location
Fig. 5

Single peak in power spectrum indicates the oscillatory behavior of the model. Parameter values are the same as in Fig. 4(2).

Grahic Jump Location
Fig. 6

Bifurcations of the system (3) under the variation of the parameter ε. Orbital diagram is constructed for the section with plane θ = 1, and calculation step ωε = 0.01, showing the asymptotic dynamics after 8 × 106 time units. At each instance, the parameters held constant are awarded the values near the equilibrium point but admitting the limit cycle, ε = 0.4 (δε = 0.4), ξ = 0.5, and γ = 0.8. Corresponding time series and phase portraits for periodic and chaotic motion are shown in Fig. 7.

Grahic Jump Location
Fig. 7

Temporal evolution of variable V and the appropriate phase portraits for: (1) ωε = 0.9 (periodic motion); (2) ωε = 0.2 (deterministic chaos). At each instance, the parameters held constant are awarded the values near the equilibrium point, but admitting the limit cycle, ε = 0.4 (δε = 0.4), ξ = 0.5, and γ = 0.8.

Grahic Jump Location
Fig. 8

(1) Fourier power spectrum of periodic motion (first peak is for fundamental frequency, other peaks represent harmonics); (2) the broadband noise in the Fourier power spectrum indicates the onset of deterministic chaos. Parameter values are the same as in Figs. 7(1) and 7(2), respectively.

Grahic Jump Location
Fig. 9

Maximal Lyapunov exponent converges well to λmax = 0.072, confirming the onset of deterministic chaos. The parameter values are identical to those in Fig. 7(2).

Grahic Jump Location
Fig. 10

Bifurcations of the system (3) under the variation of the parameter ξ. Orbital diagram is constructed for the section with plane θ = 1, and calculation step ωξ = 0.01, showing the asymptotic dynamics after 4 × 106 time units. At each instance, the parameters held constant are awarded values near the equilibrium point, but admitting the limit cycle, ε = 0.4, ξ = 0.5 (δξ = 0.4), and γ = 0.8. Corresponding time series and phase portraits for quasi-periodic motion are shown in Fig. 11.

Grahic Jump Location
Fig. 11

Temporal evolution of variable V (1) and the appropriate phase portrait (2) for ωξ = 0.65 (quasi-periodic motion). At each instance, the parameters held constant are awarded the values near the equilibrium point, but admitting the limit cycle, ε = 0.4 (δε = 0.4), ξ = 0.5, and γ = 0.8.

Grahic Jump Location
Fig. 12

Several incommensurate frequencies in Fourier power spectrum indicate quasi-periodic motion. Parameter values are the same as in Fig. 11.

Grahic Jump Location
Fig. 13

Bifurcations of the system (4), for different perturbation amplitudes δε and δξ: (1) ωε=0.2; (2) ωξ=0.3. In both cases, orbital diagram is constructed for the section with plane θ = 1, and calculation step ωε = ωξ = 0.01, showing the dynamics after 8 × 106 and 4 × 106 time units, respectively. At each instance, the parameters held constant are awarded values near the equilibrium point, but admitting the limit cycle, ε = 0.4, ξ = 0.5, and γ = 0.8.

Grahic Jump Location
Fig. 16

Temporal evolution of variable V and the appropriate phase portraits for points 1 and 2 in Fig. 15: (1) δε = 0.1, δξ = 0.05 (periodic motion); (2) δε = 0.2, δξ = 0.2 (deterministic chaos). At each instance, the parameters held constant are awarded values near the equilibrium point but admitting the limit cycle, ε = 0.4, ξ = 0.5, and γ = 0.8. Values of the angular frequency are chosen to be near the angular frequency of the block in an unperturbed state, admitting the onset of deterministic chaos for a single parameter perturbation: ωε=0.2, ωξ=0.3.

Grahic Jump Location
Fig. 17

(1) Two commensurate peaks in power spectrum (first peak for the fundamental frequency, and the second peak for the harmonic) imply the periodic motion (2). The broadband noise in the Fourier power spectrum indicates the chaotic behavior of the system. The parameter values are the same as in Figs. 16(1) and 16(2), respectively.

Grahic Jump Location
Fig. 18

Calculation of maximal Lyapunov exponent for a time series in Fig. 17(2): (1) λmax = 0.009 (method of Wolf et al. [37]); (2) λmax ≈ 0.009 (method of Rosenstein et al. [38]). Effective expansion rate S(Δn) represents the average of the logarithm of Di(Δn), defined as the average distance of all nearby trajectories to the reference trajectory as a function of the relative time Δn. The slope of dashed lines indicating the predominant slope of S(Δn) in dependence on Δndt presents a robust estimate for the maximal Lyapunov exponent. The results are determined for 1000 reference points and neighboring distance ε = 0.1–0.5.

Grahic Jump Location
Fig. 14

Bifurcations of the system (3), for different perturbation amplitudes δε and δξ: (1) ωε=0.9; (2) ωξ=0.9. In both cases, orbital diagram is constructed for the section with plane θ = 1, and calculation step ωε = ωξ = 0.01 showing the dynamics after 8 × 106 time units. At each instance, the parameters held constant are awarded values near the equilibrium point but admitting the limit cycle, ε = 0.4, ξ = 0.5, and γ = 0.8.

Grahic Jump Location
Fig. 15

Attractors of the system (3) in parameter plane δε - δξ, for the frequencies near the frequency of the block in unperturbed state, admitting chaos due to a single parameter perturbation (ωε=0.2, ωξ=0.3). Diagram is constructed for the grid 0.01 × 0.01. At each instance, the parameters held constant are awarded the values near the equilibrium point but admitting the limit cycle, ε = 0.4, ξ = 0.5, and γ = 0.8. P and C are abbreviations for periodic motion and deterministic chaos, respectively. Corresponding time series and phase portraits for points 1 and 2 are shown in Fig. 16.

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In