scholarly journals Effects of basal drag on subduction dynamics from 2D numerical models

2020 ◽  
Author(s):  
Lior Suchoy ◽  
Saskia Goes ◽  
Benjamin Maunder ◽  
Fanny Garel ◽  
Rhodri Davies
Keyword(s):  
Numerical Models ◽  
Force Balance ◽  
Plate Motion ◽  
Mantle Flow ◽  
Plate Motions ◽  
Flow Models ◽  
Plate Size ◽  
Modelling Studies ◽  
Subduction System

Abstract. Subducting slabs are an important driver of plate motions, yet the force balance governing subduction dynamics remains incompletely understood. Basal drag has been proposed to be a minor contributor to subduction forcing, because of the lack of correlation between plate size and velocity in observed and reconstructed plate motions. Furthermore, in single subduction system models, low basal drag, associated with a low ratio of asthenospheric to lithospheric viscosity, leads to subduction behaviour most consistent with the observation that trench migration velocities are generally low compared to convergence velocities. By contrast, analytical calculations and global mantle flow models indicate basal drag can be substantial. In this study, we revisit this problem by examining the drag at the base of the lithosphere, for a single subduction system, in 2D models with a free trench and composite non-linear rheology. We compare the behaviour of short and long plates for a range of asthenospheric and lithospheric rheologies. We reproduce results from previous modelling studies, including low ratios of trench over plate motions. However, we also find that any combination of asthenosphere and lithosphere viscosity that produces Earth-like subduction behaviour leads to a correlation of velocities with plate size, due to the role of basal drag. By examining Cenozoic plate motion reconstructions, we find that slab age and plate size are positively correlated: higher slab pull for older plates tends to be offset by higher basal drag below these larger plates. This, in part, explains the lack of plate velocity-size correlation in observations, despite the important role of basal drag in the subduction force-balance.

scholarly journals Effects of basal drag on subduction dynamics from 2D numerical models

Solid Earth ◽  
2021 ◽  
Vol 12 (1) ◽  
pp. 79-93
Author(s):  
Lior Suchoy ◽  
Saskia Goes ◽  
Benjamin Maunder ◽  
Fanny Garel ◽  
Rhodri Davies
Keyword(s):  
Numerical Models ◽  
Force Balance ◽  
Plate Motion ◽  
Mantle Flow ◽  
Plate Motions ◽  
Flow Models ◽  
Plate Size ◽  
Modelling Studies ◽  
Subduction System

Abstract. Subducting slabs are an important driver of plate motions, yet the relative importance of different forces in governing subduction motions and styles remains incompletely understood. Basal drag has been proposed to be a minor contributor to subduction forcing because of the lack of correlation between plate size and velocity in observed and reconstructed plate motions. Furthermore, in single subduction system models, low basal drag leads to subduction behaviour most consistent with the observation that trench migration velocities are generally low compared to convergence velocities. By contrast, analytical calculations and global mantle flow models indicate basal drag can be substantial. In this study, we revisit this problem by examining the drag at the base of the lithosphere, for a single subduction system, in 2D models with a free trench and composite non-linear rheology. We compare the behaviour of short and long plates for a range of asthenospheric and lithospheric rheologies. We reproduce results from previous modelling studies, including low ratios of trench over plate motions. However, we also find that any combination of asthenosphere and lithosphere viscosity that produces Earth-like subduction behaviour leads to a correlation of velocities with plate size, due to the role of basal drag. By examining Cenozoic plate motion reconstructions, we find that slab age and plate size are positively correlated: higher slab pull for older plates tends to be offset by higher basal drag below these larger plates. This, in part, explains the lack of plate velocity–size correlation in observations, despite the important role of basal drag in the subduction force balance.

Coupled evolution of supercontinents and lower mantle structure

