Éditeur(s) : HAL CCSD
Résumé : International audience
The strength of the lithospheric mantle controls the mode of lithosphere deformation: the presence or absence ofa high strength brittle mantle respectively leads to localized or distributed rifting. However, first order geophysicaldata question the existence of such a brittle mantle. Here we use 2-D finite-element large strain modelling toquantify the impact of a ductile localizing mantle – instead of brittle – in triggering continental rifting. As a novelty,the mantle rheology considers the effect of grain boundary sliding during strain-induced grain size reduction,which may promote a significant strength drop and subsequent strain localization at low mantle temperature (<700-800°C). Our results reveal that such ductile localizing mantle implies varying modes of continental rifting thatmainly depend on both the amount of weakening in the ductile mantle and the strength of the lower ductile crust. Amedium to strong lower crust implies coupling between the upper crust and ductile localizing mantle, yielding tonarrow continental rifting. In contrast, a weak lower crust implies decoupling between the upper crust and ductilelocalizing mantle, giving rise to a switch from distributed faulting at incipient strain to localized faulting at largestrain. Ductile strain localization in the lithospheric mantle is therefore sufficient to trigger continental rifting,although a critical amount of weakening is required. Such ductile localizing mantle provides a relevant geologicaland mechanical alternative to the brittle mantle. It moreover provides a wider variety of modes of upper crustalfaulting that are commonly observed in nature.
EGU General Assembly 2014
Vienne, Austria

insu-01565161
https://hal-insu.archives-ouvertes.fr/insu-01565161

É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
Analogue and numerical models have shown that the strength of the lithospheric mantle controls the mode of lithosphere deformation. In extension, the presence or absence of a high strength brittle mantle respectively leads to localized or distributed rifting. However, first order geophysical data question the existence of such a brittle mantle. Here we use 2-D finite-element large strain modelling to quantify the impact of a ductile localizing mantle - instead of brittle - in triggering continental rifting. As a novelty, the mantle rheology considers the effect of grain boundary sliding during strain-induced grain size reduction, which may promote a significant strength drop and subsequent strain localization at low mantle temperature (< 700-800 °C). Our results reveal that such ductile localizing mantle implies varying modes of continental rifting that mainly depend on both the amount of weakening in the ductile mantle and the strength of the lower ductile crust. A medium to strong lower crust implies coupling between the upper crust and ductile localizing mantle, yielding to narrow continental rifting. In contrast, a weak lower crust implies decoupling between the upper crust and ductile localizing mantle, giving rise to a switch from distributed faulting at incipient strain to localized faulting at large strain. Ductile strain localization in the lithospheric mantle is therefore sufficient to trigger continental rifting, although a critical amount of weakening is required. Such ductile localizing mantle provides a relevant geological and mechanical alternative to the brittle mantle. It moreover provides a wider variety of modes of upper crustal faulting that are commonly observed in nature.
ISSN: 0040-1951

