Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Laboratory models have been conducted to improve our understanding of the role that the resistance of the slab to bending and its coupling to the ambient mantle play in subduction dynamics over geological time scales. Our models are set up with a viscous plate of silicone (lithosphere) subducting under negative buoyancy in a viscous layer of glucose syrup (mantle). For our study, the lithosphere/upper mantle viscosity contrast has been systematically varied, from similar to 10 to similar to 10(5) in order to explore the parameter space between weak and strong slab dynamics. We found that subduction is characterized by a retreating mode for viscosity ratios > 10(4), by the coexistence of a retreating mode and an advancing mode for viscosity ratios between similar to 10(4) and similar to 10(2), and quasi-stationary, Rayleigh-Taylor like behaviour for ratios < 10(2). By combining our experimental results and kinematic data from current subduction zones in four reference frames which differ in the amount of net rotation, we infer that a lithosphere/upper mantle viscosity contrast of 150-500 is necessary to obtain realistic trench/subducting plate velocity ratios as well as the variability of subduction styles observed in nature.
ISSN: 0012-821X

hal-00412008
https://hal.archives-ouvertes.fr/hal-00412008
DOI : 10.1016/j.epsl.2008.04.006

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Transmission electron microscopy observations on olivine crystals deformed at moderate (≤1273 K) temperature evidenced dislocations interactions explaining the hardening observed in the experiments, but also recovery mechanisms by the absorption or emission of point defects. Thus we investigate the possibility that, at geological strain-rates, these recovery processes allow steady-state deformation by dislocation creep at low to moderate temperatures in the lithospheric mantle. We test this hypothesis using a 2.5-D dislocation dynamics (DD) model, which combines dislocation glide and recovery by climb. This model shows that diffusion-controlled recovery processes allow for steady-state deformation by dislocation creep in the lithospheric mantle at stresses <500 MPa. For stresses of 50–200 MPa, steady-state strain-rates of 10−15 s−110−15 s−1 may be attained at temperatures as low as 900 K. Fitting of the DD model produces a flow law, which represents a lower bound for the lithospheric mantle strength, since the models describe the deformation of an olivine single crystal in an easy slip orientation. Comparison of strain-rates and Moho temperatures inferred for different geodynamic environments and the predictions of this model-based flow law implies, nevertheless, that, except in incipient rifts, most of the observed deformation may be produced by stress levels ≤200 MPa, consistent with those inferred to be produced by convection. This convergence suggests that the present models, which explicitly calculate the time-dependent dislocation dynamics, may provide a correct first order estimate of the mechanical behaviour of the lithospheric mantle, which cannot be derived directly from any existing data.
ISSN: 0012-821X

hal-01277178
https://hal.archives-ouvertes.fr/hal-01277178
https://hal.archives-ouvertes.fr/hal-01277178/document
https://hal.archives-ouvertes.fr/hal-01277178/file/Boioli.pdf
DOI : 10.1016/j.epsl.2015.10.012

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
Two-dimensional simulations using a thermomechanical model based on a finite-difference method on a staggered grid and a marker in cell method are performed to study the plume-lithosphere interaction beneath moving plates. The plate and the convective mantle are modelled as a homogeneous peridotite with a Newtonian temperature- and pressure-dependent viscosity. A constant velocity, ranging from 5 to 12.5 cm yr−1, is imposed at the top of the plate. Plumes are generated by imposing a thermal anomaly of 150 to 350 K on a 50 km wide domain at the base of the model (700 km depth); the plate atop this thermal anomaly is 40 Myr old. We analyse (1) the kinematics of the plume as it impacts the moving plate, (2) the dynamics of time-dependent small-scale convection (SSC) instabilities developing in the low-viscosity layer formed by spreading of hot plume material at the base of the lithosphere and (3) the resulting thermal rejuvenation of the lithosphere. The spreading of the plume material at the base of the lithosphere, characterized by the ratio between the maximum down- and upstream horizontal (dimensionless) velocities in the plume-fed sublithospheric layer, Peup/Pedown depends on the ratio between the maximum plume upwelling velocity and the plate velocity, Peplume/Peplate. For fast plate velocities and sluggish plumes (low Peplume/Peplate), plate motion drags most plume material and downstream flow is dominant. As Peplume/Peplate increases, an increasing part of the plume material flows upstream. SSC systematically develops in the plume-fed sublithospheric layer, downstream from the plume. Onset time of SSC decreases with the Rayleigh number. For vigorous plumes, it does not depend on plate velocity. For more sluggish plumes, however, variations in the plume spreading behaviour at the base of the lithosphere result in a decrease in the onset time of SSCs with increasing plate velocity. In any case, SSC results in uplift of the isotherm 1573 K by up to 20 km relative to its initial equilibrium depth at the impact point.
ISSN: 0956-540X

