Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Upon subduction, the oceanic crust transforms to blueschists and eclogites, with seismic properties that gradually become similar to those of the surrounding mantle. In order to evaluate the anisotropy of blueschists and glaucophane-bearing eclogites, the elastic constants of glaucophane single-crystal plates from the Sesia-Lanzo Zone (Aosta Valley, Western Alps) were measured using Brillouin spectroscopy at ambient conditions. The mean P- and S-wave velocities are 7.8 and 4.6 km s−1 respectively, and the anisotropy is high (38.1% (AVP) and 27.3% (AVS)). Glaucophane develops strong LPO, characterized by the [001]-axes concentrated sub-parallel to the lineation, and the {110} poles concentrated sub-perpendicular to the foliation in both blueschist and eclogite rocks. The measured LPO is in good agreement with viscoplastic selfconsistent numerical models. Seismic properties of glaucophane-bearing blueschist and eclogite are calculated by combining measured LPO and the single-crystal elastic moduli of glaucophane with the other main mineral constituents of the rock: mostly epidote for blueschist, and omphacite and garnet for eclogite. Blueschists present stronger anisotropies (AVP=16.1% and AVS=10.3%) than eclogites (AVP=2.9% and AVS=1.7%). The shear-wave splitting and resulting delay times for a 7-km thick layer of eclogite or blueschist are low for the eclogite (b0.03 s), but signi!cant for blueschist (0.16 s). Application to the subducted oceanic crust yields a decrease of velocity contrast with the surrounding mantle and of anisotropy at depth depending on the temperature of the slab. The details of the velocity variations in subducted oceanic crust are dif!cult to reconcile with the blueschist-eclogite transition as probed by exhumed metamorphic rocks and may require additional phases such as lawsonite or chemical variations such as oxidation state.
ISSN: 0040-1951

hal-00564518
https://hal.archives-ouvertes.fr/hal-00564518
DOI : 10.1016/j.tecto.2010.09.011

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
Seismic anisotropy has been widely studied this last decade, particularly by measuring splitting of vertically propagating core shear waves. The main interest in this technique is to characterize upper mantle flow beneath seismic stations. On the other hand, the major restriction in this method is that a single station gives a single anisotropy measurement. Alternative methods have been developed in order to avoid this restriction. An accurate determination of upper mantle seismic anisotropy beneath a seismic station may allow one, by doing anisotropy correction, to characterize remote or deeper anisotropy. The Geoscope network is ideal for this purpose because it is composed of a large set (about 26) of high-quality, broadband seismometers globally distributed and because some of these stations have run for more than 10 years and most of them for more than 5 years. We selected about 100 events at each site, generally of magnitude (m b) > 6.0, and we performed systematic measurements of the splitting parameters (fast polarization direction φ and delay time δt) on SKS, SKKS, and PKS phases. Splitting on oceanic islands has been difficult to observe owing to the low quality of the signal but also perhaps owing to complex upper mantle structures beneath the stations. Station KIP (Kipapa, Hawaii) in the Pacific is the only oceanic Geoscope station with a clear anisotropy. We determined well-constrained splitting parameters for 10 of the 17 continental stations that may be explained by a single anisotropic layer. The poor correlation between fast polarization directions and the absolute plate motion together with the apparent incoherence between the plate velocities and the observed delay times suggest that a simple drag-induced asthenospheric flow alone fails to explain most of the observations. For some stations located on or near major lithospheric structures (TAM, Tamanrasset, Algeria, for instance), we observe a good correlation between fast polarization directions and regional structures. At station SCZ (Santa Cruz, California), we found clear variations of the splitting parameters as a function of the event backazimuth, compatible with two layers of anisotropy. Three stations (CAN (Canberra), HYB (Hyderabad, India) and SSB (Saint Sauveur Badole, France)) seem to be devoid of detectable anisotropy.
ISSN: 2169-9313

