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>

Controls on the present-day dynamic topography predicted from mantle flow models since 410 Ma

2021 ◽  
Author(s):  
An Yang ◽  
Ting Yang
Keyword(s):  
Rayleigh Number ◽  
Mantle Convection ◽  
Subduction Zones ◽  
Surface Heat ◽  
Heat Fluxes ◽  
Phase Changes ◽  
Plate Motion ◽  
Dynamic Topography ◽  
Convergence Region ◽  
Temperature Dependent Viscosity

Summary Mantle convection induces dynamic topography, the lithosphere's surface deflections driven by the vertical stresses from sub-lithospheric mantle convection. Dynamic topography has important influences on a range of geophysical and geological observations. Here, we studied controls on the Earth's dynamic topography through three-dimensional spherical models of mantle convection, which use reconstructed past 410 Myr global plate motion history as time-dependent surface mechanical boundary condition. The numerical model assumes the extended-Boussinesq approximation and includes strongly depth- and temperature-dependent viscosity and phase changes in the mantle. Our results show that removing the chemical layer above the CMB and including depth-dependent thermal expansivity have both a limited influence on the predicted present-day dynamic topography. Considering phase transitions in our models increases the predicted amplitude of dynamic topography, which is mainly influenced by the 410 km exothermic phase transition. The predicted dynamic topography is very sensitive to shallow temperature-induced lateral viscosity variations (LVVs) and Rayleigh number. The preservation of LVVs significantly increases the negative dynamic topography at subduction zones. A decrease (or increase) of Rayleigh number increases (or decreases) the predicted present-day dynamic topography considerably. The dynamic topography predicted from the model considering LVVs and with a Rayleigh number of 6 × 108 is most compatible with residual topography models. This Rayleigh number is consistent with the convective vigor of the Earth as supported by generating more realistic lower mantle structure, slab sinking rate, and surface and CMB heat fluxes. The evolution of the surface heat flux pattern is similar to the long-term eustatic sea-level change. Before the formation of Pangea, large negative dynamic topography formed between the plate convergence region of Gondwana and Laurussia. The predicted dynamic topography similar to that of present day has already emerged by about 262 Ma. Powers for degrees 1–3 dynamic topography at 337 Ma and 104 Ma which correspond to times of higher plate velocities and higher surface heat fluxes are larger.

scholarly journals Improving global paleogeography since the late Paleozoic using paleobiology

2017 ◽  
Author(s):  
Wenchao Cao ◽  
Sabin Zahirovic ◽  
Nicolas Flament ◽  
Simon Williams ◽  
Jan Golonka ◽  
...  
Keyword(s):  
Ocean Circulation ◽  
Mantle Convection ◽  
Growth Models ◽  
Plate Motion ◽  
Late Paleozoic ◽  
Plate Tectonic ◽  
Dynamic Topography ◽  
Tectonic Reconstructions ◽  
Marine Boundaries ◽  
Continental Growth

Abstract. Paleogeographic reconstructions are important to understand Earth's tectonic evolution, past eustatic and regional sea level change, hydrocarbon genesis, and to constrain and interpret the dynamic topography predicted by time-dependent global mantle convection models. Several global paleogeographic maps have been compiled and published but they are generally presented as static maps with varying temporal resolution and fixed spatial resolution. Existing global paleogeographic maps are also tied to a particular plate motion model, making it difficult to link them to alternative digital plate tectonic reconstructions. To address this limitation, we developed a workflow to reverse-engineer global paleogeographic maps to their present-day coordinates and enable them to be linked to any tectonic reconstruction. Published paleogeographic compilations are also tied to fixed input datasets. We used fossil data from the Paleobiology Database to identify inconsistencies between fossils paleo-environments and published paleogeographic maps, and to improve the location of inferred terrestrial-marine boundaries by resolving these inconsistencies. As a result, the overall consistency ratio between the paleogeography and fossil collections was improved from 76.9 % to 96.1 %. We estimated the surface areas of global paleogeographic features (shallow marine environments, landmasses, mountains and ice sheets), and reconstructed the global continental flooding history since the late Paleozoic based on the amended paleogeographies. Finally, we discuss the relationships between emerged land area and total continental crust area through time, continental growth models, and strontium isotope (87Sr/86Sr) signatures in ocean water. Our study highlights the flexibility of digital paleogeographic models linked to state-of-the-art plate tectonic reconstructions in order to better understand the interplay of continental growth and eustasy, with wider implications for understanding Earth's paleotopography, ocean circulation, and the role of mantle convection in shaping long-wavelength topography.

