Mineral phase transitions can either hinder or accelerate mantle flow. In the present day, the formation of the bridgmanite + ferropericlase assemblage from ringwoodite at 660 km depth has been found to cause weak and intermittent layering of mantle convection. However, for the higher temperatures in Earth’s past, different phase transitions could have controlled mantle dynamics. We investigate the potential changes in convection style during Earth’s secular cooling using a new numerical technique that reformulates the energy conservation equation in terms of specific entropy instead of temperature. This approach enables us to accurately include the latent heat effect of phase transitions for mantle temperatures different from the average geotherm, and therefore fully incorporate the thermodynamic effects of realistic phase transitions in global-scale mantle convection modeling. We set up 2-D models with the geodynamics software ASPECT, using thermodynamic properties computed by HeFESTo, while applying a viscosity profile constrained by the geoid and mineral physics data and a visco-plastic rheology to reproduce self-consistent plate tectonics and Earth-like subduction morphologies. Our model results reveal the layering of plumes induced by the wadsleyite to garnet (majorite) + ferropericlase endothermic transition (between 420–600 km depth and over the 2000–2500 K temperature range). They show that this phase transition causes a large-scale and long-lasting temperature elevation in a depth range of 500–650 km depth if the potential temperature is higher than 1800 K, indicating that mantle convection may have been partially layered in Earth’s early history.

Stochastic tomography, made possible by dense deployments of seismic sensors, is used to identify previously undetected and poorly detected changes in the composition and mineral structure of Earth’s mantle. This technique inverts the spatial coherence of amplitudes and travel times of body waves to determine the depth and lateral dependence of the spatial spectrum of seismic velocity. The inverted spectrum is interpreted using predictions from the thermodynamic stability of different compositions and mineral phases as a function of temperature and pressure, in which the metamorphic temperature derivative of seismic velocities can be used as a proxy for the effects of heterogeneity induced in a region undergoing a phase change. Peaks in the metamorphic derivative of seismic velocity are found to closely match those found from applying stochastic tomography to elements of Earthscope and FLEX arrays. Within ± 12 km, peaks in the fluctuation of P velocity at 425, 500, and 600 km depth beneath the western US agree those predicted by a mechanical mixture of harzburgite and basalt in a cooler mantle transition zone. A smaller peak at 250 km depth may be associated with chemical heterogeneity induced by dehydration of subducted oceanic sediments, and a peak at 775 km depth with a phase change in subducted basalt. Non-detection of a predicted endothermic phase change near 660 km is consistent with its width being much less than 10 km. These interpretations of the heterogeneity spectrum are consistent with the known history of plate subduction beneath North America.