Periodic Viscous Shear Heating Instability in Fine-Grained Shear Zones: Possible Mechanism for Intermediate Depth Earthquakes and Slow Earthquakes?

Localized ductile shear zones with widths of cm to m are observed in exposures of Earth's shallow mantle (e.g., Kelemen & Dick JGR 95; Vissers et al. Tectonophys 95) and dredged from oceanic fracture zones (e.g., Jaroslow et al. Tectonophys 96). These are mylonitic (grain size 10 to 100 microns) and record mineral cooling temperatures from 1100 to 600 C. Pseudotachylites in a mantle shear zone show that shear heating temperatures can exceed the mantle solidus (e.g., Obata & Karato Tectonophys 95). Simple shear, recrystallization, and grain boundary sliding all decrease the spacing between pyroxenes, so olivine grain growth at lower stress is inhibited; thus, once formed, these shear zones do not "heal" on geological time scales. Reasoning that grain-size sensitive creep will be localized within these shear zones, rather than host rocks (grain size 1 to 10 mm), and inspired by the work of Whitehead & Gans (GJRAS 74), we thought these might undergo repeated shear heating instabilities. In this view, as elastic stress increases, the shear zone weakens via shear heating; rapid deformation of the weak shear zone releases most stored elastic stress; lower stress and strain rate coupled with diffusion of heat into host rocks leads to cooling and strengthening, after which the cycle repeats. We constructed a simple numerical model incorporating olivine flow laws for dislocation creep, diffusion creep, grain boundary sliding, and low T plasticity. We assumed that viscous deformation remains localized in shear zones, surrounded by host rocks undergoing elastic deformation. We fixed the velocity along one side of an elastic half space, and calculated stress due to elastic strain. This stress drives viscous deformation in a shear zone of specified width. Shear heating and thermal diffusion control temperature evolution in the shear zone and host rocks. A maximum of 1400 C (where substantial melting of peridotite would occur) is imposed. Grain size evolves during dislocation creep and grain boundary sliding as a function of stress and strain, and undergoes diffusive growth during diffusion creep. For strain rates ca E-13 per second and initial temperatures ca 600 to 850 C, this model produces periodic viscous shear heating events with periods of 100's of years. Strain rates during these events approach 1 per second as temperatures reach 1400 C, so future models will incorporate inertial terms in the stress. Cooling between events returns the shear zone almost to its initial temperature, but ultimately shear zone temperature between events exceeds 850 C resulting in stable viscous creep. Back of the envelope calculations based on model results support the view that viscous deformation in both shear zone and host will be mainly via grain-size sensitive creep, and thus deformation will remain localized in shear zones. Similarly, we infer that inertial terms will remain small. Future models will test and quantify these inferences. The simple model described above provides an attractive explanation for intermediate-depth earthquakes, especially those in subduction zones that occur in a narrow thermal window (e.g., Hacker et al JGR 2003). We think that a "smoother"periodic instability might be produced via the same mechanism in weaker materials, which could provide a viscous mechanism for some slow earthquakes. By AGU, we will construct a second, simple model using quartz rheology to investigate this. Finally, coupling of viscous shear heating instabilities in the shallow mantle with brittle stick-slip deformation in the weaker, overlying crust may influence earthquake frequency.