hal-01388830
http://hal.univ-reunion.fr/hal-01388830
http://hal.univ-reunion.fr/hal-01388830/document
http://hal.univ-reunion.fr/hal-01388830/file/Barruol_anisoGeoscope_JGR1999.pdf
DOI : 10.1029/1999JB900033

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Mongolia represents the northernmost area affected by the India–Asia collision, and it is actively deformed along transpressive belts closely associated with large-scale strike-slip faults. The active and past mantle flow beneath this region is, however, poorly known. In order to investigate deep mantle deformation beneath central Mongolia and its relation with the surrounding major structures such as the Siberian craton, the Gobi–Altay belt and the Baikal rift, a NS-trending profile of broadband seismic stations has been deployed in the summer 2003 from the southern Siberian craton to the Gobi–Altay range, crossing the entire Hangay dome. Mantle flow is deduced from the splitting of teleseismic shear waves such as SKS phases. In eastern Mongolia, the permanent station ULN in Ulaanbaatar reveals the presence of two anisotropic layers, the upper one being oriented NE–SW, close to the trend of lithospheric structures and the lower one NW–SE, close to the trend of Eurasia absolute plate motion. Along the NS profile in central Mongolia, seismic anisotropy deduced from SKS splitting reveals a homogeneous NW–SE trending structure, fully consistent with the observations made in the Altay–Sayan in western Mongolia. The observed delay times of 1.5 to more than 2.0 s favor consistent mantle flow over large mantle thicknesses. Since the lithosphere is less than 100 km thick beneath central Mongolia and since the observed fast directions are parallel to the trend of the lithospheric structures but also close to the trend of the absolute plate motion, we propose that both the lithosphere and the asthenosphere may join their anisotropic effects beneath central Mongolia to explain the large delay times. Although GPS vectors represent the instantaneous displacement of the Earth's surface and SKS splitting the time and vertical integration of finite strain at depth, we use the opportunity of the dense geodetic measurements available in this region to discuss the anisotropy pattern in term of present-day deformation. In the Eurasia-fixed reference frame, GPS and SKS both depict a similar trend beneath central Mongolia, suggesting a lithospheric block “escaping” toward the east that could orient olivine a-axes in the upper mantle, within a transpressive tectonic regime. A different behaviour is observed in western Mongolia: the GPS vectors trend NS, close to the regional compression direction, whereas the fast SKS directions trend EW, suggesting a tectonic regime close to a mode of axial shortening, generating the development of an EW-trending foliation at depth. We therefore propose that Mongolia is a place where active and frozen lithospheric deformation may add their effects, together with the sublithospheric flow. In the three sources of anisotropy inferred, a fundamental role is played by the Siberian craton that acted as an undeformable core of the continent through time: the frozen Paleozoic and Mesozoic structures wrap around the craton, building up the fast anisotropic direction pattern; the present-day sublithospheric flow induced by the plate motion is likely deflecting around its deep roots; finally, the present-day tectonic regime appears to be controlled by the presence of the craton to the north.
ISSN: 0012-821X

hal-00407853
https://hal.archives-ouvertes.fr/hal-00407853
https://hal.archives-ouvertes.fr/hal-00407853/document
https://hal.archives-ouvertes.fr/hal-00407853/file/Barruol_SKS_Mongolia_EPSL2008_hal.pdf
DOI : 10.1016/j.epsl.2008.07.027

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : A differential effective medium (DEM) model is used to predict elastic properties for a set of porous and anisotropic aggregates, comprised of mixtures of calcite and muscovite. The DEM takes into consideration an anisotropic background medium with triclinic or higher symmetry, in which inclusions of idealized ellipsoidal shape are added incrementally. In general, the calculated elastic properties of a solid that contains inclusions representing "dry" pores/cracks are strongly dependent on the orientation and aspect ratio of the inclusions. Aspect ratios of inclusions in the synthetic aggregates, which consist of air-filled pores, are estimated from anisotropy of magnetic susceptibility (AMS) of samples whose pore space has been impregnated with a colloidal ferrofluid. The AMS derived pore shape geometry is used as an input value for inclusions in the DEM. Modeling results are compared with laboratory determined elastic properties, measured with ultrasonic waves. Calculated shear wave velocities agree in general well with laboratory measured S wave velocities, whereas calculated P wave velocities are typically 0.5-1.1 km/s higher than measured values. Differences between calculated and measured P wave velocities are attributed mainly to incomplete and biased ferrofluid saturation of pores. Spherical pores are preferably filled during imbibition, in comparison to thin cracks, which leads to overprediction of the calculated P wave velocities. The amount of ferrofluid that fills the pore space is dependent on the ratio of calcite to muscovite and the load used for compaction during sample synthesis. The permeability decreases with increasing muscovite content and increasing compaction load. Incomplete saturation of samples with high-muscovite content is confirmed by X-ray microtomography density contrast imaging of dry and ferrofluid saturated specimens.
ISSN: 0148-0227

hal-00617717
https://hal.archives-ouvertes.fr/hal-00617717
DOI : 10.1029/2010JB007750

É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 American Geophysical Union (AGU) AGU and the Geochemical Society
Résumé : International audience
To study how an impacting plume modifies the mantle lithosphere, we analyzed the microstructures and crystal preferred orientations (CPO) of 29 peridotites and 37 pyroxenites that sample the mantle root of the Ontong Java Plateau (OJP) from 60 to 120 km depth. The peridotites show a strong compositional variability, but homogeneous coarse granular to tabular microstructures, except for those equilibrated at the shallowest and deepest depths, which are porphyroclastic. All peridotites have clear olivine CPO, with dominant fiber-[010] patterns. Low intragranular misorientations and straight grain boundaries in olivine suggest that, above 100 km depth, annealing often followed deformation. Calculated density and P wave velocities of the peridotites decrease weakly with depth. S wave velocities decrease faster, resulting in increasing Vp/Vs ratio with depth. Calculated densities and seismic velocity profiles are consistent with those estimated for normal mantle compositions under a cold oceanic geotherm. Enrichment in pyroxenites may further increase seismic velocities. The calculated seismic properties cannot therefore explain the low S waves velocities predicted by Rayleigh wave tomography and ScS data in the mantle beneath the OJP. Calculated P and S waves anisotropy is variable (2–12%). It is higher on average in the deeper section of the lithosphere. Because olivine has dominantly [010]-fiber CPO patterns, if foliations are horizontal, vertically propagating S waves and Rayleigh waves will sample very weak anisotropy in the OJP mantle lithosphere. Moreover, if the orientation of the lineation changes with depth, the anisotropy-induced contrast in seismic properties might produce an intralithospheric reflector marking the stratification of the OJP mantle root.
ISSN: 1525-2027

