Éditeur(s) : HAL CCSD Elsevier
Résumé : Understanding chemical element mobility during high-pressure metamorphism is paramount to the knowledge of the transformations occurring during the course of eclogite-facies metamorphism. In particular, the role of deformation and fluid circulation appear essential. In order to better decipher this role, we investigated two eclogite-facies samples from the Lindas Nappe in the Bergen arcs (Norwegian Caledonides), which are related to a fluid/deformation metamorphic event. This study, based on detailed microtextural investigations, coupled with major and trace element analyses, either in situ EMPA and LA-ICP-MS or by conventional ICP-MS nebulisation on mineral fractions, illustrates the complexity of chemical reactions accompanying the eclogite- and amphibolite-facies metamorphism. A single eclogite-facies mineral species two eclogite-facies (i.e. phengite and epidote) depicts various chemical features depending on its textural location at thin-section scale. The correlation with the mineralogical reaction occurring in each textural site demonstrates that the compositional variations are inherited from the precursor minerals. Thus, the composition of eclogite-facies minerals is locally controlled by the composition of few hundred-micron domains. On the contrary, alteration phases of the amphibolite facies (symplectite and calcic amphibole) display clear enrichment in LILE, Pb, Sr and LREE compared to their precursor minerals (omphacite and garnet). This supports that the transport of elements at hand-sample scale was enhanced during the retrogression and elements were efficiently delocalised from one textural site to another. The contrasting behaviour of element mobility between peak eclogite-facies metamorphism and retrogression can be explained in terms of rate processes and fluid/deformation activity. During eclogitisation,deformation allows the formation of textural and mineralogical microdomains. Contemporaneously the fluid circulation favours the development of heterogeneities by transport of elements in excess towards areas protected from deformation where phengite-bearing quartz lenses crystallise. This, together with the very fast recrystallisation of eclogite-facies minerals, results in a heterogeneous redistribution of elements at sample scale and in a chemical disequilibrium between different microdomains of hundred-micron size. During retrogression in the amphibolite facies, coeval with long-term ductile deformation, the fluids act as an efficient vector to redistribute the elements from a given textural site to another and, thus, allow to reach, at hand-sample scale, a new chemical equilibrium between the different crystallising phases.
ISSN: 0024-4937

hal-00404602
https://hal.archives-ouvertes.fr/hal-00404602
DOI : 10.1016/j.lithos.2006.10.001

Éditeur(s) : HAL CCSD Mineralogical Society of America
Résumé : International audience
Accessory phases are important hosts of trace elements; allanite may contain >90% of the REE in a bulk rock. The mobility and redistribution of several trace elements, notably HREE, Th, U, and Y is thus controlled by reactions involving allanite and other REE phases, as well as several rock-forming minerals. As these elements are commonly concentrated in mature elastic sediments, a suite of impure quartzite was studied. Two eclogite facies samples from the Monometamorphic Cover Complex of the Sesia Zone (Western Italian Alps) are presented in some detail, as they reveal a remarkably rich spectrum of reaction relationships that involve REE phases.;Two allanite-forming reactions were inferred from textures and phase compositions;(1) monazite + Ca-silicate(?) + fluid -> allanite + apatite + thorite;;(2) monazite + thorite + Ca-silicate(?) + fluid -> Th-rich allanite + auerlite +/- apatite.;Petrographic observations and thermodynamic models suggest that allanite entered the HP assemblage at similar to 530 degrees C and 17-18 kbar during prograde metamorphism.;In one sample, allanite is rimmed by epidote rich in Y and HREE that grew at the expense of xenotime. Two net transfer reactions were derived;(3) xenotime + allanite + fluid -> Y-rich epidote + apatite + thorite;;(4) xenotime + allanite + fluid -> Y-rich epidote + aeschynite + thorite + (phosphate?).;Textural relationships and trace element analyses of coexisting allanite/monazite and xenotime/Y-rich epidote reveal systematic partitioning of the REE. Partition coefficients for the HREE are compatible with equilibrium fractionation, whereas those for the LREE show patterns that seem to be inherited from the precursor phases, in this case zircon with variable LREE composition.
ISSN: 0003-004X