2021 ◽  
Author(s):  
Xianzhi Cao ◽  
Nicolas Flament ◽  
Ömer Bodur ◽  
Dietmar Müller
Keyword(s):  
Lower Mantle ◽  
Reference Frames ◽  
Plate Motion ◽  
Motion Path ◽  
Mantle Flow ◽  
Plate Tectonic ◽  
Mantle Structure ◽  
Plate Motions ◽  
Flow Models ◽  
Over Time

<p>The relationships between plate motions and basal mantle structure remain poorly understood, with some models implying that the basal mantle structure has remained stable over time, while others suggest that it could be shaped by the aggregation and dispersal of supercontinents. Here we investigate the plate-basal mantle relationship through 1) building a series of end-member plate tectonic models over one billion years, and 2) creating mantle flow models assimilated by those plate models. To achieve that, we build synthetic plate tectonic models dating from 1 Ga to 250 Ma that we connect to an existing palaeogeographical plate reconstruction from 250 Ma to create a relative plate motion model for the last 1 Gyr, in which supercontinent breakup and reassembly occur via introversion. We consider three distinct reference frames that result in different net lithospheric rotation. We find that the flow models predict a dominant degree-2 lower mantle structure most of the time and that they are in first-order agreement (~70% spatial match) with tomographic models. Model thermochemical structures at the base of the mantle may split into smaller structures when slabs sink onto them, and smaller basal structures may merge into larger ones as a result of slab pushing. The basal thermochemical structure under the superocean is large and continuous, whereas the basal thermochemical structure under the supercontinent is smaller and progressively assembles during and shortly after supercontinent assembly. In the models, plumes also develop preferentially along the edge of the basal thermochemical structures and tend to migrate towards the interior of basal structures over time as they interact with the slabs. Lone plumes can also form away from the main thermochemical structures, often within a small network of sinking slabs. Lone plumes may migrate between basal structures. We analyse the relationship between imposed tectonic velocities and deep mantle flow, and find that at spherical harmonic degree 2, the maxima of lower mantle radial flow and temperature follow the motion path of the maxima of surface divergence. It may take ~160-240 Myr for lower mantle structure to reflect plate motion changes when the lower mantle is reorganised by slabs sinking onto basal thermochemical structures, and/or when slabs stagnate in the transition zone before sinking to the lower mantle. Basal thermochemical structures move at less than 0.6 °/Myr in our models with a temporal average of 0.16 °/Myr when there is no net lithospheric rotation, and between 0.20-0.23 °/Myr when net lithospheric rotation exists and is induced to the lower mantle. Our results suggest that basal thermochemical structures are not stationary, but rather linked to global plate motions and plate boundary reconfigurations, reflecting the dynamic nature of the co-evolving plate-mantle system.</p>

scholarly journals Numerical models of slab migration in continental collision zones

Solid Earth ◽  
2012 ◽  
Vol 3 (2) ◽  
pp. 293-306 ◽  
Author(s):  
V. Magni ◽  
J. van Hunen ◽  
F. Funiciello ◽  
C. Faccenna
Keyword(s):  
Subduction Zone ◽  
Plate Tectonics ◽  
Numerical Models ◽  
Continental Collision ◽  
Force Balance ◽  
Small Scale ◽  
Dip Angle ◽  
External Forces ◽  
Stress Dependent ◽  
Subduction System

Abstract. Continental collision is an intrinsic feature of plate tectonics. The closure of an oceanic basin leads to the onset of subduction of buoyant continental material, which slows down and eventually stops the subduction process. In natural cases, evidence of advancing margins has been recognized in continental collision zones such as India-Eurasia and Arabia-Eurasia. We perform a parametric study of the geometrical and rheological influence on subduction dynamics during the subduction of continental lithosphere. In our 2-D numerical models of a free subduction system with temperature and stress-dependent rheology, the trench and the overriding plate move self-consistently as a function of the dynamics of the system (i.e. no external forces are imposed). This setup enables to study how continental subduction influences the trench migration. We found that in all models the slab starts to advance once the continent enters the subduction zone and continues to migrate until few million years after the ultimate slab detachment. Our results support the idea that the advancing mode is favoured and, in part, provided by the intrinsic force balance of continental collision. We suggest that the advance is first induced by the locking of the subduction zone and the subsequent steepening of the slab, and next by the sinking of the deepest oceanic part of the slab, during stretching and break-off of the slab. These processes are responsible for the migration of the subduction zone by triggering small-scale convection cells in the mantle that, in turn, drag the plates. The amount of advance ranges from 40 to 220 km and depends on the dip angle of the slab before the onset of collision.