hal-01116428
https://hal.archives-ouvertes.fr/hal-01116428
DOI : 10.1002/2014GC005452

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
Peridotite xenoliths brought up to the surface by the volcanism of the Kerguelen Islands represent a mantle that has been affected by a high degree of partial melting followed by intense melt percolation above the Kerguelen plume. These xenoliths are therefore particularly suitable to investigate effects of melt-rock interaction on crystallographic fabrics (lattice-preferred orientation (LPO)) of peridotite minerals and on the LPO-induced seismic properties of peridotites above a mantle plume. We have studied a suite of 16 ultramafic samples representative of different degrees of partial melting and magma-rock interaction among which the protogranular harzburgites are the least metasomatised xenoliths and dunites are the ultimate stage of metasomatism. Olivine LPO is characterized by high concentration of [010] axes perpendicular to the foliation and [100] axes close to the lineation or distributed in the foliation plane in harzburgites, whereas the high concentration of [100] axes is parallel to the lineation and [010] axes is perpendicular to the assumed foliation in dunites. Olivine LPO in harzburgites is interpreted as being due to a deformation regime in axial compression or transpression. The fabric strength of olivine decreases progressively from protogranular to poikilitic harzburgites and finally to dunites, for which it remains nevertheless significant (J index 3.8). Seismic properties calculated from LPO of minerals indicate that metasomatism at higher melt/rock ratio lowers the P wave velocities. The most significant difference between harzburgites and dunites corresponds to the distribution of S wave anisotropy. Harzburgites display the maximum of anisotropy within the foliation plane and the minimum of anisotropy perpendicular to the foliation plane, whereas the lowest anisotropy is parallel to the lineation for dunites. These modifications of seismic properties as a result of metasomatic processes may induce seismic heterogeneities in the mantle above the Kerguelen plume. In addition, assuming a lithospheric mantle primarily harzburgitic and structured with a horizontal foliation, the seismic properties calculated for the Kerguelen xenoliths reconcile the rather high anisotropy evidenced by the horizontally propagating surface waves with the apparent isotropy revealed by the absence of splitting of vertically propagating teleseismic SKS waves recorded by the GEOSCOPE Kerguelen station.
ISSN: 1525-2027

hal-00366222
https://hal.archives-ouvertes.fr/hal-00366222
https://hal.archives-ouvertes.fr/hal-00366222/document
https://hal.archives-ouvertes.fr/hal-00366222/file/bascou_Kerguelen_2007GC001879_hal.pdf
DOI : 10.1029/2007GC001879

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
In the southeastern corner of the Tyrrhenian basin, in the central Mediterranean Sea, a tight alignment of earthquakes along a well-defined Benioff zone marks one of the narrowest active trenches worldwide, where one of the last fragment of the former Tethys ocean is consumed. Seismic tomography furnishes snapshot images of the present-day position of this slab, and seismic anisotropy allows to reconstruct the past kinematics of the subduction process. Using seismic anisotropy fast directions as a proxy for the present and past mantle flow, we look backward for the seismic traces of the slab motion through the western-central Mediterranean mantle, from the starting locus of subduction toward its present day position. The result of combining independent data sets provides a coherent pattern of anisotropy that illustrates an example of slab rollback from its initiation point to its present-day position.
ISSN: 0012-821X

hal-01389115
http://hal.univ-reunion.fr/hal-01389115
http://hal.univ-reunion.fr/hal-01389115/document
http://hal.univ-reunion.fr/hal-01389115/file/Lucente_medit_EPSL2006_hal.pdf
DOI : 10.1016/j.epsl.2005.10.041

Éditeur(s) : HAL CCSD Elsevier
Résumé : The depth of the oceanic lithosphere–asthenosphere boundary (LAB), as inferred from shear wave velocities, increases with lithospheric age, in agreement with models of cooling oceanic lithosphere. On the other hand, the distribution of radial anisotropy under oceanic plates is almost age-independent. In particular, radial anisotropy shows a maximum positive gradient at a depth of ∼70 km, which, if used as a proxy, indicates an age-independent LAB depth. These contrasting observations have fueled a controversy on the seismological signature of the LAB. To better understand the discrepancy between these observations, we model the development of lattice preferred orientation (LPO) in upper mantle crystal aggregates and predict the seismic anisotropy produced by plate-driven mid-ocean ridge flows. The model accounts for the progressive cooling of the lithosphere with age and can incorporate both diffusion and dislocation creep deformation mechanisms. We find that an age-independent distribution of radial anisotropy is the natural consequence of these simple flows. The depth and strength of anisotropy is further controlled by the deformation regime – dislocation or diffusion creep – experienced by crystals during their ascent towards, and subsequent motion away from, the ridge axis. Comparison to surface wave tomography models yield constraints on rheological parameters such as the activation volume. Although not excluded, additional mechanisms proposed to explain some geophysical signatures of the LAB, such as the presence of partial melt or changes in water content, are not required to explain the radial anisotropy proxy. Our prediction, that the age-independent radial anisotropy proxy marks the transition to flow-induced asthenospheric anisotropy, provides a way to reconcile thermal, mechanical and seismological views of the LAB.
ISSN: 0012-821X

