Search

In convergent systems, the interplay between tectonics, erosion, and sedimentation controls the orogenic evolution. The nature of the interactions between these factors is still elusive due to the complex feedbacks that operate across different temporal and spatial scales. Here, we investigate these feedbacks with analog models of landscape evolution designed to account for both tectonic forcing and surface processes, using a water-saturated granular material that allows to simulate contemporary brittle deformation and surface processes. The deformation is imposed by the movement of a rigid backstop, and surface processes are triggered by simulated rainfall and runoff. We vary the convergence velocity, rainfall rate, and basal angle of the box, testing how different boundary conditions affect the balance between tectonics and surface processes. We measure the competition between input fluxes (tectonics) and output fluxes (erosion) of material, showing how sedimentation strongly affects the balance between these fluxes. The results suggest that the experimental equilibrium between tectonics and erosion can be achieved, in the analog models, only for low convergence rates (about 10 mm hr−1) and/or for high basal angle (>2°, limited sedimentation). If the foreland is overfilled with sediments and/or if convergence velocity is higher, channels decrease their erosional efficiency, moving the dynamic equilibrium between tectonics and erosion toward the former.

In the last decades, seismotectonic analog models have been developed to better understand many aspects of the seismic cycle. Differently from other lab-quake experiments, seismotectonic models mimic the first order characteristics of the seismic cycle in a scaled fashion. Here we introduce Foamquake: a novel seismotectonic model with a granular frictional interface that as a whole behaves elastoplastically. The model experiences cycles of elastic loading and release via spontaneous nucleation of frictional instabilities at the base of an elastic foam wedge, hereafter called foamquakes. These analog earthquakes show source parameters (i.e., moment-duration and moment-rupture area) scaling as great interplate earthquakes and a coseismic displacement of few tens of meters when scaled to nature. Models with two asperities separated by a barrier can be performed with Foamquake given the 3D nature of the setup. Such model configuration generates sequences of full and partial ruptures with different recurrence intervals as well as rupture cascades. By tuning the normal load acting on individual asperities, Foamquake reproduces superimposed cycles rupture patterns such as those observed along natural megathrusts. The physical properties of asperities and barriers affect model seismic behavior. Asperities with similar properties and low yield strength fail preferentially in a simultaneous manner. The combination of all those characteristics suggests that Foamquake is a valuable tool for investigating megathrust seismicity and seismic processes that depend on the 3D nature of the subduction environment.

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.

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.

Geophysical Journal International, 157 (3), ISSN:0956-540X, ISSN:1365-246X

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.

Giant earthquakes (Mw≥8.5) usually occur on the boundary between subducting and overriding plates of converging margins, but it is not yet clear which (if any) subduction zones are more prone to produce such a kind of events. Here we show that subduction zones may have different capabilities to produce giant earthquakes. We analyze the frequency-magnitude distribution of the interplate earthquakes at subduction zones that occurred during 1976–2007 and calculate the propensity (defined as the average annual rate) of giant events for about half of the subduction zones. We find that the b value of interplate earthquakes is significantly different among the subduction zones, and out-of-sample giant earthquakes (before 1976 and after 2007) have occurred preferentially in high-propensity areas. Besides the importance for seismic hazard assessment and risk mitigation, our results seem to indicate that a higher seismicity rate does not necessarily imply a higher likelihood to generate giant earthquakes, and the way in which the stress is released at a subduction interface does not change significantly after the occurrence of such events. ©2016. American Geophysical Union. All Rights Reserved.

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.

The evaluation of the kinematic setting and the structural control of volcanic arcs are important in defining the tectono-magmatic processes along convergent plate boundaries. However, our knowledge is fragmented and the available data highlight different behaviours. This study analyzes the kinematic setting, the structural control and the volcanic productivity of 16 arcs. These arcs are characterized by predominant extensional, compressional, strike–slip or oblique motions. There is an overall coupling between the normal vs. parallel motions along the arc and those of the underlying slab. Therefore, the higher the trench-normal (or parallel) component of the subduction rate, the higher is the amount of arc-normal (or parallel) motion. This relation confirms that strain partitioning is, in general, feasible at many convergent settings, involving also the volcanic arc portion. The arc-normal motion may be characterized by extension or compression, as a function of the increase in the trench-normal convergence (or subduction) velocity. The lack of an evident relation between the subduction rate (or convergence rate) and the volcanic productivity of the arc is partly inconsistent with previous studies, which highlight a relation between the subduction rate and melt production below arcs. This discrepancy mainly suggests that different processes control the generation, rise and eruption of magma at different depths, varying the intrusive/extrusive ratio along arcs. The structure of the arc does not control the distribution of the volcanoes; however, it does control the volcanic output through different processes, even though regional or local extension (associated with strike–slip or compressive structures) is the ultimate requisite, in any setting. In general, the higher the amount of extension, the higher is the output rate along the arc.

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.

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.

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.

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

