A stable and 40Ar/39Ar isotope study of a major thrust in the Helvetic nappes (Swiss Alps) - Evidence for fluid-flow and constraints on nappe kinematics

Stable isotope and Ar-40/Ar-39 measurements,were made on samples associated with a major tectonic discontinuity in the Helvetic Alps, the basal thrust of the Diablerets nappe (external zone of the Alpine Belt) in order to determine both the importance of fluids in this thrust zone and the timing of thrusting. A systematic decrease in the delta(18)O values (up to 6 parts per thousand) of calcite, quartz, and white mica exists within a 10- to 70-m-wide zone over a distance of 37 km along the thrust, and they become more pronounced toward the root of the nappe. A similar decrease in the delta(13)C values of calcite is observed only in the deepest sections (up to 3 parts per thousand). The delta D-SMOW (SMOW = standard mean ocean water) values of white mica are -54 parts per thousand +/- 8 parts per thousand (n = 22) and are independent of the distance from the thrust. These variations are interpreted to reflect syntectonic solution reprecipitation during fluid passage along the thrust. The calculated delta(18)O and delta D values (versus SMOW) for the fluid in equilibrium with the analyzed minerals is 12 parts per thousand to 16 parts per thousand and -30 parts per thousand to +5 parts per thousand, respectively, for assumed temperatures of 250 to 450 degrees C. The isotopic and structural data are consistent with fluids derived from the deep-seated roots of the Helvetic nappes where large volumes of Mesozoic sediments were metamorphosed to the amphibolite facies, It is suggested that connate and metamorphic waters, overpressured by rapid tectonic burial in a subductive system escaped by upward infiltration along moderately dipping pathways until they reached the main shear zone at the base of the moving pile, where they were channeled toward the surface, This model also explains the mechanism by which large amounts of fluid were removed from the Mesozoic sediments during Alpine metamorphism. White mica Ar-49/Ar-39 ages vary from 27 Ma far from the Diablerets thrust to 15 Ma along the thrust. An older component is observed in micas far from the thrust, interpreted as a detrital signature, and indicates that regional metamorphic temperatures were less than about 350 degrees C. The;plateau and near plateau ages nearest the thrust are consistent with either neocrystallization of white mica or argon loss by recrystallization during thrusting, which may have been enhanced in the zones of highest fluid flow. The 15 Ma Ar-40/Ar-39 age plateau measured on white mica sampled exactly on the thrust surface dates the end of both fluid flow and tectonic transport.

