Francesca Funiciello

Episodic Back-arc extension during restricted mantle convection in the Central Mediterranean

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

Mantle dynamics in the Mediterranean

The Mediterranean offers a unique opportunity to study the driving forces of tectonic deformation within a complex mobile belt. Lithospheric dynamics are affected by slab rollback and collision of two large, slowly moving plates, forcing fragments of continental and oceanic lithosphere to interact. This paper reviews the rich and growing set of constraints from geological reconstructions, geodetic data, and crustal and upper mantle heterogeneity imaged by structural seismology. We proceed to discuss a conceptual and quantitative framework for the causes of surface deformation. Exploring existing and newly developed tectonic and numerical geodynamic models, we illustrate the role of mantle convection on surface geology. A coherent picture emerges which can be outlined by two, almost symmetric, upper mantle convection cells. The downwellings are found in the center of the Mediterranean and are associated with the descent of the Tyrrhenian and the Hellenic slabs. During plate convergence, these slabs migrated backward with respect to the Eurasian upper plate, inducing a return flow of the asthenosphere from the backarc regions towards the subduction zones. This flow can be found at large distance from the subduction zones, and is at present expressed in two upwellings beneath Anatolia and eastern Iberia. This convection system provides an explanation for the general pattern of seismic anisotropy in the Mediterranean, first-order Anatolia and Adria microplate kinematics, and may contribute to the high elevation of scarcely deformed areas such as Anatolia and Eastern Iberia. More generally, the Mediterranean is an illustration of how upper mantle, small-scale convection leads to intraplate deformation and complex plate boundary reconfiguration at the westernmost terminus of the Tethyan collision.

Subduction Zone Geodynamics

Subduction Zone Geodynamics.- A Review of the Role of Subduction Dynamics for Regional and Global Plate Motions.- Subduction with Variations in Slab Buoyancy: Models and Application to the Banda and Apennine Systems.- Continental Collision and the STEP-wise Evolution of Convergent Plate Boundaries: From Structure to Dynamics.- Seismic Tomography and Anisotropy.- Seismic Anisotropy of Subduction Zone Minerals-Contribution of Hydrous Phases.- Local Earthquake Tomography in the Southern Tyrrhenian Region of Italy: Geophysical and Petrological Inferences on the Subducting Lithosphere.- Great Subduction Zone Earthquakes.- Effect of Subducting Seafloor Topography on the Rupture Characteristics of Great Subduction Zone Earthquakes.- Great Earthquakes in Slow-Subduction, Low-Taper Margins.- Seismogenic Zone Characterization.- Convergent Margin Structure in High-Quality Geophysical Images and Current Kinematic and Dynamic Models.- Imaging Interseismic Locking at the Nankai Subduction Zone, Southwest Japan.- Continental and Ridge Subduction Processes.- Exhumation Processes in Oceanic and Continental Subduction Contexts: A Review.- Evolution of Subductions Indicated by Melanges in Taiwan.- Subduction of an Active Spreading Ridge Beneath Southern South America: A Review of the Cenozoic Geological Records from the Andean Foreland, Central Patagonia (46-47 S).- Configuration of the Colombian Caribbean Margin: Constraints from 2D Seismic Reflection data and Potential Fields Interpretation.

Flat subduction dynamics and deformation of the South American plate: Insights from analog modeling

Received 14 June 2007; revised 13 January 2008; accepted 12 March 2008; published 21 June 2008. [1] We present lithospheric-scale analog models, investigating how the absolute plates’ motion and subduction of buoyant oceanic plateaus can affect both the kinematics and the geometry of subduction, possibly resulting in the appearance of flat slab segments, and how it changes the overriding plate tectonic regime. Experiments suggest that flat subductions only occur if a large amount of a buoyant slab segment is forced into subduction by kinematic boundary conditions, part of the buoyant plateau being incorporated in the steep part of the slab to balance the negative buoyancy of the dense oceanic slab. Slab flattening is a long-term process (� 10 Ma), which requires the subduction of hundreds of kilometers of buoyant plateau. The overriding plate shortening rate increases if the oceanic plateau is large enough to decrease the slab pull effect. Slab flattening increases the interplate friction force and results in migration of the shortening zone within the interior of the overriding plate. The increase of the overriding plate topography close to the trench results from (1) the buoyancy of the plate subducting at trench and (2) the overriding plate shortening. Experiments are compared to the South American active margin, where two major horizontal slab segments had formed since the Pliocene. Along the South American subduction zone, flat slab segments below Peru and central Chile/NW Argentina appeared at � 7 Ma following the beginning of buoyant slab segments’ subduction. In northern Ecuador and northern Chile, the process of slab flattening resulting from the Carnegie and Iquique ridges’ subductions, respectively, seems to be active but not completed. The formation of flat slab segments below South America from the Pliocene may explain the deceleration of the Nazca plate trenchward velocity. Citation: Espurt, N., F. Funiciello, J. Martinod, B. Guillaume, V. Regard, C. Faccenna, and S. Brusset (2008), Flat subduction dynamics and deformation of the South American plate: Insights

