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 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 The evolution of basal mantle structure in response to supercontinent aggregation and dispersal

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Xianzhi Cao ◽  
Nicolas Flament ◽  
Ömer F. Bodur ◽  
R. Dietmar Müller
Keyword(s):  
Flow Model ◽  
Shear Velocity ◽  
Plate Motion ◽  
Mantle Flow ◽  
Mantle Structure ◽  
Exchange Material ◽  
Plate Reconstruction ◽  
Tectonic Terranes ◽  
Over Time ◽  
Low Shear

AbstractSeismic 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.

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.

scholarly journals Feedbacks between a non-Newtonian upper mantle, mantle viscosity structure and mantle dynamics

2020 ◽  
Vol 224 (2) ◽  
pp. 961-972
Author(s):  
A G Semple ◽  
A Lenardic
Keyword(s):  
Upper Mantle ◽  
Lower Mantle ◽  
Mantle Convection ◽  
Mantle Flow ◽  
Mantle Dynamics ◽  
Plate Motions ◽  
Low Viscosity ◽  
Mantle Viscosity ◽  
Flow Feature ◽  
Mantle Rheology

SUMMARY Previous studies have shown that a low viscosity upper mantle can impact the wavelength of mantle flow and the balance of plate driving to resisting forces. Those studies assumed that mantle viscosity is independent of mantle flow. We explore the potential that mantle flow is not only influenced by viscosity but can also feedback and alter mantle viscosity structure owing to a non-Newtonian upper-mantle rheology. Our results indicate that the average viscosity of the upper mantle, and viscosity variations within it, are affected by the depth to which a non-Newtonian rheology holds. Changes in the wavelength of mantle flow, that occur when upper-mantle viscosity drops below a critical value, alter flow velocities which, in turn, alter mantle viscosity. Those changes also affect flow profiles in the mantle and the degree to which mantle flow drives the motion of a plate analogue above it. Enhanced upper-mantle flow, due to an increasing degree of non-Newtonian behaviour, decreases the ratio of upper- to lower-mantle viscosity. Whole layer mantle convection is maintained but upper- and lower-mantle flow take on different dynamic forms: fast and concentrated upper-mantle flow; slow and diffuse lower-mantle flow. Collectively, mantle viscosity, mantle flow wavelengths, upper- to lower-mantle velocities and the degree to which the mantle can drive plate motions become connected to one another through coupled feedback loops. Under this view of mantle dynamics, depth-variable mantle viscosity is an emergent flow feature that both affects and is affected by the configuration of mantle and plate flow.

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>

Mantle convection beneath East Asia under global mantle convection spherical framework

2020 ◽  
Author(s):  
Qunfan Zheng ◽  
Huai Zhang
Keyword(s):  
Tibetan Plateau ◽  
East Asia ◽  
Upper Mantle ◽  
Mantle Convection ◽  
Mantle Flow ◽  
The Tibetan Plateau ◽  
Mantle Structure ◽  
Lithospheric Deformation ◽  
Mantle Dynamics ◽  
Flow Models

<p>East Asia is a tectonically active area on earth and has a complicated lithospheric deformation due to the western Indo-Asian continental collision and the eastern oceanic subduction mainly from Pacific plate. Till now, mantle dynamics beneath this area is not well understood due to its complex mantle structure, especially in the framework of global spherical mantle convection. Hence, a series of numerical models are conducted in this study to reveal the key controlling parameters in shaping the present-day observed mantle structure beneath East Asia under 3-D global mantle flow models. Global mantle flow models with coarse mesh are firstly applied to give a rough constraint on global mantle convection. The detailed description of upper mantle dynamics of East Asia is left with regional refined mesh. A power-law rheology and absolute plate field are applied subsequently to get a better constraint on the related regional mantle rheological structure and surficial boundary conditions. Thus, the refined and reasonable velocity and stress distributions of upper mantle beneath East Asia at different depths are retrieved based on our 3-D global mantle flow simulations. The derived large shallow mantle flow beneath the Tibetan Plateau causes significant lithospheric shear drag and dynamic topography that result in prominent tectonic evolution of this area. And the Indo–Asian collision may have induced mantle flow beneath the Indian plate and the different velocity structures between the asthenosphere and lithosphere indicate the shear drag of asthenospheric mantle. That may explain the reason that Indo–Asian collision has occurred for 50 Ma, and this collision can still continue to accelerate uplift in the Tibetan plateau. Finally, we also consider the possible implementations of 3-D numerical simulations combined with global lithosphere and deep mantle dynamics so as to discuss the relevant influences.</p>

Dynamics of upper and lower mantle subduction and its effects on the amplitude and pattern of mantle convection.

2021 ◽  
Author(s):  
Erik van der Wiel ◽  
Cedric Thieulot ◽  
Wim Spakman ◽  
Douwe van Hinsbergen
Keyword(s):  
Tectonic Evolution ◽  
Lower Mantle ◽  
Mantle Convection ◽  
Subduction Zones ◽  
Morphological Evolution ◽  
Phase Changes ◽  
Plate Motion ◽  
Plate Tectonic ◽  
Model Setup ◽  
Modeling Strategy

<p>Long-lived, Mesozoic-Cenozoic subduction zones such as the Pacific slab under the Americas and the Tethyan slab under Eurasia consumed thousands of kms of lithosphere of which remnants are detected in today’s mantle by seismic tomography. Major differences, however, in subduction zone evolution occurred between these systems which include strong variations in subduction rate, slab morphological evolution, and trench motion, which all appear mostly to be accommodated in the upper 1000 km of the mantle (van der Meer et al. 2018). Furthermore, sinking rates of slabs below this zone tend to be similar for different subduction systems and an order of magnitude smaller than their plate/subduction velocities. Working from the premise that the mantle rheology that accommodated these subduction systems is basically similar, although still poorly constrained, we test the hypothesis that the contrasting evolution of these subduction systems is primarily tied in with the global plate tectonic forcing of subduction.</p><p>It is generally accepted that plate motion is primarily driven by slab pull with contributions from ridge push, rather than the drag of the underlying mantle. If correct, numerical subduction models should be able to obtain upper as well as lower mantle subduction velocities and sinking rates similar to those reconstructed from geological records. We are at the start of this investigation and will present the numerical model setup, modeling strategy, and preliminary results of a 2-D subduction modelling experiment. We implement a 2D-cylindrical model setup for solving the conservation of momentum, mass and energy with the open-source geodynamics code ASPECT (Kronbichler et al. 2012) using a nonlinear visco-plastic rheology and including the major phase changes. Our focus is on the possible role of the absolute motion of the subducting and overriding plates in concert with slab pull variation reconstructed from plate tectonic evolution models, while in both subduction cases the same (partly nonlinear) mantle rheological processes are required to accommodate slab morphology change and slab sinking. Kinematic modelling constraints are derived from global plate tectonic evolution models, while the tomographically inferred present-day stage provides the end-stage geometry of slabs.</p><p>van der Meer, D. G., Van Hinsbergen, D. J., & Spakman, W. (2018). Atlas of the underworld: Slab remnants in the mantle, their sinking history, and a new outlook on lower mantle viscosity. Tectonophysics, 723, 309-448.</p><p>Kronbichler, M., Heister, T., & Bangerth, W. (2012). High accuracy mantle convection simulation through modern numerical methods. Geophysical Journal International, 191(1), 12-29.</p>

Sign in / Sign up

Export Citation Format

Share Document