International audience; A variety of atypical plume-like structures and focused upwellings that are not rooted in the lower mantle have recently been discussed, and seismological imaging has shown ubiquitous small-scale convection in the uppermost mantle in regions such as the Mediterranean region, the western US, and around the western Pacific. We argue that the three-dimensional return flow and slab fragmentation associated with complex oceanic subduction trajectories within the upper mantle can generate focused upwellings and that these may play a significant role in regional tectonics. The testable surface expressions of this process are the outside-arc alkaline volcanism, topographic swell, and low-velocity seismic anomalies associated with partial melt. Using three-dimensional, simplified numerical subduction models, we show that focused upwellings can be generated both ahead of the slab in the back-arc region (though -five times further inward from the trench than arc-volcanism) and around the lateral edges of the slab (in the order of 100 km away from slab edges). Vertical mass transport, and by inference the associated decompression melting, in these regions appears strongly correlated with the interplay between relative trench motion and subduction velocities. The upward flux of material from the depths is expected to be most pronounced during the first phase of slab descent into the upper mantle or during slab fragmentation. We discuss representative case histories from the Pacific and the Mediterranean where we find possible evidence for such slab-related volcanism.

International audience; Recent global kinematic studies reveal that most of the trenches roll back but a significant number of them advance toward the upper plate. Those advancing trenches are mostly located in the Western Pacific and correspond to the subduction of very old, Mesozoic oceanic lithosphere. While retreating trenches are commonly explained by the slab pull action of the descending lithosphere, the origin of advancing trenches is still debated. Since this relationship is dependent upon the adopted reference frame, we select region where geological studies show the variability of trench migration style with time. The Izu-Bonin-Mariana (IBM) region represents a key example. The detailed reconstruction of the trench migration of the IBM subducting system reveals that after a long episode of asymmetric rollback, the IBM trench recently started advancing. We propose that this change from retreating to advancing trench mode results from the subduction of progressively older and stiffer lithospheric material. We test this hypothesis by means of two-dimensional (2-D) numerical models, reproducing the effects of the lithospheric aging during subduction. The result of our numerical tests shows that the entrance of old and stiff lithosphere forces the trench to advance because the increasing stiffness of the slab prevents the slab to unbend once it has subducted. We adapt this physical result to the IBM evolution, showing that the switch from trench retreat to trench advance can reconcile the shape of IBM slab as previously suggested. Finally we discuss the possibility that the subduction of the Cretaceous lithosphere might have triggered the change of the motion of the Pacific plate around 5 Ma.

International audience; Laboratory models have been conducted to improve our understanding of the role that the resistance of the slab to bending and its coupling to the ambient mantle play in subduction dynamics over geological time scales. Our models are set up with a viscous plate of silicone (lithosphere) subducting under negative buoyancy in a viscous layer of glucose syrup (mantle). For our study, the lithosphere/upper mantle viscosity contrast has been systematically varied, from similar to 10 to similar to 10(5) in order to explore the parameter space between weak and strong slab dynamics. We found that subduction is characterized by a retreating mode for viscosity ratios > 10(4), by the coexistence of a retreating mode and an advancing mode for viscosity ratios between similar to 10(4) and similar to 10(2), and quasi-stationary, Rayleigh-Taylor like behaviour for ratios < 10(2). By combining our experimental results and kinematic data from current subduction zones in four reference frames which differ in the amount of net rotation, we infer that a lithosphere/upper mantle viscosity contrast of 150-500 is necessary to obtain realistic trench/subducting plate velocity ratios as well as the variability of subduction styles observed in nature.

A combination of statistical studies on present-day subduction zones and three-dimensional (3D) laboratory models is performed with the aim to clarify the way that plate kinematics control the geometry of the slab and the overriding plate deformation in subduction zones. In 3D laboratory models, the analogue of a two layer linearly viscous lithosphere–upper mantle system is achieved by means of silicon putty glucose syrup tank experiment. The subducting and overriding plate velocities are systematically changed by exploring the variability field of natural plate kinematics. Both statistical and modelling approaches recognize the importance of overriding plate motion on subduction process behavior: (1) trenches migrate at a rate close to the overriding plate motion, but always move slower than the overriding plates. The mechanism at work is a direct consequence of “slab anchoring” opposed by both lithosphere and mantle viscous resistance and is responsible for overriding plate deformation and slab geometry variability. (2) An overriding plate shortens when the overriding plate moves toward the trench and conditions that are favourable for overriding plate extension are created when the overriding plate moves away from the trench. (3) Shallow and steep dips are found if the overriding plate moves toward and away from the trench, respectively.

