Subduction controls the distribution and fragmentation of Earth’s tectonic plates

Spreading ridge migration enabled by plume-ridge de-anchoring

It has long been recognised that spreading ridges are kept in place by competing subduction forces that drive plate motions. Asymmetric strain rates pull spreading ridges in the direction of the strongest slab pull force, which partially explains why spreading ridges can migrate vast distances. However, the interaction between mantle plumes and spreading ridges plays a relatively unknown role on the evolution of plate boundaries. Using a numerical model of mantle convection, we show that plumes with high buoyancy flux (>3000 kg/s) can capture spreading ridges within a 1000 km radius and anchor them in place. Exceptionally high buoyancy fluxes may fragment the overriding plate into smaller plates to accommodate more efficient plate motion. If the plume buoyancy flux wanes below 1000 kg/s the ridge may be de-anchored, leading to rapid ridge migration rates when combined with asymmetric plate boundary forces. Our results show that plume-ridge de-anchoring may have contributed to the rapid migration of the SE Indian Ridge from 43 million years ago (Ma) due to waning buoyancy flux from the Kerguelen plume, supported by magma flux estimates and radiogenic isotope geochemistry of eruption products. The plume-ridge de-anchoring mechanism we have identified has global implications for the evolution of plate boundaries near mantle plumes.

Artificial intelligence for geoscience: Progress, challenges, and perspectives

This paper explores the evolution of geoscientific inquiry, tracing the progression from traditional physics-based models to modern data-driven approaches facilitated by significant advancements in artificial intelligence (AI) and data collection techniques. Traditional models, which are grounded in physical and numerical frameworks, provide robust explanations by explicitly reconstructing underlying physical processes. However, their limitations in comprehensively capturing Earth's complexities and uncertainties pose challenges in optimization and real-world applicability. In contrast, contemporary data-driven models, particularly those utilizing machine learning (ML) and deep learning (DL), leverage extensive geoscience data to glean insights without requiring exhaustive theoretical knowledge. ML techniques have shown promise in addressing Earth science-related questions. Nevertheless, challenges such as data scarcity, computational demands, data privacy concerns, and the "black-box" nature of AI models hinder their seamless integration into geoscience. The integration of physics-based and data-driven methodologies into hybrid models presents an alternative paradigm. These models, which incorporate domain knowledge to guide AI methodologies, demonstrate enhanced efficiency and performance with reduced training data requirements. This review provides a comprehensive overview of geoscientific research paradigms, emphasizing untapped opportunities at the intersection of advanced AI techniques and geoscience. It examines major methodologies, showcases advances in large-scale models, and discusses the challenges and prospects that will shape the future landscape of AI in geoscience. The paper outlines a dynamic field ripe with possibilities, poised to unlock new understandings of Earth's complexities and further advance geoscience exploration.

Coupled fates of Earth's mantle and core: Early sluggish-lid tectonics and a long-lived geodynamo

Conventional Earth evolution models are unable to simultaneously reproduce two fundamental observations: the mantle's secular temperature record and a long-lived geodynamo before inner core nucleation. Today, plate tectonics efficiently cools the mantle, but if assumed to operate throughout Earth's history, past mantle temperature and plate motion become unrealistically high. Through coupled core-mantle modeling that self-consistently predicts multiple mantle convection regimes, we show that over most of the Precambrian, Earth likely operated in a distinct "sluggish-lid" tectonic mode, characterized by partial decoupling between the lithosphere and mantle. This dominant early regime is due to a hotter Earth and the presence of the asthenosphere. This mode regulates the core-mantle boundary heat flow, which powers the geodynamo before inner core nucleation. Both sluggish-lid tectonics and a long-lived dynamo demonstrate the inextricably connected paths of the core-mantle system. Moreover, our simulations simultaneously satisfy diverse geological observations and are consistent with emerging interpretations of such records.

Visibility Graph Analysis of the Seismic Activity of Three Areas of the Cocos Plate Mexican Subduction Where the Last Three Large Earthquakes (M > 7) Occurred in 2017 and 2022

