Journal cover Journal topic
Solid Earth An interactive open-access journal of the European Geosciences Union
Journal topic

Journal metrics

IF value: 2.921
IF 5-year value: 3.087
IF 5-year
CiteScore value: 4.8
SNIP value: 1.314
IPP value: 2.87
SJR value: 0.993
Scimago H <br class='widget-line-break'>index value: 38
Scimago H
h5-index value: 36
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.

  26 Oct 2020

26 Oct 2020

Review status
This preprint is currently under review for the journal SE.

Analytical solution for residual stress and strain preserved in anisotropic inclusion entrapped in isotropic host

Xin Zhong1,2, Marcin Dabrowski1,3, and Bjørn Jamtveit2 Xin Zhong et al.
  • 1Institut für Geologische Wissenschaften, Freie Universität Berlin, Malteserstrasse 74–100, 12449 Berlin, Germany
  • 2Physics of Geological Processes, The Njord Center, University of Oslo, Norway
  • 3Computational Geology Laboratory, Polish Geological Institute - NRI, Wrocław, Poland

Abstract. Raman elastic thermobarometry has recently been applied in many petrological studies to recover the pressure-temperature (P-T) conditions of mineral inclusion entrapment. Existing modelling methods in petrology either adopt an assumption of a spherical, isotropic inclusion embedded in an isotropic, infinite host, or use numerical techniques such as finite element method to simulate the residual stress and strain state preserved in the non-spherical anisotropic inclusion. Here, we use the Eshelby solution to develop an analytical framework for calculating the residual stress and strain state of an elastically anisotropic, ellipsoidal inclusion in an infinite, isotropic host. The analytical solution is applicable to any class of inclusion symmetry and an arbitrary inclusion aspect ratio. Explicit expressions are derived for some symmetry classes including e.g. tetragonal, hexagonal and trigonal.

The effect of changing the aspect ratio on residual stress is investigated including quartz, zircon, rutile, apatite and diamond inclusions in garnet host. Quartz is demonstrated to be the least affected, while rutile is the most affected. For prolate quartz inclusion (c-axis longer than a-axis), the effect of varying the aspect ratio on Raman shift is demonstrated to be insignificant. When c/a = 5, only ca. 0.3 cm−1 wavenumber variation is induced as compared to the spherical inclusion shape. For oblate quartz inclusions, the effect is more significant, when c/a = 0.5 ca. 0.8 cm−1 wavenumber variation for the 464 cm−1 band is induced compared to the reference spherical inclusion case. We also show that it is possible to fit an effective ellipsoid to obtain a proxy for the averaged residual stress/strain within faceted inclusion. The difference between the volumetrically averaged stress of a faceted inclusion and the analytically calculated stress from the best-fitted effective ellipsoid is calculated to obtain the root mean square deviation (RMSD) for quartz, zircon, rutile, apatite and diamond inclusions in garnet host. Based on the results of 500 randomly generated (a wide range of aspect ratio and random crystallographic orientation) faceted inclusion, we show that the volumetrically averaged stress serves as an excellent stress measure and the associated RMSD is less than 2 %, except for diamond with a systematically higher RMSD (ca. 8 %). This expands the applicability of the analytical solution for any arbitrary inclusion shape in practical Raman measurements.

Xin Zhong et al.

Interactive discussion

Status: open (until 07 Dec 2020)
Status: open (until 07 Dec 2020)
AC: Author comment | RC: Referee comment | SC: Short comment | EC: Editor comment
[Subscribe to comment alert] Printer-friendly Version - Printer-friendly version Supplement - Supplement

Xin Zhong et al.


Total article views: 315 (including HTML, PDF, and XML)
HTML PDF XML Total Supplement BibTeX EndNote
280 32 3 315 80 0 1
  • HTML: 280
  • PDF: 32
  • XML: 3
  • Total: 315
  • Supplement: 80
  • BibTeX: 0
  • EndNote: 1
Views and downloads (calculated since 26 Oct 2020)
Cumulative views and downloads (calculated since 26 Oct 2020)

Viewed (geographical distribution)

Total article views: 281 (including HTML, PDF, and XML) Thereof 278 with geography defined and 3 with unknown origin.
Country # Views %
  • 1



No saved metrics found.


No discussed metrics found.
Latest update: 04 Dec 2020
Publications Copernicus
Short summary
Elastic thermobarometry is an useful tool to recover the paleo-pressure and temperature. Here, we provide an analytical model based on the Eshelby solution to calculate the residual stress/strain preserved in a mineral inclusion exhumed from depth. The method applies to ellipsoidal, anisotropic inclusion in infinite isotropic host. Finite element method is also used for facet effect. It is found that the volumetrically avereged stress is a good proxy for the overall heterogeneous stress stage.
Elastic thermobarometry is an useful tool to recover the paleo-pressure and temperature. Here,...