Quantifying the tectono-metamorphic evolution of pelitic rocks from a wide range of tectonic settings: mineral compositions in equilibrium

Commonly used thermometer and barometer calibrations are sensitive to mineral assemblage and, thus, bulk-rock composition. Calculated mineral stabilities for an average pelitic rock over a pressure–temperature (PT) range appropriate for normal, thickened, heated and shallowly subducted continental crust (400–900°C at 0.1–3.0 GPa) reveal more than one hundred possible assemblages. Individual phase compositions are dependent on the assemblage in which they belong and combining isopleth sets to represent and reveals several PT-ranges where commonly used mineral thermobarometers are less effective. For example, the garnet-biotite thermometer becomes increasingly P dependent in the absence of muscovite in high T melt-bearing assemblages, and biotite and plagioclase are not stable at pressures appropriate for lower thickened continental crust. Compositional thermobarometers involving equilibration between alternative phases (namely garnet, phengite and omphacite) are presented. Although the equilibrium compositions of phases at any P and T may change significantly as a function of bulk-rock composition, compositional-ratio thermobarometers are typically insensitive to this, unless a pseudo-univariant reaction is crossed and the buffering assemblage is altered. Quantification of the limits of efficacy of various thermobarometers allows the mineralogy of metapelites to be used to precisely determine segments of PT paths and infer their likely tectonic controls. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

---

Forward calculation of zoned garnet growth with limited diffusion transport in the matrix

Summary A forward model is proposed to reproduce the formation of garnet under conditions of sluggish diffusion transport in the matrix. Starting from a matrix consisting of chlorite and quartz, the amount of garnet growth and the chemical composition was calculated at each P–T increment in the system MnO–FeO–MgO–Al₂O₃–SiO₂–H₂O. Sluggish diffusion transport was introduced considering the local equilibrium between garnet surface and the matrix within a given diffusion distance (equilibration volume). Varying the diffusion distance, calculations were performed along the prograde P–T path of the Sambagawa metamorphic belt, Japan. The final size of the garnet grains was largely proportional to the diffusion distance. In contrast to the model without diffusion limitations, a shorter diffusion distance resulted in a rise of the Mg/(Mg + Fe) ratio in garnet before Mn approached zero. These results indicate that the chemical composition trend in zoned garnet from the Sambagawa belt is consistent with growth under sluggish material transport. The calculated amount of garnet growth increases dramatically with temperature. The amount of newly grown natural garnet in the Sambagawa metamorphic rocks was plotted against temperatures, where chemical compositions of garnet were calibrated against temperatures with the Gibbs’ method. This trend was also consistent with the modelled garnet behaviour.     

Calcium diffusivity in alumino-silicate garnets: an experimental and ATEM study

Concentration gradients in calcium are common in metamorphic or magmatic garnets and can be used to determine the timescales of geological processes. However, the kinetics of Ca diffusion in garnet is poorly constrained and experimental studies have to date yielded widely varying diffusion coefficients. In this paper, we describe a new method for generating diffusion profiles in garnet. We incorporated polished and compositionally homogeneous garnet seeds in a finely ground powder of clinopyroxene and garnet. During the experiments (1.3 GPa, 1,050–1,250°C, and ƒO₂ ≤ the graphite-O₂ buffer), the mineral powder partially melted, recrystallized, and formed a 10–50 μm wide overgrowth zone of compositionally distinct garnet around the seeds. Long duration experiments generated measurable relaxation profiles at these seed/overgrowth interfaces. We performed analytical transmission electron microscope traverses across the interfaces in each experiment. Thirteen usable compositional profiles were obtained with characteristic distances of diffusion ranging from 300 to 1,000 nm. From these profiles, Ca–(Fe, Mg) interdiffusion coefficients were retrieved using an analytical solution for the diffusion equation and the data were cast in an Arrhenius relation. Linear regression of the data yields an activation energy Q _{Ca–(Fe, Mg)} equal to 188 ± 48 kJ mol⁻¹ and a frequency factor D ₀ equal to 6.6 × 10⁻¹⁴ m² s⁻¹. Within the compositional range studied, the composition of garnet has no major effect on the Ca–(Fe, Mg) interdiffusion coefficient. The very slow diffusion rate of Ca is in agreement with natural observations indicating that Ca diffuses more slowly than Fe and Mg. The Ca diffusion coefficients derived from this study are not model-dependent and can be used to determine the durations of geological events from Ca relaxation profiles in natural garnets.