An analytical model concerning the possible initiation of subduction between the India and Africa plates, caused by the Morondova plume

2020 ◽  
Author(s):  
Bernhard Steinberger ◽  
Douwe van Hinsbergen
Keyword(s):  
Subduction Zones ◽  
Plate Tectonics ◽  
Eastern Mediterranean ◽  
Plate Boundary ◽  
Plate Motion ◽  
Mantle Flow ◽  
Plate Motions ◽  
Subduction Initiation ◽  
Euler Pole ◽  
Geodynamic Processes

<p>Identifying the geodynamic processes that trigger the formation of new subduction zones is key to understand what keeps the plate tectonic cycle going, and how plate tectonics once started. Here we discuss the possibility of plume-induced subduction initiation. Previously, our numerical modeling revealed that mantle upwelling and radial push induced by plume rise may trigger plate motion change, and plate divergence as much as 15-20 My prior to LIP eruption. Here we show that, depending on the geometry of plates, the distribution of cratonic keels and where the plume rises, it may also cause a plate rotation around a pole that is located close to the same plate boundary where the plume head impinges: If that occurs near one end of the plate boundary, an Euler pole of the rotation may form along that plate boundary, with extension on one side, and convergence on the other.  This concept is applied to the India-Africa plate boundary and the Morondova plume, which erupted around 90 Ma, but may have influenced plate motions as early as 105-110 Ma. If there is negligible friction, i.e. there is a pre-existing weak plate boundary, we estimate that the total amount of convergence generated in the northern part of the India-Africa plate boundary can exceed 100 km, which is widely thought to be sufficient to initiate forced, self-sustaining subduction. This may especially occur if the India continental craton acts like an “anchor” causing a comparatively southern location of the rotation pole of the India plate. Geology and paleomagnetism-based reconstructions of subduction initiation below ophiolites from Pakistan, through Oman, to the eastern Mediterranean reveal that E-W convergence around 105 Ma caused forced subduction initiation, and we tentatively postulate that this is triggered by Morondova plume head rise. Whether the timing of this convergence is appropriate to match observations on subduction initiation as early as 105 Ma depends on the timing of plume head arrival, which may predate eruption of the earliest volcanics. It also depends on whether a plume head already can exert substantial torque on the plate while it is still rising – for example, if the plate is coupled to the induced mantle flow by a thick craton.</p>

scholarly journals The role of the margins in the dynamics of an active ice stream

1994 ◽  
Vol 40 (136) ◽  
pp. 527-538 ◽  
Author(s):  
K. A. Echelmeyer ◽  
W. D. Harrison ◽  
C. Larsen ◽  
J. E. Mitchell
Keyword(s):  
Numerical Models ◽  
Force Balance ◽  
Shear Stresses ◽  
List Type ◽  
Ice Stream ◽  
Strain Heating ◽  
Ice Stream B ◽  
Basal Shear ◽  
Marginal Zones