The understanding of the dynamical behavior of seismic phenomena is currently an open problem, mainly because seismic series can be considered to be produced by phenomena exhibiting dynamic phase transitions; that is, with some complexity. For this purpose, the Middle America Trench in central Mexico is considered a natural laboratory for examining subduction because of its heterogenous natural structure. In this study, the Visibility Graph method was applied to study the seismic activity of three regions within the Cocos plate: the Tehuantepec Isthmus, the Flat slab and Michoacan, each one with a different level of seismicity. The method maps time series into graphs, and it is possible to connect the topological properties of the graph with the dynamical features underlying the time series. The seismicity analyzed was monitored in the three areas studied between 2010 and 2022. At the Flat Slab and Tehuantepec Isthmus, two intense earthquakes occurred on 7 and 19 September 2017, respectively, and, on 19 September 2022, another one occurred at Michoacan. The aim of this study was to determine the dynamical features and the possible differences between the three areas by applying the following method. First, the time evolution of the a- and b-values in the Gutenberg-Richter law was analyzed, followed by the relationship between the seismic properties and topological features using the VG method, the k-M slope and the characterization of the temporal correlations from the γ-exponent of the power law distribution, P(k) ∼ k-γ, and its relationship with the Hurst parameter, which allowed us to identify the correlation and persistence of each zone.

Why supercontinents became shorter lived as the Earth evolved

Periodic assembly and break-up of supercontinents since at least two billion years ago (Ga), dubbed the supercontinent cycle, provides the first-order tectonic control on the evolution of the Earth System including episodic orogenic events, mineralization, the formation and closure of oceans and superoceans, and even the evolution of life. However, the lifespan of the supercontinents appears to decrease with time, from ∼300 million years (Myr) for Nuna/Columbia, to 200-250 Myr for Rodinia and ∼150 Myr for the youngest supercontinent Pangaea. To understand what caused such a secular decrease in supercontinental lifespan, we conduct 3-D geodynamic modeling using realistic tectonic settings. The results show that the yield stress of newly formed orogens during the assembly of a supercontinent provides the dominant control on the lifespan of the supercontinent, implying that the yield stress of young orogens becomes lower with time. We hypothesize that the decreasing mantle temperature due to Earth's secular cooling might have caused new orogens to become weaker.

Back-Arc Spreading Centers and Superfast Subduction: The Case of the Northern Lau Basin (SW Pacific Ocean)

The Lau Basin is a back-arc region formed by the subduction of the Pacific plate below the Australian plate. We studied the regional morphology of the back-arc spreading centers of the Northern Lau basin, and we compared it to their relative spreading rates. We obtained a value of 60.2 mm/year along the Northwest Lau Spreading Centers based on magnetic data, improving on the spreading rate literature data. Furthermore, we carried out numerical models including visco-plastic rheologies and prescribed surface velocities, in an upper plate-fixed reference frame. Although our thermal model points to a high temperature only near the Tonga trench, the model of the second invariant of the strain rate shows active deformation in the mantle from the Tonga trench to ~800 km along the overriding plate. This explains the anomalous magmatic production along all the volcanic centers in the Northern Lau Back-Arc Basin.

The Deep-Time Digital Earth program: data-driven discovery in geosciences

Current barriers hindering data-driven discoveries in deep-time Earth (DE) include: substantial volumes of DE data are not digitized; many DE databases do not adhere to FAIR (findable, accessible, interoperable and reusable) principles; we lack a systematic knowledge graph for DE; existing DE databases are geographically heterogeneous; a significant fraction of DE data is not in open-access formats; tailored tools are needed. These challenges motivate the Deep-Time Digital Earth (DDE) program initiated by the International Union of Geological Sciences and developed in cooperation with national geological surveys, professional associations, academic institutions and scientists around the world. DDE’s mission is to build on previous research to develop a systematic DE knowledge graph, a FAIR data infrastructure that links existing databases and makes dark data visible, and tailored tools for DE data, which are universally accessible. DDE aims to harmonize DE data, share global geoscience knowledge and facilitate data-driven discovery in the understanding of Earth's evolution.

Explosive fragmentation of Prince Rupert's drops leads to well-defined fragment sizes