hal-00856763
https://hal.archives-ouvertes.fr/hal-00856763
DOI : 10.1093/gji/ggt128

Éditeur(s) : HAL CCSD Geological Society of America
Résumé : International audience
We present the results of dynamic modeling of the western Mediterranean that accounts for observed global positioning system (GPS) surface deformation of the Alboran Sea and surrounding Rif and Betic Mountains as the result of the combined effects of relative motion of the Eurasian and Nubian plates, low strength in the Alboran Sea region and sub-lithospheric processes occurring beneath the External Rif domain. Assuming that the lithosphere behaves elastically over the short time period of the GPS observations, an elastic plate model is considered in our study, including an area of weak lithosphere (factor of 10) centered on the Alboran Sea and in which lateral boundary conditions consist of the Nubia-Eurasia oblique convergence. Sub-crustal processes are modeled by application of a horizontal traction on a small area (patch) at the base of the elastic plate. Our modeling studies demonstrate the need for sub-crustal or sub-lithospheric, southwestward-directed forcing to account for observed southwestward motion of the Rif and Betic domains. Based on the location, orientation, and small area of the traction patch, we hypothesize that forcing is associated with delamination and rollback of the subducted African continental lithospheric mantle beneath the External Rif zone, due to the pull of the oceanic part of the Western Mediterranean slab, a dynamic process that may be similar to that where the over-riding plate is driven toward the subduction zone during slab rollback.
ISSN: 0091-7613

hal-00524546
https://hal.archives-ouvertes.fr/hal-00524546
DOI : 10.1130/G30963.1

É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 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
As subducting plates reach the base of the upper mantle, some appear to flatten and stagnate, while others seemingly go through unimpeded. This variable resistance to slab sinking has been proposed to affect long-term thermal and chemical mantle circulation. A review of observational constraints and dynamic models highlights that neither the increase in viscosity between upper and lower mantle (likely by a factor 20–50) nor the coincident endothermic phase transition in the main mantle silicates (with a likely Clapeyron slope of –1 to –2 MPa/K) suffice to stagnate slabs. However, together the two provide enough resistance to temporarily stagnate subducting plates, if they subduct accompanied by significant trench retreat. Older, stronger plates are more capable of inducing trench retreat, explaining why backarc spreading and flat slabs tend to be associated with old-plate subduction. Slab viscosities that are ∼2 orders of magnitude higher than background mantle (effective yield stresses of 100–300 MPa) lead to similar styles of deformation as those revealed by seismic tomography and slab earthquakes. None of the current transition-zone slabs seem to have stagnated there more than 60 m.y. Since modeled slab destabilization takes more than 100 m.y., lower-mantle entry is apparently usually triggered (e.g., by changes in plate buoyancy). Many of the complex morphologies of lower-mantle slabs can be the result of sinking and subsequent deformation of originally stagnated slabs, which can retain flat morphologies in the top of the lower mantle, fold as they sink deeper, and eventually form bulky shapes in the deep mantle.
ISSN: 1553-040X

hal-01571965
https://hal.archives-ouvertes.fr/hal-01571965
DOI : 10.1130/GES01476.1

