Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
The association of experimental data showing that the plastic deformation of olivine, the main constituent of the upper mantle, is highly anisotropic and the ubiquitous seismic anisotropy in the upper mantle, which indicates that olivine crystals show coherent orientations over scales of tens to hundreds of kilometers, implies that the long-term deformation in the upper mantle is anisotropic. We propose a multiscale approach, based on a combination of finite element and homogenization techniques, to model the deformation of a lithospheric plate while fully considering the mechanical anisotropy stemming from a strain-induced orientation of olivine crystals in the mantle. This multiscale model explicitly takes into account the evolution of crystal preferred orientations (CPO) of olivine and of the mechanical anisotropy during the deformation. We performed a series of numerical experiments simulating the uniaxial extension of a homogeneous (100% olivine) but anisotropic plate to test the role of the olivine CPO on the plate mechanical behavior and the link between CPO and mechanical anisotropy evolution. Even for this simple solicitation, different orientations and intensity of the initial olivine CPO result in variable plate strengths and deformation regimes. A plate with an initial CPO where the olivine [100] and [010] axes are concentrated at 45 degrees to the extension direction has high resolved shear stresses on the easy (010)[100] and (001)[100] slip systems of olivine. This results in low strength and in deformation by transtension. Plates with an initial CPO where the maximum of [100] axes is parallel or normal to the extension direction show a high initial strength. Isotropic plates have an intermediate behavior. The progressive rotation of olivine [100] axes toward the imposed stretching direction results in hardening in all models, except in those characterized by an initial concentration of olivine [100] axes normal to the imposed extension, in which softening is followed by hardening.
ISSN: 1525-2027

hal-00420063
https://hal.archives-ouvertes.fr/hal-00420063
DOI : 10.1029/2009GC002423

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
1] Recent statistical analysis by Lallemand et al. (2008) of subduction zone parameters revealed that the back-arc deformation mode depends on the combination between the subducting (nu(sub)) and upper (nu(up)) plate velocities. No significant strain is recorded in the arc area if plate kinematics verifies nu(up) = 0.5 vsub - 2.3 (cm/a) in the HS3 reference frame. Arc spreading ( shortening) occurs if nu(up) is greater ( lower) than the preceding relationship. We test this statistical law with numerical models of subduction, by applying constant plate velocities far away from the subduction zone. The subducting lithosphere is free to deform at all depths. We quantify the force applied on the two converging plates to sustain constant surface velocities. The simulated rheology combined viscous (non-Newtonian) and brittle behaviors, and depends on water content. The influence of subduction rate vs is first studied for a fixed upper plate. After 950 km of convergence ( steady state slab pull), the transition from extensional to compressive stresses in the upper plate occurs for vs similar to 1.4 cm/a. The effect of upper plate velocity is then tested at constant subduction rate. Upper plate retreat ( advance) with respect to the trench increases extension ( compression) in the arc lithosphere and increases ( decreases) the subducting plate dip. Our modeling confirms the statistical kinematic relationship between vsub and nu(up) that describes the transition from extensional to compressive stresses in the arc lithosphere, even if the modeled law is shifted toward higher rates of upper plate retreat, using our set of physical parameters ( e. g., 100 km thick subducting oceanic plate) and short- term simulations. Our results make valid the choice of the HS3 reference frame for assessing plate velocity influence on arc tectonic regime. The subduction model suggests that friction along the interplate contact and the mantle Stokes reaction could be the two main forces competing against slab pull for upper mantle subductions. Besides, our simulations show that the arc deformation mode is strongly time dependent.
ISSN: 1525-2027