hal-01553243
https://hal.archives-ouvertes.fr/hal-01553243
DOI : 10.1016/j.epsl.2016.12.004

Éditeur(s) : HAL CCSD
Résumé : To understand the effect of inherited lithospheric structures on the continental plates deformation process, we have developed a multi scale model that explicitly takes into account an anisotropic viscosity due to the Crystal Preferred Orientation (CPO) of olivine in the mantle by coupling a viscoplastic self-consistent simulation of the deformation at the polycrystalline aggregate to a 3D Finite Element code. Using this tool, we calculate the mechanical behavior of a continental lithosphere submitted to extension for various initial structures. The first set of experiments investigates the deformation of a plate with a homogeneous initial olivine CPO. Results show that the macroscopic anisotropy is strongly dependent on the initial orientation of the olivine CPO relative to the solicitation and to its evolution. The second set of experiments aims to analyze the deformation of a continental plate containing pre-existing lithospheric scale shear zones. These multi-domain models show that the olivine mechanical anisotropy may induce the reactivation of pre-existing lithospheric scale shear zones in a plate submitted to rifting. In both cases, the results shows that the macroscopic mechanic anisotropy induced by the evolution of olivine CPO in the upper mantle is a key parameter to understand the deformation of the continental plates and to explain the observed reactivation of inherited structure from past tectonics events
Afin de comprendre l'effet des structures préexistantes sur la déformation des plaques continentales, nous avons développé une méthode de modélisation mécanique multi-échelle utilisant une viscosité anisotrope induite par les Orientations Préférentielles du Réseau (OPR) des cristaux d'olivine présents dans le manteau. Cette méthode en couple le calcul viscoplastique autocohérent de la déformation d'un agrégat d'olivine à un code éléments finis 3D. Elle a été utilisée pour le calcul du comportement mécanique d'une plaque continentale soumise à une extension. Une première série d'expériences numériques sur l'extension de plaque lithosphérique présentant des OPR initiales d'olivine homogènes a permis de montrer que l'anisotropie mécanique macroscopique est fortement dépendante de l'orientation initiale de l'OPR. La seconde série de simulations analyse l'extension d'une plaque continentale présentant une zone présentant des OPR d'olivine héritées d'un événement tectonique précédant. Les résultats de ces calculs multi-domaines que l'anisotropie mécanique de l'olivine peut expliquer la réactivation des zones de cisaillement héritées lors du rifting. Les résultats obtenus dans les deux séries de modèles montrent que l'anisotropie mécanique macroscopique induite par l'évolution des OPR d'olivine dans le manteau supérieur est un paramètre clé pour la compréhension de la déformation des plaques continentales et pour expliquer la réactivation observée de structures tectoniques, tels que les grands décrochements, lors des épisodes tectoniques ultérieurs.
https://tel.archives-ouvertes.fr/tel-00369463

tel-00369463
https://tel.archives-ouvertes.fr/tel-00369463
https://tel.archives-ouvertes.fr/tel-00369463/document
https://tel.archives-ouvertes.fr/tel-00369463/file/KnollThese09.pdf

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We investigated the deformation mechanisms and the seismic properties of 10 eclogite samples from different localities (Alps, Norway, Mali and eastern China) through the analysis of their microstructures and lattice-preferred orientations (LPO). These samples are representative of various types and intensity of deformation under eclogitic metamorphic conditions. Omphacite and garnet LPO were determined from electron backscatter diffraction (EBSD) technique. Garnet appears to be almost randomly oriented whereas omphacite develops strong LPO, characterized by the [001]-axes concentrated sub-parallel to the lineation, and the (010)-poles concentrated sub-perpendicular to the foliation. In order to analyze the deformation mechanisms that produced such omphacite LPO, we compare our observations to LPO simulated by viscoplastic self-consistent numerical models. A good fit to the measured LPO is obtained for models in which the dominant slip systems are 1/2h110i{11 ¯ 0}, [001] {110} and [001] (100). Dominant activation of these slip systems is in agreement with TEM studies of naturally deformed omphacite. Seismic properties of eclogite are calculated by combining the measured LPO and the single crystal elastic constants of omphacite and garnet. Although eclogite seismic anisotropies are very weak (less than 3% for both P-and S-wave), they are generally characterized by a maximum P-wave velocity sub-parallel to the lineation and by a minimum velocity approximately normal to foliation. The mean P-and S-wave velocities are high (respectively, 8.6 and 4.9 km/s). The S-wave anisotropy pattern displays complex relationships with the structural frame but the fast polarization plane generally tends to be parallel to the foliation. Calculated reflection coefficients show that an eclogite/crust interface is generally a good reflector (Rc > 0.1), whereas an eclogite body embedded in the upper mantle would be hardly detectable.
ISSN: 0040-1951

