Subduction dynamics as revealed by trench migration

International audienceNew estimates of trench migration rates allow us to address the dynamics of trench migration and back-arc strain. We show that trench migration is primarily controlled by the subducting plate velocity V-sub, which largely depends on its age at the trench. Using the hot and weak arc to back-arc region as a strain sensor, we define neutral arcs characterized by the absence of significant strain, meaning places where the forces (slab pull, bending, and anchoring) almost balance along the interface between the plates. We show that neutral subduction zones satisfy the kinematic relation between trench and subducting plate absolute motions: V-t = 0.5V(sub) - 2.3 (in cm a(-1)) in the HS3 reference frame. Deformation occurs when the velocity combination deviates from kinematic equilibrium. Balancing the torque components of the forces acting at the trench indicates that stiff (old) subducting plates facilitate trench advance by resisting bending

Neogene to Quaternary evolution of the Calabrian Subduction System, (Central Mediterranean)

We construct an ESE striking to WNW geological cross-section across the Calabrian Subduction System (Central Mediterranean) using seismic near vertical profiles and field data. The interpreted profiles were time-to-depth converted, merged and translated in a geological section stretching from the Marsili Oceanic Basin (Southern Tyrrhenian Sea) to the Ionian accretionary complex . Moving toward the east, the resulting section through the Paola, Amantea, and Crati basins, the Coastal Chain and Sila Massif and Crotone basin. The maximum elongation of these basins change progressively moving toward the east: from NNW in the Paola to NS in the Crati to the NNE in the Crotone basins. Data we present suggest that: Across the Calabria Tyrrhenian Continental Margin (CTCM), top of Kabilian-Calabrian Unit (KCU) is laterally variable in depth forming basins, which are separated by major structures with contractional or transcurrent kinematics, filled by Oligo-Miocene clastic to evaporitic deposits up to 1500m thick. Plio-Quaternary deposits display a remarkable variation in thickness from 4.5 km in the Paola Basin to less than 400m in the central sector of the margin. Plio-Quaternary sediments are internally sub-divisible into four sub-units (namely D1-D4) separated by tectonics enhanced angular unconformities. W-ward vergent reverse faults with limited vertical displacement offset the top of KCU as well as the Oligo-Miocene sedimentary and evaporitic units in the eastern side of the Paola basin. On land (Amantea \u2013 Crati) and farther to the east (Crotone basin) below a Messinian-Pleistocene deposits the top of KCU is variable structured and covered by a Oligo-Miocene clastic deposits with different thickness. The Plio-Quaternary deposits, unconformably overlay the Messinian and older deposits, show the maximum thickness in the Crotone basin. Two main tectonic unconformities within the Plio-Quaternary deposits have been recognised allowing the separation of this unit into three sub-units. In the offshore portion of the Crotone basin, SE-ward reverse faults dissect the KCU and the Oligo-Miocene up to the Messinian deposits. While the pre-Messinian tectonic history across the Calabrian Subduction System seems to be quite similar, a main reorganization of the system occurring during the (?) early and (?) middle-Pliocene. Geometrical and stratigraphic relationship show that several W-ward and E-ward vergent reverse faults in the Paola and Crotone basins, respectively, cut and offset Messinian evaporites and older sedimentary units, controlling the geometry of the basins. In the Paola Basin the amount of subsidence gradually increase during deposition of subunits D2 and D3, which are probably Pliocene in age. On land, the evidence of the unconformities in the Crotone basin indicate that Pliocene deposition occurring during the uplift of the Sila Massif. Therefore uplift of the Sila range occurred during the strong subsidence of the Paola and Crotone basins. The evolution of the overall structure can be then divided in two different steps: 1)the onset of subsidence started in the Late Miocene and covered a large areas presently occupied by the Paola and Crotone basins. This basin, was probably already separated into sub-basins but evolved in a slowly subsiding and poorly deformed area located between the active accretionary prism and the volcanic arc. Therefore in the Middle-Upper Miocene this basin could be defined as forearc basin. 2)In the Pliocene the structure of this large basin was fragmented due to the uplift of a central range (Sila Massif) with an overall pop-up like structure. 3)Uplift of the belt producing subsidence along the flanks and simultaneously formation of two distinct basins: the Paola and Crotone basins. This process probably occurred during episodes of fast roll-back of the subducting slab, as attested by the opening of two ocean floor basins in the back-arc region