hal-00745674
https://hal.archives-ouvertes.fr/hal-00745674
DOI : 10.2138/am.2012.3832

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : Geochemistry, Geophysics, Geosystems, v. 9, p. Q12001, 2008. http://dx.doi.org/10.1029/2008GC002098
International audience
Rb/Sr, 40Ar/39Ar, and Sm/Nd isotopic data are reported in Caledonian eclogites from the Lindås Nappe, Bergen Arcs, Norway, in order to investigate processes controlling isotopic equilibrium at mineral scale in polymetamorphic rocks. The Bergen Arcs exposes Sveconorwegian ∼950–930 Ma granulites, partially overprinted by Caledonian eclogite-facies metamorphism at ∼425 Ma and amphibolite-facies metamorphism at ∼410 Ma. Geochemical and Rb-Sr data from more than 10 phengite fractions separated from one sample reflect the composition of the microdomain (a few hundred microns in size) in which phengite crystallized. Phengite crystallized after garnet or plagioclase by dissolution-precipitation processes yield apparent age between 700 and 600 Ma. At the time of their crystallization, these phengites inherited the isotopic composition of their precursor minerals, at a microdomain scale. Phengite from quartz veins, which crystallized from elements mobilized by the circulating fluid, yield an age closer to the eclogite-facies metamorphic age. The closed system evolution of the eclogitizing fluid, the segregation of textural and reactional microdomains, the high Sr content of the studied phengite, and the short duration of the recrystallization processes (<1Ma) are interpreted as the main factors responsible for the lack of a Rb/Sr isotopic equilibrium, at the scale of hand samples. Such equilibrium is nevertheless reached in quartz veins where the crystallization of minerals implies that the fluid circulation acted as a factor of isotopic homogenization. The in situ single mineral Ar/Ar datings revealed that both eclogite- and amphibolite-facies minerals are characterized by an excess of radiogenic argon (ages between 425 and 520 Ma). The excess of argon has been inherited from the previous granulite and has been only partially evacuated from the system by the circulating fluid. The apparent Sm/Nd ages from garnets inherited from the granulite-facies metamorphism (c. 930 Ma) are in agreement with previous estimates. This result confirms that the diffusion of Nd in garnet does not occur at temperatures lower than 700°C. This study highlights the complexity of radiochronometer behavior during HP metamorphism and demonstrates that coupling different radiochronometers, such as Ar/Ar and Rb/Sr, does not always guarantee the validity of the geochronological results.
ISSN: 1525-2027

hal-00407885
https://hal.archives-ouvertes.fr/hal-00407885
DOI : 10.1029/2008GC002098

É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
Before attaining the mantle wedge, where they trigger partial melting, volatiles released from dehydration reactions in the slab have to migrate across a relatively cold (<750 degrees C), peridotite-layer above the incoming slab. In order to unravel the mechanisms allowing for this initial stage of fluid transport, we performed a detailed field and microstructural study of metamorphic prograde peridotites in the Cerro del Almirez ultramafic massif (Betic Cordillera, Spain), where evidences of one of the most important dehydration reactions in subduction zones, the high-pressure antigorite breakdown (P = 1.6-1.9 GPa and T approximate to 680 degrees C), can be mapped in the field. This reaction led to arborescent growth of centimeter-size olivine and orthopyroxene, producing a chlorite-harzburgite with a spinifex-like texture. Microstructural observations and crystal preferred (CPO) mapping show no evidences of solid-state deformation during the prograde growth of olivine and orthopyroxene at the expenses of antigorite. However, a few tens to a hundred meters away from the reaction front, the metamorphic texture is partially obliterated by grain-size reduction in roughly planar conjugate zones, a few mm to meters wide. Grain size reduction zones (GSRZ) are characterized by (1) sharp contacts with undeformed spinifex-like texture domains, (2) important reduction of the olivine grain size (60-250 mu m), (3) olivine color change from brownish to colorless, (4) decrease in the modal amount of orthopyroxene, and (5) at the mm- to cm-scale, irregular shapes and abrupt terminations. Field and microstructural observations exclude that relative displacement took place across these GSRZ. Changes in modal composition imply reactions with fluids undersaturated in silica. Analysis of olivine crystal-preferred orientations (CPO) in GSRZ shows patterns similar, but more dispersed, than those in neighboring spinifex-like domains. It also reveals mm- to cm-scale discrete domains with rather homogeneous crystallographic orientations suggesting inheritance from the preexisting spinifex-like olivines in the host peridotite. Misorientation angles between neighboring grains in the GSRZ show peaks at similar to 5-10 degrees and similar to 20 degrees, but rotations are not crystallographically controlled. Based on these observations, we rule out the formation of the GSRZ by dynamic recrystallization during dislocation creep and propose that they record brittle deformation (microcraking) of the spinifex-like chlorite-harzburgite, probably induced by hydrofracturing at high pressure and relative low temperature conditions (680-710 degrees C). High-pressure hydrofracturing can, thus, be invoked as an efficient mechanism for fluid flow across the cold top-slab mantle layer, hence allowing the slab-derived fluids to ingress in the wedge.
ISSN: 0012-821X