---

Metabasites with eclogite facies relics from Variscides in Sardinia,Italy: a review

Metabasites with eclogite facies relics occur in northern Sardinia as massive to strongly foliated lenses or boudins embedded within low- to medium-grade rocks (Anglona) and migmatites (NE Sardinia). U–Pb zircon dating yielded 453 ± 14, 457 ± 2 and 460 ± 5 Ma as the protolith ages; 400 ± 10 and 403 ± 4 Ma have been interpreted as the ages of the HP event and 352 ± 3 and 327 ± 7 Ma as the ages of the main Variscan retrograde events. A pre-eclogite stage is documented by the occurrence of tschermakite, zoisite relics within garnet porphyroblasts (Punta de li Tulchi) and an edenite–andesine inclusion within a relict kyanite porphyroblast (Golfo Aranci). Four main metamorphic stages have been distinguished in the eclogite evolution: (1) eclogite stage, revealed by the occurrence of armoured omphacite relics within garnet porphyroblasts. The Golfo Aranci eclogites also include kyanite, Mg-rich garnet and pargasite; (2) granulite stage, producing orthopyroxene and clinopyroxene–plagioclase symplectites replacing omphacite. At Golfo Aranci, the symplectitic rims around relict kyanite consist of sapphirine, anorthite, corundum and spinel; (3) amphibolite stage, leading to the formation of amphibole–plagioclase kelyphites between garnet porphyroblasts and pyroxene–plagioclase symplectites and to the growth of cummingtonite on orthopyroxene. Tschermakite to Mg-hornblende, plagioclase, cummingtonite, ilmenite, titanite and biotite are coexisting phases; (4) greenschist to sub-greenschist stage, defined by the appearance of actinolite, chlorite, epidote ss, titanite, sericite and prehnite. The following P–T ranges have been estimated for the different stages. Eclogite stage 550–700°C; 1.3–1.7 GPa; granulite stage 650–900°C; 0.8–1.2 GPa, clustering in the range 1.0–1.2 GPa; amphibolite stage 550–740°C; 0.3–0.7 GPa; greenschist stage 300–400°C; 0.2–0.3 GPa. Comparable ranges characterise the other Variscan massifs in Europe; eclogite stage: T = 530–800°C; P from 0.7–1.1 to 1.7 ± 0.3 GPa; granulite stage T = 760–870°C and P from 1.1–1.4 to 7.2–9.9 GPa, clustering around 1.0–1.2 GPa. Whole-rock chemistry: Sardinian eclogites are N- to T-MORB; European ones N- to E-MORB or calc-alkaline.     

Long length scales of element transport during reaction texture development in orthoamphibole-cordierite gneiss: Thor-Odin dome,British Columbia,Canada

First-order factors controlling the textural and chemical evolution of metamorphic rocks are bulk composition and pressure–temperature–time (P–T–t) path. Although it is common to assume that major element bulk composition does not change during regional metamorphism, rocks with reaction textures such as corona structures record evidence for major changes in effective bulk composition (EBC) and therefore provide significant insight into the scale, pathways, and mechanisms of element transport during metamorphism. Quantifying changes in EBC is essential for petrologic applications such as calculation of phase diagrams (pseudosections). The progressive growth of complex corona structures on garnet and Al₂SiO₅ porphyroblasts in orthoamphibole-cordierite gneiss Thor-Odin dome (British Columbia, Canada) reduced the EBC volume of the rock during metamorphism and therefore had a dramatic effect on the evolution of the stable mineral assemblage. These rocks contain a chemical and textural record of metamorphic reactions and preserve 3D networks (reaction pathways) connecting corona structures. These coronal networks record long (>cm) length scales of localized element transport during metamorphism. P–T, T–X, and P–X pseudosections are used to investigate the control of effective bulk composition on phase assemblage evolution. Despite textural complexity and evidence for disequilibrium, mineral assemblages and compositions were successfully modeled and peak metamorphic conditions estimated at 750°C and 9 kbar. These results illustrate how textural and chemical changes during metamorphism can be evaluated using an integrated petrographic and pseudosection approach, highlight the importance of effective bulk composition choice for application of phase equilibria methods in metamorphic rocks, and show how corona structures can be used to understand the scale of compositional change and element transport during metamorphism.     

Chloritoid-Bearing Mineral Assemblages in High-Pressure Metapelites from the Bughea Complex, Leaota Massif (South Carpathians)

