Earth's interior dynamics drive marine fossil diversity cycles of tens of millions of years.

The fossil record reveals that biotic diversity has fluctuated quasi-cyclically through geological time. However, the causal mechanisms of biotic diversity cycles remain unexplained. Here, we highlight a common, correlatable 36 ± 1 Myr (million years) cycle in the diversity of marine genera as well as in tectonic, sea-level, and macrostratigraphic data over the past 250 Myr of Earth history. The prominence of the 36 ± 1 Myr cycle in tectonic data favors a common-cause mechanism, wherein geological forcing mechanisms drive patterns in both biological diversity and the preserved rock record. In particular, our results suggest that a 36 ± 1 Myr tectono-eustatically driven sea-level cycle may originate from the interaction between the convecting mantle and subducting slabs, thereby pacing mantle-lithospheric deep-water recycling. The 36 ± 1 Myr tectono-eustatic driver of biodiversity is likely related to cyclic continental inundations, with expanding and contracting ecological niches on shelves and in epeiric seas.

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.

Post-extinction recovery of the Phanerozoic oceans and biodiversity hotspots.

The fossil record of marine invertebrates has long fuelled the debate as to whether or not there are limits to global diversity in the sea1-5. Ecological theory states that, as diversity grows and ecological niches are filled, the strengthening of biological interactions imposes limits on diversity6,7. However, the extent to which biological interactions have constrained the growth of diversity over evolutionary time remains an open question1-5,8-11. Here we present a regional diversification model that reproduces the main Phanerozoic eon trends in the global diversity of marine invertebrates after imposing mass extinctions. We find that the dynamics of global diversity are best described by a diversification model that operates widely within the exponential growth regime of a logistic function. A spatially resolved analysis of the ratio of diversity to carrying capacity reveals that less than 2% of the global flooded continental area throughout the Phanerozoic exhibits diversity levels approaching ecological saturation. We attribute the overall increase in global diversity during the Late Mesozoic and Cenozoic eras to the development of diversity hotspots under prolonged conditions of Earth system stability and maximum continental fragmentation. We call this the 'diversity hotspots hypothesis', which we propose as a non-mutually exclusive alternative to the hypothesis that the Mesozoic marine revolution led this macroevolutionary trend12,13.

Evolution of Earth's tectonic carbon conveyor belt.

Concealed deep beneath the oceans is a carbon conveyor belt, propelled by plate tectonics. Our understanding of its modern functioning is underpinned by direct observations, but its variability through time has been poorly quantified. Here we reconstruct oceanic plate carbon reservoirs and track the fate of subducted carbon using thermodynamic modelling. In the Mesozoic era, 250 to 66 million years ago, plate tectonic processes had a pivotal role in driving climate change. Triassic-Jurassic period cooling correlates with a reduction in solid Earth outgassing, whereas Cretaceous period greenhouse conditions can be linked to a doubling in outgassing, driven by high-speed plate tectonics. The associated 'carbon subduction superflux' into the subcontinental mantle may have sparked North American diamond formation. In the Cenozoic era, continental collisions slowed seafloor spreading, reducing tectonically driven outgassing, while deep-sea carbonate sediments emerged as the Earth's largest carbon sink. Subduction and devolatilization of this reservoir beneath volcanic arcs led to a Cenozoic increase in carbon outgassing, surpassing mid-ocean ridges as the dominant source of carbon emissions 20 million years ago. An increase in solid Earth carbon emissions during Cenozoic cooling requires an increase in continental silicate weathering flux to draw down atmospheric carbon dioxide, challenging previous views and providing boundary conditions for future carbon cycle models.

Constraining the response of continental-scale groundwater flow to climate change.