hal-00534035
https://hal.archives-ouvertes.fr/hal-00534035
DOI : 10.1016/j.epsl.2010.06.029

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
In the Tahalgha district (southwestern Hoggar, Algeria), the Cenozoic volcanism has sampled subcontinental mantle beneath two crustal terranes that collided during the Pan-African orogeny: the “Polycyclic Central Hoggar” to the east and the “Western Hoggar” to the west. Two major lithospheric shear zones separate these terranes: the “4°35” and the “4°50” faults. Mantle xenoliths were collected between the two faults and across the 4°35 fault. In addition to a range in equilibrium temperatures and chemical compositions reported elsewhere, the samples show variations in their microstructures and crystallographic preferred orientations. Equilibrium temperatures and geochemical characteristics allow dividing them into low — (LT; 700–900 °C), intermediate — (IT; 900–1000 °C), and high-temperature (HT; 1000–1100 °C) xenoliths. The LT and IT peridotites occur on both sides of the 4°35 fault; they are usually coarse-grained. HT xenoliths are present only east of the 4°35 fault, in the narrow domain stuck between the two faults; they are fine-grained and extensively affected by annealing and melt–rock reactions. Microstructures and crystallographic textures indicate that deformation in the LT- and IT-xenoliths occurred through dislocation creep under relatively high-temperature, low-pressure conditions, followed by post-kinematic cooling. The fine-grained HT-xenoliths were deformed under relatively high-stress conditions before being annealed.Combining microstructural and CPO data with petrological and geochemical informations suggests that: (1) the LT xenoliths are remnants of the Neoproterozoic lithospheric mantle that preserved microstructural and chemical characteristics inherited from the Pan-African orogeny, and (2) the HT xenoliths record localized Cenozoic deformation associated with melt channeling through feed-back processes that culminated in the formation of high-permeability porous-flow conduits. Limited grain-growth in HT xenoliths suggests that advective heating of melt conduits was transient and rapidly followed by thermal relaxation due to conductive heat loss into wall-rock peridotites represented by the IT xenoliths, then by exhumation due to volcanic activity.
ISSN: 0040-1951

hal-01171761
https://hal.archives-ouvertes.fr/hal-01171761
DOI : 10.1016/j.tecto.2014.11.012

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Highly siderophile elements (Platinum-group elements, Au and Re) are currently assumed to reside inside base metal sulfides (BMS) in the convecting upper mantle. However, fertile lherzolites sampled by Pyrenean orogenic peridotite massifs are unexpectedly rich in 0.5-3 mu m large micronuggets of platinum-group minerals (PGM). Among those, sulfides from the laurite-erlichmanite series (Ru, Os(Ir)S(As)2), Pt-Ir-Os alloys and Pt-Pd-Te-Bi phases (moncheite-merenskyite) are predominant. Not only the BMS phases but also the PGM micronuggets must be taken into account in calculation of the PGE budget of orogenic fertile lherzolites. Laurite is a good candidate for equilibrating the whole-rock budget of Os, it and Ru while accounting for supra-chondritic Ru/Ir-N. Textural relationships between PGMs and BMS highlight heterogeneous mixing between refractory PGMs (laurite/Pt-Ir-Os alloys) inherited from ancient refractory lithospheric mantle and late-magmatic metasomatic sulfides precipitated from tholeiitic melts. "Low-temperature" PGMs, especially Pt-Pd bismuthotellurides should be added to the list of mineral indicators of lithosphere refertilization process. Now disseminated within fertile lherzolites, "lithospheric" PGMs likely account for local preservation of ancient Os model ages (up to 2 Ga) detected in BMS by in-situ isotopic analyses. These PGMs also question the reliability of orogenic lherzolites for estimating the PGE signature of the Primitive Silicate Earth.
ISSN: 0012-821X

hal-00486109
https://hal.archives-ouvertes.fr/hal-00486109
DOI : 10.1016/j.epsl.2009.11.017

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
ISSN: 0012-821X