Assessing the volcanic hazard for Rome. 40Ar/39Ar and In-SAR constraints on the most recent eruptive activity and present-day uplift at Colli Albani Volcanic District

We present new 40Ar/39Ar data which allow us to refine the recurrence time for the most recent eruptive activity occurred at Colli Albani Volcanic District (CAVD) and constrain its geographic area. Time elapsed since the last eruption (36 kyr) overruns the recurrence time (31 kyr) in the last 100 kyr. New interferometric synthetic aperture radar data, covering the years 1993–2010, reveal ongoing inflation with maximum uplift rates (>2 mm/yr) in the area hosting the most recent (<200 ka) vents, suggesting that the observed uplift might be caused by magma injection within the youngest plumbing system. Finally, we frame the present deformation within the structural pattern of the area of Rome, characterized by 50 m of regional uplift since 200 ka and by geologic evidence for a recent (<2000 years) switch of the local stress-field, highlighting that the precursors of a new phase of volcanic activity are likely occurring at the CAVD

The Dynamics of Forearc – Back‐Arc Basin Subsidence: Numerical Models and Observations From Mediterranean Subduction Zones

The subsidence history of forearc and back-arc basins reflects the relationship between subduction kinematics, mantle dynamics, magmatism, crustal tectonics, and surface processes. The distinct contributions of these processes to the topography variations of active margins during subduction initiation, oceanic subduction, and collision are less understood. We ran 2D elasto-visco-plastic numerical models including surface and hydration processes. The models show the evolution of wedge-top and retro-forearc basins on the continental overriding plate, separated by a forearc high. They are affected by repeated compression and extension phases. Compression-induced subsidence is recorded in the syncline structure of the retro-forearc basin from the onset of subduction. The 2-4 km upper plate negative residual topography is produced by the gradually steepening slab, which drags down the upper plate. Trench retreat leads to slab unbending and decreasing slab dip angle that leads to upper plate trench-ward tilting. Back-arc basins are either formed along inherited weak zones at a large distance from the arc or are created above the hydrated mantle wedge originating from arc rifting. Back-arc subsidence is primarily governed by crustal thinning that is controlled by slab roll-back and supported by the underlying mantle convection. High subduction and mantle convection velocities result in large wavelength negative dynamic topography. Collision and continental subduction are linked to the uplift of the forearc basins; however, the back-arc records ongoing extension during a soft collision. During the hard collision, both the forearc and back-arc basins are ultimately affected by the compression. Our modeling results provide insights into the evolution of Mediterranean subduction zones

Recent extension driven by mantle upwelling beneath the Admiralty Mountains (East Antarctica)

Northern Victoria Land is located at the boundary between an extended, presumably hot, region (West Antarctic Rift System) and the thick, possibly cold, East Antarctic craton. The style and timing of Tertiary deformation along with relationships with the magmatic activity are still unclear, and contrasting models have been proposed. We performed structural and morphotectonic analyses at the NE termination of northern Victoria Land in the Admiralty Mountains area, where the relationship between topography, tectonics, and magmatism is expected to be well pronounced. We found evidence of two subsequent episodes of faulting, occurring concurrently with the Neogene McMurdo volcanism. The first episode is associated with dextral transtension, and it is overprinted by extensional tectonics during the emplacement of large shield alkaline volcanoes. Upper mantle seismic tomography shows that the extensional regime is limited to regions overlying a low-velocity anomaly. We interpret this anomaly to be of thermal origin, and have tested the role of largescale upwelling on lithosphere deformation in the area. The results of this integrated analysis suggest that the morphotectonic setting of the region and the magmatism is likely the result of upwelling flow at the boundary between the cold cratonic and the hot stretched province (WARS), at work until recent time in this portion of the northern Victoria Land