Numerical models of groundwater flow play a critical role for water management scenarios under climate extremes. Large-scale models play a key role in determining long range flow pathways from continental interiors to the oceans, yet struggle to simulate the local flow patterns offered by small-scale models. We have developed a highly scalable numerical framework to model continental groundwater flow which capture the intricate flow pathways between deep aquifers and the near-surface. The coupled thermal-hydraulic basin structure is inferred from hydraulic head measurements, recharge estimates from geochemical proxies, and borehole temperature data using a Bayesian framework. We use it to model the deep groundwater flow beneath the Sydney-Gunnedah-Bowen Basin, part of Australia's largest aquifer system. Coastal aquifers have flow rates of up to 0.3 m/day, and a corresponding groundwater residence time of just 2,000 years. In contrast, our model predicts slow flow rates of 0.005 m/day for inland aquifers, resulting in a groundwater residence time of [Formula: see text] 400,000 years. Perturbing the model to account for a drop in borehole water levels since 2000, we find that lengthened inland flow pathways depart significantly from pre-2000 streamlines as groundwater is drawn further from recharge zones in a drying climate. Our results illustrate that progressively increasing water extraction from inland aquifers may permanently alter long-range flow pathways. Our open-source modelling approach can be extended to any basin and may help inform policies on the sustainable management of groundwater.

The evolution of basal mantle structure in response to supercontinent aggregation and dispersal.

Seismic studies have revealed two Large Low-Shear Velocity Provinces (LLSVPs) in the lowermost mantle. Whether these structures remain stable over time or evolve through supercontinent cycles is debated. Here we analyze a recently published mantle flow model constrained by a synthetic plate motion model extending back to one billion years ago, to investigate how the mantle evolves in response to changing plate configurations. Our model predicts that sinking slabs segment the basal thermochemical structure below an assembling supercontinent, and that this structure eventually becomes unified due to slab push from circum-supercontinental subduction. In contrast, the basal thermochemical structure below the superocean is generally coherent due to the persistence of a superocean in our imposed plate reconstruction. The two antipodal basal thermochemical structures exchange material several times when part of one of the structures is carved out and merged with the other one, similarly to "exotic" tectonic terranes. Plumes mostly rise from thick basal thermochemical structures and in some instances migrate from the edges towards the interior of basal thermochemical structures due to slab push. Our results suggest that the topography of basal structures and distribution of plumes change over time due to the changing subduction network over supercontinent cycles.

Intraplate volcanism triggered by bursts in slab flux.

Long-lived, widespread intraplate volcanism without age progression is one of the most controversial features of plate tectonics. Previously proposed edge-driven convection, asthenospheric shear, and lithospheric detachment fail to explain the ~5000-km-wide intraplate volcanic province from eastern Australia to Zealandia. We model the subducted slab volume over 100 million years and find that slab flux drives volcanic eruption frequency, indicating stimulation of an enriched mantle transition zone reservoir. Volcanic isotope geochemistry allows us to distinguish a high-µ (HIMU) reservoir [>1 billion years (Ga) old] in the slab-poor south, from a northern EM1/EM2 reservoir, reflecting a more recent voluminous influx of oceanic lithosphere into the mantle transition zone. We provide a unified theory linking plate boundary and slab volume reconstructions to upper mantle reservoirs and intraplate volcano geochemistry.

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

Atmospheric carbon dioxide (CO2) 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 CO2 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.

Abrupt plate accelerations shape rifted continental margins.

Rifted margins are formed by persistent stretching of continental lithosphere until breakup is achieved. It is well known that strain-rate-dependent processes control rift evolution, yet quantified extension histories of Earth's major passive margins have become available only recently. Here we investigate rift kinematics globally by applying a new geotectonic analysis technique to revised global plate reconstructions. We find that rifted margins feature an initial, slow rift phase (less than ten millimetres per year, full rate) and that an abrupt increase of plate divergence introduces a fast rift phase. Plate acceleration takes place before continental rupture and considerable margin area is created during each phase. We reproduce the rapid transition from slow to fast extension using analytical and numerical modelling with constant force boundary conditions. The extension models suggest that the two-phase velocity behaviour is caused by a rift-intrinsic strength--velocity feedback, which can be robustly inferred for diverse lithosphere configurations and rheologies. Our results explain differences between proximal and distal margin areas and demonstrate that abrupt plate acceleration during continental rifting is controlled by the nonlinear decay of the resistive rift strength force. This mechanism provides an explanation for several previously unexplained rapid absolute plate motion changes, offering new insights into the balance of plate driving forces through time.

