Search

The formation of mountain belts or rift zones is commonly attributed to interactions between plates along their boundaries, but the widely distributed deformation of Asia from Himalaya to the Japan Sea and other back‐arc basins is difficult to reconcile with this notion. Through comparison of the tectonic and kinematic records of the last 50 Ma with seismic tomography and anisotropy models, we show that the closure of the former Tethys Ocean and the extensional deformation of East Asia can be best explained if the asthenospheric mantle transporting India northward, forming the Himalaya and the Tibetan Plateau, reaches East Asia where it overrides the westward flowing Pacific mantle and contributes to subduction dynamics, distributing extensional deformation over a 3,000‐km wide region. This deep asthenospheric flow partly controls the compressional stresses transmitted through the continent‐continent collision, driving crustal thickening below the Himalayas and Tibet and the propagation of strike‐slip faults across Asian lithosphere further north and east, as well as with the lithospheric and crustal flow powered by slab retreat east of the collision zone below East and SE Asia. The main shortening direction in the deforming continent between the collision zone and the Pacific subduction zones may in this case be a proxy for the direction of flow in the asthenosphere underneath, which may become a useful tool for studying mantle flow in the distant past. Our model of the India‐Asia collision emphasizes the role of asthenospheric flow underneath continents and may offer alternative ways of understanding tectonic processes., Key Points Coeval India‐Asia collision and Pacific back‐arc extension explained by asthenospheric flowMantle flow correlates with long‐term kinematics from India to Pacific borderIndian mantle reaches Pacific border and Asian lithosphere spreads above ponding slabs

Tectonics and erosion–sedimentation are the main processes responsible for shaping the Earth's surface. The link between these processes has a strong influence on the evolution of landscapes. One of the tools we have for investigating coupled process models is analog modeling. Here we contribute to the utility of this tool by presenting laboratory-scaled analog models of erosion. We explore the erosional response of different materials to imposed boundary conditions, trying to find the composite material that best mimics the behavior of the natural prototype. The models recreate conditions in which tectonic uplift is no longer active, but there is an imposed fixed slope. On this slope the erosion is triggered by precipitation and gravity, with the formation of channels in valleys and diffusion on hillslope that are functions of the analog material. Using digital elevation models (DEMs) and a laser scan correlation technique, we show model evolution and measure sediment discharge rates. We propose three main components of our analog material (silica powder, glass microbeads and PVC powder; PVC: polyvinyl chloride), and we investigate how different proportions of these components affect the model evolution and the development of landscapes. We find that silica powder is mainly responsible for creating a realistic landscape in the laboratory. Furthermore, we find that varying the concentration of silica powder between 40 wt % and 50 wt % (with glass microbeads and PVC powder in the range 35 wt %–40 wt % and 15 wt %–20 wt %, respectively) results in metrics and morphologies that are comparable with those from natural prototypes., Earth Surface Dynamics, 8 (4), ISSN:2196-632X, ISSN:2196-6311

The study of slab-slab interactions has come to the front of geodynamics researches to explain geological and geophysical observations from tectonically complex areas. Here we aim to better understand the geodynamics of the Central Mediterranean, where the Adria plate subducts on its two opposite sides. Additionally, the slab below the Central South Apennines has been progressively breaking off during the last 3 Myr. The role of a slab window in a single slab or in an outward dipping double-sided subduction system is addressed by analog models at the scale of the upper mantle, realized using glucose syrup and silicone putty, to model the interaction between the Earth's mantle and the lithosphere. Our results show that the presence of a slab window modifies the pattern of mantle circulation, as well as the trench geometry and kinematics. In particular, the opening of the slab window induces the formation of two arcs flanking the window, while the mantle flows through it and turns toward the arcs, creating a small-scale toroidal flow. The effect of a slab window is more pronounced on double subduction systems, as the outflow through the window is amplified, while internal deformation is induced in the plate by the opposite slab pull force. These experimental results suggest that the origin of the Apenninic and the Calabrian arcs is the result of the formation of a slab window, providing a new interpretation of the surface deformation and the SKS shear wave splitting pattern of the Adria microplate.Plain Language Summary The Adria microplate in the Central Mediterranean is an important part of the tectonic puzzle between the converging Africa and Eurasia plates. The heavy lithosphere sinks into the underlying mantle on two sides of the Adria plate, under the Apennines in the west and under the Dinarides in the east. During the most recent few millions of years, the lithosphere started to detach from the surface in the Central South Apennines while sinking into the mantle. This opened a slab window, an additional escape route for the mantle, which is squeezed by the two subducting lithospheres. We scaled down and simplified this setting in order to model and investigate the role of the dual subduction and the opening slab window in the evolution of the Central Mediterranean by means of analog models. Our results show that the two opposite subduction zones have a regional tectonics effect via the complex subduction-induced mantle circulation. Furthermore, the detaching lithosphere modifies the geometry of the subduction zone, which results in a double-arc geometry resembling the Northern Apennines and Calabria.

The importance of slab–slab interactions is manifested in the kinematics and geometry of the Philippine Sea Plate and western Pacific subduction zones, and such interactions offer a dynamic basis for the first-order observations in this complex subduction setting. The westward subduction of the Pacific Sea Plate changes, along-strike, from single slab subduction beneath Japan, to a double-subduction setting where Pacific subduction beneath the Philippine Sea Plate occurs in tandem with westward subduction of the Philippine Sea Plate beneath Eurasia. Our 3-D numerical models show that there are fundamental differences between single slab systems and double slab systems where both subduction systems have the same vergence. We find that the observed kinematics and slab geometry of the Pacific–Philippine subduction can be understood by considering an along-strike transition from single to double subduction, and is largely independent from the detailed geometry of the Philippine Sea Plate. Important first order features include the relatively shallow slab dip, retreating/stationary trenches, and rapid subduction for single slab systems (Pacific Plate subducting under Japan), and front slabs within a double slab system (Philippine Sea Plate subducting at Ryukyu). In contrast, steep to overturned slab dips, advancing trench motion, and slower subduction occurs for rear slabs in a double slab setting (Pacific subducting at the Izu–Bonin–Mariana). This happens because of a relative build-up of pressure in the asthenosphere beneath the Philippine Sea Plate, where the asthenosphere is constrained between the converging Ryukyu and Izu–Bonin–Mariana slabs. When weak back-arc regions are included, slab–slab convergence rates slow and the middle (Philippine) plate extends, which leads to reduced pressure build up and reduced slab–slab coupling. Models without back-arcs, or with back-arc viscosities that are reduced by a factor of five, produce kinematics compatible with present-day observations.

Subducted sediments are thought to lubricate the subduction interface and promote faster plate speeds. However, global observations are not clear-cut on the relationship between the amount of sediments and plate motion. Sediments are also thought to influence slab dip, but variations in subduction geometry depend on multiple factors. Here we use 2D thermomechanical models to explore how sediments can influence subduction dynamics and geometry. We find that thick sediments can lead to slower subduction due to an increase of the megathrust shear stress as the accretionary wedge gets wider, and a decrease in slab pull as buoyant sediments are subducted. Our results also show that larger slab buoyancy and megathrust stress due to thick sediments increase the slab bending radius. This offers a new perspective on the role of sediments, suggesting that sediment buoyancy and wedge geometry also play an important role on large-scale subduction dynamics.

The central Apennines, an accretionary wedge overlying an area of slab detachment, are characterized by prominent topography, active normal faulting, and high uplift rates. However, previous studies have failed to resolve the surface uplift history, complicating efforts to link the topographic evolution with underlying geodynamic processes. We aim to better quantify orographic changes by using stable oxygen isotope paleoaltimetry. Modern surface water δ 18 O are 5‰ lower at high elevation than at sea level, reflecting orographic rainout over the Apennines. We present 262 new lacustrine and paleosol carbonate δ 18 O measurements collected from ten extensional intermontane basins—spanning both high and low elevations—and combine these with 1,166 published δ 18 O data, permitting us to constrain changes in δ 18 O both spatially and temporally. Since the Pliocene, δ 18 O in present-day high-elevation basins has continuously decreased, even as δ 18 O in lowland basins has remained constant over time. We attribute this continuous 5‰ shift to increased orographic rainout as the central Apennines were uplifted. We estimate an increase in mean elevation of approximately 1–2 km since the late Pliocene, and these estimates match the suggested timing and expected amplitude of slab break-off related uplift. This supports the hypothesis that the opening of the Adriatic slab window and associated mantle flow contributed significantly to building topography in the central Apennines.