Plate tectonic chain reaction constrained from noise in the Cretaceous Quiet Zone

2020 ◽  
Author(s):  
Derya Gürer ◽  
Roi Granot ◽  
Douwe J.J. van Hinsbergen
Keyword(s):  
Subduction Zones ◽  
Kinematic Model ◽  
Plate Boundary ◽  
Western Mediterranean ◽  
Plate Motion ◽  
Plate Boundaries ◽  
Plate Tectonic ◽  
Test Bed ◽  
Chain Reaction ◽  
Tectonic Plates

<p>The relative motions of the tectonic plates show remarkable variation throughout Earth’s history. Major changes in relative motion between the tectonic plates are traditionally viewed as spatially and temporally isolated events linked to forces acting on plate boundaries (i.e., formation of same-dip double subduction zones, changes in the strength of the boundary), or thought to be associated with mantle dynamics. A Cretaceous global plate reorganization event has been postulated to have affected all major plates. The Cretaceous ‘swing’ in Africa-Eurasia relative plate motion provides an ideal test-bed for assessing the temporal and spatial evolution of both relative plate motions and surrounding geological markers. Here we show a novel plate kinematic model for the closure of the Tethys Ocean by implementing intra-Cretaceous Quiet Zone time markers and combine the results with the geological constraints found along the convergent plate boundary. Our results allow to assess the order, causes and consequences of geological events and unravel a chain of tectonic events that set off with the onset of horizontally-forced double subduction ~105 Myr ago, followed by a 40 Myr long period of acceleration of the Africa relative to Eurasia that peaked at 80 Myr ago (at rates four times as high as previously predicted). This acceleration, which was likely caused by the pull of two same-dip subduction zones was followed by a sharp decrease in plate velocity, when double subduction terminated with ophiolite obduction onto the African margin. These tectonic forces acted on the eastern half of the Africa-Eurasia plate boundary, which led to counterclockwise rotation of Africa and sparked new subduction zones in the western Mediterranean region. Our analysis identifies the Cretaceous double subduction episode between Africa and Eurasia as a link in the global plate tectonic chain reaction and provides a dynamic view on plate reorganizations.</p>

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>

Using geodynamic modeling to test plate tectonic scenarios of the Mediterranean-Alpine area

2020 ◽  
Author(s):  
Boris Kaus ◽  
Eline Le Breton ◽  
Georg Reuber ◽  
Christian Schuler
Keyword(s):  
Subduction Zones ◽  
Thermal Structure ◽  
Shear Zones ◽  
Western Mediterranean ◽  
Plate Motion ◽  
Plate Tectonic ◽  
Test Plate ◽  
Alpine Area ◽  
Tectonic Reconstruction ◽  
The Mediterranean

