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: On the Role of Rheological Memory for Convection‐Driven Plate Reorganizations
Abstract Understanding the temporal variability of plate tectonics is key to unraveling how mantle convection transports heat, and one critical factor for the formation and evolution of plate boundaries is rheological “memory,” that is, the persistence of weak zones. Here, we analyze the impact of such damage memory in global, oceanic‐lithosphere‐only models of visco‐plastic mantle convection. Self‐consistently‐formed weak zones are found to be reactivated in distinct ways, and convection preferentially selects such damaged zones for new plate boundaries. Reactivation of damage zones increases the frequency of plate reorganizations, and hence reduces the dominant periods of surface heat loss. The inheritance of distributed lithospheric damage thus dominates global surface dynamics over any local boundary stabilizing effects of weakening. In nature, progressive generation of weak zones may thus counteract and perhaps overcome any effects of reduced convective vigor throughout planetary cooling, with implications for the frequency of orogeny and convective transport throughout Wilson cycles.  more » « less
Award ID(s):
1853856 1927216
PAR ID:
10373410
Author(s) / Creator(s):
 ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
49
Issue:
18
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Abstract Mantle convection models based on geophysical constraints have provided us with a basic understanding of the forces driving and resisting plate motions on Earth. However, existing studies computing the balance of underlying forces are contradicting, and the impact of plate boundary geometry on surface deformation remains unknown. We address these issues by developing global instantaneous 3‐D mantle convection models with a heterogeneous density and viscosity distribution and weak plate boundaries prescribed using different geometries. We find that the plate boundary geometry of the Global Earthquake Model (GEM, Pagani et al., 2018,https://doi.org/10.1177/8755293020931866), featuring open plate boundaries with discrete lithospheric‐depth weak zones in the oceans and distributed crustal faults within continents, achieves the best fit to the observed GPS data with a directional correlation of 95.1% and a global point‐wise velocity residual of 1.87 cm/year. A good fit also requires plate boundaries being 3 to 4 orders of magnitude weaker than the surrounding lithosphere and low asthenospheric viscosities between 5 × 1017and 5 × 1018 Pa s. Models without asthenospheric and lower mantle heterogeneities retain on average 30% and 70% of the plate speeds, respectively. Our results show that Earth's plate boundaries are not uniform and better described by more discrete plate boundaries within the oceans and distributed faults within continents. Furthermore, they emphasize the impact of plate boundary geometry on the direction and speed of plate motions and reaffirm the importance of slab pull in the uppermost mantle as a major plate driving force. 
    more » « less
  3. Abstract Mantle plumes are typically considered secondary features of mantle convection, yet their surface effects over Earth's evolution may have been significant. We use 2‐D convection models to show that mantle plumes can in fact cause the termination of a subduction zone. This extreme case of plume‐slab interaction is found when the slab is readily weakened, for example, by damage‐type rheology, and the subducting slab is young. We posit that this mechanism may be relevant, particularly for the early Earth, and a subdued version of these plume‐slab interactions may remain relevant for modern subduction zones. Such core‐mantle boundary–surface interactions may be behind some of the complexity of tomographically imaged mantle structures, for example, in South America. More generally, plume “talk back” to subduction zones may make plate tectonics more episodic. 
    more » « less
  4. 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
  5. Abstract The Earth's long‐ and intermediate‐wavelength geoid anomalies are surface expressions of mantle convection and are sensitive to mantle viscosity. While previous studies of the geoid provide important constraints on the mantle radial viscosity variations, the mantle buoyancy in these studies, as derived from either seismic tomography or slab density models, may suffer significant uncertainties. In this study, we formulate 3‐D spherical mantle convection models with plate motion history since the Cretaceous that generate dynamically self‐consistent mantle thermal and buoyancy structures, and for the first time, use the dynamically generated slab structures and the observed geoid to place important constraints on the mantle viscosity. We found that non‐uniform weak plate margins and strong plate interiors are critical in reproducing the observed geoid and surface plate motion, especially the net lithosphere rotation (i.e., degree‐1 toroidal plate motion). In the best‐fit model, which leads to correlation of 0.61 between the modeled and observed geoid at degrees 4–12, the lower mantle viscosity is ∼1.3–2.5 × 1022 Pa⋅s and is ∼30 and ∼600–1,000 times higher than that in the transition zone and asthenosphere, respectively. Slab structures and the geoid are also strongly affected by slab strength, and the observations prefer moderately strong slabs that are ∼10–100 times stronger than the ambient mantle. Finally, a thin weak layer below the 670‐km phase change on a regional scale only in subduction zones produces stagnant slabs in the mantle transition zone as effectively as a weak layer on a global scale. 
    more » « less