É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
Résumé : This work provides new constraints on the development and on the distribution of the deformation in the upper mantle and particularly beneath transform plate boundaries. USArray experiment and the remarkable increase of the dataset in California for the past ten years allowed us to scrutinize the lateral variations of the anisotropy in the vicinity of the San Andreas Fault zone. We have confirmed and increased the detection of two layers of anisotropy beneath this plate boundary. The first layer, located in the lithosphere, is related to the deformation induced at the fault, and the other one, located in the asthenosphere, is coherent with the anisotropy observed far from it, its origin is however less clear. We show that the deformation zone associated both to the San Andreas, Calaveras and Hayward Faults, is likely 40 km wide at 70 km depth. We then performed numerical thermomechanical modeling (ADELI) of the displacement of a transform plate boundary associated with the computation of the development of crystallographic fabrics using a viscoplastic selfconsistent approach (VPSC). We analyzed the distribution of the deformation in the model ant looked after the possible interactions at depth between deformation caused at surface by the strike-slip dynamic of the fault and the shearing at the base of the lithosphere caused by the horizontal displacement of the plates. Elastic properties derived from the crystallographic fabrics modeled, show that such interactions exist and induce, beneath the fault zone, a progressive rotation of the crystallographic fabrics with depth. Seismological signature of these smooth rotations is however not relevant with the presence of two anisotropic layers as proposed beneath California. We thus consider that a decoupling zone exists between the lithosphere and the asthenosphere beneath the California to account for the sharp separation between a lithospheric and an asthenospheric deformation. We furthermore estimate that anisotropy observed far form the San Andreas Fault in California cannot be explained only by the drag of the asthenosphere by the North America lithosphere as proposed in our article. Indeed, we can only expect few tenths of second of splitting delay from the anisotropic properties derived from the numerical modeling of a plate moving in the same direction and in the same velocity than the North American lithosphere only for 10 Ma of displacement. As delays observed in California rather reach 1.5 s, anisotropy in this region thus requires the existence of an active asthenospheric flow to be explained.
Le travail réalisé pendant cette thèse a permis d'apporter de nouvelles contraintes sur le développement et la distribution de la déformation dans le manteau supérieur et plus particulièrement au niveau des grandes limites de plaques décrochantes. Grâce à l'apport de l'expérience USArray et d'une dizaine d'années d'enregistrements sismologiques supplémentaires, nous avons pu étudier, de manière précise, les variations d'anisotropie dans le voisinage de la Faille de San Andreas. Nous avons confirmé et étendu l'observation de deux couches anisotropes sous cette limite de plaque. On y observe une première couche localisée dans la lithosphère marquant la déformation induite à la limite de plaque, et une autre, asthénosphérique, cohérente avec l'anisotropie observée loin de la faille et d'origine plus discutée. Nous avons montré que la zone de déformation associée aux failles de San Andreas, Calaveras et d'Hayward a, vraisemblablement, une largeur d'au moins 40 kilomètres en base de lithosphère, sous chacune de ces failles. Nous avons ensuite procédé à la modélisation thermomécanique (ADELI) de la migration d'une limite de plaques décrochante couplée à une modélisation du développement de fabriques cristallographiques par une approche viscoplastique auto-cohérente (VPSC). Ceci nous a permis d'y observer le développement de la déformation et les conséquences des possibles interactions entre la déformation décrochante en surface et le cisaillement en base de lithosphère dû au déplacement horizontal des plaques. Les propriétés élastiques déduites des fabriques cristallographiques modélisées montrent que de telles interactions existent et provoquent, sous la limite de plaques, une rotation des orientations cristallographiques avec la profondeur. Le signal associé à ces rotations progressives n'est toutefois pas cohérent avec la présence de deux couches d'anisotropie comme proposée sous la faille de San Andreas. Nous pensons par conséquent qu'il existe, sous la Californie, une zone de découplage entre la lithosphère et l'asthénosphère, permettant d'individualiser une déformation lithosphérique d'une déformation asthénosphérique. Nous estimons, en outre, que l'anisotropie observée dans l'asthénosphère sous la Californie ne peut être expliquée seulement par le cisaillement induit par le déplacement de la lithosphère Nord Amérique. En effet, les propriétés anisotropes obtenues par modélisation à partir d'une plaque se déplaçant dans une direction et une vitesse proche de celle de la plaque Amérique du Nord montrent qu'on ne peut espérer guère plus que quelques dixièmes de seconde de délai au bout de 10 Ma de déplacement. Les déphasages mesurés en Californie étant de l'ordre de 1,5 s, il est donc nécessaire d'invoquer la présence d'écoulements mantelliques actifs sous cette région.
https://tel.archives-ouvertes.fr/tel-00767124

