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.


Title: An entropy method for geodynamic modelling of phase transitions: capturing sharp and broad transitions in a multiphase assemblage
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
Award ID(s):
1925677 2054605 1853388
PAR ID:
10370703
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Geophysical Journal International
Volume:
231
Issue:
3
ISSN:
0956-540X
Page Range / eLocation ID:
p. 1833-1849
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. The formation and preservation of compositional heterogeneities inside the Earth affect mantle convection patterns globally and control the long-term evolution of geochemical reservoirs. However, the distribution, nature, and size of reservoirs in the Earth’s mantle are poorly constrained. Here, we invert measurements of travel times and amplitudes of seismic waves interacting with mineralogical phase transitions at 400–700-km depth to obtain global probabilistic maps of temperature and bulk composition. We find large basalt-rich pools (up to 60% basalt fraction) surrounding the Pacific Ocean, which we relate to the segregation of oceanic crust from slabs that have been subducted since the Mesozoic. Segregation of oceanic crust from initially cold and stiff slabs may be facilitated by the presence of a weak hydrated layer in the slab or by weakening upon mineralogical transition due to grain-size reduction. 
    more » « less
  3. For plate tectonics to operate on a planet, mantle convective forces must be capable of forming weak, localized shear zones in the lithosphere that act as plate boundaries. Otherwise, a planet's mantle will convect in a stagnant lid regime, where subduction and plate motions are absent. Thus, when and how plate tectonics initiated on Earth is intrinsically tied to the ability of mantle convection to form plate boundaries; however, the physics behind this process are still uncertain. Most mantle convection models have employed a simple pseudoplastic model of the lithosphere, where the lithosphere "fails" and develops a mobile lid when stresses in the lithosphere reach the prescribed yield stress. With pseudoplasticity high mantle temperatures and high rates of internal heating, conditions relevant for the early Earth, impede plate boundary formation by decreasing lithospheric stresses, and hence favor a stagnant lid for the early Earth. However, when a model for shear zone formation based on grain size reduction is used, early Earth thermal conditions do not favor a stagnant lid. While lithosphere stress drops with increasing mantle temperature or heat production rate, the deformational work, which drives grain size reduction, increases. Thus the ability of convection to form weak plate boundaries is not impeded by early Earth thermal conditions. However, mantle thermal state does change the style of subduction and lithosphere mobility; high mantle temperatures lead to a more sluggish, drip-like style of subduction. This "sluggish lid" convection may be able to explain many of the key observations of early Earth crust formation processes preserved in the geologic record. Moreover, this work highlights the importance of understanding the microphysics of plate boundary formation for assessing early Earth tectonics, as different plate boundary formation mechanisms are influenced by mantle thermal state in fundamentally different ways. 
    more » « less
  4. Abstract Mantle convection drives changes in Earth's ellipsoidal figure and corresponding moment of inertia, causing shifts in the planet's rotation axis known as true polar wander (TPW). Using seismic tomography and the Back-and-Forth Nudging (BFN) method, we developed a time-dependent convection model that reconstructs the evolution of mantle density distribution and Earth’s moment of inertia over the past 70 million years. This modelling approach provides a close match with independent paleomagnetic constraints on the Cenozoic shifts of Earth’s rotation pole, specifically resolving the previously unexplained U-turn in TPW trajectory at approximately 50 million years ago. Our findings reveal TPW shifts exceeding 5 degrees, despite the temporal stability imposed by high lower-mantle viscosity and the stabilizing effect of Earth's remnant rotational bulge. Verification of our predicted changes in Earth’s ellipsoidal figure through independent paleomagnetic data provides a robust calibration for new predictions of convection-induced dynamic flattening variations. Over the past 70 million years, we find convection-induced changes of flattening that shift from -0.2% to +0.1 % during the Paleogene. Our predictions of Earth's axial precession frequency in the Paleogene align with recent independent cyclostratigraphic studies, thus validating our model's accuracy and supporting the hypothesis of reduced luni-solar tidal dissipation during this period. 
    more » « less
  5. Abstract The existence of a thin, weak asthenospheric layer beneath Earth’s lithospheric plates is consistent with existing geological and geophysical constraints, including Pleistocene glacio-isostatic adjustment, modeling of gravity anomalies, studies of seismic anisotropy, and post-seismic rebound. Mantle convection models suggest that a pronounced weak zone beneath the upper thermal boundary layer (lithosphere) may be essential to the plate tectonic style of convection found on Earth. The asthenosphere is likely related to partial melting and the presence of water in the sub-lithospheric mantle, further implying that the long-term evolution of the Earth may be controlled by thermal regulation and volatile recycling that maintain a geotherm that approaches the wet mantle solidus at asthenospheric depths. 
    more » « less