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.