Multistage progressive evolution of rare osumilite‐bearing assemblages preserved in ultrahigh‐temperature granulites from In Ouzzal (Hoggar,Algeria)

Osumilite is reported in Palaeoproterozoic Al–Mg‐rich granulites from the Khanfous area (Tekhamalt, In Ouzzal, Hoggar, Algeria). The main peak assemblages are osumilite + sapphirine + biotite + orthopyroxene + sillimanite and osumilite + orthopyroxene + sillimanite + quartz ± biotite (±K‐feldspar) in silica‐deficient and silica‐saturated granulites respectively. Osumilite coexists with F‐rich biotite (X_F ≈ 0.6). The observed microstructures, the mass balance of metamorphic reactions and P–T pseudosections modelled for bulk‐rock and reaction‐microdomain compositions indicate a clockwise P–T metamorphic evolution at ultrahigh temperatures, without substantial post‐peak deformation. The peak P–T conditions recorded by the osumilite‐bearing assemblages are 8.5–9.0 kbar and 930–980 °C. During retrogression, osumilite was partially or totally replaced by fine‐grained pseudomorphs of cordierite + orthopyroxene + K‐feldspar + quartz at ~7 kbar and ~850 °C. This study confirms that osumilite can occur only in Mg‐rich metamorphic rocks that experienced ultrahigh‐temperature metamorphism under anhydrous conditions. In the presence of a hydrous fluid, it is replaced, even at high temperatures, by cordierite‐bearing assemblages. This important feature explains the rarity of osumilite in granulite facies rocks and its common replacement by cordierite + orthopyroxene + K‐feldspar + quartz pseudomorphs. The peak conditions suggest that a delamination of the lithospheric mantle underneath the In Ouzzal crust brought the asthenosphere close to the Mohorovi?i? discontinuity.     

Two successive phases of ultrahigh temperature metamorphism in Rogaland,S. Norway: Evidence from Y‐in‐monazite thermometry

In Rogaland, South Norway, a polycyclic granulite facies metamorphic domain surrounds the late‐Sveconorwegian anorthosite–mangerite–charnockite (AMC) plutonic complex. Integrated petrology, phase equilibria modelling, monazite microchemistry, Y‐in‐monazite thermometry, and monazite U–Th–Pb geochronology in eight samples, distributed across the apparent metamorphic field gradient, imply a sequence of two successive phases of ultrahigh temperature (UHT) metamorphism in the time window between 1,050 and 910 Ma. A first long‐lived metamorphic cycle (M1) between 1,045 ± 8 and 992 ± 11 Ma is recorded by monazite in all samples. This cycle is interpreted to represent prograde clockwise P–T path involving melt production in fertile protoliths and culminating in UHT conditions of ~6 kbar and 920°C. Y‐in‐monazite thermometry, in a residual garnet‐absent sapphirine–orthopyroxene granulite, provides critical evidence for average temperature of 931 and 917°C between 1,029 ± 9 and 1,006 ± 8 Ma. Metamorphism peaked after c. 20 Ma of crustal melting and melt extraction, probably supported by a protracted asthenospheric heat source following lithospheric mantle delamination. Between 990 and 940 Ma, slow conductive cooling to 750–800°C is characterized by monazite reactivity as opposed to silicate metastability. A second incursion (M2) to UHT conditions of ~3.5–5 kbar and 900–950°C, is recorded by Y‐rich monazite at 930 ± 6 Ma in an orthopyroxene–cordierite–hercynite gneiss and by an osumilite gneiss. This M2 metamorphism, typified by osumilite paragenesis, is related to the intrusion of the AMC plutonic complex at 931 ± 2 Ma. Thermal preconditioning of the crust during the first UHT metamorphism may explain the width of the aureole of contact metamorphism c. 75 Ma later, and also the rarity of osumilite‐bearing assemblages in general.     

Ultrahigh‐temperature osumilite gneisses in southern Madagascar record combined heat advection and high rates of radiogenic heat production in a long‐lived high‐T orogen

We report the discovery of osumilite in ultrahigh‐temperature (UHT) metapelites of the Anosyen domain, southern Madagascar. The gneisses equilibrated at ~930°C/0.6 GPa. Monazite and zircon U–Pb dates record 80 Ma of metamorphism. Monazite compositional trends reflect the transition from prograde to retrograde metamorphism at 550 Ma. Eu anomalies in monazite reflect changes in fO₂ relative to quartz–fayalite–magnetite related to the growth and breakdown of spinel. The ratio Gd/Yb in monazite records the growth and breakdown of garnet. High rates of radiogenic heat production were the primary control on metamorphic grade at the regional scale. The short duration of prograde metamorphism in the osumilite gneisses (<29 ± 8 Ma) suggests that a thin mantle lithosphere (<80 km) or advective heating may have also been important in the formation of this high‐T, low‐P terrane.     

Neoproterozoic evolution and Cambrian reworking of ultrahigh temperature granulites in the Eastern Ghats Province,India

