Overriding Lithosphere Research Articles

Numerical Modeling of Melting Processes During Slab Break-off: Insights Into Tectonic Setting for Massif-Type Anorthosites

Abstract The concept that lithosphere detachment or break-off has long been conceived as a viable mechanism to explain prominent geological phenomena in Earth’s crust and the surface. One of the strengths of slab delamination mechanism is that it can account for the extensive magmatism in active orogenic belts due to the upwelling of the asthenosphere after the slab break-off. However, in the last 20 years, geodynamic simulations show that the inflow of the asthenosphere upon slab break-off is insufficient to cause significant melting of the overriding lithosphere adjacent to the slab window. The primary reasons include the occurrence of slab break-off at a location that is too deep to effectively heat the overriding lithospheric mantle. Another factor is the presence of a thin film of crustal material that is retained during the slab break- off, inhibiting a significant thermal perturbation within the lithosphere. In this work, we couple petrological–thermomechanical simulations with magmatic melting processes to examine the lithospheric melting and surface lithological expression associated with slab break-off. Our work shows that in the early Earth when the mantle temperature is relatively higher, shallow slab break-off can give rise to significant lithospheric melting during the development of slab break-off. Moreover, because the slab becomes weaker in the earlier hotter mantle, it may break-off prior to the stage of continental collision, thus the magmatism it induced may not give a direct constraint on the time of continental collision. Our study has implications for the interpretation of geological and tomography studies in orogenic belts. It also provides insights into reconciling conflicts between geodynamic and geological studies regarding slab break off-induced melting and magmatism.

Toroidal flow around the Tonga slab moved the Samoan plume during the Pliocene

Abstract Age-progressive seamount tracks generated by lithospheric motion over a stationary mantle plume have long been used to reconstruct absolute plate motion (APM) models. However, the basis of these models requires the plumes to move significantly slower than the overriding lithosphere. When a plume interacts with a convergent or divergent plate boundary, it is often deflected within the strong local mantle flow fields associated with such regimes. Here, we examined the age progression and geometry of the Samoa hotspot track, focusing on lava flow samples dredged from the deep flanks of seamounts in order to best reconstruct when a given seamount was overlying the mantle plume (i.e., during the shield-building stage). The Samoan seamounts display an apparent local plate velocity of 7.8 cm/yr from 0 to 9 Ma, 11.1 cm/yr from 9 to 14 Ma, and 5.6 cm/yr from 14 to 24 Ma. Current fixed and mobile hotspot Pacific APM models cannot reproduce the geometry of the Samoa seamount track if a long-term fixed hotspot location, currently beneath the active Vailulu’u Seamount, is assumed. Rather, reconstruction of the eruptive locations of the Samoan seamounts using APM models indicates that the surface expression of the plume migrated ~2° northward in the Pliocene. Large-scale mantle flow beneath the Pacific Ocean Basin cannot explain this plume migration. Instead, the best explanation is that toroidal flow fields—generated by westward migration of the Tonga Trench and associated slab rollback—have deflected the conduit northward over the past 2–3 m.y. These observations provide novel constraints on the ways in which plume-trench interactions can alter hotspot track geometries.

Neotethyan Ankara Melange, central Turkey: Formation by accretion of seamounts and supra-subduction zone ophiolites in an oceanic fore-arc setting

The Neotethyan Ankara Melange is a global reference for oceanic accretionary processes preserved in a collisional orogen. Here we show, based on eight representative outcrops around Ankara, that the melange encompasses stratigraphically coherent volcanic-sedimentary successions of mainly Late Jurassic-Early Cretaceous age. Our results indicate that the melange has a relatively organised tectonostratigraphy, in contrast to previous interpretations that emphasised its chaotic nature. Coherent volcanic-sedimentary successions of Late Jurassic-Early Cretaceous age are interpreted as flank facies of a large oceanic seamount (probably plume related) and its capping carbonate platform (atoll), including reef-and back-reef facies. Remnants of other seamounts are dated as Late Triassic-Early Jurassic and Late Cretaceous. Variably dismembered supra-subduction zone ophiolites of latest Early Jurassic-Middle Jurassic age, including boninites, formed in an oceanic fore-arc setting, similar to that of the Cenozoic Izu-Bonin fore arc, NW Pacific. The melange in the Ankara region partly accreted during the Early Cretaceous (pre-Albian) in response to collision of the inferred large seamount with the oceanic (ophiolitic) fore arc. During the collision which involved subduction-erosion slices of proximal–distal seamount flank lithologies accreted sequentially, whereas the seamount core subducted. Also during the collision, the distal outer edge of the ophiolitic fore arc, mainly mantle harzburgite, detached from the over-riding lithosphere and was incorporated into the melange. Structural evidence emphasises the importance of tectonic processes (rather than olistostromes) in the melange accretion. The accretionary prism in the outcrops studied was covered by arc-related sediments and volcanics during the Late Cretaceous and deformed by collision-related processes during latest Cretaceous-Paleogene.