Slab stiffness control of trench motion: Insights from numerical models

Subduction zones are not static features, but trenches retreat (roll back) or advance. Here, we investigate the dominant dynamic controls on trench migration by means of two- and three-dimensional numerical modeling of subduction. This investigation has been carried out by systematically varying the geometrical and rheological model parameters. Our viscoplastic models illustrate that advancing style subduction is promoted by a thick plate, a large viscosity ratio between plate and mantle, and a small density contrast between plate and mantle or an intermediate width (w ∼ 1300 km). Advancing slabs dissipate ∼45% to ∼50% of the energy in the system. Thin plates with relatively low viscosity or relatively high density, or wide slabs (w ∼ 2300 km), on the other hand, promote subduction in the retreating style (i.e., slab roll-back). The energy dissipated by a retreating slab is ∼35% to ∼40% of the total dissipated energy. Most of the energy dissipation occurs in the mantle to accommodate the slab motion, whereas the lithosphere dissipates the remaining part to bend and “unbend.” With a simple scaling law we illustrate that this complex combination of model parameters influencing trench migration can be reduced to a single one: plate stiffness. Stiffer slabs cause the trench to advance, whereas more flexible slabs lead to trench retreat. The reason for this is that all slabs will bend into the subduction zone because of their low plastic strength near the surface, but stiff slabs have more difficulty “unbending” at depth, when arriving at the 660-km discontinuity. Those bent slabs tend to cause the trench to advance. In a similar way, variation of the viscoplasticity parameters in the plate may change the style of subduction: a low value of friction coefficient weakens the plate and results in a retreating style, while higher values strengthen the plate and promote the advancing subduction style. Given the fact that also on Earth the oldest (and therefore probably stiffest) plates have the fastest advancing trenches, we hypothesize that the ability of slabs to unbend after subduction forms the dominant control on trench migration

Tectonics and seismicity of the Tindari Fault System, southern Italy: Crustal deformations at the transition between ongoing contractional and extensional domains located above the edge of a subducting slab

The Tindari Fault System (southern Tyrrhenian Sea, Italy) is a regional zone of brittle deformation located at the transition between ongoing contractional and extensional crustal compartments and lying above the western edge of a narrow subducting slab. Onshore structural data, an offshore seismic reflection profile, and earthquake data are analyzed to constrain the present geometry of the Tindari Fault System and its tectonic evolution since Neogene, including the present seismicity. Results show that this zone of deformation consists of a broad NNW trending system of faults including sets of right-lateral, left-lateral, and extensional faults as well as early strike-slip faults reworked under late extension. Earthquakes and other neotectonic data provide evidence that the Tindari Fault System is still active in the central and northern sectors and mostly accommodates extensional or rightlateral transtensional displacements on a diffuse array of faults. From these data, a multiphase tectonic history is inferred, including an early phase as a right-lateral strike-slip fault and a late extensional reworking under the influence of the subductionrelated processes, which have led to the formation of the Tyrrhenian back-arc basin. Within the present, regional, geodynamic context, the Tindari Fault System is interpreted as an ongoing accommodation zone between the adjacent contractional and extensional crustal compartments, these tectonic compartments relating to the complex processes of plate convergence occurring in the region. The Tindari Fault System might also be included in an incipient, oblique-extensional, transfer zone linking the ongoing contractional belts in the Calabrian-Ionian and southern Tyrrhenian compartments

Slab disruption, mantle circulation, and the opening of the Tyrrhenian basins