A fundamental scientific question is, what controls the Earth's topography? Although the theoretical principles of isostasy, flexure, and dynamic topography are widely discussed, the parameters needed to apply these principles are frequently not available. Isostatic factors controlling lithospheric buoyancy are frequently uncertain and non-isostatic factors, such as lithospheric bending towards subduction zones and dynamic topography, are hard to distinguish. The question discussed here is whether a set of simple rules that relate topography to lithospheric structure in various tectonic environments can be deduced in a way that missing parameters can be approximated; or does each area behave differently, making generalizations problematic. We contribute to this issue analyzing the Asia–Africa–Arabia–Europe domain following a top-down strategy. We compile a new crustal thickness map and remove the contribution of the crust from the observed elevation. Then, the challenge is to interpret the residual topography in terms of mantle lithosphere buoyancy and dynamics. Based on systematic relationships between tectonic environments and factors controlling topography, we argue that crustal buoyancy and mantle lithospheric density can be approximated from available geological data and that regions near mantle upwelling or downwelling are easily identified by their extreme residual topography. Yet, even for other areas, calculating lithospheric thickness from residual topography is problematic, because distinguishing variations in mantle lithosphere thickness from sub-lithospheric dynamics is difficult. Fortunately, the area studied here provides an opportunity to examine this issue. Based on the conjunction between the Afar Plume and the mid-ocean ridge in the nearby Gulf of Aden and southern Red Sea, we constrain the maximal amplitude of dynamic topography to ~ 1 km. This estimate is based on a narrow definition of dynamic topography that only includes sub-lithospheric processes and using mid-ocean ridges as reference, where mantle lithosphere buoyancy is zero.

Seismic evidence and thermal and topographic transients have led to the interpretation of lithospheric removal beneath the Southeast Carpathians region. A series of numerical geodynamic experiments in the context of the tectonic evolution of the region are conducted to test the surface-crustal response to lithosphere delamination and slab break-off. The results show that a delamination-type removal (“plate-like” migrating instability) causes a characteristic pattern of surface uplift/subsidence and crustal extension/shortening to occur due to the lithospheric deformation and dynamic/thermal forcing of the sublithospheric mantle. These features migrate with the progressive removal of the underlying lithosphere. Model results for delamination are comparable with observables related to the geodynamic evolution of the Southeast Carpathians since 10 Ma: the mantle structure inferred by seismic tomography, migrating patterns of uplift (>1.5 km) and subsidence (>2 km) in the region, crustal thinning in the Carpathian hinterland and thickening at the Focsani depression, and regional extension in the Carpathian corner (e.g., opening of Brasov basin) correlating with volcanism (e.g., Harghita and Persani volcanics) in the last 3 Myr.

International audience; Unravelling the evolution of supradetachment basins developed in the hanging wall of low‐angle detachment faults may be an invaluable tool in reconstructing the tectonic evolution of highly extended terrains. These basins may record major regional tectonic events related to the exhumation of metamorphic core complexes, and the reconstruction of their evolution helps to quantify the amount of extension accommodated by such processes. Here we present stratigraphic and structural field evidence and micropaleontological constraints to the Neogene‐to‐Quaternary evolution of the supradetachment Gediz Graben that developed on top of the exhuming Central Menderes Massif (SW Turkey). This basin displays three different structural styles during its evolution: (i) it initiated as a ramp basin following the activation of the Gediz Detachment in the Middle Miocene; (ii) evolved as a half graben during the Late Miocene following the activation of high‐angle brittle faults at its southern margin; and (iii) reached its final symmetric graben configuration in Late Pliocene (?) – Quaternary times following the activation of its northern margin. New micropaleontological data document a short‐lived upper Tortonian marine episode in the basin, and major along‐strike variations in exhumation are documented on its southern margin. Our reconstruction shows how sedimentary basins originally formed in the hanging‐wall of detachment faults may eventually end up in tectonic contact with the mylonitic footwall. Finally, we highlight the importance of magmatism in localizing the deformation in highly extended terrains and in controlling the evolution of supradetachment systems.

International audience; Trench motions represent the surface expression of the interaction between subducting plates and the underlyingmantle, but the inherent dynamics are not fully understood. One interesting case is the migration of the Izu-Bonin Marianas trench (IBM) that accommodates the subduction of the Pacific beneath the Philippine Sea Plate (PSP), which is in turn subducting beneath the Eurasian plate along the Ryukyu trench. The history of the IBM trench is dominated by fast, episodic retreat from 40 to 15 Ma. However, around 10–5 Ma, the IBM trenchreversed its motion from retreat to advance. The switch in trench motion occurred soon after the breakoff of thePSP slab along the Ryukyu trench and the onset of new subduction, and represents a fundamental change in thedynamics of the western Pacific subduction zones. Here, Following the modelling study of Čížková and Bina(2015), which suggested a link between IBM trench advance and Ryukyu subduction, we run 2-D numericalexperiments to test the influence of a newly formed Ryukyu slab on the established IBM subduction zone, we runtwo-dimensional numerical experiments to test the influence of this newly formed Ryukyu slab on the IBMsubduction zone. The results from our geodynamic model compare favourably with the reconstructed trend, indicating that the switch in trench motion along the IBM trench may indeed be related to the onset of a new subduction zone along the Ryukyu trench. Our analysis substantiates the idea that advancing trench motions in the western Pacific are due to the establishment of a double subduction system. Further analysis of such dynamics provides insights for the mechanisms controlling subducting plate and trench motions and mantle force transmission.

Active and fossil subduction systems consisting of two adjacent plates with opposite retreating directions occur in several areas on Earth, as the Mediterranean or Western Pacific. The goal of this work is to better understand the first-order plate dynamics of these systems using the results of experimental models. The laboratory model is composed of two separate plates made of silicon putty representing the lithosphere, on top of a tank filled with glucose syrup representing the mantle. The set of experiments is designed to test the influence of the width of plates and the initial separation between them on the resulting trench velocities, deformation of plates, and mantle flow. Results show that the mantle flow induced by both plates is asymmetric relative to the axis of each plate causing a progressive merging of the toroidal cells that prevents a steady state phase of the subduction process and generates a net outward drag perpendicular to the plates. Trench velocities increase when trenches approach each other and decrease when they separate after their intersection. The trench curvature of both plates increases linearly with time during the entire evolution of the process regardless their width and initial separation. The interaction between the return flows associated with each retreating plate, particularly in the interplate region, is stronger for near plate configurations and correlates with variations of rollback velocities. We propose that the inferred first-order dynamics of the presented analog models can provide relevant clues to understand natural complex subduction systems. ©2018. American Geophysical Union. All Rights Reserved., The work of M. P., M. F., and J. V. is part of the projects ALPIMED (PIE-CSIC- 201530E082) and MITE (CGL2014- 59516-P) and partly financed by the Generalitat de Catalunya (2014-SGR-1471 and 2017-SGR-847). The work of S. Z. is funded by the Spanish Ministry of Economy and Competitiveness through grant CICYT DPI2017-85139-C2-2-R and Generalitat de Catalunya (2017-SGR-1278).