tel-00767124
https://tel.archives-ouvertes.fr/tel-00767124
https://tel.archives-ouvertes.fr/tel-00767124/document
https://tel.archives-ouvertes.fr/tel-00767124/file/theseBonnin.pdf

Éditeur(s) : HAL CCSD Taylor & Francis
Résumé : Modern geochronology has moved beyond the acquisition of dates: the goal is to understand the significance of these numbers for the geodynamic evolution of Earth at all scales. The coupling of the laser-ablation microprobe (LAM) to inductively coupled plasma mass spectrometers (ICPMS, multicollector (MC)-ICPMS) has revolutionised geochronology and geochemistry over the last 10 years. These systems enable the rapid and precise in situ analysis of trace-element patterns and isotopic systems, while adding information related to microstructural context and major-element composition. The integration of these multiple sources of data is crucial in constraining the origin of the sample and the processes leading to its formation, so that we can understand the meaning of a date in terms of geological events. LAM-ICPMS measurement of U-Pb ages and trace-element patterns in zircon, coupled with LAM-MC-ICPMS analysis of Hf isotopes in the same grains, gives new insights into the processes of magma genesis. Applied to detrital zircons from modern drainages or sedimentary rocks (the TerraneChron approach), it becomes a powerful tool to investigate problems of crustal evolution on scales ranging from single terranes to continents. The in situ analysis (LAM-MC-ICPMS) of Re-Os systematics in single grains of sulfides in mantle-derived peridotites has demonstrated that most mantle rocks contain several generations of Os-bearing sulfides; whole-rock analyses are mixtures reflecting multiple melting and metasomatic events in the lithospheric mantle. These deep-seated events are commonly mirrored in the crust; Os model-age spectra from xenolith suites show age 'peaks' that correspond to the ages of thermal/tectonic events in the overlying crust, suggesting strong linkages between crust and mantle. Integrated studies of the timing and nature of crustal and mantle events, using these techniques, will be important for understanding the large-scale dynamics of the Earth.
ISSN: 0812-0099

hal-00411328
https://hal.archives-ouvertes.fr/hal-00411328
DOI : 10.1080/08120090802097450

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
We image the lithospheric and upper asthenospheric structure of western continental Yemen with 24 broadband stations to evaluate the role of the Afar plume on the evolution of the continental margin and its extent eastwards along the Gulf of Aden. We use teleseismic tomography to compute relative P-wave velocity variations in southwestern Yemen down to 300 km depth. Published receiver function analysis suggest a dramatic and localized thinning of the crust in the vicinity of the Red Sea and the Gulf of Aden, consistent with the velocity structure that we retrieve in our model. The mantle part of the model is dominated by the presence of a low velocity anomaly in which we infer partial melting just below thick Oligocene flood basalts and recent off-axis volcanic events (from 15 Ma to present). This low velocity anomaly could correspond to an abnormally hot mantle and could be responsible for dynamic topography and recent magmatism in western Yemen. Our new P-wave velocity model beneath western Yemen suggests the young rift flank volcanoes beneath margins and on the flanks of the Red Sea rift are caused by focused small-scale diapiric upwelling from a broad region of hot mantle beneath the area. Our work shows that relatively hot mantle, along with partial melting of the mantle, can persist beneath rifted margins after breakup has occurred.
ISSN: 1525-2027

