A Global Plate Model Including Lithospheric Deformation Along Major Rifts and Orogens Since the Triassic

Global deep‐time plate motion models have traditionally followed a classical rigid plate approach, even though plate deformation is known to be significant. Here we present a global Mesozoic–Cenozoic deforming plate motion model that captures the progressive extension of all continental margins since the initiation of rifting within Pangea at ~240 Ma. The model also includes major failed continental rifts and compressional deformation along collision zones. The outlines and timing of regional deformation episodes are reconstructed from a wealth of published regional tectonic models and associated geological and geophysical data. We reconstruct absolute plate motions in a mantle reference frame with a joint global inversion using hot spot tracks for the last 80 million years and minimizing global trench migration velocities and net lithospheric rotation. In our optimized model, net rotation is consistently below 0.2°/Myr, and trench migration scatter is substantially reduced. Distributed plate deformation reaches a Mesozoic peak of 30 × 106 km2 in the Late Jurassic (~160–155 Ma), driven by a vast network of rift systems. After a mid‐Cretaceous drop in deformation, it reaches a high of 48 x 106 km2 in the Late Eocene (~35 Ma), driven by the progressive growth of plate collisions and the formation of new rift systems. About a third of the continental crustal area has been deformed since 240 Ma, partitioned roughly into 65% extension and 35% compression. This community plate model provides a framework for building detailed regional deforming plate networks and form a constraint for models of basin evolution and the plate‐mantle system

Decoding Earth’s plate tectonic history using sparse geochemical data

Accurately mapping plate boundary types and locations through time is essential for understanding the evolution of the plate-mantle system and the exchange of material between the solid Earth and surface environments. However, the complexity of the Earth system and the cryptic nature of the geological record make it difficult to discriminate tectonic environments through deep time. Here we present a new method for identifying tectonic paleo-environments on Earth through a data mining approach using global geochemical data. We first fingerprint a variety of present-day tectonic environments utilising up to 136 geochemical data attributes in any available combination. A total of 38301 geochemical analyses from basalts aged from 5–0 Ma together with a well-established plate reconstruction model are used to construct a suite of discriminatory models for the first order tectonic environments of subduction and mid-ocean ridge as distinct from intraplate hotspot oceanic environments, identifying 41, 35, and 39 key discriminatory geochemical attributes, respectively. After training and validation, our model is applied to a global geochemical database of 1547 basalt samples of unknown tectonic origin aged between 1000–410 Ma, a relatively ill-constrained period of Earth's evolution following the breakup of the Rodinia supercontinent, producing 56 unique global tectonic environment predictions throughout the Neoproterozoic and Early Paleozoic. Predictions are used to discriminate between three alternative published Rodinia configuration models, identifying the model demonstrating the closest spatio-temporal consistency with the basalt record, and emphasizing the importance of integrating geochemical data into plate reconstructions. Our approach offers an extensible framework for constructing full-plate, deep-time reconstructions capable of assimilating a broad range of geochemical and geological observations, enabling next generation Earth system models

Constraining Absolute Plate Motions Since the Triassic

The absolute motion of tectonic plates since Pangea can be derived from observations of hotspot trails, paleomagnetism, or seismic tomography. However, fitting observations is typically carried out in isolation without consideration for the fit to unused data or whether the resulting plate motions are geodynamically plausible. Through the joint evaluation of global hotspot track observations (for times <80 Ma), first‐order estimates of net lithospheric rotation (NLR), and parameter estimation for paleo–trench migration (TM), we present a suite of geodynamically consistent, data‐optimized global absolute reference frames from 220 Ma to the present. Each absolute plate motion (APM) model was evaluated against six published APM models, together incorporating the full range of primary data constraints. Model performance for published and new models was quantified through a standard statistical analyses using three key diagnostic global metrics: root‐mean square plate velocities, NLR characteristics, and TM behavior. Additionally, models were assessed for consistency with published global paleomagnetic data and for ages <80 Ma for predicted relative hotspot motion, track geometry, and time dependence. Optimized APM models demonstrated significantly improved global fit with geological and geophysical observations while performing consistently with geodynamic constraints. Critically, APM models derived by limiting average rates of NLR to ~0.05°/Myr and absolute TM velocities to ~27‐mm/year fit geological observations including hotspot tracks. This suggests that this range of NLR and TM estimates may be appropriate for Earth over the last 220 Myr, providing a key step toward the practical integration of numerical geodynamics into plate tectonic reconstructions

