Cenozoic India-Asia collision driven by mantle dragging the cratonic root

The driving force behind the Cenozoic India-Asia collision remains elusive. Using global-scale geodynamic modeling, we find that the continuous motion of the Indian plate is driven by a prominent upper-mantle flow pushing the thick Indian lithospheric root, originated from the northward rollover of the detached Neo-Tethyan slab and sinking slabs below East Asia. The maximum mantle drag occurs within the strong Indian lithosphere and is comparable in magnitude to that of slab pull (1013 N m-1). The thick cratonic root enhances both lithosphere-asthenosphere coupling and upper-plate compressional stress, thereby sustaining the topography of Tibetan Plateau. We show that the calculated resistant force from the India-Asia plate boundary is also close to that due to the gravitational potential energy of Tibetan Plateau. Here, we demonstrate that this mantle flow is key for the formation of the Tibetan Plateau and represents part of a hemispheric convergent flow pattern centered on central Asia.

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.

Kimberlite eruptions driven by slab flux and subduction angle

Kimberlites are sourced from thermochemical upwellings which can transport diamonds to the surface of the crust. The majority of kimberlites preserved at the Earth's surface erupted between 250 and 50 million years ago, and have been attributed to changes in plate velocity or mantle plumes. However, these mechanisms fail to explain the presence of strong subduction signatures observed in some Cretaceous kimberlites. This raises the question whether there is a subduction process that unifies our understanding of the timing of kimberlite eruptions. We develop a novel formulation for calculating subduction angle based on trench migration, convergence rate, slab thickness and density to connect the influx of slab material into the mantle with the timing of kimberlite eruptions. We find that subduction angles combined with peaks in slab flux predict pulses of kimberlite eruptions. High rates of subducting slab material trigger mantle return flow that stimulates fertile reservoirs in the mantle. These convective instabilities transport slab-influenced melt to the surface at a distance inbound from the trench corresponding to the subduction angle. Our deep-time slab dip formulation has numerous potential applications including modelling the deep carbon and water cycles, and an improved understanding of subduction-related mineral deposits.

Oceanic crustal carbon cycle drives 26-million-year atmospheric carbon dioxide periodicities

Atmospheric carbon dioxide (CO₂) data for the last 420 million years (My) show long-term fluctuations related to supercontinent cycles as well as shorter cycles at 26 to 32 My whose origin is unknown. Periodicities of 26 to 30 My occur in diverse geological phenomena including mass extinctions, flood basalt volcanism, ocean anoxic events, deposition of massive evaporites, sequence boundaries, and orogenic events and have previously been linked to an extraterrestrial mechanism. The vast oceanic crustal carbon reservoir is an alternative potential driving force of climate fluctuations at these time scales, with hydrothermal crustal carbon uptake occurring mostly in young crust with a strong dependence on ocean bottom water temperature. We combine a global plate model and oceanic paleo-age grids with estimates of paleo-ocean bottom water temperatures to track the evolution of the oceanic crustal carbon reservoir over the past 230 My. We show that seafloor spreading rates as well as the storage, subduction, and emission of oceanic crustal and mantle CO₂ fluctuate with a period of 26 My. A connection with seafloor spreading rates and equivalent cycles in subduction zone rollback suggests that these periodicities are driven by the dynamics of subduction zone migration. The oceanic crust-mantle carbon cycle is thus a previously overlooked mechanism that connects plate tectonic pulsing with fluctuations in atmospheric carbon and surface environments.