AbstractA transverse profile of velocity was measured across Ice Stream B, West Antarctica, in order to determine the role of the margins in the force balance of an active ice stream. The profile extended from near the ice-stream center line, through a marginal shear zone and on to the slow-moving ice sheet. The velocity profile exhibits a high degree of shear deformation within a marginal zone, where intense, chaotic crevassing occurs. Detailed analysis of the profile, using analytical and numerical models of ice flow, leads to the following conclusions regarding the roles of the bed and the margins in ice-stream dynamics:(i)The overall resistive drag on the ice stream is partitioned nearly equally between the margins and the bed and, thus, both are important in the force balance of the ice stream.(ii)The ice within the chaotic zone must be about 10 times softer than the ice in the central part of the ice stream.(iii)The average basal shear stress is 0.06 × 105Pa. This implies that the entire bed cannot be blanketed by the weak, deformable till observed by Engelhardt and others (1990) near the center of the ice stream — there must be regions of increased basal drag.(iv)High strain rates and shear stresses in the marginal zones indicate that strain heating in the margins may be significant.While the exact quantitative values leading to these conclusions are somewhat model and location-dependent, the overall conclusions are robust. As such, they are likely to have importance for ice-stream dynamics in general.

scholarly journals Numerical models of trench migration in continental collision zones

2012 ◽  
Vol 4 (1) ◽  
pp. 429-458 ◽  
Author(s):  
V. Magni ◽  
J. van Hunen ◽  
F. Funiciello ◽  
C. Faccenna
Keyword(s):  
Subduction Zone ◽  
Plate Tectonics ◽  
Numerical Models ◽  
Continental Collision ◽  
Force Balance ◽  
Dip Angle ◽  
External Forces ◽  
Stress Dependent ◽  
Subduction System ◽  
Intrinsic Feature

Abstract. Continental collision is an intrinsic feature of plate tectonics. The closure of an oceanic basin leads to the onset of subduction of buoyant continental material, which slows down and eventually stops the subduction process. We perform a parametric study of the geometrical and rheological influence on subduction dynamics during the subduction of continental lithosphere. In 2-D numerical models of a free subduction system with temperature and stress-dependent rheology, the trench and the overriding plate move self-consistently as a function of the dynamics of the system (i.e. no external forces are imposed). This setup enables to study how continental subduction influences the trench migration. We found that in all models the trench starts to advance once the continent enters the subduction zone and continues to migrate until few million years after the ultimate slab detachment. Our results support the idea that the trench advancing is favoured and, in part provided by, the intrinsic force balance of continental collision. We suggest that the trench advance is first induced by the locking of the subduction zone and the subsequent steepening of the slab, and next by the sinking of the deepest oceanic part of the slab, during stretching and break-off of the slab. The amount of trench advancing ranges from 40 to 220 km and depends on the dip angle of the slab before the onset of collision.

Subduction zone seismo-dynamics: how to bridge the gap between long-term subduction dynamics and megathrust seismicity?

2021 ◽  
Author(s):  
Adam Beall ◽  
Fabio A. Capitanio ◽  
Ake Fagereng ◽  
Ylona van Dinther
Keyword(s):  
Stress State ◽  
Subduction Zone ◽  
Subduction Zones ◽  
Large Scale ◽  
Large Earthquake ◽  
Force Balance ◽  
Mantle Flow ◽  
Plate Interface ◽  
Plate Motions

<p>The largest and most devastating earthquakes on Earth occur along subduction zones. Here, long-term plate motions are accommodated in cycles of strain accumulation and release. Episodic strain release occurs by mechanisms ranging from rapid earthquakes to slow-slip and quasi-static creep along the plate interface. Slip styles can vary between and within subduction zones, though it is unclear what controls margin-scale variability. Current approaches to seismo-tectonics primarily relate the stress state and seismogenesis at subduction margins to interface material properties and plate kinematics, constrained by recorded seismic slip, GPS motions and integrated strain. At larger spatio-temporal scales, significant progress has been made towards the understanding of subduction dynamics and emerging self-consistent plate motions, tectonics and stress coupling at plate margins. The margin stress state is ultimately linked to the force balance arising from interactions between the slab, mantle flow and upper plate. These mantle and lithosphere dynamics are thus expected to govern the tectonic regimes under which seismicity occurs. It remains unclear how these longer- and shorter-term perspectives can be reconciled. We review the aspects of large-scale subduction dynamics that control tectonic loading at plate margins, discuss possible influences on the stress state of the plate interface, and summarise recent advances in integrating the earthquake cycle and large-scale dynamics. It is plausible that variations in large-scale subduction dynamics could systematically influence seismicity, though it remains unclear to what degree this interplay occurs directly through the plate interface stress state and/or indirectly, corresponding to variation of other subduction zone characteristics. While further constraints of the geodynamic controls on the nature of the plate interface and their incorporation into probabilistic earthquake models is required, their ongoing development holds promise for an improved understanding of the global variation of large earthquake occurrence and their associated risk.</p>

