We present a revised global plate motion model with continuously closing plate boundaries ranging from the Triassic at 230 Ma to the present day, assess differences among alternative absolute plate motion models, and review global tectonic events. Relatively high mean absolute plate motion rates of approximately 9–10 cm yr−1 between 140 and 120 Ma may be related to transient plate motion accelerations driven by the successive emplacement of a sequence of large igneous provinces during that time. An event at ∼100 Ma is most clearly expressed in the Indian Ocean and may reflect the initiation of Andean-style subduction along southern continental Eurasia, whereas an acceleration at ∼80 Ma of mean rates from 6 to 8 cm yr−1 reflects the initial northward acceleration of India and simultaneous speedups of plates in the Pacific. An event at ∼50 Ma expressed in relative, and some absolute, plate motion changes around the globe and in a reduction of global mean plate speeds from about 6 to 4–5 cm yr−1 indicates that an increase in collisional forces (such as the India–Eurasia collision) and ridge subduction events in the Pacific (such as the Izanagi–Pacific Ridge) play a significant role in modulating plate velocities.

No abstract

The breakup of Pangea in the Jurassic saw the opening of major ocean basins at the expense of older Tethyan and Pacific oceanic plates. Although the Tethyan seafloor spreading history has been lost to subduction, proxy indicators from multiple generations of Tethyan ribbon terranes, as well as the active margin geological histories of volcanism and ophiolite obduction events can be used to reconstruct these ancient oceanic plates. The plate reconstructions presented in this study reconcile observations from ocean basins and the onshore geological record to provide a regional synthesis, embedded in a global plate motion model, of the IndiaEurasia convergence history, the accretionary growth of Southeast Asia and the Tethyan-Pacific tectonic link through the New Guinea margin. The global plate motion model presented in this study captures the timedependent evolution of plates and their tectonic boundaries since 160 Ma, which are assimilated as surface boundary conditions for numerical experiments of mantle convection. We evaluate subducted slab locations and geometries predicted by forward mantle flow models against P-and S-wave seismic tomography models. This approach harnesses modern plate reconstruction techniques, mantle convection models with imposed one-sided subduction, and constraints from the surface geology to address a number of unresolved Tethyan geodynamic controversies. Our synthesis reveals that north-dipping subduction beneath Eurasia in the latest Jurassic consumed the Meso-Tethys, and suggests that northward slab pull opened the younger Neo-Tethyan ocean basin from ~ 155 Ma. We model the rifting of ' Argoland', representing the East Java and West Sulawesi continental fragments, as a northward transfer of continental terranes in the latest Jurassic from the northwest Australian shelf -likely colliding first with parts of the Woyla intra-oceanic arc in the mid-Cretaceous, and accreting to the Borneo (Sundaland) core by ~ 80 Ma. The Neo-Tethyan ridge was likely consumed along an intra-oceanic subduction zone south of Eurasia from ~ 105 Ma, leading to a major change in the motion of the Indian Plate by ~ 100 Ma, as observed in the Wharton Basin fracture zone bends. We investigate the geodynamic consequences of long-lived intra-oceanic subduction within the Neo-Tethys, requiring a twostage India-Eurasia collision involving first contact between Greater India and the Kohistan-Ladakh Arc sometime between ~ 60 and 50 Ma, followed by continent-continent collision from ~ 47 Ma. Our models suggest that the Sunda slab kink beneath northwest Sumatra in the mantle transition zone results from the rotation and extrusion of Indochina from ~ 30 Ma. Our results are also the first to reproduce the enigmatic Proto South China Sea slab beneath northern Borneo, as well as the Tethyan/Woyla slab that is predicted at mid-mantle depths south of Sumatra. Further east, our revised reconstructions of the New Guinea margin, notably the evolution of the Sepik composite terrane and the Maramuni subduction zone, produce...

