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.

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 10:00 PM ET on Friday, February 6 until 10:00 AM ET on Saturday, February 7 due to maintenance. We apologize for the inconvenience.


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. 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 −1 to seismic 10 − 10 3  s −1 strain 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 Whole rock compositions at Buldir Volcano, western Aleutian arc, record a strong, continuous trend of iron depletion with decreasing MgO, classically interpreted as a calc-alkaline liquid line of descent. In contrast, olivine-hosted melt inclusions have higher total iron (FeO*) than whole rocks and show little change in FeO* with decreasing MgO. To investigate this discrepancy and determine the conditions required for strong iron depletion, we conducted oxygen fugacity (ƒO2) buffered, water-saturated crystallization experiments at 800 MPa and ƒO2 = QFM + 1.6 ± 0.4 (1$$\sigma$$) (where QFM refers to the quartz-fayalite-magnetite buffer) on a high-Al, basaltic starting material modeled after a Buldir lava. Experimental conditions were informed by olivine-hosted melt inclusions that record minimum entrapment pressures as high as 570 MPa, >6 wt % H2O, and ƒO2 of QFM + 1.4 (±0.2), making Buldir one of the most oxidized and wettest arc volcanoes documented globally. The experiments produce melts with Si-enrichment and Fe-depletion signatures characteristic of evolved, calc-alkaline magmas at the lowest MgO, although FeO* remains roughly constant over most of the experimental temperature range. Experiments saturate CrAl-spinel and olivine at 1160°C, followed by clinopyroxene and Al-spinel at 1085°C, hornblende at 1060°C, and, finally, plagioclase and magnetite between 1040°C and 960°C. Hornblende crystallization, not magnetite, generates the largest increase in SiO2 and largest decrease in FeO* in coexisting melts. Compositions of melt inclusions are consistent with experimental melts and reflect crystallization of a basaltic parent magma at high PH2O. In contrast, the whole rock compositional trends are influenced by magma mixing and phenocryst redistribution and accumulation. The crystallization experiments and natural liquids (melt inclusions and groundmass glass) from Buldir suggest that for an oxidized, hydrous primary basalt starting composition, significant Fe depletion from the melt will not occur until intermediate to late stages of magma crystallization (< ~4.5 wt % MgO). We conclude that the Buldir whole rock trend cannot be reproduced by crystallization at arc-relevant oxygen fugacities and is not a true liquid line of descent, warranting caution when interpreting volcanic trends globally. 
    more » « less