hal-00380865
https://hal.archives-ouvertes.fr/hal-00380865
DOI : 10.1016/j.epsl.2008.04.020

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
We report on microstructural data obtained by optical microscopy and transmission electron microscopy concerning the crystallographic relationships of serpentine minerals with their host olivine in two contrasting situations. In the first case, mesh-textured lizardite (liz) is developed in a standard 60% serpentinized oceanic harzburgite from the Oman ophiolite where olivine converts to columnar lizardite. The joined columns are perpendicular to the basal plane (001)(liz), corresponding to the pseudofibres observed optically. The plane (001)(liz) is locally parallel to the narrow boundary ol-liz; thus column orientations register the interface of serpentinization. The ol-liz relationships are not strictly topotactic, but reflect preferred cracking orientations in olivine, parallel to (010)(ol). In the second case, antigorite (atg) develops in a rare sample of antigorite schist in a kimberlite from Moses Rock (Colorado Plateau), representative of a suprasubduction-zone mantle wedge. High-resolution transmission electron microscopy (HRTEM) images along [010](atg) show domains of very regular modulation with a 43 center dot 5 A wavelength (m = 17, where m is the number of silicate tetrahedra along the wave), with few defects, indicative of HP-HT antigorite, and also heavily kinked regions as fingerprints of strong tectonic shear. TEM imaging and electron diffraction patterns reveal two topotactic relationships between antigorite and olivine: [100](atg)//[010](ol) and < 100 >(atg)//< 100 >(ol); the planes in contact are (001)(atg)//(100)(ol) and (001)(atg)//(010)(ol), respectively. The [010](atg)//[001](ol) and antigorite lamellae are parallel to the forsterite b-axis. In both cases, the topography of olivine-serpentine interfaces is controlled by open fluid pathways along microcracks oriented according to the anisotropy of the olivine aggregate. In the cases studied, the serpentine aggregate exhibits a preferred orientation inherited from that of the peridotite. These results have some relevance to the seismic anisotropy of serpentinized mantle. Anisotropy of propagation of seismic waves as a result of the olivine fabric is maintained and reinforced with the development of lizardite. Conversely, the development of antigorite produces a trench-parallel fast S-wave polarization and an anisotropy that is lowered at low degrees of serpentinization and then increased with increasing serpentinization.
ISSN: 0022-3530

hal-00456061
https://hal.archives-ouvertes.fr/hal-00456061
DOI : 10.1093/petrology/egp049

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Microstructures, U-Pb SHRIMP zircon ages and temperature estimates were used to constrain the fabric evolution and timing of mylonitisation in the Patos shear zone (NE Brazil). Melt-bearing HT-mylonitic orthogneiss displays solid-state fabrics with coarse quartz ribbons and sutured grain boundaries. K-feldspar is often fractured and shows peripheral myrmekite. The textures typically become magmatic at the contacts with the Espinho Branco anatexite. Zircon crystals from two samples of the leucosome show Paleoproterozoic (ca. 2.2 Ga) inherited cores that are enveloped by recrystallised Neoproterozoic rims. These zircon grains define a discordia with lower intercept ages of 558 and 562 Ma but with large analytical errors. A leucogranite with transitional contacts with host diatexites provides a well-constrained zircon mean age of 566 ± 6 Ma (n = 15, 95% conf.), which is considered the best estimate for the peak metamorphic conditions. The transition from HT-mylonites to fine-grained mylonites and ultramylonites is marked by a progressive grain size reduction assisted by recrystallisation. Quartz crystallographic fabrics display [0 0 0 1] concentrations between Z and Y. K-feldspar and plagioclase fabrics record the activity of the (0 1 0)[0 0 1] and (0 1 0)[1 0 0] slip systems and, in the lower-temperature tectonites, the (1 0 0)[0 1 0] slip system. These results suggest a diachronic evolution in which an initial HT northern domain was deformed by intracrystalline slip and melt-assisted granular flow, followed by medium- to low-temperature solid-state creep in the southern shear zone domain. The mylonites therefore record the rheological heterogeneities responsible for the formation of a late continental shear zone that nucleated under partially molten conditions and was later reactivated through melt-absent strain localization following exhumation.
ISSN: 0301-9268

hal-01053497
https://hal.archives-ouvertes.fr/hal-01053497
DOI : 10.1016/j.precamres.2013.12.020