Droits : info:eu-repo/semantics/OpenAccess
insu-00924254
https://hal-insu.archives-ouvertes.fr/insu-00924254
https://hal-insu.archives-ouvertes.fr/insu-00924254/document
https://hal-insu.archives-ouvertes.fr/insu-00924254/file/Gueydan_precigout_2014_Tectonophysics.pdf
DOI : 10.1016/j.tecto.2013.11.029

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
This paper aims to present an overview on the influence of rheological heterogeneity and mechanical anisotropy on the deformation of continents. After briefly recapping the concept of rheological stratification of the lithosphere, we discuss two specific issues: (1) as supported by a growing body of geophysical and geological observations, crust=mantle mechanical coupling is usually efficient, especially beneath major transcurrent faults which probably crosscut the lithosphere and root within the sublithospheric mantle; and (2) in most geodynamic environments, mechanical properties of the mantle govern the tectonic behaviour of the lithosphere. Lateral rheological heterogeneity of the continental lithosphere may result from various sources, with variations in geothermal gradient being the principal one. The oldest domains of continents, the cratonic nuclei, are characterized by a relatively cold, thick, and consequently stiff lithosphere. On the other hand, rifting may also modify the thermal structure of the lithosphere. Depending on the relative stretching of the crust and upper mantle, a stiff or a weak heterogeneity may develop. Observations from rift domains suggest that rifting usually results in a larger thinning of the lithospheric mantle than of the crust, and therefore tends to generate a weak heterogeneity. Numerical models show that during continental collision, the presence of both stiff and weak rheological heterogeneities significantly influences the large-scale deformation of the continental lithosphere. They especially favour the development of lithospheric-scale strike-slip faults, which allow strain to be transferred between the heterogeneities. An heterogeneous strain partition occurs: cratons largely escape deformation, and strain tends to localize within or at the boundary of the rift basins provided compressional deformation starts before the thermal heterogeneity induced by rifting are compensated. Seismic and electrical conductivity anisotropies consistently point towards the existence of a coherent fabric in the lithospheric mantle beneath continental domains. Analysis of naturally deformed peridotites, experimental deformations and numerical simulations suggest that this fabric is developed during orogenic events and subsequently frozen in the lithospheric mantle. Because the mechanical properties of single-crystal olivine are anisotropic, i.e. dependent on the orientation of the applied forces relative to the dominant slip systems, a pervasive fabric frozen in the mantle may induce a significant mechanical anisotropy of the whole lithospheric mantle. It is suggested that this mechanical anisotropy is the source of the so-called tectonic inheritance, i.e. the systematic reactivation of ancient tectonic directions; it may especially explain preferential rift propagation and continental break-up along pre-existing orogenic belts. Thus, the deformation of continents during orogenic events results from a trade-off between tectonic forces applied at plate boundaries, plate geometry, and the intrinsic properties (rheological heterogeneity and mechanical anisotropy) of the continental plates.
ISSN: 0040-1951

hal-01389719
http://hal.univ-reunion.fr/hal-01389719
http://hal.univ-reunion.fr/hal-01389719/document
http://hal.univ-reunion.fr/hal-01389719/file/vauchez_etal_rheology_1998_hal.pdf
DOI : 10.1016/S0040-1951(98)00137-1

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The origin of plate tectonics is one of the most fundamental issues in earth and planetary sciences. Laboratory experiments indicate that the viscosity of silicate rocks is so strongly temperature-dependent that the entire surface of the Earth should be one immobile rigid plate. The rheology of oceanic lithosphere is, however, still poorly understood, and there exist few constraints on the temperature dependency of viscosity on the field scale. Here we report a new kind of observational constraint based on the geoid along oceanic fracture zones. We identify a large number of conspicuous small-scale geoid anomalies, which cannot be explained by the standard evolution model of oceanic lithosphere, and estimate their source density perturbations using a new Bayesian inversion method. Our results suggest that they are caused most likely by small-scale convection involving temperature perturbations of ∼300 K±100 K∼300 K±100 K. Such thermal contrast requires the activation energy of mantle viscosity to be as low as View the MathML source100±50 kJmol−1 in case of diffusion creep, and View the MathML source225±112 kJmol−1 in case of dislocation creep, substantially reducing the thickness of the stiffest part of oceanic lithosphere. Oceanic lithosphere may thus be broken and bent much more easily than previously thought, facilitating the operation of plate tectonics.
ISSN: 0012-821X