The time‐scales and P–T conditions recorded by granulite facies metamorphic rocks permit inferences about the geodynamic regime in which they formed. Two compositionally heterogeneous cordierite–spinel‐bearing granulites from Vizianagaram, Eastern Ghats Province (EGP), India, were investigated to provide P–T–time constraints using petrography, phase equilibrium modelling, U–Pb geochronology, the rare earth element composition of zircon and monazite, and Ti‐in‐zircon thermometry. These ultrahigh temperature (UHT) granulites preserve discrete compositional layering in which different inferred peak assemblages are developed, including layers bearing garnet–sillimanite–spinel, and others bearing orthopyroxene–sillimanite–spinel. These mineral associations cannot be reproduced by phase equilibrium modelling of whole‐rock compositions, indicating that the samples became domainal on a scale less than that of a thin section, even at UHT conditions. Calculation of the P–T stability fields for six compositional domains within which the main rock‐forming minerals are considered to have attained equilibrium suggests peak metamorphic conditions of ~6.8–8.3 kbar at ~1,000°C. In most of these domains, the subsequent evolution resulted in the growth of cordierite and final crystallization of melt at an elevated (residual) H₂O‐undersaturated solidus, consistent with <1 kbar of decompression. Concordant U–Pb ages obtained by SHRIMP from zircon (spread 1,050–800 Ma) and monazite (spread 950–800 Ma) demonstrate that crystallization of these minerals occurred during an interval of c. 250 Ma. By combining LA‐ICP‐MS U–Pb zircon ages with Ti‐in‐zircon temperatures from the same analysis sites, we show that the crust may have remained above 900°C for a minimum of c. 120 Ma between c. 1,000 and c. 880 Ma. Overall, the results suggest that, in the interval 1,050 to 800 Ma, the evolution of the Vizianagaram granulites culminated with UHT conditions from c. 1,000 Ma to c. 880 Ma, associated with minor decompression, before further zircon crystallization at c. 880–800 Ma during cooling to the solidus. However, these rocks are adjacent to the Paderu–Anantagiri–Salur crustal block to the NW that experienced counterclockwise P–T–t paths, and records similar UHT peak metamorphic conditions (7–8 kbar, ~950°C) followed by near‐isobaric cooling, and has a similar chronology during the Neoproterozoic. The limited decompression inferred at Vizianagaram may be explained by partial exhumation due to thrusting of this crustal block over the adjacent Paderu–Anantagiri–Salur crustal block. The residual granulites in both blocks have high concentrations of heat‐producing elements and likely remained hot at mid‐crustal depths throughout a period of relative tectonic quiescence in the interval 800–550 Ma. During the Cambrian Period, the EGP was located in the hinterland of the Denman–Pinjarra–Prydz orogen. A later concordant population of zircon dated at 511 ± 6 Ma records crystallization at temperatures of ~810°C. This age may record a low‐degree of melting due to limited influx of fluid into hot, weak crust in response to convergence of the Crohn craton with a composite orogenic hinterland comprising the Rayner terrane, EGP, and cratonic India.     

Metamorphism of high‐P metagreywacke from the Eastern Himalayan syntaxis: phase equilibria and P–T path

High‐pressure (HP) metagreywacke from the Namche Barwa Complex, Eastern Himalayan Syntaxis (EHS), consists of garnet, biotite, plagioclase, quartz, rutile and ilmenite with or without K‐feldspar, sillimanite, cordierite, spinel and orthopyroxene. Two types of metagreywacke are recognized: medium‐temperature (MT) and high‐temperature (HT) types. Garnet in the MT metagreywacke shows significant growth zoning and contains lower MgO than the weakly zoned garnet in the HT metagreywacke. Petrographic observations and phase equilibria modelling for four representative samples indicate that both types of metagreywacke experienced clockwise P–T paths subdivided into three stages: stage I is the pre‐peak prograde to pressure peak (P_max) stage characterized by progressive increase in P–T conditions. The P_max conditions are estimated using the garnet composition with maximum CaO, being 12.5–13.5 kbar and 685–725 °C for the MT metagreywacke, and 15–16 kbar and 825–835 °C for the HT one. Stage II is the post‐P_max decompression with heating or near‐isothermal to T_max stage and the T_max conditions, constrained using the garnet compositions with maximum MgO, are 11 kbar and 760 °C for the MT metagreywacke, and ~12 kbar and 830–845 °C for the HT one. The modelled mineral assemblages at T_max are garnet + biotite + K‐feldspar + rutile + plagioclase ± ilmenite in the presence of melt for both types of metagreywacke, consistent with the petrographic observations. Stage III is the post‐T_max retrograde metamorphism, characterized by decompression and cooling. The modelling suggests that the melts with high Na/K ratios (1.7–5.2) have been produced during stages I and II, which could be responsible for the formation of sodium‐rich leucogranites. This study and previous results indicate that the Higher Himalayan Crystallines in the EHS consist of MT–HP and HT–HP metamorphic units separated by a speculated tectonic contact. Petrological and structural discontinuities within the EHS cannot be easily interpreted with ‘tectonic aneurysm’ model.     

Significance of post‐peak metamorphic reaction microstructures in the ultrahigh temperature Eastern Ghats Province,India

