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: High‐Resolution Mantle Flow Models Reveal Importance of Plate Boundary Geometry and Slab Pull Forces on Generating Tectonic Plate Motions
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
Award ID(s):
2015848 2028346 1925575 1902308 1925677
PAR ID:
10441512
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
128
Issue:
8
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Subduction initiation remains one of the least understood global processes of plate tectonics. Prominent models have been cast in terms of two broad classes: spontaneous cases due to lithospheric gravitational instabilities and induced cases due to forced plate convergence. Yet gravitationally unstable lithosphere is old, strong, and difficult to begin to bend into a subduction zone and convergent forces necessary to begin subduction are often too large given the plates involved. These models also consider the asthenospheric mantle as passive, even though relative motion between slabs and the asthenosphere has long been regarded as a strong control on subduction dynamics. Here I propose that subduction-transform edge propagator (STEP) faults can initiate subduction depending on the absolute motion of the STEP fault with respect to the asthenosphere. STEP faults form where subduction zones end and the subducting plate tears forming a down flexed transcurrent plate boundary at the surface shearing against the adjacent rear arc lithospheric plate. However, STEP faults are not simple transcurrent boundaries. Absolute motion of the down flexed STEP fault edge with respect to the surrounding asthenosphere can produce a strong sea anchor force that either continues to bend the edge downward, initiating subduction, or opposes slab bending, inhibiting subduction. In the south Pacific, the southern end of the New Hebrides Trench and the northern end of the Tonga Trench are type-example STEP faults with opposite senses of dip but both moving northward with respect to the asthenosphere. The northward dipping New Hebrides STEP fault moves northward in a mantle reference frame creating a strong asthenospheric flow against the STEP fault edge, inducing active subduction at the Matthew-Hunter trench. In contrast, the Tonga STEP fault dips southward but also has a northward component of motion with respect to the mantle. Asthenosphere thus flows southward beneath the down flexed Tonga STEP fault edge opposing further bending.  Subduction does not initiate at the Tonga STEP fault despite a ~100 Myr age contrast between the Pacific and north Fiji and Lau basin lithospheres. Since absolute plate motions reflect the sum of all forces acting on the entire lithospheric plate, a strong sea anchor mantle force may be generated at a STEP fault edge, initiating subduction (or inhibiting it), even where lithosphere is old, strong, and resists bending and without requiring large convergent forces between plates, overcoming these objections to previous models. 
    more » « less
  2. Abstract Plate reconstruction models are constructed to fit constraints such as magnetic anomalies, fracture zones, paleomagnetic poles, geological observations and seismic tomography. However, these models do not consider the physical equations of plate driving forces when reconstructing plate motion. This can potentially result in geodynamically-implausible plate motions, which has implications for a range of work based on plate reconstruction models. We present a new algorithm that calculates time-dependent slab pull, ridge push (GPE force) and mantle drag resistance for any topologically closed reconstruction, and evaluates the residuals—or missing components—required for torques to balance given our assumed plate driving force relationships. In all analyzed models, residual torques for the present-day are three orders of magnitude smaller than the typical driving torques for oceanic plates, but can be of the same order of magnitude back in time—particularly from 90 to 50 Ma. Using the Pacific plate as an example, we show how our algorithm can be used to identify areas and times with high residual torques, where either plate reconstructions have a high degree of geodynamic implausibility or our understanding of the underlying geodynamic forces is incomplete. We suggest strategies for plate model improvements and also identify times when other forces such as active mantle flow were likely important contributors. Our algorithm is intended as a tool to help assess and improve plate reconstruction models based on a transparent and expandable set of a priori dynamic constraints. 
    more » « less
  3. Variable viscosity in Earth’s mantle exerts a fundamental control on mantle convection and plate tectonics, yet rigorously constraining the underlying parameters has remained a challenge. Inverse methods have not been sufficiently robust to handle the severe viscosity gradients and nonlinearities (arising from dislocation creep and plastic failure) while simultaneously resolving the megathrust and bending slabs globally. Using global plate motions as constraints, we overcome these challenges by combining a scalable nonlinear Stokes solver that resolves the key tectonic features with an adjoint-based Bayesian approach. Assuming plate cooling, variations in the thickness of continental lithosphere, slabs, and broad scale lower mantle structure as well as a constant grain size through the bulk of the upper mantle, a good fit to global plate motions is found with a nonlinear upper mantle stress exponent of 2.43 ± 0.25 (mean ± SD). A relatively low yield stress of 151 ± 19 MPa is required for slabs to bend during subduction and transmit a slab pull that generates asymmetrical subduction. The recovered long-term strength of megathrusts (plate interfaces) varies between different subduction zones, with South America having a larger strength and Vanuatu and Central America having lower values with important implications for the stresses driving megathrust earthquakes. 
    more » « less
  4. 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
  5. 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