<p>Using geological and geophysical data, it is possible to reconstruct the past motion of tectonic plates involved in the Alpine orogeny and propose possible scenarios for their geological evolution. However, those scenarios have not yet been tested for geodynamic consistency.</p><p>Here, we perform 3D thermomechanical geodynamic simulations of the Mediterranean-Alpine area starting with a plate tectonic reconstruction at 20 Ma based on the work of Le Breton et al. (2017). The models include viscoelastoplastic rheologies and a free surface, and thus simulate the spontaneous occurrence of shear zones as well as the development of topography. Whereas some aspects of the tectonic reconstruction are well constrained (i.e. past position of the plates and subduction-collision fronts), many details such as the dip and length of the subducted plates, their thermal structure as well as their rheology, are unknown. The models are run forward in time to see to which extent they are consistent with the kinematic reconstructions. Perhaps unsurprisingly, our initial modelling attempts show a wide variety of behavior, including slab break off events and slab rollbacks in the wrong directions. Yet, in all cases tested so far, the model evolution does not reproduce the present-day geological setting, with Adria frequently moving too far towards the east and breaking apart internally, frequently no Alpine chain forming and in some cases new subduction zones developing within the Western Mediterranean that swallow Sardinia and Corsica.</p><p>Reproducing geological scenarios with thermomechanical geodynamic modelling thus requires substantial additional work, both from the modelling side (testing the effect of uncertain parameters on the behaviour of plates and subduction zones), as well as from the plate reconstruction side (assessing which parameters are well constrained and need to be reproduced).  Nevertheless, interesting insights can already be obtained from our models, and in our presentation, we will highlight some of the links between interacting subducting plates and plate motion.</p><p>Le Breton E, Handy MR, Molli G, Ustaszewski K (2017) Post-20 Ma Motion of the Adriatic Plate: New Constraints From Surrounding Orogens and Implications for Crust-Mantle Decoupling. Tectonics 36:3135–3154. doi: 10.1002/2016TC004443</p><div> <div> <div> </div> </div> </div>

Evolution of Orogenic Plateaus at Subduction Zones - Sinking and raising the southern margin of the Central Anatolian Plateau

2018 ◽  
Author(s):  
David Fernández-Blanco
Keyword(s):  
Tectonic Evolution ◽  
Subduction Zones ◽  
Inversion Tectonics ◽  
High Discharge ◽  
Anatolian Plateau ◽  
Surface Uplift ◽  
Southern Margin ◽  
Wide Range ◽  
Forearc Basins ◽  
Orogenic Plateaus

Orogenic plateaus have raised abundant attention amongst geoscientists during the last decades, offering unique opportunities to better understand the relationships between tectonics and climate, and their expression on the Earth’s surface.Orogenic plateau margins are key areas for understanding the mechanisms behind plateau (de)formation. Plateau margins are transitional areas between domains with contrasting relief and characteristics; the roughly flat elevated plateau interior, often with internally drained endorheic basins, and the external steep areas, deeply incised by high-discharge rivers. This thesis uses a wide range of structural and tectonic approaches to investigate the evolution of the southern margin of the Central Anatolian Plateau (CAP), studying an area between the plateau interior and the Cyprus arc. Several findings are presented here that constrain the evolution, timing and possible causes behind the development of this area, and thus that of the CAP. After peneplanation of the regional orogeny, abroad regional subsidence took place in Miocene times in the absence of major extensional faults, which led to the formation of a large basin in the northeast Mediterranean. Late Tortonian and younger contractional structures developed in the interior of the plateau, in its margin and offshore, and forced the inversion tectonics that fragmented the early Miocene basin into the different present-day domains. The tectonic evolution of the southern margin of the CAP can be explained based on the initiation of subduction in south Cyprus and subsequent thermo-mechanical behavior of this subduction zone and the evolving rheology of the Anatolian plate. The Cyprus slab retreat and posterior pull drove subsidence first by relatively minor stretching of the crust and then by its flexure. The growth by accretion and thickening of the upper plate, and that of the associated forearc basins system, caused by accreting sediments, led to rheological changes at the base of the crust that allowed thermal weakening, viscous deformation, driving subsequent surface uplift and raising the modern Taurus Mountains. This mechanism could be responsible for the uplifted plateau-like areas seen in other accretionary margins. ISBN: 978-90-9028673-0

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.

scholarly journals On the enigmatic birth of the Pacific Plate within the Panthalassa Ocean