hal-00412328
https://hal.archives-ouvertes.fr/hal-00412328
https://hal.archives-ouvertes.fr/hal-00412328/document
https://hal.archives-ouvertes.fr/hal-00412328/file/arcayG32008.pdf
DOI : 10.1029/2007GC001875

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We use first first-principle methods and the Peierls-Nabarro model to evaluate the resistance to glide, characterized by the Peierls stress, of glide systems for end-member MgSiO3 Perovskite at mantle pressures. [010](100) is the easiest glide system in Mg-Perovskite at all pressures. Peierls stresses increase systematically with pressure for all systems except [001](010), indicating the importance of lattice friction at lower mantle pressures. The ratio of the maximum Peierls stress for each system relative to the [010](100) value defines their critical resolved shear stress (CRSS). These CRSS are used in a visco-plastic self-consistent homogenization model to predict the evolution of crystal preferred orientations (CPO) during deformation of polycrystalline Mg-Perovskite. In axial compression, [100] tends to align with the compression direction, in agreement with in situ observations in axial compression experiments. In simple shear, [010] concentrates near the shear direction and (100), although more dispersed, tends to align near the shear plane, consistent with the dominant activity of the easier [010](100) system. The calculated seismic anisotropy for a 100% Mg-Perovskite aggregate using the CPO in simple shear and the elastic constants Of MgSiO3 perovskite at lower mantle pressures and temperatures is weak (>3% for P-waves with and >2% for S-waves) and decreases with increasing temperature and pressure. P-waves show the fastest propagation parallel to the lineation and S-waves fast polarization is in the foliation at 38 GPa and normal to the lineation at 88 GPa. This weak anisotropy is consistent with global seismological observations of a nearly isotropic lower mantle. There are however two regions where strain-induced Mg-Perovskite CPO could contribute to anisotropy; a) low temperature regions in the uppermost lower mantle, where the predicted S-wave polarization anisotropy may attain 1.6% with a fast polarization parallel to the foliation, b) in high high-temperature domains in the D '' layer, where Mg-Perovskite may be the major stable phase, leading to polarization of fast S-waves normal to the lineation for propagation directions at high angle to the lineation and an apparent isotropy for all other propagation directions.
ISSN: 0012-821X

hal-00411460
https://hal.archives-ouvertes.fr/hal-00411460
DOI : 10.1016/j.epsl.2008.03.058

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Although subducting slabs undergo a bending deformation that resists tectonic plate motions, the magnitude of this resistance is not known because of poor constraints on slab strength. However, because slab bending slows the relatively rapid motions of oceanic plates, observed plate motions constrain the importance of bending. We estimated the slab pull force and the bending resistance globally for 207 subduction zone transects using new measurements of the bending curvature determined from slab seismicity. Predicting plate motions Using a global mantle flow model, we constrain the viscosity of the bending slab to be at most similar to 300 times more viscous than the upper mantle: stronger slabs are intolerably slowed by the bending deformation. Weaker slabs, however, cannot transmit a pull force sufficient to explain rapid trenchward plate motions unless slabs stretch faster than seismically observed rates of similar to 10(-15) s(-1). The constrained bending viscosity (similar to 2 x 10(23) Pa s) is larger than previous estimates that yielded similar or larger bending resistance (here similar to 25% of forces). This apparent discrepancy occurs because slabs bend more gently than previously thought, with an average radius of curvature of 390 km that permits subduction Of Strong slabs. This gentle bending may ultimately permit plate tectonics on Earth.
ISSN: 0012-821X

hal-00412383
https://hal.archives-ouvertes.fr/hal-00412383
DOI : 10.1016/j.epsl.2008.05.009

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
In this study, we propose a theoretical approach to test the validity of the isomechanical analogues for post-perovskite structures. Intrinsic plastic properties are evaluated for three materials exhibiting a post-perovskite phase: MgSiO3, MgGeO3 and CaIrO3. Dislocation properties of each structure are determined using the Peierls-Nabarro model based on first-principles calculations of generalised stacking fault and the plastic properties are extended to crystal-preferred orientations using a visco-plastic self-consistent method. This study provides intrinsic parameters of plastic deformation such as dislocation structures and Peierls stresses that can be directly compared between the three materials. It appears that it is very difficult to draw any simple conclusions on polycrystalline deformation simply by comparing single crystal properties. In particular, contrasting single crystal properties of MgSiO3 and CaIrO3 lead to similar crystal-preferred orientation of the polycrystal aggregates.
ISSN: 0031-9201