High-Mg chloritoid (X_Mg = 0·40–0·47) andrelatively high-Mg staurolite (X_Mg = 0·25–0·28)coexisting with kyanite and garnet were identified in a mica–garnet-richrock associated with very high-pressure eclogites in the BugheaComplex of the Leaota Massif (South Carpathians). Major andtrace element geochemical data for both fresh eclogites andassociated rocks which represent a metasomatic or retrogradealteration rind of the eclogites, indicate a pelitic precursor.Magnesian chloritoid was found as inclusions in garnet as partof a chloritoid–kyanite–garnet assemblage whichis indicative of high-pressure conditions. The host garnet showsa typically prograde chemical zoning pattern. The chloritoid-bearingassemblage is confined to the inner part of the garnet porphyroblasts,whereas the matrix assemblage in equilibrium with Mg-rich garnetrims has exceeded the thermal stability limit of chloritoid.Pressure–temperature pseudosections for simplified compositionsapproaching the rock bulk-chemistry show a high-pressure fieldfor the identified chloritoid-bearing assemblage in good agreementwith pressure–temperature estimates in the CFMASH andKCFMASH chemical subsystems using analysed mineral compositions.The derived pressure–temperature path is clockwise, indicatingoverprinting during exhumation from 1·8 GPa and 580°Cto 1·15 GPa and 620°C, at a water activity approachinga_H2O = 1. These conditions were attained in a subduction mélangeindicating transient thermal perturbations of a subduction channel. KEY WORDS: high-pressure metapelite; Mg-rich chloritoid; P–T path; P–T pseudosection; very high-pressure eclogite     

Metamorphic <Emphasis Type="Italic">P</Emphasis>–<Emphasis Type="Italic">T</Emphasis> paths of the Zanhuang amphibolites and metapelites: constraints on the tectonic evolution of the Paleoproterozoic Trans-North China Orogen

Garnet-bearing metapelites and amphibolites are exposed in the south and middle parts of the Zanhuang complex, which is located in the central segment of the nearly NS-striking Trans-North China Orogen. These rocks preserve three metamorphic mineral assemblages forming at the prograde, peak and post-peak decompression stages. The prograde metamorphic stage (M₁) is represented by mineral inclusions within garnet porphyroblasts, the peak metamorphic stage (M₂) is represented by garnet rims and matrix minerals, whereas the retrograde stage (M₃) is represented by amphibole + plagioclase symplectite rimming garnet porphyroblasts in the amphibolites and biotite + plagioclase symplectite rimming garnet porphyroblasts in the metapelites. All garnet porphyroblasts in the metapelites preserve prograde chemical zoning except for the ubiquitous, quite narrow zones from the underwent post-peak decompression. It has been determined through thermobarometric computation that the metamorphic conditions are 650–710°C at 8.2−9.2 kbar for the M₁ (inclusion) assemblages, >810°C at >12.5 kbar for the metamorphic peak M₂ (matrix) assemblages, and 660–680°C at 4.4–4.5 kbar for the retrograde M₃ (symplectite) assemblages. These rocks are thus determined to have undergone metamorphism with clockwise P–T paths involving nearly isothermal decompression (ITD) segments, which is inferred to be related to the amalgamation of the Eastern and Western Blocks to form the coherent basement of the North China Craton along the Trans-North China Orogen in the late Paleoproterozoic (1.88–1.85 Ga).     

Prograde garnet growth along complex <Emphasis Type="Italic">P–T–t</Emphasis> paths: results from numerical experiments on polyphase garnet from the Wölz Complex (Austroalpine basement)

Garnet in metapelites from the Wölz Complex of the Austroalpine crystalline basement east of the Tauern Window characteristically consists of two growth phases, which preserve a comprehensive record of the geothermal history during polymetamorphism. From numerical modelling of garnet formation, detailed information on the pressure–temperature–time (P–T–t) evolution during prograde metamorphism is obtained. In that respect, the combined influences of chemical fractionation associated with garnet growth, modification of the original growth zoning through intragranular diffusion and the nucleation history on the chemical zoning of garnet as P and T change during growth are considered. The concentric chemical zoning observed in garnet and the homogenous rock matrix, which is devoid of chemical segregation, render the simulation of garnet growth through successive equilibrium states reliable. Whereas the first growth phase of garnet was formed at isobaric conditions of ～3.8 kbar at low heating/cooling rates, the second growth phase grew along a Barrovian P–T path marked with a thermal peak of ～625°C at ～10 kbar and a maximum in P of ～10.4 kbar at ～610°C. For the heating rate during the growth of the second phase of garnet, average rates faster than 50°C Ma^?1 are obtained. From geochronological investigations the first growth phase of garnet from the Wölz Complex pertains to the Permian metamorphic event. The second growth phase grew in the course of Eo-Alpine metamorphism during the Cretaceous.     

