skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on February 1, 2026

Title: Numerical Investigation Into Mechanical Behavior of Metastable Olivine During Phase Transformation: Implications for Deep‐Focus Earthquakes
Abstract One hypothesized mechanism that triggers deep‐focus earthquakes in oceanic subducting slabs below ∼300 km depth is transformational faulting due to the olivine‐to‐spinel phase transition. This study uses finite element modeling to investigate phase transformation‐induced stress redistribution and material weakening in olivine. A thermodynamically consistent constitutive model is developed to capture the evolution of phase transformation in olivine under different pressure and temperature conditions. The overall numerical model enables considering multiscale material features, including the polycrystalline structure, mesoscale heterogeneity, and various phases or variants of phases at the microscopic level, and accounts for viscoplastic behaviors with thermo‐mechanical coupling effects. The model is validated with several benchmarks, including a phase diagram of phase transformation from olivine to spinel. The validated model is used to study the interactive behaviors between defects (heterogeneity) and phase transformation. The simulation results reveal that spinel formation under pressure initiates near inclusions and along the grain boundaries, consistent with experimental observations. At lower temperatures, the transformation leads to the formation of thin conjugate bands of spinel diagonal to the compression loading direction. Local stress analysis along these bands also suggests the initiation of faulting. In contrast, the numerical results at higher transformation rates show that significant spinel formation occurs over a larger area at elevated temperatures, leading to ductile behavior, which agrees with experimental findings. Numerical simulation of multiple inclusions under confined pressure also shows the formation of a network of spinel bands resembling phase‐transformation patterns observed in the laboratory experiments. Additionally, stress softening patterns due to phase transformation are similar to experimental observations.  more » « less
Award ID(s):
1925920 1926011
PAR ID:
10569913
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
American Geophysical Union
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
130
Issue:
2
ISSN:
2169-9313
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. Abstract Deep‐focus earthquakes at 350–660 km are presumably caused by olivine‐spinel phase transformation (PT). This cannot, however, explain the observed high seismic strain rate, which requires PT to complete within seconds, while metastable olivine does not transform for over a million years. Recent theory quantitatively describes how severe plastic deformations (SPD) can solve this dilemma but lacking experimental proof. Here, we introduce dynamic rotational diamond anvil cell with rough diamond anvils to impose SPD on San Carlos olivine. While olivine never transformed to spinel at room temperature, we obtained reversible olivine‐ringwoodite PT under SPD at 15–28 GPa within tens of seconds. The PT pressure reduces with increasing dislocation density, microstrain, plastic strain, and decreasing crystallite size. Results demonstrate a new strain‐induced PT mechanism compared to a pressure/temperature‐induced one. Combined with SPD during olivine subduction, this mechanism can accelerate olivine‐ringwoodite PT from millions of years to timescales relevant to earthquakes. 
    more » « less
  3. Composite mantle xenoliths from the Cima Volcanic Field (CA, USA) contain a variety of melt (now glassy) inclusions hosted within mantle phases. The compositions and textures of these melt inclusions have the po- tential to constrain their trapping processes, melt sources, and the rates of ascent of their parent xenoliths. Here we focus on unusual spinel-hosted melt inclusions from one composite xenolith, reporting glass and daughter mineral compositions and textures and attempting to reconstruct inclusion bulk compositions. The xenolith contains spinel-hosted melt inclusions in its harzburgite, olivine-websterite and lherzolite layers; there are none in its orthopyroxenite layer. The glass compositions and reconstructed bulk compositions of the partly-crystallized inclusions correspond to alkaline intermediate melts, mostly trachyandesites. Such melts are most likely to be generated and trapped by vapor-undersaturated phlogopite or amphibole dehydration melting to an assemblage of liquid + spinel + olivine ± pyroxenes. We modeled the near-liquidus phase relations of the inclusion bulk compositions and noted the closest approach of each inclusion to simultaneous saturation with spinel and either phlogopite or amphibole, resulting in estimated trapping pressures of ~0.5–1.5 GPa and temperatures of ~1000–1100 ◦C. The large size of the hosting spinel grains suggests a slow process associated with these breakdown reactions, probably thinning of the lithosphere and steepening of the geotherm during regional extension. A linear correlation between the vesicle area and inclusion area (as proxies for volume) suggests an in-situ exsolution process from melts of relatively uniform volatile initial contents, consistent with trapping of vapor- undersaturated melts that later exsolve vapor during cooling and daughter crystal growth. A negative correla- tion between the glass content in melt inclusions and the size of the inclusion itself suggests a control on the degree of crystallinity with the size. There appears to be a two-stage cooling history captured by the inclusions, forming first prismatic daughter crystals and large round vesicles at the wall of the inclusion, followed by quenching to form a mat of fine crystallites and small vesicles in most inclusions. We connect the final quench to rapid ascent of the xenolith in its host melt, which also triggered partial breakdown of remaining amphibole to fine glassy symplectites. 
    more » « less
  4. 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
  5. Abstract The application of melt inclusions (MI) to infer magmatic processes assumes the MI have remained as constant mass, constant volume systems since the time of trapping. Understanding the effects of both compositional and volumetric re‐equilibration is key for the interpretation of MI data. Although the re‐equilibration behavior MI in quartz and olivine has been studied in some detail, the process is less understood for other MI host phases such as plagioclase, a common phase in igneous rocks. A MI can re‐equilibrate when it experiences pressure and temperature (PT) conditions that differ from formation PT conditions. During laboratory heating, irreversible MI expansion may occur. As a result, the internal pressure within the MI decreases, resulting in chemical and structural changes to the MI and host. We present results of heating experiments on plagioclase‐hosted MI designed to induce volumetric re‐equilibration. The experiments consisted of incrementally heating the MI to temperatures above the homogenization temperatures. At ∼40°C above, the temperature at which the daughter minerals melted, irreversible volume expansion lowered the pressure in the MI, and led to exsolution of CO2into vapor bubbles. With each additional few degrees of heating, additional episodes of CO2exsolution, bubble nucleation and expansion of the vapor bubblesoccurred. Re‐equilibration of MI in plagioclase occurred through a combination of ductile and brittle deformation of the host surrounding the MI, whereas previous studies have shown that MI in olivine re‐equilibrate dominantly through ductile deformation associated with movement along dislocations. This behavior is consistent with the differing rheological properties of these phases. 
    more » « less