hal-00420104
https://hal.archives-ouvertes.fr/hal-00420104
DOI : 10.1016/j.pepi.2008.04.003

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
Upper-mantle flow beneath the South Pacific is investigated by analysing shear wave splitting parameters at eight permanent long-period and broad-band seismic stations and 10 broad-band stations deployed in French Polynesia from 2001 to 2005 in the framework of the Polynesian Lithosphere and Upper Mantle Experiment (PLUME). Despite the small number of eventsand the rather poor backazimuthal coverage due to the geographical distribution of the natural seismicity, upper-mantle seismic anisotropy has been detected at all stations except at Tahiti where two permanent stations with 15 yr of data show an apparent isotropy. The median value of fast polarization azimuths (N67.5°W) is parallel to the present Pacific absolute platemotion direction in French Polynesia (APM: N67°W). This suggests that the observed SKS fast polarization directions result mainly from olivine crystal preferred orientations produced by deformation in the sublithospheric mantle due to viscous entrainment by the moving Pacific Plate and preserved in the lithosphere as the plate cools. However, analysis of individualmeasurements highlights variations of splitting parameters with event backazimuth that imply an actual upper-mantle structure more complex than a single anisotropic layer with horizontal fast axis. A forward approach shows that a two-layer structure of anisotropy beneath French Polynesia better explains the splitting observations than a single anisotropic layer. Second ordervariations in the measurements may also indicate the presence of small-scale lateral heterogeneities. The influence of plumes or fracture zones within the studied area does not appear to dominate the large-scale anisotropy pattern but may explain these second-order splitting variations across the network.
ISSN: 0956-540X

hal-01249275
https://hal.archives-ouvertes.fr/hal-01249275
https://hal.archives-ouvertes.fr/hal-01249275/document
https://hal.archives-ouvertes.fr/hal-01249275/file/GJI2007SKS.pdf
DOI : 10.1111/j.1365-246X.2007.03475.x

Éditeur(s) : HAL CCSD American Geophysical Union (AGU)
Résumé : International audience
Nowadays, technological advances in satellite imagery measurements as well as the development of dense geodetic and seismologic networks allow for a detailed analysis of surface deformation associated with active fault seismic cycle. However, the study of earthquake dynamics faces several limiting factors related to the difficulty to access the deep source of earthquake and to integrate the characteristic time scales of deformation processes that extend from seconds to thousands of years. To overcome part of these limitations and better constrain the role and couplings between kinematic and mechanical parameters, we have developed a new experimental approach allowing for the simulation of strike-slip fault earthquakes and analyze in detail hundreds of successive seismic cycle. Model rheology is made of multilayered visco-elasto-plastic analog materials to account for the mechanical behavior of the upper and lower crust and to allow simulating brittle/ductile coupling, postseismic deformation phase and far-field stress transfers. The kinematic evolution of the model surface is monitored using an optical system, based on subpixel spectral correlation of high-resolution digital images. First, results show that the model succeed in reproducing the deformation mechanisms and surface kinematics associated to the main phases of the seismic cycle indicating that model scaling is satisfactory. These results are comforted by using numerical algorithms to study the strain and stress distribution at the surface and at depth, along the fault plane. Our analog modeling approach appears, then, as an efficient complementary approach to investigate earthquake dynamics.
ISSN: 0278-7407

hal-01174176
https://hal.archives-ouvertes.fr/hal-01174176
DOI : 10.1002/2014TC003701

Éditeur(s) : HAL CCSD American Geophysical Union (AGU)
Résumé : International audience
The Ronda peridotite supplies one of the best objects to document subcontinental mantle deformation, but its internal deformation and exhumation mechanisms remain controversial. Here we provide new structural data and numerical results that constrain the Oligocene-Miocene deformation history of the Ronda massif. We first describe a mantle shear zone in the northern massif that developed subcrustal strain localization during decompression and related partial melting. The deformation regime of this mantle shear zone evolved from penetrative NE-SW stretching to sinistral shear highlighted by discrete shear bands. We then show structural observations that document a viscous deformation in the southern massif occurring prior to the thrust-assisted emplacement of the peridotite during large decompression. Finally, we performed numerical investigations that quantify a high temperature of 980°C for the basal peridotite lens at the time of its crustal emplacement. Our numerical results constrain the timing of ductile deformation of the peridotite just before 22 Ma, probably between 30 and 22 Ma. Altogether, these features led us to conclude that the deformation and exhumation of the Ronda peridotite results from lithosphere thinning subsequently inverted in the course of the Oligocene-Miocene. Among available models, our findings support the hypothesis of peridotite exhumation by the inversion of a thinned back-arc continental lithosphere during westward slab rollback through the Alboran region.
ISSN: 0278-7407

