skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on February 1, 2026

Title: How Phase Transitions Impact Changes in Mantle Convection Style Throughout Earth's History: From Stalled Plumes to Surface Dynamics
Abstract 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 softwareAspect, 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 plate‐like behavior and Earth‐like subduction morphologies. Our model results reveal the layering of plumes induced by the wadsleyite to garnet (majorite) + ferropericlase endothermic transition (between 450 and 590 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 of the mantle is higher than 1800 K, indicating that mantle convection may have been partially layered in Earth's early history.  more » « less
Award ID(s):
2223935 1925677
PAR ID:
10615197
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
American Geophysical Union
Date Published:
Journal Name:
Geochemistry, Geophysics, Geosystems
Volume:
26
Issue:
2
ISSN:
1525-2027
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. SUMMARY Phase transitions play an important role for the style of mantle convection. While observations and theory agree that a substantial fraction of subducted slabs and rising plumes can move through the whole mantle at present day conditions, this behaviour may have been different throughout Earth’s history. Higher temperatures, such as in the early Earth, cause different phase transitions to be dominant, and also reduce mantle viscosity, favouring a more layered style of convection induced by phase transitions. A period of layered mantle convection in Earth’s past would have significant implications for the secular evolution of the mantle temperature and the mixing of mantle heterogeneities. The transition from layered to whole mantle convection could lead to a period of mantle avalanches associated with a dramatic increase in magmatic activity. Consequently, it is important to accurately model the influence of phase transitions on mantle convection. However, existing numerical methods generally preclude modelling phase transitions that are only present in a particular range of pressures, temperatures or compositions, and they impose an artificial lower limit on the thickness of phase transitions. To overcome these limitations, we have developed a new numerical method that solves the energy equation for entropy instead of temperature. This technique allows for robust coupling between thermodynamic and geodynamic models and makes it possible to model realistically sharp phase transitions with a wide range of properties and dynamic effects on mantle processes. We demonstrate the utility of our method by applying it in regional and global convection models, investigating the effect of individual phase transitions in the Earth’s mantle with regard to their potential for layering flow. We find that the thickness of the phase transition has a bigger influence on the style of convection than previously thought: with all other parameters being the same, a thin phase transition can induce fully layered convection where a broad phase transition would lead to whole-mantle convection. Our application of the method to convection in the early Earth illustrates that endothermic phase transitions may have induced layering for higher mantle temperatures in the Earth’s past. 
    more » « less
  2. SUMMARY We derive exact expressions for the thermal expansivity, heat capacity and bulk modulus for assemblages with arbitrarily large numbers of components and phases, including the influence of phase transformations and chemical exchange. We illustrate results in simple two-component, two-phase systems, including Mg–Fe olivine-wadsleyite and Ca–Mg clinopyroxene-orthopyroxene and for a multicompontent model of mantle composition in the form of pyrolite. For the latter we show results for the thermal expansivity and heat capacity over the entire mantle pressure–temperature regime to 40 GPa, or a depth of 1000 km. From the thermal expansivity, we derive a new expression for the phase buoyancy parameter that is valid for arbitrarily large numbers of phases and components and which is defined at every point in pressure–temperature space. Results reveal regions of the mantle where the magnitude of the phase buoyancy parameter is larger in magnitude than for those phase transitions that are most commonly included in mantle convection simulations. These regions include the wadsleyite to garnet and ferropericlase transition, which is encountered along hot isentropes (e.g. 2000 K potential temperature) in the transition zone, and the ferropericlase and stishovite to bridgmanite transition, which is encountered along cold isentropes (e.g. 1000 K potential temperature) in the shallow lower mantle. We also show the bulk modulus along a typical mantle isentrope and relate it to the Bullen inhomogeneity parameter. All results are computed with our code HeFESTo, updates and improvements to which we discuss, including the implementation of the exact expressions for the thermal expansivity, heat capacity and bulk modulus, generalization to allow for pressure dependence of non-ideal solution parameters and an improved numerical scheme for minimizing the Gibbs free energy. Finally, we present the results of a new global inversion of parameters updated to incorporate more recent results from experiment and first principles theory, as well as a new phase (nal phase), and new species: Na-majorite and the NaAlO2 end-member of ferropericlase. 
    more » « less
  3. Abstract Phase transitions in the mantle control its internal dynamics and structure. The post-spinel transition marks the upper–lower mantle boundary, where ringwoodite dissociates into bridgmanite plus ferropericlase, and its Clapeyron slope regulates mantle flow across it. This interaction has previously been assumed to have no lateral spatial variations, based on the assumption of a linear post-spinel boundary in pressure and temperature. Here we present laser-heated diamond anvil cell experiments with synchrotron X-ray diffraction to better constrain this boundary, especially at higher temperatures. Combining our data with results from the literature, and using a global analysis based on machine learning, we find a pronounced nonlinearity in the post-spinel boundary, with its slope ranging from –4 MPa/K at 2100 K, to –2 MPa/K at 1950 K, and to 0 MPa/K at 1600 K. Changes in temperature over time and space can therefore cause the post-spinel transition to have variable effects on mantle convection and the movement of subducting slabs and upwelling plumes. 
    more » « less
  4. Abstract Understanding the mineralogy of the Earth's interior is a prerequisite for unravelling the evolution and dynamics of our planet. Here, we conducted high pressure-temperature experiments mimicking the conditions of the deep lower mantle (DLM, 1800–2890 km in depth) and observed surprising mineralogical transformations in the presence of water. Ferropericlase, (Mg, Fe)O, which is the most abundant oxide mineral in Earth, reacts with H2O to form a previously unknown (Mg, Fe)O2Hx (x≤1) phase. The (Mg, Fe)O2Hx has the pyrite structure and it coexists with the dominant silicate phases, bridgmanite and post-perovskite. Depending on Mg content and geotherm temperatures, the transformation may occur at 1800 km for (Mg0.6Fe0.4)O or beyond 2300 km for (Mg0.7Fe0.3)O. The (Mg, Fe)O2Hx is an oxygen excess phase that stores an excessive amount of oxygen beyond the charge balance of maximum cation valences (Mg2+, Fe3+, and H+). This important phase has a number of far-reaching implications including the extreme redox inhomogeneity, deep-oxygen reservoirs in the DLM, and an internal source for modulating oxygen in the atmosphere. 
    more » « less
  5. Sublithospheric diamonds and the inclusions they may carry crystallize in the asthenosphere, transition zone, or uppermost lower mantle (from 300 to ∼800 km), and are the deepest minerals so far recognized to form by plate tectonics. These diamonds are distinctive in their deformation features, low nitrogen content, and inclusions of these major mantle minerals: majoritic garnet, clinopyroxene, ringwoodite, CaSi perovskite, ferropericlase, and bridgmanite or their retrograde equivalents. The stable isotopic compositions of elements within these diamonds (δ11B, δ13C, δ15N) and their inclusions (δ18O, δ56Fe) are typically well outside normal mantle ranges, showing that these elements were either organic (C) or modified by seawater alteration (B, O, Fe) at relatively low temperatures. Metamorphic minerals in cold slabs are effective hosts that transport C as CO3 and H as H2O, OH, or CH4 below the island arc and mantle wedge. Warming of the slab generates carbonatitic melts, supercritical aqueous fluids, or metallic liquids, forming three types of sublithospheric diamonds. Diamond crystallization occurs by movement and reduction of mobile fluids as they pass through host mantle via fractures—a process that creates chemical heterogeneity and may promote deep focus earthquakes. Geobarometry of majoritic garnet inclusions and diamond ages suggest upward transport, perhaps to the base of mantle lithosphere. From there, diamonds are carried to Earth's surface by eruptions of kimberlite magma. Mineral assemblages in sublithospheric diamonds directly trace Earth's deep volatile cycle, demonstrating how the hydrosphere of a rocky planet can connect to its solid interior.▪Sublithospheric diamonds from the deep upper mantle, transition zone, and lower mantle host Earth's deepest obtainable mineral samples.▪Low-temperature seawater alteration of the ocean floor captures organic and inorganic carbon at the surface eventually to become some of the most precious gem diamonds.▪Subduction transports fluids in metamorphic minerals to great depth. Fluids released by slab heating migrate, react with host mantle to induce diamond crystallization, and may trigger earthquakes.▪Sublithospheric diamonds are powerful tracers of subduction—a plate tectonic process that deeply recycles part of Earth's planetary volatile budget. 
    more » « less