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: Scaling laws for stagnant-lid convection with a buoyant crust
SUMMARY Stagnant-lid convection, where subduction and surface plate motion is absent, is common among the rocky planets and moons in our solar system, and likely among rocky exoplanets as well. How stagnant-lid planets thermally evolve is an important issue, dictating not just their interior evolution but also the evolution of their atmospheres via volcanic degassing. On stagnant-lid planets, the crust is not recycled by subduction and can potentially grow thick enough to significantly impact convection beneath the stagnant lid. We perform numerical models of stagnant-lid convection to determine new scaling laws for convective heat flux that specifically account for the presence of a buoyant crustal layer. We systematically vary the crustal layer thickness, crustal layer density, Rayleigh number and Frank–Kamenetskii parameter for viscosity to map out system behaviour and determine the new scaling laws. We find two end-member regimes of behaviour: a ‘thin crust limit’, where convection is largely unaffected by the presence of the crust, and the thickness of the lithosphere is approximately the same as it would be if the crust were absent; and a ‘thick crust limit’, where the crustal thickness itself determines the lithospheric thickness and heat flux. Scaling laws for both limits are developed and fit the numerical model results well. Applying these scaling laws to rocky stagnant-lid planets, we find that the crustal thickness needed for convection to enter the thick crust limit decreases with increasing mantle temperature and decreasing mantle reference viscosity. Moreover, if crustal thickness is limited by the formation of dense eclogite, and foundering of this dense lower crust, then smaller planets are more likely to enter the thick crust limit because their crusts can grow thicker before reaching the pressure where eclogite forms. When convection is in the thick crust limit, mantle heat flux is suppressed. As a result, mantle temperatures can be elevated by 100 s of degrees K for up to a few Gyr in comparison to a planet with a thin crust. Whether convection enters the thick crust limit during a planet’s thermal evolution also depends on the initial mantle temperature, so a thick, buoyant crust additionally acts to preserve the influence of initial conditions on stagnant-lid planets for far longer than previous thermal evolution models, which ignore the effects of a thick crust, have found.  more » « less
Award ID(s):
1723057
PAR ID:
10389156
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Geophysical Journal International
Volume:
228
Issue:
1
ISSN:
0956-540X
Page Range / eLocation ID:
631 to 663
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The trondhjemite-tonalite-granodiorite (TTG) suite of rocks prominent in Earth’s Archean continents is thought to form by melting of hydrated basalt, but the specific tectonic settings of formation are unclear. Models for TTG genesis range from melting of downgoing mafic crust during subduction into a hotter mantle to melting at the base of a thick crustal plateau; while neither uniquely defines a global tectonic regime, the former is consistent with mobile lid tectonics and the latter a stagnant lid. One major problem for a subduction model is slabs sinking too quickly and steeply in a hotter mantle to melt downgoing crust. I show, however, that grain size reduction in the lithosphere leads to relatively strong plate boundaries on the early Earth, which slow slab sinking. During this “sluggish subduction,” sinking plates can heat up enough to melt when the mantle temperature is ≳1600 °C. Crustal melting via sluggish subduction can thus explain TTG formation during the Archean due to elevated mantle temperatures and the paucity of TTG production since due to mantle cooling. 
    more » « less
  2. 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
  3. Abstract Terrestrial planets in the habitable zone (HZ) of Sun-like stars are priority targets for detection and observation by the next generation of space telescopes. Earth's long-term habitability may have been tied to the geological carbon cycle, a process critically facilitated by plate tectonics. In the modern Earth, plate motion corresponds to a mantle convection regime called mobile lid. The alternate, stagnant-lid regime is found on Mars and Venus, which may have lacked strong enough weathering feedback to sustain surface liquid water over geological timescales if initially present. Constraining observational strategies able to infer the most common regime in terrestrial exoplanets requires quantitative predictions of the atmospheric composition of planets in either regime. We use end-member models of volcanic outgassing and crust weathering for the stagnant- and mobile-lid convection regimes, which we couple to models of atmospheric chemistry and climate and ocean chemistry to simulate the atmospheric evolution of these worlds in the HZ. In our simulations under the two alternate regimes, we find that the fraction of planets possessing climates consistent with surface liquid water is virtually the same. Despite this unexpected similarity, we predict that a mission capable of detecting atmospheric CO2abundance above 0.1 bar in 25 terrestrial exoplanets is extremely likely (≥95% of samples) to infer the dominant interior convection regime in that sample with strong evidence (10:1 odds). This offers guidance for the specifications of the Habitable Worlds Observatory NASA concept mission and other future missions capable of probing samples of habitable exoplanets. 
    more » « less
  4. Abstract Seismic observations indicate accumulation of subducted slabs in the mantle transition zone in many subduction zones. By systematically conducting 2‐D numerical experiments, we demonstrate that a weak layer or zone beneath the spinel‐to‐post‐spinel phase transition leads to horizontally deflected (stagnant) slab structures in the mantle transition zone, which is consistent with recent studies of 3‐D global mantle convection models. Trench retreat velocity, Clapeyron slope and the viscosity contrast between the lower mantle and mantle transition zone also affect horizontally deflected slab formation. By considering grain size dependent viscosity and grainsize evolution for slabs going through the phase change in the lower mantle, our models with a dynamically generated weak zone beneath the phase boundary indicate that the geometry and viscosity reduction of the weak zone is strongly affected by grain growth rate. A smaller grain growth rate results in a thicker and wider weak zone that promotes deflected slab formation. 
    more » « less
  5. Abstract Despite significant study, when and how plate tectonics initiated on Earth remains contentious. Geologic evidence from some of Earth's earliest cratons has been interpreted as reflecting the formation of initial continental blocks by non‐subduction processes, which then trigger subduction initiation at their margins. Numerical models of mantle convection with a plastic yield stress rheology have shown this scenario is plausible. However, whether continents can trigger subduction initiation has not been tested with other rheologies. We, therefore, use numerical models of mantle convection with an imposed continental block to test whether continents facilitate subduction initiation with a grain‐damage mechanism, where weak shear zones form by grain size reduction. Our results show that continents modestly enhance stresses in the lithosphere, but not enough to significantly impact lithospheric damage or subduction initiation: continents have minimal influence on lithospheric damage or plate speed, nor does subduction preferentially initiate at the continental margin. A new regime diagram that includes continental blocks shows only a small shift in the boundary between the mobile‐lid and stagnant‐lid regimes when continents are added. However, as we do find that stresses are modestly enhanced at the continental margin in our models, we develop a scaling law for this stress enhancement to more fully test whether continents could trigger subduction initiation on early Earth. We find that lithospheric stresses supplied by continents are not sufficient to initiate subduction on the early Earth on their own with grain‐damage rheology; instead, additional factors would be required. 
    more » « less