hal-01467827
https://hal.archives-ouvertes.fr/hal-01467827
DOI : 10.1016/j.epsl.2016.05.027

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
Rifting in a cratonic lithosphere is strongly controlled by several interacting processes including crust/mantle rheology, magmatism, inherited structure and stress regime. In order to better understand how these physical parameters interact, a 2 yr long seismological experiment has been carried out in the North Tanzanian Divergence (NTD), at the southern tip of the eastern magmatic branch of the East African rift, where the southward-propagating continental rift is at its earliest stage. We analyse teleseismic data from 38 broad-band stations ca. 25 km spaced and present here results from their receiver function (RF) analysis. The crustal thickness and Vp/Vs ratio are retrieved over a ca. 200 × 200 km2 area encompassing the South Kenya magmatic rift, the NTD and the Ngorongoro-Kilimanjaro transverse volcanic chain. Cratonic nature of the lithosphere is clearly evinced through thick (up to ca. 40 km) homogeneous crust beneath the rift shoulders. Where rifting is present, Moho rises up to 27 km depth and the crust is strongly layered with clear velocity contrasts in the RF signal. The Vp/Vs ratio reaches its highest values (ca. 1.9) beneath volcanic edifices location and thinner crust, advocating for melting within the crust. We also clearly identify two major low-velocity zones (LVZs) within the NTD, one in the lower crust and the second in the upper part of the mantle. The first one starts at 15–18 km depth and correlates well with recent tomographic models. This LVZ does not always coexist with high Vp/Vs ratio, pleading for a supplementary source of velocity decrease, such as temperature or composition. At a greater depth of ca. 60 km, a mid-lithospheric discontinuity roughly mimics the step-like and symmetrically outward-dipping geometry of the Moho but with a more slanting direction (NE–SW) compared to the NS rift. By comparison with synthetic RF, we estimate the associated velocity reduction to be 8–9 per cent. We relate this interface to melt ponding, possibly favouring here deformation process such as grain-boundary sliding (EAGBS) due to lithospheric strain. Its geometry might have been controlled by inherited lithospheric fabrics and heterogeneous upper mantle structure. We evidence that crustal and mantle magmatic processes represent first order mechanisms to ease and locate the deformation during the first stage of a cratonic lithospheric breakup.
ISSN: 0956-540X

hal-01622078
https://hal.archives-ouvertes.fr/hal-01622078
https://hal.archives-ouvertes.fr/hal-01622078/document
https://hal.archives-ouvertes.fr/hal-01622078/file/plasmanGJI2017.pdf
DOI : 10.1093/gji/ggx177

É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é : 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
This study reports petrostructural observations in the pargasite and phlogopite-bearing Finero peridotite massif (Italian Western Alps), which suggest that the pervasive foliation in this massif was formed by deformation concomitant with percolation of hydrous Si-rich melts: (1) diffuse contacts, but systematic parallelism between the pyroxenitic layers and the foliation of the peridotite (2) strong shape and crystal preferred orientations (SPO and CPO), but subhedral or interstitial shapes and weak intracrystalline deformation of the hydrous phases, (3) CPO, but interstitial shapes of the pyroxenes, (4) very coarse olivine grain sizes, which are correlated to the olivine abundance, and (5) elongated shapes, but weak intracrystalline deformation, and extremely weak and highly variable CPO of olivine. The pervasive deformation of the Finero peridotite occurred therefore under conditions that allowed coexistence of H2O–CO2-bearing melts, pargasite, and spinel, that is, temperatures of 980–1080 °C and pressures <2 GPa. The petrostructural observations suggest that the presence of hydrous melts results in accommodation of large amounts of deformation by stress-controlled dissolution–precipitation and advective transport of matter by the melts and in fast grain boundary migration in olivine. By consequence, it produces significant rheological weakening. Water contents in olivine are <4 ppm wt., implying limited contribution of hydration of olivine to weakening. In addition, the analysis of protomylonites composing the external domains of the shear zones that overprint the pervasive foliation indicates that the transition to melt-free conditions results in enhanced contribution of dislocation creep to the deformation. The associated increase of the peridotites' strength leads to onset of strain localization. The latter is not correlated to the local abundance in pargasite or phlogopite, implying that crystallization of amphiboles or phlogopite, even at concentrations of 25 vol.%, does not produce rheological weakening in the upper mantle.
ISSN: 0012-821X

hal-01667879
https://hal.archives-ouvertes.fr/hal-01667879
DOI : 10.1016/j.epsl.2017.08.015