hal-01015496
https://hal.archives-ouvertes.fr/hal-01015496
https://hal.archives-ouvertes.fr/hal-01015496/document
https://hal.archives-ouvertes.fr/hal-01015496/file/Korostelev_etal_GGG_2014.pdf
DOI : 10.1002/2014GC005316

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
The nature and extent of the regional lithosphere-asthenosphere interaction beneath Ireland and Britain remains unclear. Although it has been established that ancient Caledonian signatures pervade the lithosphere, tertiary structure related to the Iceland plume has been inferred to dominate the asthenosphere. To address this apparent contradiction in the literature, we image the 3-D lithospheric and deeper upper-mantle structure beneath Ireland via non-linear, iterative joint teleseismic-gravity inversion using data from the ISLE (Irish Seismic Lithospheric Experiment), ISUME (Irish Seismic Upper Mantle Experiment) and GRACE (Gravity Recovery and Climate Experiment) experiments. The inversion combines teleseismic relative arrival time residuals with the GRACE long wavelength satellite derived gravity anomaly by assuming a depth-dependent quasilinear velocity-density relationship. We argue that anomalies imaged at lithospheric depths probably reflect compositional contrasts, either due to terrane accretion associated with Iapetus Ocean closure, frozen decompressional melt that was generated by plate stretching during the opening of the north Atlantic Ocean, frozen Iceland plume related magmatic intrusions, or a combination thereof. The continuation of the anomalous structure across the lithosphere-asthenosphere boundary is interpreted as possibly reflecting sub-lithospheric small-scale convection initiated by the lithospheric compositional contrasts. Our hypothesis thus reconciles the disparity which exists between lithospheric and asthenospheric structure beneath this region of the north Atlantic rifted margin.
ISSN: 0956-540X

hal-00617815
https://hal.archives-ouvertes.fr/hal-00617815
DOI : 10.1111/j.1365-246X.2011.04921.x

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
This study presents microstructural evidence for deformation-controlled melt organization and for changes in olivine deformation associated with the presence of melt in an 80 m vertical section of the crust-mantle transition zone in the Oman ophiolite. This zone represents an 'end member' case for analyzing feedbacks between deformation and melt distribution in the upper mantle, since it experienced strong shear strains in presence of large melt fractions. It is characterized by a subhorizontal compositional layering at the mm to meter scale, from weakly impregnated dunites to olivine-rich gabbros, which parallels a pervasive foliation containing a strong stretching lineation. The parallelism between the compositional layering and the foliation, the diffuse limits of the layers, the alignment of elongated plagioclase-rich aggregates devoid of internal deformation structures with the elongation of olivine crystals in the dunitic layers, and the sharp compositional changes across some, but not all layer limits suggest deformation plays an essential role on the development of the layering. The variation on a mm-scale of the olivine crystal preferred orientation (CPO) symmetry as a function of the modal content: from axial-[100] symmetry in layers with <70% modal olivine to axial-[010] in more gabbroic levels (<40% olivine), which is repeated over the entire section, implies deformation in presence of variable melt fractions. Axial-[100] olivine CPO in olivine-rich layers is consistent with deformation by dislocation creep under high temperature, low pressure, dry conditions. Axial-[010] olivine CPO patterns imply additional sliding along preferentially wetted (010) grain boundaries, increase in the activity of [001] glide, or transpression localized in the melt-rich layers. Since the change in CPO symmetry is not accompanied by dispersion, instantaneous melt fractions must have remained <30-40% in all layers. The continuous variation in olivine CPO symmetry with decreasing olivine content implies therefore that the former depends on the cumulated strain in presence of melt rather than on the instantaneous melt fraction.
ISSN: 0012-821X

hal-00795508
https://hal.archives-ouvertes.fr/hal-00795508
DOI : 10.1016/j.epsl.2012.10.003

