Probabilistic tomography maps chemical heterogeneities throughout the lower mantle.

We obtained likelihoods in the lower mantle for long-wavelength models of bulk sound and shear wave speed, density, and boundary topography, compatible with gravity constraints, from normal mode splitting functions and surface wave data. Taking into account the large uncertainties in Earth's thermodynamic reference state and the published range of mineral physics data, we converted the tomographic likelihoods into probability density functions for temperature, perovskite, and iron variations. Temperature and composition can be separated, showing that chemical variations contribute to the overall buoyancy and are dominant in the lower 1000 kilometers of the mantle.