skip to main content


Title: Resolving puzzles of the phase-transformation-based mechanism of the strong deep-focus earthquake
Abstract

Deep-focus earthquakes that occur at 350–660 km are assumed to be caused by olivine → spinel phase transformation (PT). However, there are many existing puzzles: (a) What are the mechanisms for jump from geological 10−17 − 10−15 s−1to seismic 10 − 103 s−1strain rates? Is it possible without PT? (b) How does metastable olivine, which does not completely transform to spinel for over a million years, suddenly transform during seconds? (c) How to connect shear-dominated seismic signals with volume-change-dominated PT strain? Here, we introduce a combination of several novel concepts that resolve the above puzzles quantitatively. We treat the transformation in olivine like plastic strain-induced (instead of pressure/stress-induced) and find an analytical 3D solution for coupled deformation-transformation-heating in a shear band. This solution predicts conditions for severe (singular) transformation-induced plasticity (TRIP) and self-blown-up deformation-transformation-heating process due to positive thermomechanochemical feedback between TRIP and strain-induced transformation. This process leads to temperature in a band, above which the self-blown-up shear-heating process in the shear band occurs after finishing the PT. Our findings change the main concepts in studying the initiation of the deep-focus earthquakes and PTs during plastic flow in geophysics in general.

 
more » « less
Award ID(s):
1943710
NSF-PAR ID:
10376569
Author(s) / Creator(s):
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Nature Communications
Volume:
13
Issue:
1
ISSN:
2041-1723
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Deep-focus earthquakes that occur at 350–660 km, where pressures p =12-23 GPa and temperature T =1800-2000 K, are generally assumed to be caused by olivine→spinel phase transformation, see pioneering works [1–10]. However, there are many existing puzzles: (a) What are the mechanisms for jump from geological 10−17−10−15 s−1 to seismic 10−103s−1(see [3]) strain rates? Is it possible without phase transformation? (b) How does metastable olivine, which does not completely transform to spinel at high temperature and deeply in the region of stability of spinel for over the million years, suddenly transforms during seconds and generates seismic strain rates 10−103s−1 that produce strong seismic waves? (c) How to connect deviatorically dominated seismic signals with volume-change dominated transformation strain during phase transformations [9,11]? Here we introduce a combination of several novel concepts that allow us to resolve the above puzzles quantitatively. We treat the transformation in olivine like plastic strain-induced (instead of pressure/stress-induced) and find an analytical 3D solution for coupled deformation-transformation-heating processes in a shear band. This solution predicts conditions for severe (singular) transformation-induced plasticity (TRIP) and self-blown-up deformation-transformation-heating process due to positive thermomechanochemical feedback between TRIP and strain-induced transformation. In nature, this process leads to temperature in a band exceeding the unstable stationary temperature, above which the self-blown-up shear-heating process in the shear band occurs after finishing the phase transformation. Without phase transformation and TRIP, significant temperature and strain rate increase is impossible. Due to the much smaller band thickness in the laboratory, heating within the band does not occur, and plastic flow after the transformation is very limited. Our findings change the main concepts in studying the initiation of the deep-focus earthquakes and phase transformations during plastic flow in geophysics in general. The latter may change the interpretation of different geological phenomena, e.g., the possibility of the appearance of microdiamond directly in the cold Earth crust within shear-bands [12] during tectonic activities without subduction to the mantle and uplifting. Developed theory of the self-blown-up transformation-TRIP-heating process is applicable outside geophysics for various processes in materials under pressure and shear, e.g., for new routes of material synthesis [12,13], friction and wear, surface treatment, penetration of the projectiles and meteorites, and severe plastic deformation and mechanochemical technologies. 
    more » « less
  2. The first in situ quantitative synchrotron X-ray diffraction (XRD) study of plastic strain-induced phase transformation (PT) has been performed on $\alpha-\omega$ PT in ultra-pure, strongly plastically predeformed Zr as an example, under different compression-shear pathways in rotational diamond anvil cell (RDAC). Radial distributions of pressure in each phase and in the mixture, and concentration of $\omega$-Zr, all averaged over the sample thickness, as well as thickness profile were measured. The minimum pressure for the strain-induced $\alpha-\omega$ PT, $p^d_{\varepsilon}$=1.2 GPa, is smaller than under hydrostatic loading by a factor of 4.5 and smaller than the phase equilibrium pressure by a factor of 3; it is independent of the compression-shear straining path. The theoretically predicted plastic strain-controlled kinetic equation was verified and quantified; it is independent of the pressure-plastic strain loading path and plastic deformation at pressures below $p^d_{\varepsilon}$. Thus, strain-induced PTs under compression in DAC and torsion in RDAC do not fundamentally differ. The yield strength of both phases is estimated using hardness and x-ray peak broadening; the yield strength in shear is not reached by the contact friction stress and cannot be evaluated using the pressure gradient. Obtained results open a new opportunity for quantitative study of strain-induced PTs and reactions with applications to material synthesis and processing, mechanochemistry, and geophysics. 
    more » « less
  3. Abstract

    Oceanic transform faults represent abundant yet relatively unexplored components of the hydrologic cycle in the mantle lithosphere. Current models limit fluid circulation to 600 °C, the thermal limit of earthquakes recorded by teleseismic surveys. However, recent ocean‐bottom seismic surveys have located earthquakes at depths corresponding to >1000 °C in modeled thermal structure. To constrain the depth extent of brittle deformation and fluid infiltration, we analyzed peridotite mylonites dredged from the Shaka Transform Fault, Southwest Indian Ridge. Samples range from high strain mylonites that preserve ductile microstructures to lower strain mylonites that are fractured and overprinted by hydrothermal alteration. Microstructural analysis of the high strain samples reveals brittle deformation of pyroxene concomitant with ductile deformation of olivine and growth of amphibole. Porphyroclasts preserve healed fractures filled with fluid inclusions, implying repeated episodes of fracture, fluid infiltration, and healing. The association of hydration features with brittle structures points to seawater, rather than melt, as the fluid source. Textural analysis indicates that strain localization was initiated by grain boundary pinning and that olivine grain size was reduced to ~1 μm in the presence of amphibole. Comparing the amphibole stability field to thermometry estimates for the limit of recrystallization suggests that fluid flow extended to ~650–850 °C. Our results indicate that the hydrologic cycle extends past the brittle‐ductile transition and promotes strain localization via hydrolytic weakening and hydration reactions. We propose that seawater infiltration on oceanic transform faults is driven by the seismic cycle and represents a first order control on the rheology of the oceanic lithosphere.

     
    more » « less
  4. Abstract

    Experimentally quantifying the viscoplastic rheology of olivine at the high stresses and low temperatures of the shallow lithosphere is challenging due to olivine's propensity to deform by brittle mechanisms at these conditions. In this study, we use microscale uniaxial compression tests to investigate the rheology of an olivine single crystal at room pressure and temperature. Pillars with nominal diameters of 1.25 μm were prepared using a focused ion beam milling technique and were subjected to sustained axial stresses of several gigapascal. The majority of the pillars failed after dwell times ranging from several seconds to a few hours. However, several pillars exhibited clear evidence of plastic deformation without failure after 4–8 hr under load. The corresponding creep strain rates are estimated to be on the order of 10−6to 10−7 s−1. The uniaxial stresses required to achieve this deformation (4.1–4.4 GPa) are in excellent agreement with complementary data obtained using nanoindentation techniques. Scanning transmission electron microscopy observations indicate that deformation occurred along amorphous shear bands within the deformed pillars. Electron energy loss spectroscopy measurements revealed that the bands are enriched in Fe and depleted in Mg. We propose that inhomogeneities in the cation distribution in olivine concentrate stress and promote the amorphization of the Fe‐rich regions. The time dependence of catastrophic failure events suggests that the amorphous bands must grow to some critical length scale to generate an unstable defect, such as a shear crack.

     
    more » « less
  5. Abstract

    Flexure occurs on intermediate geologic timescales (∼1 Myr) due to volcanic‐island building at the Island of Hawaii, and the deformational response of the lithosphere is simultaneously elastic, plastic, and ductile. At shallow depths and low temperatures, elastic deformation transitions to frictional failure on faults where stresses exceed a threshold value, and this complex rheology controls the rate of deformation manifested by earthquakes. In this study, we estimate the seismic strain rate based on earthquakes recorded between 1960 and 2019 at Hawaii, and the estimated strain rate with 10−18–10−15s−1in magnitude exhibits a local minimum or neutral bending plane at 15 km depth within the lithosphere. In comparison, flexure and internal deformation of the lithosphere are modeled in 3D viscoelastic loading models where deformation at shallow depths is accommodated by frictional sliding on faults and limited by the frictional coefficient (μf), and at larger depths by low‐temperature plasticity and high‐temperature creep. Observations of flexure and the seismic strain rate are best‐reproduced by models withμf = 0.3 ± 0.1 and modified laboratory‐derived low‐temperature plasticity. Results also suggest strong lateral variations in the frictional strength of faults beneath Hawaii. Our models predict a radial pattern of compressive stress axes relative to central Hawaii consistent with observations of earthquake pressure (P) axes. We demonstrate that the dip angle of this radial axis is essential to discerning a change in the curvature of flexure, and therefore has implications for constraining lateral variations in lithospheric strength.

     
    more » « less