Since several decades, the processes allowing for the subduction of the continental lithosphere less dense than the mantle in a collision context have been widely explored, but models that are based upon the premise that slab pull is the prominent driver of plate tectonics fail. The India–Asia collision, where several episodes of continental subduction have been documented, constitute a case study for alternative views. One of these episodes occurred in the early collision time within the Asian plate where continental lithosphere not attached to any oceanic lithosphere subducted southward in front of the Indian lithosphere during its northward subduction that followed the oceanic subduction of the Tethys ocean. This process, known as collisional subduction, has a counter-intuitive behavior since the subduction is not driven by slab pull. It has been speculated that the mantle circulation can play an important role in triggering collisional subduction but a detailed, qualitative analysis of it is not available, yet. In this work we explore the southward subduction dynamics of the Asian lithosphere below Tibet by means of analogue experiments with the aim to highlight how the mantle circulation induces or responds to collisional subduction. We found that during the northward oceanic subduction (analogue of Tethys subduction) attached to the indenter (Indian analogue), the main component of slab motion is driven vertically by its negative buoyancy, while the trench rolls back. In the mantle the convective pattern consists in a pair of wide convective cells on both sides of the slab. But when the indenter starts to bend and plunge in the mantle, trench motion reverses. Its advance transmits the far field forces to two upper plates (Asian analogues). The more viscous frontal plate thickens, and the less viscous hinterland plate, which is attached to the back wall of the box, subducts. During this transition, a pair of sub-lithospheric convective cells is observed on both sides of the Asian analogue slab, driven by the shortening of the frontal plate. It favors the initiation of the backwall plate subduction. Such subduction is maintained during the entire collision by a wide cell with a mostly horizontal mantle flow below Tibet, passively advecting the Asian analogue slab. Experimental results suggest that once the tectonic far-field force related to the forward horizontal motion becomes dominant upon the buoyancy forces, trench advancing and the transmission of the tectonic force to the upper and backwall plates are promoted. This peculiar condition triggers the subduction of the backwall plate, despite it is light and buoyant.

Slab‐slab interaction is a characteristic feature of tectonically complex areas. Outward dipping double‐sided subduction is one of these complex cases, which has several examples on Earth, most notably the Molucca Sea and Adriatic Sea. This study focuses on developing a framework for linking plate kinematics and slab interactions in an outward dipping subduction geometry. We used analog and numerical models to better understand the underlying subduction dynamics. Compared to a single subduction model, double‐sided subduction exhibits more time‐dependent and vigorous toroidal flow cells that are elongated (i.e., not circular). Because both the Molucca and Adriatic Sea exhibit an asymmetric subduction configuration, we also examine the role that asymmetry plays in the dynamics of outward dipping double‐sided subduction. We introduce asymmetry in two ways; with variable initial depths for the two slabs (“geometric” asymmetry), and with variable buoyancy within the subducting plate (“mechanical” asymmetry). Relative to the symmetric case, we probe how asymmetry affects the overall slab kinematics, whether asymmetric behavior intensifies or equilibrates as subduction proceeds. While initial geometric asymmetry disappears once the slabs are anchored to the 660 km discontinuity, the mechanical asymmetry can cause more permanent differences between the two subduction zones. In the most extreme case, the partly continental slab stops subducting due to the unequal slab pull force. The results show that the slab‐slab interaction is most effective when the two trenches are closer than 10–8 cm in the laboratory, which is 600–480 km when scaled to the Earth.

The elevation of an orogenic belt is commonly related to crustal/lithosphere thickening. Here, we discuss the Apennines as an example to show that topography at a plate margin may be controlled not only by isostatic adjustment but also by dynamic, mantle-driven processes. Using recent structural constraints for the crust and mantle we find that the expected crustal isostatic component explains only a fraction of the topography of the belt, indicating positive residual topography in the central Apennines and negative residual topography in the northern Apennines and Calabria. The trend of the residual topography matches the mantle flow induced dynamic topography estimated from regional tomography models. We infer that a large fraction of the Apennines topography is related to mantle dynamics, producing relative upwellings in the central Apennines and downwellings in the northern Apennines and Calabria where subduction is still ongoing. Comparison between geodetic and geological data on vertical motions indicates that this dynamic process started in the early Pleistocene and the resulting uplift appears related to the formation and enlargement of a slab window below the central Apennines. The case of the Apennines shows that at convergent margins the elevation of a mountain belt may be significantly different from that predicted solely by crustal isostasy and that a large fraction of the elevation and its rate of change are dynamically controlled by mantle convection.

Isostatic and dynamic models of Earth's surface topography can provide important insights into the driving processes of tectonic deformation. We analyze these two estimates for the tectonically-active western United States using refined structural models derived from EarthScope USArray. For the crust, use of recent Moho depth measurements and crustal density anomalies inferred from passive source seismology improve isostatic models. However, seismically determined lithospheric thickness variations from “lithosphere–asthenosphere boundary” (LAB) maps, and lithospheric and mantle density anomalies derived from heat flow or uppermost mantle tomography, do not improve isostatic models substantially. Perhaps this is a consequence of compositional heterogeneity, a mismatch between thermal and seismological LAB, and structural complexity caused by smaller-scale dynamics. The remaining, non-isostatic (“dynamic”) component of topography is large. Topography anomalies include negative residuals likely due to active subduction of the Juan de Fuca plate, and perhaps remnants of formerly active convergence further south along the margin. Our finding of broad-scale, positive residual topography in the Basin and Range substantiates previous results, implying the presence of anomalous buoyancy there which we cannot fully explain. The Colorado Plateau does not appear dynamically anomalous at present, except at its edges. Many of the residual topography features are consistent with predictions from mantle flow computations. This suggests a convective origin, and important interactions between vigorous upper mantle convection and intraplate deformation.

The Middle East region represents a key site within the Tethyan domain where continental break-up, collision, backarc extension and escape tectonics are kinematically linked together. We perform global mantle circulation computations to test the role of slab pull and mantle upwellings as driving forces for the kinematics of the Arabia-Anatolia-Aegean (AAA) system, evaluating different boundary conditions and mantle density distributions as inferred from seismic tomography or slab models. Model results are compared with geodetically inferred crustal motions, residual topography, and shear wave splitting measurements. The AAA velocity field with respect to Eurasia shows an anti-clockwise toroidal pattern, with increasing velocities toward the Aegean trench. The best match to these crustal motions can be obtained by combining the effect of slab pull exerted in the Aegean with a mantle upwelling underneath Afar and, more generally, with the large-scale flow associated with a whole mantle, Tethyan convection cell. Neogene volcanism for AAA is widespread, not only in the extensional or subduction settings, but also within plates, such as in Syria-Jordan-Israel and in Turkey, with geochemical fingerprints similar of those of the Afar lava. In addition, morphological features show large uplifting domains far from plate boundaries. We speculate that the tectonic evolution of AAA is related to the progressive northward entrainment of upwelling mantle material, which is itself associated with the establishment of the downwelling part of a convection cell through the segmented Tethyan slab below the northern Zagros and Bitlis collision zone. The recently established westward flow dragged Anatolia and pushed the Aegean slab south-westward, thus accelerating backarc extension. Our model reconciles Afar plume volcanism, the collision in the Bitlis mountains and northern Zagros, and the rapid increase of Aegean trench rollback in a single coherent frame of large scale mantle convection, initiated during the last similar to 40 Ma. (C) 2013 Elsevier B.V. All rights reserved.

International audience; We review the geodynamic evolution of the Aegean-Anatolia region and discuss strain localisation there over geological times. From Late Eocene to Present, crustal deformation in the Aegean backarc has localised progressively during slab retreat. Extension started with the formation of the Rhodope Metamorphic Core Complex (Eocene) and migrated to the Cyclades and the northern Menderes Massif (Oligocene and Miocene), accommodated by crustal-scale detachments and a first series of core complexes (MCCs). Extension then localised in Western Turkey, the Corinth Rift and the external Hellenic arc after Messinian times, while the North Anatolian Fault penetrated the Aegean Sea. Through time the direction and style of extension have not changed significantly except in terms of localisation. The contributions of progressive slab retreat and tearing, basal drag, extrusion tectonics and tectonic inheritance are discussed and we favour a model (1) where slab retreat is the main driving engine, (2) successive slab tearing episodes are the main causes of this stepwise strain localisation and (3) the inherited heterogeneity of the crust is a major factor for localising detachments. The continental crust has an inherited strong heterogeneity and crustal-scale contacts such as major thrust planes act as weak zones or as zones of contrast of resistance and viscosity that can localise later deformation. The dynamics of slabs at depth and the asthenospheric flow due to slab retreat also have influence strain localisation in the upper plate. Successive slab ruptures from the Middle Miocene to the Late Miocene have isolated a narrow strip of lithosphere, still attached to the African lithosphere below Crete. The formation of the North Anatolian Fault is partly a consequence of this evolution. The extrusion of Anatolia and the Aegean extension are partly driven from below (asthenospheric flow) and from above (extrusion of a lid of rigid crust).

