skip to main content


Title: Natural and Experimental Constraints on a Flow Law for Dislocation‐Dominated Creep in Wet Quartz
Abstract

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: = Aσne−(Q+VP)/RT. From this initial analysis, we identify differenteffectivestress 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 MPanr s−1;n = 3.5 ± 0.2;r = 0.49 ± 0.13;Q = 118 ± 5 kJ mol−1; andV = 2.59 ± 2.45 cm3 mol−1. The “high‐pressure” data set produces a different set of parameters: log(A) = −7.90 ± 0.34 MPanr s−1;n = 2.0 ± 0.1;r = 0.49 ± 0.13;Q = 77 ± 8 kJ mol−1; andV = 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 innat high confining pressure may result from an increase in the activity of grain size sensitive creep.

 
more » « less
Award ID(s):
1650173
NSF-PAR ID:
10361488
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
126
Issue:
5
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The effect of pressure on grain‐growth kinetics of olivine was investigated up to 10 GPa at 1773 K under relatively water‐poor conditions. The results are interpreted using a relationto obtain the activation volume = 5.0 ± 1.1 cm3/mol forn = 2 or = 5.2 ± 1.1 cm3/mol forn = 3. The small activation volume means that grain‐growth kinetics in pure olivine aggregates is fast even in the dry deep upper mantle, implying that grain‐size is controlled by the pinning by other phases or by dynamic recrystallization except for the early stage after the phase transformation from wadsleyite in upwelling materials. The present results are applied to seismic wave attenuation that is likely controlled by grain‐boundary processes. The inferred peak in attenuation just below the oceanic lithosphere‐asthenosphere boundary from the NoMelt array is difficult to be explained by the pressure effects assuming the absorption band behavior because such a model requires a much larger activation volume than determined in this work and it also fails to explain high attenuation in the deep asthenosphere. This suggests that either melt accumulation or other processes such as elastically accommodated grain‐boundary sliding (EAGBS) is responsible for the peak in attenuation. The present results are also applied to EAGBS. We suggest that the deep upper mantle is likely to be relaxed by EAGBS, which implies that the shear velocity of the deep upper mantle can be several percent smaller than that inferred from single crystal elasticity.

     
    more » « less
  2. Abstract

    Mineral inclusions are ubiquitous in metamorphic rocks and elastic models for host‐inclusion pairs have become frequently used tools for investigating pressure–temperature (P–T) conditions of mineral entrapment. Inclusions can retain remnant pressures () that are relatable to their entrapmentP–Tconditions using an isotropic elastic model andP–T–Vequations of state for host and inclusion minerals. Elastic models are used to constrainP–Tcurves, known as isomekes, which represent the possible inclusion entrapment conditions. However, isomekes require a temperature estimate for use as a thermobarometer. Previous studies obtained temperature estimates from thermometric methods external of the host‐inclusion system. In this study, we present the firstP–Testimates of quartz inclusion entrapment by integrating the quartz‐in‐garnet elastic model with titanium concentration measurements of inclusions and a Ti‐in‐quartz solubility model (QuiG‐TiQ). QuiG‐TiQ was used to determine entrapmentP–Tconditions of quartz inclusions in garnet from a quartzofeldspathic gneiss from Goodenough Island, part of the (ultra)high‐pressure terrane of Papua New Guinea. Raman spectroscopic measurements of the 128, 206, and 464 cm−1bands of quartz were used to calculate inclusion pressures using hydrostatic pressure calibrations (), a volume strain calculation (), and elastic tensor calculation (), that account for deviatoric stress.values calculated from the 128, 206, and 464 cm−1bands’ hydrostatic calibrations are significantly different from one another with values of 1.8 ± 0.1, 2.0 ± 0.1, and 2.5 ± 0.1 kbar, respectively. We quantified elastic anisotropy using the 128, 206 and 464 cm−1Raman band frequencies of quartz inclusions and stRAinMAN software (Angel, Murri, Mihailova, & Alvaro, 2019, 234:129–140). The amount of elastic anisotropy in quartz inclusions varied by ~230%. A subset of inclusions with nearly isotropic strains gives an averageandof 2.5 ± 0.2 and 2.6 ± 0.2 kbar, respectively. Depending on the sign and magnitude, inclusions with large anisotropic strains respectively overestimate or underestimate inclusion pressures and are significantly different (<3.8 kbar) from the inclusions that have nearly isotropic strains. Titanium concentrations were measured in quartz inclusions exposed at the surface of the garnet. The average Ti‐in‐quartz isopleth (19 ± 1 ppm [2σ]) intersects the average QuiG isomeke at 10.2 ± 0.3 kbar and 601 ± 6°C, which are interpreted as theP–Tconditions of quartzofeldspathic gneiss garnet growth and entrapment of quartz inclusions. TheP–Tintersection point of QuiG and Ti‐in‐quartz univariant curves represents mechanical and chemical equilibrium during crystallization of garnet, quartz, and rutile. These three minerals are common in many bulk rock compositions that crystallize over a wide range ofP–Tconditions thus permitting application of QuiG‐TiQ to many metamorphic rocks.

     
    more » « less
  3. Abstract

    We performed deformation and grain growth experiments on natural olivine aggregates with olivine water contents (COH = 600 ± 300 H/106 Si) similar to upper mantle olivine, at 1000–1200°C and 1,400 ± 100 MPa confining pressure. Our experiments differ from published grain growth studies in that most were (1) conducted on natural olivine cores rather than hot‐pressed aggregates and (2) dynamically recrystallized prior to or during grain growth. We combine our results with similar experiments performed at 1200–1300°C and fit the data to a grain growth relationship, yielding a growth exponent (p) of 3.2, activation energy (EG) 620 ± 145 kJ mol1(570 ± 145 kJ mol1when accounting for the role of temperature on water content), activation volume (VG) ~5 × 10−6 m3mol1, and rate constant (k0) 1.8 × 103 mp s−1. OurEGis within uncertainty of that predicted for dislocation creep of wet olivine (E* = 480 ± 40 kJ mol−1). Grain size in strain rate‐stepping samples adjusted to the olivine piezometer within 1.3–7.9% strain. The active grain boundary migration processes during deformation and dynamic recrystallization affect the kinetics of postdeformation grain growth, as grain boundary migration driven by strain energy density (ρGBM) may delay the onset of grain growth driven by interfacial energy (γGBM). We compared our postdeformation grain growth rates with data from previously published hydrostatic annealing experiments on synthetic olivine. At geologic timescales, the growth rates are much slower than predicted by the existing wet olivine grain growth law.

     
    more » « less
  4. Abstract

    Small repeating earthquakes are thought to represent rupture of isolated asperities loaded by surrounding creep. The observed scaling between recurrence interval and seismic moment,TrM1/6, contrasts with expectation assuming constant stress drop and no aseismic slip (TrM1/3). Here we demonstrate that simple crack models of velocity‐weakening asperities in a velocity‐strengthening fault predict theM1/6scaling; however, the mechanism depends on asperity radius,R. For small asperities (, whereis the nucleation radius) numerical simulations with rate‐state friction show interseismic creep penetrating inward from the edge, and earthquakes nucleate in the center and rupture the entire asperity. Creep penetration accounts for ∼25% of the slip budget, the nucleation phase takes up a larger fraction of slip. Stress drop increases with increasingR; the lack of self‐similarity being due to the finite nucleation dimension. Forsimulations exhibit simple cycles with ruptures nucleating from the edge. Asperities withexhibit complex cycles of partial and full ruptures. HereTris explained by an energy criterion: full rupture requires that the energy release rate everywhere on the asperity at least equals the fracture energy, leading to the scalingTrM1/6. Remarkably, in spite of the variability in behavior with source dimension, the scaling ofTrwith stress drop Δτ, nucleation length and creep ratevplis the same across all regimes:. This supports the use of repeating earthquakes as creepmeters and provides a physical interpretation for the scaling observed in nature.

     
    more » « less
  5. Summary

    Chemotrophic microorganisms gain energy for cellular functions by catalyzing oxidation–reduction (redox) reactions that are out of equilibrium. Calculations of the Gibbs energy (ΔGr) can identify whether a reaction is thermodynamically favourable and quantify the accompanying energy yield at the temperature, pressure and chemical composition in the system of interest. Based on carefully calculated values ofΔGr, we predict a novel microbial metabolism – sulfur comproportionation (3H2S ++ 2H+4S0+ 4H2O). We show that at elevated concentrations of sulfide and sulfate in acidic environments over a broad temperature range, this putative metabolism can be exergonic (ΔGr<0), yielding ~30–50 kJ mol−1. We suggest that this may be sufficient energy to support a chemolithotrophic metabolism currently missing from the literature. Other versions of this metabolism, comproportionation to thiosulfate (H2S ++ H2O) and to sulfite (H2S + 34+ 2H+), are only moderately exergonic or endergonic even at ideal geochemical conditions. Natural and impacted environments, including sulfidic karst systems, shallow‐sea hydrothermal vents, sites of acid mine drainage, and acid–sulfate crater lakes, may be ideal hunting grounds for finding microbial sulfur comproportionators.

     
    more » « less