Éditeur(s) : HAL CCSD
Résumé : International audience
Understanding subduction rheology in both space and time has been a challenge since the advent of plate tectonics. We herein focus on "subduction infancy", that is the first ~1-5 My immediately following subduction nucleation, when a newly born slab penetrates into the upper plate mantle and heats up. The only remnants of this critical yet elusive geodynamic step are thin metamorphic soles, commonly found beneath pristine, 100-1000 km long portions of oceanic lithosphere emplaced on continents (i.e., ophiolites).Through the (i) worldwide compilation of pressure-temperature conditions of metamorphic sole formation augmented by pseudosection thermodynamic modeling, (ii) calculations of the viscosity of materials along the plate interface and (iii) generic numerical thermal models, we provide a conceptual model of dynamic plate interface processes during subduction infancy (and initiation s.l.).We show in particular how major rheological switches across the subduction interface control slab penetration, and the formation of metamorphic soles. Due to the downward progression of hydration and weakening of the mantle wedge with cooling, the lower plate (basalt, sediment) and the upper plate (mantle wedge) rheologies equalize and switch over a restricted temperature-time-depth interval (e.g., at ~800°C and ~1 GPa, during 0.1-2 My, for high-temperature metamorphic sole formation). These switches result in episodes of maximum interplate mechanical coupling, thereby slicing the top of the slab and welding pieces (basalt, sediment) to the base of the mantle wedge. Similar mechanical processes likely apply for the later, deeper accretion and exhumation of high-temperature oceanic eclogites in serpentinite mélanges, or for the accretion of larger tectonic slices.This model provides constraints on the effective rheologies of the crust and mantle and general understanding, at both rock and plate scale, for accretion processes and early slab dynamics.
AGU Fall Meeting 2015
San Francisco, United States

insu-01242103
https://hal-insu.archives-ouvertes.fr/insu-01242103

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
We investigate flexural geometry and rheology of the India plate beneath the eastern Himalaya from a new gravity data set acquired in Bhutan. Compared to the well‒studied Nepal Himalaya, the obtained Bouguer anomaly profiles across the range show shorter wavelength flexure of the lithosphere with a narrower and shallower foreland basin. This new data set and seismic Moho constraints are interpreted together in terms of lithospheric flexure using a 2‒D thermomechanical model. Our results suggest that the strongest layer of the continental lithosphere beneath Bhutan is the upper mantle, as it is beneath Nepal. The observed west‒to‒east decrease in flexural wavelength is associated with weakening mantle rheology. The simulations show that this decrease can be related to ductile mantle behavior: either hydrated wet dunite or dry and hot olivine rheology. Both models display decoupled lithospheric layers leading to an eastward decrease of flexural rigidity from ∼1024 to ∼5*1022 N m in Nepal and Bhutan, respectively.
ISSN: 0094-8276

hal-00903564
https://hal.archives-ouvertes.fr/hal-00903564
DOI : 10.1002/grl.50793

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Subduction infancy corresponds to the first few million years following subduction initiation, when slabs start their descent into the mantle. It coincides with the transient (yet systematic) transfer of material from the top of the slab to the upper plate, as witnessed by metamorphic soles welded beneath obducted ophiolites. Combining structure–lithology–pressure–temperature–time data from metamorphic soles with flow laws derived from experimental rock mechanics, this study highlights two main successive rheological switches across the subduction interface (mantle wedge vs. basalts, then mantle wedge vs. sediments; at ∼800 °C and ∼600 °C, respectively), during which interplate mechanical coupling is maximized by the existence of transiently similar rheologies across the plate contact. We propose that these rheological switches hinder slab penetration and are responsible for slicing the top of the slab and welding crustal pieces (high- then low-temperature metamorphic soles) to the base of the mantle wedge during subduction infancy. This mechanism has implications for the rheological properties of the crust and mantle (and for transient episodes of accretion/exhumation of HP-LT rocks in mature subduction systems) and highlights the role of fluids in enabling subduction to overcome the early resistance to slab penetration.
ISSN: 0012-821X

