An official website of the United States government
The strength of lithospheric plates is a central component of plate tectonics, governed by brittle processes in the shallow portion of the plate and ductile behavior in the deeper portion. We review experimental constraints on ductile deformation of olivine, the main mineral in the upper mantle and thus the lithosphere. Olivine deforms by four major mechanisms: low-temperature plasticity, dislocation creep, dislocation-accommodated grain-boundary sliding (GBS), and diffusion-accommodated grain-boundary sliding (diffusion creep). Deformation in most of the lithosphere is dominated by GBS, except in shear zones—in which diffusion creep dominates—and in the brittle-ductile transition—in which low-temperature plasticity may dominate. We find that observations from naturally deformed rocks are consistent with extrapolation of the experimentally constrained olivine flow laws to geological conditions but that geophysical observations predict a weaker lithosphere. The causes of this discrepancy are unresolved but likely reside in the uncertainty surrounding processes in the brittle-ductile transition, at which the lithosphere is strongest. ▪ Ductile deformation of the lithospheric mantle is constrained by experimental data for olivine. ▪ Olivine deforms by four major mechanisms: low-temperature plasticity, dislocation creep, dislocation-accommodated grain-boundary sliding, and diffusion creep. ▪ Observations of naturally deformed rocks are consistent with extrapolation of olivine flow laws from experimental conditions. ▪ Experiments predict stronger lithosphere than geophysical observations, likely due to gaps in constraints on deformation in the brittle-ductile transition.
more »
« less
Rheology of the lower mantle: a review
Abstract We review our current understanding of the rheological properties of the lower mantle based both on materials science and geophysics points of view. We assume a simple model of the lower mantle that is made of only two minerals: bridgmanite (Br) (Mg,Fe)SiO3and ferropericlase (Fp) (Mg,Fe)O, and address a question of (i) which mineral is weaker (lower viscosity), (ii) how does lower mantle viscosity change with depth and location, and (iii) discuss implications for shear localization. We first review plausible mechanisms of deformation based on the deformation mechanism map on the normalized stress and temperature space. We conclude that likely mechanism of deformation in the lower mantle is either diffusion creep or power-law dislocation creep. Based on this review, we discuss recently proposed models by Cordier and his group (Cordier in Nature 481:177–181, 2012; Cordier in Nature 613:303–306 , 2023) where either asthermal creep (i.e., low-temperature plasticity) or pure climb creep (not power-law dislocation creep) would play an important role. We conclude that these models are not acceptable because (1) many aspects of their models are incompatible with experimental observations and theoretical models of deformation of most materials including oxides and metals and (2) these models are not consistent with the distribution of seismic anisotropy. Hence, we focus on power-law dislocation creep and diffusion creep. We review previously published results on deformation (by dislocation creep) and diffusion, we conclude that Fp is weaker than Br. The radial (depth) depth and lateral variation of viscosity is discussed based on the estimated activation volume and estimated variation of grain-size. Geophysical studies suggest only modest depth variation of viscosity that demands relatively small activation volume (V* (< 3$$\times$$ 10–6m3/mol)). Plausible models to explain small activation volume are discussed including the role of extrinsic diffusion. Grain-size also controls viscosity if deformation is by diffusion creep. Okamoto and Hiraga (J Geophys Res, 2024. 10.1029/2023JB027803), Solomatov et al. (Phys Earth Planet Inter 129:265–282, 2002) estimated the grain-size evolution in the lower mantle based on the kinetics of grain-growth and the role of a phase transformation. In contrast, there are other papers (e.g., Paul et al. in Prog Earth Planet Sci 11:64, 2024; Rozel in Geochem Geophys Geosyst, 2012. 10.1029/2012GC004282) where grain-size distribution is estimated assuming that grain-size is controlled by dynamic recrystallization. The validity of assumption is questionable because dynamic recrystallization occurs due to deformation by dislocation creep but not by diffusion creep and the absence of seismic anisotropy indicates that diffusion creep dominates in most of the lower mantle. Finally, we review the published models of shear localization that would explain the long-term preservation of geochemical reservoirs in the lower mantle. Accepting that two minerals (Fp and Br) in the lower mantle have largely different viscosity, Ballmer et al. (Nat Geosci 10:236–240, 2017) proposed that the presence of regions of compositional difference (difference in Fp/Br ratio) leads to localized deformation (deformation mainly in the weaker regions). However, in addition to the ad hoc nature of this model, there is no strong evidence for the presence of large variation in Fp/Br in the lower mantle that makes the validity of this model questionable. There are some papers where processes of shear localization are explored without invoking the presence of regions of large rheological contrast. Thielmann et al. (Geochem Geophys Geosyst, 2020. 10.1029/2019GC008688) presented the results of theoretical study of deformation of initially homogeneous two-phase mixture (Fp and Br) and showed that deformation causes the elongation of a weak Fp that promotes shear localization. In this model, the rheological contrast between Fp and Br was assumed to be independent of strain. However, Cho and Karato (J Geophys Res 2022. 10.1029/2021JB022673 ; Phys Earth Planet Inter, 2024. 10.1016/j.pepi.2024 ) showed that when deformation is by diffusion creep, the rheological contrast increases with strain due to the evolution of stress concentration caused by grain elongation. They showed that this will promote strain weakening particularly in simple shear that would lead to shear localization. Consequently, the tendency for shear localization is stronger in their model than a model where rheological contrast is assumed to be independent of strain.
more »
« less
- Award ID(s):
- 2322719
- PAR ID:
- 10601098
- Editor(s):
- Kawahata, Hotaka
- Publisher / Repository:
- Progress in Earth and Planetary Sciences
- Date Published:
- Journal Name:
- Progress in Earth and Planetary Science
- Volume:
- 12
- Issue:
- 1
- ISSN:
- 2197-4284
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Transport of heat from the interior of the Earth drives convection in the mantle, which involves the deformation of solid rocks over billions of years. The lower mantle of the Earth is mostly composed of iron-bearing bridgmanite MgSiO 3 and approximately 25% volume periclase MgO (also with some iron). It is commonly accepted that ferropericlase is weaker than bridgmanite 1 . Considerable progress has been made in recent years to study assemblages representative of the lower mantle under the relevant pressure and temperature conditions 2,3 . However, the natural strain rates are 8 to 10 orders of magnitude lower than in the laboratory, and are still inaccessible to us. Once the deformation mechanisms of rocks and their constituent minerals have been identified, it is possible to overcome this limitation thanks to multiscale numerical modelling, and to determine rheological properties for inaccessible strain rates. In this work we use 2.5-dimensional dislocation dynamics to model the low-stress creep of MgO periclase at lower mantle pressures and temperatures. We show that periclase deforms very slowly under these conditions, in particular, much more slowly than bridgmanite deforming by pure climb creep. This is due to slow diffusion of oxygen in periclase under pressure. In the assemblage, this secondary phase hardly participates in the deformation, so that the rheology of the lower mantle is very well described by that of bridgmanite. Our results show that drastic changes in deformation mechanisms can occur as a function of the strain rate.more » « less
-
Olivine lattice preferred orientation (LPO), or texture, forms in relation to deformation mechanisms such as dislocation creep and can be observed in the upper mantle as seismic anisotropy. Olivine is also mechanically anisotropic, meaning that it responds to stresses differently depending on the direction of the stress. Understanding the interplay between anisotropic viscosity (AV) and LPO, and their role in deformation, is necessary for relating seismic anisotropy to mantle flow patterns. In this study, we employ three methods to predict olivine texture (D-Rex, MDM, and MDM+AV) in a shear box model and a subduction model. D-Rex and MDM are two representative texture development methods that have been compared before, and our results are in line with previous studies showing that textures computed by D-Rex develop faster and are stronger and more point-like than textures calculated with MDM. MDM+AV uses the same isotropic mantle stresses and particle paths as D-Rex and MDM but includes the effect of AV for texture predictions. MDM+AV predicts a texture similar to MDM with a distinct girdle-like orientation for simple shear deformation or at low strain in the subduction model. At larger strains, MDM+AV’s textures are more point-like and stronger compared to the other two methods. The effective viscosity for MDM+AV drops by up to 60% in a shear box model and can be either strengthened or weakened relative to isotropic viscosity for different regions of the subduction model experiencing different patterns of deformation. Our results emphasize the significant role of AV in olivine texture development, which could substantially affect geodynamic processes in the upper mantle.more » « less
-
null (Ed.)We present a flow law for dislocation-dominated creep in wet quartz derived from compiled experimental and field-based rheological data. By integrating the field-based data, including independently calculated strain rates, deformation temperatures, pressures, and differential stresses, we add constraints for dislocation-dominated creep at conditions unattainable in quartz deformation experiments. A Markov Chain Monte Carlo (MCMC) statistical analysis computes internally consistent parameters for the generalized flow law: urn:x-wiley:21699313:media:jgrb54871:jgrb54871-math-0001 = Aσnurn:x-wiley:21699313:media:jgrb54871:jgrb54871-math-0002e−(Q+VP)/RT. From this initial analysis, we identify different effective stress exponents for quartz deformed at confining pressures above and below ∼700 MPa. To minimize the possible effect of confining pressure, compiled data are separated into “low-pressure” (<560 MPa) and “high-pressure” (700–1,600 MPa) groups and reanalyzed using the MCMC approach. The “low-pressure” data set, which is most applicable at midcrustal to lower-crustal confining pressures, yields the following parameters: log(A) = −9.30 ± 0.66 MPa−n−r s−1; n = 3.5 ± 0.2; r = 0.49 ± 0.13; Q = 118 ± 5 kJ mol−1; and V = 2.59 ± 2.45 cm3 mol−1. The “high-pressure” data set produces a different set of parameters: log(A) = −7.90 ± 0.34 MPa−n−r s−1; n = 2.0 ± 0.1; r = 0.49 ± 0.13; Q = 77 ± 8 kJ mol−1; and V = 2.59 ± 2.45 cm3 mol−1. Predicted quartz rheology is compared to other flow laws for dislocation creep; the calibrations presented in this study predict faster strain rates under geological conditions by more than 1 order of magnitude. The change in n at high confining pressure may result from an increase in the activity of grain size sensitive creep.more » « less
-
Abstract. Viscous flow in ice is often described by the Glen flow law – anon-Newtonian, power-law relationship between stress and strain rate with astress exponent n ∼ 3. The Glen law is attributed tograin-size-insensitive dislocation creep; however, laboratory and fieldstudies demonstrate that deformation in ice can be strongly dependent ongrain size. This has led to the hypothesis that at sufficiently lowstresses, ice flow is controlled by grain boundary sliding, which explicitly incorporates the grain size dependence of ice rheology. Experimental studiesfind that neither dislocation creep (n ∼ 4) nor grain boundarysliding (n ∼ 1.8) have stress exponents that match the value ofn ∼ 3 in the Glen law. Thus, although the Glen law provides anapproximate description of ice flow in glaciers and ice sheets, itsfunctional form is not explained by a single deformation mechanism. Here weseek to understand the origin of the n ∼ 3 dependence of theGlen law by using the “wattmeter” to model grain size evolution in ice.The wattmeter posits that grain size is controlled by a balance between themechanical work required for grain growth and dynamic grain size reduction.Using the wattmeter, we calculate grain size evolution in two end-membercases: (1) a 1-D shear zone and (2) as a function of depth within anice sheet. Calculated grain sizes match both laboratory data and ice coreobservations for the interior of ice sheets. Finally, we show thatvariations in grain size with deformation conditions result in an effectivestress exponent intermediate between grain boundary sliding and dislocationcreep, which is consistent with a value of n = 3 ± 0.5 over the rangeof strain rates found in most natural systems.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1186/s40645-025-00695-6
