Physics of the Earth and Planetary Interiors 270 (2017) 195–212
Contents lists available at ScienceDirect
Physics of the Earth and Planetary Interiors
journal homepage: www.elsevier.com/locate/pepi
Formation of lithospheric shear zones: Effect of temperature on two-phase grain damage ⇑ Elvira Mulyukova a, , David Bercovici a a Yale University, Department of Geology and Geophysics, New Haven, CT, USA article info abstract
Article history: Shear localization in the lithosphere is a characteristic feature of plate tectonic boundaries, and is evident Received 17 February 2017 in the presence of small grain mylonites. Localization and mylonitization in the ductile portion of the Received in revised form 15 June 2017 lithosphere can arise when its polymineralic material deforms by a grain-size sensitive rheology in com- Accepted 24 July 2017 bination with Zener pinning, which can impede, or possibly even reverse, grain growth and thus pro- Available online 1 August 2017 motes a self-softening feedback mechanism. However, the efficacy of this mechanism is not ubiquitous and depends on lithospheric conditions such as temperature and stress. Therefore, we explore the con- Keywords: ditions under which self-weakening takes place, and, in particular, the effect of temperature and defor- Shear zones mation state (stress or strain-rate) on these conditions. In our model, the lithosphere-like polymineralic Grain damage mechanics Mylonites material is deformed in a two-dimensional simple shear driven by constant stress or strain rate. The min- eral grains evolve to a stable size, which is obtained when the rate of coarsening by normal grain growth and the rate of grain size reduction by damage are in balance. Damage involves processes by which some of the deformational energy gets transferred into surface energy. This can happen by (i) dynamic recrys- tallization (grain damage) and (ii) stretching, deforming and stirring the material interface (interface damage). The influence of temperature enters through rheological laws (which govern the rate of work and damage), grain growth kinetics, and the damage partitioning fraction, which is the fraction of defor- mational work that goes into creating new surface energy. We demonstrate that a two-phase damage model, in which the partitioning fraction depends on both temperature and roughness of the interface between the phases, can successfully match the field data, including the reported correlation of grain size and temperature, the increasing dominance of dislocation creep at higher temperatures and a large range of grain sizes observed across the depth of a single shear zone. Ó 2017 Elsevier B.V. All rights reserved.
1. Introduction is, for example, geologically evident in the presence of small grain mylonites and ultramylonites (White et al., 1980). Such subdivi- The history of deformation of the Earth’s crust and lithosphere sion of the Earth’s surface into largely undeformed tectonic plates, is recorded in the microstructure of its constitutive rocks, which separated by narrow, strongly deformed plate boundaries, suggests deform in response to the differential stress imposed by forces in that its constitutive materials experience shear localization, as can the lithosphere and underlying convecting mantle. The heat and for example be caused by some self-weakening mechanism (Kaula, the mechanical work supplied by the mantle drives deformation 1980; Bercovici, 1993; Bercovici, 1995; Tackley, 2000; Bercovici, on a multitude of scales: from the movement of the crystalline 2003; Montési, 2013; Bercovici et al., 2015). To understand how defects on the subgrain scale, associated with the solid state diffu- Earth evolved to have such an exotic style of surface dynamics, sion and dislocation creep, to the formation of mylonitic shear we need to understand the mechanisms responsible for strain zones, and further to the planetary scale deformation of the Earth’s localization across the crust and the lithosphere, including their lithosphere, manifested as plate tectonics. brittle and ductile regions. Localization in the brittle portion of The relative motion of the plates is almost entirely accommo- the Earth’s outer shell is an expected feature at low crustal and dated by highly localized deformation at the plate boundaries, as lithospheric pressures and temperatures (Kohlstedt et al., 1995). However, the ductile portion of the lithosphere is where the peak
⇑ Corresponding author. lithospheric strength occurs and is thus a bottle-neck to localiza- E-mail addresses: [email protected] (E. Mulyukova), david.bercovici@ tion; in this region, deformation is also enigmatic since it involves yale.edu (D. Bercovici). http://dx.doi.org/10.1016/j.pepi.2017.07.011 0031-9201/Ó 2017 Elsevier B.V. All rights reserved. 196 E. Mulyukova, D. Bercovici / Physics of the Earth and Planetary Interiors 270 (2017) 195–212 combined brittle-ductile behavior as well as mylonitization (1998) and Linckens et al. (2015). The reported microstructure, (Kohlstedt et al., 1995; Karato, 2008). including grain size and crystallographic orientation, together with A crystalline rock deforming viscously can experience self- the independent temperature measurements, are used to constrain weakening and localization due to a number of processes, such the dominant creep mechanism. In combination with the olivine as shear heating, presence of multiple phases (including melt) flow laws and grain growth kinematics obtained from previous and grain size reduction (via phase transformation or dynamic experimental studies (Hirth and Kohlstedt, 2003), we are able to recrystallization), to name a few (e.g., Hobbs et al., 1990; Hobbs deduce the deformation conditions, such as stress or strain rate, et al., 1990; Montési and Zuber, 2002; Montési and Zuber, 2002; for a given tectonic environment. Some of the physical parameters Karato, 2008; Karato, 2008, pp. 288–301; Bercovici et al., 2015; that enter the model are less well constrained, such as the grain Bercovici et al., 2015, and references therein). The microstructure growth activation energy and the partitioning fraction of the that develops due to the action of a given strain-localizing mecha- mechanical work that goes into creating new grain boundaries. nism can reveal whether and to what degree it is active during Here, we provide further constraints, by comparing the grain size lithospheric deformation. Such interpretations require a thorough dependence on temperature predicted by our model with that quantitative understanding of these mechanisms at the conditions reported in the observational data. representative of the lithospheric crust and mantle. The physical Field observations demonstrate a correlation between grain conditions, characterized by temperature, composition and the size and temperature (with the grain size increasing by over magnitude of the mechanical forcing (e.g., stress), vary spatially two orders of magnitude across a few hundred degrees increase and through time, so that the dependence of shear localizing mech- in temperature), increasing dominance of dislocation creep at anisms and their associated microstructure on these properties higher temperatures, and a coexistence of different deformation must also be understood. regimes across small spatial scales (Jin et al., 1998; Linckens A common geometrical quantity used to characterize the et al., 2015). These observations are used as constraints in the microstructure of a crystalline material is its grain size distribu- model presented below. tion. It can be readily measured in laboratory experiments, as well as in the natural samples collected from the field. From the physi- cal point of view, the size of the grains reflects the amount of inter- 2. A simple shear zone model nal energy that is stored as surface energy on the grain boundaries. In particular, when the material deforms by dislocation creep and A detailed mathematical description of the two phase grain undergoes dynamic recrystallization, by mobilization of a pre- damage model was presented previously (e.g.Bercovici and existing grain boundary or by the formation of subgrain bound- Ricard, 2012; Bercovici et al., 2015; Bercovici and Ricard, 2016). aries, its grain boundary area increases (Karato et al., 1980; Van We provide a brief overview of the equations here, to introduce der Wal et al., 1993). This reduces the grain size and thereby stores the governing physical processes and to highlight the main the deformational work in the form of increased grain boundary assumptions made along the way. area (Bercovici and Ricard, 2005; Austin and Evans, 2007; Ricard We consider the deformation of a polycrystalline material sub- and Bercovici, 2009; Rozel et al., 2011). In polyphase materials, ject to simple shear. The material consists of a non-dilute mixture the relative movement of grains of different phases can also lead of two immiscible fluids, or phases (e.g., a peridotite mixture with to distortion of the inter-phase boundaries, which further increases volume fractions of 40% and 60% pyroxene and olivine, respec- the surface energy (Bercovici and Ricard, 2012). tively). In the general formulation of the model, the phases have For shear localization to occur, there needs to be a mechanism different densities, viscosities and other properties, although in by which deformation can cause weakening, which would increase this study, for simplicity, the phases are assumed to have identical deformation, and which would further enhances weakening. This material properties. is, for example, the case when a material deforms by dislocation- The two-phase material is characterized by a mean grain size accommodated grain boundary sliding and undergoes dynamic within each phase, as well as the geometry (i.e., roughness or, recrystallization (Hirth and Kohlstedt, 1995; Warren and Hirth, specifically, the characteristic radius of curvature) of the interface 2006; Hansen et al., 2011; Hansen et al., 2013; Kohlstedt and that separates the two phases. The grain size and the interface Hansen, 2015; Skemer and Hansen, 2016). Another possible mech- roughness directly or indirectly influence the mechanical proper- anism, which is the focus of this study, is when a polyphase mate- ties of the material, but also evolve in response to the applied rial deforms by grain size sensitive diffusion creep and undergoes mechanical or kinematic forcing. This introduces a non-linear cou- grain size reduction due to interface-damage and pinning, which pling between the applied deformation and the effective rheology are induced by the presence of multiple phases (Bercovici and of the material. Ricard, 2012; Cross and Skemer, 2017). The grain size distribution within each phase is assumed to fol- There is abundant support from field (Jin et al., 1998; Montési low a lognormal distribution, characterized by a mean grain size. and Hirth, 2003; Warren and Hirth, 2006; Skemer et al., 2010; Moreover, the grain size evolution is assumed to be slaved to that Herwegh et al., 2011; Linckens et al., 2011; Linckens et al., 2015; of the interface roughness r, as the grain growth within each phase Hansen and Warren, 2015), experimental (Evans et al., 2001; is blocked by the presence of the other phase. This is known as the Hiraga et al., 2010; Farla et al., 2013; Platt, 2015; Cross and pinned state limit (Bercovici and Ricard, 2012), in which the grain Skemer, 2017) and theoretical studies (Bercovici and Ricard, size is proportional to the interface roughness. For the grain-size 2012; Bercovici and Ricard, 2016) that secondary phases play an distribution and phase volume fraction (e.g., 60% olivine, 40% important role in the mechanical properties of lithospheric materi- pyroxene) assumed in this study (following Bercovici and Ricard als. The goal of this study is to use this observational evidence to (2013, 2014, 2016, 2015)), the mean grain size is approximately test and further constrain the theoretical models of grain damage, p=2 times the interface roughness (the limitations of the pinned- with a particular focus on the effect of temperature, and effectively state assumption are discussed in Appendix A.). depth, on the grain-size dependent shear-localizing mechanisms The model has a two-dimensional plane geometry, with shear (see alsoLanduyt et al., 2008; Behn et al., 2009; Gueydan et al., acting parallel to the x-axis, and assumed to be infinite and uni- 2014). form in the x-direction. The temperature varies with depth, To constrain the models presented here, we use the data from increasing in the negative z-direction, while the applied stress or four different upper mantle shear zones, published in Jin et al. strain rate are constant and uniform in all directions (Fig. 1). E. Mulyukova, D. Bercovici / Physics of the Earth and Planetary Interiors 270 (2017) 195–212 197
Fig. 1. A sketch of the shear zone model used in this study. The three inserts are cartoons of some of the possible grain size distributions, with grain-colors indicating grain- size (blue is small and red is large). (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.)
In a natural shear zone, both stress and strain rate are likely to grain size by the interface roughness in the pinned state limit, is vary in space and time. Constant stress and strain-rate conditions included in the prefactor B0 entering the diffusion creep compliance approximate the end-member tectonic states for how driving B. forces or loads are supported by the lithosphere and/or mantle. The evolution of the interface roughness can be expressed as When the lithosphere is strong and provides the primary resis- (Bercovici and Ricard, 2012; Bercovici and Ricard, 2016): tance to plate-driving forces like slab-pull or effective ridge push, it is approximately deforming under constant stress conditions. dr gG f r2 ¼ I I W ð2Þ Once the lithosphere undergoes localization and weakens, the dt qrq 1 cg resistance to the driving forces is taken up by the viscous drag from the mantle. This force balance dictates the plate velocity, which where thus imposes a constant mean strain-rate on a deformation zone R m of a given width. W ¼ 2e_s ¼ 2Asnþ1 1 þ F ð3Þ The spatial uniformity of stress or strain-rate is also not guaran- r teed for field or experimental conditions. However, for the basic and case of simple shear, stress is uniform. While it is difficult to 1=m impose a uniform strain rate in simple shear when viscosity is vari- B R ¼ ð4Þ able, the mean strain rate can be controlled. F Asn 1 When deformation is driven by an applied stress, the variations is the field boundary roughness (or equivalent grain-size) at which in temperature and grain size with depth may induce vertical the strain-rates from diffusion and dislocation creep are equal. The strain rate variations and thus vertical shear e_ @v =@z. How- xz x compliances A and B are dependent on temperature, and thus the ever, our assumed width of the shear zone in the y-direction field boundary relationship between stress and interface roughness (e.g., 1kmin Linckens et al., 2015) is much smaller than the varies with temperature. width of the lithosphere in the z-direction (of order 100 km), hence The diffusion-driven process of grain growth (first term on the we assumed e_ e_ . xy xz right hand side of (2)) is counteracted by the grain size reduction In the two-dimensional simple shear set up, the only non-zero driven by deformational work (second term on the right hand side component of the symmetric strain rate tensor is @v of (2)). A more detailed mathematical and physical description of e_ ¼ e_ ¼ e_ ¼ 1 x, where v is the rate of displacement in x- xy yx 2 @y x these processes is presented in Bercovici and Ricard (2012). The direction. Similarly, there is only one non-zero component of the grain size reaches a steady state when the grain growth and grain stress-tensor, which we denote as s. The stress - strain rate rela- size reduction processes are in balance. tionship for our model with composite (diffusion and dislocation The general formulation for the fraction of deformational work creep) rheology is then given by: that goes towards distorting the interface in multiphase materials (f ) is: B I e_ ¼ Asn þ s ð1Þ rm b R ¼ m ð Þ f I f 5 where A and B are the dislocation and diffusion creep compliances I b þ b Rm r (assumed equal for both phases), respectively (see Table 1 for the list of all the parameters featured in our model, together with their where f I is given by: ranges of values), and r is the interface roughness (or effective grain f ¼ F expð C ðT=T Þkf Þ¼F f ðTÞð6Þ size). The scaling factor ðp=2Þm, which arises when we represent I 0 f ref 0 0 198 E. Mulyukova, D. Bercovici / Physics of the Earth and Planetary Interiors 270 (2017) 195–212
Table 1 Material and model properties.
Property Symbol Value/definition Dimension Rate of deformational work W 2se_ Pa s 1 Surface tension c 1Jm 2 ¼ : ; ¼ : Phase volume fraction /i /1 0 4 /2 0 6 1 : Phase distribution function g 3/1/2 0 72 Reference temperature Tref 1000 K 2 Damage partitioning fraction Rb 6 f ¼ f m f I 1 I I b þ b Rm r Temperature-dependent prefactor ð ð = Þkf Þ f I F0 exp Cf T Tref ; P Experimentally constrained constants Cf kf Cf 0
Constant prefactor F0 0 < F0 < 1 Mixing transition roughness Rm RmRF m 6 < Mixing transition prefactor Rm 0 Rm 1 Mixing transition exponent b 4
3 _ n Dislocation creep edisl ¼ As