Éditeur(s) : HAL CCSD Springer Verlag (Germany) Springer Verlag
Résumé : Dehydration reactions result in minerals’ replacement and a transient fluid-filled porosity. These reactions involve interface-coupled dissolution–precipitation and might therefore lead to fixed crystallographic orientation relations between reactant (protolith) and product phases (i.e. topotaxy). We investigate these two phenomena in the dehydration of a foliated antigorite (atg) serpentinite by comparing the crystallographic preferred orientation (CPO) developed by olivine (ol), orthopyroxene (opx) and chlorite (chl) during high-pressure antigorite and chlorite dehydration in piston-cylinder experiments and in natural samples recording the dehydration of antigorite (Cerro del Almirez, Betic Cordillera, Spain). Experiments were performed under undrained conditions resulting in fluid-filled porosity and in strong CPO of the prograde minerals, controlled by the pre-existing antigorite CPO in the reactant foliated serpentinite. The orientation of a ol,opx and c∗chl is parallel to c∗atg from the protolith. The Cerro del Almirez samples show similar, locally well-developed topotactic relations between orthopyroxene, chlorite and antigorite, but the product CPOs are weaker and more complex at the thin section scale. In contrast to the experiments, olivine from natural samples shows a weak correlation between b ol and the former c∗atg. We relate the strengthening of local topotactic relations and the weakening of the inherited CPO at a larger scale in natural samples to compaction and associated fluid migration. Microstructural features that might be related to compaction in the natural samples include: (1) smooth bending of the former foliation, (2) gradual crystallographic misorientation (up to 16°) of prismatic orthopyroxene due to buckling by dislocation creep, (3) inversion of enstatite to low clinoenstatite (P21/c) along lamellae and (4) brittle fracturing of prismatic orthopyroxene enclosed by plastically deformed chlorite. The coexistence of orthopyroxene buckling and clinoenstatite lamellae enables estimating the local strain rates and shear stresses generated during compaction. An lower bound for the strain rates in the order of 10−12 to 10−13 s−1 and shear stresses of 60–70 MPa are estimated based on creep data. Lower shear stresses (20–40 MPa) are retrieved using a theoretical approach. These data point to slow compaction (and fluid extraction) in nature if the system is not perturbed by external forces, with rates only marginally higher than the viscoplastic deformation of the solid matrix.
ISSN: 0010-7999

hal-01174158
https://hal.archives-ouvertes.fr/hal-01174158
DOI : 10.1007/s00410-015-1129-4

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Synkinematic phyllosilicates in fault zones can be used to deduce the deformation mechanisms and the conditions of fault activity, as their chemical composition, crystal structure and texture can record the different stages of deformation and fluid-rock interactions. The Pic de Port Vieux, a second-order thrust related to the major Gavarnie thrust in the southern central part of the Pyrenees Axial Zone, juxtaposes Triassic pelites of the hanging wall and Cretaceous limestones of the footwall. In order to investigate the mineralogical and geochemical changes and constrain the deformation conditions of thrusting, characterization of phyllosilicates was performed along a transect in the fault hanging wall pelites. The Triassic pelites are mainly composed of quartz, calcite, phyllosilicates (chlorite and K-white mica) and ± hematite. The core fault zone thickness is estimated to be about one meter of intensively foliated green pelites, whereas the damage zone is composed of several meters of red pelites. XRD data demonstrated that the difference in color is related to hematite which is only present in the damage zone. Phyllosilicates of the damage zone are mainly inherited/diagenetic K-white mica and chlorite. In the core zone, newly formed chlorite is abundant and preferentially located in veins. It is enriched in Fe compared to the chlorite of the red pelites. The well-defined foliation in the core zone is overlaid by preferentially oriented muscovite grains which have homogeneous compositions, with less Na and relatively more Fe than mica from the red pelites. Newly formed chlorite and muscovite in the core zone are synkinematic to the fault activity. They are both related to deformation processes and fluid rock interactions. Kübler index measurements and chlorite thermometry show that synkinematic phyllosilicates have registered a temperature of 270°C ± 23°C for the damage zone and 285°C ± 28°C for the core zone which are corresponding to lower-anchizone/epizone grade conditions. Numerical modelling with the geochemical modelling program PhreeqC was performed to determine the favorable conditions for the mineralogical change between red and green pelites. It suggests that the main critical parameter favoring hematite dissolution and chlorite precipitation in the core zone, is the redox conditions. According to the model, the chemical changes in the core zone occurred due to interactions with highly reductive fluids. In addition, the hematite dissolution maybe the source of iron for newly formed phyllosilicates in the core zone that are more iron rich compared to those from the red pelites.
ISSN: 0169-1317

hal-01536006
https://hal.archives-ouvertes.fr/hal-01536006
DOI : 10.1016/j.clay.2017.05.008