Plate tectonic history, geological, geochemical (element and isotope ratios), and seismological (P-wave tomography and SKS splitting) data are combined with laboratory modeling to present a three-dimensional reconstruction of the subduction history of the central Mediterranean subduction. We fi nd that the dynamic evolution of the Calabrian slab is characterized by a strong episodicity revealed also by the discrete opening of the Tyrrhenian Sea. The Calabrian slab has been progressively disrupted by means of mechanical and thermal erosion leading to the formation of large windows, both in the southern Tyrrhenian Sea and in the southern Apennines. Windows at lateral slab edges have caused a dramatic reorganization of mantle convection, permitting infl ow of subslab mantle material and causing a complicated pattern of magmatism in the Tyrrhenian region, with coexisting K- and Na-alkaline igneous rocks. Rapid, intermittent avalanches of large amounts of lithospheric material at slab edges progressively reduced the lateral length of the Calabrian slab to a narrow (200 km) slab plunging down into the mantle and enhancing the end of the subduction process

Shear-Velocity Structure and Dynamics Beneath the Sicily Channel and Surrounding Regions of the Central Mediterranean Inferred From Seismic Surface Waves

The evolution of the Sicily Channel Rift Zone (SCRZ) is thought to accommodate the regional tectonic stresses of the Calabrian subduction system. Much of the observations we have today are either limited to the surface or to the upper crust or deeper from regional seismic tomography, missing important details about the lithospheric structure and dynamics. It is unclear whether the rifting is passive from far-field extensional stresses or active from mantle upwelling beneath. We measure Rayleigh-and Love-wave phase velocities from ambient seismic noise and invert for 3-D shear-velocity and radial anisotropic models. Variations in crustal S-velocities coincide with topographic and tectonic features. The Tyrrhenian Sea has a ∼10 km thin crust, followed by the SCRZ (∼20 km). The thickest crust is beneath the Apennine-Maghrebian Mountains (∼55 km). Areas experiencing extension and intraplate volcanism have positive crustal radial anisotropy (VSH > VSV); areas experiencing compression and subduction-related volcanism have negative anisotropy. The crustal anisotropy across the Channel shows the extent of the extension. Beneath the Tyrrhenian Sea, we find very low sub-Moho S-velocities. In contrast, the SCRZ has a thin mantle lithosphere underlain by a low-velocity zone. The lithosphere-asthenosphere boundary rises from 60 km depth beneath Tunisia to ∼33 km beneath the SCRZ. Negative radial anisotropy in the upper mantle beneath the SCRZ is consistent with vertical mantle flow. We hypothesize a more active mantle upwelling beneath the rift than previously thought from an interplay between poloidal and toroidal fluxes related to the Calabrian slab, which in turn produces uplift at the surface and induces volcanism

Slab Driven Quaternary Rock‐Uplift and Topographic Evolution in the Northern‐Central Apennines From Linear Inversion of the Drainage System

Investigating rock-uplift variations in time and space provides insights into the processes driving mountain-belt evolution. The Apennine Mountains of Italy underwent substantial Quaternary rock uplift that shaped the present-day topography. Here, we present linear river-profile inversions for 28 catchments draining the eastern flank of the Northern-Central Apennines to reconstruct rock-uplift histories. We calibrated these results by estimating an erodibility coefficient (K) from incision rates and catchment-averaged erosion rates obtained from cosmogenic-nuclide data, and we tested whether a uniform or variable K produces a rock-uplift model that satisfactorily fits independent geochronological constraints. We employ a landscape-evolution model to demonstrate that our inversion results are reliable despite substantial seaward lengthening of the catchments during uplift. Our findings suggest that a rock-uplift pulse started around 3.0-2.5 Ma, coinciding with the onset of extension in the Apennines, and migrated southward at a rate of similar to 90 km/Myr. The highest reconstructed rock-uplift rates (>1 km/Myr) occur in the region encompassing the highest Apennine massifs. These results are consistent with numerical models and field evidence from other regions exhibiting rapid rock-uplift pulses and uplift migration related to slab break-off. Our results support the hypothesis of break-off of the Adria slab under the central Apennines and its southward propagation during the Quaternary. Moreover, the results suggest a renewed increase in rock-uplift rates after the Middle Pleistocene along the Adriatic coast, coeval with recent uplift acceleration along the eastern coast of southern Italy in the Apulian foreland