Subduction controls the distribution and fragmentation of Earth's tectonic plates.

The theory of plate tectonics describes how the surface of Earth is split into an organized jigsaw of seven large plates of similar sizes and a population of smaller plates whose areas follow a fractal distribution. The reconstruction of global tectonics during the past 200 million years suggests that this layout is probably a long-term feature of Earth, but the forces governing it are unknown. Previous studies, primarily based on the statistical properties of plate distributions, were unable to resolve how the size of the plates is determined by the properties of the lithosphere and the underlying mantle convection. Here we demonstrate that the plate layout of Earth is produced by a dynamic feedback between mantle convection and the strength of the lithosphere. Using three-dimensional spherical models of mantle convection that self-consistently produce the plate size­frequency distribution observed for Earth, we show that subduction geometry drives the tectonic fragmentation that generates plates. The spacing between the slabs controls the layout of large plates, and the stresses caused by the bending of trenches break plates into smaller fragments. Our results explain why the fast evolution in small back-arc plates reflects the marked changes in plate motions during times of major reorganizations. Our study opens the way to using convection simulations with plate-like behaviour to unravel how global tectonics and mantle convection are dynamically connected.

A rapid burst in hotspot motion through the interaction of tectonics and deep mantle flow.

Volcanic hotspot tracks featuring linear progressions in the age of volcanism are typical surface expressions of plate tectonic movement on top of narrow plumes of hot material within Earth's mantle. Seismic imaging reveals that these plumes can be of deep origin--probably rooted on thermochemical structures in the lower mantle. Although palaeomagnetic and radiometric age data suggest that mantle flow can advect plume conduits laterally, the flow dynamics underlying the formation of the sharp bend occurring only in the Hawaiian-Emperor hotspot track in the Pacific Ocean remains enigmatic. Here we present palaeogeographically constrained numerical models of thermochemical convection and demonstrate that flow in the deep lower mantle under the north Pacific was anomalously vigorous between 100 million years ago and 50 million years ago as a consequence of long-lasting subduction systems, unlike those in the south Pacific. These models show a sharp bend in the Hawaiian-Emperor hotspot track arising from the interplay of plume tilt and the lateral advection of plume sources. The different trajectories of the Hawaiian and Louisville hotspot tracks arise from asymmetric deformation of thermochemical structures under the Pacific between 100 million years ago and 50 million years ago. This asymmetric deformation waned just before the Hawaiian-Emperor bend developed, owing to flow in the deepest lower mantle associated with slab descent in the north and south Pacific.

The GPlates Portal: Cloud-Based Interactive 3D Visualization of Global Geophysical and Geological Data in a Web Browser.

The pace of scientific discovery is being transformed by the availability of 'big data' and open access, open source software tools. These innovations open up new avenues for how scientists communicate and share data and ideas with each other and with the general public. Here, we describe our efforts to bring to life our studies of the Earth system, both at present day and through deep geological time. The GPlates Portal (portal.gplates.org) is a gateway to a series of virtual globes based on the Cesium Javascript library. The portal allows fast interactive visualization of global geophysical and geological data sets, draped over digital terrain models. The globes use WebGL for hardware-accelerated graphics and are cross-platform and cross-browser compatible with complete camera control. The globes include a visualization of a high-resolution global digital elevation model and the vertical gradient of the global gravity field, highlighting small-scale seafloor fabric such as abyssal hills, fracture zones and seamounts in unprecedented detail. The portal also features globes portraying seafloor geology and a global data set of marine magnetic anomaly identifications. The portal is specifically designed to visualize models of the Earth through geological time. These space-time globes include tectonic reconstructions of the Earth's gravity and magnetic fields, and several models of long-wavelength surface dynamic topography through time, including the interactive plotting of vertical motion histories at selected locations. The globes put the on-the-fly visualization of massive data sets at the fingertips of end-users to stimulate teaching and learning and novel avenues of inquiry.