2016 ◽  
Vol 2 (7) ◽  
pp. e1600022 ◽  
Author(s):  
Lydian M. Boschman ◽  
Douwe J. J. van Hinsbergen
Keyword(s):  
Triple Junction ◽  
Tectonic Evolution ◽  
Plate Boundary ◽  
Pacific Plate ◽  
Plate Tectonic ◽  
Point Location ◽  
Marine Magnetic Anomalies ◽  
The Pacific ◽  
History Of ◽  
Ridge System

The oceanic Pacific Plate started forming in Early Jurassic time within the vast Panthalassa Ocean that surrounded the supercontinent Pangea, and contains the oldest lithosphere that can directly constrain the geodynamic history of the circum-Pangean Earth. We show that the geometry of the oldest marine magnetic anomalies of the Pacific Plate attests to a unique plate kinematic event that sparked the plate’s birth at virtually a point location, surrounded by the Izanagi, Farallon, and Phoenix Plates. We reconstruct the unstable triple junction that caused the plate reorganization, which led to the birth of the Pacific Plate, and present a model of the plate tectonic configuration that preconditioned this event. We show that a stable but migrating triple junction involving the gradual cessation of intraoceanic Panthalassa subduction culminated in the formation of an unstable transform-transform-transform triple junction. The consequent plate boundary reorganization resulted in the formation of a stable triangular three-ridge system from which the nascent Pacific Plate expanded. We link the birth of the Pacific Plate to the regional termination of intra-Panthalassa subduction. Remnants thereof have been identified in the deep lower mantle of which the locations may provide paleolongitudinal control on the absolute location of the early Pacific Plate. Our results constitute an essential step in unraveling the plate tectonic evolution of “Thalassa Incognita” that comprises the comprehensive Panthalassa Ocean surrounding Pangea.

The effect of temperature-dependent thermal parameters in thermal models of subduction zones

2021 ◽  
Author(s):  
Iris van Zelst ◽  
Timothy J. Craig ◽  
Cedric Thieulot
Keyword(s):  
Subduction Zones ◽  
Thermal Parameters ◽  
Dislocation Creep ◽  
Seismogenic Zone ◽  
Effect Of Temperature ◽  
Temperature Dependent ◽  
Thermal Models ◽  
Intermediate Depth ◽  
Dehydration Reactions ◽  
Model Setup

<p>The thermal structure of subduction zones plays an important role in the seismicity that occurs there with e.g., the downdip limit of the seismogenic zone associated with particular isotherms (350 °C - 450 °C) and intermediate-depth seismicity linked to dehydration reactions that occur at specific temperatures and pressures. Therefore, accurate thermal models of subduction zones that include the complexities found in laboratory studies are necessary. One of the often-ignored effects in models is the temperature-dependence of the thermal parameters such as the thermal conductivity, heat capacity, and density.<span> </span></p><p>Here, we build upon the model setup presented by Van Keken et al., 2008 by including temperature-dependent thermal parameters to an otherwise clearly constrained, simple model setup of a subducting plate. We consider a fixed kinematic slab dipping at 45° and a stationary overriding plate with a dynamic mantle wedge. Such a simple setup allows us to isolate the effect of temperature-dependent thermal parameters. We add a more complex plate cooling model for the oceanic plate for consistency with the thermal parameters.<span> </span></p><p>We test the effect of temperature-dependent thermal parameters on models with different rheologies, such as an isoviscous wedge, diffusion and dislocation creep. We find that slab temperatures can change by up to 65 °C which affects the location of isotherm depths. The downdip limit of the seismogenic zone defined by e.g., the 350 °C isotherm shifts by approximately 4 km, thereby increasing the maximum possible rupture area of the seismogenic zone. Similarly, the 600 °C isotherm is shifted approximately 30 km deeper, affecting the depth at which dehydration reactions and hence intermediate-depth seismicity occurs. Our results therefore show that temperature-dependent thermal parameters in thermal models of subduction zones cannot be ignored when studying subduction-related seismicity.<span> </span></p>

Sign in / Sign up

Export Citation Format

Share Document