Ultrahigh temperature (UHT) granulites in the Eastern Ghats Province (EGP) have a complex P–T–t history. We review the P–T histories of UHT metamorphism in the EGP and use that as a framework for investigating the P–T–t history of Mg–Al‐rich granulites from Anakapalle, with the express purpose of trying to reconcile the down‐pressure‐dominated P–T path with other UHT localities in the EGP. Mafic granulite that is host to Mg–Al‐rich metasedimentary granulites at Anakapalle has a protolith age of c. 1,580 Ma. Mg–Al‐rich metasedimentary granulites within the mafic granulite at Anakapalle were metamorphosed at UHT conditions during tectonism at 960–875 Ma, meaning that the UHT metamorphism was not the result of contact metamorphism from emplacement of the host mafic rock. Reworking occurred during the Pan‐African (c. 600–500 Ma) event, and is interpreted to have produced hydrous assemblages that overprint the post‐peak high‐T retrograde assemblages. In contrast to rocks elsewhere in the EGP that developed post‐peak cordierite, the metasedimentary granulites at Anakapalle developed post‐peak, generation ‘2’ reaction products that are cordierite‐absent and nominally anhydrous. Therefore, rocks at Anakapalle offer the unique opportunity to quantify the pressure drop that occurred during so‐called M2 that affected the EGP. We argue that M2 is either a continuation of M1 and that the overall P–T path shape is a complex counter‐clockwise loop, or that M1 is an up‐temperature counter‐clockwise deviation superimposed on the M2 path. Therefore, rather than the rocks at Anakapalle having a metamorphic history that is apparently anomalous from the rest of the EGP, we interpret that other previously studied localities in the EGP record a different part of the same P–T path history as Anakapalle, but do not preserve a significant record of pressure decrease. This is due either to the inability of refractory rocks to extensively react to produce a rich mineralogical record of pressure decrease, or because the earlier high‐P part of the rocks history was erased by the M1 loop. Irrespective of the specific scenario, models for the tectonic evolution of the EGP must take the substantial pressure decrease during M2 into account, as it is probable the P–T record at Anakapalle is a reflection of tectonics affecting the entire province.     

Refining the P–T records of UHT crustal metamorphism

Ultra‐high‐temperature (UHT) metamorphism occurs when the continental crust is subjected to temperatures of greater than 900 °C at depths of 20–40 km. UHT metamorphism provides evidence that major tectonic processes may operate under thermal conditions more extreme than those generally produced in numerical models of orogenesis. Evidence for UHT metamorphism is recorded in mineral assemblages formed in magnesian pelites, supported by high‐temperature indicators including mesoperthitic feldspar, aluminous orthopyroxene and high Zr contents in rutile. Recent theoretical, experimental and thermodynamic data set constraints on metamorphic phase equilibria in FMAS, KFMASH and more complex chemical systems have greatly improved quantification of the P–T conditions and paths of UHT metamorphic belts. However, despite these advances key issues that remain to be addressed include improving experimental constraints on the thermodynamic properties of sapphirine, quantifying the effects of oxidation state on sapphirine, orthopyroxene and spinel stabilities and quantifying the effects of H₂O–CO₂ in cordierite on phase equilibria and reaction texture analysis. These areas of uncertainty mean that UHT mineral assemblages must still be examined using theoretical and semi‐quantitative approaches, such as P(–T)–μ sections, and conventional thermobarometry in concert with calculated phase equilibrium methods. In the cases of UHT terranes that preserve microtextural and mineral assemblage evidence for steep or ‘near‐isothermal’ decompression P–T paths, the presence of H₂O and CO₂ in cordierite is critical to estimates of the P–T path slopes, the pressures at which reaction textures have formed and the impact of fluid infiltration. Many UHT terranes have evolved from peak P–T conditions of 8–11 kbar and 900–1030 °C to lower pressure conditions of 8 to 6 kbar whilst still at temperature in the range of 950 to 800 °C. These decompressional P–T paths, with characteristic dP/dT gradients of ～25 ± 10 bar °C^?1, are similar in broad shape to those generated in deep‐crustal channel flow models for the later stages of orogenic collapse, but lie at significantly higher temperatures for any specified pressure. This thermal gap presents a key challenge in the tectonic modelling of UHT metamorphism, with implications for the evolution of the crust, sub‐crustal lithosphere and asthenospheric mantle during the development of hot orogens.     

Phase equilibria modelling of kyanite‐bearing anatectic paragneisses from the central Grenville Province