Why did Arabia separate from Africa? Insights from 3-D laboratory experiments

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

Dynamics of retreating slabs: 1. Insights from two‐dimensional numerical experiments

from the effects of mantle flux (part 2). Therefore, in this paper, we apply forces to the slab using simple analytical functions related to buoyancy and viscous forces in order to isolate the role of rheology on slab dynamics. We analyze parameters for simplified elastic, viscous, and nonlinear viscoelastoplastic single-layer models of slabs and compare them with a stratified thermomechanical viscoelastoplastic slab embedded in a thermal solution. The near-surface behavior of slabs is summarized by assessing the amplitude and wavelength of forebulge uplift for each rheology. In the complete thermomechanical solutions, vastly contrasting styles of slab dynamics and force balance are observed at top and bottom bends. However, we find that slab subduction can be modeled using simplified rheologies characterized by a narrow range of selected benchmark parameters. The best fit linear viscosity ranges between 5 � 10 22 Pa s and 5 � 10 23 Pa s. The closeness of the numerical solution to nature can be characterized by a Deborah number >0.5, indicating that elasticity is an important ingredient in subduction. INDEX TERMS: 8120 Tectonophysics: Dynamics of lithosphere and mantle—general; 8159 Tectonophysics: Rheology—crust and lithosphere; 8160 Tectonophysics: Rheology—general; KEYWORDS: subduction, numerical models, lithospheric rheology

Dynamical effects of subducting ridges: insights from 3-D laboratory models

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

Spreading pulses of the Tyrrhenian Sea during the narrowing of the Calabrian slab

The opening of the Tyrrhenian Sea has been punctuated by short-lived episodes of oceanic accretion on separate small backarc basins during early Pliocene (Vavilov basin) and early Pleistocene (Marsili basin) time. These spreading pulses are related to slab rollback and are synchronous with the reduction of the subduction zone width during the formation of the narrow Calabrian arc. Using laboratory models, we investigated the long-term and transient effects of the reduction of slab width on the subduction kinematics. We found that the abrupt reduction in slab width results in a pulse of acceleration of the trench retreat velocity, as the balance between driving and resisting forces acting on the slab is temporarily modified. Our findings also show that the time scale and amplitude of spreading observed in the Tyrrhenian Sea can be experimentally fitted if the scaled viscosity of the uppermost part of the mantle ranges between 10 19 and 10 20 Pa s.

The seismic cycle at subduction thrusts: 1. Insights from laboratory models

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

On the relation between trench migration, seafloor age, and the strength of the subducting lithosphere

Oceanic lithosphere thickens and strengthens as it grows older. Worldwide databases reveal that the age of the slab to a certain extent controls the subduction style. Old and thick (and consequently strong) slabs show a trench “advance,” while younger, thin (weak) slabs are migrating in retreating style (trench “rollback”). We performed numerical models to show that this configuration could be the result of the dynamic equilibrium between gravity and resisting forces. In particular we show that energy dissipation caused by bending and unbending of the slab, although less important than other resisting forces, could be a primary control on trench migration. Our results fit well with global compilations of kinematic data from modern subduction zones in two reference frames with different amounts of net rotations. Based on the age at which the transition from retreating to advancing style occurs, we propose an effective lithosphere/mantle viscosity of ~200 during bending of the lithosphere into the subduction zone.