Éditeur(s) : HAL CCSD
Résumé : Volatile-rich, silica-undersaturated alkaline lavas record the important role of fluids during fluid-rock and magma-rock interactions in the mantle, which are key processes to understand the dynamics of the convective mantle and lithosphere-asthenosphere interactions in intracontinental settings. The aim of this thesis is to bring new constraints on the genesis of alkaline magmas by characterizing the crystallization conditions, the source and the partial melting processes taking part in the genesis of olivine nephelinites, pyroxene nephelinites and basanites from the Jbel Saghro volcanic field in the Moroccan Anti Atlas.The petrological and geochemical study of rocks and minerals coupled with the analysis of fluid inclusions constrains the pre-eruptive conditions of Saghro nephelinites to 1.7–2.2 GPa and ~1350 °C. Minerals show that nephelinitic magmas are rich in volatile elements (Cl, F, S), and fluid inclusions indicate that magmas were saturated with a CO2-rich fluid at pressures > 590 MPa. The various mineralogical assemblages and the presence of peridotite xenoliths suggest a rapid ascent for olivine nephelinites and more complex processes at depth for pyroxene nephelinites.Fractional crystallization and partial melting modelling of Saghro mafic lavas indicate that they are low-degree melts (0.6–2.5 %) of an amphibole-bearing carbonated peridotite enriched in incompatible elements, at the garnet-spinel transition (~80–85 km). Saghro nephelinites display a temporal evolution with a slight increase of the degree of melting and a decrease of the amount of residual amphibole from the oldest (olivine nephelinites, 9.6 Ma) to the most recent (pyroxene nephelinites, 2.9 Ma). Basanites form a system that is independent from nephelinites and are slightly higher-degree melts. Important variations in their chemical composition suggest variable amounts of fractional crystallization during ascent.The peculiar characteristics of Saghro nephelinites and basanites (enrichment in incompatible elements, negative anomalies in K, Zr, Hf and Ti, high Ca/Al and Zr/Hf ratios) indicate that their source was affected by carbonatitic metasomatism. The influence of this metasomatism is stronger for pyroxene nephelinites than for olivine nephelinites. These results suggest fluid-rock interactions beneath the Northwest African Craton, leading to the formation of a metasomatized mantle by CO2-rich carbonatitic components at the lithosphere-asthenosphere transition. The origin of the metasomatism inducing source enrichment and the formation of amphibole veins could be attributed to the melting of relict subducted oceanic lithosphere. The relatively low melting temperatures (< 1350 °C) suggest the absence of a thermal anomaly beneath the Jbel Saghro, and thus support a lithosphere delamination model as precursor of Saghro volca0,3nism. However, the increasing degree of partial melting over time, also observed in the Middle Atlas, together with the isotopic and geochemical similarities with Canary Islands alkaline lavas does not allow us to discard the influence of a deviation of the Canary mantle plume beneath northwest Africa.
Les laves alcalines sous-saturées riches en éléments volatils sont les marqueurs du rôle important des fluides dans le manteau et des interactions fluide-roche et magma-roche, processus clés pour comprendre la dynamique du manteau convectif et les interactions asthénosphère-lithosphère en domaine intracontinental. L’objectif de cette thèse est d’apporter de nouvelles contraintes sur la genèse des magmas alcalins en caractérisant les conditions de cristallisation, la source et les processus de fusion partielle à l’origine des néphélinites à olivine, des néphélinites à pyroxène et des basanites du champ volcanique du Jbel Saghro dans l’Anti-Atlas marocain.L’étude pétrologique et géochimique des roches et des minéraux, couplée à l’analyse des inclusions fluides a permis de contraindre les conditions pré-éruptives des néphélinites de Saghro à 1.7–2.2 GPa et ~1350 °C. Les minéraux montrent que les magmas néphélinitiques sont riches en éléments volatils (Cl, F, S), et les inclusions fluides indiquent que les magmas étaient saturés en fluide riche en CO2 à des pressions > 590 MPa. Les différents assemblages minéralogiques des néphélinites et la présence de xénolites péridotitiques suggèrent une ascension rapide des néphélinites à olivine et des processus plus complexes en profondeur pour les néphélinites à pyroxène.La modélisation des processus de cristallisation fractionnée et de fusion partielle des laves mafiques de Saghro a permis de déterminer qu’elles sont issues de faibles taux de fusion partielle (0.6–2.5 %) d’une péridotite carbonatée enrichie en éléments incompatibles, au niveau de la transition grenat–spinelle (~80–85 km) et en présence d’amphibole. Les néphélinites de Saghro montrent une évolution temporelle avec une légère augmentation du taux de fusion et une diminution de la quantité d’amphibole au résidu des plus anciennes (néphélinites à olivine, 9.6 Ma) aux plus récentes (néphélinites à pyroxène, 2.9 Ma). Les basanites forment un système indépendant des néphélinites et sont issues de taux de fusion plus élevés. Les fortes variations dans leur composition chimique suggèrent qu’elles ont subi de la cristallisation fractionnée lors de leur ascension.Les caractéristiques particulières des néphélinites et basanites de Saghro (enrichissement en éléments incompatibles, anomalies négatives en K, Zr, Hf et Ti, rapports Ca/Al et Zr/Hf élevés) indiquent que leur source a subi un métasomatisme principalement carbonatitique. L’influence de ce métasomatisme est plus forte pour les néphélinites à pyroxène que pour les néphélinites à olivine, impliquant une évolution temporelle de l’intensité du métasomatisme. Ces résultats suggèrent des interactions fluide-roche sous le craton Nord-Ouest Africain, entraînant la formation d'un manteau métasomatisé par des composants carbonatitiques riches en CO2 au niveau de la transition lithosphère-asthénosphère. L’origine du métasomatisme provoquant l’enrichissement de la source et la formation de veines d’amphibole pourrait être liée à la fusion de reliquats de croûte océanique subductée. Les températures de fusion relativement faibles (< 1350 °C) suggèrent l’absence d’anomalie thermique sous le Jbel Saghro, et favorisent donc un modèle de délamination de la lithosphère comme initiateur du volcanisme. Cependant, l’augmentation du taux de fusion partielle au cours du temps, également observée dans le Moyen Atlas, et les similitudes isotopiques et géochimiques avec les laves alcalines des îles Canaries ne permettent pas d’exclure une influence du panache des Canaries sur la source du volcanisme alcalin du Jbel Saghro.
https://tel.archives-ouvertes.fr/tel-01615888