The Cordillera of the Andes is a double-vergent orogenic belt built up by thickening of South American plate crust. Several models provide plausible explanations for the evolution of the Andes, but the reason why shortening started at ∼50 Ma is still unclear. We explore the evolution of the subduction zone through time by restoring the position of the Nazca trench in an absolute reference frame, comparing its position with seismic tomography models and balancing the evolution of the subducting slab. Reconstructions show that the slab enters into the lower mantle at ∼ 50 ± 10 Ma , and then progressed, moving horizontally at shallow lower mantle depth while thickening and folding in the transition zone. We test this evolutionary scenario by numerical models, which illustrate that compression in the upper plate intensifies once the slab is anchored in the lower mantle. We conclude that onset of significant shortening and crustal thickening in the Andes and its sustained action over tens of million years is related to the penetration of the slab into the lower mantle, producing a slowdown of lateral slab migration, and dragging the upper plate against the subduction zone by large-scale return flow.

Density anomalies beneath the lithosphere are expected to generate dynamic topography at the Earth's surface due to the induced mantle flow stresses which scale linearly with density anomalies, while the viscosity of the upper mantle is expected to control uplift rates. However, limited attention has been given to the role of the lithosphere. Here we present results from analogue modeling of the interactions between a density anomaly rising in the mantle and the lithosphere in a Newtonian system. We find that, for instabilities with wavelengths of the same order of magnitude as lithosphere thickness, the uplift rate and the geometry of the surface bulge are inversely correlated to the lithosphere thickness. We also show that a layered lithosphere may modulate the topographic signal. With respect to previous approaches our models represent a novel attempt to unravel the way normal stresses generated by mantle flow are transmitted through a rheologically stratified lithosphere and the resulting topographic signal. © 2017. American Geophysical Union. All Rights Reserved., T.W.B. was partially supported through NASA OSP 201601412-001.

Subduction-induced mantle circulation plays an important role in the dynamics of convergent margins. Different components of the flow, i.e. toroidal and poloidal, provide relevant driving forces for back-arc basin formation, overriding plate deformation, curvature of subduction zones and volcanic activity. Here, we investigate on the emergence and controls on the toroidal component of the subduction–induced mantle flow by means of numerical modeling. To characterize the toroidal cell's three-dimensional flow, size and length-scales and its disposing factors, we test separately a series of lithospheric and mantle parameters, such as the density difference and viscosity ratio between the slab and the mantle, the width of the slab, as opposed to the size, the stratification and the rheology of the mantle. Out of the tested parameters, the numerical results show that the strength of the flow depends on the mantle viscosity and the magnitude of the slab pull force, that is slab-mantle density difference and the mantle thickness, however the characteristic length, axis and the shape of the toroidal cell are almost independent of the slab's properties and mainly depend on the thickness of the convecting mantle.

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.

Exhumation of high-pressure rocks has long remained a controversial issue in the Earth sciences. In this article, we analyze the tectono-metamorphic, stratigraphic and plate-motion constraints from the Western Alps region, providing new insights on exhumation mechanisms and tectonic evolution during the earliest orogenic stages. Eocene eclogites of the Western Alps form a 20–25 km wide belt on the upper-plate side of the orogen (Eclogite belt), exposed beneath extensional shear zones at the rear of a lower-pressure accretionary wedge. Units of the Eclogite belt show the youngest peak-pressure assemblages within the subduction zone, and experienced superfast tectonic exhumation since 45–40 Ma. The role of erosion was negligible during the whole of this stage. Eocene foreland basins remained starved, and the massive arrival of axial-belt detritus began well after exhumation was completed. Tectonic reconstructions based on fixed-boundaries exhumation models (e.g. channel flow), and/or implying fast erosion at the surface (e.g. slab breakoff), are thus not consistent with geological evidence. In the lack of erosion, exhumation through the overburden requires divergence within the subduction zone. We demonstrate that this was not attained by rollback of the lower plate (Europe), but was instead attained by NNEward motion of the upper plate (Adria-Africa) away from the Western Alps trench. Such motion induced localized extension within the weak portion of the upper plate, at the rear of the accretionary wedge, and allowed tectonic emplacement of the Eclogite belt in the upper crust at rates much faster than subduction rates. Tectonic exhumation ceased in the Oligocene, when oblique-divergence along the Western Alps traverse changed into oblique-convergence. The onset of slow erosional unroofing was synchronously recorded by pressure–temperature paths in all major tectonic units of the Western Alps, and by arrival of orogenic detritus in sedimentary basins. This work demonstrates that divergence between upper plate and trench is a viable mechanism to exhume large and coherent eclogite units in continental subduction zones. Our exhumation model can be applied to other eclogite belts showing a similar exhumational record, including the Western Gneiss Region, the Dabie-Sulu, and eastern Papua New Guinea.

International audience; The tectonic and topographic evolution of southeast Asia is attributed to the indentation of India into Eurasia, gravitational collapse of the uplifted terrains and the dynamics of the Sunda and other western Pacific subduction zones, but their relative contributions remain elusive. Here, we analyse 3D numerical geodynamic modeling results involving a collision-subduction system and show that vigorous asthenospheric flow due to differential along-strike slab kinematics may contribute to the surface strain and elevations at collision-subduction transition zones. We argue that protracted northward migration of the collisional front and Indian slab during south to south-westward rollback subduction along the Sunda margin might have produced a similar asthenospheric flow. This flow could have contributed to the southeast Asia extrusion tectonics and uplift of the terrains around the eastern Himalayan syntaxis and protruding from southeast Tibet. Therefore, we suggest that the tectonics and topographic growth east and southeast of Tibet are controlled not only by crustal and lithospheric deformation but also by asthenospheric dynamics.

