Significance This study tackles two long-standing questions: 1) the morphology and physical state of the subducting Indian slab under Tibet, along with their impact on growth dynamics of the plateau, and 2) the driving forces behind the ongoing northward motion of the Indian Plate during its postcollision with Asia. We present evidence that the Indian slab subducts beneath almost the entire plateau, with a subhorizontal geometry at ∼300 km depth. The mechanism by which the plateau maintains its high-elevation, low-relief topography is potentially attributed to buoyancy of the low-velocity zone above the flat slab. Mantle drag, exerted on the base of the Indian Plate by plume-associated northward mantle convective flows, likely serves as the primary force driving the India–Asia collision.

Abstract Mantle convection models based on geophysical constraints have provided us with a basic understanding of the forces driving and resisting plate motions on Earth. However, existing studies computing the balance of underlying forces are contradicting, and the impact of plate boundary geometry on surface deformation remains unknown. We address these issues by developing global instantaneous 3‐D mantle convection models with a heterogeneous density and viscosity distribution and weak plate boundaries prescribed using different geometries. We find that the plate boundary geometry of the Global Earthquake Model (GEM, Pagani et al., 2018, https://doi.org/10.1177/8755293020931866 ), featuring open plate boundaries with discrete lithospheric‐depth weak zones in the oceans and distributed crustal faults within continents, achieves the best fit to the observed GPS data with a directional correlation of 95.1% and a global point‐wise velocity residual of 1.87 cm/year. A good fit also requires plate boundaries being 3 to 4 orders of magnitude weaker than the surrounding lithosphere and low asthenospheric viscosities between 5 × 10 17 and 5 × 10 18  Pa s. Models without asthenospheric and lower mantle heterogeneities retain on average 30% and 70% of the plate speeds, respectively. Our results show that Earth's plate boundaries are not uniform and better described by more discrete plate boundaries within the oceans and distributed faults within continents. Furthermore, they emphasize the impact of plate boundary geometry on the direction and speed of plate motions and reaffirm the importance of slab pull in the uppermost mantle as a major plate driving force.

Earth's tectonic plates are generally considered to be driven largely by negative buoyancy associated with subduction of oceanic lithosphere. In this context, mid-ocean ridges (MORs) are passive plate boundaries whose divergence accommodates flow driven by subduction of oceanic slabs at trenches. We show that over the past 80 million years (My), the East Pacific Rise (EPR), Earth's dominant MOR, has been characterized by limited ridge-perpendicular migration and persistent, asymmetric ridge accretion that are anomalous relative to other MORs. We reconstruct the subduction-related buoyancy fluxes of plates on either side of the EPR. The general expectation is that greater slab pull should correlate with faster plate motion and faster spreading at the EPR. Moreover, asymmetry in slab pull on either side of the EPR should correlate with either ridge migration or enhanced plate velocity in the direction of greater slab pull. Based on our analysis, none of the expected correlations are evident. This implies that other forces significantly contribute to EPR behavior. We explain these observations using mantle flow calculations based on globally integrated buoyancy distributions that require core-mantle boundary heat flux of up to 20 TW. The time-dependent mantle flow predictions yield a long-lived deep-seated upwelling that has its highest radial velocity under the EPR and is inferred to control its observed kinematics. The mantle-wide upwelling beneath the EPR drives horizontal components of asthenospheric flows beneath the plates that are similarly asymmetric but faster than the overlying surface plates, thereby contributing to plate motions through viscous tractions in the Pacific region.

The development of plate tectonics from a pre-plate tectonics regime requires both the initiation of subduction and the development of nascent subduction zones into long-lived contiguous features. Subduction itself has been shown to be sensitive to system parameters such as thermal state and the specific rheology. While generally it has been shown that cold-interior high-Rayleigh-number convection (such as on the Earth today) favours plates and subduction, due to the ability of the interior stresses to couple with the lid, a given system may or may not have plate tectonics depending on its initial conditions. This has led to the idea that there is a strong history dependence to tectonic evolution-and the details of tectonic transitions, including whether they even occur, may depend on the early history of a planet. However, intrinsic convective stresses are not the only dynamic drivers of early planetary evolution. Early planetary geological evolution is dominated by volcanic processes and impacting. These have rarely been considered in thermal evolution models. Recent models exploring the details of plate tectonic initiation have explored the effect of strong thermal plumes or large impacts on surface tectonism, and found that these 'primary drivers' can initiate subduction, and, in some cases, over-ride the initial state of the planet. The corollary of this, of course, is that, in the absence of such ongoing drivers, existing or incipient subduction systems under early Earth conditions might fail. The only detailed planetary record we have of this development comes from Earth, and is restricted by the limited geological record of its earliest history. Many recent estimates have suggested an origin of plate tectonics at approximately 3.0 Ga, inferring a monotonically increasing transition from pre-plates, through subduction initiation, to continuous subduction and a modern plate tectonic regime around that time. However, both numerical modelling and the geological record itself suggest a strong nonlinearity in the dynamics of the transition, and it has been noted that the early history of Archaean greenstone belts and trondhjemite-tonalite-granodiorite record many instances of failed subduction. Here, we explore the history of subduction failure on the early Earth, and couple these with insights from numerical models of the geodynamic regime at the time.This article is part of a discussion meeting issue 'Earth dynamics and the development of plate tectonics'.

Sub-plate mantle flow traction (MFT) has been considered as a major driving force for plate motion; however, the force acting on the overlying plate is difficult to constrain. One of the reasons lies in the variable rheological flow laws of mantle rocks, e.g. linear versus power-law rheology, applied in previous studies. Here, systematic numerical models are conducted to evaluate MFT under variable rheological, geometrical and kinematic conditions. The results indicate that MFT with power-law rheology is much lower than that with linear rheology under the same mantle/plate velocity contrast. In addition, existence of a lithospheric root in the overlying plate could enhance MFT, where integrated normal force acting on the walls of lithospheric root is much lower than the shear force in a large-scale domain. In the acting domain of several thousand kilometers, MFT with power-law rheology is comparable to the ridge push of about 3×1012 N/m, whereas that with linear rheology is comparable to the slab pull of about 3×1013 N/m. The roles of MFT in driving plate motion are further analyzed for the Tethyan evolution. It indicates that MFT with power-law rheology could partially support the Wilson cycles experienced in the Tethyan system, whereas that with linear rheology could easily dominate any kinds of plate tectonic evolutions. The quantitative evaluation of MFT in this study clarifies the roles of rheological flow laws on MFT and could help to better understand the contrasting results in previous numerical studies.

Plate tectonics is primarily driven by the constant gravitational pull of slabs where oceanic lithosphere sinks into the mantle at subduction zones. Under stable plate boundary configurations, changes in plate motion are then thought to occur gradually. Surprisingly, recent high-resolution Indian plate reconstructions revealed rapid (2–3 Million-year) plate velocity oscillations of ±50 %. We show, through numerical experiments, that the buckling of slabs in the mantle transition zone causes such oscillations. This buckling results from the deceleration of slabs as they sink into the lower mantle. The amplitude and period of buckling-associated oscillations depend on average subduction velocity and the available space in the mantle transition zone. The oscillations also affect the upper plate which may explain enigmatic observations of episodic deformation and fluid flow in subduction-related orogens. We infer that the slab pull that drives plate tectonics is generated in just the top few hundred kilometers of the mantle. Slab buckling in the mantle transition zone can explain 2-3 million-year oscillations in the Indian plate velocity during rapid subduction, based on two-dimensional numerical models.

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 (10¹³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.

In his paper “Über das Bewegungsbild von Faltengebirgen” [On the movement pattern of folded mountains], published in the almanac of the Austrian Geological Survey in Vienna, Otto Ampferer from Innsbruck (Austria) presented a series of geotectonic considerations and interpretations, which today are summarized under the term “theory of undercurrent”. The interpretation of these processes occurring in the deep crust of the Earth and in the upper mantle was mainly kinematic. For a long time, the tectonic passivity of the magma being anorogenic was dogma until Ampferer’s undercurrent theory changed this in 1906, according to which folds and thrusts on the Earth’s surface portray motions of the deeper magmatic substratum. In these undercurrents, Ampferer recognized the crucial forces that lead to the formation of ocean basins and high mountains on the edges of the drifting continents. In his paper on the history of the formation of the Atlantic region, he presented already, in 1941, a process anticipating what is now known as seafloor spreading.

Forces associated with subduction are an integral aspect of plate tectonics and mantle processes. However due to the non-linear overlap of physical processes and forces driving subduction their tectonic implications remain difficult to resolve. The aim of this thesis is to examine insufficiently understood drivers of subduction and relate them to their explicit tectonic expressions at regional and global scales. To systematically study the parameters controlling subduction this investigation utilises the Stokes flow solver, BEM-Earth, and the global mantle convection code, Terra. I investigate the role of the plates overriding the subduction zone and the feedbacks between the down-going slab interacting with the surface plates. Particularly, I distinguish the force partitioning between the trench-orthogonal and trench-perpendicular mantle flow induced by a subducting slab’s overriding plate. It is determined that overriding plate rheology, boundary configuration, and kinematics interact and partially control subducting slab morphology and velocity. This is important for explaining surface plate motions and apparent structure of positive wave-speed anomalies in the mantle imaged with seismic tomography. I further analyse subduction driven forces acting on the Pacific Plate, during the Late Cretaceous and early Cenozoic, and spatially and temporally correlate resulting plate motions and subsequent lithospheric deformation to geophysical observations. The modeled buoyancy driven

Plate tectonics can explain several geological and geophysical phenomena on Earth, and a number of mantle flow models have been developed to investigate the underlying plate tectonic forces. However, these models have come to contradictory conclusions on the balance between the resisting and driving forces. Additionally, they have used the same simplified model to represent the geometry of the plates, and therefore the impact of plate boundary geometry on surface deformation remains unknown.To address these issues, we have developed high-resolution global instantaneous mantle convection models based on recent geophysical constraints with a heterogeneous density and viscosity distribution and weak plate boundaries prescribed using different plate boundary configurations. We find a good fit to the observed GPS data for models with plate boundaries that are 3 to 4 orders of magnitude weaker than the surrounding lithosphere and low asthenospheric viscosities between 5×1017 and 5×1018 Pa s for all plate boundary configurations. We also find that the model with plate boundaries defined by the Global Earthquake Model (GEM, Pagani et al., 2018)—featuring open plate boundaries with discrete lithospheric-depth weak zones in the oceans and distributed crustal faults within continents—achieves the best fit to the observed GPS data with a directional correlation of 95.1% and a global point-wise velocity residual of 1.87 cm/year. These results show that Earth’s plate boundaries are not uniform and better described by more discrete plate boundaries within the oceans and distributed faults within continents.Our models also quantify the contributions to the plate driving forces originating from heterogeneities in the upper mantle and the lower mantle, respectively, finding that the slab-pull in the top 300 km alone contributes ~70% of the total plate speeds. Noting the importance of slab pull as a major plate driving force, we further investigate the influence of subduction zone and slab geometry on surface plate motions and their fit to GPS data. Specifically, our models compare a simplified slab structure to a more detailed representation of slabs based on the Slab2 database (Hayes et al., 2018), and reaffirm that a realistic slab geometry is a crucial factor in the transmission of slab pull forces to the plate.

sink into the mantle can cause viscous mantle forces driving plates through slab suction. In the theory of plume tectonics followed by numerous researchers Plate tectonics (from Latin tectonicus, from Ancient Greek τεκτονικός (tektonikós) 'pertaining to building') is the scientific theory that Earth's lithosphere comprises a number of large tectonic plates, which have been slowly moving since 3–4 billion years ago. The model builds on the concept of continental drift, an idea developed during the first decades of the 20th century. Plate tectonics c For much of the first quarter of the 20th century, the leading theory of the driving force behind tectonic plate motions envisaged large scale convection currents in the upper mantle, which can be transmitted through the asthenosphere. This theory was launched by Arthur Holmes and some forerunners in the 1930s and was immediately recognized as the solution for the acceptance of the theory as originally discussed in the papers of Alfred Wegener in the early years of the 20th century. However, despite its acceptance, it was long debated in the scientific community because the leading theory still envisaged a static Earth without moving continents up until the major breakthroughs of the early sixties. Two- and three-dimensional imaging of Earth's interior (seismic tomography) shows a varying lateral density distribution throughout the mantle. Such density variations can be material (from rock chemistry), mineral (from variations in mineral structures), or thermal (through thermal expansion and contraction from heat energy). The manifestation of this varying lateral density is mantle convection from buoyancy forces. How mantle convection directly and indirectly relates to plate motion is a matter of ongoing study and discussion in geodynamics. Somehow, this energy must be transferred to the lithosphere for tectonic plates to move. There are essentially two main types of mechanisms that are thought to exist related to the dynamics of the mantle that influence plate motion which are primary (through the large scale convection cells) or secondary. The secondary mechanisms view plate motion driven by friction between the convection currents in the asthenosphere and the more rigid overlying lithosphere. This is due to the inflow of mantle material related to the downward pull on plates in subduction zones at ocean trenches. Slab pull may occur in a geodynamic setting where basal tractions continue to act on the plate as it dives into the mantle (although perhaps to a greater extent acting on both the under and upper side of the slab). Furthermore, slabs that are broken off and sink into the mantle can cause viscous mantle forces driving plates through slab suction. In…

Abstract Large shallow intraslab earthquakes (LSIEs) of over magnitude-7 rarely occur at shallower depths (≤ 60 km) in the subducting plates of the circum-Pacific and northern margin of the Australian Plate off Indonesia. Previous studies have suggested that most LSIEs occur under surface tectonic conditions with a lateral stress gradient across the back-arc to the fore-arc on the overriding plate based on seismological and geological evidence. In this study, dynamic conditions for the occurrence of LSIEs were studied using the intraslab stress state, stress state of the overriding plate, motion speed of overriding and subducting plates, and trench migration speed determined by plate motion models. LSIEs prefer tectonic conditions under (1) an extensional-to-neutral stress state in the back-arc, (2) a compressional stress state in the fore-arc, and (3) down-dip tension in the shallower part of the subducting plate. The results suggest that the Earth’s subduction zones can be classified into four categories. This categorization is informed by the behavior of underlying mantle flow, i.e., the magnitude of slab suction flow in the mantle under the subducting plates, which is related to the degree of plate-slab coupling, and the scale of return flow in the mantle wedge under the overriding plates, which varies the combination of stress states in the fore-arc and back-arc.

Abstract The paper considers the history of the spreading of the Eurasian basin. The sharp deceleration of the spreading rate in the Eocene about 46 million years ago, which is fixed by the distribution of linear magnetic anomalies, is noted. That jump in velocity is clarified from the perspective of the geodynamic model but shouldnt be explained by the northern motion of Greenland. The geodynamic processes of the Pacific subduction zone generate an upper mantle convective cell with return flow dragging the Arctic continental lithosphere in the direction of the Pacific subduction zone. The geodynamic mechanism is confirmed by seismic tomographic mantle sections of the northeastern margin of Asia and the numerical model of the upper mantle convection of the active continental margin. It is the activity of the upper mantle convective return cell, which is determined by the runoff volume and, ultimately, the speed and direction of the Kula plate and Pacific plate subduction vectors in the subduction zone, affects tectonics and kinematics of the plates of the Eurasian basin. In the Middle CretaceousMiddle Eocene and for about 73 Ma the return cell has been active, since the Kula and Pacific plates move north and submerged orthogonally beneath the Central Arctic. After the Middle Eocene geodynamic reorganization about 47.5 million years ago, oceanic plates in the Pacific Ocean begin to move to the northwest. As a result, the transport of the oceanic Pacific Ocean lithospheric subs

crust is formed through volcanic activity and then gradually moves away from the ridge Slab pull – That part of the motion of a tectonic plate that is caused The following outline is provided as an overview of and introduction to Oceanography. Oceanography (from Ancient Greek ὠκεανός (ōkeanós) 'ocean' and γραφή (graphḗ) 'writing'), also known as oceanology, sea science, ocean science, and marine science, is the scientific study of the ocean, including its physics, chemistry, biology, and geology. It is an Earth science, which covers a wide range o Ast…

Other volcanoes have unusual shapes, such as the “pancake domes” illustrated in Figure 10.8. All of the volcanism is the result of eruption of lava onto the surface of the planet. But the hot lava rising from the interior of a planet does not always make it to the surface. On both Earth and Venus, this upwelling lava can collect to produce bulges in the crust. Many of the granite mountain ranges on Earth, such as the Sierra Nevada in California, involve such subsurface volcanism. These bulges are common on Venus, where they produce large circular or oval features called coronae (singular: corona) (Figure 10.9). Convection currents of molten material in the mantle of Venus push and stretch the crust. Such forces are called tectonic, and the geological features that result from these forces are called tectonic features. On Venus’ lowland plains, tectonic forces have broken the lava surface to create remarkable patterns of ridges and cracks (Figure 10.10). In a few places, the crust has even torn apart to generate rift valleys. The circular features associated with coronae are tectonic ridges and cracks, and most of the mountains of Venus also owe their existence to tectonic forces.

Flat subduction, an intermittent phenomenon along active margins, arises from well-known causes, yet the mechanisms driving its expansion remain poorly understood. The prevailing view suggests that trenchward continental motion drives slab overthrusting, causing the flat slab to expand oceanward. Here, we explore an alternative mechanism: underthrusting of the subducting plate through forward propagation of the flat-slab hinge. We directly evaluate both hypotheses through a kinematic analysis of trench and flat-slab motions using a global flat subduction database cast into multiple absolute plate motion models. Our results indicate that flat-slab expansion reflects distinct end-member processes, with forward propagation emerging as the dominant mode. We present a framework for flat-slab propagation that emphasizes the dynamic interaction between lower-plate motion and slab pull from adjacent subduction zones, an interplay validated through numerical modeling. These findings challenge conventional assumptions and underscore the need to reconsider the role of lower-plate kinematics in flat-slab dynamics.