NNT : 2015MONTS070
tel-01615888
https://tel.archives-ouvertes.fr/tel-01615888
https://tel.archives-ouvertes.fr/tel-01615888/document
https://tel.archives-ouvertes.fr/tel-01615888/file/53752_CHAMBOREDON_2015_archivage_cor.pdf

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : We study the Africa-Iberia plate boundary in the vicinity of Gibraltar. Numerous models have been proposed for that region throughout the last decades, proposing mechanisms that range widely from continental delamination, convective removal, to subduction of oceanic lithosphere. To better constrain upper-mantle structure under the region, we study waveforms of P-waves that traverse the Alboran Sea region between Spain and Morocco. These show dispersive behavior, which, together with early arrival times, confirms the presence of an anomalous upper mantle structure under the Alboran Sea. The dispersion is consistent with that expected from subducted lithosphere. Waveforms of body waves therefore provide a way to better constrain the elusive mantle structure and dynamics of the Alboran Sea region
ISSN: 0094-8276

hal-00407670
https://hal.archives-ouvertes.fr/hal-00407670
DOI : 10.1029/2007GL030964

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We discuss mantle structure and dynamics under the Gibraltar arc. For that region, a large variety of geodynamic models has been proposed, and there was, up to recently, no consensus on which of these models is the best one to describe the geodynamic situation of the region. The key to distinguishing the different models probably lies in the upper mantle, and we propose two specially adapted seismological techniques that help to distinguish each of the two models that are more frequently invoked. These are, on one hand, the subduction-rollback (SR) model, and on the other, continental delamination or convective removal (CDCR). The first technique that we use is based on observations of dispersed P-wavetrains that are expected for certain ray directions. Observations made on the Gibraltar arc suggest the presence of subducted oceanic lithosphere under the Alboran Sea. The second technique is based on shear-wave splitting of SKS phases, and it allows us to look for the radially-oriented flow that must be present in the mantle beneath the Alboran Sea, for a CDCR mechanism. There are however no indications of any radial flow up to this point. Both approaches therefore suggest that subduction-rollback is the best model to explain the complex geodynamic situation under the region.
ISSN: 0040-1951

hal-01388041
http://hal.univ-reunion.fr/hal-01388041
http://hal.univ-reunion.fr/hal-01388041/document
http://hal.univ-reunion.fr/hal-01388041/file/Bokelmann_Gibraltar_Tphy2011_hal.pdf
DOI : 10.1016/j.tecto.2010.08.004