hal-01389715
http://hal.univ-reunion.fr/hal-01389715
http://hal.univ-reunion.fr/hal-01389715/document
http://hal.univ-reunion.fr/hal-01389715/file/Bascou_eclo-sismo_tphy2001_hal.pdf
DOI : 10.1016/S0040-1951(01)00156-1

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Seismic anisotropy originating within the continental crust is used to determine kinematic flow lines within active mountain belts and is widely attributed to regionally aligned mica. However, naturally deformed micaceous rocks commonly show composite (e.g. S-C) fabrics. It is necessary therefore to understand how both varying mica content and differing intensities of multiple foliations impact on seismic interpretations in terms of deformation fields. An outcrop analogue for granitic mid-crustal deformed zones is used here to calibrate the seismic response against both parameters. Seismic responses are modelled using crystallographic preferred orientations for polymineralic, micaceous granitic gneisses, measured using Electron Back-Scatter Diffraction. The sample results are generalised by modelling the effects of variations in modal composition and the relative importance of deformation fabrics of variable orientation, so-called rock and fabric recipes. The maximum P- and S-anisotropy are calculated at 16.6% and 23.9% for single-foliation gneisses but for mixed (i.e. S-C) foliation gneisses these values reduce to 5.8% and 7.5% respectively. Furthermore, mixtures of multiple foliations generate significant variations in the geometry of the seismic anisotropy. This effect, coupled with the geographical orientation of fabrics in nature. can generate substantial variations in the orientation and magnitude of seismic anisotropy (especially for shear waves) as measured for the continental crust using existing receiver function and teleseismic near-vertical incidence methods. Thus, maps of seismic anisotropy varying with depth in deforming continents need not imply necessarily depth-varying deformation kinematics and tectonic decoupling.
ISSN: 0012-821X

hal-00445247
https://hal.archives-ouvertes.fr/hal-00445247
DOI : 10.1016/j.epsl.2009.09.035

Éditeur(s) : HAL CCSD Springer Verlag (Germany)
Résumé : International audience
The influence of magnetic interactions on the anisotropy of magnetic susceptibility (AMS) have been largely studied by several theoretical models or experiments. Numerical models have shown that when magnetostatic interactions occur, the distributions of particles over the volume rather than their individual orientations control the AMS. We have shown recently from a comprehensive rock magnetic study and from a theoretical 2-dimensional (2-D) model that single domain particles closely packed in globule aggregates could produce strong local random interaction magnetic fields which could influence the magnetic susceptibility and decrease the degree of anisotropy. In this paper, we first present in detail this 2-D theoretical model and then we extend it to the 3-D case. The possible distribution function of the magnetostatic interaction fields comprises two extreme states: it is either isotropic or ordered. The former case corresponds to the thermal-demagnetized state while the second case corresponds to the alternating field (AF) demagnetized state. We show that when easy axes of magnetization are not uniformly distributed, the degree of anisotropy decreases as the interaction field increases in both AF- and thermal-demagnetized states in 2-D and 3-D geometry. Thus we conclude that random magnetic fields generated by a random arrangement of magnetic particles over the sample volume decrease the degree of anisotropy of AMS and may alter the magnetic fabric.
ISSN: 0039-3169

hal-00758113
https://hal.archives-ouvertes.fr/hal-00758113
https://hal.archives-ouvertes.fr/hal-00758113/document
https://hal.archives-ouvertes.fr/hal-00758113/file/SGEG-2012-HAL.pdf
DOI : 10.1007/s11200-010-9097-7

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The seismic anisotropy of the deep Earth is reviewed as a profile from the upper mantle to the solid inner core at the centre of the Earth. The upper mantle is by far the most anisotropic region of the Earth, followed the D” layer above the core mantle boundary. In contrast it is shown that many minerals that are present in the upper mantle, transition zone, lower mantle, D” layer and solid inner core are elastically anisotropic to different degrees. The anisotropy of hydrous phases present in subduction zones is briefly introduced. The basic concepts of single crystal and polycrystalline elasticity and the extrapolation of elastic properties to high temperature and pressure are presented. The specific features of elastic wave propagation in an anisotropic medium of any arbitrary symmetry are illustrated using mantle and inner core phases. The roles of crystal preferred orientation, water and melt in producing seismic anisotropy in the upper mantle are discussed.
Treatese of Geophysics, vol.2