International audience; Since the Mesozoic, Africa has been under extension with shorter periods of compression associated with obductionof ophiolites on its northern margin. Less frequent than “normal” subduction, obduction is a first order process that remainsenigmatic. The closure of the Neo-Tethys Ocean, by the Upper Cretaceous, is characterized by a major obduction event, from theMediterranean region to the Himalayas, best represented around the Arabian Plate, from Cyprus to Oman. These ophiolites wereall emplaced in a short time window in the Late Cretaceous, from 100 to 75 Ma, on the northern margin of Africa, in a contextof compression over large parts of Africa and Europe, across the convergence zone. The scale of this process requires anexplanation at the scale of several thousands of kilometres along strike, thus probably involving a large part of the convectingmantle. We suggest that alternating extension and compression in Africa could be explained by switching convection regimes.The extensional situation would correspond to steady-state whole-mantle convection, Africa being carried northward by alarge-scale conveyor belt, while compression and obduction would occur when the African slab penetrates the upper–lowermantle transition zone and the African plate accelerates due to increasing plume activity, until full penetration of the Tethys slabin the lower mantle across the 660 km transition zone during a 25 Myr long period. The long-term geological archives on whichsuch scenarios are founded can provide independent time constraints for testing numerical models of mantle convection andslab–plume interactions.; Depuis le Mésozoïque, l’Afrique a été en extension avec de courtes périodes de compression associées a`l’obduction d’ophiolites sur sa marge nord. Moins fréquente que la subduction, l’obduction est néanmoins un phénomènede premier ordre qui reste énigmatique. La fermeture de la Neo-Téthys au Crétacé supérieur est caractérisée par un épisodemajeur d’obduction, depuis la Méditerranée jusqu’a` l’Himalaya, en particulier sur la marge de l’Arabie, de Chypre a` l’Oman.Ces ophiolites furent toutes mises en place dans un court laps de temps pendant le Crétacé supérieur, de 100 a` 75 Ma, dansun contexte de compression enregistré en de larges portions de l’Afrique et de l’Europe, au travers de la zone de convergence.L’échelle de ce processus requiert une explication a` l’échelle de plusieurs milliers de kilomètres et donc impliquantvraisemblablement l’ensemble du manteau convectif. Nous suggérons que l’alternance de périodes extensives et compressivesen Afrique résulte de changements du régime convectif. Les périodes extensives correspondraient a` la convectionimpliquant tout le manteau, l’Afrique étant portée par une grande cellule de type tapis-roulant, tandis que la compressionet l’obduction se produiraient quand le panneau plongeant africain pénètre la transition entre le manteau supérieur et lemanteau inférieur et quand la plaque Afrique accélère en conséquence d’une plus grande activité du panache, jusqu’a`pénétration complète, durant une période d’environ 25 Myr. Les archives géologiques sur lesquelles ce type de scenarios estfondé peuvent fournir des contraintes temporelles indépendantes pour tester les modèles numériques de convectionmantellique et les interactions panache–panneau plongeant.

Using silicone slabs as a model analogue for lithospheric plates subducting into a box of glucose syrup, as an analogue of the mantle, we explore the subduction of continental lithosphere in a context of intercontinental collision. The continental indenter pushed by a piston, reproducing the collision, attached to a dense oceanic plate, subducts to two-thirds of the depth of the mantle box. We show that, surprisingly, the continental plate attached to the back wall of the box subducts, even if not attached to a dense oceanic slab. The engine of this subduction is not the weight of the slab, because the slab is lighter than the mantle, but the motion of the piston, which generates horizontal tectonic forces. These are transmitted to the back wall plate through the indenter and the upper plate at the surface, and by the advancing indenter slab through the mantle at shallow depth. We define this process as collisional subduction occurring in a compressional context. The collisional subduction absorbs between 14% and 20% of the convergence, and represents an unexplored component of collisional mass balance. The transmission of tectonic forces far from the collision front favors the formation of a wide plateau. Our experiments reproduce adequately the amount and geometry of the Asian lithosphere subduction episodes inferred during the collision, leading us to conclude that it reproduces adequately the physics of such process.

The combination of analytical and analogue models gives new opportunities to better understand the kinematic parameters controlling the evolution of transpression zones. In this work, we carried out a set of analogue models using the kinematic parameters of transpressional deformation obtained by applying a general triclinic transpression analytical model to a tabular-shaped shear zone in the external Betic Chain (Torcal de Antequera massif). According to the results of the analytical model, we used two oblique convergence angles to reproduce the main structural and kinematic features of structural domains observed within the Torcal de Antequera massif (α = 15° for the outer domains and α = 30° for the inner domain). Two parallel inclined backstops (one fixed and the other mobile) reproduce the geometry of the shear zone walls of the natural case. Additionally, we applied digital particle image velocimetry (PIV) method to calculate the velocity field of the incremental deformation. Our results suggest that the spatial distribution of the main structures observed in the Torcal de Antequera massif reflects different modes of strain partitioning and strain localization between two domain types, which are related to the variation in the oblique convergence angle and the presence of steep planar velocity – and rheological – discontinuities (the shear zone walls in the natural case). In the 15° model, strain partitioning is simple and strain localization is high: a single narrow shear zone is developed close and parallel to the fixed backstop, bounded by strike-slip faults and internally deformed by R and P shears. In the 30° model, strain partitioning is strong, generating regularly spaced oblique-to-the backstops thrusts and strike-slip faults. At final stages of the 30° experiment, deformation affects the entire model box. Our results show that the application of analytical modelling to natural transpressive zones related to upper crustal deformation facilitates to constrain the geometrical parameters of analogue models.

Ethiopia is a key site to investigate the interactions between mantle dynamics and surface processes because of the presence of the Main Ethiopian Rift (MER), Cenozoic continental flood basalt volcanism, and plateau uplift. The role of mantle plumes in causing Ethiopia's flood basalts and tectonics has been commonly accepted. However, the location and number of plumes and their impact on surface uplift are still uncertain. Here we present new constraints on the geological and topographic evolution of the Ethiopian Plateau (NW Ethiopia) prior to and after the emplacement of the large flood basalts (40–20 Ma). Using geological information and topographic reconstructions, we show that the large topographic dome that we see today is a long-term feature, already present prior to the emplacement of the flood basalts. We also infer that large-scale doming operated even after the emplacement of the flood basalts. Using a comparison with the present-day topographic setting, we show that an important component of the topography has been and is presently represented by a residual, nonisostatic, dynamic contribution. We conclude that the growth of the Ethiopian Plateau is a long-term, probably still active, dynamically supported process. Our arguments provide constraints on the processes leading to the formation of one of the largest igneous plateaus on Earth.

Highlights • Plate boundary re-organization in the central Mediterranean Sea • Segmentation of the subduction complex along lithospheric transverse faults • STEP faults in the Ionian Sea • Pleistocene active faulting and Mt. Etna formation Abstract The Calabrian Arc is a narrow subduction-rollback system resulting from Africa/Eurasia plate convergence. While crustal shortening is taken up in the accretionary wedge, transtensive deformation accounts for margin segmentation along transverse lithospheric faults. One of these structures is the NNW-SSE transtensive fault system connecting the Alfeo seamount and the Etna volcano (Alfeo-Etna Fault, AEF). A second, NW-SE crustal discontinuity, the Ionian Fault (IF), separates two lobes of the CA subduction complex (Western and Eastern Lobes) and impinges on the Sicilian coasts south of the Messina Straits. Analysis of multichannel seismic reflection profiles shows that: 1) the IF and the AEF are transfer crustal tectonic features bounding a complex deformation zone, which produces the downthrown of the Western lobe along a set of transtensive fault strands; 2) during Pleistocene times, transtensive faulting reactivated structural boundaries inherited from the Mesozoic Tethyan domain which acted as thrust faults during the Messinian and Pliocene; 3) the IF and the AEF, and locally the Malta escarpment, accommodate a recent tectonic event coeval and possibly linked to the Mt. Etna formation. Regional geodynamic models show that, whereas AEF and IF are neighboring fault systems, their individual roles are different. Faulting primarily resulting from the ESE retreat of the Ionian slab is expressed in the northwestern part of the IF. The AEF, on the other hand, is part of the overall dextral shear deformation, resulting from differences in Africa-Eurasia motion between the western and eastern sectors of the Tyrrhenian margin of northern Sicily, and accommodating diverging motions in the adjacent compartments, which results in rifting processes within the Western Lobe of the Calabrian Arc accretionary wedge. As such, it is primarily associated with Africa-Eurasia relative motion.

Seismicity at the northern terminus of the Nazca subduction is diffused over a wide area containing the puzzling seismic feature known as the Bucaramanga nest. We relocate about 5000 earthquakes recorded by the Colombian national seismic network and produce the first 3-D velocity model of the area to define the geometry of the lithosphere subducting below the Colombian Andes. We found lateral velocity heterogeneities and an abrupt offset of the Wadati-Benioff zone at 5°N indicating that the Nazca plate is segmented by an E-W slab tear, that separates a steeper Nazca segment to the south from a flat subduction to the north. The flat Nazca slab extends eastward for about 400 km, before dip increases to ∼50° beneath the Eastern Cordillera, where it yields the Bucaramanga nest. We explain this puzzling locus of intermediate-depth seismicity located beneath the Eastern Cordillera of Colombia as due to a massive dehydration and eclogitization of a thickened oceanic crust. We relate the flat subducting geometry to the entrance at the trench at ca. 10 Ma of a thick - buoyant oceanic crust, likely a volcanic ridge, producing a high coupling with the overriding plate. Sub-horizontal plate subduction is consistent with the abrupt disappearance of volcanism in the Andes of South America at latitudes > 5°N.

9 figures, 1 table., New paleomagnetic results from Neogene sedimentary sequences from the Betic chain (Spain) are here presented. Sedimentary basins located in different areas were selected in order to obtain paleomagnetic data from structural domains that experienced different tectonic evolution during the Neogene. Whereas no rotations have been evidenced in the Late Tortonian sediments in the Guadalquivir foreland basin, clockwise vertical axis rotations have been measured in sedimentary basins located in the central part of the Betics: the Aquitanian to Messinian sediments in the Alcalà la Real basin and the Tortonian and Messinian sediments in the Granada basin. Moreover, counterclockwise vertical axis rotations, associated to left lateral strike-slip faults have been locally measured from sedimetary basins in the eastern Betics: the Middle Miocene to Lower Pliocene sites from the Lorca and Vera basins and, locally, the Tortonian units of the Huercal-Overa basin. Our results show that, conversely from what was believed up to now, paleomagnetic rotations continued in the Betics after Late Miocene, enhancing the role of vertical axis rotations in the recent tectonic evolution of the Gibraltar Arc., Financial support for this work was provided by Projects: a) 01-LECEMA22F [WESTMED] by the European Science Foundation under the EUROCORES Programme EUROMARGINS, through contract No. ERAS-CT-2003-980409 of the European Commission, DG Research, FP6 (M.M., F.C., M.P. and C.F.); b) BTE2003-05057-C02 (A.C–B) and c) REN2001-3868-C03-01 (M.C.) from the MEC and FEDER founding (Spain). We thank G. Booth-Rea for his help in taking samples in the eastern Betics. M.M. and C.F. are particularly grateful to Victor Garcia Duenas who introduce them to the geology of the Betics.

Recent modeling shows that surface processes, such as erosion and deposition, may drive the deformation of the Earth's surface, interfering with deeper crustal and mantle signals. To investigate the coupling between the surface and deep process, we designed a three-dimensional laboratory apparatus, to analyze the role of erosion and sedimentation, triggered by deep mantle instability. The setup is constituted and scaled down to natural gravity field using a thin viscous sheet model, with mantle and lithosphere simulated by Newtonian viscous glucose syrup and silicon putty, respectively. The surface process is simulated assuming a simple erosion law producing the downhill flow of a thin viscous material away from high topography. The deep mantle upwelling is triggered by the rise of a buoyant sphere. The results of these models along with the parametric analysis show how surface processes influence uplift velocity and topography signals.

In the Central Mediterranean, during the last 80 Myr Africa has been slowly converging toward Eurasia, with an average relative velocity of less than 1 cm yr−1. For 50 Myr this slow convergence led to the build-up of the Alpine orogenic belt and to the initiation of consumption of the former Tethys ocean; in the last 30 Myr, however, the whole Central-Western Mediterranean has been swept by rapid episodes of trench migration (up to 6 cm yr−1) and back-arc opening, which consumed the whole former Tethys ocean. We combine plate tectonic history, geological timing, and laboratory modelling to reconstruct the subduction history and the evolution of the Central Mediterranean over the last 80 Myr. We find that the dynamic evolution of the subducting oceanic lithosphere can reconcile the rapid, episodic back-arc migration with the slow African convergence, but only if the subduction process is restricted to the upper mantle.

The effects of convergence rate variations on the evolution of mature orogenic wedges are investigated using temperature-dependent variations in the viscosity of paraffin as an analogue for the natural strength distribution in the crust. Experiments were performed in a thermomechanical shortening apparatus for strain rates ranging between 10−5 and 5×10−5 s−1. Experimental wedges accreted at fast convergence rates show steep surface tapers and narrow deformed areas, whereas at low strain rates they show lower taper topographies and wider deformed areas during self-similar growth. A decrease in convergence rate during deformation corresponds to a change from localised to delocalised deformation in concomitance with the outward migration of the deformation front. The model results provide insights into deformation partitioning, geometry and accretionary mechanisms within orogenic wedges.

The Abruzzi–Molise sector in the Central Apennines is a part of a fold and thrust belt that has been deforming since the Late Cretaceous as a result of collision tectonics between the European and Adriatic plates. The superposition of different deformational styles highly reworked the originally complex palaeogeography of this portion of the southern Tethyan margin. Analogue modelling has been performed on thrusting mechanisms in the Abruzzi–Molise area in order to (1) reduce the number of admissible hypotheses regarding palaeogeographic setting, and (2) define thrusting mechanics. Both of these goals are crucial for hydrocarbon exploration purposes. The sandbox apparatus used to simulate the undeformed passive margin consisted of a thin compartment juxtaposed against a thick one along a linear boundary having variable geometries and mechanical stratigraphy; a rigid but mobile backstop was used to deform the stratigraphy in a Coulomb thrust wedge. Results from six experiments show that the geometric relationships between different structural units depend on the distribution of palaeogeographic domains. These domains are defined by mechanical and/or geometrical parameters, such as the orientation between the maximum compression direction and the palaeogeographic boundary, the mechanical stratigraphy and the thickness of the successions reproduced in the models. The present-day tectonic styles and Meso–Cenozoic palaeogeography of the Abruzzi–Molise area are discussed in terms of the mechanisms and structures analysed through the models.

Obduction is an odd geodynamic process characterized by the emplacement of dense oceanic “ophiolites” atop light continental plates in convergent settings. We herein present analog models specifically designed to explore the conditions (i.e., sharp increase of plate velocities — herein coined as ‘acceleration’, slab interaction with the 660 km discontinuity, ridge subduction) under which obduction may develop as a result of subduction initiation. The experimental setup comprises an upper mantle modeled as a low-viscosity transparent Newtonian glucose syrup filling a rigid Plexiglas tank and high-viscosity silicone plates. Convergence is simulated by pushing a piston with plate tectonics like velocities (1–10 cm/yr) onto a model comprising a continental margin, a weakness zone with variable resistance and dip ( W ), an oceanic plate (with or without a spreading ridge), a preexisting subduction zone ( S ) dipping away from the piston and an upper active continental margin, below which the oceanic plate is being subducted at the start of the model (as for the Neotethyan natural example). Several configurations were tested over thirty-five parametric models, with special emphasis on comparing different types of weakness zone and the degree of mechanical coupling across them. Measurements of displacements and internal deformation allow for a precise and reproducible tracking of deformation. Models consistently demonstrate that once conditions to initiate subduction are reached, obduction may develop further depending on the effective strength of W .Results(1)constrainthe range of physical conditions required for obduction to develop/nucleate and (2) underline the key role of such perturbations for triggering obduction, particularly plate ‘acceleration’. They provide an explanation to the short-lived Peri-Arabic obduction, which took place along thousands of km almost synchronously (within ∼ 50–10 Myr), from Turkey to Oman, while the subduction zone beneath Eurasia became temporarily jammed. They also demonstrate that the emplacement of dense, oceanic material on continental lithosphere is not a mysterious process requiring extraordinary boundary conditions but results from large-scale, normal (oceanic then continental) subduction processes.

We present analogue models simulating the subduction of a buoyant ridge oriented perpendicularly or obliquely with respect to the trench, beneath an advancing overriding plate. The convergence velocity is imposed by lateral boundary conditions in this experimental set. We analyze the three-dimensional geometry of the slab, the deformation and topography of the overriding plate. Experiments suggest that ridge subduction diminishes the dip of the slab, eventually leading to the appearance of a horizontal slab segment in case boundary conditions impose a rapid convergence. This result contrasts with that obtained in free subduction experiments, in which ridge subduction diminishes the convergence velocity which, in turn, increases the dip of the slab beneath the ridge. The slab dip decrease is accompanied by the indentation of the overriding plate by the ridge, resulting in arc curvature. Experiments suggest that indentation is larger for small convergence velocity and large slab dip. Ridge subduction also uplifts the overriding plate. Uplift first occurs close to the trench (~ fore-arc area) and is accompanied by the flexural subsidence of the overriding plate behind the uplifted area (~ back-arc subsidence). The uplifted area migrates within the overriding plate interiors following the appearance of a horizontal slab segment. These results are compared with natural examples of ridge subduction in the circum-Pacific area. They explain why ridge subduction may have contrasted effects on the overriding plate dynamics depending on the global conditions that constrain the converging system.

Crustal earthquake focal mechanisms are investigated in the Calabrian Arc region, where the western Mediterranean subduction process is close to ending and the residual Ionian subducting slab affected by gravitational roll-back produces a strong variation of faulting regimes at shallow depth along the local section of the convergent margin. An updated database of earthquake focal mechanisms has been compiled by selecting the best-quality solutions available in the literature and in catalogs, and by adding 17 new solutions estimated in the present work. A total of 164 mechanisms are included in this database, 142 computed by waveform inversion and 22 by analysis of P-wave first motions from earthquakes with good network coverage and no less than 14 records. 60% of the solutions included in the database have never been used for regional-scale geodynamic investigations before the present study, and this makes the compiled database substantially new for our application. Focal mechanisms have been inverted for stress tensor orientations to obtain the principal stress axes over the study region. Results are compatible with three major tectonic domains subject to markedly different regional stresses along the arc. These three domains are separated by two transitional domains, which are located on top of the Ionian subducting slab edges and are likely forced in their horizontal transfer kinematics by the different tectonic regimes occurring in the adjacent major domains rather than by the regional tectonics. This along-arc differential tectonics is at least in part interpreted as the surface expression of the different deep mechanisms occurring in correspondence of the narrow Ionian slab and their lateral edges. Open tectonic questions are emphasized and proposed for future studies.

The Neogene tectonic evolution of the western Mediterranean region is accompanied and outlasted by diffuse magmatism. This study describes the tectonic setting, the petrography and geochemistry of the Early Miocene granitic dyke swarm that occurs in the Oued Amter area, in the core of the Alboran Domain of the Moroccan Rif. The structural setting indicates dyke intrusion assisted and controlled by strike-slip tectonics that operated through conjugate, NW–SE left-lateral and NE–SW right-lateral fault strands. The overall composition of the dykes (the high SiO2 contents (69–77 wt.%), coupled with low concentrations of TiO2, MgO, FeO, CaO), the high 87Sr/86Sr (0.719–0.722) and the low 143Nd/144Nd (ca. 0.5120) values point to a predominantly crustal origin of these magmatic bodies, compatible with a process involving muscovite–(biotite) dehydration melting of fertile metasedimentary sources. The isotopic signature is similar to that of the basement rocks of the Alboran Sea, suggesting for a similar crustal source for the Early Miocene felsic magmatism of the Betic–Rif realm. When framed within the regional setting, these data are used to propose a synthetic geodynamic model for the Early Miocene tectonic and magmatic evolution of the western Mediterranean region.

Subduction megathrust earthquakes occur at the interface between the subducting and overriding plates. These hazardous phenomena are only partially understood because of the absence of direct observations, the restriction of the instrumental seismic record to the past century, and the limited resolution/completeness of historical to geological archives. To overcome these restrictions, modeling has become a key-tool to study megathrust earthquakes. We present a novel model to investigate the seismic cycle at subduction thrusts using complementary analog (paper 1) and numerical (paper 2) approaches. Here we introduce a simple scaled gelatin-on-sandpaper setup including realistic tectonic loading, spontaneous rupture nucleation, and viscoelastic response of the lithosphere. Particle image velocimetry allows to derive model deformation and earthquake source parameters. Analog earthquakes are characterized by "quasi-periodic" recurrence. Consistent with elastic theory, the interseismic stage shows rearward motion, subsidence in the outer wedge and uplift of the "coastal area" as a response of locked plate interface at shallow depth. The coseismic stage exhibits order of magnitude higher velocities and reversal of the interseismic deformation pattern in the seaward direction, subsidence of the coastal area, and uplift in the outer wedge. Like natural earthquakes, analog earthquakes generally nucleate in the deeper portion of the rupture area and preferentially propagate upward in a crack-like fashion. Scaled rupture width-slip proportionality and seismic moment-duration scaling verifies dynamic similarities with earthquakes. Experimental repeatability is statistically verified. Comparing analog results with natural observations, we conclude that this technique is suitable for investigating the parameter space influencing the subduction interplate seismic cycle. © 2012. American Geophysical Union. All Rights Reserved.

We designed three-dimensional dynamically self-consistent laboratory models of subduction to analyse the relationships between overriding plate deformation and subduction dynamics in the upper mantle. We investigated the effects of the subduction of a lithosphere of laterally variable buoyancy on the temporal evolution of trench kinematics and shape, horizontal flow at the top of the asthenosphere, dynamic topography and deformation of the overriding plate. Two subducting units, which correspond to a negatively buoyant oceanic plate and positively buoyant continental one, are juxtaposed via a trench-perpendicular interface (analogue to a tear fault) that is either fully-coupled or shear-stress free. Differential rates of trench retreat, in excess of 6 cm yr−1 between the two units, trigger a more vigorous mantle flow above the oceanic slab unit than above the continental slab unit. The resulting asymmetrical sublithospheric flow shears the overriding plate in front of the tear fault, and deformation gradually switches from extension to transtension through time. The consistency between our models results and geological observations suggests that the Late Cenozoic deformation of the Aegean domain, including the formation of the North Aegean Trough and Central Hellenic Shear zone, results from the spatial variations in the buoyancy of the subducting lithosphere. In particular, the lateral changes of the subduction regime caused by the Early Pliocene subduction of the old oceanic Ionian plate redesigned mantle flow and excited an increasingly vigorous dextral shear underneath the overriding plate. The models suggest that it is the inception of the Kefalonia Fault that caused the transition between an extension dominated tectonic regime to transtension, in the North Aegean, Mainland Greece and Peloponnese. The subduction of the tear fault may also have helped the propagation of the North Anatolian Fault into the Aegean domain.

""Convergent margins often exhibit spatial and temporal correlations between trench curvature, overriding plate shortening and topography uplift that provide insights into the dynamics of subduction. The Andean system, where the Nazca plate plunges beneath continental South America, is commonly regarded as the archetype of this class of tectonics systems. There is distinctive evidence that the degree of mechanical coupling between converging plates, i.e. the amount of resistive force mutually transmitted in the direction opposite to their motions, may be at the present-day significantly higher along the central Andean margin compared to the northern and southern limbs. However quantitative estimates of such resistance are still missing and would be desirable. Here we present laboratory models of subduction performed to investigate quantitatively how strong lateral coupling variations need to be to result in trench curvature, tectonic shortening and distribution of topography comparable to estimates from the Andean margin. The analogue of a two-layers Newtonian lithosphere\\\/upper mantle system is established in a silicone putty\\\/glucose syrup tank-model where lateral coupling variations along the interface between subducting and overriding plates are pre-imposed. Despite the simplicity of our setup, we estimate that coupling in the central margin as large as 20% of the driving force is sufficient to significantly inhibit the ability of the experimental overriding plate to slide above the subducting one. As a consequence, the central margin deforms and shortens more than elsewhere while the trench remains stationary, as opposed to the advancing lateral limbs. This causes the margin to evolve into a peculiar shape similar to the present-day trench of the Andean system. (C) 2011 Elsevier B.V. All rights reserved.""

International audience; The geometric and kinematic evolution of the Andes provides insight onto the nature of the force balance beneath the South American plate. While the Andean load is opposed on its western edge by the force induced by subduction of the Nazca plate, its more elusive eastern counterpart, which we explore herein, requires some contribution from the mantle beneath the South Atlantic. Using a mantle flow model, we show that the Andes owe their existence to basal drag beneath South America caused by a cylindrical convection cell under the South Atlantic. We find that the observed Andean uplift requires both westward push fromactive upwelling beneath Africa andwestward drag toward the downgoing Nazca slab. These mutually-reinforcing downwellings and upwellings amount to 38% and 23% of the total driving force, respectively. Further decomposition reveals that the South Atlantic cell is most vigorous near its center, rendering the net drag force higher where the Andes also reach their highest elevation. Kinematic reconstructions suggest that the South Atlantic cell could have grown owing to the migration of the Nazca slab until ~50 Ma. We propose that from 50 Ma onwards, the cell may have ceased growing westward because (i) it had reached an optimal aspect ratio and (ii) the Nazca slab became anchored into the lower mantle. Continued westward motion of the plates, however, moved the surface expressions of spreading and convergence away from the upwelling and downwelling arms of this cell. Evidence for this scenario comes from the coeval tectonic, morphologic, and magmatic events in Africa and South America during the Tertiary.

""The relief around subduction zones results from the interplay of dynamic processes that may locally exceed the (iso)static contributions. The viscous dissipation of the energy in and around subduction zones is capable of generating kilometer scale vertical ground movements. In order to evaluate dynamic topography in a self-consistent subduction system, we carried out a set of laboratory experiments, wherein the lithosphere and mantle are simulated by means of Newtonian viscous materials, namely silicone putty and glucose syrup. Models are kept in their most simple form and are made of negative buoyancy plates, of variable width and thickness, freely plunging into the syrup. The surface of the model and the top of the slab are scanned in three dimensions. A forebulge systematically emerges from the bending of the viscous plate, adjacent to the trench. With a large wavelength, dynamic pressure offsets the foreside and backside of the slab by similar to 500 m on average. The suction, that accompanies the vertical descent of the slab depresses the surface on both sides. At a distance equal to the half-width of the slab, the topographic depression amounts to similar to 500 m on average and becomes negligible at a distance that equals the width of the slab. In order to explore the impact of slab rollback on the topography, the trailing edge of the plates is alternatively fixed to (fixed mode) and freed from (free mode) the end wall of the tank. Both the pressure and suction components of the topography are similar to 30% lower in the free mode, indicating that slab rollback fosters the dynamic subsidence of upper plates. Our models are compatible with first order observations of the topography around the East Scotia, Tonga, Kermadec and Banda subduction zones, which exhibit anomalous depths of nearly 1 km as compared to adjacent sea floor of comparable age. (C) 2011 Elsevier B.V. All rights reserved. RI Guillaume, Benjamin\\\/D-6570-2012""

[1] The fate of a convergent continental margin is investigated. We perform a set of 2D numerical models to study how and why continental collision can evolve in different scenarios. Since the rheology of continental lithosphere has a major control on the dynamics of subduction, we explore a range of different lithosphere and lower crust viscosity values to understand their sensitivity on the possible scenarios. We find that with a rheologically layered crust both delamination and break-off are feasible. We identify three modes: (1) slab detachment, in which the lithospheric mantle and the crust are strongly coupled, subduction slows down and the slab eventually breaks; (2) delamination of the lithospheric mantle that separates from the crust and continue to subduct and (3) an intermediate mode where the lithospheric mantle and the crust remain partially coupled, resulting in an initial stage of delamination followed by the slow down and cessation of subduction.

Alpine orogens in the central Mediterranean region have revealed the concomitance of crustal extension in back-arc domain and crustal shortening in frontal domain. Quantitative data of deformation in the frontal orogenic wedges are necessary to understand how the shortening-extension pair evolves in terms of structures, orogenic transport, timing, and exhumation rate. This paper deals with kinematics and ages of the frontal thrust systems of the Calabria-Peloritani Arc (Italy) exposed in the eastern Sila Massif. We first present structural fieldwork, onshore and offshore well log data, and new apatite fission-track (AFT) thermochronology. Then, we describe the structural architecture of the studied area as an ENE-verging stacking of thrust sheets involving basement units and syn-orogenic sediments. The AFT study documents that thrust sheets entered the partial annealing zone from 18. Ma to 13. Ma. This Early-Middle Miocene thrusting phase was coeval with exhumation of high-pressure/low temperature metamorphic rocks in the hinterland of the orogen (Coastal Chain area), mainly driven by top-to-the-W extensional tectonics. Opposite kinematic shear senses (contractional top-to-the-E and extensional top-to-the-W) and different exhumation rates (slow in the frontal, more rapid in the hinterland) are framed in a tectonic scenario of a critically tapered orogenic wedge during the eastward retreating of the Apennine slab. © 2012 Elsevier B.V.

[1] Laboratory models are realized to investigate the role of interface roughness, driving rate, and pressure on friction dynamics. The setup consists of a gelatin block driven at constant velocity over sand paper. The interface roughness is quantified in terms of amplitude and wavelength of protrusions, jointly expressed by a reference roughness parameter obtained by their product. Frictional behavior shows a systematic dependence on system parameters. Both stick slip and stable sliding occur, depending on driving rate and interface roughness. Stress drop and frequency of slip episodes vary directly and inversely, respectively, with the reference roughness parameter, reflecting the fundamental role for the amplitude of protrusions. An increase in pressure tends to favor stick slip. Static friction is a steeply decreasing function of the reference roughness parameter. The velocity strengthening/weakening parameter in the state- and rate-dependent dynamic friction law becomes negative for specific values of the reference roughness parameter which are intermediate with respect to the explored range. Despite the simplifications of the adopted setup, which does not address the problem of off-fault fracturing, a comparison of the experimental results with the depth distribution of seismic energy release along subduction thrust faults leads to the hypothesis that their behavior is primarily controlled by the depth- and time-dependent distribution of protrusions. A rough subduction fault at shallow depths, unable to produce significant seismicity because of low lithostatic pressure, evolves into a moderately rough, velocity-weakening fault at intermediate depths. The magnitude of events in this range is calibrated by the interplay between surface roughness and subduction rate. At larger depths, the roughness further decreases and stable sliding becomes gradually more predominant. Thus, although interplate seismicity is ultimately controlled by tectonic parameters (velocity of the plates/trench and the thermal regime), the direct control is exercised by the resulting frictional properties of the plate interface.

Mobile belts are long-lived deformation zones composed of an ensemble of crustal fragments, distributed over hundreds of kilometres inside continental convergent margins. The Mediterranean represents a remarkable example of this tectonic setting: the region hosts a diffuse boundary between the Nubia and Eurasia plates comprised of a mosaic of microplates that move and deform independently from the overall plate convergence. Surface expressions of Mediterranean tectonics include deep, subsiding backarc basins, intraplate plateaux and uplifting orogenic belts. Although the kinematics of the area are now fairly well defined, the dynamical origins of many of these active features are controversial and usually attributed to crustal and lithospheric interactions. However, the effects of mantle convection, well established for continental interiors, should be particularly relevant in a mobile belt, and modelling may constrain important parameters such as slab coherence and lithospheric strength. Here we compute global mantle flow on the basis of recent, high-resolution seismic tomography to investigate the role of buoyancy-driven and plate-motion-induced mantle circulation for the Mediterranean. We show that mantle flow provides an explanation for much of the observed dynamic topography and microplate motion in the region. More generally, vigorous small-scale convection in the uppermost mantle may also underpin other complex mobile belts such as the North American Cordillera or the Himalayan-Tibetan collision zone.

We present dynamically self-consistent mantle-scale laboratory models that have been conducted to improve our understanding of the influence of slab window opening on subduction dynamics, mantle flow and associated dynamic topography over geological time scales. The adopted setup consists of a two-layer linearly viscous system simulating the subduction of a fixed plate of silicone (lithosphere) under negative buoyancy in a viscous layer of glucose syrup (mantle). Our experimental setting is also characterized by a constant-width rectangular window located at the center of a laterally confined slab, modeling the case of the interaction of a trench-parallel spreading ridge with a wide subduction zone. We found that the opening of a slab window does not produce consistent changes of the geometry and the kinematics of the slab. On the contrary, slab-induced mantle circulation, quantified both in the vertical and horizontal sections using the Feature Tracking image analysis technique, is strongly modified. In particular, rollback subduction and the opening of the slab window generate a complex mantle circulation pattern characterized by the presence of poloidal and toroidal components, with the importance of each evolving according to kinematic stages. Mantle coming from the oceanic domain floods through the slab window, indenting the supra-slab mantle zone and producing its deformation without any mixing between mantle portions. The opening of the slab window and the upwelling of sub-slab mantle produce a regional-scale non-isostatic topographic uplift of the overriding plate that would correspond to values ranging between ca. 1 and 5 km in nature. Assuming that our modeling results can be representative of the natural behavior of subduction zones, we compared them to the tectonics and volcanism of the Patagonian subduction zone. We found that the anomalous backarc volcanism that has been developing since the middle Miocene could result from the lateral flow of sub-slab mantle and that part of the Patagonian uplift could be dynamically supported.

No abstract has been provided

International audience; The origin of forces driving the deformation of the continental crust near subduction zones and especially in backarc regions is debated. This work is based on a compilation of SKS fast splitting directions that give an image of flow lines in the mantle around the Mediterranean subduction zones and a comparison with stretching and shear directions in metamorphic core complexes that show the pattern of deformation at the scale of the middle and lower crusts. We find that : (1) the two sets of directions are parallel in the three main backarc regions, namely the Alboran Sea, the Tyrrhenian Sea and the Aegean Sea showing that the lithosphere deforms with the same direction of stretching in the crust and the mantle, suggesting that (2) crustal deformation is mainly driven from below by slab retreat, and (3) the lithospheric fabric is reset within a few millions of years in backarc environments.