You are here: Home / November 2025
206.189.94.68

November 2025

Phase Transitions Control Plume Layering During Earth’s Secular Cooling

Earth’s mantle convects, cooling the planet and shaping both its internal structure and surface tectonic regimes. While models and observations provide reasonable constraints on present-day mantle flow, it remains an open question how patterns of mantle convection have evolved throughout Earth’s history. A key factor is temperature. Petrological studies (e.g., Herzberg et al., 2010) estimate mantle cooling rates of up to ~100 K/Gyr, suggesting that the Archean mantle could have been about 300 K hotter than today. A hotter mantle would likely have been less viscous and hosted different mineral phase assemblages and transitions (Figure 1, comparing the different lines). Because endothermic phase transitions can hinder mantle flow and promote layering, it is important to understand how variations in these phase transitions affect mantle dynamics, including the degree of convective layering, the behavior of plumes and subducted slabs, and the material exchange between the upper and lower mantle.

In our recent study (Li et al., 2025), we investigated how varying phase transitions affect convection styles using numerical models. We present a series of 2D cylindrical annulus simulations covering a wide range of core–mantle boundary temperatures and initial mantle adiabats, representing Earth at different stages of cooling. The models evolve forward in time and include plastic deformation to generate plate-like subduction. We use the community code ASPECT (Heister et al., 2017) with the entropy method (Dannberg et al., 2022), which captures the full dynamic effects of phase transitions. To incorporate realistic material behavior, we use a look-up table generated by the thermodynamic software HeFESTo (Stixrude and Lithgow-Bertelloni, 2011), which includes complex phase relations through their variations in physical properties. By solving the energy equation for entropy instead of temperature, we fully include latent heat and avoid numerical instabilities, allowing us to model realistically sharp transitions and their impact on mantle dynamics.

Our models show that when the mantle potential temperature exceeds 1800 K (following the orange or yellow adiabat in Figure 1), the endothermic wadsleyite ⇔ garnet (majorite) + ferropericlase transition (boxed by a dot-dashed line in the right panel) occurs and significantly affects plume behavior. Unlike the widely-studied ringwoodite ⇔ bridgmanite + ferropericlase transition at 660 km, this wadsleyite transition (occurring between 450 and 590 km depth and over the 2000–2500 K temperature range) has received little attention before this work. In our hotter-mantle models (bottom panel in Figure 2), plumes are often trapped by this phase transition near transition-zone depths (500–650 km), forming a global hot layer and producing a variety of stalled plume shapes. This layering effect can reduce vertical plume transport between the lower and upper mantle by nearly 100%. Although slabs continue to sink through the transition and thus convection is not fully layered, the strong impedance of plumes may still influence Earth’s cooling history and the return of deep material to the surface through mantle upwellings. In addition, numerous secondary plumes rise from this hot layer, potentially affecting surface volcanism and plume-driven tectonic motion.

Summarizing the story (Figure 3), the wadsleyite ⇔ garnet (majorite) + ferropericlase transition may have dominated mantle dynamics in early Earth, causing strong plume layering until the mantle cooled below ~1750 K. This was likely followed by a transitional regime in which the effect weakened gradually and the hot layer diffused away. Although today’s mantle is too cool for this transition to occur, it may have played an important role in the past, and could still be dynamically significant on hotter terrestrial planets.

Contributed by: Ranpeng Li, Juliane Dannberg, Rene Gassmöller (GEOMAR Helmholtz Centre for Ocean Research Kiel); Carolina Lithgow-Bertelloni, Lars Stixrude (Earth, Planetary, and Space Sciences, University of California, Los Angeles)

Acknowledgements

This work was supported by NSF award EAR-1925677 and by Helmholtz funding EBP-01-08. University of Florida Research Computing has provided computational resources and support that have contributed to the research results reported in this publication.

References

Dannberg, J., Gassmöller, R., Li, R., Lithgow-Bertelloni, C., & Stixrude, L. (2022). An entropy method for geodynamic modelling of phase transitions: Capturing sharp and broad transitions in a multiphase assemblage. Geophysical Journal International, 231(3), 1833–1849. https://doi.org/10.1093/gji/ggac293

Heister, T., Dannberg, J., Gassmöller, R., & Bangerth, W. (2017). High accuracy mantle convection simulation through modern numerical methods. II: Realistic models and problems. Geophysical Journal International, 210(2), 833–851. https://doi.org/10.1093/gji/ggx195

Herzberg, C., Condie, K., & Korenaga, J. (2010). Thermal history of the Earth and its petrological expression. Earth and Planetary Science Letters, 292(1–2), 79–88. https://doi.org/10.1016/j.epsl.2010.01.022

Li, R., Dannberg, J., Gassmöller, R., Lithgow-Bertelloni, C., & Stixrude, L. (2025). How phase transitions impact changes in mantle convection style throughout Earth's history: From stalled plumes to surface dynamics. Geochemistry, Geophysics, Geosystems, 26, e2024GC011600. https://doi.org/10.1029/2024GC011600

Stixrude, L., & Lithgow-Bertelloni, C. (2011). Thermodynamics of mantle minerals-II. Phase equilibria. Geophysical Journal International, 184(3), 1180–1213. https://doi.org/10.1111/j.1365-246x.2010.04890.x

Xu, W., Lithgow-Bertelloni, C., Stixrude, L., & Ritsema, J. (2008). The effect of bulk composition and temperature on mantle seismic structure. Earth and Planetary Science Letters, 275(1–2), 70–79. https://doi.org/10.1016/j.epsl.2008.08.012

Figure 1. Material properties of pyrolite, a modeled bulk composition of Earth’s mantle, approximately 18% basalt and 82% harzburgite (Xu et al., 2008). Left: thermal expansivity plotted over pressure-temperature space; spikes indicate phase transitions. Some sharp transitions, like ringwoodite to bridgmanite + ferropericlase, appear only as small dots due to limited resolution. Right: density differences between two points at the same pressure but with a 200 K temperature difference, showing plume buoyancy relative to the mantle. Red areas mark endothermic phase transitions where plumes are denser and impeded. Solid colored lines show isentropes for various mantle potential temperatures. The inset highlights the wadsleyite → garnet (majorite) + ferropericlase transition region. This phase diagram is calculated by HeFESTo. For detailed mineral labeling, see Stixrude and Lithgow-Bertelloni (2011), Figure 16.

Figure 2. Snapshots of three models that have different dominated phase transitions, resulting in distinct convection styles. Yellow lines mark transitions from positive to negative effective thermal expansivity, showing the regions of phase transitions with negative Clapeyron slope. Click here to see videos of full model evolution:

  • Top panel: https://youtu.be/wD2iinksTj8
  • Middle panel: https://youtu.be/-hAFfHk7yD4
  • Bottom panel: https://youtu.be/R7cNCqxnjCc

Figure 3. A schematic summary. Time evolves in a clockwise direction. The dashed line marks a depth of 500 km, and the dotted shading indicates the pressure–temperature range where the wadsleyite ⇔ garnet (majorite) + ferropericlase phase transformation occurs. Purple shading marks areas of partial melting. The timeline calculation assumes a mantle potential temperature of 1900 K at ~3 Ga, cooling by about 100 K per Gyr.