The motion of oceanic plates is commonly related to the subduction of cold and dense oceanic material into the mantle. These models predict plate velocities from subduction velocities but the trench motion is not directly included in the computation. Here, using a recent compilation of a global data set, we found that the motion of trenches (either advancing or retreating with respect to upper plates) scales with their corresponding subducting plates motion. Based on simple experimental tests, we found that subduction of strong slabs inside the upper mantle correctly predicts these kinematic relationships. We deduce that the motion of the trenches represent the surface manifestation of the resistance encountered by the subducting lithosphere to bend and penetrate within the upper mantle.

International audience; We model using analogue experiments the subduction of buoyant ridges and plateaus to study their effect on slab dynamics. Experiments show that simple local (1-D) isostatic considerations are not appropriate to predict slab behaviour during the subduction of a buoyant ridge perpendicular to the trench, because the rigidity of the plate forces the ridge to subduct with the dense oceanic lithosphere. Oceanic ridges parallel to the trench have a stronger effect on the process of subduction because they simultaneously affect a longer trench segment. Large buoyant slab segments sink more slowly into the asthenosphere, and their subduction result in a diminution of the velocity of subduction of the plate. We observe a steeping of the slab below those buoyant anomalies, resulting in smaller radius of curvature of the slab that augments the energy dissipated in folding the plate and further diminishes the velocity of subduction. When the 3-D geometry of a buoyant plateau is modelled, the dip of the slab above the plateau decreases, as a result of the larger velocity of subduction of the dense `normal' oceanic plate on both sides of the plateau. Such a perturbation of the dip of the slab maintains long time after the plateau has been entirely incorporated into the subduction zone. We compare experiments with the present-day subduction zone below South America. Experiments suggest that a modest ridge perpendicular to the trench such as the present-day Juan Fernandez ridge is not buoyant enough to modify the slab geometry. Already subducted buoyant anomalies within the oceanic plate, in contrast, may be responsible for some aspects of the present-day geometry of the Nazca slab at depth.

International audience; Three-dimensional laboratory experiments have been designed to investigate the way slab-bearing plates move during subduction inside the mantle. In our experiments a viscous plate of silicone (lithosphere) subducts under its negative buoyancy in a viscous layer of pure honey (mantle). Varying thickness, width, viscosity, and density of the plate and mantle, three characteristic modes of subduction are observed: a retreating trench mode (mode I), a retreating trench mode following a transient period of advancing trench (mode II), and an advancing trench mode (mode III). These modes are characterized by different partitioning of the amount of subduction into plate and trench motion. Our experiments show that the velocity of subduction can be modeled by the dynamic interaction between acting and resisting forces, where lithospheric bending represents 75-95% of the total resisting forces. However, our experimental results also show the impossibility to predict a priori the plate velocity only from the velocity of subduction without considering trench migrations. We find experimentally that the lithospheric radius of curvature, which depends upon plate characteristics (stiffness and thickness) and the mantle thickness, exerts a primary control on the trench behavior. Our results suggest that the complexity of the style of subduction could be controlled by geometrical rules of a plate bending inside a stratified mantle. The Earth system is in the crucial range for the interplay between the rigidity of the plate and the mantle stratification: this setting may be the responsible for the complexity of the past and present tectonic styles.

We have performed 3-D scaled lithospheric experiments to investigate the role of the gravitational force exerted by a subducting slab on the deformation of the subducting plate itself. Experiments have been constructed using a dense silicone putty plate, to simulate a thin viscous lithosphere, floating in the middle of a large box filled with glucose syrup, simulating the upper mantle. We examine three different plate configurations: (i) subduction of a uniform oceanic plate, (ii) subduction of an oceanic–continental plate system and, (iii) subduction of a more complex oceanic–continental system simulating the asymmetric Africa–Eurasia system. Each model has been performed with and without the presence of a circular weak zone inside the subducting plate to test the near-surface weakening effect of a plume activity. Our results show that a subducting plate can deform in its interior only if the force distribution varies laterally along the subduction zone, i.e. by the asymmetrical entrance of continental material along the trench. In particular, extensional deformation of the plate occurs when a portion of the subduction zone is locked by the collisional process. The results of this study can be used to analyze the formation of the Arabian plate. We found that intraplate stresses, similar to those that generated the Africa–Arabia break-up, can be related to the Neogene evolution of the northern convergent margin of the African plate, where a lateral change from collision (Mediterranean and Bitlis) to active subduction (Makran) has been described. Second, intraplate stress and strain localization are favored by the presence of a weakness zone, such as the one generated by the Afar plume, producing a pattern of extensional deformation belts resembling the Red Sea–Gulf of Aden rift system.

[1] The physics governing the seismic cycle at seismically active subduction zones remains poorly understood due to restricted direct observations in time and space. To investigate subduction zone dynamics and associated interplate seismicity, we validate a continuum, visco-elasto-plastic numerical model with a new laboratory approach (Paper 1). The analogous laboratory setup includes a visco-elastic gelatin wedge underthrusted by a rigid plate with defined velocity-weakening and -strengthening regions. Our geodynamic simulation approach includes velocity-weakening friction to spontaneously generate a series of fast frictional instabilities that correspond to analog earthquakes. A match between numerical and laboratory source parameters is obtained when velocity-strengthening is applied in the aseismic regions to stabilize the rupture. Spontaneous evolution of absolute stresses leads to nucleation by coalescence of neighboring patches, mainly occurring at evolving asperities near the seismogenic zone limits. Consequently, a crack-, or occasionally even pulse-like, rupture propagates toward the opposite side of the seismogenic zone by increasing stresses ahead of its rupture front, until it arrests on a barrier. The resulting surface displacements qualitatively agree with geodetic observations and show landward and, from near the downdip limit, upward interseismic motions. These are rebound and reversed coseismically. This slip increases adjacent stresses, which are relaxed postseismically by afterslip and thereby produce persistent seaward motions. The wide range of observed physical phenomena, including back-propagation and repeated slip, and the agreement with laboratory results demonstrate that visco-elasto-plastic geodynamic models with rate-dependent friction form a new tool that can greatly contribute to our understanding of the seismic cycle at subduction zones.

International audience; The role of seafloor roughness on the seismogenic behavior of subduction zones has been increasingly addressed over the past years, although their exact relationship remains unclear. Do subducting features like seamounts, fracture zones, or submarine ridges act as barriers, preventing ruptures from propagating, or do they initiate megathrust earthquakes instead? We address this question using a global approach, taking into account all oceanic subduction zones and a 117-year time window of megathrust earthquake recording. We first compile a global database, SubQuake, that provides the location of a rupture epicenter, the overall rupture area, and the region where the largest displacement occurs (the seismic asperity) for M W ≥ 7.5 subduction interplate earthquakes. With these data, we made a quantitative comparison with the seafloor roughness seaward of the trench, which is assumed to be a reasonable proxy for the subduction interface roughness. We compare the spatial occurrence of megathrust ruptures, seismic asperities, and epicenters, with two roughness parameters: the short-wavelength roughness R SW (12-20 km) and the long-wavelength roughness R LW (80-100 km). We observe that ruptures with M W ≥ 7.5 tend to occur preferentially on smooth subducting seafloor at long wavelengths, which is especially clear for the M W > 8.5 events. At both short and long wavelengths, seismic asperities show a more amplified relation with smooth seafloor than rupture segments in general. For the epicenter correlation, we see a slight difference in roughness signal, which suggests that there might be a physical relationship between rupture nucleation and subduction interface roughness. Plain Language Summary Subduction zones are regions on Earth where an oceanic plate dives below another plate. Earthquakes that occur along the contact between plates in such regions are among the largest and most destructive on Earth. To better understand where these large earthquakes are most likely to occur, we look at the effect of seafloor roughness. A rough seafloor is often characterized by many topographic features, such as seamounts or ridges, while a smooth seafloor is generally more flat. On a global scale, we compared the roughness of the incoming seafloor of the downgoing plate, with the occurrence of large earthquakes in each subduction zone. We find that the seafloor in front of large earthquakes is generally smoother than in areas where no large earthquakes have occurred. This is the clearest for very large earthquakes, with magnitudes larger than 8.5. Investigating which parameters play a role in the location of earthquakes helps us to understand where future earthquakes are more likely to occur.

High-density macroseismic surveys were carried out in Rome, following both the 14 October 1997 ( M w = 5.6; I 0 = VIII MCS), and the 26 March 1998 ( M w = 5.3; I 0 = VI-VII MCS), Umbria-Marche (Central Italy) earthquakes. The aim of the work was to investigate the relationship between intensity and local geology for the modern urban area of Rome, an issue yet to be examined in contemporary times. Intensity data for both earthquakes were correlated with near-surface geology. The results show a statistically significant amplification in Holocene alluvial deposits. In particular, we found one degree of difference in MCS intensity between hard rock and recent alluvium formations. Furthermore, the distribution of the earthquake effects in alluvium-filled valleys was studied as a function of the shape ratio ( h/l ) of the alluvial basins. In the Tiber River valley, higher effects were found to be concentrated within a 150-m-wide band along the bedrock-alluvium surface contact, suggesting the occurrence of an edge effect. Surprising results were observed in the minor alluvial valleys, which displayed the strongest effects in terms of shaking. The narrowest valleys (width