Fluid-assisted retrogression of garnet and P–T history of metapelites from HP/UHP metamorphic terrane (Pohorje Mountains, Eastern Alps)

Fluid inclusions in garnet combined with element X-ray mapping, phase equilibrium modelling and conventional thermobarometry have been used to constrain the metamorphic evolution of metapelitic gneiss from the HP/UHP metamorphic terrane of Pohorje Mountains in the Eastern Alps, Slovenia. Retrograde P–T trajectory from ~2.75 GPa and 780°C is constrained by the composition of matrix phengite (6.66 apfu Si) coexisting with garnet cores, kyanite and quartz. The intersection of the X _Prp = 0.25 isopleth for the garnet with the upper stability boundary for K-feldspar in the matrix indicates near-isothermal decompression to ~0.9 GPa at 720°C. Temperatures over 650°C during this stage are corroborated by the high degree of ordering of graphite inclusions associated with Zn, Mg-rich staurolite and phlogopite in the Mg-rich (X _Prp = 0.22–0.25) garnet cores. Majority of garnet porphyroblasts are depleted in Mg (down to X _Prp = 0.09) and enriched in Mn (up to X _Sps = 0.12) along cracks and at their margins. The associated retrograde mineral assemblage comprises Zn, Mg-poor staurolite, muscovite, biotite–siderophyllite, sillimanite and quartz. The onset of the retrogression and the compositional modification of the garnet porphyroblasts were accompanied by the addition of fluid-deposited graphite around older graphite inclusions, probably due to removal of water from a graphite-buffered COH fluid by dissolution in partial silicic melt. Instantaneous expulsion of water near the melt solidus (640°C, max. 0.45 GPa) caused dissolution of the graphite at redox conditions corresponding to 0.25–1.25 logfO₂ units below the QFM buffer, giving rise to a H₂O–CO₂–CH₄ fluid trapped in primary inclusions in Mn-rich, Mg-poor, almandine garnet that reprecipitated within the retrogressed domains. The absence of re-equilibration textures and consistent densities of the fluid inclusions reflect a near-isochoric cooling postdating the near-isothermal decompression. Bulk water content in the metapelite attained 2 wt% during this stage. The low-degree partial melting and extensive hydration due to the release of the internally derived, low-pressure aqueous fluids led to the reset of peak-pressure mineral assemblage.     

Decoding the complex internal chemical structure of garnet porphyroblasts from the Zermatt area,Western Alps

Garnet is a prototypical mineral in metamorphic rocks because it commonly preserves chemical and textural features that can be used for untangling its metamorphic development. Large garnet porphyroblasts may show extremely complex internal structures as a result of a polycyclic growth history, deformation, and modification of growth structures by intra‐ and intercrystalline diffusion. The complex internal structure of garnet porphyroblasts from garnet–phengite schists (GPS) of the Zermatt area (Western Alps) has been successfully decoded. The centimetre‐sized garnet porphyroblasts are composed of granulite facies garnet fragments overgrown by a younger generation of grossular‐rich eclogite facies garnet. The early granulite facies garnet (G‐Grt) formed from low‐P, high‐T metamorphism during a pre‐Alpine orogenic event. The late garnet (E‐Grt) is typical of high‐pressure, low‐temperature (HPLT) metamorphism and can be related to Alpine subduction of the schists. Thus, the garnet of the GPS are polycyclic (polymetamorphic). G‐Grt formation occurred at ~670 MPa and 780°C, E‐Grt formed at ~1.7 GPa and 530°C. The G‐Grt is relatively rich in Prp and poor in Grs, while E‐Grt is rich in Grs and poor in Prp. The Alm content (mol.%) of G‐Grt is 68 of E‐Grt 55. After formation of E‐Grt between and around fragmented G‐Grt at 530°C, the GPS have been further subducted and reached a maximum temperature of 580°C before exhumation started. Garnet composition profiles indicate that the initially very sharp contacts between the granulite facies fragments of G‐Grt and fracture seals of HPLT garnet (E‐Grt) have been modified by cation diffusion. The profiles suggest that Ca did not exchange at the scale of 1 µm, whereas Fe and Mg did efficiently diffuse at the derived maximum temperature of 580°C for the GPS at the scale of 7–8 µm. The Grt–Grt diffusion profiles resulted from spending c. 10 Ma at 530–580°C along the P–T–t path. The measured Grt composition profiles are consistent with diffusivities of log D_MgFe = ?25.8 m²/s from modelled diffusion profiles. Mg loss by diffusion from G‐Grt is compensated by Fe gain by diffusion from E‐Grt to maintain charge balance. This leads to a distinctive Fe concentration profile typical of uphill diffusion.     