hal-00408321
https://hal.archives-ouvertes.fr/hal-00408321
DOI : 10.1016/B978-044452748-6.00045-6

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Homogenization models are widely used to predict the evolution of texture (crystal preferred orientations) and resulting anisotropy of physical properties in metals, rocks, and ice. They fail, however, in predicting two main features of texture evolution in simple shear (the dominant deformation regime on Earth) for highly anisotropic crystals, like olivine: (1) the fast rotation of the CPO towards a stable position characterized by parallelism of the dominant slip system and the macroscopic shear and (2) the asymptotical evolution towards a constant intensity. To better predict CPO-induced anisotropy in the mantle, but limiting computational costs and use of poorly-constrained physical parameters, we modified a viscoplastic self-consistent code to simulate the effects of subgrain rotation recrystallization. To each crystal is associated a finite number of fragments (possible subgrains). Formation of a subgrain corresponds to introduction of a disorientation (relative to the parent) and resetting of the fragment strain and internal energy. The probability of formation of a subgrain is controlled by comparison between the local internal energy and the average value in the polycrystal. A two-level mechanical interaction scheme is applied for simulating the intracrystalline strain heterogeneity allowed by the formation of low-angle grain boundaries. Within a crystal, interactions between subgrains follow a constant stress scheme. The interactions between grains are simulated by a tangent viscoplastic self-consistent approach. This two-level approach better reproduces the evolution of olivine CPO in simple shear in experiments and nature. It also predicts a marked weakening at low shear strains, consistently with experimental data.
ISSN: 0012-821X

hal-01243503
https://hal.archives-ouvertes.fr/hal-01243503
DOI : 10.1016/j.epsl.2015.08.018

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
This study analyses the interactions between deformation and reactive fluid and melt percolation, and their effects on sub-arc mantle seismic properties based on microstructural observations on mantle xenoliths extracted by the Tubaf and Edison seamounts close to the Lihir Island, in the Papua New Guinea archipelago. These xenoliths sample an oceanic lithosphere, which has experienced high-temperature deformation in the presence of fluids or melts. This was followed by metasomatism under static conditions. Syn-kinematic percolation of reactive Si-rich melts or fluids in peridotites has produced pyroxene-enrichment, grain size reduction, and dispersion of olivine crystal preferred orientation (CPO). Fourier transform infrared spectroscopy analyses show that olivine has very low water contents (1-4 wt. ppm H2O), similar to spinel peridotites from other subduction zones. These low values may record both low water solubility in olivine at low pressure and dehydration during transport and exhumation. Water contents in pyroxenes are highly variable and likely result from spatially heterogeneous melt or fluid percolation. Analysis of olivine CPO indicates dominant activation of both (010)[100] and (001)[100] slip systems, which are characteristic of deformation under high temperature, low stress, low pressure and low to moderate hydrous conditions. Fast S-wave polarization and P- and Rayleigh propagation directions are thus parallel to the mantle flow direction. The pyroxene enrichment by melt-rock reactions is accompanied by dispersion of olivine CPO and induces a significant decrease of the maximum S-wave and P-wave anisotropy in the peridotites. The calculated seismic properties also show that the lowest Vp/Vs ratios (< 1.7) mapped in fore-arc mantle may only be explained by taking in consideration the CPO-induced elastic anisotropy of the peridotites.
ISSN: 0040-1951

hal-00950058
https://hal.archives-ouvertes.fr/hal-00950058
DOI : 10.1016/j.tecto.2013.09.024

