skip to main content

Title: Global Analysis of Experimental Data on the Rheology of Olivine Aggregates

Following the reanalysis of individual experimental runs of some widely cited studies (Jain et al., 2018,, we revisit the global data analysis of Korenaga and Karato (2008, with a significantly improved version of their Markov chain Monte Carlo inversion. Their algorithm, previously corrected by Mullet et al. () to minimize potential parameter bias, is further modified here to estimate more efficiently interrun biases in global data sets. Using the refined Markov chain Monte Carlo inversion technique, we simultaneously analyze experimental data on the deformation of olivine aggregates compiled from different studies. Realistic composite rheological models, including both diffusion and dislocation creep, are adopted, and the role of dislocation‐accommodated grain boundary sliding is also investigated. Furthermore, the influence of interrun biases on inversion results is studied using experimental and synthetic data. Our analysis shows that existing data can tightly constrain the grain‐size exponent for diffusion creep at ∼2, which is different from the value commonly assumed (p= 3). Different data sets and model assumptions, however, yield nonoverlapping estimates on other flow‐law parameters, and the flow‐law parameters for grain boundary sliding are poorly resolved in most cases. We thus provide a few plausible candidate flow‐law models for olivine rheology to facilitate future geodynamic modeling. The availability of more data that explore a wider range of experimental conditions, especially higher pressures, is essential to improve our understanding of upper mantle rheology.

more » « less
Award ID(s):
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Page Range / eLocation ID:
p. 310-334
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The rheology of the upper mantle impacts a variety of geodynamic processes, including postseismic deformation following great earthquakes and post-glacial rebound. The deformation of upper mantle rocks is controlled by the rheology of olivine, the most abundant upper mantle mineral. The mechanical properties of olivine at steady state are well constrained. However, the physical mechanism underlying transient creep, an evolutionary, hardening phase converging to steady state asymptotically, is still poorly understood. Here, we constrain a constitutive framework that captures transient creep and steady state creep consistently using the mechanical data from laboratory experiments on natural dunites containing at least 94% olivine under both hydrous and anhydrous conditions. The constitutive framework represents a Burgers assembly with a thermally activated nonlinear stress-versus-strain-rate relationship for the dashpots. Work hardening is obtained by the evolution of a state variable that represents internal stress. We determine the flow law parameters for dunites using a Markov chain Monte Carlo method. We find the activation energy $$430\pm 20$$ 430 ± 20   and $$250\pm 10$$ 250 ± 10  kJ/mol for dry and wet conditions, respectively, and the stress exponent $$2.0\pm 0.1$$ 2.0 ± 0.1 for both the dry and wet cases for transient creep, consistently lower than those of steady-state creep, suggesting a separate physical mechanism. For wet dunites in the grain-boundary sliding regime, the grain-size dependence is similar for transient creep and steady-state creep. The lower activation energy of transient creep could be due to a higher jog density of the corresponding soft-slip system. More experimental data are required to estimate the activation volume and water content exponent of transient creep. The constitutive relation used and its associated flow law parameters provide useful constraints for geodynamics applications. Graphical Abstract 
    more » « less
  2. Abstract

    Relationships between the recrystallized grain size and stress are investigated for experimentally deformed water‐added quartz aggregates. For stresses ≥100 MPa there is a variation in the measured recrystallized grain size for a given stress. This variation correlates with a change in thec‐axis fabric in general shear experiments, where samples with larger recrystallized grain sizes for a given stress have dominantly prism c‐axis fabrics and samples with smaller recrystallized grain sizes for a given stress have dominantly basal c‐axis fabrics. The dislocation creep flow law also changes at conditions where these twoc‐axis fabrics form (Tokle et al., 2019, Using the wattmeter model (Austin & Evans, 2007,, different piezometric relationships are quantified for samples that develop prism and basal c‐axis fabrics, respectively. The wattmeter model is sensitive to grain growth kinetics; a new grain growth law for quartz is formulated based on reanalysis of microstructures in samples from previous work. The activation enthalpies and water fugacity exponents for our grain growth law and dislocation creep flow laws are the same within error, suggesting the recrystallized grain size versus stress relationships are nearly independent of temperature and water fugacity, consistent with laboratory observations. The wattmeters successfully predict the recrystallized grain size versus stress relationships of all quartzite samples from experiments with added water. These results support the use and extrapolation of the wattmeter model for both experimental and geologic conditions to investigate the stress state and grain size evolution of quartz rich rocks.

    more » « less
  3. Abstract

    Apparent polar wander paths (APWPs) synthesized from paleomagnetic poles provide the most direct data for reconstructing past paleogeography and plate motions for times earlier than ca. 200 Ma. In this contribution, we describe a new method for APWP synthesis that extends the paleomagnetic Euler pole analysis of Gordon et al. (1984, by placing it within the framework of a Bayesian inverse problem. This approach incorporates uncertainties in pole positions and age that are often ignored in standard treatments. The paleomagnetic Euler poles resulting from the inversions provide estimates for full‐vector plate motion (both latitude and longitude) and associated uncertainty. The method allows for inverting for one or more Euler poles with the timing of changepoints being solved as part of the inversion. In addition, the method allows the incorporation of true polar wander rotations, thus providing an avenue for probabilistic partitioning of plate tectonic motion and true polar wander based on paleomagnetic poles. We show example inversions on synthetic data to demonstrate the method's capabilities. We illustrate application of the method to Cenozoic Australia paleomagnetic poles which can be compared to independent plate reconstructions. A two‐Euler pole inversion for the Australian record recovers northward acceleration of Australia in the Eocene with rates that are consistent with plate reconstructions. We also apply the method to constrain rapid rates of motion for cratonic North America associated with the Keweenawan Track of late Mesoproterozoic paleomagnetic poles. The application of Markov chain Monte Carlo methods to estimate paleomagnetic Euler poles can open new directions in quantitative paleogeography.

    more » « less
  4. 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
  5. Abstract

    In the D‐region, the ionization rate cannot be detected directly with any known measurement technique, therefore it must be estimated. Starting from space‐based measurements of precipitating particle flux, we estimate the ionization rate in the atmosphere using the Electron Precipitation Monte Carlo transport method. This ionization rate is used to calculate the expected electron density in the D‐region with the Glukhov‐Pasko‐Inan five species (GPI5) atmospheric chemistry model. We then compare the simulated electron density with that measured by the Poker Flat Incoherent Scatter Radar (PFISR). From ground‐based radar measurements of electron density enhancements due to sub‐relativistic and relativistic electron precipitation, we present a method to extract the ionization rate altitude profiles using inverse theory. We use this estimation of ionization rate to find the energy distribution of the precipitating particles. With this inverse method, we are able to link ground measurements of electron density to the precipitating flux in a time dependent manner and with uncertainty in the inverted parameters. The method was tested on synthetic data and applied to specific PFISR data sets. The method is able to retrieve the ionization rate altitude profiles that, when forward modeled, return the expected electron densities within ∼7% error as compared to the PFISR data. For the case presented here, the arbitrary energy distribution inversion results are comparable in magnitude and shape to those presented in Turunen et al. (2016, for the inversion of a single event of pulsating aurora observed by EISCAT.

    more » « less