Petrophysical constraints on the seismic properties of the Kaapvaal craton mantle root
Abstract. We calculated the seismic properties of 47 mantle xenoliths from 9 kimberlitic pipes in the Kaapvaal craton based on their modal composition, the crystal preferred orientations (CPO) of olivine, ortho- and clinopyroxene, and garnet, the Fe content of olivine, and the pressures and temperatures at which the rocks were equilibrated. These data allow constraining the variation of seismic anisotropy and velocities with depth. The fastest P wave and fast split shear wave (S1) polarization direction is always close to olivine [100] maximum. Changes in olivine CPO symmetry result in minor variations in the seismic anisotropy patterns. Seismic anisotropy is higher for high olivine contents and stronger CPO. Maximum P waves azimuthal anisotropy (AVp) ranges between 2.5 and 10.2% and S waves polarization anisotropy (AVs) between 2.7 and 8%. Seismic properties averaged in 20 km thick intervals depth are, however, very homogeneous. Based on these data, we predict the anisotropy that would be measured by SKS, Rayleigh (SV) and Love (SH) waves for 5 end-member orientations of the foliation and lineation. Comparison to seismic anisotropy data in the Kaapvaal shows that the coherent fast directions, but low delay times imaged by SKS studies and the low azimuthal anisotropy and SH faster than SV measured using surface waves may only be consistently explained by dipping foliations and lineations. The strong compositional heterogeneity of the Kaapvaal peridotite xenoliths results in up to 3% variation in density and in up to 2.3% of variation Vp, Vs and the Vp/Vs ratio. Fe depletion by melt extraction increases Vp and Vs, but decreases the Vp/Vs ratio and density. Orthopyroxene enrichment decreases the density and Vp, but increases Vs, strongly reducing the Vp/Vs ratio. Garnet enrichment increases the density, and in a lesser manner Vp and the Vp/Vs ratio, but it has little to no effect on Vs. These compositionally-induced variations are slightly higher than the velocity perturbations imaged by body-wave tomography, but cannot explain the strong velocity anomalies reported by surface wave studies. Comparison of density and seismic velocity profiles calculated using the xenoliths' compositions and equilibrium conditions to seismological data in the Kaapvaal highlights that: (i) the thickness of the craton is underestimated in some seismic studies and reaches at least 180 km, (ii) the deep sheared peridotites represent very local modifications caused and oversampled by kimberlites, and (iii) seismological models probably underestimate the compositional heterogeneity in the Kaapvaal mantle root, which occurs at a scale much smaller than the one that may be sampled seismologically.