Éditeur(s) : HAL CCSD
Résumé : This work provides new constraints on the interactions between deformation and percolation of melt or fluids, as well as on the implications of these processes for the seismic properties of the mantle wedge. It is based on the analysis of spinel peridotites from the massif of Ronda (Spain) and two xenolith suites from active subduction zones (Kamchatka, Papua New Guinea). A coupled structural, petrological and geochemical study of these samples shows that they underwent a reactive percolation of melts or Si-rich fluids, which was synchronous to a deformation event occurring under high temperature and low stress conditions, consistent with the PT conditions at the base of the lithosphere or in the asthenosphere. This reactive percolation is responsible for enrichment in pyroxenes, mainly orthopyroxene, which is often localized in bands parallel to the high-temperature foliation. This enrichment is associated with grain size reduction and dispersion of the crystallographic orientation of olivine. The dominant slip system in olivine is {0kl}[100], which results in fast S-wave polarization parallel to the flow direction in the mantle. The decrease in the intensity of the olivine crystal preferred orientations associated with the enrichment in pyroxene results in significant decrease of the anisotropy that may induce an underestimation of the thickness of the anisotropic layer by up to 33%. The observed orthopyroxene enrichment also lowers the Vp/Vs ratio, but cannot explain Vp/Vs < 1.7 mapped locally in the fore-arc mantle in Japan and the Andes. Such low Vp/Vs ratios may however be explained by considering the intrinsic anisotropy of the peridotites, which is generally ignored in large-scale Vp/Vs ratio mapping of the mantle wedge. Infrared analyses show that olivine from the both xenolith collections contains less water than the theoretical saturation calculated for their estimated equilibrium temperature in the spinel stability field. These low water content are similar to those observed in spinel peridotites from other subduction zones and probably record both the low solubility of water in olivine at relatively low pressure and dehydration during exhumation of the xenoliths. These measured water contents as well as theoretical saturation estimations are not sufficient to change the dominant slip direction of in olivine from [100] to [001].
Ce travail apporte de nouvelles contraintes sur les interactions entre déformation et processus d'hydratation et de percolation de magma ou de fluides dans le manteau à l'aplomb d'une zone de subduction et leurs implications sur les propriétés sismiques dans le coin mantellique. Il se base sur l'analyse de péridotites à spinelle provenant du massif de Ronda (Espagne) et deux séries de xénolites issues de zones de subduction actives (Kamchatka, Papouasie-Nouvelle-Guinée). L'étude structurale, pétrologique et géochimique de ces échantillons montrent qu'ils ont subit une percolation réactive de magma ou de fluide synchrone d'une déformation de haute température basse contrainte cohérente avec les condition PT de la base de la lithosphère ou de l'asthénosphère. Cette percolation réactive est responsable d'un enrichissement en pyroxènes, qui est souvent localisé dans des bandes parallèles à la foliation. Cet enrichissement est associé à une dispersion de l'orientation cristallographique de l'olivine. Les systèmes de glissement dominant dans l'olivine sont dans tous les cas {0kl}[100], ce qui implique que la direction de polarisation rapide des ondes S dans la partie supérieure du coin mantellique est parallèle à la direction d'écoulement dans le manteau. La décroissance de l'intensité des OPR de l'olivine associée à l'enrichissement en pyroxènes a pour conséquence une réduction non négligeable de l'anisotropie sismique qui peut induire jusqu'à 33% d'erreur sur l'interprétation de l'épaisseur de la couche anisotrope. Un enrichissement en orthopyroxène peut aussi entraîner une baisse du rapport Vp/Vs, mais ne peut expliquer les rapports Vp/Vs <1,7 cartés dans certains avant-arcs. Cependant de telles valeurs peuvent être expliquées si l'anisotropie des péridotites du coin mantellique est prise en compte. Les analyses par spectroscopie infra-rouge montrent que les olivines des deux séries de xénolites étudiées contiennent moins d'eau que la saturation théorique calculée pour leurs températures d'équilibre dans le domaine du spinelle. Ces faibles concentrations en eau sont similaires à celles observées dans les olivines des péridotites à spinelle des autres zones de subduction. Elles enregistrent probablement la faible solubilité de l'eau dans l'olivine à des pressions relativement basses et la déshydratation au cours de l'exhumation des xénolites. Ces teneurs en eau mesurées ainsi que les estimations de la saturation théorique ne sont pas suffisantes pour changer la direction de glissement dominante dans l'olivine de [100] à [001].
https://tel.archives-ouvertes.fr/tel-00689832

tel-00689832
https://tel.archives-ouvertes.fr/tel-00689832
https://tel.archives-ouvertes.fr/tel-00689832/document
https://tel.archives-ouvertes.fr/tel-00689832/file/These-Soustelle2010.pdf

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
We investigate upper mantle anisotropy beneath the Pyrenean range along three N-S profiles across the mountain belt. The results of a first profile that operated in 1993 in the central part of the belt have been presented elsewhere. We present the results of two other profiles that ran in 1995-1996 and 1996-1997 in the eastern and western part of the belt, respectively and propose an interpretation of the whole results. Teleseismic shear waves (SKS, SKKS, and PKS) are used to determine splitting parameters: the fast polarization direction φ and the delay time δt. Teleseismic shear wave splitting in the eastern Pyrenees displays homogeneous φ values trending N100°E and δt values in the range 1.1 to 1.5 s. A station located in the southern Massif Central, 100 km north of the range, is characterized by different splitting parameters (φ = N70°E, δt = 0.7 s). In the western part of the belt, anisotropy parameters are similar across the whole belt (φ = N110°E and δt = 1.3 to 1.5 s). Most of the measured delay times, including those obtained in the central part of the range, are above the global average of the SKS splitting (around 1 s). At the belt scale, φ is generally poorly correlated with recent estimations of the absolute plate motion, which predicts a fast direction ranging between N50°E and N80°E. Instead, the orientation of φ (N100°E) is parallel to the trend of the Pyrenean belt but also to Hercynian preexisting structures. This parallelism supports an anisotropy primarily related to frozen or active lithospheric structures. We show that a signature related to the Pyrenean orogeny is likely for the stations located in the internal domains of the belt. By contrast, the anisotropy measured at the stations located on the external parts of the belt could reflect a pre-Pyrenean (Hercynian) deformation. We suggest that a late Hercynian strike-slip deformation is responsible for this frozen upper mantle anisotropy and that the Pyrenean tectonic fabric developped parallel to this preexisting fabric. Finally, no particularly strong splitting is related to the North Pyrenean Fault, commonly believed to represent the plate boundary between Iberia and Eurasia.
ISSN: 2169-9313