Kyanite‐bearing paragneisses from the Manicouagan Imbricate Zone and its footwall (high‐P belt of the central Grenville Province) preserve evidence of partial melting with development of metamorphic textures involving biotite–garnet ± kyanite ± plagioclase ± K‐feldspar–quartz. Garnet in these rocks displays a variety of zoning patterns with respect to Ca. Pseudosection modelling in the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O (NCKFMASHTO) system using measured bulk rock compositions accounts for the textural evolution of two aluminous and two sub‐aluminous samples from the presumed thermal peak to conditions at which retained melt solidified. The prograde features are best explained by pseudosections calculated with compositions to account for melt loss. The intersection of isopleths of grossular content and Fe/(Fe + Mg) relating to large porphyroblasts of garnet provide constraints on the P–T conditions of the metamorphic peak. These P–T estimates are considered to be minima because of the potential for diffusional modification of the composition of garnet at high‐T and during the early stages of cooling. However, they are consistent with textural observations and pseudosection topology, with peak assemblages best preserved in rocks for which the calculated pseudosections predict only small changes in mineral proportions in the P–T interval, in which retrograde reactions are inferred to have occurred between the thermal peak and the solidus. Maximum P–T conditions (14.5–15.5 kbar and 840–890 °C) and steep retrograde P–T paths inferred for rocks from the Manicouagan Imbricate Zone are comparable with those determined for mafic rocks from the same area. In contrast, maximum P–T conditions of 12.5–13 kbar and 815–830 °C and flatter P–T paths are inferred for the rocks of the footwall to the Manicouagan Imbricate Zone. The general consistency between textures, mineral compositions and the topologies of the calculated pseudosections suggests that the pseudosection approach is an appropriate tool for inferring the P–T evolution of high‐P anatectic quartzo‐feldspathic rocks.     

Integrated pressure–temperature–time constraints for the Tso Morari dome (Northwest India): implications for the burial and exhumation path of UHP units in the western Himalaya

Northward subduction of the leading edge of the Indian continental margin to depths greater than 100 km during the early Eocene resulted in high‐pressure (HP) quartz‐eclogite to ultrahigh‐pressure (UHP) coesite–eclogite metamorphism at Tso Morari, Ladakh Himalaya, India. Integrated pressure–temperature–time determinations within petrographically well‐constrained settings for zircon‐ and/or monazite‐bearing assemblages in mafic eclogite boudins and host aluminous gneisses at Tso Morari uniquely document segments of both the prograde burial and retrograde exhumation path for HP/UHP units in this portion of the western Himalaya. Poikiloblastic cores and inclusion‐poor rims of compositionally zoned garnet in mafic eclogite were utilized with entrapped inclusions and matrix minerals for thermobarometric calculations and isochemical phase diagram construction, the latter thermodynamic modelling performed with and without the consideration of cation fractionation into garnet during prograde metamorphism. Analysis of the garnet cores document (M1) conditions of 21.5 ± 1.5 kbar and 535 ± 15 °C during early garnet growth and re‐equilibration. Sensitive high resolution ion microprobe (SHRIMP) U–Pb analysis of zircon inclusions in garnet cores yields a maximum age determination of 58.0 ± 2.2 Ma for M1. Peak HP/UHP (M2) conditions are constrained at 25.5–27.5 kbar and 630–645 °C using the assemblage garnet rim–omphacite–rutile–phengite–lawsonite–talc–quartz (coesite), with mineral compositional data and regional considerations consistent with the upper P–T bracket. A SHRIMP U–Pb age determination of 50.8 ± 1.4 Ma for HP/UHP metamorphism is given by M2 zircons analysed in the eclogitic matrix and that are encased in the garnet rim. Two garnet‐bearing assemblages from the Puga gneiss (host to the mafic eclogites) were utilized to constrain the subsequent decompression path. A non‐fractionated isochemical phase diagram for the assemblage phengite–garnet–biotite–plagioclase–quartz–melt documents a restricted (M3) P–T stability field centred on 12.5 ± 0.5 kbar and 690 ± 25 °C. A second non‐fractionated isochemical phase diagram calculated for the lower pressure assemblage garnet–cordierite–sillimanite–biotite–plagioclase–quartz–melt (M4) documents a narrow P–T stability field ranging between 7–8.4 kbar and 705–755 °C, which is consistent with independent multiequilibria P–T determinations. Th–Pb SHRIMP dating of monazite cores surrounded by allanite rims is interpreted to constrain the timing of the M4 equilibration to 45.3 ± 1.1 Ma. Coherently linking metamorphic conditions with petrographically constrained ages at Tso Morari provides an integrated context within which previously published petrological or geochronological results can be evaluated. The new composite path is similar to those published for the Kaghan UHP locality in northern Pakistan, although the calculated 12‐mm a^?1 rate of post‐pressure peak decompression at Tso Morari would appear less extreme.     

A new thermodynamic model for sapphirine: calculated phase equilibria in K2O–FeO–MgO–Al2O3–SiO2–H2O–TiO2–Fe2O3

The activity–composition (a–x) relations of sapphirine are re‐evaluated in the light of a recent new internally‐consistent data set of phase end‐members for use in phase equilibria modelling, particularly of ultra‐high‐temperature (UHT) rocks. This is achieved with the aid of relatively oxidized sapphirine+quartz‐bearing granulites from Wilson Lake, Canada. Calculated P–T projections and compatibility diagrams in the K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–Fe₂O₃ (KFMASHTO) system are used to illustrate sapphirine+quartz‐bearing phase equilibria in the context of UHT metamorphism. These new a–x relations for sapphirine should allow pseudosection thermobarometry in NCKFMASHTO for estimating peak P–T conditions of sapphirine‐bearing rocks.     

Closed system behaviour of argon in osumilite records protracted high‐T metamorphism within the Rogaland–Vest Agder Sector,Norway