Marine geophysics. New global marine gravity model from CryoSat-2 and Jason-1 reveals buried tectonic structure.

Gravity models are powerful tools for mapping tectonic structures, especially in the deep ocean basins where the topography remains unmapped by ships or is buried by thick sediment. We combined new radar altimeter measurements from satellites CryoSat-2 and Jason-1 with existing data to construct a global marine gravity model that is two times more accurate than previous models. We found an extinct spreading ridge in the Gulf of Mexico, a major propagating rift in the South Atlantic Ocean, abyssal hill fabric on slow-spreading ridges, and thousands of previously uncharted seamounts. These discoveries allow us to understand regional tectonic processes and highlight the importance of satellite-derived gravity models as one of the primary tools for the investigation of remote ocean basins.

How supercontinents and superoceans affect seafloor roughness.

Seafloor roughness varies considerably across the world's ocean basins and is fundamental to controlling the circulation and mixing of heat in the ocean and dissipating eddy kinetic energy. Models derived from analyses of active mid-ocean ridges suggest that ocean floor roughness depends on seafloor spreading rates, with rougher basement forming below a half-spreading rate threshold of 30-35 mm yr(-1) (refs 4, 5), as well as on the local interaction of mid-ocean ridges with mantle plumes or cold-spots. Here we present a global analysis of marine gravity-derived roughness, sediment thickness, seafloor isochrons and palaeo-spreading rates of Cretaceous to Cenozoic ridge flanks. Our analysis reveals that, after eliminating effects related to spreading rate and sediment thickness, residual roughness anomalies of 5-20 mGal remain over large swaths of ocean floor. We found that the roughness as a function of palaeo-spreading directions and isochron orientations indicates that most of the observed excess roughness is not related to spreading obliquity, as this effect is restricted to relatively rare occurrences of very high obliquity angles (>45 degrees ). Cretaceous Atlantic ocean floor, formed over mantle previously overlain by the Pangaea supercontinent, displays anomalously low roughness away from mantle plumes and is independent of spreading rates. We attribute this observation to a sub-Pangaean supercontinental mantle temperature anomaly leading to slightly thicker than normal Late Jurassic and Cretaceous Atlantic crust, reduced brittle fracturing and smoother basement relief. In contrast, ocean crust formed above Pacific superswells, probably reflecting metasomatized lithosphere underlain by mantle at only slightly elevated temperatures, is not associated with basement roughness anomalies. These results highlight a fundamental difference in the nature of large-scale mantle upwellings below supercontinents and superoceans, and their impact on oceanic crustal accretion.

Long-term sea-level fluctuations driven by ocean basin dynamics.

Earth's long-term sea-level history is characterized by widespread continental flooding in the Cretaceous period (approximately 145 to 65 million years ago), followed by gradual regression of inland seas. However, published estimates of the Late Cretaceous sea-level high differ by half an order of magnitude, from approximately 40 to approximately 250 meters above the present level. The low estimate is based on the stratigraphy of the New Jersey margin. By assimilating marine geophysical data into reconstructions of ancient ocean basins, we model a Late Cretaceous sea level that is 170 (85 to 270) meters higher than it is today. We use a mantle convection model to suggest that New Jersey subsided by 105 to 180 meters in the past 70 million years because of North America's westward passage over the subducted Farallon plate. This mechanism reconciles New Jersey margin-based sea-level estimates with ocean basin reconstructions.