International audienceThe development of the Alpine mountain belt has been governed by the convergence of the African and European plates since the Late Cretaceous. During the Cenozoic, this orogeny was accompanied with two major kinds of intraplate deformation in the NW-European foreland: (1) the European Cenozoic Rift System (ECRIS), a left-lateral transtensional wrench zone striking NNE-SSW between the western Mediterranean Sea and the Bohemian Massif; (2) long-wavelength lithospheric folds striking NE and located between the Alpine front and the North Sea. The present-day geometry of the European crust comprises the signatures of these two events superimposed on all preceding ones. In order to better define the processes and causes of each event, we identify and separate their respective geometrical signatures on depth maps of the pre-Mesozoic basement and of the Moho. We derive the respective timing of rifting and folding from sedimentary accumulation curves computed for selected locations of the Upper Rhine Graben. From this geometrical and chronological separation, we infer that the ECRIS developed mostly from 37 to 17 Ma, in response to north-directed impingement of Adria into the European plate. Lithospheric folds developed between 17 and 0 Ma, after the azimuth of relative displacement between Adria and Europe turned counter-clockwise to NW SE. The geometry of these folds (wavelength = 270 km; amplitude = 1,500 m) is consistent with the geometry, as predicted by analogue and numerical models, of buckle folds produced by horizontal shortening of the whole lithosphere. The development of the folds resulted in ca. 1,000 m of rock uplift along the hinge lines of the anticlines (Burgundy Swabian Jura and Normandy Vogelsberg) and ca. 500 m of rock subsidence along the hinge line of the intervening syncline (Sologne Franconian Basin). The grabens of the ECRIS were tilted by the development of the folds, and their rift-related sedimentary infill was reduced on anticlines, while sedimentary accumulation was enhanced in synclines. We interpret the occurrence of Miocene volcanic activity and of topographic highs, and the basement and Moho configurations in the Vosges Black Forest area and in the Rhenish Massif as interference patterns between linear lithospheric anticlines and linear grabens, rather than as signatures of asthenospheric plumes

Our study compares the seismic properties between the flat and normal subduction regions in central Chile, to better understand the links between the slab geometry, surface deformation and the deeper structures. ln comparison with previous studies, we show the most complete 3-D regional seismic tomography images for this region, in which we use (1) a larger seis mic data set compiled from several short-term seismic catalogues, (2) a denser seismic array allowing a better resolution of the subduction zone from the trench to the backarc and into the upper ~30km of the slab and (3) a starting 1-D background velocity mode] specifically calculated for this region and refined over the years. We assess and discuss our tomography results using regional seismic attenuation models and estimating rock types on the basis of pressure and temperature conditions computed from thermomechanical models. Our results show significant seismic differences between the flat and normal subduction zones. As ex pected, the faster seismic velocities and increased seismicity within the flat slab and overriding lithosphere are generally consistent with a cooler thermal state. Our results are also consistent with dehydration of the mantle above the subducted Juan Fernandez Ridge at the eastern tip of the flat slab segment, indicating that the latter retains some fluids during subduction. However, fluids in the upper portion of the flat slab segment are not seismically detected, since we report instead fast slab seismic velocities which contradict the argument of its buoyancy being the cause of horizontal subduction. The forearc region, above the flat slab, exhibits high Vs and very low Vp/Vs ratios, uncorrelated with typical rock compositions, increased density or reduced temperature; this feature is possibly linked with the aftershock effects of the Mw7 .1 1997 Punitaqui earthquake, the flat slab geometry and/or seismic anisotropy. At the surface, the seismic variations correlate with the geological terranes. The Andean crust is strongly reduced in seismic velocities along the La Ramada-Aconcagua deformation belt, suggesting structural damage. Slow seismic velocities along the Andean Moho match non-eclogitized hydrated rocks, consistent with a previous delamination event or a felsic composition, which in turn supports the extent of the Chilenia terrane at these depths. We confirm previous studies that suggest that the Cuyania terrane in the backarc region is mafic and contains an eclogitized lower crust below 50-km depth. We also hypothesize major Andean basement detachment faults (or shear zones) to extend towards the plate interface and canalize slab-derived fluids into the continental crust.