Here, we present results of the first ⁴⁰Ar/³⁹Ar dating of osumilite, a high‐T mineral that occurs in some volcanic and high‐grade metamorphic rocks. The metamorphic osumilite studied here is from a metapelitic rock within the Rogaland–Vest Agder Sector, Norway, an area that experienced regional granulite facies metamorphism and subsequent contact metamorphism between 1,100 Ma and 850 Ma. The large grain size (~1 cm) of osumilite in the studied rock, which preserves a nominally anhydrous assemblage, increases the potential for large portions of individual grains to have remained essentially unaffected by the effects of diffusive argon loss, potentially preserving prograde ages. Step‐heating diffusion experiments yielded a maximum activation energy of ~461 kJ/mol and a pre‐exponential factor of ~8.34 × 10⁸ cm²/s for Ar diffusion in osumilite. These parameters correspond to a relatively high closure temperature of ~620°C for a cooling rate of 10°C/Ma in an osumilite crystal with a 175 µm radius. Fragments of osumilite separated from the sample preserve a range of ages between c. 1,070 and 860 Ma. The oldest ages are inferred to date the growth of coarse‐grained osumilite during prograde granulite facies regional metamorphism, which pre‐date contact metamorphism that has historically been ascribed to the growth of osumilite in the region. The majority of fragments record ages between c. 920 and 860 Ma, inferred to reflect the growth of osumilite and/or diffusive argon loss during contact metamorphism. The retention of old ⁴⁰Ar/³⁹Ar dates was facilitated by the low diffusivity of Ar in osumilite (i.e. a closed system), large grain sizes, and anhydrous metamorphic conditions. The ability to date osumilite with the ⁴⁰Ar/³⁹Ar method provides a valuable new thermochronometer that may constrain the timing and duration of high‐T magmatic and metamorphic events.     

Shooting at a moving target: phase equilibria modelling of high‐temperature metamorphism

Thermodynamic modelling and calculation of P–T pseudosections are commonly employed for quantifying the P–T evolution of metamorphic rocks. A key assumption involved in interpreting a P–T pseudosection is that the bulk‐rock composition used is representative of the effective bulk composition (EBC) from which apparently equilibrated mineral assemblages grew. Choosing an EBC can be difficult in cases where the rock has evolved significantly throughout the P–T history and has become domainal for whatever reason (e.g. loss of fluid and/or melt), particularly at suprasolidus conditions. During partial melting, melt migration may not only change the bulk composition by melt loss but also may generate local variations due to the variable consumption/loss of melt from domain to domain to create volumes of rock that were once internally equilibrated in the presence of a grain boundary melt, but which departed from equilibrium as inter‐granular mobility was slowed by local reductions in melt volume. As well as careful consideration of an EBC, the results of thermodynamic modelling are highly dependent on the specific thermodynamic data set and solution models used, as updates to these data sets may lead to substantially different calculated phase equilibria. This contribution addresses: (1) how consideration of evolving EBCs at multiple scales of observation can be used to resolve the history of complex high‐grade rocks, and (2) how use of different thermodynamic data sets and a–x models (i.e. thermocalc ds5.5 v. ds6) can result in different interpretations of metamorphic evolution. This study investigates the evolution of a mineralogically heterogeneous and texturally complex hand sample of granulite from the Gruf Complex (Central Alps). At the hand‐specimen scale, an EBC can be identified and used to constrain the P–T conditions at which the ‘whole rock’ was last in mutual equilibrium, in the presence of intergranular melt that has subsequently been lost or consumed. Smaller macrodomains (~cm scale) and microdomains (~mm scale) can be identified that represent subsequent evolution during and after melt channelization and loss, and P–T pseudosections can be calculated for the compositions of these domains. Using this approach reveals that the sample experienced a clockwise P–T path marked by near‐isothermal decompression following attainment of peak UHT conditions (~960 °C, 8.5 kbar). The approach enables construction of a P–T history of a rock for which P–T pseudosections are otherwise difficult to interpret. Thermodynamic modelling using ds6 yields similar results to those stated above, but suggests: (1) near‐isothermal decompression occurred over a wider pressure range (~0.5 kbar v. 1.5 kbar), and (2) that not all microdomains record this part of the P–T evolution.     

Ultrahigh‐temperature metamorphism in the Tuguiwula area,Khondalite Belt,North China Craton

In this study, sapphirine‐bearing granulites and sapphirine‐absent garnet–sillimanite gneisses from the Tuguiwula area in the eastern segment of the Khondalite Belt, North China Craton (NCC) are interpreted to show a P–T evolution involving cooling at pressures of 8–9 kbar from >960°C to the solidus (~820°C) and late subsolidus decompression. This interpretation is based on the sequence of mineral appearance and thermodynamic modelling of phase equilibria. Sapphirine is observed to coexist with spinel within the peak assemblages. This observation conflicts with the traditional view that spinel generally appears prior to sapphirine and thus indicates pre‐T_max compression. For ultrahigh‐temperature (UHT) metapelites at Tuguiwula, a clockwise P–T path may be more likely, which would be consistent with the clockwise P–T evolution of the extensive “normal” granulites (T_max <900°C) and UHT granulites at other localities in the eastern segment of the Khondalite Belt. At Tuguiwula, for UHT metapelites with low bulk‐rock Mg/(Mg+Fe^T), the oxidation state/Fe³⁺ content is interpreted to be a significant factor in controlling the mineral assemblages. We find that these compositions tend to contain sapphirine under oxidized conditions but spinel (without sapphirine) under reduced conditions. This difference may account for the simultaneous presence of both sapphirine‐bearing UHT granulites and sapphirine‐absent garnet–sillimanite UHT gneisses at Tuguiwula. LA‐ICP‐MS U–Pb dating of metamorphic zircon in the UHT metapelites yields mean ²⁰⁷Pb/²⁰⁶Pb ages of c. 1.92 Ga (two samples), which are interpreted to record the timing of cooling of the UHT rocks to the solidus. The UHT metamorphism is interpreted to have been generated by mantle upwelling and emplacement of mafic magmas within a post‐orogenic setting.     