[1] We present a global community data set of fracture zones (FZs), discordant zones, propagating ridges, V-shaped structures and extinct ridges, digitized from vertical gravity gradient (VGG) maps. We use a new semi-automatic FZ tracking program to test the precision of our hand-digitized traces and find a Mean Absolute Deviation of less than 3.4 km from the raw VGG minima that most clearly delineate each feature, and less than 5.4 km from the FZ location predicted by fitting model profiles to the VGG data that represent the morphology of the individual FZs. These offsets are small considering gravity data only provide an approximation for the underlying basement morphology. We further investigate the origin of non-FZ seafloor fabric by combining published abyssal hill heights computed from gravity anomalies with global half-spreading rates. A residual abyssal hill height grid, with spreading rate effects removed, combined with our interpreted tectonic fabric reveals several types of seafloor fabric distinct from typical abyssal hills. Where discordant zones do not overprint abyssal hill signals, residual abyssal hill height anomalies correspond to seafloor that accreted near mantle thermal anomalies or zones of melt-depletion. Our analysis reveals several areas where residual abyssal hill height anomalies reflect pseudo-faults and extinct ridges associated with ridge propagation and/or microplate formation in the southern Pacific Ocean.

We present an updated oceanic crustal age grid and a set of complementary grids including spreading rate, asymmetry, direction, and obliquity. Our data set is based on a selected set of magnetic anomaly identifications and the plate tectonic model of Müller et al. (2019, 10.1029/2018TC005462). We find the mean age of oceanic crust is 64.2 Myr, slightly older than previous estimates, mainly due to the inclusion of pockets of Mesozoic aged crust in the Atlantic and Mediterranean and improvements to the Jurassic Pacific triangle. This older crust is partly compensated by additional Cenozoic-aged back-arc basin crust not included in previous models. The distribution of spreading modes based on area of preserved crust is relatively equal between slow (20-55 mm/yr) and fast (75-180 mm/yr) spreading systems at 33% and 39%, respectively. Crust transitional between fast and slow, or intermediate systems (55-75 mm/yr), cover 20% of the preserved ocean floor with much smaller proportions of crust formed at ultraslow (5%) and super-fast (3%) spreading systems. Slow and intermediate spreading systems exhibit the most stable behavior in terms of spreading asymmetry and obliquity, with the widest distribution of obliquities occurring at ultraslow spreading systems, consistent with present-day observations. Our confidence grid provides a complementary resource for nonexperts to identify those parts of the age grid that are least well constrained. Our grids in 6, 2, and 1 arc min resolution as well as our python workflow, isopolate, used to compute our data sets are freely available in online repositories and on the GPlates data portal. Plain Language Summary The age of the ocean floor is one of the fundamental data sets for understanding plate tectonic processes. We have generated an updated present-day age grid of oceanic crust by incorporating recent improvements to the underlying plate motion model, magnetic anomaly identifications, and fracture zones. The resultant mean age of the ocean floor is~64 million years, slightly older than previous studies. We also present data sets showing various seafloor spreading parameters such as spreading rate, asymmetry, direction, and obliquity. Slow and intermediate seafloor spreading rates produce the most stable spreading behavior in terms of asymmtery and obliquity. Our data sets and workflows are freely available and can be used for multidisciplinary scientific studies. Our maps are an effective visual tool to understand plate tectonic concepts and can be used for education and outreach.

No abstract

No abstract

The geoscience community is increasingly utilizing seismic tomography to interpret mantle heterogeneity and its links to past tectonic and geodynamic processes. To assess the robustness and distribution of positive seismic anomalies, inferred as subducted slabs, we create a set of vote maps for the lower mantle with 14 global P-wave or S-wave tomography models. Based on a depth-dependent threshold metric, an average of 20% of any given tomography model depth is identified as a potential slab. However, upon combining the 14 models, the most consistent positive wavespeed features are identified by an increasing vote count. An overall peak in the most robust anomalies is found between 1000–1400 km depth, followed by a decline to a minimum around 2000 km. While this trend could reflect reduced tomographic resolution in the middle mantle, we show that it may alternatively relate to real changes in the time-dependent subduction flux and/or a mid-lower mantle viscosity increase. An apparent secondary peak in agreement below 2500 km depth may reflect the degree-two lower mantle slow seismic structures. Vote maps illustrate the potential shortcomings of using a limited number or type of tomography models and slab threshold criteria.