Recrystallisation rims in zircon (Valle d'Arbedo, Switzerland): an integrated cathodoluminescence, LA-ICP-MS, SHRIMP, and TEM study

Recrystallization rims are a common feature of zircon crystals that underwent metamorphism. We present a microstructural and microchemical study of partially recrystallized zircon grains collected in polymetamorphic migmatites (Valle d'Arbedo, Ticino, Switzerland). The rims are bright in cathodo-luminescence (CL), with sharp and convex contacts characterized by inward-penetrating embayments transgressing igneous zircon cores. Laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) data and transmission electron microscopy (TEM) imaging indicate that the rims are chemically and microstructurally different from the cores. The rims are strongly depleted in REE, with concentrations up to two orders of magnitude lower than in the cores, indicating a significant loss of REE during zircon recrystallization. Enrichment in non-formula elements, such as Ca, has not been observed in the rims. The microstructure of zircon cores shows a dappled intensity at and below the 100 nm scale, possibly due to radiation damage. Other defects such as pores and dislocations are absent in the core except at healed cracks. Zircon rims are mostly dapple-free, but contain nanoscale pores and strain centers, interpreted as fluid inclusions and chemical residues, respectively. Sensitive high-resolution ion microprobe (SHRIMP) U-Pb ages show that the recrystallization of the rims took place >200 Ma ago when the parent igneous zircon was not metamict. The chemical composition and the low-Ti content of the rims indicate that they form at sub-solidus temperatures (550-650 degrees C). Recrystallization rims in Valle d'Arbedo zircon are interpreted as the result of the migration of chemical reaction fronts in which fluid triggered in situ and contemporaneous interface-coupled dissolution-reprecipitation mechanisms. This study indicates that strong lattice strain resulting from the incorporation of a large amount of impurities and structural defects is not a necessary condition for zircon to recrystallize. Our observations suggest that the early formation of recrystallization rims played a major role in preserving zircon from the more recent Alpine metamorphic overprint.

Early Carboniferous age of the Versoyen ophiolites and consequences: non-existence of a "Valais ocean'' (Lower Penninic, western Alps)

Ophiolites occur at several places in the Lower Penninic of the W and Central Alps. They are generally ascribed to oceanic crust of a so-called ``Valais ocean'' of Cretaceous age which plays a fundamental role in many models of Alpine paleogeography and geodynamics. The type locality and only observational base for the definition of a ``Valais ocean'' in the W Alps is the Versoyen ophiolitic complex, on the French-Italian boundary W of the Petit St-Bernard col. The idea of a "Valais ocean'' is based on two propositions that are since 40 years the basis for most reconstructions of the Lower Penninic: (1) The Versoyen forms the (overturned) stratigraphic base of the Cretaceous-Tertiary Valais-Tarentaise series; and (2) it has a Cretaceous age. We present new field and isotopic data that severely challenge both propositions. (1) The base of the Versoyen ophiolite is a thrust. It overlies a wildflysch with blocks of Versoyen rocks, named the Mechandeur Formation. This ``supra-Tarentaise'' wildflysch has been confused with an (overturned) stratigraphic transition from the Versoyen to the Valais-Tarentaise series. Thus the contact Versoyen/Tarentaise is not stratigraphic but tectonic, and the Versoyen ophiolite has no link with the Valais basin. This thrust corresponds to an inverse metamorphic discontinuity and to an abrupt change in tectonic style. (2) The contact of the Versoyen complex with the overlying Triassic-Jurassic Petit St-Bernard (PSB) series is stratigraphic (and not tectonic as admitted by all authors since 50 years). Several types of sedimentary structures polarize it and show that the PSB series is younger than the Versoyen. Consequently the Versoyen ophiolitic complex is Paleozoic and forms the basement of the PSB Mesozoic sediments. They both belong to a single tectonic unit, named the Versoyen-Petit St-Bernard nappe. (3) Ion microprobe U-Pb isotopic data on zircons from the main gabbroic intrusion in the Versoyen complex give a crystallization age of 337.0 +/- 4.1 Ma (Visean, Early Carboniferous). These zircons show typical oscillatory zoning and no overgrowth or corrosion. and are interpreted to date the Versoyen magmatism. These U-Pb data are in excellent agreement with our field observations and confirm the Paleozoic age of the Versoyen ophiolite. The existence of a ``Valais ocean'' of Cretaceous age in the W Alps becomes very improbable. The eclogite facies metamorphism of the Versoyen-Petit St-Bernard nappe results from an Alpine intra-continental subduction, guided by a Paleozoic oceanic suture. This is an example of the lone term influence of inherited deep-seated structures on a Much younger orogeny. This might well be a major cause of of the inherent complexity of the Alps.

Devolatilization-induced pressure build-up: Implications for reaction front movement and breccia pipe formation

Generation of fluids during metamorphism can significantly influence the fluid overpressure, and thus the fluid flow in metamorphic terrains. There is currently a large focus on developing numerical reactive transport models, and with it follows the need for analytical solutions to ensure correct numerical implementation. In this study, we derive both analytical and numerical solutions to reaction-induced fluid overpressure, coupled to temperature and fluid flow out of the reacting front. All equations are derived from basic principles of conservation of mass, energy and momentum. We focus on contact metamorphism, where devolatilization reactions are particularly important owing to high thermal fluxes allowing large volumes of fluids to be rapidly generated. The analytical solutions reveal three key factors involved in the pressure build-up: (i) The efficiency of the devolatilizing reaction front (pressure build-up) relative to fluid flow (pressure relaxation), (ii) the reaction temperature relative to the available heat in the system and (iii) the feedback of overpressure on the reaction temperature as a function of the Clapeyron slope. Finally, we apply the model to two geological case scenarios. In the first case, we investigate the influence of fluid overpressure on the movement of the reaction front and show that it can slow down significantly and may even be terminated owing to increased effective reaction temperature. In the second case, the model is applied to constrain the conditions for fracturing and inferred breccia pipe formation in organic-rich shales owing to methane generation in the contact aureole.

The Tertiary structural and thermal evolution of the central Alps - Compressional and extensional structures in an orogenic belt

The Western Alpine Are has been created during the Cretaceous and the Tertiary orogenies. The interference patterns of the Tertiary structures suggest their formation during continental collision of the European and the Adriatic Plates, with an accompanying anticlockwise rotation of the Adriatic indenter. Extensional structures are mainly related to ductile deformation by simple shear. These structures developed at a deep tectonic level, in granitic crustal rocks, at depths in excess of 10 km. In the early Palaeogene period of the Tertiary Orogeny, the main Tertiary nappe emplacement resulted from a NW-thrusting of the Austroalpine, Penninic and Helvetic nappes. Heating of the deep zone of the Upper Cretaceous and Tertiary nappe stack by geothermal heat flow is responsible for the Tertiary regional metamorphism, reaching amphibolite-facies conditions in the Lepontine Gneiss Dome (geothermal gradient 25 degrees C/ km). The Tertiary thrusting occurred mainly during prograde metamorphic conditions with creation of a penetrative NW-SE-oriented stretching lineation, X(1) (finite extension), parallel to the direction of simple shear. Earliest cooling after the culmination of the Tertiary metamorphism, some 38 Ma ago, is recorded by the cooling curves of the Monte Rosa and Mischabel nappes to the west and the Suretta Nappe to the east of the Lepontine Gneiss Dome. The onset of dextral transpression, with a strong extension parallel to the mountain belt, and the oldest S-vergent `'backfolding'' took place some 35 to 30 Ma ago during retrograde amphibolite-facies conditions and before the intrusion of the Oligocene dikes north of the Periadriatic Line. The main updoming of the Lepontine Gneiss Dome started some 32-30 Ma ago with the intrusion of the Bergell tonalites and granodiorites, concomitant with S-vergent backfolding and backthrusting and dextral strike-slip movements along the Tonale and Canavese Lines (Argand's Insubric phase). Subsequently, the center of main updoming migrated slowly to the west, reaching the Simplon region some 20 Ma ago. This was contemporaneous with the westward migration of the Adriatic indenter. Between 20 Ma and the present, the Western Aar Massif-Toce culmination was the center of strong uplift. The youngest S-vergent backfolds, the Glishorn anticline and the Berisal syncline fold the 12 Ma Rb/Sr biotite isochron and are cut by the 11 Ma old Rhone-Simplon Line. The discrete Rhone-Simplon Line represents a late retrograde manifestation in the preexisting ductile Simplon Shear Zone. This fault zone is still active today. The Oligocene-Neogene dextral transpression and extension in the Simplon area were concurrent with thrusting to the northwest of the Helvetic nappes, the Prealpes (35-15 Ma) and with the Jura thin-skinned thrust (11-3 Ma). It was also contemporaneous with thrusting to the south of the Bergamasc (> 35-5 Ma) and Milan thrusts (16-5 Ma).

Modelling changes in stable isotope compositions of minerals during net transfer reactions in a contact aureole: Wollastonite growth at the northern Hunter Mountain Batholith (Death Valley National Park, USA)

One of the world's largest wollastonite deposits was formed at the contact of the northern Hunter Mountain Batholith (California, USA) in Paleozoic sediments. Wollastonite occurs as zones of variable thickness surrounding layers or nodules of quartzite in limestones. A minimum formation temperature of 650 degrees C is estimated from isolated periclase-bearing lenses in that area. Contact metamorphism of siliceous carbonates has produced mineral assemblages that are consistent with heterogeneous, and partly limited infiltration of water-rich fluids, compatible with O-18/O-16 and C-13/C-12 isotopic patterns recorded in carbonates. Oxygen isotope compositions of wollastonites in the study area may also not require infiltration of large quantities of externally-derived fluids that were out of equilibrium with the rocks. 8180 values of wollastonite are high (14.8 parts per thousand to 25.0 parts per thousand; median: 19.7 parts per thousand) and close to those of the host limestone (19.7 parts per thousand to 28 parts per thousand; median: 24.9 parts per thousand) and quartz (18.0 parts per thousand. to 29.1 parts per thousand; median: 22.6 parts per thousand). Isotopic disequilibrium exists at quartz/wollastonite and wollastonite/calcite boundaries. Therefore, classical batch/Rayleigh fractionation models based on reactant and product equilibrium are not applicable to the wollastonite rims. An approach that relies on local instantaneous mass balance for the reactants, based on the wollastonite-forming reaction is suggested as an alternative way to model wollastonite reaction rims. This model reproduces many of the measured delta O-18 values of wollastonite reaction rims of the current study to within +/- 1 parts per thousand, even though the wollastonite compositions vary by almost 10 parts per thousand. (C) 2011 Elsevier B.V. All rights reserved.

Structural and metamorphic evolution of the northern Himachal Himalaya, NW India - (Spiti-eastern Lahul-Parvati valley traverse)

The Himalayan orogen is the result of the collision between the Indian and Asian continents that began 55-50 Ma ago, causing intracontinental thrusting and nappe formation. Detailed mapping as well as structural and microfabric analyses on a traverse from the Tethyan Himalaya southwestward through the High Himalayan Crystalline and the Main Central Thrust zone (MCT zone) to the Lesser Himalayan Sequence in the Spiti-eastern Lahul-Parvati valley area reveal eight main phases of deformation, a series of late stage phases and five stages of metamorphic crystallization. This sequence of events is integrated into a reconstruction of the tectonometamorphic evolution of the Himalayan orogen in northern Himachal Pradesh. The oldest phase D-1 is preserved as relies in the High Himalayan Crystalline. Its deformational conditions are poorly known, but the metamorphic evolution is well documented by a prograde metamorphism reaching peak conditions within the upper amphibolite facies. This indicates that D-1 was an important tectonometamorphic event including considerable crustal thickening. The structural, metamorphic and sedimentary record suggest that D-1 most probably represents an early stage of continental collision. The first event clearly attributed to the collision between India and Asia is documented by two converging nappe systems, the NE-verging Shikar Beh Nappe and the SW-verging north Himalayan nappes. The D-2 Shikar Beh Nappe is characterized by isoclinal folding and top-to-the NE shearing, representing the main deformation in the High Himalayan Crystalline. D-2 also caused the main metamorphism in the High Himalayan Crystalline that was of a Barrovian-type, reaching upper amphibolite facies peak conditions. The Shikar Beh Nappe is interpreted to have formed within the Indian crust SW of the subduction zone. Simultaneously with NE-directed nappe formation, incipient subduction of India below Asia caused stacking of the SW-verging north Himalayan Nappes, that were thrust from the northern edge of the subducted continent toward the front of the Shikar Beh Nappe. As a result, the SW-verging folds of the D-3 Main Fold Zone formed in the Tethyan Himalaya below the front of the north Himalayan nappes. D-3 represents the main deformation in the Tethyan Himalaya, associated with a greenschist facies metamorphism. Folding within the Main Fold Zone subsequently propagated toward SW into the High Himalayan Crystalline, where it overprinted the preexisting D-2 structures. After subduction at the base of the north Himalayan nappes, the subduction zone stepped to the base of the High Himalayan Crystalline, where D-3 folds were crosscut by SW-directed D-4 thrusting. During D-4, the Crystalline Nappe, comprising the Main Fold Zone and relies of the Shikar Beh Nappe was thrust toward SW over the Lesser Himalayan Sequence along the 4 to 5 kms thick Main Central Thrust zone. Thrusting was related to a retrograde greenschist facies overprint at the base of the Crystalline Nappe and to pro-grade greenschist facies conditions in the Lesser Himalayan Sequence. Simultaneously with thrusting at the base of the Crystalline Nappe, higher crustal levels were affected by NE-directed D-5 normal extensional shearing and by dextral strike-slip motion, indicating that the high-grade metamorphic Crystalline Nappe was extruded between the low-grade metamorphic Lesser Himalayan Sequence at the base and the north Himalayan nappes at the top. The upper boundary of the Crystalline Nappe is not clearly delimited and passes gradually into the low-grade rocks at the front of the north Himalayan nappes. Extrusion of the Crystalline Nappe was followed by the phase D-6, characterized by large-scale, upright to steeply inclined, NE-verging folds and by another series of normal and extensional structures D-7+D-8 that may be related to ongoing extrusion of the Crystalline Nappe. The late stage evolution is represented by the phases D-A and D-B that indicate shortening parallel to the axis of the mountain chain and by D-C that is interpreted to account for the formation of large-scale domes with NNW-SSE-trending axes, an example of which is exposed in the Larji-Kullu-Rampur tectonic window.

Implications of the serpentine phase transition on the behaviour of beryllium and lithium-boron of subducted ultramafic rocks

The Totalp-Platta-Malenco ophiolites in the Eastern Central Alps offer a unique opportunity to study the behaviour of Li, Be and B in ultramafic rocks in response to serpentinization and to progressive Alpine metamorphism. These units represent the remnants of a former ocean-continent transition that was intensely serpentinized during exposure on the Jurassic seafloor of the Ligurian Tethys. From north to the south, three isograd reactions (lizardite double right arrow antigorite + brucite; lizardite + talc double right arrow antigorite; lizardite + tremolite double right arrow antigorite + diopside) have been used to quantify the evolution of the light element content of metamorphic minerals. We determined the Li, Be and B concentrations in major silicate minerals from the ultramafic bodies of Totalp, Platta and Malenco by secondary ion mass spectrometry. Mantle minerals have Be concentrations (e.g. <0.001-0.009 mu g/g in olivine) similar to the metamorphic minerals that replace them (e.g. <0.001-0.016 mu g/g in serpentine). The mantle signature of Be is thus neither erased during seafloor alteration nor by progressive metamorphism from prehnite-pumpellyite to epidote-amphibolite facies. In contrast, the Li and B inventories of metamorphic minerals are related to the lizardite-to-antigorite transition. Both elements display higher concentrations in the low-temperature serpentine polymorph lizardite (max. 156 mu/g Li, max. 318 mu g/g B) than in antigorite (max. 0.11 mu g/g Li, max. 12 mu g/g B). Calculated average B/Li ratios for lizardite (similar to 1395) and antigorite (similar to 115) indicate that Li fractionates from B during the lizardite-to-antigorite transition during prograde metamorphism in ultramafic rocks. In subduction zones, this signature is likely to be recorded in the B-rich nature of forearc fluids. Relative to oceanic mantle the Be content of mantle clinopyroxene is much higher, but similar to Be values from mantle xenoliths and subduction-related peridotite massifs. These data support previous hypothesis that the mantle rocks from the Eastern Central Alps have a subcontinental origin. We conclude that Be behaves conservatively during subduction metamorphism of ultramafic rocks, at least at low-temperature, and thus retains the fingerprint of ancient subduction-related igneous events in mantle peridotites. (C) 2010 Elsevier Ltd. All rights reserved.

The duration of prograde garnet crystallization in the UHP eclogites at Lago di Cignana, Italy

The distinct core-to-rim zonation of different REEs in garnet in metamorphic rocks, specifically Sm relative to Lu, suggests that Sm-Nd and Lu-Hf isochron ages will record different times along a prograde garnet growth history. Therefore, REE zonations in garnet must be measured in order to correctly interpret the isochron ages in terms of the garnet growth interval, which could span several m.y. New REE profiles, garnet crystal size distributions, and garnet growth modeling, combined with previously published Sm-Nd and Lu-Hf geochronology on a UHP eclogite of the Zermatt-Saas Fee (ZSF) ophiolite, Lago di Cignana (Italy), demonstrate that prograde garnet growth of this sample occurred over a similar to 30 to 40 m.y. interval. Relative to peak metamorphism at 38 to 40 Ma, garnet growth is estimated to have begun at similar to 11 to 14 kbar pressure at similar to 70 to 80 Ma. Although such a protracted garnet growth interval is surprising, this is supported by plate tectonic reconstructions which suggest that subduction of the Liguro-Piemont ocean occurred through slow and oblique convergence. These results demonstrate that REE zonations in garnet, coupled to crystal size distributions, provide a powerful means for understanding prograde metamorphic paths when combined with Sm-Nd and Lu-Hf geochronology. (C) 2009 Elsevier B.V. All rights reserved.

Kinetics of crystal growth and equilibrium domains in eclogite of the Sesia Zone, Western Alps

Kinetics of crystal growth and equilibrium domains in eclogite of the Sesia Zone, Western Alps Darbellay Bastien Institut de Minéralogie et Géochimie Résumé grand public Comme toute matière, la roche est sensible à son environnement et cherche à s'adapter pour acquérir un état stable (état d'équilibre). Les changements des conditions physiques (température et pression) vont ainsi impliquer des modifications dans la roche. Le métamorphisme est l'étude de ces changements. Les minéraux qui constituent la roche peuvent modifier, leur structure, leur chimie ou être remplacer par d'autres minéraux plus stables. Il est ainsi crucial de déterminer les processus responsables et limitant de la croissance minérale. Trois processus permettent la croissance ; (1) la dissolution des éléments du réactant, (2) le transport de ces éléments vers le site de croissance, (3) l'incorporation de ces éléments dans la nouvelle structure. Cette thèse se focalise sur les structures des minéraux de haute pression (forme, zonation chimique, structure interne) pour essayer de déterminer les facteurs importants à l'origine de leur état final. Les zones d'étude se situent dans la zone de Sésia. La première partie traite de la problématique liée à l'incorporation d'un élément dans une structure minérale. A l'image de la croissance humaine, les irrégularités minéralogiques permettent de mettre en lumière un dysfonctionnement de la croissance due à un excès ou à une carence d'un élément. Bien dosé, cet élément est cependant essentiel à la croissance. Les zoisites (épidotes) des métabasites de la région de Cima di Bonze montrent une zonation chimique en sablier. Dans cette zonation la teneur en fer excède la capacité maximum que peut contenir la structure orthorhombique de la zoisite. Des défauts de structure permettent l'accommodation de cet excès. La zoisite peut ainsi adapter sa structure pour permettre l'incorporation d'une relativement grande quantité de fer. Les études précédentes montraient, pour des conditions similaires, la formation de deux épidotes distinctes. La deuxième partie se penche sur la compétition entre le minéral qui fait sa croissance et les minéraux (réactants) qui l'entourent. Les métapélites de la région du Monte Mucrone contiennent des grenats atollaires. Des études détaillées de la texture et de la zonation chimique du grenat ainsi qu'une modélisation thermodynamique ont permis de mieux cerner les facteurs importants responsables de la forme atollaire. Cette structure est obtenue par un changement du comportement de la croissance du grenat le long d'un chemin P-T hercynien. Dans un premier stade, le grenat croît rapidement et consume peu le quartz de la matrice. La croissance se fait ainsi le long des jointures des grains de quartz. Dans un second temps, les changements de conditions PT donnent une croissance lente du grenat et une forte consommation du quartz. Le grenat peut ainsi développer sa forme dodécaédrale classique. La troisième partie s'intéresse aux distances de transport par diffusion d'un élément (ici l'argon) durant la haute pression. Pour ce faire, un profile d'âges 40Ar/39Ar sur biotite a été mesuré depuis un veine de haute pression riche en argon jusque dans son encaissant (granitoïd du Monte Mucrone). Le profile montre une répartition des âges suivant une courbe de diffusion. Le transport se fait sur une longueur de deux centimètre avec l'aide d'un fluide. Il est réduit à une échelle millimétrique quand la phase fluide disparaît. Cette étude montre ainsi les difficultés de transport des éléments durant la haute pression ne permettant pas un rééquilibrage de la roche à grande échelle. Kinetics of crystal growth and equilibrium domains in eclogite of the Sesia Zone, Western Alps Darbellay Bastien Institut de Minéralogie et Géochimie Résumé de thèse Les processus de croissance (diffusion des éléments et les réactions d'interface) et les conditions dans lesquelles les minéraux grandissent (température, pression, fluide, composition chimique de la roche), déterminent la texture ainsi que la zonation des minéraux. Cette thèse se focalise, par le biais de textures peu communes, sur trois différents processus impliqués dans la croissance minérale à haute pression (Zone de Sésia, Alpes de l'Ouest, Italie). L'incorporation d'un élément dans une structure minérale ne peut se faire que dans des sites en accord avec la taille et la charge ionique de l'élément. De plus, la balance de charge doit être maintenue dans le minéral. La régularité de la structure cristalline fixe ainsi une limite maximum de concentration d'un élément donné. Les zoisites provenant des métabasites de la région de Cima di Bonze montrent des zonations en sablier caractérisées par une concentration anormale en fer. La zonation se marque par une différente teinte de biréfringence et par un plus grand angle d'extinction que le reste de la zoisite. Une inter-croissance de clinozoisite à l'intérieur de la structure orthorhombique de la zoisite peut ainsi être suspectée. Les analyses XRD (diffraction des rayons x) ainsi que les analyses Raman ne confirment pas cette suspicion. Seules les analyses TEM (microscope à électrons transmis) montrent des défauts de structure pouvant être interprétés comme des modules de clinozoisite. Ils ne peuvent cependant pas être considérés comme une phase thermodynamique. Un nouveau trou d'immiscibilité entre deux zoisite (X ep= 0.1 and Xep = 0.15) a ainsi pu être établi. Dans les métapélites la région du Monte Mucrone, des grenats fortement zonés montrent une évolution texturale singulière. Ils présentent une forme initiale de `champignon' qui se développe pour former une structure atollaire finale. L'étude conjuguée de la structure 3D et des zonations, ainsi que l'établissement d'un model thermodynamique, indiquent que ces structures proviennent de deux épisodes de croissances : (1) La croissance du grenat durant un chemin prograde hercynien (de 525 °C et 6.2 kbar à 640 °C et 9 kbar) permet la formation des textures atollaires. Elles sont le résultat d'une croissance poecilitique initiale suivie d'une croissance idiomorphique du grenat. (2) La structure est rendue plus complexe par la cristallisation d'un grenat homogène tout autour ainsi qu'à l'intérieur du grenat hercynien durant la haute pression alpine (550 °C and 20 kbar). L'arrivée de l'eau durant la haute pression facilite le transport d'éléments et permet une cristallisation rapide du grenat. La diffusion peut être un facteur limitant de la croissance minéralogique. Elle a aussi une grande importance pour la géochronologie. Une veine de haute pression à l'intérieur du granitoïde du Monte Mucrone a été étudiée dans le but de déterminer la distance de diffusion de l'argon. Le profile d'âges 40Ar/39Ar sur biotites, établi de la veine vers le métagranitoïde, suit une courbe de diffusion. Les âges sont élevés proche de la veine (800 Ma) puis décroissent jusqu'à des âges homogènes (170-150 Ma) à deux centimètres de la veine. La présence de fluide, marqué par de hautes concentrations en chlore, permet une diffusion centimétrique. Cependant, la distance est réduite à une échelle millimétrique quand le fluide est absent. Les très faibles distances de diffusion préservent les âges pré-alpins et impliquent un événement géologique pour les âges de 170-150 Ma. Kinetics of crystal growth and equilibrium domains in eclogite of the Sesia Zone, Western Alps Darbellay Bastien Institut de Minéralogie et Géochimie Thesis abstract Rock textures and zonings are the consequence of growth processes (element diffusion and interface reaction) steered by the environment in which they grew (pressure, differential stress, temperature, fluid and rock composition). The thesis presented here focuses on three different topics, each of it dealing with aspects of mineral growth processes during subduction, in a high-pressure environment. All studies were conducted in the Sesia Zone of the Western European Alps, Italy. The first study addresses the crystallography and geochemistry of element incorporation in zoisite, one of the major hydrous minerals found in subduction zone rocks. Elements can be incorporated into a mineral structure only on crystallographic sites that offer enough space for the ion and the overall charge balance has to be maintained. Element concentrations are hence limited. Incorporation of some elements produces complex zoning, including hourglass like patterns, which are the focus of the first contribution. Zoisites from Cima di Bonze (Sesia Zone) show spectacular hourglass zoning defined by Fe-content variations. The hourglass zones have a distinct birefringence and a different extinction angle than the regular part of the zoisite. We show by detailed XRD (X-ray diffraction) and confocal Raman analyses that the high Fe-zones are nevertheless zoisite, and not clinozoisite as one might expect. High resolution TEM (transmission electron microscopy) analyses show planar defects on (100) that can be interpreted as small-scale clinozoisite modules. However, these clinozoisites cannot be interpreted as a distinctive thermodynamic phase and the entire mineral has to be considered as zoisite. The miscibility gap between two zoisites (Xep = 0.1 and Xep = 0.15) can be then definite at 550 ± 50°C and 14 to 20 Kbar. Strongly zoned garnets in quartz rich metapelite from the Monte Mucrone area (Sesia Zone) show evolution form 3D mushroom to atoll structure. The second contribution presents textural investigations, garnet zoning and thermodynamic modeling that demonstrate that atoll garnets are the result of two distinctive growth events. (1) Garnet atoll structure is already formed during a prograde Hercynian path from 525 °C and 6.2 kbar to 640 °C and 9 kbar. It results in an initial poikilitic growth followed by a final idiomorphic growth event. (2) Alpine HP garnet are homogenous (550 °C and 20 kbar) and grew around and also inside the Hercynian garnet. Lack of prograde Alpine garnet and fast growth of the HP garnet is explained by the absence of water during much of the prograde path. Water saturation was only observed towards the end, close towards the peak metamorphic conditions. Diffusion could be a limiting factor for crystal growth. It has also a great importance in geochronology. HP vein inside the metagranitoide of the Monte Mucrone (~300 Ma) was investigated to determine argon diffusion scales during high-pressure metamorphism. 40Ar/39Ar biotite ages profile from the vein toward the metagranodiorite show a diffusion curve: old ages (800 Ma) located close to the vein decrease until homogenous 170-150 Ma ages are obtained, two centimeter away from the vein. Centimeter-scale diffusion occurs with help of a fluid phase marked by high chlorine concentrations. Argon diffusion is reduced to a millimeter scale when free fluid is absent. Very short diffusion distance permits to preserve pre-Alpine ages. The 170-150 Ma ages are considered to be geologic meaningful, probably resulting from the extensional tectonics linked to opening of the Tethian ocean.

THE ADULA NAPPE: STRATIGRAPHY, STRUCTURE AND KINEMATICS OF AN EXHUMED HIGH-PRESSURE NAPPE

This study analyses the stratigraphy, structure and kinematics of the northern part of the Adula nappe of the Central Alps. The Adula nappe is one of the highest basement nappes in the Lower Penninic nappe stack of the Lepontine Dome. This structural position makes possible the investigation of the transition between the Helvetic and North Penninic paleogeographic domains. The Adula nappe is principally composed of crystalline basement rocks. The investigation of the pre-Triassic basement shows that it contains several Palaeozoic detrital metasedimentary formations dated from the Cambrian to the Ordovician. These formations contain also some volcanic or intrusive magmatic rocks. Ordovician metagranites dated at ~450 Ma are also a common rock-type of the Adula basement. These formations underwent Alpine and Variscan deformation and metamorphism. Permian granites (Zervreila orthogneiss, dated at ~290 Ma) have intruded this pre-structured basement in a post-orogenic geodynamic context. Due to their age, the Zervreila orthogneiss are good markers for alpine deformation. The stratigraphy of the Mesozoic and Paleogene sedimentary cover of the Adula nappe is essential to unraveling its pre- orogenic history. The autochthonous cover is assigned to a North Penninic Triassic series that testifies for a transition between the Helvetic and Briançonnais Triassic domains. The Adula domain goes through an emersion during the Middle Jurassic, and is part of a topographic high during the first phase of the Alpine rift. The sediments of the late Middle Jurassic show a drowning phase associated with a tectonic activity and a breccia formation. In the neighbouring domains, coeval with the drowning phase in the Adula domain, a strong extensional crustal delamination and a scattered magmatic activity is associated with the main opening of the North Penninic domain. The Upper Jurassic of the Adula nappe is characterized by a carbonate formation comparable with those in the Helvetic or Subbriaçonnais domains. Flysch s.l. deposition starts probably at the end of the Cretaceous. These sediments are deposited on a large unconformity testifying for a Cretaceous sedimentary gap. The Adula nappe exhibits a very complex structure. This structure is formed by several deformation phases. Two ductile deformations are responsible for the nappe emplacement. The first deformation phase is associated with a folding compatible with a top-to-south movement at the top of the nappe. The second phase is dominant and pervasive throughout the whole nappe. It goes with a strong north vergent folding and the main nappe emplacement. These two phases cause the exhumation and emplacement of a coherent, although pre-structured, piece of continental crust. Two further deformation phases postdate the nappe emplacement. - Ce travail concerne l'étude géologique de la partie nord de la nappe de l'Adula dans les Alpes centrales. La nappe de l'Adula est l'une des nappes cristallines la plus élevée dans la pile des nappes du Pennique inférieur des Alpes lepontines. Cette position particulière permet d'étudier la transition entre les nappes des domaines helvétique et pennique inférieur. La nappe de l'Adula est principalement composée de socle cristallin : l'étude de l'histoire géologique du socle est donc l'un des thèmes de cette recherche. Ce socle contient plusieurs formations métasédimentaires paléozoïques du Cambrien à I'Ordovicien. Ces métasédiments sont issus de formations clastiques comprenant souvent des roches magmatiques volcaniques et intrusives. Ces métasédiments ont subi les cycles orogéniques varisque et alpin. La nappe de l'Adula contient plusieurs corps magmatiques granitiques métamorphisés. Les premiers métagranites sont Ordovicien et témoignent d'un environnement de marge active. Ces granites sont aussi polymétamorphiques. Les deuxièmes métagranites sont représentés par les orthogneiss de type Zervreila. Ce métagranite est d'âge permien (-290 Ma). Il est mis en place dans un contexte tectonique post-orogénique. Ce granite est un maqueur de la déformation alpine car il n'est pas affecté par les orogenèses précédentes, flippy Le contenu stratigraphique des roches mésozoïques et cénozoiques de la couverture sédimentaire de la nappe de l'Adula est'important pour en étudier son histoire pré-alpine. La couverture autochtone est composée d'une série d'âge triasique d'affinité nord-pennique, un faciès qui marque la transition entre les domaines helvétiques et briançonnais au Trias. Le domaine paléogéographique représenté dans la nappe de l'Adula connaît une émersion pendant le Jurassique moyen. Cette émersion marque le commencement du rift dans le domaine alpin. La sédimentation de la fin du Jurassique moyen est marquée par une transgression marine accompagnée par des mouvements tectoniques et la formation d'une brèche. Cette transgression est contemporaine des importants mouvements tectoniques et des manifestations magmatiques dans les unités voisines qui marquent la phase principale d'ouverture du bassin nord-pennique. Le Jurassique supérieur est caractérisé par l'instauration d'une sédimentation carbonatée comparable à celle du domaine helvétique ou subbriançonnais. Une sédimentation flyschoïde, probablement du Crétacé à Tertiaire, est déposée sur une importante discordance qui témoigne d'une lacune au Crétacé. La structure complexe de la nappe de l'Adula témoigne de nombreuses phases de déformation. Ces phases de déformation sont en partie issues de la mise en place de la nappe et de déformations plus tardives. La mise en place de la nappe produit deux phases de déformation ductile : la première produit un plissement compatible avec un cisaillement top-vers-le sud dans la partie supérieure de la nappe; la deuxième produit un intense plissement qui accompagne la mise en place de la nappe vers le nord. Ces deux phases de déformation témoignent d'un mécanisme d'exhumation par déformation ductile d'un bloc cohérent.

Geological transect across the Northwestern Himalaya in eastern Ladakh and Lahul (A model for the continental collision of India and Asia)

The detailed geological mapping and structural study of a complete transect across the northwestern Himalaya allow to describe the tectonic evolution of the north Indian continental margin during the Tethys ocean opening and the Himalayan Orogeny. The Late Paleozoic Tethys rifting is associated with several tectonomagmatic events. In Upper Lahul and SE Zanskar, this extensional phase is recorded by Lower Carboniferous synsedimentary transtensional faults, a Lower Permian stratigraphic unconformity, a Lower Permian granitic intrusion and middle Permian basaltic extrusions (Panjal Traps). In eastern Ladakh, a Permian listric normal fault is also related to this phase. The scarcity of synsedimentary faults and the gradual increase of the Permian syn-rift sediment thickness towards the NE suggest a flexural type margin. The collision of India and Asia is characterized by a succession of contrasting orogenic phases. South of the Suture Zone, the initiation of the SW vergent Nyimaling-Tsarap Nappe corresponds to an early phase of continental underthrusting. To the S, in Lahul, an opposite underthrusting within the Indian plate is recorded by the NE vergent Tandi Syncline. This structure is associated with the newly defined Shikar Beh Nappe, now partly eroded, which is responsible for the high grade (amphibolite facies) regional metamorphism of South Lahul. The main thrusting of the Nyimaling-Tsarap Nappe followed the formation of the Shikar Beh Nappe. The Nyimaling-Tsarap Nappe developed by ductile shear of the upper part of the subducted Indian continental margin and is responsible for the progressive regional metamorphism of SE Zanskar, reaching amphibolite facies below the frontal part of the nappe, near Sarchu. In Upper Lahul, the frontal parts of the Nyimaling-Tsarap and Shikar Beh nappes are separated by a zone of low grade metamorphic rocks (pumpellyite-actinolite facies to lower greenschist facies). At high structural level, the Nyimaling-Tsarap Nappe is characterized by imbricate structures, which grade into a large ductile shear zone with depth. The related crustal shortening is about 87 km. The root zone and the frontal part of this nappe have been subsequently affected by two zones of dextral transpression and underthrusting: the Nyimaling Shear Zone and the Sarchu Shear Zone. These shear zones are interpreted as consequences of the counterclockwise rotation of the continental underthrusting direction of India relative to Asia, which occurred some 45 and 36 Ma ago, according to plate tectonic models. Later, a phase of NE vergent `'backfolding'' developed on these two zones of dextral transpression, creating isoclinal folds in SE Zanskar and more open folds in the Nyimaling Dome and in the Indus Molasse sediments. During a late stage of the Himalayan Orogeny, the frontal part of the Nyimaling-Tsarap Nappe underwent an extension of about 15 km. This phase is represented by two types of structures, responsible for the tectonic unroofing of the amphibolite facies rocks of the Sarchu area: the Sarchu high angle Normal Fault, cutting a first set of low angle normal faults, which have been created by reactivation of older thrust planes related to the Nyimaling-Tsarap Nappe.

TECTONICS OF THE HELVETIC NAPPE SYSTEM (W SWITZERLAND): CONTROL OF STRAIN LOCALIZATION AND BASEMENT-COVER DEFORMATION. Insights from models based on continuum mechanics

The Helvetic nappe system in Western Switzerland is a stack of fold nappes and thrust sheets em-placed at low grade metamorphism. Fold nappes and thrust sheets are also some of the most common features in orogens. Fold nappes are kilometer scaled recumbent folds which feature a weakly deformed normal limb and an intensely deformed overturned limb. Thrust sheets on the other hand are characterized by the absence of overturned limb and can be defined as almost rigid blocks of crust that are displaced sub-horizontally over up to several tens of kilometers. The Morcles and Doldenhom nappe are classic examples of fold nappes and constitute the so-called infra-Helvetic complex in Western and Central Switzerland, respectively. This complex is overridden by thrust sheets such as the Diablerets and Wildhörn nappes in Western Switzerland. One of the most famous example of thrust sheets worldwide is the Glariis thrust sheet in Central Switzerland which features over 35 kilometers of thrusting which are accommodated by a ~1 m thick shear zone. Since the works of the early Alpine geologist such as Heim and Lugeon, the knowledge of these nappes has been steadily refined and today the geometry and kinematics of the Helvetic nappe system is generally agreed upon. However, despite the extensive knowledge we have today of the kinematics of fold nappes and thrust sheets, the mechanical process leading to the emplacement of these nappe is still poorly understood. For a long time geologist were facing the so-called 'mechanical paradox' which arises from the fact that a block of rock several kilometers high and tens of kilometers long (i.e. nappe) would break internally rather than start moving on a low angle plane. Several solutions were proposed to solve this apparent paradox. Certainly the most successful is the theory of critical wedges (e.g. Chappie 1978; Dahlen, 1984). In this theory the orogen is considered as a whole and this change of scale allows thrust sheet like structures to form while being consistent with mechanics. However this theoiy is intricately linked to brittle rheology and fold nappes, which are inherently ductile structures, cannot be created in these models. When considering the problem of nappe emplacement from the perspective of ductile rheology the problem of strain localization arises. The aim of this thesis was to develop and apply models based on continuum mechanics and integrating heat transfer to understand the emplacement of nappes. Models were solved either analytically or numerically. In the first two papers of this thesis we derived a simple model which describes channel flow in a homogeneous material with temperature dependent viscosity. We applied this model to the Morcles fold nappe and to several kilometer-scale shear zones worldwide. In the last paper we zoomed out and studied the tectonics of (i) ductile and (ii) visco-elasto-plastic and temperature dependent wedges. In this last paper we focused on the relationship between basement and cover deformation. We demonstrated that during the compression of a ductile passive margin both fold nappes and thrust sheets can develop and that these apparently different structures constitute two end-members of a single structure (i.e. nappe). The transition from fold nappe to thrust sheet is to first order controlled by the deformation of the basement. -- Le système des nappes helvétiques en Suisse occidentale est un empilement de nappes de plis et de nappes de charriage qui se sont mis en place à faible grade métamorphique. Les nappes de plis et les nappes de charriage sont parmi les objets géologiques les plus communs dans les orogènes. Les nappes de plis sont des plis couchés d'échelle kilométrique caractérisés par un flanc normal faiblement défor-mé, au contraire de leur flanc inverse, intensément déformé. Les nappes de charriage, à l'inverse se caractérisent par l'absence d'un flanc inverse bien défini. Elles peuvent être définies comme des blocs de croûte terrestre qui se déplacent de manière presque rigide qui sont déplacés sub-horizontalement jusqu'à plusieurs dizaines de kilomètres. La nappe de Mordes et la nappe du Doldenhorn sont des exemples classiques de nappes de plis et constitue le complexe infra-helvétique en Suisse occidentale et centrale, respectivement. Ce complexe repose sous des nappes de charriages telles les nappes des Diablerets et du Widlhörn en Suisse occidentale. La nappe du Glariis en Suisse centrale se distingue par un déplacement de plus de 35 kilomètres qui s'est effectué à la faveur d'une zone de cisaillement basale épaisse de seulement 1 mètre. Aujourd'hui la géométrie et la cinématique des nappes alpines fait l'objet d'un consensus général. Malgré cela, les processus mécaniques par lesquels ces nappes se sont mises en place restent mal compris. Pendant toute la première moitié du vingtième siècle les géologues les géologues ont été confrontés au «paradoxe mécanique». Celui-ci survient du fait qu'un bloc de roche haut de plusieurs kilomètres et long de plusieurs dizaines de kilomètres (i.e., une nappe) se fracturera de l'intérieur plutôt que de se déplacer sur une surface frictionnelle. Plusieurs solutions ont été proposées pour contourner cet apparent paradoxe. La solution la plus populaire est la théorie des prismes d'accrétion critiques (par exemple Chappie, 1978 ; Dahlen, 1984). Dans le cadre de cette théorie l'orogène est considéré dans son ensemble et ce simple changement d'échelle solutionne le paradoxe mécanique (la fracturation interne de l'orogène correspond aux nappes). Cette théorie est étroitement lié à la rhéologie cassante et par conséquent des nappes de plis ne peuvent pas créer au sein d'un prisme critique. Le but de cette thèse était de développer et d'appliquer des modèles basés sur la théorie de la méca-nique des milieux continus et sur les transferts de chaleur pour comprendre l'emplacement des nappes. Ces modèles ont été solutionnés de manière analytique ou numérique. Dans les deux premiers articles présentés dans ce mémoire nous avons dérivé un modèle d'écoulement dans un chenal d'un matériel homogène dont la viscosité dépend de la température. Nous avons appliqué ce modèle à la nappe de Mordes et à plusieurs zone de cisaillement d'échelle kilométrique provenant de différents orogènes a travers le monde. Dans le dernier article nous avons considéré le problème à l'échelle de l'orogène et avons étudié la tectonique de prismes (i) ductiles, et (ii) visco-élasto-plastiques en considérant les transferts de chaleur. Nous avons démontré que durant la compression d'une marge passive ductile, a la fois des nappes de plis et des nappes de charriages peuvent se développer. Nous avons aussi démontré que nappes de plis et de charriages sont deux cas extrêmes d'une même structure (i.e. nappe) La transition entre le développement d'une nappe de pli ou d'une nappe de charriage est contrôlé au premier ordre par la déformation du socle. -- Le système des nappes helvétiques en Suisse occidentale est un emblement de nappes de plis et de nappes de chaînage qui se sont mis en place à faible grade métamoiphique. Les nappes de plis et les nappes de charriage sont parmi les objets géologiques les plus communs dans les orogènes. Les nappes de plis sont des plis couchés d'échelle kilométrique caractérisés par un flanc normal faiblement déformé, au contraire de leur flanc inverse, intensément déformé. Les nappes de charriage, à l'inverse se caractérisent par l'absence d'un flanc inverse bien défini. Elles peuvent être définies comme des blocs de croûte terrestre qui se déplacent de manière presque rigide qui sont déplacés sub-horizontalement jusqu'à plusieurs dizaines de kilomètres. La nappe de Morcles and la nappe du Doldenhorn sont des exemples classiques de nappes de plis et constitue le complexe infra-helvétique en Suisse occidentale et centrale, respectivement. Ce complexe repose sous des nappes de charriages telles les nappes des Diablerets et du Widlhörn en Suisse occidentale. La nappe du Glarüs en Suisse centrale est certainement l'exemple de nappe de charriage le plus célèbre au monde. Elle se distingue par un déplacement de plus de 35 kilomètres qui s'est effectué à la faveur d'une zone de cisaillement basale épaisse de seulement 1 mètre. La géométrie et la cinématique des nappes alpines fait l'objet d'un consensus général parmi les géologues. Au contraire les processus physiques par lesquels ces nappes sont mises en place reste mal compris. Les sédiments qui forment les nappes alpines se sont déposés à l'ère secondaire et à l'ère tertiaire sur le socle de la marge européenne qui a été étiré durant l'ouverture de l'océan Téthys. Lors de la fermeture de la Téthys, qui donnera naissance aux Alpes, le socle et les sédiments de la marge européenne ont été déformés pour former les nappes alpines. Le but de cette thèse était de développer et d'appliquer des modèles basés sur la théorie de la mécanique des milieux continus et sur les transferts de chaleur pour comprendre l'emplacement des nappes. Ces modèles ont été solutionnés de manière analytique ou numérique. Dans les deux premiers articles présentés dans ce mémoire nous nous sommes intéressés à la localisation de la déformation à l'échelle d'une nappe. Nous avons appliqué le modèle développé à la nappe de Morcles et à plusieurs zones de cisaillement provenant de différents orogènes à travers le monde. Dans le dernier article nous avons étudié la relation entre la déformation du socle et la défonnation des sédiments. Nous avons démontré que nappe de plis et nappes de charriages constituent les cas extrêmes d'un continuum. La transition entre nappe de pli et nappe de charriage est intrinsèquement lié à la déformation du socle sur lequel les sédiments reposent.

Volatile transport during the crystallization of anatectic melts: Oxygen, boron and hydrogen stable isotope study on the metamorphic complex of Naxos, Greece

Crystallization of anatectic melts in high-temperature metamorphic terrains releases volatile-rich magmas that can be transported into adjacent lithologies. This study addresses the variations in the oxygen, boron and hydrogen isotopic composition of aplite-pegmatite dikes that formed during the crystallization of anatectic melts in regional high-temperature metamorphism on the island of Naxos, Greece, and propagated upward into the overlying sequences of metamorphic schist. The transport distance of these dikes was increased through a significant horizontal component of travel that was imposed by contemporaneous low-angle extensional shearing. Laser fluorination oxygen isotope analyses of quartz, tourmaline, garnet, and biotite mineral separates from the aplite-pegmatite dikes show a progressive rise in delta(18)O values with increasing distance from the core. Oxygen isotope fractionations among quartz, tourmaline, and garnet show temperature variations from > 700degreesC down to similar to400degreesC. This range is considered to reflect isotopic fractionation beginning with crystallization at high temperatures in water-undersaturated conditions and then evolving through lower temperature crystallization and retrograde sub-solidus exchange. Two processes are examined for the cause of the progressive increase in delta(18)O values: (1) heterogeneous delta(18)O sources and (2) fluid-rock exchange between the aplite/pegmatite magmas and their host rock. Although the former process cannot be ruled out, there is as yet no evidence in the exposed sequences on Naxos for the presence of a suitable high delta(18)O magma source. In contrast, a tendency for the delta(18)O of quartz in the aplite/pegmatite dikes to approach that of the quartz in the metamorphic rock suggests that fluid-rock exchange with the host rock may potentially be an important process. Advection of fluid into the magma is examined based on Darcian fluid flow into an initially water-undersaturated buoyantly propagating aplitic dike magma. It is shown that such advective flow could only account for part of the O-18-enrichment, unless it were amplified by repeated injection of magma pulses, fluid recycling, and deformation-assisted post-crystallization exchange. The mechanism is, however, adequate to account for hydrogen isotope equilibration between dike and host rock. In contrast, variations in the delta(11)B values of tourmalines suggest that B-11/B-10 fractionation during crystallization and/or magma degassing was the major control of boron geochemistry rather than fluid-rock interaction and that the boron isotopic system was decoupled from that of oxygen. Copyright (C) 2003 Elsevier Ltd.