Age and duration of ultrahigh-temperature metamorphism in the Anápolis–Itauçu Complex, Southern Brasília Belt, central Brazil – constraints from U–Pb geochronology, mineral rare earth element chemistry and trace-element thermometry

Integrated petrographic and chemical analysis of zircon, garnet and rutile from ultrahigh‐temperature (UHT) granulites in the Anápolis–Itauçu Complex, Brazil, is used to constrain the significance of zircon ID‐TIMS U–Pb geochronological data. Chondrite‐normalized rare earth element (REE) profiles of zircon cores have positive‐sloping heavy‐REE patterns, commonly inferred to be magmatic, whereas unambiguous metamorphic grains and overgrowths have flat to slightly negatively sloping heavy‐REE patterns. However, in one sample, a core of zircon interpreted as having formed prior to garnet crystallization and a metamorphic zircon formed within microstructures involving garnet breakdown both display elevated heavy‐REE (and Y) with positive‐sloping patterns. D_REE(zrc/grt) partition coefficients suggest an approximation to equilibrium between zircon and garnet cores, although progressive enrichment in heavy REE towards garnet rims occurs in two of the samples investigated. Titanium‐in‐zircon thermometry indicates zircon growth during both the prograde and post‐peak evolution, but not at peak temperatures of the UHT metamorphism. By contrast, zirconium‐in‐rutile thermometry of inclusions armoured by garnet records crystallization temperatures, based on the upper end of the interquartile range of the data, of ～890 to 870 °C and maximum temperature around 980 °C, indicating prograde and retrograde growth close to and after peak conditions. Rutile located in retrograde microstructures records crystallization temperatures of 850 to 820 °C. Rutile intergrown with ilmenite and included within orthopyroxene, which is associated with exsolved zircon, records temperatures ～760 °C, consistent with Ti‐in‐zircon crystallization temperatures. ID‐TIMS U–Pb geochronological data from two of the four samples investigated define upper intercept ages of 641.3 ± 8.4 Ma (MSWD 0.91) and 638.8 ± 2.5 Ma (MSWD 1.03) that correlate with periods of zircon growth along the prograde segment of the P–T path. Individual zircon U–Pb dates retrieved from all samples range from 649 to 634 Ma, indicating a maximum duration of c. 15 Myr for the UHT event. This period is interpreted as recording modest thickening of hot backarc lithosphere located behind the Arenópolis Arc at the edge of the São Francisco Craton consequent upon terminal collision of the Parána Block with the arc during the amalgamation of West Gondwana.     

A curious case of agreement between conventional thermobarometry and phase equilibria modelling in granulites: New constraints on P–T estimates in the Antarctica segment of the Musgrave–Albany–Fraser–Wilkes Orogen

The Windmill Islands region in Wilkes Land, east Antarctica, preserves granulite facies metamorphic mineral assemblages that yield seemingly comparable P–T estimates from conventional thermobarometry and mineral equilibria modelling. This is uncommon in granulite facies terranes, where conventional thermobarometry and phase equilibria modelling generally produce conflicting P–T estimates because peak mineral compositions tend to be modified by retrograde diffusion processes. In situ U–Pb monazite geochronology and calculated metamorphic phase diagrams show that the Windmill Islands experienced two phases of high thermal gradient metamorphism during the Mesoproterozoic. The first phase of metamorphism is recorded by monazite ages in two widely separated samples and occurred at c. 1,305 Ma. This event was regional in extent, involved crustally derived magmatism and reached conditions of ~3.2–5 kbar and 690–770°C corresponding to very high thermal gradients of >150°C/kbar. The elevated thermal regime is interpreted to reflect a period of extension or increased extension in a back‐arc setting that existed prior to c. 1,330 Ma. The first metamorphic event was overprinted by granulite facies metamorphism at c. 1,180 Ma that was coeval with the intrusion of charnockite. This event involved peak temperatures of ~840–850°C and pressures of ~4–5 kbar. A phase of granitic magmatism at c. 1,250–1,210 Ma, prior to the intrusion of the charnockite, is interpreted to reflect a phase of compression within an overall back‐arc setting. Existing conventional thermobarometry suggests conditions of ~4 kbar and 750°C for M₁ and 4–7 kbar and 750–900°C for M₂. The apparent similarities between the phase equilibria modelling and existing conventional thermobarometry may suggest either that the terrane cooled relatively quickly, or that the P–T ranges obtained from conventional thermobarometry are sufficiently imprecise that they cover the range of P–T conditions obtained in this study. However, without phase equilibria modelling, the veracity of existing conventional P–T estimates cannot be evaluated. The calculated phase diagrams from this study allow the direct comparison of P–T conditions in the Windmill Islands with phase equilibria models from other regions in the Musgrave–Albany–Fraser–Wilkes Orogen. This shows that the metamorphic evolution of the Wilkes Land region is very similar to that of the eastern Albany–Fraser Orogen and Musgrave Province in Australia, and further demonstrates the remarkable consistency in the timing of metamorphism and the thermal gradients along the ~5,000 km strike length of this system.     