hal-01383368
http://hal.upmc.fr/hal-01383368
http://hal.upmc.fr/hal-01383368/document
http://hal.upmc.fr/hal-01383368/file/Agard_2016_Plate_interface.pdf
DOI : 10.1016/j.epsl.2016.06.054

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
Numerical models show that small-scale convection (SSC) occurring atop a mantle plume is a plausible mechanism to rejuvenate the lithosphere. The triggering of SSC depends on the density contrast and on the rheology of the unstable layer underlying the stagnant upper part of the thermal boundary layer (TBL). Partial melting may change both properties. We analyze, using 2-D numerical simulations, how partial melting influences the dynamics of time-dependent SSC instabilities and the resulting thermo-mechanical rejuvenation of an oceanic plate moving atop of a plume. Our simulations show a complex behavior, with acceleration, no change, or delay of the SSC onset, due to competing effects of the latent heat of partial melting, which cools the plume material, and of the buoyancy increase associated with both melt retention and depletion of residue following melt extraction. The melt-induced viscosity reduction is too localized to affect significantly SSC dynamics. Faster SSC triggering is promoted for low melting degrees (low plume temperature anomalies, thick lithosphere, or fast moving plates), which limit both the temperature reduction due to latent heat of melting and the accumulation of depleted buoyant residue in the upper part of the unstable layer. In contrast, high partial melting degrees lead to a strong temperate decrease due to latent heat of melting and development of a thick depleted layer within the sublithospheric convecting layer, which delay the development of gravitational instabilities. Despite differences in SSC dynamics, the thinning of the lithosphere is not significantly enhanced relatively to simulations that neglect partial melting.
ISSN: 1525-2027

hal-01277175
https://hal.archives-ouvertes.fr/hal-01277175
DOI : 10.1002/2015GC005967

É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 Springer
Résumé : International audience
In this study, we used data recorded by two consecutive passive broadband deployments on the Gulf of Aden northern margin, Dhofar region, Sultanate of Oman. The objective of these deployments is to map the young eastern Gulf of Aden passive continental margin crust and upper mantle structure and rheology. In this study, we use shear-wave splitting analysis to map lateral variations of upper mantle anisotropy beneath the study area. In this study, we found splitting magnitudes to vary between 0.33 and 1.0 s delay times, averaging about 0.6 s for a total of 17 stations from both deployment periods. Results show distinct abrupt lateral anisotropy variation along the study area. Three anisotropy zones are identified: a western zone dominated by NW-SE anisotropy orientations, an eastern zone dominat- ed with NE-SW anisotropy orientations, and central zone with mixed anisotropy orientations similar to the east and west zones. We interpret these shorter wavelength anisotropy zones to possibly represent fossil lithospheric mantle anisotropy. We postulate that the central anisotropy zone may be representing a Proterozoic suture zone that separates two terranes to the east and west of it. The anisotropy zones west and east were being used indicative of different terranes with different upper mantle anisotropy signatures.
ISSN: 1866-7411

hal-00647478
https://hal.archives-ouvertes.fr/hal-00647478
DOI : 10.1007/s12517-011-0477-2

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The present study aims at better deciphering the different mechanisms involved in the functioning of the subduction interplate. A 2D thermo-mechanical model is used to simulate a subduction channel, made of oceanic crust, free to evolve. Convergence at constant rate is imposed under a 100 km thick upper plate. Pseudo-brittle and non-Newtonian behaviours are modelled. The influence of the subduction channel strength, parameterized by the difference in activation energy between crust and mantle () is investigated to examine in detail the variations in depth of the subduction plane down-dip extent, . First, simulations show that numerical resolution may be responsible for an artificial and significant shallowing of if the weak crustal layer is not correctly resolved. Second, if the age of the subducting plate is 100 Myr, subduction occurs for any . The stiffer the crust is, that is, the lower is, the shallower is (60 km depth if kJ/mol) and the hotter the fore-arc base is. Conversely, imposing a very weak subduction channel ( J/mol) leads there to an extreme mantle wedge cooling and inhibits mantle melting in wet conditions. Partial kinematic coupling at the fore-arc base occurs if kJ/mol. If the incoming plate is 20 Myr old, subduction can occur under the conditions that the crust is either stiff and denser than the mantle, or weak and buoyant. In the latter condition, cold crust plumes rise from the subduction channel and ascend through the upper lithosphere, triggering (1) partial kinematic coupling under the fore-arc, (2) fore-arc lithosphere cooling, and (3) partial or complete hindrance of wet mantle melting. then ranges from 50 to more than 250 km depth and is time-dependent if crust plumes form. Finally, subduction plane dynamics is intimately linked to the regime of subduction-induced corner flow. Two different intervals of are underlined: 80–120 kJ/mol to reproduce the range of slab surface temperature inferred from geothermometry, and 10–40 kJ/mol to reproduce the shallow hot mantle wedge core inferred from conditions of last equilibration of near-primary arc magmas and seismic tomographies. Therefore, an extra process controlling mantle wedge dynamics is needed to satisfy simultaneously the aforementioned observations. A mantle viscosity reduction, by a factor 4–20, caused by metasomatism in the mantle wedge is proposed. From these results, I conclude that the subduction channel down-dip extent, , should depend on the subduction setting, to be consistent with the observed variability of sub-arc depths of the subducting plate surface.
ISSN: 0031-9201