Experimental evidence of decompression melting during exhumation of subducted continental crust

Experiments have been carried out on a metagreywacke at 800, 850 and 900°C, in the pressure range 0.5–5 GPa to locate the solidus and the eclogite/amphibolite facies transition in felsic rocks, identify the nature of the reactions responsible for major mineralogical changes, and determine the proportions of phases as a function of pressure. The mineral assemblage phengite + clinopyroxene + garnet + quartz/coesite is stable above 2.3 GPa while biotite + plagioclase + garnet + quartz is stable below 2 GPa. The model reaction for the eclogite/amphibolite facies transition in metagreywackes is:

Simultaneous operation of opposing reaction mechanisms: The influence of matrix heterogeneity on post-kinematic garnet crystallization in an inverted metamorphic sequence

Metapelites from the inverted Barrovian sequence in the Sikkim Himalaya (northeast India) are shown to be largely continuous with respect to their bulk rock compositions, microstructures and pressure–temperature–time–deformation (P–T–t–D) histories. However, the upper garnet–lower staurolite zone demarcates a region of microstructurally anomalous post-kinematic garnet populations contained within strongly segregated matrices. The different microstructures within samples from this region cannot be attributed to differences in their thermobarometric histories or bulk compositions, but are instead interpreted to represent an otherwise unexposed level of the Daling Group that is now exposed along a post-metamorphic thrust splay. These heterogeneous samples contain several discrete garnet populations that progressively crystallized with increasing P–T. Garnet populations that experienced the most protracted growth now form complex polycrystals that exhibit crystallographically controlled and morphologically irregular interfaces adjacent to micaceous and quartzofeldspathic domains respectively. Electron backscatter diffraction indicates that these polycrystalline garnet structures contain numerous coalesced porphyroblasts that are structurally uncorrelated across their grain boundaries. However, a crystallographically preferred orientation at the polycrystal scale is interpreted to derive from epitaxial crystallization of early-formed garnet porphyroblasts on precursor mica. Later-nucleated porphyroblasts within polycrystals preferentially concentrated towards quartzofeldspathic domains, with the overall nucleation distribution likely controlled by a complex interplay between chemical heterogeneities, strain partitioning and epitaxial crystallization. The subsequent growth of these polycrystals was equally spatially heterogeneous; it was moderated by differences in the efficiency of grain boundary transfer between quartzofeldspathic and micaceous domains that precluded thin section-scale chemical equilibration. In contrast to samples from Sikkim containing more typical porphyroblastic populations in continuous and disseminated matrices, heterogeneous availability of garnet-forming components within this strongly layered matrix is shown to have resulted in grain-scale variations in growth rates and the spatial juxtapositioning of interface-controlled microstructures and locally equilibrated chemical compositions with those that were transport controlled.     

High-Pressure (~2000 MPa) Kyanite- and Glaucophane-bearing Pelitic Schist and Eclogite from Cordillera de la Costa Belt, Venezuela

Cretaceous melange of the Cordillera de la Costa belt, north–centralVenezuela, there are knockers of eclogite, barroisite-bearingeclogite, and pelitic glaucophane schist. These occur in a metamorphicmelange matrix that locally consists of marble, serpentinite,amphibolite, actinolite schist, feldspathic schist and gneiss,graphitic schist, chloritoid schist, and garnet-bearing micaschist. The protoliths for these various rock types exhibita wide age range (Cambrian to Early Cretaceous?). Recently discoveredknockers of pelitic glaucophane schist contain Mg-glaucophane+ paragonite + kyanite + garnet + talc + graphite + rutile +quartz. The coexistence of kyanite and Mg-glaucophane suggestsminimum P 2000 MPa at T > 600°C. Eclogite knockers fromthe same outcrop contain garnet and clinopyroxene which yield500°C for cores, 700°C for rims, and P 1200 MPa. Theassemblage garnet–biotite–phengite–albitewithin schists of the melange matrix of this locality indicatesmetamorphic conditions of T = 450–520°C at P = 1800MPa. Because all lithologies in this outcrop record high-P conditions,this metamorphic melange formed before or during peak metamorphismin a mid-Cretaceous subduction zone. KEY WORDS: geothermobarometry; high-P pelitic schist; eclogite; Puerto Cabello; Venezuela     