Long‐lived high‐T,low‐P granulite facies metamorphism in the Arunta Region,central Australia

In situ LA–ICP–MS U–Pb monazite geochronology from the Boothby Hills in the Aileron Province, central Australia, indicates that the region records more than 80 Ma of high‐T, low‐P (HTLP) anatectic conditions during the Early Mesoproterozoic. Monazite ages from granulite facies rocks and leucosomes span the interval 1576–1542 Ma. Pegmatites that overprint the regional gneissic fabric and are interpreted to record the last vestiges of melt crystallization give ages between 1523 and 1513 Ma. Calculated P–T pseudosections suggest peak metamorphic conditions in excess of 850 °C at 0.65–0.75 GPa. The retrograde evolution was characterized by a P–T path that involved minor decompression and then cooling, culminating with the development of andalusite. Integration of the geochronological data set with the inferred P–T path trajectory suggests that suprasolidus cooling must have been slow, in the order of 2.5–4 °C Ma^?1. In addition, the retrograde P–T path trajectory suggests that HTLP conditions were generated within crust of relatively normal thickness. Despite the long duration over which anatectic conditions occurred, there is no evidence for external magmatic inputs or evidence that HTLP conditions were associated with long‐lived extension. Instead, it seems probable that the long‐lived HTLP metamorphism was driven to a significant extent by long‐lived conductive heating provided by high crustal heat production in voluminous pre‐metamorphic granitic rocks.     

P–T evolution of Precambrian eclogite in the Sveconorwegian orogen,SW Sweden

Conditions of the prograde, peak‐pressure and part of the decompressional P–T path of two Precambrian eclogites in the eastern Sveconorwegian orogen have been determined using the pseudosection approach. Cores of garnet from a Fe–Ti‐rich eclogite record a first prograde and syn‐deformational stage along a Barrovian gradient from ~670 °C and 7 kbar to 710 °C and 8.5 kbar. Garnet rims grew during further burial to 16.5–19 kbar at ~850–900 °C, along a steep dP/dT gradient. The pseudosection model of a kyanite‐bearing eclogite sample of more magnesian bulk composition confirms the peak conditions. Matrix reequilibration associated with subsequent near‐isothermal decompression and partial exhumation produced plagioclase‐bearing symplectites replacing kyanite and clinopyroxene at an estimated 850–870 °C and 10–11 kbar. The validity of the pseudosections is discussed in detail. It is shown that in pseudosection modelling the fractionation of FeO in accessory sulphides may cause a significant shift of field boundaries (here displaced by up to 1.5 kbar and 70 °C) and must not be neglected. Fast burial, exhumation and subsequent cooling are supported by the steepness of both the prograde and the decompressional P–T paths as well as the preservation of garnet growth zoning and the symplectitic reaction textures. These features are compatible with deep tectonic burial of the eclogite‐bearing continental crust as part of the underthrusting plate (Eastern Segment, continent Baltica) in a collisional setting that led to an effectively doubled crustal thickness and subsequent exhumation of the eclogites through tectonic extrusion. Our results are in accordance with regional structural and petrological relationships, which demonstrate foreland‐vergent partial exhumation of the eclogite‐bearing nappe along a basal thrust zone and support a major collisional stage at c. 1 Ga. We argue that the similarities between Sveconorwegian and Himalayan eclogite occurrences emphasize the modern style of Grenvillian‐aged tectonics.     

New P–T constraints on the Tamayen glaucophane‐bearing rocks,eastern Taiwan: Perple_X modelling results and geodynamic implications

New pseudosection modelling was applied to better constrain the P–T conditions and evolution of glaucophane‐bearing rocks in the Tamayen block of the Yuli belt, recognized as the world's youngest known blueschist complex. Based on the predominant clinoamphibole, textural relationships, and mineral compositions, these glaucophane‐bearing high‐P rocks can be divided into four types. We focused on the three containing garnet. The chief phase assemblages are (in decreasing mode): amphibole + quartz + epidote + garnet + chlorite + rutile/titanite (Type‐I), phengite + amphibole + quartz + garnet + chlorite + epidote + titanite + biotite + magnetite (Type‐II), and amphibole + quartz + albite + epidote + garnet + rutile + hematite + titanite (Type‐III). Amphibole exhibits compositional zoning from core to rim as follows: glaucophane → pargasitic amphibole → actinolite (Type‐I), barroisite → Mg‐katophorite/taramite → Fe‐glaucophane (Type‐II), glaucophane → winchite (Type‐III). Using petrographic data, mineral compositions and Perple_X modelling (pseudosections and superimposed isopleths), peak P–T conditions were determined as 13 ± 1 kbar and 550 ± 40 °C for Type‐I, 10.5 ± 0.5 kbar and 560 ± 30 °C for Type‐II (thermal peak) and 11 ± 1 kbar and 530 ± 30 °C for Type‐III. The calculations yield higher pressures and temperatures than previously thought; the difference is ~1–6 kbar and 50–200 °C. The three rock types record similar P–T retrograde paths with clockwise trajectories; all rocks followed trajectories with substantial pressure decrease under near‐isothermal conditions (Type‐I and Type‐III), with the probable exception of Type‐II where decompression followed colder geotherms. The P–T paths suggest a tectonic environment in which the rocks were exhumed from maximum depths of ~45 km within a subduction channel along a relative cold geothermal gradient of ~11–14 °C km^?1.     