É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 Oxford University Press (OUP)
Résumé : International audience
Most kimberlites contain abundant dunitic nodules. These are centimetre-sized, rounded and multi-grained assemblages of xenocrystic olivine with a wide range of compositions (Fo83 to Fo94). The absence of orthopyroxene and other mantle minerals and the range of olivine compositions have been attributed to reaction between mantle peridotite and (proto)kimberlitic fluid or melt, but the timing of the reaction is a subject of debate. In a kimberlite from the Kangamiut region of Greenland, nodule cores are surrounded by fine-grained outer margins with near-constant Fo contents (∼Fo88) but highly variable minor element contents (e.g. 500–2500 ppm Ni). These margins crystallized from the kimberlite melt and we show that their compositions can be explained by crystallization of olivine alone, if a high partition coefficient for Ni between melt and olivine (DNi > 20) is assumed. Orthopyroxene assimilation is not required, removing the constraint that its dissolution occurred during ascent of the kimberlite magma. Within some nodules, in addition to the usual core-to-margin gradients, we observe asymmetric compositional changes (variable Fo but near-constant minor element contents) across mobile grain boundaries. These changes document fluid percolation at the grain scale that occurred during dynamic recrystallization in the deforming lithospheric mantle. We note that chemical gradients associated with mobile grain boundaries are observed in olivines that cover the entire compositional range of the nodules, and propose that fluid-assisted dynamic recrystallization took place in dunite that was already compositionally heterogeneous. Reaction between peridotite and protokimberlitic melt or fluid and dissolution of orthopyroxene thus occurred within the lithospheric mantle, immediately (a few days) prior to the ascent of the kimberlite melt and the entrainment of the dunite nodules. We propose that the grain boundary zones probably mimic, at a fine scale, the fluid–peridotite interaction that caused, at a larger scale, orthopyroxene dissolution and formation of compositionally diverse olivine in kimberlites.
ISSN: 0022-3530

hal-01277963
https://hal.archives-ouvertes.fr/hal-01277963
DOI : 10.1093/petrology/egv054

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
In order to constrain the vertical and lateral extent of deformation and the interactions between lithosphere and asthenosphere in a context of a transpressional plate boundary, we performed teleseismic shear wave splitting measurements for 65 permanent and temporary broadband stations in central California. We present evidence for the presence of two anisotropic domains: (1) one with clear E–W trending fast directions and delay times in the range 1.5 to 2.0 s and (2) the other closely associated with the San Andreas Fault system with large azimuthal variations of the splitting parameters that can be modeled by two anisotropic layers. The upper of the two layers provides fast directions close to the strike of the main Californian faults and averaged delay times of 0.7 s; the lower layers show E–W directions and delay times in the range 1.5 to 2.5 s and thus can be compared to what is observed in stations that require a single layer. We propose the E–W trending anisotropic layer to be a 150 to 200 km thick asthenospheric layer explained by the shearing associated with the absolute plate motion of the North American lithosphere. The shallower anisotropic layer ought to be related to the dynamics of the San Andreas Fault system and thus characterized by a vertical foliation with lineation parallel to the strike of the faults localized in the lithosphere. We also propose that the anisotropic layer associated with each fault of the San Andreas Fault system is about 40 km wide at the base of the lithosphere.
ISSN: 0148-0227

hal-01236161
http://hal.univ-reunion.fr/hal-01236161
http://hal.univ-reunion.fr/hal-01236161/document
http://hal.univ-reunion.fr/hal-01236161/file/Bonnin_SKS_California_JGR2010.pdf
DOI : 10.1029/2009JB006438

Éditeur(s) : HAL CCSD Springer
Résumé : Subduction is a major process that plays a first-order role in the dynamics of the Earth. The sinking of cold lithosphere into the mantle is thought by many authors to be the most important source of energy for plates driving forces. It also deeply modifies the thermal and chemical structure of the mantle, producing arc volcanism and is responsible for the release of most of the seismic energy on Earth. There has been considerable achievements done during the past decades regarding the complex interactions between the various processes acting in subduction zones. This volume contains a collection of contributions that were presented in June 2007 in Montpellier (France) during a conference that gave a state of the art panorama and discussed the perspectives about "Subduction Zone Geodynamics". The papers included in this special volume offer a unique multidisciplinary picture of the recent research on subduction zones geodynamics. They are organized into five main topics: Subduction zone geodynamics, Seismic tomography and anisotropy, Great subduction zone earthquakes, Seismogenic zone characterization, Continental and ridge subduction processes. Each of the 13 papers collected in the present volume is primarily concerned with one of these topics. However, it is important to highlight that papers always treat more than one topic so that all are related lighting on different aspects of the complex and fascinating subduction zones geodynamics.
https://hal.archives-ouvertes.fr/hal-00808864

ISBN : 978-3-540-87974
hal-00808864
https://hal.archives-ouvertes.fr/hal-00808864