Anyone who has ever broken a dish or a glass knows that the resulting fragments range from roughly the size of the object all the way down to indiscernibly small pieces: typical fragment size distributions of broken brittle materials follow a power law, and therefore lack a characteristic length scale. The origin of this power-law behavior is still unclear, especially why it is such an universal feature. Here we study the explosive fragmentation of glass Prince Rupert’s drops, and uncover a fundamentally different breakup mechanism. The Prince Rupert’s drops explode due to their large internal stresses resulting in an exponential fragment size distribution with a well-defined fragment size. We demonstrate that generically two distinct breakup processes exist, random and hierarchical, that allows us to fully explain why fragment size distributions are power-law in most cases but exponential in others. We show experimentally that one can even break the same material in different ways to obtain either random or hierarchical breakup, giving exponential and power-law distributed fragment sizes respectively. That a random breakup process leads to well-defined fragment sizes is surprising and is potentially useful to control fragmentation of brittle solids.

Fragmentation of breaking glass as a brittle solid is a problem of equal practical and theoretical importance. Kooij et al. demonstrate that the fragment size distribution can surprisingly be both, either power-law or exponential, depending on how a particular specimen is broken.

Plato's cube and the natural geometry of fragmentation

Plato envisioned Earth's building blocks as cubes, a shape rarely found in nature. The solar system is littered, however, with distorted polyhedra-shards of rock and ice produced by ubiquitous fragmentation. We apply the theory of convex mosaics to show that the average geometry of natural two-dimensional (2D) fragments, from mud cracks to Earth's tectonic plates, has two attractors: "Platonic" quadrangles and "Voronoi" hexagons. In three dimensions (3D), the Platonic attractor is dominant: Remarkably, the average shape of natural rock fragments is cuboid. When viewed through the lens of convex mosaics, natural fragments are indeed geometric shadows of Plato's forms. Simulations show that generic binary breakup drives all mosaics toward the Platonic attractor, explaining the ubiquity of cuboid averages. Deviations from binary fracture produce more exotic patterns that are genetically linked to the formative stress field. We compute the universal pattern generator establishing this link, for 2D and 3D fragmentation.

What drives tectonic plates?

Does Earth's mantle drive plates, or do plates drive mantle flow? This long-standing question may be ill posed, however, as both the lithosphere and mantle belong to a single self-organizing system. Alternatively, this question is better recast as follows: Does the dynamic balance between plates and mantle change over long-term tectonic reorganizations, and at what spatial wavelengths are those processes operating? A hurdle in answering this question is in designing dynamic models of mantle convection with realistic tectonic behavior evolving over supercontinent cycles. By devising these models, we find that slabs pull plates at rapid rates and tear continents apart, with keels of continents only slowing down their drift when they are not attached to a subducting plate. Our models show that the tectonic tessellation varies at a higher degree than mantle flow, which partly unlocks the conceptualization of plate tectonics and mantle convection as a unique, self-consistent system.

Crustal evolution and mantle dynamics through Earth history

Resolving the modes of mantle convection through Earth history, i.e. when plate tectonics started and what kind of mantle dynamics reigned before, is essential to the understanding of the evolution of the whole Earth system, because plate tectonics influences almost all aspects of modern geological processes. This is a challenging problem because plate tectonics continuously rejuvenates Earth's surface on a time scale of about 100 Myr, destroying evidence for its past operation. It thus becomes essential to exploit indirect evidence preserved in the buoyant continental crust, part of which has survived over billions of years. This contribution starts with an in-depth review of existing models for continental growth. Growth models proposed so far can be categorized into three types: crust-based, mantle-based and other less direct inferences, and the first two types are particularly important as their difference reflects the extent of crustal recycling, which can be related to subduction. Then, a theoretical basis for a change in the mode of mantle convection in the Precambrian is reviewed, along with a critical appraisal of some popular notions for early Earth dynamics. By combining available geological and geochemical observations with geodynamical considerations, a tentative hypothesis is presented for the evolution of mantle dynamics and its relation to surface environment; the early onset of plate tectonics and gradual mantle hydration are responsible not only for the formation of continental crust but also for its preservation as well as its emergence above sea level. Our current understanding of various material properties and elementary processes is still too premature to build a testable, quantitative model for this hypothesis, but such modelling efforts could potentially transform the nature of the data-starved early Earth research by quantifying the extent of preservation bias.This article is part of a discussion meeting issue 'Earth dynamics and the development of plate tectonics'.