Quantifying Intrinsic and Extrinsic Contributions to Seismic Anisotropy
in Tomographic Models
Abstract
Seismic anisotropy in the Earth’s mantle inferred from seismic
observations is usually interpreted either in terms of intrinsic
anisotropy due to Crystallographic Preferred Orientation (CPO) of
minerals, or extrinsic anisotropy due to rock-scale Shape Preferred
Orientation (SPO). The coexistence of both contributions misconstrues
the origins of seismic anisotropy observed in seismic tomography models.
It is thus essential to discriminate CPO from SPO.
Homogenization/upscaling theory provides means to achieve this goal.
This theory enables to compute the effective elastic properties of a
heterogeneous medium, as seen by long-period waves. In this work, we
investigate the effects of upscaling an intrinsically anisotropic and
highly heterogeneous Earth’s mantle. We show analytically in 1-D that
the full effective radial anisotropy ξ * is approximately the product of
the effective intrinsic radial anisotropy ξ * CPO and the extrinsic
radial anisotropy ξ * SPO : ξ * ≈ ξ * CPO x ξ * SPO. This law is
verified numerically in the case of a 2-D marble cake model of the
mantle with a binary composition, and in the presence of CPO obtained
from a micro-mechanical model of olivine deformation. We compute the
long-wavelength effective equivalent of this mantle model using the 3-D
non-periodic elastic homogenization technique. Our numerical findings
predict that for wavelenghts smaller than the scale of deformation
patterns, tomography may overestimate the true anisotropy (i.e.
intrinsic anisotropy due to CPO) due to significant SPO-induced
extrinsic anisotropy. However, at wavelenghts larger than deformation
patterns, intrinsic anisotropy is always underestimated in tomographic
models due to the spatial averaging of the preferred orientation of
anisotropic minerals. Thus, we show that it is imperative to homogenize
a CPO evolution model first before drawing comparisons with tomographic
models. As a demonstration, we use our composite law with a homogenized
CPO model of a plate-driven flow underneath a mid-ocean ridge, to
estimate the SPO contibution to an existing tomographic model of radial
anisotropy.