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   
                    
                            
                            The dependence of planetary tectonics on mantle thermal state: Applications to early Earth evolution
                        
                    
    
            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   
        
    
                            - Award ID(s):
- 1723057
- PAR ID:
- 10065516
- Date Published:
- Journal Name:
- Philosophical transactions - Royal Society. Mathematical, physical and engineering sciences
- ISSN:
- 1471-2962
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
- 
            
- 
            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
- 
            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
- 
            Earth's style of planetary heat transport is characterized by plate tectonics which requires rock strength to be reduced plastically in order to break an otherwise stagnant lithospheric lid, and for rocks to have a memory of past deformation to account for strain localization and the hysteresis implied by geological sutures. Here, we explore ∼107Rayleigh number, visco‐plastic, 3‐D global mantle convection with damage. We show that oceanic lithosphere‐only models generate strong toroidal‐poloidal power ratios and features such as a mix of long‐wavelength tectonic motions and smaller‐scale, back‐arc tectonics driven by downwellings. Undulating divergent plate boundaries can evolve to form overlapping spreading centers and microplates, promoted and perhaps stabilized by the effects of damage with long memory. The inclusion of continental rafts enhances heat flux variability and toroidal flow, including net rotation of the lithosphere, to a level seen in plate reconstructions for the Cenozoic. Both the super‐continental cycle and local rheological descriptions affect heat transport and tectonic deformation across a range of scales, and we showcase both general tectonic dynamics and regionally applied continental breakup scenarios. Our work points toward avenues for renewed analysis of the typical, mean behavior as well as the evolution of fluctuations in geological and model plate boundary evolution scenarios.more » « less
- 
            Subduction zones play a pivotal role in the mechanics of plate tectonics by providing the driving force through slab pull and weak megathrusts that facilitate the relative motion between tectonic plates. The initiation of subduction zones is intricately linked to the accumulation of slab pull and development of weakness at plate boundaries and, by consequence, the largest changes in the energetics of mantle convection. However, the transient nature of subduction initiation accompanied by intense subsequent tectonic activity, leaves critical evidence poorly preserved and making subduction initiation difficult to constrain. We overcome these limitations through a comprehensive analysis focused on Puysegur, a well-constrained extant example of subduction initiation offshore South Island, New Zealand. Through time-dependent, three-dimensional thermo-mechanical computations and quantitative comparison to new geophysical and geological observations, including topography, stratigraphy, and seismicity, we demonstrate that subduction initiation develops with a fast strain weakening described with a small characteristic displacement ( 4 to 8 km). Potential physical mechanisms contributing to the strain weakening are explored and we find that the observed fast weakening may arise through a combination of grain-size reduction within the lower lithosphere and fluid pressurization at shallower depths. With the shared commonality in the underlying physics of tectonic processes, the rapid strain weakening constrained at Puysegur offers insights into the formation of the first subduction during early Earth and the onset of plate tectonics.more » « less
 An official website of the United States government
An official website of the United States government 
				
			 
					 
					
 
                                    