Strain and stress gradients through the backarc regions of Miocene western Japan: A new type of arc-parallel extension

The overriding lithosphere is occasionally stretched parallel to the consuming plate boundary. Such arc-parallel extension manifests itself mostly in forearc regions. Here, we report unique arc-parallel extension, which affected only the backarc regions with arc-parallel strain and stress gradients. The extension took place in the backarc regions of the northern Ryukyu and western SW Japan arcs in the Middle Miocene. The extension was evidenced by NW-SE trending normal faults in the extensional zone and paleostresses that were determined from dike orientations. To determine the age of and spatial extent of this zone, we analyzed fault-slip data from the Amakusa region, northern Ryukyu arc, where early middle Miocene intrusions allowed us to judge the magmatism to be older than the faults that affected Cretaceous and Eocene formations. The slip senses of the faults have been controversial in Amakusa. As a result, the region was found to have belonged to the extensional zone. The extensional zone is 300–400 km long and ∼100 km wide. We also found that the arc-parallel extension resulted in crustal strain with a southwestward increasing trend along the zone. The strain magnitudes showed a negative correlation with distance from the Okinawa Trough. This correlation is explained by the increasing stress magnitude toward the trough. We suggest that the stalled rift propagation of the northern Okinawa Trough caused the unique arc-parallel extension in the Middle Miocene.

Dynamics of the Subducted Izanagi-Pacific Plates Since the Mesozoic and Its Implications for the Formation of Big Mantle Wedge Beneath Eastern Asia

The slab dynamics of the subducted Izanagi-Pacific plate is still a subject of controversy and its relationship with the tectonic evolution of Eastern Asia remains not well explored. Here, we perform 3-D global convection models to investigate the slab dynamics of the Izanagi-Pacific plate beneath Eastern Asia since the Mesozoic time. We introduce a tracking technique in numerical models to explicitly distinguish the Izanagi slab and the Pacific slab during their subduction processes. We find that all subducted Izanagi slabs have completely fallen into the lower mantle until the late Cenozoic and the stagnant slabs currently observed at the mantle transition zone depth beneath Eastern Asia are entirely from the Pacific plate. We also find that multiple slab stagnation events have occurred during the subduction of the Izanagi plate in the Mesozoic time (∼150–120 Ma, 90–70 Ma) with a timescale of tens of million years. The stagnation of the subducted slabs facilitates the formation of a big mantle wedge beneath the overriding lithosphere and the time periods of the mantle wedge are consistent with the episodes of magmatic activities in Eastern Asia.

The delamination of lower crust in continental back-arc basin: Evidence from Sr isotope and elemental compositions of plagioclase and clinopyroxene in andesites from Kueishantao, north of Taiwan, China