A tectonic-rules-based mantle reference frame since 1 billion years ago - implications for supercontinent cycles and plate-mantle system evolution

Understanding the long-term evolution of Earth\u27s plate-mantle system is reliant on absolute plate motion models in a mantle reference frame, but such models are both difficult to construct and controversial. We present a tectonic-rules-based optimization approach to construct a plate motion model in a mantle reference frame covering the last billion years and use it as a constraint for mantle flow models. Our plate motion model results in net lithospheric rotation consistently below 0.25g g€Myr-1, in agreement with mantle flow models, while trench motions are confined to a relatively narrow range of -2 to +2g€cmg€yr-1 since 320g€Ma, during Pangea stability and dispersal. In contrast, the period from 600 to 320g€Ma, nicknamed the zippy tricentenary here, displays twice the trench motion scatter compared to more recent times, reflecting a predominance of short and highly mobile subduction zones. Our model supports an orthoversion evolution from Rodinia to Pangea with Pangea offset approximately 90g eastwards relative to Rodinia - this is the opposite sense of motion compared to a previous orthoversion hypothesis based on paleomagnetic data. In our coupled plate-mantle model a broad network of basal mantle ridges forms between 1000 and 600g€Ma, reflecting widely distributed subduction zones. Between 600 and 500g€Ma a short-lived degree-2 basal mantle structure forms in response to a band of subduction zones confined to low latitudes, generating extensive antipodal lower mantle upwellings centred at the poles. Subsequently, the northern basal structure migrates southward and evolves into a Pacific-centred upwelling, while the southern structure is dissected by subducting slabs, disintegrating into a network of ridges between 500 and 400g€Ma. From 400 to 200g€Ma, a stable Pacific-centred degree-1 convective planform emerges. It lacks an antipodal counterpart due to the closure of the Iapetus and Rheic oceans between Laurussia and Gondwana as well as due to coeval subduction between Baltica and Laurentia and around Siberia, populating the mantle with slabs until 320g€Ma when Pangea is assembled. A basal degree-2 structure forms subsequent to Pangea breakup, after the influence of previously subducted slabs in the African hemisphere on the lowermost mantle structure has faded away. This succession of mantle states is distinct from previously proposed mantle convection models. We show that the history of plume-related volcanism is consistent with deep plumes associated with evolving basal mantle structures. This Solid Earth Evolution Model for the last 1000 million years (SEEM1000) forms the foundation for a multitude of spatio-temporal data analysis approaches

Decoding Earth’s plate tectonic history using sparse geochemical data

Accurately mapping plate boundary types and locations through time is essential for understanding the evolution of the plate-mantle system and the exchange of material between the solid Earth and surface environments. However, the complexity of the Earth system and the cryptic nature of the geological record make it difficult to discriminate tectonic environments through deep time. Here we present a new method for identifying tectonic paleo-environments on Earth through a data mining approach using global geochemical data. We first fingerprint a variety of present-day tectonic environments utilising up to 136 geochemical data attributes in any available combination. A total of 38301 geochemical analyses from basalts aged from 5–0 Ma together with a well-established plate reconstruction model are used to construct a suite of discriminatory models for the first order tectonic environments of subduction and mid-ocean ridge as distinct from intraplate hotspot oceanic environments, identifying 41, 35, and 39 key discriminatory geochemical attributes, respectively. After training and validation, our model is applied to a global geochemical database of 1547 basalt samples of unknown tectonic origin aged between 1000–410 Ma, a relatively ill-constrained period of Earth's evolution following the breakup of the Rodinia supercontinent, producing 56 unique global tectonic environment predictions throughout the Neoproterozoic and Early Paleozoic. Predictions are used to discriminate between three alternative published Rodinia configuration models, identifying the model demonstrating the closest spatio-temporal consistency with the basalt record, and emphasizing the importance of integrating geochemical data into plate reconstructions. Our approach offers an extensible framework for constructing full-plate, deep-time reconstructions capable of assimilating a broad range of geochemical and geological observations, enabling next generation Earth system models

A global plate model including lithospheric deformation along major rifts and orogens since the Triassic