insu-00924180
https://hal-insu.archives-ouvertes.fr/insu-00924180
https://hal-insu.archives-ouvertes.fr/insu-00924180/document
https://hal-insu.archives-ouvertes.fr/insu-00924180/file/tect20062.pdf
DOI : 10.1002/tect.20062

Éditeur(s) : HAL CCSD Nature Publishing Group
Résumé : International audience
Large earthquakes nucleate at tectonic plate boundaries, and their occurrence within a plate's interior remains rare and poorly documented, especially offshore. The two large earthquakes that struck the northeastern Indian Ocean on 11 April 2012 are an exception: they are the largest strike-slip events reported in historical times1, 2 and triggered large aftershocks worldwide3. Yet they occurred within an intra-oceanic setting along the fossil fabric of the extinct Wharton basin, rather than on a discrete plate boundary4, 5, 6, 7, 8. Here we show that the 11 April 2012 twin earthquakes are part of a continuing boost of the intraplate deformation between India and Australia that followed the Aceh 2004 and Nias 2005 megathrust earthquakes, subsequent to a stress transfer process recognized at other subduction zones9, 10. Using Coulomb stress change calculations, we show that the coseismic slips of the Aceh and Nias earthquakes can promote oceanic left-lateral strike-slip earthquakes on pre-existing meridian-aligned fault planes. We further show that persistent viscous relaxation in the asthenospheric mantle several years after the Aceh megathrust explains the time lag between the 2004 megathrust and the 2012 intraplate events. On a short timescale, the 2012 events provide new evidence for the interplay between megathrusts at the subduction interface and intraplate deformation offshore. On a longer geological timescale, the Australian plate, driven by slab-pull forces at the Sunda trench, is detaching from the Indian plate, which is subjected to resisting forces at the Himalayan front6, 8, 11.
ISSN: 0028-0836

hal-00767569
https://hal.archives-ouvertes.fr/hal-00767569
DOI : 10.1038/nature11520

Éditeur(s) : HAL CCSD European Geosciences Union
Résumé : The properties of the subduction interplate domain are likely to affect not only the seismogenic potential of the subduction area but also the overall subduction process, as it influences its viability. Numerical simulations are performed to model the long-term equilibrium state of the subduction interplate when the diving lithosphere interacts with both the overriding plate and the surrounding convective mantle. The thermomechanical model combines a non-Newtonian viscous rheology and a pseudo-brittle rheology. Rock strength here depends on depth, temperature and stress, for both oceanic crust and mantle rocks. I study the evolution through time of, on one hand, the brittle-ductile transition (BDT) depth, zBDT, and, on the other hand, of the kinematic decoupling depth, zdec, simulated along the subduction interplate. The results show that both a high friction and a low ductile strength at the asthenospheric wedge tip shallow zBDT. The influence of the weak material activation energy is of second order but not negligible. zBDT becomes dependent on the ductile strength increase with depth (activation volume) if the BDT occurs at the interplate decoupling depth. Regarding the interplate decoupling depth, it is shallowed (1) significantly if mantle viscosity at asthenospheric wedge tip is low, (2) if the difference in mantle and interplate activation energy is weak, and (3) if the activation volume is increased. Very low friction coefficients and/or low asthenospheric viscosities promote zBDT = zdec. I then present how the subducting lithosphere age affects the brittle-ductile transition depth and the kinematic decoupling depth in this model. Simulations show that a rheological model in which the respective activation energies of mantle and interplate material are too close hinders the mechanical decoupling at the down-dip extent of the interplate, and eventually jams the subduction process during incipient subduction of a young (20-Myr-old) and soft lithosphere under a thick upper plate. Finally, both the BDT depth and the decoupling depth are a function of the subducting plate age, but are not influenced in the same fashion: cool and old subducting plates deepen the BDT but shallow the interplate decoupling depth. Even if BDT and kinematic decoupling are intrinsically related to different mechanisms of deformation, this work shows that they are able to interact closely. Comparison between modelling results and observations suggests a minimum friction coefficient of 0.045 for the interplate plane, even 0.069 in some cases, to model realistic BDT depths. The modelled zdec is a bit deeper than suggested by geophysical observations. Eventually, the better way to improve the adjustment to observations may rely on a moderate to strong asthenosphere viscosity reduction in the metasomatised mantle wedge.
ISSN: 1869-9510

hal-00795559
https://hal.archives-ouvertes.fr/hal-00795559
DOI : 10.5194/se-3-467-2012, 2012