hal-01623155
https://hal.archives-ouvertes.fr/hal-01623155
DOI : 10.1016/j.pepi.2017.05.008

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The minerals constituting the Earth's upper mantle are nominally anhydrous silicates (NAMs). However they do contain hydrogen as a trace element, decorating point defects in their crystalline structure. Experimental petrology and mineralogy have quantified the maximum concentration under several compositional and thermodynamic conditions, but systematic studies on the hydrogen concentration in minerals from mantle-derived rocks have only recently been carried out. Here, we have compiled the distribution of hydrogen in upper mantle peridotite xenoliths, from which several conclusions can be drawn. NAMs from peridotite xenoliths contain a few ppm wt. H2O in their structure. From the current database, the hydrogen concentrations in olivine regularly increase with increasing depth. The amount of hydrogen in NAMs from peridotite xenoliths from subduction contexts is not higher than in other geological context for similar temperature and pressure conditions. The highest hydrogen concentrations is found in peridotitic olivines from cratonic mantle, and are likely due to the depth of origin. The increasing hydrogen concentration in olivine with increasing depth is likely controlled by the increase of H partitioning into olivine at the expense of orthopyroxene as imposed by a decrease in Al content in opx with depth. However, the sparse data could also indicate that the bulk hydrogen concentration slightly increases with depth > 150 km. In this case, it would suggest, locally (Udachnaya for example), a possible increase in water fugacity due to fluid saturation. Even if the most abundant mineral in mantle rocks is olivine, the bulk hydrogen concentration in peridotites is controlled by the amount of hydrogen stored in pyroxenes. However, hydrogen concentration in olivine remains crucial for consequences on physical properties such as rheology and electrical conductivity. Kinetics of hydrogen transport is reviewed and hydrous melt/fluid percolation appears necessary to homogenize the hydrogen distribution at km-scale. Sampling of natural rock specimens is currently biased (e.g., in favor of optically attractive samples) and needs improvement, which will be achieved by increasing sample diversity (all type of grain sizes, lithologies, geological settings) and size of rock samples, as well as advances in analytical techniques. Acquisition of high quality data will be achieved by studying the co-existing minerals in mantle specimens, exploring each sample by linear and mapping measurements, and using appropriate FTIR calibrations for polarized and unpolarized radiation.
ISSN: 0024-4937

hal-01306676
https://hal.archives-ouvertes.fr/hal-01306676
DOI : 10.1016/j.lithos.2015.11.012

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Deep crustal and mantle rocks are exhumed in core complex mode of extension in three types of structures: metamorphic core complexes, oceanic core complexes and magma poor passive margins. Using available analogue and numerical models and their comparison with natural examples, the present paper reviews the mechanical processes involved in these different types of extensional setting. Three main aspects are considered: i) the primary role of lithosphere rheology, ii) the lithosphere-scale patterns of progressive deformation that lead to the exhumation of deep metamorphic or mantle rocks and iii) the initiation and development of detachment zones. Crustal core complexes develop in continental lithospheres whose Moho temperature is higher than 750 °C with “upper crust-dominated” strength profiles. Contrary to what is commonly believed, it is argued from analogue and numerical models that detachments that accommodate exhumation of core complexes do not initiate at the onset of extension but in the course of progressive extension when the exhuming ductile crust reaches the surface. In models, convex upward detachments result from a rolling hinge process. Mantle core complexes develop in either the oceanic lithosphere, at slow and ultra-slow spreading ridges, or in continental lithospheres, whose initial Moho temperature is lower than 750 °C, with “sub-Moho mantle-dominated” strength profiles. It is argued that the mechanism of mantle exhumation at passive margins is a nearly symmetrical necking process at lithosphere scale without major and permanent detachment, except if strong strain localization could occur in the lithosphere mantle. Distributed crustal extension, by upper crust faulting above a décollement along the ductile crust increases toward the rift axis up to crustal breakup. Mantle rocks exhume in the zone of crustal breakup accommodated by conjugate mantle shear zones that migrate with the rift axis, during increasing extension.
ISSN: 0040-1951