P–T–t evolution of spinel–cordierite–garnet gneisses from the Sauwald Zone (Southern Bohemian Massif, Upper Austria): is there evidence for two independent late-Variscan low-P/high-T events in the Moldanubian Unit?

The Sauwald Zone, located at the southern rim of the Bohemian Massif in Upper Austria, belongs to the Moldanubian Unit. It exposes uniform biotite + plagioclase ± cordierite paragneisses that formed during the post-collisional high-T/low-P stage of the Variscan orogeny. Rare metapelitic inlayers contain the mineral assemblage garnet + cordierite + green spinel + sillimanite + K-feldspar + plagioclase + biotite + quartz. Mineral chemical and textural data indicate four stages of mineral growth: (1) peak assemblage as inclusions in garnet (stage 1): garnet core + cordierite + green spinel + sillimanite + plagioclase (An_35–65); (2) post-peak assemblages in the matrix (stages 2, 3): cordierite + spinel (brown-green and brown) ± sillimanite ± garnet rim + plagioclase (An_10–45); and (3) late-stage growth of fibrolite, muscovite and albite (An_0–15) during stage 4. Calculation of the P–T conditions of the peak assemblage (stage 1) yields 750–840°C, 0.29–0.53 GPa and for the stage 2 matrix assemblage garnet + cordierite + green spinel + sillimanite + plagioclase 620–730°C, 0.27–0.36 GPa. The observed phase relations indicate a clockwise P–T path, which terminates below 0.38 GPa. The P–T evolution of the Sauwald Zone and the Monotonous Unit are very similar, however, monazite ages of the former are younger (321 ± 9 Ma vs. 334 ± 1 Ma). This indicates that high-T/low-P metamorphism in the Sauwald Zone was either of longer duration or there were two independent phases of late-Variscan low-P/high-T metamorphism in the Moldanubian Unit.     

Zr-in-rutile thermometry in HP/UHP eclogites from Western China

Four Zr-in-rutile thermometry calibrations are applied to eclogites from Western China. Here, we show that if rutile grows in equilibrium with Qtz and Zrn, and is isolated inside garnet, it preserves its Zr composition and does not undergo compositional change due to cation exchange with the host garnet. It thus preserves the composition for the P–T conditions of its formation and the growth zoning of the host garnet. For the HP/UHP metamorphic temperature, the Tomkins et al. (J Metamorph Geol 25:703–713, 2007) calibration yields temperatures that agree well with previous studies, whereas the other three calibrations (Zack et al. in Contrib Mineral Petrol 148:471–488, 2004; Watson et al. in Contrib Mineral Petrol 151:413–433, 2006; Ferry and Watson in Contrib Mineral Petrol in 154:429–437, 2007), which do not include a pressure correction, give systematically lower temperatures. Zr contents of rutile inclusions within garnet show systematic decrease from garnet core to rim. The rutile inclusions in garnet rims contain the lowest Zr content, similar to that in the matrix. Analyses confirm that the pressure plays a significant role in modifying the primary temperature dependence of the Zr content of rutile. Rutiles trapped in garnets are unable to re-equilibrate easily during retrogression, but those in the matrix can do so, providing retrograde P–T path information.     

Subduction-related P–T path for eclogites and garnet glaucophanites from the Samaná Peninsula basement complex, northern Hispaniola