Éditeur(s) : HAL CCSD Elsevier
Résumé : 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.
ISSN: 0012-821X

hal-00406598
https://hal.archives-ouvertes.fr/hal-00406598
DOI : 10.1016/j.epsl.2007.02.004

Éditeur(s) : HAL CCSD
Résumé : International audience
Two opposing visions of the Aegean backarc tectonics implicitly contain contrasting images of the rheological behaviour of the lihosphere, supported by different sets of observations. The propagation of the NAF and extension in the Corinth Rift suggest a strongly localising rheology whereas the formation of MCC in the Cyclades and the Rhodope suggests instead a more viscous behaviour. This paradox was the seminal question addressed by the ANR-EGEO programme (2007-2010) that leads to the following conclusions: (1) The exhumation of the Cycladic MCC's is accommodated by the N-dipping North Cycladic Detachment System (NCDS) that partly reworks the Vardar suture and at the base by a series of thrusts, including the basal contact of the Cycladic Blueschists (CBS) over the Cycladic Basement (CB). The activity of LANF is due to the reactivation of pre-existing discontinuities such as thrusts or earlier detachments and much less to the interaction with granitic plutons. (2) Exhumation in the Cyclades has proceeded in two stages: (a) Eocene syn-orogenic exhumation within the subduction lower/upperplate interface while the post-orogenic Rhodope MCC formed further north, (b) Oligo-Miocene post-orogenic extension in an MCC mode coeval with syn-orogenic exhumation in Crete and the Peloponnese. (3) Both stages were associated with slab retreat, as early as the Eocene with an acceleration at 30-35 Ma. (4) The localisation of the presently active steeply-dipping normal faults on the southern margin of the Corinth Rift may have been preceded by a partial reworking of thrusts and syn-orogenic detachments by shallow-dipping decollements. (5) The localisation of the NAF in the Marmara Sea region and the Northern Aegean Sea is now accurately dated between 5.3 and 5 Ma Ma using the erosion and deposition surfaces that mark the Messinian salinity Crisis and nannofossils ages. (6) The formation of MCCs in the Cyclades resulted in the draining of the low viscosity lower crust from the northern Aegean Sea and led to an easier coupling between upper mantle and upper crust, easing the localisation and propagation of the NAF. (7) A comparison between crustal and mantle finite strain suggests that mantle flow due to slab retreat controls the stretching of the crust. Our model involves progressive slab retreat and mantle flow below the overriding plate inducing the observed succession of tectonic episodes with a progressive localisation of extension due to several slab tearing events. There is no paradox in the tectonic history and rheological behaviour of the Aegean lithosphere but an evolution of the kinematic boundary conditions due to slab deformation at depth.
American Geophysical Union
San Francisco, United States

insu-00860518
https://hal-insu.archives-ouvertes.fr/insu-00860518
BIBCODE : 2010AGUFM.T13G..04J

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
In order to estimate the seismic anisotropy of subducted crust, polycrystalline samples of KA1Si(3)O(8) K-lingunite (25% of the total subducted transformed sediments), have been synthesized and deformed under the temperature and pressure conditions of the subducted slabs. Transmission electron microscopy (TEM) characterizations of the recovered samples reveal that the microstructures are clearly dominated by [001] glide involving screw dislocations. For this reason, only {100} could be identified as glide planes, the question of [001] slip on {110} remains open. Few 1/2 < 111 > dislocations were observed gliding on {110} planes, which implies that 1/2 < 111 > {110} is a harder slip system than those involving [001] slip. The occurrence of sub-grain boundaries suggests that diffusion and climb might be active under these conditions.;To assess the texture of polycrystalline K-lingunite, the crystal preferred orientations (CP05) were calculated using visco-plastic self-consistent (VPSC) polycrystalline plasticity model in simple shear using the slip systems identified by TEM. Finally, the seismic properties of K-lingunite aggregates were calculated from the CPO and single crystal elasticity tensor. K-lingunite is predicted to have a high seismic anisotropy, which could combine constructively with one of the stishovite (same proportion as K-lingunite at the transition zone depth ranges).
ISSN: 0031-9201

hal-00534050
https://hal.archives-ouvertes.fr/hal-00534050
DOI : 10.1016/j.pepi.2010.06.013