Subduction zones (arcs and back-arcs) are major sites for elemental cycling via slab dehydration and subsequent mantle metasomatism and melting; many models have been suggested by researchers to explain these processes. However, the influence of the overriding lithosphere, especially the lower crust, on the generation of back-arc magmas during the early spreading stage remains largely unknown. Kueishantao (KST), located at the tectonic junction of the southernmost part of the Okinawa Trough (SPOT) and the Ryukyu arc, is an ideal setting to investigate these issues because of the obvious differences between the subducted slab and the overriding continental lithosphere. To better address the contribution of the overriding lithosphere to a nascent back-arc basin, we examined the in situ major and trace element compositions of plagioclase (Pl) and clinopyroxene (Cpx) phenocrysts and the in situ 87Sr/86Sr of Pl phenocrysts in andesites from KST.The Pl and Cpx in KST andesites show obvious evolution (An contents of Pl, 39–91; Mg# of Cpx, 73–92). The valid thermobarometers of Pl and Cpx imply crystallization temperatures and depths of 1059–1117 °C and 1–14 km, respectively, similar to those of the SPOT. The compositions of Cpx suggest that the KST andesitic magma belonged to the subalkaline, tholeiitic and calc-alkaline magma series. The δEu values of Pl are 9.27–25.73, and those of Cpx are 0.50–0.92, indicating that the magma experienced Pl crystallization and separation. The 87Sr/86Sr values of Pl range from 0.705630 to 0.706985 and remain constant with varying An contents (52–85), indicating an enriched magma source not resulting from assimilation of continental crust. Taking into consideration all geochemical, trace elemental, 87Sr/86Sr isotope and geophysical constraints, the principal magma source of KST andesites was probably the lower continental crust, which foundered into the convecting mantle and subsequently melted and interacted with asthenospheric mantle materials. The over-thickened crust, mantle corner flow and addition of subducted fluids in the KST area led to delamination. The structure and geochemistry of Pl and Cpx imply a simple two- or multilayer magma chamber system in the continental crust and a mixing magma process between a relatively silicic magma (H2O-saturated and high fO2) and a primitive mafic magma at shallow depths.

Plate corner subduction and rapid localized exhumation: Insights from3Dcoupled geodynamic and geomorphological modelling

AbstractRapid, localized exhumation has been reported at many plate corners between adjacent subduction/collision segments. Here we use a fully‐coupled geodynamic and geomorphological modelling approach to investigate overriding plate deformation and resulting rock uplift patterns in these narrow, cuspate regions. In this study, we focus on the effects of internal deformation within a subducting convex‐upward‐shaped indenter and the strength of the interface between the upper and downgoing plate. The strongest localization of high rock uplift rates in the region above the indenter apex is predicted in experiments with a deformable lower plate, a weak interface layer and lateral shortening accommodated only by subduction (i.e., without an upper plate advance component). Our results suggest that bull’s eye shaped structures characterized by young thermochronological ages can, in principle, be reproduced numerically when taking into account a non‐rigid subducting plate together with complex brittle‐ductile rheology and stratification of the overriding lithosphere and realistically implemented fluvial erosion at its surface.

Melt Percolation, Concentration and Dyking in the Hawaiian Mantle Plume and Overriding Lithosphere: Links to the Evolution of Lava Composition along the Volcanic Chain

Abstract Oceanic island basalts and related magmatic rocks from Hawaii are derived from a compositionally heterogeneous mantle plume. Here we describe how this heterogeneity results from the transport of filaments of a specific composition in the plume, representing a relatively small volume of rocks (~15 %) interbedded inside a dry peridotite mantle. Four types of filaments are considered: sub-primitive mantle, ultralow-velocity zone, fertilized-harzburgite and eclogite type filaments. We present a model that describes the flow within a plume and the stress field in the overriding viscoelastic lithosphere and that can determine, from depth to the surface, the melting rate, composition and trajectory of melts produced within each type of filament. Our model shows that (1) the filaments melt at a depth corresponding to >5 GPa, where the temperature gap between the solidus and liquidus is narrow (~40–80 °C), and (2) the volume of filaments is small relative to the total volume of mantle, which therefore allows the latent heat required for the partial melting to be provided via conduction inside the hot plume. The primitive melts produced inside the filaments, occasionally mixed with the melt derived from an eclogite filament, represent a volume comparable with that expected in a plume composed only of dry peridotite that partially melts to a degree of ~10 % at the interface between the spinel and garnet fields (60–70 km depth). In particular, in the centre of the plume, sub-primitive mantle filaments produce up to 30 % tholeiite–picrite melts, whereas in fertilized-harzburgite filaments, the mantle melts completely to produce a melt having a meimechite-like composition. A key finding is that the fractional crystallization of these melts probably forms the so-called ‘primary mantle-derived alkaline magmas’ along with dunites and olivine-rich cumulates. Our plume model shows that the mantle flow divides into two parts. The first corresponds to hot flowlines that originate at a depth of ~200 km and at a distance of less than 25 km from the plume axis. Along these flowlines, when the mantle reaches a pressure of 5 GPa, the partially molten horizon in filaments is sufficiently thick for the interstitial melt to be squeezed out via dykes. This melt eventually ponds as sills in a subrectangular zone that is located inside the overlying lithosphere, between 70 and 50 km depth and centred over a distance of less than 40 km on either side of the axis. This zone is designated as the shield magmatic reservoir. The volatile-rich melt inside the sills infiltrates the surrounding mantle lithosphere and partially melts it. After ~0·1 Myr, the melt resumes its vertical ascent via dykes and eventually ponds and differentiates within subcrustal magma chambers located below active shield volcanoes. This sequence of processes matches the expected volume, petrology and geochemistry recorded for shield volcanoes. The second part of the melt flow does not pond within the shield magmatic reservoir. Rather, the mantle cold flowlines, originating at ~200 km depth and at 25–35 km from the plume axis, discharge their interstitial melt through dykes that were initially generated deeper, at ~5 GPa. The melt reaches the Moho at 100–150 km from the plume axis, where it forms magmatic bodies within which the melt differentiates. This melt probably represents that observed in pre- and postshield volcanoes. Finally, at ~70 km from the plume axis and at a depth greater than 200 km, the flowlines are subvertical. They then deflect at ~180 km depth and rotate toward the horizontal and eventually transit at 10–20° to the horizontal across an ~200 km distance from the axis and reach ~140 km depth. The fertilized-harzburgite and sub-primitive mantle/ultralow-velocity zone filaments that flow along these elbows partially melt by a few to several per cent. The resulting interstitial melt has a kimberlite-like composition. Thereafter, the excess pressure at the top of the filament at ~200 km from the axis overcomes the threshold for dyking and thus allows the escape of the interstitial melt via dykes ponding in subcrustal magma chambers or emerging directly at the surface. These melts have a composition similar to that associated with rejuvenated volcanism. We use the nature and the composition of whole erupted magmas and the seismic structure along the Hawaiian chain to validate this model.

Contrasting collision-induced far-field orogenesis controlled by thermo-rheological properties of the composite terrane

The continental collision-induced far-field intraplate deformation may be topographically separated from or coherent with the collision front, the mechanism of which is poorly constrained. Here, we build a series of 2D thermo–mechanical models to explore the controlling factors of distal orogenesis in the framework of a composite terrane sandwiched between two strong continental blocks. The model results suggest that the development of distal orogenesis requires a strong rheological heterogeneity in the overriding lithosphere. Two types of far-field collisional effects with distinct patterns of topographic evolution are identified, characterized by (I) a delayed and discrete distal orogenesis separated from the collision front by a low-relief, less-deformed terrane interior, and (II) an instant distal orogenesis subsequently integrated into a progressively growing plateau. These two types of far-field orogenesis are revealed to be analogous to the development of the Kopet Dagh and northern Tibet. The comparisons between model results and geological observations provide significant constraints on the dynamics of crustal-scale deformation propagation and the formation of distinct lithospheric architectures.

Oceanic Plateau and Subduction Zone Jump: Two‐Dimensional Thermo‐Mechanical Modeling

AbstractOceanic plateaus are characterized by overthickened crust and rheologically strong and chemically depleted mantle root with respect to normal oceanic lithosphere. The buoyant nature makes them have a great potential to resist subduction and initiate a new subduction zone behind the plateaus, yet the geodynamic process remains enigmatic. Here, we use 2D thermo‐mechanical modeling method to investigate the role of an oceanic plateau in subduction zone reorganization. We systematically examine the rheological strength and compositional density of oceanic plateau mantle and the size of oceanic plateau. Three distinct subduction styles are observed: (a) steep subduction without subduction zone reorganization, favored by a small‐sized (≤150 km) or lowly depleted oceanic plateau; (b) breakoff of subducting slab followed by flat‐slab subduction, which is conditional on a moderate‐sized (≥300 km) oceanic plateau with high degree depletion but relatively low strength; and (c) initiation of a new subduction zone behind the plateau (subduction zone jump), facilitated by a moderate‐sized oceanic plateau with a strongly depleted mantle root. The geodynamic results demonstrate that the overriding lithosphere undergoes a tectonic force switch from compression to tension in ∼5 Myr as the subduction zone is reorganized. This can be reflected in the crustal deformation of the overriding plate and used as an indicator to identify past subduction zone jumps. The spatial relationship predicted by the subduction zone jump models explains the Jurassic closure event of the accretion on southern Qiangtang at ∼175 Ma and the intraoceanic infant arc magmatism in the Meso‐Tethys at ∼165 Ma.

3He evidence for fluid transfer and continental hydration above a flat slab

The transfer of large volumes of fluid to the overriding lithosphere during flat-slab subduction should drastically alter the physical and chemical properties of continental margins. However, this process is poorly understood and without active magmatism, direct evidence for fluid transfer has remained elusive in modern systems. New helium isotope ratios (3He/4He) and gas species abundances in thermal springs above the Peruvian flat slab in South America demonstrate widespread transit of mantle-derived volatiles through thick (>45 km) continental crust. An unambiguous mantle signature is present in all 52 springs investigated over a 200,000 km2 region, with 3He/4He ratios ranging from 0.26 to 2.52 RC/RA (3He/4He ratio relative to the air ratio, RA, of 1.39×10−6, corrected for air-derived helium). In the absence of recent magmatism, the observed 3He/4He ratios are best explained by slab-derived fluids mobilizing helium from the subcontinental lithospheric mantle (SCLM). As slab- and mantle-derived volatiles migrate through the continental crust, 3He/4He ratios are reduced from SCLM values (6.1±2.1 RA) by ingrowth of radiogenic 4He and mixing with crustal fluids. The highest 3He/4He ratios in the flat slab region are observed above a probable tear in the Nazca Plate and along the Cordillera Blanca Detachment fault, implying that these features promote faster fluid transport through the crust. Our observations require a persistent flux of slab-derived fluids to move lithospheric-mantle-derived volatiles into and through the continental crust, providing evidence for active continental hydration in a modern flat-slab subduction system. Similar results from springs located in the backarc region to the south of the flat slab indicate that the mantle lithosphere below the northern Altiplano is experiencing dehydration and/or partial melting resulting from a mid-Cenozoic episode of flat-slab subduction.

Flat Subduction Versus Big Mantle Wedge: Contrasting Modes for Deep Hydration and Overriding Craton Modification

AbstractSubduction‐induced deep hydration and water cycling may strongly control the modification of overriding cratonic lithosphere. Two contrasting modes are generally proposed: (1) flat subduction (FS) regime with slab subducting subhorizontally beneath the overriding lithosphere and (2) big mantle wedge (BMW) regime with slab flattening at the bottom of mantle transition zone. Here, systematic numerical models are built to study the subduction‐induced deep hydration processes in the contrasting FS and BMW regimes as well as their effects on the modification of overriding lithosphere. The model results indicate that the dehydration process in the FS regime can significantly modify the overriding lithosphere for a region of about 600 km from the trench. During the progressive flat subduction, the partial melting and magmatism migrate toward the inner land of the overriding plate, which will be reversed and backward to the trench during the transition from flat to steep slab subduction. On the other hand, the deep hydration in the BMW regime is strongly dependent on the subcrustal serpentinite layer in the subducting slab, whereas the oceanic crust cannot carry water to the transition zone. The modification of the overriding lithosphere in the BMW regime occurs in a larger region of >1,000 km from trench, which is, however, generally slower and weaker. Thus, the flat subduction of Izanagi plate may play a more significant role in the modification of North China Craton in the late Jurassic to early Cretaceous, which could be accompanied by the effects of deep water cycling in the BMW regime.

The effects of the thermal state of overriding continental plate on subduction dynamics: Two-dimensional thermal-mechanical modeling

The dynamic process of ocean-continent subduction depends on not only the properties of the subducting oceanic plate, but also the characteristics and state of the overriding continental plate. Numerical models conducted to date have mostly focused on the oceanic lithosphere in this regard; research on the properties of overriding continental lithosphere remains relatively limited, especially the influence of its thermal state on subduction dynamics. Here we explored the performance of continental lithosphere with different thermal states during the subduction process using two-dimensional thermal-mechanical modeling and systematically investigated the effects of the thermal state of overriding continental plate, the age of subducting oceanic plate, and relative convergence rate on subduction dynamics. Modeling results show that: (1) When the geothermal gradient of continental crust is low (between 10 and 15 °C km−1), the oceanic plate first subducts at a low angle. As subduction continues, the slab dip gradually increases and the slab begins to retreat rapidly driven by its negative buoyancy, opening an ocean basin ranging from 600 to 1100 km in width. This leads to the decoupling between the overriding continental plate and oceanic plate. As the trench retreat continues, the horizontal deviatoric stress inside the overriding continental crust alternates between being positive and negative in a local area. Thinning of the overriding lithosphere mainly occurs at the region adjacent to the subduction zone, where the surface experiences significant subsidence. (2) When the geothermal gradient of continental crust is higher (greater than 15 °C km−1), oceanic plate retreat causes the overriding continental plate to be strongly stretched. In this case, the trench retreat distance decreases and the width of the ocean basin also reduces by between 100 and 1000 km. The horizontal deviatoric stress inside the whole overriding continental crust first manifests as compression and then changes into extension, which causes the surface to first uplift and then slowly subside. (3) Increasing the age of oceanic lithosphere accelerates trench retreat and promotes overriding plate thinning. (4) An advancing overriding continental plate slows down trench retreat. In cases where the geothermal gradient of continental crust is greater than 17.5 °C km−1, the hot continental crust experiences gravitational collapse and is overthrusted onto oceanic lithosphere, resulting in slow trench retreat. We analyzed the subduction process of the western Paleo-Pacific Plate in the Early Cretaceous based on our modeling results and discussed its possible control on the tectonic evolution of the rift basins in east Asia. We suggest that the development of a wide rift basin system on the Amurian Superterrane in the Early Cretaceous was likely related to slow trench retreat and the collapse of the hot crust, and the formation of a series of passive rift basins in the North China Craton was likely caused by the relatively cold thermal state of the lithosphere and the rapid retreat of the Paleo-Pacific Plate.

Influence of the Thickness of the Overriding Plate on Convergence Zone Dynamics.

The important role played by the upper plate in convergence zones dynamics has long been underestimated but is now more and more emphasized. However, the influence of its thickness and/or strength on orogenic systems evolution remains largely unknown. Here we present results from 3D thermo‐mechanical numerical simulations of convergence zones (including oceanic subduction followed by continental subduction/collision), in which we vary the rheological profile of the overriding plate (OP). For this, we systematically modify the crustal thickness of the overriding lithosphere and the temperature at the Moho to obtain a thermal thickness of the overriding lithosphere ranging from 80 to 180 km. While all models share a common global evolution (i.e., slab sinking, interaction between slab and the 660 km discontinuity, continental subduction/collision, and slab breakoff), they also highlight first‐order differences arising from the variations in the OP strength (thermal thickness). With a thin/weak OP, slab rollback is favored, the slab dip is low, the mantle flow above the slab is vigorous, and the trench migrates at a high rate compared to a thick/strong OP. In addition, slab breakoff and back‐arc basin formation events occur significantly earlier than in models involving a thick OP. Our models therefore highlight the major role played by the thickness/strength of the OP on convergence zone dynamics and illustrate its influence in a quantitative way.

Detailed tectonic reconstructions of the Western Mediterranean region for the last 35 Ma, insights on driving mechanisms

Slab retreat, slab tearing and interactions of slabs are first-order drivers of the deformation of the overriding lithosphere. An independent description of the tectonic evolution of the back-arc and peripheral regions is a pre-requisite to test the proposed conceptual, analogue and numerical models of these complex dynamics in 3-D. We propose here a new series of detailed kinematics and tectonic reconstructions from 35 Ma to the Present shedding light on the driving mechanisms of back-arc rifting in the Mediterranean where several back-arc basins all started to form in the Oligocene. The step-by-step backward reconstructions lead to an initial situation 35 Ma ago with two subduction zones with opposite direction, below the AlKaPeCa block (i.e.belonging to the Alboran, Kabylies, Peloritani, Calabrian internal zones). Extension directions are quite variable and extension rates in these basins are high compared to the Africa-Eurasia convergence velocity. The highest rates are found in the Western Mediterranean, the Liguro-Provençal, Alboran and Tyrrhenian basins. These reconstructions are based on shortening rates in the peripheral mountain belts, extension rates in the basins, paleomagnetic rotations, pressure-temperature-time paths of metamorphic complexes within the internal zones of orogens, and kinematics of the large bounding plates. Results allow visualizing the interactions between the Alps, Apennines, Pyrenean-Cantabrian belt, Betic Cordillera and Rif, as well as back-arc basins. These back-arc basins formed at the emplacement of mountain belts with superimposed volcanic arcs, thus with thick, hot and weak crusts explaining the formation of metamorphic core complexes and the exhumation of large portions of lower crustal domains during rifting. They emphasize the role of transfer faults zones accommodating differential rates of retreat above slab tears and their relations with magmatism. Several transfer zones are identified, separating four different kinematic domains, the largest one being the Catalan-Balearic-Sicily Transfer Zone. Their integration in the wider Mediterranean realm and a comparison of motion paths calculated in several kinematic frameworks with mantle fabric shows that fast slab retreat was the main driver of back-arc extension in this region and that large-scale convection was a subsidiary driver for the pre-8 Ma period, though it became dominant afterward. Slab retreat and back-arc extension was mostly NW-SE until ∼ 20 Ma and the docking of the AlKaPeCa continental blocks along the northern margin of Africa induced a slab detachment that propagated eastward and westward, thus inducing a change in the direction of extension from NW-SE to E-W. Fast slab retreat between 32 and 8 Ma and induced asthenospheric flow have prevented the transmission of the horizontal compression due to Africa-Eurasia convergence from Africa to Eurasia and favored instead upper-plate extension driven by slab retreat. Once slab retreat had slowed down in the Late Miocene, this N-S compression was felt and recorded again from the High Atlas to the Paris Basin.

Continental lithospheric-scale subduction versus crustal-scale underthrusting in the collision zone: Numerical modeling

Two distinct styles of continental collision are observed in nature, i.e. lithospheric-scale subduction versus crustal-scale underthrusting. The former refers to a collision scenario where the lower plate subducts below the lithospheric mantle of the overriding plate, whereas the latter is characterized by the underthrusting of the subducting plate beneath the crust of the overriding plate. Systematic numerical models were constructed to study the dynamics of these two different collision modes. The model results indicate that low convergence rate and cold overriding plate contribute to the lithospheric-scale subduction, which incorporates limited shortening of the overriding lithosphere during collision. In contrast, the crustal-scale underthrusting requires relatively high convergence rate and hot/weak overriding plate, in which the thickened lithospheric mantle of the overriding plate is progressively delaminated from the overlying crust and is replaced by the sub-horizontally underthrusted plate. The model results were systematically compared with four ongoing collisional orogens along the Alpine–Himalayan belt, i.e. the central-western Alps, the East Anatolian plateau, the Zagros–Iranian plateau, and the southern Tibetan plateau. It indicates that the lithospheric-scale subduction is widely observed and generally consistent with the natural orogens, for example, the Alps, Anatolia, and Zagros. On the other hand, the crustal-scale underthrusting is reconciled with the India–Asia collision zone, with respect to the geodynamic conditions, crustal-lithospheric structure, as well as the uplift history.

P Wave Azimuthal Anisotropic Tomography in Northern Chile: Insight Into Deformation in the Subduction Zone

AbstractBased on a large data set of local body wave travel times, we determined the first 3‐D model of azimuthal P wave anisotropic tomography of northern Chile in order to study the deformation in the subduction zone. Our results indicate different deformation patterns in the overriding lithosphere, mantle wedge, and the subducting Nazca slab as well as significant along‐arc variations. Radiating fast velocity directions (FVDs) around the rupture zone of the 2014 Iquique earthquake (Mw 8.2) are notable in the crust, which may reflect the specific crustal extension induced by a point‐like stress source in the form of a fully locked asperity at the plate interface. In most of the study region, FVDs in the mantle wedge are trench normal, which we interpret to have resulted from mantle wedge flow driven by the oceanic plate subduction. However, trench‐parallel FVDs are found beneath the northern segment, which may imply the occurrence of B‐type olivine fabrics in the cold forearc mantle. Strong along‐arc variations are also observed in the subducting slab, which may reflect intraslab heterogeneity and change in the slab geometry. In the northern segment, slab bending induces generally trench‐normal extension and leads to trench‐normal FVDs in the upper part of the slab. In contrast, in the southern segment, flat slab subduction and thus slab unbending dominate, producing trench‐normal compression and trench‐parallel FVDs instead. Along the subducting slab interface, FVDs change from trench‐parallel to trench‐normal in the downdip direction, which may indicate the first‐order transition from compression in the locking area to dominant shear in the decoupling area along the slab interface.

Fossilized Melts in Mantle Wedge Peridotites

The shallow oxidized asthenosphere may contain a small fraction of potassic silicate melts that are enriched in incompatible trace elements and volatiles. Here, to determine the chemical composition of such melt, we analysed fossilized melt inclusions, preserved as multiphase solid inclusions, from an orogenic garnet peridotite in the Bohemian Massif. Garnet-poor (2 vol.%) peridotite preserves inclusions of carbonated potassic silicate melt within Zn-poor chromite (<400 ppm) in the clinopyroxene-free harzburgite assemblage that equilibrated within the hot mantle wedge (Stage 1, > 1180 °C at 3 GPa). The carbonated potassic silicate melt, which has a major element oxide chemical composition of K2O = 5.2 wt.%, CaO = 17 wt.%, MgO = 18 wt.%, CO2 = 22 wt.%, and SiO2 = 20 wt.%, contains extremely high concentrations of large ion lithophile elements, similar to kimberlite melts. Peridotites cooled down to ≅800 °C during Stage 2, resulted in the growth of garnet relatively poor in pyrope content, molar Mg/(Mg + Fe + Ca + Mn), (ca. 67 mol.%). This garnet displays a sinusoidal REE pattern that formed in equilibrium with carbonatitic fluid. Subsequently, subduction of the peridotite resulted in the formation of garnet with a slightly higher pyrope content (70 mol.%) during the Variscan subduction Stage 3 (950 °C, 2.9 GPa). These data suggest the following scenario for the generation of melt in the mantle wedge. Primarily, infiltration of sediment-derived potassic carbonatite melt into the deep mantle wedge resulted in the growth of phlogopite and carbonate/diamond. Formation of volatile-bearing minerals lowered the density and strength of the peridotite. Finally, phlogopite-bearing carbonated peridotite rose as diapirs in the mantle wedge to form carbonated potassic silicate melts at the base of the overriding lithosphere.

Dynamics of thinning and destruction of the continental cratonic lithosphere: Numerical modeling

Thinning of the cratonic lithosphere is common in nature, but its destruction is not. In either case, the mechanisms for both thinning and destruction are still widely under debate. In this study, we have made a review on the processes and mechanisms of thinning and destruction of cratonic lithosphere according to previous studies of geological/geophysical observations and numerical simulations, with specific application to the North China Craton (NCC). Two main models are suggested for the thinning and destruction of the NCC, both of which are related to subduction of the oceanic lithosphere. One is the “bottom-up” model, in which the deeply subducting slab perturbs and induces upwelling from the hydrous mantle transition zone (MTZ). The upwelling produces mantle convection and erodes the bottom of the overriding lithosphere by the fluid-melt-peridotite reaction. Mineral compositions and rheological properties of the overriding lithospheric mantle are changed, allowing downward dripping of lithospheric components into the asthenosphere. Consequently, lithospheric thinning or even destruction occurs. The other is the “top-down” model, characterized by the flat subduction of oceanic slab beneath the overriding cratonic lithosphere. Dehydration reactions from the subducting slab would significantly hydrate the lithospheric mantle and decrease its rheological strength. Then the subduction angle may be changed from shallow to steep, inducing lateral upwelling of the asthenosphere. This upwelling would heat and weaken the overriding lithospheric mantle, which led to the weakened lithospheric mantle dripping into the asthenosphere. These two models have some similarities, in that both take the subducting oceanic slab and relevant fluid migration as the major driving mechanism for thinning or destruction of the overriding cratonic lithosphere. The key difference between the two models is the effective depth of the subducting oceanic slab. One is stagnation and flattening in the MTZ, whereas the other is flat subduction at the bottom of the cratonic lithosphere. In the NCC, the eastern lithosphere was likely affected by subduction of the Izanagi slab during the Mesozoic, which would have perturbed the asthenosphere and the MTZ, and induced fluid migration beneath the NCC lithosphere. The upwelling fluid may largely have controlled the reworking of the NCC lithosphere. In order to discuss and analyze these two models further, it is crucial to understand the role of fluids in the subduction zone and the MTZ. Here, we systematically discuss phase transformations of hydrous minerals and the transport processes of water in the subduction system. Furthermore, we analyze possible modes of fluid activity and the problems to explore the applied feasibility of each model. In order to achieve a comprehensive understanding of the mechanisms for thinning and destruction of cratonic lithosphere, we also consider four additional possible dynamic models: extension-induced lithospheric thinning, compression-induced lithospheric thickening and delamination, large-scale mantle convection and thermal erosion, and mantle plume erosion. Compared to the subduction-related models presented here, these four models are primarily controlled by the relatively simple and single process and mechanism (extension, compression, convection, and mantle plume, respectively), which could be the secondary driving mechanisms for the thinning and destruction of lithosphere.

Constraining quantitatively the timing and process of continent-continent collision using magmatic record: Method and examples

Based on the main driving force of plate motion (the slab pull force generated by the descent of the oceanic plate in subduction zones) and the three primary mechanisms for magma generation (adding fluid, increasing temperature, and decreasing pressure), the continent-continent collisional process has been divided into three stages, including initial collision, ongoing collision, and tectonic transition. These stages are characterized by normal calc-alkaline andesitic magma (dehydration of the oceanic crust to release fluids), the migration of calc-alkaline magma toward the trench (dehydration of the oceanic crust or an increase in temperature) or small-scale crust-derived peraluminous magma (heat from intra-crustal shearing), and extensive magmatism with compositional diversity induced by slab break-off (increasing temperature and decreasing pressure), respectively. On the basis of the obtained age of slab break-off, the timing of the initial continent-continent collision can be quantitatively back-dated using the convergence rate, depth of slab break-off, and subduction angle. The spatio-temporal migration of the magmatic activity of the Gangdese Batholith, the onset of magmatic flare-up, and the increase of magma temperature at 52–51 Ma documented by the volcanic rocks of the Linzizong Pana Formation were most likely the result of the break-off of the Yarlung-Zangbo Neo-Tethyan oceanic lithosphere at approximately 53 Ma. This proposed age of slab break-off suggests that the initial India-Asia collision likely occurred at approximately 55–54 Ma, which is consistent with the collision ages constrained by other abundant geological data (60–55 Ma). This magmatic method has been applied to the Bitlis orogenic belt in southern Turkey in the Arabia-Eurasia continental collision zone, yielding an age range of approximately 29–22 Ma for the initial Arabia-Asia continental collision that is close to the collision ages recently obtained by apatite fission-track dating (approximately 20 Ma) and regional tectonic shortening (approximately 27 Ma). The intense folding of the Upper Cretaceous and the angular unconformity between the overlying Linzizong volcanic rocks in the southern Lhasa Terrane (90−69 Ma) are not related to the initial continental collision between India and Asia, but can be interpreted as the consequences of the strong coupling between the hot and young subducting oceanic crust immediately south of the spreading ridge and the overriding lithosphere or the subduction of the Neo-Tethys oceanic plateaux or seamounts. The tectonic event documented by the angular unconformity between the Linzizong Dianzhong Formation and the Nianbo Formation lasted approximately 3 Ma and likely marks the initial India-Asia collision. The significant deceleration of the Indian continent at approximately 51 Ma can be attributed to the disappearance of the slab pull force in the subduction zone due to the break-off of the Yarlung-Zangbo Neo-Tethyan oceanic lithosphere. The descent of the eclogitized lower crust of the northern Indian continent provides the main driving force for the current northward motion of Indian plate. The weak deformation of the lithospheric plate in the overriding plate of the India-Asia collisional zone between 60 and 40 Ma can be attributed to the high-angle subduction related to the rollback of the Yarlung-Zangbo Neo-Tethyan oceanic lithosphere after the initial India-Asia continental collision, the presence of the thick crust and high elevation on the southern margin of the Lhasa Terrane, and the decoupling between the mid-upper and lower crust and between the lower crust and lithospheric mantle of the Indian continent.