Eclogite and garnet glaucophanite lenses from the Punta Balandra unit of the Samaná basement complex (northern Hispaniola) preserve information of the early metamorphic and tectonic history of subduction in the Caribbean island-arc and its collision with the North America plate. For this reason, P–T paths were reconstructed from the interpretation of meso- and microfabrics, mineral assemblages and chemistry, with the aid of equilibrium phase diagrams calculated for specific bulk compositions in the CKNFMASH system and isopleths for selected solution end-members. The obtained results suggest that the subduction-related prograde path evolved from garnet-free and garnet-bearing lawsonite-blueschist facies, to phengite eclogite facies conditions at P=22–24 kbar and T=610–625°C, with a probable intermediate stage of low-P lawsonite eclogite facies. The subsequent retrograde P–T path entered the epidote-blueschist (garnet-free) facies and ended within the greenschist facies field, similar to the prograde evolution at low-P. Eclogites and garnet glaucophanites formed in a subduction zone in which oceanic lithosphere was subducted WSW/W beneath the Caribbean plate.     

Thermochemistry of sulfide liquids III: Ni-bearing liquids at 1 bar

Experiments were performed in the system O–S–Fe–Ni designed to extend our understanding of the chemistry of sulfide liquids. Results indicate that adding nickel to Fe-rich sulfide liquids in equilibrium with silicate liquids extends their stability field to much higher oxygen fugacities and lower sulfur fugacities. Increasing Ni/Fe at a given temperature and sulfur and oxygen fugacity is accompanied by a significant decrease in the oxygen content of the sulfide liquid. Results of these experiments are combined with data from the literature to calibrate an associated regular solution model for O–S–Fe–Ni liquids. This model represents a complete refit of the associated regular solution model of Kress (Contrib Mineral Petrol 139:316–325, 2000). The resulting model is combined with the olivine solution model of Hirschmann (Am Mineral 76:1232–1248, 1991) to explore the effect of variations in oxygen and sulfur fugacities on the distribution of Fe and Ni between olivine and sulfide liquid. Predicted olivine–sulfide distribution trends parallel those observed by Gaetani and Grove (Geochim Cosmochim Acta 61:1829–1846, 1997), Gaetani and Grove (Earth Planet Sci Lett 169:147–163, 1999), Brenan and Caciagli (Geochim Cosmochim Acta 64:307–320, 2000) and Brenan (Geochim Cosmochim Acta 67:2663–2681, 2003), but are systematically offset toward lower predicted Ni in the sulfide. Nevertheless our results are consistent with the assertion that low K _D ^os values in magmatic ore deposits such as the J-M Reef reflect high iron contents in the sulfides combined with relatively high oxygen fugacities.

---

Population-wide garnet growth zoning revealed by LA-ICP-MS mapping: implications for trace element equilibration and syn-kinematic deformation during crystallisation

The strong partitioning of many trace elements into garnet and their slow diffusivities in both garnet and the rock matrix means that their distribution may record valuable petrogenetic information not documented by major elements in metamorphic rocks. Complex trace element growth zoning in garnet porphyroblasts from a garnet-grade metapelite from the Barrovian sequence of the Sikkim Himalaya is assessed using quantified LA-ICP-MS raster mapping coupled with X-ray micro-computed tomography. The data document systematic changes in the zoning patterns from early- to late-nucleated crystals, and also suggest that the REE+Y chemistry incorporated into garnet is dependent on persistent disequilibrium in the rock volume. There is evidence for HREE+Y diffusion haloes surrounding growing garnets and a heterogeneous HREE+Y distribution in the rock matrix. Annuli superimposed on oscillatory zoning are not consistent with formation during some rock-wide event, but are dependent on the spatial disposition of the garnet. Annuli are interpreted to reflect an integrated history of varying growth rates and the incorporation of pre-existing heterogeneities due to relatively slow matrix diffusivities. Conversely, smooth zoning of many transition metals indicate that their distribution in garnet may be controlled by equilibrium partitioning between garnet and the rock matrix. Significant rotation of garnet porphyroblasts during growth is revealed due to immobility of Cr over the duration of the crystallisation interval and overprinting of the heterogenous precursor Cr distribution. Strain rate estimates derived from this zoning are on the order of \(10^{-11}\)–\(10^{-12}\, \hbox {s}^{-1}\).     

Pressure-temperature paths from P-T pseudosections and zoned garnets: potential, limitations and examples from the Zanskar Himalaya, NW India