scholarly journals The convergence history of India-Eurasia records multiple subduction dynamics processes

2020 ◽  
Vol 6 (19) ◽  
pp. eaaz8681 ◽  
Author(s):  
Adina E. Pusok ◽  
Dave R. Stegman
Keyword(s):  
Indian Ocean ◽  
Numerical Models ◽  
Plate Motion ◽  
Indian Plate ◽  
Geological Evidence ◽  
Sequence Of Events ◽  
History Of ◽  
The Indian Ocean ◽  
Fast Motion ◽  
Subduction System

During the Cretaceous, the Indian plate moved towards Eurasia at the fastest rates ever recorded. The details of this journey are preserved in the Indian Ocean seafloor, which document two distinct pulses of fast motion, separated by a noticeable slowdown. The nature of this rapid acceleration, followed by a rapid slowdown and then succeeded by a second speedup, is puzzling to explain. Using an extensive observation dataset and numerical models of subduction, we show that the arrival of the Reunion mantle plume started a sequence of events that can explain this history of plate motion. The forces applied by the plume initiate an intra-oceanic subduction zone, which eventually adds enough additional force to drive the plates at the anomalously fast speeds. The two-stage closure of a double subduction system, including accretion of an island arc at 50 million years ago, may help reconcile geological evidence for a protracted India-Eurasia collision.

scholarly journals The role of the margins in the dynamics of an active ice stream

1994 ◽  
Vol 40 (136) ◽  
pp. 527-538 ◽  
Author(s):  
K. A. Echelmeyer ◽  
W. D. Harrison ◽  
C. Larsen ◽  
J. E. Mitchell
Keyword(s):  
Numerical Models ◽  
Force Balance ◽  
Shear Stresses ◽  
List Type ◽  
Ice Stream ◽  
Strain Heating ◽  
Ice Stream B ◽  
Basal Shear ◽  
Marginal Zones

AbstractA transverse profile of velocity was measured across Ice Stream B, West Antarctica, in order to determine the role of the margins in the force balance of an active ice stream. The profile extended from near the ice-stream center line, through a marginal shear zone and on to the slow-moving ice sheet. The velocity profile exhibits a high degree of shear deformation within a marginal zone, where intense, chaotic crevassing occurs. Detailed analysis of the profile, using analytical and numerical models of ice flow, leads to the following conclusions regarding the roles of the bed and the margins in ice-stream dynamics: (i)The overall resistive drag on the ice stream is partitioned nearly equally between the margins and the bed and, thus, both are important in the force balance of the ice stream.(ii)The ice within the chaotic zone must be about 10 times softer than the ice in the central part of the ice stream.(iii)The average basal shear stress is 0.06 × 105 Pa. This implies that the entire bed cannot be blanketed by the weak, deformable till observed by Engelhardt and others (1990) near the center of the ice stream — there must be regions of increased basal drag.(iv)High strain rates and shear stresses in the marginal zones indicate that strain heating in the margins may be significant.While the exact quantitative values leading to these conclusions are somewhat model and location-dependent, the overall conclusions are robust. As such, they are likely to have importance for ice-stream dynamics in general.

Sign in / Sign up

Export Citation Format

Share Document