hal-01388841
http://hal.univ-reunion.fr/hal-01388841
http://hal.univ-reunion.fr/hal-01388841/document
http://hal.univ-reunion.fr/hal-01388841/file/Barruol_Aniso_Pyr_JGR1998.pdf
DOI : 10.1029/98JB02790

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We present an integrated study of geochemistry, petrofabrics and seismic properties of strongly sheared eclogites from the Chinese Continental Scientific Drilling (CCSD) project in the Sulu ultrahigh-pressure (UHP) metamorphic terrane, eastern China. First, geochemical data characterize diverse protoliths of the studied eclogites. The positive Eu- and Sr-anomalies, negative Nb anomaly and flat portion of heavy rare earth elements in coarse-grained rutile eclogites (samples B270 and B295) suggest a cumulate origin in the continental crust, whereas the negative Nb anomaly and enrichment of light rare earth elements in retrograde eclogites (samples B504, B15 and B19) imply an origin of continental basalts or island arc basalts. Second, P-wave velocities (V-p) of three typical eclogite samples were measured under confining pressures up to 500 MPa and temperatures to 700 degrees C. At 500 MPa and room temperature, the mean V-p reaches 8.50-8.53 km/s in samples B270 and B295 but drops to 7.86 km/s in sample B504, and the P-wave anisotropy changes from 1.7-2.7% to 5.5%, respectively. The pressure and temperature derivatives of V-p are larger in the retrograde eclogite than in fresh ones. Third, the electron backscatter diffraction (EBSD) measurements of the eclogites reveal random crystal preferred orientation (CPO) of garnet and pronounced CPO of omphacite, which is characterized by a strong concentration of [001]-axes sub-parallel to the lineation and of (010)-poles perpendicular to the foliation. The asymmetric CPO of omphacite in sample B270 recorded a top-to-the-south shear event during subduction of the Yangtze plate. The calculated fastest V-p is generally sub-parallel to the lineation, but a different deformation environment during exhumation could form second-order variations in omphacite CPO and affect the V-p distribution in eclogites (e.g., the fastest V-p is at similar to 35 degrees from the foliation in sample 13295). Comparison between measured and calculated seismic properties indicates that the CPO of omphacite controls the seismic anisotropy of eclogites at high pressure, and compositional layering and retrograde minerals will increase the anisotropy. Calculated P-wave velocities agree well with velocities measured at 500 MPa and room temperature for fresh eclogites, but much higher than those of retrograde eclogite. As a case study, the laboratory-derived V-p-P and V-p-T relationships were used to estimate P-wave velocities of eclogites and peridotites beneath the Western Superior Province, Canada. The results indicate that besides the fabric-induced anisotropy, the direction dependence of pressure and temperature derivatives of V-p can significantly increase seismic anisotropy of eclogites with depth, which results in eclogites being an important candidate for the seismic anisotropy in the upper mantle. Due to their very high density and velocity, garnet-rich eclogites within peridotite could be detected in seismic reflections in subduction zones.
ISSN: 0040-1951

hal-00463512
https://hal.archives-ouvertes.fr/hal-00463512
DOI : 10.1016/j.tecto.2008.09.027

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
A detailed rock-magnetic and archeointensity study was carried out on materials baked by a western Mexican artisan following traditional techniques to produce faithful reproductions of archeological pieces of the Michoacan region (Western Mesoamerica). The field strength at the site (41.0 +/- 0.51 mu T) was measured with a fluxgate magnetometer and the temperature of the furnace during the baking process was monitored continually by means of a thermocouple placed in the middle of the baking cavity.;Rock-magnetic experiments performed on the raw material (clay and paste) and on in situ prepared baked ceramics and bricks included measurement of thermomagnetic curves (susceptibility and strong-field magnetization versus temperature), first-order reversal curves (FORC), anisotropy of magnetic susceptibility (AMS) and anisotropy of thermoremanent magnetization (A-TRM). Magnetite and probably hematite are present in the samples as carriers of the remanence. Hysteresis ratios suggest that the samples fall in the pseudo-single-domain grain size region, which may indicate a mixture of multidomain and a significant amount of single-domain grains.;Ceramic pieces and brick fragments were subjected to the Thellier-Coe archeointensity method and to an alternative paleointensity experiment, with a TRIAXE magnetometer, in order to check whether they are faithful recorders of the local geomagnetic field strength.;Mean raw-intensity of sample M1 (pottery) overestimates a 7% the expected site intensity, while those corresponding to the brick samples (LQ1 and LQ2) underestimate it 15%. Brick sample LNQ shows a slightly lower intensity (7%), but agrees with the expected site intensity within the experimental uncertainty. The intensity retrieved from the volcanic fragment also included closely reproduces the expected intensity. After A-TRM and cooling-rate corrections, all mean raw values move closer to the expected intensity.;Measurement of temperatures at different parts inside the kiln (bottom and upper parts of both central and peripheral parts) revealed the existence of significant thermal gradients, similar to those observed in ovens from other localities. Different cooling rates are then expected in a single oven.;The scatter in the intensity determinations observed in this study, retrieved from pieces elaborated together in the same oven, could arise from this differentiated cooling rate within the oven and thus, to an inappropriate cooling rate correction in the archeointensity protocol. As this situation was probably reproduced in the baking of ancient ceramic artifacts, a better knowledge of the temperature distribution inside these types of kiln would be desirable in order to choose the appropriate cooling rate correction.
ISSN: 0031-9201

hal-00644559
https://hal.archives-ouvertes.fr/hal-00644559
DOI : 10.1016/j.pepi.2011.06.007