The extraction of P-T histories from metamorphic rocks provides a valuable dataset for the elucidation of the tectonic mechanisms for orogeny. While continued re-equilibration frequently obliterates early information, garnet zonation and inclusion assemblages can often surmount this problem. The task is more difficult in high variance assemblages or if inclusions are not preserved, but one approach is to use pseudosections that are specific to the bulk composition of a given rock. In the latter case, the compositions and abundances of all the minerals are fixed at a given P-T point such that, if the effective bulk composition is known, the garnet composition alone can be used to reconstruct the history. Here, we explore this approach using examples from the Zanskar Himalaya, NW India. Pseudosections have been calculated for four pelitic to semipelitic rocks from the Zanskar Himalaya and have been contoured for garnet composition. The calculations adequately model the mineral assemblages in the rocks and predict the presence of chlorite in the early assemblage where chlorite is found as inclusions within garnet. Moreover, the pseudosections successfully model the garnet core compositions, with all three independent compositional contours overlapping at a single pressure and temperature. This occurs at ∼550 °C and at pressures varying from 3–7 kbar for the four rocks studied. We have been less successful, however, at modelling garnet compositions beyond the cores because fractionation of the effective bulk composition is caused by garnet growth itself. However, in this case, a combination of the␣pseudosection and conventional thermobarometry using␣Fe-Ti inclusions and matrix phases allows us to reconstruct␣the entire P-T history. The resulting P-T paths record burial of 3–5 kbar without significant temperature increase followed by isobaric heating of 50–100 °C. This evolution is consistent with Himalayan collision in the early Tertiary but a combination of the P-T data presented here and published geochronological data suggests renewed thrusting south of the suture zone in the Oligocene. In addition, the data demonstrate that no extra heat source is required to cause melting of the Himalayan crust in the Miocene. While melting could have occurred both by dehydration during decompression or in the presence of a fluid, the lack of garnet resorption does suggest decompression was rapid and followed quickly by cooling. This scenario favours melting by decompression. Received: 17 July 1997 / Accepted: 6 April 1998     

Geochemical syntheses among the cratonic,off-cratonic and orogenic garnet peridotites and their tectonic implications

Garnet-bearing mantle peridotites, occurring as either xenoliths in volcanic rocks or lenses/massifs in high-pressure and ultrahigh-pressure terrenes within orogens, preserve a record of deep lithospheric mantle processes. The garnet peridotite xenoliths record chemical equilibrium conditions of garnet-bearing mineral assemblage at temperatures (T) ranging from ~700 to 1,400°C and pressures (P) > 1.6–8.9 GPa, corresponding to depths of ~52–270 km. A characteristic mineral paragenesis includes Cr-bearing pyropic garnet (64–86 mol% pyrope; 0–10 wt% Cr₂O₃), Cr-rich diopside (0.5–3.5 wt% Cr₂O₃), Al-poor orthopyroxene (0–5 wt% Al₂O₃), high-Cr spinel (Cr/(Cr + Al) × 100 atomic ratio = 2–86) and olivine (88–94 mol% forsterite). In some cases, partial melting, re-equilibration involving garnet-breakdown, deformation, and mantle metasomatism by kimberlitic and/or carbonatitic melt percolations are documented. Isotope model ages of Archean and Proterozoic are ubiquitous, but Phanerozoic model ages are less common. In contrast, the orogenic peridotites were subjected to ultrahigh-pressure (UHP) metamorphism at temperature ranging from ~700 to 950°C and pressure >3.5–5.0 GPa, corresponding to depths of >110–150 km. The petrologic comparisons between 231 garnet peridotite xenoliths and 198 orogenic garnet peridotites revealed that (1) bulk-rock REE (rare earth element) concentrations in xenoliths are relatively high, (2) clinopyroxene and garnet in orogenic garnet peridotites show a highly fractionated REE pattern and Ce-negative anomaly, respectively, (3) Fo contents of olivines for off-cratonic xenolith are in turn lower than those of orogenic garnet and cratonic xenolith but mg-number of garnet for orogenic is less than that of off-cratonic and on-cratonic xenolith, (4) Al₂O₃, Cr₂O₃, CaO and Cr# of pyroxenes and chemical compositions of whole rocks are very different between these garnet peridotites, (5) orogenic garnet peridotites are characterized by low T and high P, off-cratonic by high T and low P, and cratonic by medium T and high P and (6) garnet peridotite xenoliths are of Archean or Proterozoic origin, whereas most of orogenic garnet peridotites are of Phanerozoic origin. Taking account of tectonic settings, a new orogenic garnet peridotite exhumation model, crust-mantle material mixing process, is proposed. The composition of lithospheric mantle is additionally constrained by comparisons and compiling of the off-cratonic, on-cratonic and orogenic garnet peridotite.