Phase equilibria modelling constraints on P–T conditions during fluid catalysed conversion of granulite to eclogite in the Bergen Arcs,Norway

Exhumed eclogitic crust is rare and exposures that preserve both protoliths and altered domains are limited around the world. Nominally anhydrous Neoproterozoic anorthositic granulites exposed on the island of Holsnøy, in the Bergen Arcs in western Norway, preserve different stages of progressive prograde deformation, together with the corresponding fluid‐assisted metamorphism, which record the conversion to eclogite during the Ordovician–Silurian Caledonian Orogeny. Four stages of deformation can be identified: (1) brittle deformation resulting in the formation of fractures and the generation of pseudotachylites in the granulite; (2) development of mesoscale shear zones associated with increased fluid–rock interaction; (3) the complete large‐scale replacement of granulite by hydrous eclogite with blocks of granulite sitting in an eclogitic “matrix”; and (4) the break‐up of completely eclogitized granulite by continued fluid influx, resulting in the formation of coarse‐grained phengite‐dominated mineral assemblages. P–T constraints derived from phase equilibria forward modelling of mineral assemblages of the early and later stages of the conversion to eclogite document burial and partial exhumation path with peak metamorphic conditions of ~21–22 kbar and 670–690°C. The P–T models, in combination with existing T–t constraints, imply that the Lindås Nappe underwent extremely rapid retrogressive pressure change. Fluid infiltration began on the prograde burial path and continued throughout the recorded P–T evolution, implying a fluid source that underwent progressive dehydration during subduction of the granulites. However, in places limited fluid availability on the prograde path resulted in assemblages largely consuming the available fluid, essentially freezing in snapshots of the prograde evolution. These were carried metastably deeper into the mantle with strain and metamorphic recrystallization partitioned into areas where ongoing fluid infiltration was concentrated.     

Retrograde processes in migmatites and granulites revisited

Many migmatites and granulites preserve evidence of a clockwise P–T evolution involving decompression (decrease in P) while close to the thermal peak. The extent of post‐thermal peak reaction is influenced by several factors, including: (1) the P–T path in relation to invariants in the system and the Clapeyron slopes of the equilibria; (2) the rate of cooling; and (3) the availability of fluid (H₂O‐rich volatile phase or melt) for fluid‐consuming reactions. Reaction may occur between products of a prograde (increasing T) fluid‐generating reaction as the same equilibrium is re‐crossed in the retrograde (decreasing T) sense. In general, reaction reversal or ‘back reaction’ requires the P–T path to approximate isobaric heating and cooling, without significant decompression, and evolved fluid to remain within the equilibration volume. The larger the decompression segment in the P–T evolution, the more chance there is of crossing different reactions along the retrograde segment from those crossed along the prograde segment. For common pelite compositions, we may generalize by considering three pressure regimes separated by the [Spl, Ms, H₂O] invariant in KFMASH (approximately 9 kbar) and the intersection of muscovite breakdown with the H₂O‐rich volatile phase‐saturated solidus (approximately 4 kbar). Reaction reversal cannot occur along P–T paths that traverse around one of these points, but may occur along P–T paths confined to one of the three regimes in between. Additionally, above the solidus, melt segregation and loss potentially change the composition of the equilibration volume; and, the size of the equilibration volume shrinks with decreasing T. Since the proportion of melt to residue in the equilibration volume may change with decreasing size, the composition of the equilibration volume may change throughout the supra‐solidus part of the retrograde segment of the P–T evolution. If melt has been lost from the equilibration volume, reaction reversal may not be possible or may be only partial; indeed, the common preservation of close‐to‐peak mineral assemblages in migmatite and granulite demonstrates that extensive reaction with melt is uncommon, which implies melt isolation or loss prior to crossing potential melt‐consuming reactions. Water dissolved in melt is transported through the crust to be exsolved on crystallization at the solidus appropriate to the intrinsic a(H₂O). This recycled water causes retrogression at subsolidus conditions. Consideration of the evidence for supra‐solidus decompression‐dehydration reactions, and review of microstructures that have proven controversial, such as corona and related microstructures, selvage microstructures and ‘late’ muscovite, leads to the conclusion that there is more than one way for these microstructures to form and reminds us that we should always consider multiple working hypotheses!