SUMMARY The aim of this study is to quantify the relationship between the kinematics of subduction, deformation in the overriding plate and the evolution of slab geometry. A 2‐D finite element numerical code is used, and a first objective consists in benchmarking previously published analogue models. Far‐field plate velocities are applied, and once the subducting plate reaches the 660 km discontinuity, modelled as a rigid base, we obtain two different forms or styles of subduction that depend on the overriding plate velocity vop: if vop > 0, the slab lies forwards on the 660 km discontinuity (style 1), and if vop≤ 0, the slab lies backwards on the discontinuity (style 2). We also obtain a cyclic pattern with the slab folding on itself repeatedly when vsp > 0 and 2vop+vsp > 0 (where vsp is the subducting plate velocity). These conditions result from the analysis of several simulations in which the subduction velocities and plate viscosities are varied. When the slab periodically folds on the 660 km discontinuity, periods of shallow slab dip and compression in the overriding plate are followed by periods of slab steepening and relative extension in the overriding plate. Folding periodicity is controlled by the slab viscosity and subduction velocity. When a low‐viscosity zone is included in the overriding plate, the trench motion is effectively decoupled from the overriding plate velocity, therefore allowing it to be directly controlled by the deep dynamics of the slab. For the cyclic style 2 corresponding to forward folding of the slab, the low‐viscosity region in the overriding plate increases the stress amplitudes oscillations, the trench motion and the folding periodicity with time. Therefore the strength of the entire overriding plate is shown to directly control the dynamics of subduction. Using the Nazca and South American plate velocities we produce models of cyclic folding with a period of ca. 22 Ma and a minimal dip angle of ca. 10°. Episodic folding of the slab on the 660 km discontinuity would produce the necessary changes in slab dip and overriding plate deformation that explain episodes of volcanic quiescence alternating with greater rates of shortening along the Andes.

In this study, we support the buckling hypothesis by direct numerical modeling. We compare our results with the data on three most proeminent cases of the oceanic and continental folding-like deformation (Indian Ocean, Western Gobi (Central Asia) and Central Australia). We demonstrate that under reasonable tectonic stresses, folds can develop from brittle faults cutting through the brittle parts of a lithosphere. The predicted wavelengths and finite growth rates are in agreement with observations. We also show that within a continental lithosphere with thermal age greater than 400 My, either a bi-harmonic mode (two superimposed wavelengths, crustal and mantle one) or a coupled mode (mono-layer deformation) of inelastic folding can develop, depending on the strength and thickness of the lower crust.

[1] This paper presents an interdisciplinary study of the northern Chile double seismic zone. First, a highresolution velocity structure of the subducting Nazca plate has been obtained by the tomoDD doubledifference tomography method. The double seismic zone (DSZ) is observed between 80 and 140 km depth, and the two seismic planes is 20 km apart. Then, the chemical and petrologic characteristics of the oceanic lithosphere associated with this DSZ are deduced by using current thermal-petrological-seismological models and are compared to pressure-temperature conditions provided by a numerical thermomechanical model. Our results agree with the common hypothesis that seismicity in both upper and lower planes is related to fluid releases associated with metamorphic dehydration reactions. In the seismic upper plane located within the upper crust, these reactions would affect material of basaltic (MORB) composition and document different metamorphic reactions occurring within high-P (>2.4 GPa) and low-T (<570°C) jadeite-lawsonite blueschists and, at greater depth (>130 km), lawsonite-amphibole eclogite conditions. The lower plane lying in the oceanic mantle can be associated with serpentinite dehydration reactions. The Vp and Vs characteristics of the region in between both planes are consistent with a partially ($25-30 vol % antigorite, $0-10% vol % brucite, and $4-10 vol % chlorite) hydrated harzburgitic material. Discrepancies persist that we attribute to complexities inherent to heterogeneous structural compositions. While various geophysical indicators evidence particularly cold conditions in both the descending Nazca plate and the continental fore arc, thermomechanical models indicate that both seismic planes delimit the inner slab compressional zone around the 400°C (±50°C) isotherm. Lower plane earthquakes are predicted to occur in the slab's flexural neutral plane, where fluids released from surrounding metamorphic reactions could accumulate and trigger seismicity. Fluids migrating upward from the tensile zone below could be blocked in their ascension by the compressive zone above this plane, thus producing a sheeted layer of free fluids, or a serpentinized layer. Therefore earthquakes may present either downdip compression and downdip tensile characteristics. Numerical tests indicate that the slab's thermal structure is not the only factor that controls the occurrence of inner slab compression. (1) A weak ductile subduction channel and (2) a cold mantle fore arc both favor inner slab compression by facilitating transmission of compressional stresses from the continental lithosphere into the slab. (3) Decreasing the radius of curvature of the slab broadens the depth of inner slab compression, whereas (4) decreasing upper plate convergence diminishes its intensity. All these factors indicate that if DSZs indeed contour inner slab compression, they cannot be linked only to slab unbending, but also to the transmission of high compressional stresses from the upper plate into the slab.