Global deep‐time plate motion models have traditionally followed a classical rigid plate approach, even though plate deformation is known to be significant. Here we present a global Mesozoic–Cenozoic deforming plate motion model that captures the progressive extension of all continental margins since the initiation of rifting within Pangea at ~240 Ma. The model also includes major failed continental rifts and compressional deformation along collision zones. The outlines and timing of regional deformation episodes are reconstructed from a wealth of published regional tectonic models and associated geological and geophysical data. We reconstruct absolute plate motions in a mantle reference frame with a joint global inversion using hot spot tracks for the last 80 million years and minimizing global trench migration velocities and net lithospheric rotation. In our optimized model, net rotation is consistently below 0.2°/Myr, and trench migration scatter is substantially reduced. Distributed plate deformation reaches a Mesozoic peak of 30 × 10^6 km^2 in the Late Jurassic (~160–155 Ma), driven by a vast network of rift systems. After a mid‐Cretaceous drop in deformation, it reaches a high of 48 x 10^6 km^2 in the Late Eocene (~35 Ma), driven by the progressive growth of plate collisions and the formation of new rift systems. About a third of the continental crustal area has been deformed since 240 Ma, partitioned roughly into 65% extension and 35% compression. This community plate model provides a framework for building detailed regional deforming plate networks and form a constraint for models of basin evolution and the plate‐mantle system

A global plate model including lithospheric deformation along major rifts and orogens since the Triassic

Global deep‐time plate motion models have traditionally followed a classical rigid plate approach, even though plate deformation is known to be significant. Here we present a global Mesozoic–Cenozoic deforming plate motion model that captures the progressive extension of all continental margins since the initiation of rifting within Pangea at ~240 Ma. The model also includes major failed continental rifts and compressional deformation along collision zones. The outlines and timing of regional deformation episodes are reconstructed from a wealth of published regional tectonic models and associated geological and geophysical data. We reconstruct absolute plate motions in a mantle reference frame with a joint global inversion using hot spot tracks for the last 80 million years and minimizing global trench migration velocities and net lithospheric rotation. In our optimized model, net rotation is consistently below 0.2°/Myr, and trench migration scatter is substantially reduced. Distributed plate deformation reaches a Mesozoic peak of 30 × 10^6 km^2 in the Late Jurassic (~160–155 Ma), driven by a vast network of rift systems. After a mid‐Cretaceous drop in deformation, it reaches a high of 48 x 10^6 km^2 in the Late Eocene (~35 Ma), driven by the progressive growth of plate collisions and the formation of new rift systems. About a third of the continental crustal area has been deformed since 240 Ma, partitioned roughly into 65% extension and 35% compression. This community plate model provides a framework for building detailed regional deforming plate networks and form a constraint for models of basin evolution and the plate‐mantle system

Assessing plate reconstruction models using plate driving force consistency tests

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

Assessing plate reconstruction models using plate driving force consistency tests

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

A Global Plate Model Including Lithospheric Deformation Along Major Rifts and Orogens Since the Triassic

Global deep‐time plate motion models have traditionally followed a classical rigid plate approach, even though plate deformation is known to be significant. Here we present a global Mesozoic–Cenozoic deforming plate motion model that captures the progressive extension of all continental margins since the initiation of rifting within Pangea at ~240 Ma. The model also includes major failed continental rifts and compressional deformation along collision zones. The outlines and timing of regional deformation episodes are reconstructed from a wealth of published regional tectonic models and associated geological and geophysical data. We reconstruct absolute plate motions in a mantle reference frame with a joint global inversion using hot spot tracks for the last 80 million years and minimizing global trench migration velocities and net lithospheric rotation. In our optimized model, net rotation is consistently below 0.2°/Myr, and trench migration scatter is substantially reduced. Distributed plate deformation reaches a Mesozoic peak of 30 × 10⁶ km² in the Late Jurassic (~160–155 Ma), driven by a vast network of rift systems. After a mid‐Cretaceous drop in deformation, it reaches a high of 48 x 10⁶ km² in the Late Eocene (~35 Ma), driven by the progressive growth of plate collisions and the formation of new rift systems. About a third of the continental crustal area has been deformed since 240 Ma, partitioned roughly into 65% extension and 35% compression. This community plate model provides a framework for building detailed regional deforming plate networks and form a constraint for models of basin evolution and the plate‐mantle system