insu-01610186
https://hal-insu.archives-ouvertes.fr/insu-01610186
https://hal-insu.archives-ouvertes.fr/insu-01610186/document
https://hal-insu.archives-ouvertes.fr/insu-01610186/file/brun-tect-2017.pdf
DOI : 10.1016/j.tecto.2017.09.017

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
We performed forward modelling of seismic anisotropy beneath a migrating strike-slip plate boundary to: (1) test if such a geodynamic context might explain teleseismic shear wave splitting data in the vicinity of the central part of the San Andreas Fault and (2) constrain the power of such data to unravel vertical and lateral variations in deformation patterns in the upper mantle. The modelling involves five steps: (1) thermomechanical modelling, using a finite-element code, of the deformation field, (2) viscoplastic self-consistent modelling of the resulting olivine and pyroxene crystal preferred orientations, (3) calculation of the elasticity tensors for different domains of the finite elements (FEs) model, (4) forward modelling of seismic wave propagation through the model using ray theory, finite-frequency theory and a full-wave approach and (5) performing splitting measurements on the synthetic seismograms. SKS splitting data in central California are best fitted by a model with a hotter geotherm within 60 km of the plate boundary accounting for the opening of an asthenospheric window due to the northward migration of the Mendocino Triple Junction. The westward motion of the plate boundary cannot however explain the rotation of fast polarization directions east of the San Andreas Fault in central California. Comparison between modelled and measured individual shear wave splitting also implies that the homogeneity of the two-layer models accounting for the observations in the vicinity of the San Andreas Fault indicates a sharp transition between lithospheric and asthenospheric deformations beneath this plate boundary. The ability of different synthetic approaches to localize horizontally and vertically the plate boundary-related deformation differs significantly. Splitting in ray theory synthetics closely follow variations in olivine crystallographic preferred orientations in the model. In contrast, splitting analysis on full-wave synthetics, which should be more representative of actual long-period SKS waves, results in smooth apparent lateral variations of anisotropy; the exact location and width of the plate boundary may only be retrieved by comparing fast-polarization profiles obtained using a multichannel analysis on waves with different periods.
ISSN: 0956-540X

hal-01236150
http://hal.univ-reunion.fr/hal-01236150
http://hal.univ-reunion.fr/hal-01236150/document
http://hal.univ-reunion.fr/hal-01236150/file/Bonnin_San_Andreas_modelling_GJI2012.pdf
DOI : 10.1111/j.1365-246X.2012.05650.x

Éditeur(s) : HAL CCSD Nature Publishing Group
Résumé : International audience
Mantle flow involves large strains of polymineral aggregates. The strongly anisotropic plastic response of each individual grain in the aggregate results from the interactions between neighbouring grains and the continuity of material displacement across the grain boundaries. Orthorhombic olivine, which is the dominant mineral phase of the Earth's upper mantle, does not exhibit enough slip systems to accommodate a general deformation state by intracrystalline slip without inducing damage. Here we show that a more general description of the deformation process that includes the motion of rotational defects referred to as disclinations can solve the olivine deformation paradox. We use high-resolution electron backscattering diffraction (EBSD) maps of deformed olivine aggregates to resolve the disclinations. The disclinations are found to decorate grain boundaries in olivine samples deformed experimentally and in nature. We present a disclination-based model of a high-angle tilt boundary in olivine, which demonstrates that an applied shear induces grain-boundary migration through disclination motion. This new approach clarifies grain-boundary-mediated plasticity in polycrystalline aggregates. By providing the missing mechanism for describing plastic flow in olivine, this work will permit multiscale modelling of the rheology of the upper mantle, from the atomic scale to the scale of the flow.
ISSN: 0028-0836

hal-01054315
https://hal.archives-ouvertes.fr/hal-01054315
DOI : 10.1038/nature13043