We reevaluate the effects of slab age, speed, and dip on slab temperature with numerical models. The thermal parameter Φ = t v sin θ, where t is age, v is speed, and θ is angle, is traditionally used ...as an indicator of slab temperature. However, we find that an empirically derived quantity, in which slab temperature T ∝ log (t−av−b) , is more accurate at depths <120 km, with the constants a and b depending on position within the slab. Shallower than the decoupling depth (~70–80 km), a~1 and b~0, that is, temperature is dependent on slab age alone. This has important implications for the early devolatilization of slabs in the hottest (youngest) cases and for shallow slab seismicity. At subarc depths (~100 km), within the slab mantle, a~1 and b~0 again. However, for the slab crust, now a~0.5 and b~1, that is, speed has the dominant effect. This is important when considering the generation of arc magmatism, and in particular, slab melting and the generation of slab‐derived melange diapirs. Moving deeper into the Earth, the original thermal parameter performs well as a temperature indicator, initially in the core of the slab (the region of interest for deep water cycling). Finally, varying the decoupling depth between 40 and 100 km has a dominant effect on slab temperatures down to 140‐km depth, but only within the slab crust. Slab mantle temperature remains primarily dependent on age.
Key Points
New slab temperature indicator is more accurate/detailed than the thermal parameter for depths <120 km
Slab mantle temperature is controlled by slab age; subarc slab crust temperature is controlled primarily by subduction speed
A variable plate decoupling depth has a larger effect on slab crust temperature, in the subarc region, than other slab parameters
Subduction zones are not static features, but trenches retreat (roll back) or advance. Here, we investigate the dominant dynamic controls on trench migration by means of two‐ and three‐dimensional ...numerical modeling of subduction. This investigation has been carried out by systematically varying the geometrical and rheological model parameters. Our viscoplastic models illustrate that advancing style subduction is promoted by a thick plate, a large viscosity ratio between plate and mantle, and a small density contrast between plate and mantle or an intermediate width (w ∼ 1300 km). Advancing slabs dissipate ∼45% to ∼50% of the energy in the system. Thin plates with relatively low viscosity or relatively high density, or wide slabs (w ∼ 2300 km), on the other hand, promote subduction in the retreating style (i.e., slab roll‐back). The energy dissipated by a retreating slab is ∼35% to ∼40% of the total dissipated energy. Most of the energy dissipation occurs in the mantle to accommodate the slab motion, whereas the lithosphere dissipates the remaining part to bend and “unbend.” With a simple scaling law we illustrate that this complex combination of model parameters influencing trench migration can be reduced to a single one: plate stiffness. Stiffer slabs cause the trench to advance, whereas more flexible slabs lead to trench retreat. The reason for this is that all slabs will bend into the subduction zone because of their low plastic strength near the surface, but stiff slabs have more difficulty “unbending” at depth, when arriving at the 660‐km discontinuity. Those bent slabs tend to cause the trench to advance. In a similar way, variation of the viscoplasticity parameters in the plate may change the style of subduction: a low value of friction coefficient weakens the plate and results in a retreating style, while higher values strengthen the plate and promote the advancing subduction style. Given the fact that also on Earth the oldest (and therefore probably stiffest) plates have the fastest advancing trenches, we hypothesize that the ability of slabs to unbend after subduction forms the dominant control on trench migration.
Some oceanic volcano chains violate the predictions of the hotspot hypothesis for geographic age progressions. One mechanism invoked to explain these observations is small‐scale sublithospheric ...convection (SSC). In this study, we explore this concept in thermo‐chemical, 3D‐numerical models. Melting due to SSC is shown to emerge in elongated features (∼750 km) parallel to plate motion and not just at a fixed spot; therefore volcanism occurs in chains but not with hotspot‐like linear age progressions. The seafloor age at which volcanism first occurs is sensitive to mantle temperature, as higher temperatures increase the onset age of SSC because of the stabilizing influence of thicker residue from previous mid‐ocean ridge melting. Mantle viscosity controls the rate of melt production with decreasing viscosities leading to more vigorous convection and volcanism. Calculations predict many of the key observations of the Pukapuka ridges, and the volcano groups associated with the Line, Cook‐Austral, and Marshall Islands.
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.
Mesozoic‐Cenozoic rifting between Greenland and North America created the Labrador Sea and Baffin Bay, while leaving preserved continental lithosphere in the Davis Strait, which lies between them. ...Inherited crustal structures from a Palaeoproterozoic collision have been hypothesized to account for the tectonic features of this rift system. However, the role of mantle lithosphere heterogeneities in continental suturing has not been fully explored. Our study uses 3‐D numerical models to analyze the role of crustal and subcrustal heterogeneities in controlling deformation. We implement continental extension in the presence of mantle lithosphere suture zones and deformed crustal structures and present a suite of models analyzing the role of local inheritance related to the region. In particular, we investigate the respective roles of crust and mantle lithospheric scarring during an evolving stress regime in keeping with plate tectonic reconstructions of the Davis Strait. Numerical simulations, for the first time, can reproduce first‐order features that resemble the Labrador Sea, Davis Strait, Baffin Bay continental margins, and ocean basins. The positioning of a mantle lithosphere suture, hypothesized to exist from ancient orogenic activity, produces a more appropriate tectonic evolution of the region than the previously proposed crustal inheritance. Indeed, the obliquity of the continental mantle suture with respect to extension direction is shown here to be important in the preservation of the Davis Strait. Mantle lithosphere heterogeneities are often overlooked as a control of crustal‐scale deformation. Here, we highlight the subcrust as an avenue of exploration in the understanding of rift system evolution.
Key Points
The role of mantle sutures during Mesozoic rifting and Cenozoic ocean basin development in Laurentia has not been studied
We present 3‐D models that reproduce first‐order tectonics of Labrador Sea, Davis Strait, and Baffin Bay through mantle suture reactivation
The obliquity of the mantle suture to the extension preserves continental lithosphere, interpreted as the creation of the Davis Strait
Many volcano chains in the Pacific do not follow the most fundamental predictions of hot spot theory in terms of geographic age progressions. One possible explanation for non–hot spot intraplate ...volcanism is small‐scale sublithospheric convection (SSC), and we explore this concept using 3‐D numerical models that simulate melting with rheology laws that account for the effects of dehydration. SSC spontaneously self‐organizes beneath relatively mature oceanic lithosphere. Whenever this lithosphere is sufficiently young and thin, SSC replaces the shallow layer of harzburgite, which was formed by partial melting at the mid‐ocean ridge, with fresh peridotite. This mechanism enables magma generation without any preexisting thermochemical anomalies. However, the additional effect of melting‐induced dehydration to stiffen the harzburgite requires lower background viscosities to allow for vigorous SSC, overturn of the compositional stratification, and related magmatism. The intrinsic stiffness of the dehydrated harzburgite furthermore restricts penetration of SSC into very shallow and cooler levels. On the one hand, such a restriction precludes high degrees of melting, but on the other hand, it slows asthenospheric cooling and thus prolongs the duration of melting (to ∼25 Ma). Volcanism over such an elongated melting anomaly continues for at least 10–20 Ma and occurs on seafloor ages of ∼20 to ∼60 Ma. These seafloor ages increase with increasing mantle temperature due to the effect of forming a thicker harzburgite layer from more extensive mid‐ocean ridge melting. The long durations of volcanism predicted reconcile observations of extended activity of individual seamounts and synchronous activity over great distances along some volcanic chains. SSC thus gives an explanation for previously enigmatic volcano ages along the Line Islands and the Gilbert and Pukapuka ridges, as well as along the individual subchains of the Wakes, Marshalls, and Cook‐Australs.
Dynamic models of subduction and continental collision are used to predict dynamic topography changes on the overriding plate. The modelling results show a distinct evolution of topography on the ...overriding plate, during subduction, continental collision and slab break-off. A prominent topographic feature is a temporary (few Myrs) basin on the overriding plate after initial collision. This "collisional mantle dynamic basin" (CMDB) is caused by slab steepening drawing, material away from the base of the overriding plate. Also, during this initial collision phase, surface uplift is predicted on the overriding plate between the suture zone and the CMDB, due to the subduction of buoyant continental material and its isostatic compensation. After slab detachment, redistribution of stresses and underplating of the overriding plate cause the uplift to spread further into the overriding plate. This topographic evolution fits the stratigraphy found on the overriding plate of the Arabia-Eurasia collision zone in Iran and south east Turkey. The sedimentary record from the overriding plate contains Upper Oligocene-Lower Miocene marine carbonates deposited between terrestrial clastic sedimentary rocks, in units such as the Qom Formation and its lateral equivalents. This stratigraphy shows that during the Late Oligocene-Early Miocene the surface of the overriding plate sank below sea level before rising back above sea level, without major compressional deformation recorded in the same area. Our modelled topography changes fit well with this observed uplift and subsidence.
Orogenic crustal anatexis is a still poorly understood process due to the complexity of the thermal and geodynamical interaction between mantle and crustal processes during and after continental ...collision. Here we present a novel conceptual model for the formation of granite‐migmatite belts: we propose that convective thinning of the lithosphere results in minor amounts of partial melts within the lowermost crust that trigger further instabilities. This will lead to positive feedback effects between melt weakening, mantle upwelling, and wholesale mantle lithosphere removal, causing a strong pulse of mantle and crustal melting. We test this model numerically, and results show that this process, taking between 20 and 50 Myr in total, can explain the temporal evolution of melting in granite‐migmatite zones and associated mantle‐derived mafic rocks and provides a heat source for crustal melting without the need for other processes, such as slab break‐off or increased radiogenic heating. Furthermore, the generation of a refractory residue after mantle and crustal melting is also shown to control the progress of the lithospheric mantle removal, providing another feedback mechanism between melting and lithospheric reequilibration.
Key Points
Partial lithosphere melting leads to positive feedback with convective lithosphere thinning
Postcollisional lithosphere thermal relaxation initiates lower crustal anatexis
Process explains temporal order of migmatites, granitic plutons, and intervening mantle melts