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.     

Hot lherzolite exhumation,UHT migmatite formation,and acid volcanism driven by Miocene rollback of the Banda Arc,eastern Indonesia

The northern Banda Arc, eastern Indonesia, exposes upper mantle/lower crustal complexes comprising lherzolites and granulite facies migmatites of the ‘Kobipoto Complex’. Residual garnet–sillimanite granulites, which contain spinel + quartz inclusions within garnet, experienced ultrahigh-temperature (UHT; > 900 °C) conditions at 16 Ma due to heat supplied by lherzolites exhumed during slab rollback in the Banda Arc. Here, we present U–Pb zircon ages and new whole-rock geochemical analyses that document a protracted history of high-T metamorphism, melting, and acid magmatism of a common sedimentary protolith. Detrital zircons from the Kobipoto Complex migmatites, with ages between 3.4 Ga and 216 Ma, show that their protolith was derived from both West Papua and the Archean of Western Australia, and that metamorphism of these rocks on Seram could not have occurred until the Late Triassic. Zircons within the granulites then experienced three subsequent episodes of growth – at 215–173 Ma, 25–20 Ma, and at c. 16 Ma. The population of zircon rims with ages between 215 and 173 Ma document significant metamorphic (± partial melting) events that we attribute to subduction beneath the Bird's Head peninsula and Sula Spur, which occurred until the Banda and Argo continental blocks were rifted from the NW Australian margin of Gondwana in the Late Jurassic (from c. 160 Ma). Late Oligocene-Early Miocene collision between Australia (the Sula Spur) and SE Asia (northern Sulawesi) was then recorded by crystallisation of several 25–20 Ma zircon rims. Thereafter, a large population of c. 16 Ma zircon rims grew during subsequent and extensive Middle Miocene metamorphism and melting of the Kobipoto complex rocks beneath Seram under high- to ultrahigh-temperature (HT–UHT) conditions. Lherzolites located adjacent to the granulite-facies migmatites in central Seram equilibrated at 1280–1300 °C upon their exhumation to 1 GPa (~ 37 km) depth, whereupon they supplied sufficient heat to have metamorphosed adjacent Kobipoto Complex migmatites under UHT conditions at 16 Ma. Calculations suggesting slight (~ 10 vol%) mantle melting are consistent with observations of minor gabbroic intrusions and scarce harzburgites. Subsequent extension during continued slab rollback exhumed both the lherzolites and adjacent granulite-facies migmatites beneath extensional detachment faults in western Seram at 6.0–5.5 Ma, and on Ambon at 3.5 Ma, as recorded by subsequent zircon growth and ⁴⁰Ar/³⁹Ar ages in these regions. Ambonites, cordierite- and garnet-bearing dacites sourced predominantly from melts generated in the Kobipoto Complex migmatites, were later erupted on Ambon from 3.0 to 1.9 Ma.     

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.     

Ultra high temperature metamorphism of mafic granulites from South Altyn Orogen,West China: A result from the rapid exhumation of deeply subducted continental crust

The South Altyn orogen in West China contains ultra high pressure (UHP) terranes formed by ultra‐deep (>150–300 km) subduction of continental crust. Mafic granulites which together with ultramafic interlayers occur as blocks in massive felsic granulites in the Bashiwake UHP terrane, are mainly composed of garnet, clinopyroxene, plagioclase, amphibole, rutile/ilmenite, and quartz with or without kyanite and sapphirine. The kyanite/sapphirine‐bearing granulites are interpreted to have experienced decompression‐dominated evolution from eclogite facies conditions with peak pressures of 4–7 GPa to high pressure (HP)–ultra high temperature (UHT) granulite facies conditions and further to low pressure (LP)–UHT facies conditions based on petrographic observations, phase equilibria modelling, and thermobarometry. The HP–UHT granulite facies conditions are constrained to be 2.3–1.6 GPa/1,000–1,070°C based on the observed mineral assemblages of garnet+clinopyroxene+rutile+plagioclase+amphibole±quartz and measured mineral compositions including the core–rim increasing anorthite in plagioclase (X_An = 0.52–0.58), core–rim decreasing jadeite in clinopyroxene (X_Jd = 0.20–0.15), and TiO₂ in amphibole (Ti^M2/2 = 0.14–0.18). The LP–UHT granulite facies conditions are identified from the symplectites of sapphirine+plagioclase+spinel, formed by the metastable reaction between garnet and kyanite at <0.6–0.7 GPa/940–1,030°C based on the calculated stability of the symplectite assemblages and sapphirine–spinel thermometer results. The common granulites without kyanite/sapphirine are identified to record a similar decompression evolution, including eclogite, HP–UHT granulite, and LP–UHT granulite facies conditions, and a subsequent isobaric cooling stage. The decompression under HP–UHT granulite facies is estimated to be from 2.3 to 1.3 GPa at ~1,040°C on the basis of textural records, anorthite content in plagioclase (X_An = 0.25–0.32), and grossular content in garnet (X_Grs = 0.22–0.19). The further decompression to LP–UHT facies is defined to be >0.2–0.3 GPa based on the calculated stability for hematite‐bearing ilmenite. The isobaric cooling evolution is inferred mainly from the amphibole (Ti^M2/2 = 0.14–0.08) growth due to the crystallization of residual melts, consistent with a temperature decrease from >1,000°C to ~800°C at ~0.4 GPa. Zircon U–Pb dating for the two types of mafic granulite yields similar protolith and metamorphic ages of c. 900 Ma and c. 500 Ma respectively. However, the metamorphic age is interpreted to represent the HP–UHT granulite stage for the kyanite/sapphirine‐bearing granulites, but the isobaric cooling stage for the common granulites on the basis of phase equilibria modelling results. The two types of mafic granulite should share the same metamorphic evolution, but show contrasting features in petrography, details of metamorphic reactions in each stage, thermobarometric results, and also the meaning of zircon ages as a result of their different bulk‐rock compositions. Moreover, the UHT metamorphism in UHP terranes is revealed to represent the lower pressure overprinting over early UHP assemblages during the rapid exhumation of ultra‐deep subducted continental slabs, in contrast to the cause of traditional UHT metamorphism by voluminous heat addition from the mantle.     

Jadeite–garnet glaucophane schists in the Bizan area,Sambagawa metamorphic belt,eastern Shikoku,Japan: significance and extent of eclogite facies metamorphism

Eclogite facies metamorphic rocks have been discovered from the Bizan area of eastern Shikoku, Sambagawa metamorphic belt. The eclogitic jadeite–garnet glaucophane schists occur as lenticular or sheet‐like bodies in the pelitic schist matrix, with the peak mineral assemblage of garnet + glaucophane + jadeite + phengite + quartz. The jadeitic clinopyroxene (X_Jd 0.46–0.75) is found exclusively as inclusions in porphyroblastic garnet. The eclogite metamorphism is characterized by prograde development from epidote–blueschist to eclogite facies. Metamorphic P–T conditions estimated using pseudosection modelling are 580–600 °C and 18–20 kbar for eclogite facies. Compared with common mafic eclogites, the jadeite–garnet glaucophane schists have low CaO (4.4–4.5 wt%) and MgO (2.1–2.3 wt%) bulk‐rock compositions. The P–T– pseudosections show that low X_Ca bulk‐rock compositions favour the appearance of jadeite instead of omphacite under eclogite facies conditions. This is a unique example of low X_Ca bulk‐rock composition triggered to form jadeite at eclogite facies conditions. Two significant types of eclogitic metamorphism have been distinguished in the Sambagawa metamorphic belt, that is, a low‐T type and subsequent high‐T type eclogitic metamorphic events. The jadeite–garnet glaucophane schists experienced low‐T type eclogite facies metamorphism, and the P–T path is similar to lawsonite‐bearing eclogites recently reported from the Kotsu area in eastern Shikoku. During subduction of the oceanic plate (Izanagi plate), the hangingwall cooled gradually, and the geothermal gradient along the subduction zone progressively decreased and formed low‐T type eclogitic metamorphic rocks. A subsequent warm subduction event associated with an approaching spreading ridge caused the high‐T type eclogitic metamorphism within a single subduction zone.     

Paleoproterozoic UHT metamorphism with isobaric cooling (IBC) followed by decompression–heating in the Khondalite Belt (North China Craton): New evidence from two sapphirine formation processes

Interpretation of reaction microstructures may provide constraints on the P–T path followed by rocks, with implications for the geodynamic evolution. Sapphirine generally occurs in diverse microstructures in ultrahigh-temperature (UHT) Mg–Al-rich granulites. Understanding multi-stage sapphirine formation processes and the resultant P–T path may provide insights into the cause of UHT metamorphism, which is otherwise under broad debate. Here, we investigate samples of UHT granulite containing two distinct types of sapphirine from the Dongpo locality in the Khondalite Belt, North China Craton, with the aim of understanding the processes of sapphirine formation and the metamorphic evolution of the host rocks. Petrographic observations show that early sapphirine, which occurs as coronas on spinel and as single porphyroblasts, formed together with biotite, sillimanite, and inclusion-rich garnet. Late symplectitic sapphirine along with fine-grained plagioclase and spinel plus plagioclase symplectites, formed by consumption of sillimanite, biotite, and garnet. Three pseudosections based on the bulk compositions of microdomains inferred to reflect spatially restricted equilibrium suggest that the rocks record near isobaric cooling (IBC) from ~980 to 830ºC at ~0.9 GPa for early sapphirine formation, and decompression and heating to ≤0.7 GPa and ~900ºC for late sapphirine formation. Our study in combination with other metamorphic P–T and age information reveals the common occurrence of IBC paths and long duration (c. 1.93 to 1.86 Ga) regional UHT metamorphism in the Khondalite Belt, North China Craton. Locally, this is followed by decompression–heating paths at c. 1.86 Ga. The Palaeoproterozoic UHT metamorphism with long-lived IBC path in the Khondalite Belt, North China Craton supports large hot orogen model in the amalgamation of this part in the supercontinent Nuna.     

Phase equilibria modelling and zircon age dating of pelitic granulites in Zhaojiayao,from the Jining Group of the Khondalite Belt,North China Craton

The Jining Group occurs as the eastern segment of the Khondalite Belt, North China Craton and is dominated by a series of granulite facies rocks involving ‘normal’ pelitic granulites recording peak temperatures of ~850 °C and ultrahigh‐temperature (UHT) pelitic granulites recording peak temperatures of 950–1100 °C. The P–T paths and ages of these two types of granulites are controversial. Three pelitic granulite samples in the Jining Group comprising two sillimanite–garnet gneiss samples (J1208 and J1210) and one spinel–garnet gneiss sample (J1303) were collected from Zhaojiayao, where ‘normal’ pelitic granulites occur, for determination of their metamorphic evolution and ages. Samples J1208 and J1210 are interpreted to record cooling paths from the T_max stages with P–T conditions respectively of ~870–890 °C/7–8 kbar and >840 °C/>7.5 kbar constrained from the stability fields of the observed mineral assemblages and the isopleths of plagioclase, garnet and biotite compositions in pseudosections. Sample J1303 is interpreted to record pre‐T_max decompression from the kyanite‐stability fields to the T_max stage of 950–1020 °C/8–9 kbar and a post‐T_max cooling path revealed mainly from the stability field of the observed mineral assemblage, the plagioclase zoning and the biotite composition isopleth in pseudosections. The post‐T_max cooling stage can be divided into suprasolidus and subsolidus stages. The suprasolidus cooling may not result in an equilibrium state at the solidus in a rock. Therefore, different minerals may record different P–T conditions along the cooling path; the inferred maximum temperature is commonly higher than the solidus as well as different solidi being recorded for different samples from the same outcrop but experiencing different degrees of melt loss. Plagioclase compositions, especially its zoning in plagioclase‐rich granulites, are predicted to be useful for recording the higher temperature conditions of a granulite's thermal history. The three samples studied seem to record the temperature range covering those of the ‘normal’ and UHT pelitic granulites in the Jining Group, suggesting that UHT conditions may be reached in ‘normal’ granulites without diagnostic UHT indicators. LA‐ICP‐MS zircon U–Pb data provide a continuous trend of concordant ²⁰⁷Pb/²⁰⁶Pb ages from 1.89 to 1.79 Ga for sample J1210, and from 1.94 to 1.80 Ga for sample J1303. These continuous and long age spectrums are interpreted to represent a slow cooling and exhumation process corresponding to the post‐T_max cooling P–T paths recorded by the pelitic granulites, which may have followed the exhumation of deeply buried rocks in a thickened crust region resulted from a collision event at c. 1.95 Ga as suggested by the pre‐T_max decompression P–T path.     

Constraints from monazite and xenotime growth modelling in the MnCKFMASH‐PYCe system on the P–T path of a metapelite from Shillong‐Meghalaya Plateau: implications for the Indian shield assembly

The Palaeo‐Mesoproterozoic metapelite granulites from northern Garo Hills, western Shillong‐Meghalaya Gneissic Complex (SMGC), northeast India, consist of resorbed garnet, cordierite and K‐feldspar porphyroblasts in a matrix comprising shape‐preferred aggregates of biotite±sillimanite+quartz that define the penetrative gneissic fabric. An earlier assemblage including biotite and sillimanite occurs as inclusions within the garnet and cordierite porphyroblasts. Staurolite within cordierite in samples without matrix sillimanite is interpreted to have formed by a reaction between the sillimanite inclusion and the host cordierite during retrogression. Accessory monazite occurs as inclusions within garnet as well as in the matrix, whereas accessory xenotime occurs only in the matrix. The monazite inclusions in garnet contain higher Ca, and lower Y and Th/U than the matrix monazite outside resorbed garnet rims. On the other hand, matrix monazite away from garnet contains low Ca and Y, and shows very high Th/U ratios. The low Th/U ratios (<10) of the Y‐poor garnet‐hosted monazite indicate subsolidus formation during an early stage of prograde metamorphism. A calculated P–T pseudosection in the MnCKFMASH‐PYCe system indicates that the garnet‐hosted monazite formed at <3 kbar/600 °C (Stage A). These P–T estimates extend backward the previously inferred prograde P–T path from peak anatectic conditions of 7–8 kbar/850 °C based on major mineral equilibria. Furthermore, the calculated P–T pseudosections indicate that cordierite–staurolite equilibrated at ~5.5 kbar/630 °C during retrograde metamorphism. Thus, the P–T path was counterclockwise. The Y‐rich matrix monazite outside garnet rims formed between ~3.2 kbar/650 °C and ~5 kbar/775 °C (Stage B) during prograde metamorphism. If the effect of bulk composition change due to open system behaviour during anatexis is considered, the P–T conditions may be lower for Stage A (<2 kbar/525 °C) and Stage B (~3 kbar/600 °C to ~3.5 kbar/660 °C). Prograde garnet growth occurred over the entire temperature range (550–850 °C), and Stage‐B monazite was perhaps initially entrapped in garnet. During post‐peak cooling, the Stage‐B monazite grains were released in the matrix by garnet dissolution. Furthermore, new matrix monazite (low Y and very high Th/U ≤80, ~8 kbar/850–800 °C, Stage C), some monazite outside garnet rims (high Y and intermediate Th/U ≤30, ~8 kbar/800–785 °C, Stage D), and matrix xenotime (<785 °C) formed through post‐peak crystallization of melt. Regardless of textural setting, all monazite populations show identical chemical ages (1630–1578 Ma, ±43 Ma). The lithological association (metapelite and mafic granulites), and metamorphic age and P–T path of the northern Garo Hills metapelites and those from the southern domain of the Central Indian Tectonic Zone (CITZ) are similar. The SMGC was initially aligned with the southern parts of CITZ and Chotanagpur Gneissic Complex of central/eastern India in an ENE direction, but was displaced ~350 km northward by sinistral movement along the north‐trending Eastern Indian Tectonic Zone in Neoproterozoic. The southern CITZ metapelites supposedly originated in a back‐arc associated with subducting oceanic lithosphere below the Southern Indian Block at c. 1.6 Ga during the initial stage of Indian shield assembly. It is inferred that the SMGC metapelites may also have originated contemporaneously with the southern CITZ metapelites in a similar back‐arc setting.     

Modelling of mineral equilibria in ultrahigh-temperature metamorphic rocks from the Anápolis–Itauçu Complex, central Brazil

A new quantitative approach to constraining mineral equilibria in sapphirine‐bearing ultrahigh‐temperature (UHT) granulites through the use of pseudosections and compatibility diagrams is presented, using a recently published thermodynamic model for sapphirine. The approach is illustrated with an example from an UHT locality in the Anápolis–Itauçu Complex, central Brazil, where modelling of mineral equilibria indicates peak metamorphic conditions of about 9 kbar and 1000 °C. The early formed, coarse‐grained assemblage is garnet–orthopyroxene–sillimanite–quartz, which was subsequently modified following peak conditions. The retrograde pressure–temperature (P–T) path of this locality involves decompression across the FeO–MgO–Al₂O₃–SiO₂ (FMAS) univariant reaction orthopyroxene + sillimanite = garnet + sapphirine + quartz, resulting in the growth of sapphirine–quartz, followed by cooling and recrossing of this reaction. The resulting microstructures are modelled using compatibility diagrams, and pseudosections calculated for specific grain boundaries considered as chemical domains. The sequence of microstructures preserved in the rocks constrains a two‐stage isothermal decompression–isobaric cooling path. The stability of cordierite along the retrograde path is examined using a domainal approach and pseudosections for orthopyroxene–quartz and garnet–quartz grain boundaries. This analysis indicates that the presence or absence of cordierite may be explained by local variation in a_H2O. This study has important implications for thermobarometric studies of UHT granulites, mainly through showing that traditional FMAS petrogenetic grids based on experiments alone may overestimate P–T conditions. Such grids are effectively constant a_H2O sections in FMAS‐H₂O (FMASH), for which the corresponding a_H2O is commonly higher than that experienced by UHT granulites. A corollary of this dependence of mineral equilibria on a_H2O is that local variations in a_H2O may explain the formation of cordierite without significant changes in P–T conditions, particularly without marked decompression.     

The interpretation of reaction textures in Fe‐rich metapelitic granulites of the Musgrave Block,central Australia: constraints from mineral equilibria calculations in the system K2O–FeO–MgO–Al2O3–SiO2–H2O–TiO2–Fe2O3

Fe‐rich metapelitic granulites of the Musgrave Block, central Australia, contain several symplectic and coronal reaction textures that post‐date a peak S2 metamorphic assemblage involving garnet, sillimanite, spinel, ilmenite, K‐feldspar and quartz. The earliest reaction textures involve spinel‐ and quartz‐bearing symplectites that enclose garnet and to a lesser extent sillimanite. The symplectic spinel and quartz are in places separated by later garnet and/or sillimanite coronas. The metamorphic effects of a later, D3, event are restricted to zones of moderate to high strain where a metamorphic assemblage of garnet, sillimanite, K‐feldspar, magnetite, ilmenite, quartz and biotite is preserved. Quantitative mineral equilibria calculations in the system K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–Fe₂O₃ (KFMASHTO) using Thermocalc 3.0 and the accompanying internally consistent dataset provide important constraints on the influence of TiO₂ and Fe₂O₃ on biotite‐bearing and spinel‐bearing equilibria, respectively. Biotite‐bearing equilibria are shifted to higher temperatures and spinel‐bearing equilibria to higher pressures and lower temperatures in comparison to the equivalent equilibria in K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O (KFMASH). The sequence of reaction textures involving spinel is consistent with a D2 P–T path that involved a small amount of decompression followed predominantly by cooling within a single mineral assemblage stability field. Thus, the reaction textures reflect changes in modal proportions within an equilibrium assemblage rather than the crossing of a univariant reaction. The D3 metamorphic assemblage is consistent with lower temperatures than those inferred for D2.     

A unique Paleoproterozoic HP–UHT metamorphic event recorded by the Bengbu mafic granulites in the southwestern Jiao–Liao–Ji Belt,North China Craton

The Wuhe Complex in the Bengbu area of the Jiao–Liao–Ji Belt, southeast North China Craton, contains garnet-bearing mafic granulites that have undergone high-pressure (HP) and ultrahigh-temperature (UHT) metamorphism. These granulites also experienced partial melting and occur as lenses within marbles. Petrographic observations and quantitative phase equilibria modeling reveal clockwise P–T paths, involving an inferred HP stage followed by decompressional, medium-pressure, granulite-facies metamorphism and subsequent cooling. The HP assemblage of garnet + clinopyroxene + plagioclase + K-feldspar ± amphibole ± quartz ± rutile indicates P–T conditions of 840–980 °C and 12–17 kbar. This was followed by post-peak, near-isothermal decompression with the development of orthopyroxene + clinopyroxene + plagioclase + K-feldspar + garnet + amphibole + ilmenite at 850–960 °C and 7–10 kbar, resulting in the development of orthopyroxene rims on resorbed garnet. Pyroxene and ternary feldspar thermometry yielded high temperatures of ~1150 °C and 1055–1087 °C at 10 kbar, respectively, which constrain the minimum crystallization temperatures of the igneous protoliths. The host and lamellae of the pyroxene and ternary feldspar are relict magmatic minerals/textures that survived metamorphism due to the silica-undersaturated bulk-rock conditions. Zr-in-rutile thermometry yielded temperatures of ~935 °C and 800 °C, with the former being consistent with the predicted peak metamorphic temperatures. Small amounts of melts (up to 5%) were generated during decompression of the Bengbu mafic granulites. The generated partial melts were mainly (quartz) monzonite at 900–920 °C, and the silica contents of the melts were controlled by the quartz stability field in P–T pseudosections. The partial melts were enriched in Na and strongly depleted in Fe–Mg at the peak pressure of ~14 kbar and 920 °C, and later evolved to Fe–Mg-rich and high-K compositions during decompression. The melt compositions in the studied rocks are similar when the pressures reached ~9 kbar. The modal proportion of amphibole increased as the melt H₂O content decreased at lower pressures, indicating that the limited H₂O remaining in the host rocks was consumed to produce amphibole. U–Pb geochronology of zircon containing inclusions of clinopyroxene, plagioclase, and apatite constrains the timing of metamorphism to 1930–1840 Ma, as is the case for HP granulites from Shandong, Liaoning, and southern Jilin in the central and northeastern Jiao–Liao–Ji Belt. The Wuhe HP–UHT mafic granulites were ultimately sourced from upwelling asthenosphere-derived magma at ~2.1 Ga, which intruded and crystallized at shallower depths. The igneous protoliths were then buried to middle–lower crustal levels and experienced HP–UHT granulite-facies metamorphism and partial melting at 1.95–1.90 Ga related to continental subduction and overthickening. The HP–UHT mafic granulites were rapidly exhumed at ~1.85 Ga and generated small volumes of (quartz) monzonite during decompression. The newly discovered Paleoproterozoic HP–UHT mafic granulites associated with partial melting suggest that the continent materials were deeply subducted to the lower crustal levels and that additional heating was not involved. The finding of the HP–UHT granulites, together with the widespread distributions of the granulite-facies metamorphic rocks and the determination of the clockwise P–T–t paths, reveal that the Paleoproterozoic Jiao–Liao–Ji orogenic belt extends at least 1000 km, starting from southern Jilin, passing through the southeastern Liaoning and Jiaobei terranes, and elongating to the Bengbu area in Anhui.     

Late‐stage orogenic loading revealed by contact metamorphism in the northern Appalachians,New York

Pelitic schists from contact aureoles surrounding mafic–ultramafic plutons in Westchester County, NY record a high‐P (~0.8 GPa) high‐T (~790 °C) contact overprint on a Taconic regional metamorphic assemblage (~0.5 GPa). The contact metamorphic assemblage of a pelitic sample in the innermost aureole of the Croton Falls pluton, a small (<10 km²) gabbroic body, consists of quartz–plagioclase–biotite–garnet–sillimanite–ilmenite–graphite–Zn‐rich Al‐spinel. Both K‐feldspar and muscovite are absent, and abundant biotite, plagioclase, sillimanite, quartz and ilmenite inclusions are found within subhedral garnet crystals. Unusually low bulk‐rock Na and K contents imply depletion of alkalic components and silica through anatexis and melt extraction during contact heating relative to typical metapelites outside the aureole. Thermobarometry on nearby samples lacking a contact overprint yields 620–640 °C and 0.5–0.6 GPa. In the aureole sample, WDS X‐ray chemical maps show distinct Ca‐enriched rims on both garnet and matrix plagioclase. Furthermore, biotite inclusions within garnet have significantly higher Mg concentration than matrix biotite. Thermobarometry using GASP and garnet–biotite Mg–Fe exchange equilibria on inclusions and adjacent garnet host interior to the high‐Ca rim zone yield ~0.5 ± 0.1 GPa and ~620 ± 50 °C. Pairs in the modified garnet rim zone yield ~0.9 ± 0.1 GPa and ~790 ± 50 °C. Thermocalc average P–T calculations yield similar results for core (~0.5 ± ~0.1 GPa, ~640 ± ~80 °C) and rim (~0.9 ± ~0.1 GPa, ~800 ± ~90 °C) equilibria. The core assemblages are interpreted to record the P–T conditions of peak metamorphism during the Taconic regional event whereas the rim compositions and matrix assemblages are interpreted to record the P–T conditions during the contact event. The high pressures deduced for this later event are interpreted to reflect loading due to the emplacement of Taconic allochthons in the northern Appalachians during the waning stages of regional metamorphism (after c. 465 Ma) and before contact metamorphism (c. 435 Ma). In the absence of contact metamorphism‐induced recrystallization, it is likely that this regional‐scale loading would remain cryptic or unrecorded.     

Phase equilibria modelling and LASS monazite petrochronology: P–T–t constraints on the evolution of the Priest River core complex,northern Idaho

Phase equilibria modelling, laser‐ablation split‐stream (LASS)‐ICP‐MS petrochronology and garnet trace‐element geochemistry are integrated to constrain the P–T–t history of the footwall of the Priest River metamorphic core complex, northern Idaho. Metapelitic, migmatitic gneisses of the Hauser Lake Gneiss contain the peak assemblage garnet + sillimanite + biotite ± muscovite + plagioclase + K‐feldspar ± rutile ± ilmenite + quartz. Interpreted P–T paths predict maximum pressures and peak metamorphic temperatures of ~9.6–10.3 kbar and ~785–790 °C. Monazite and xenotime ²⁰⁸Pb/²³²Th dates from porphyroblast inclusions indicate that metamorphism occurred at c. 74–54 Ma. Dates from HREE‐depleted monazite formed during prograde growth constrain peak metamorphism at c. 64 Ma near the centre of the complex, while dates from HREE‐enriched monazite constrain the timing of garnet breakdown during near‐isothermal decompression at c. 60–57 Ma. Near‐isothermal decompression to ~5.0–4.4 kbar was followed by cooling and further decompression. The youngest, HREE‐enriched monazite records leucosome crystallization at mid‐crustal levels c. 54–44 Ma. The northernmost sample records regional metamorphism during the emplacement of the Selkirk igneous complex (c. 94–81 Ma), Cretaceous–Tertiary metamorphism and limited Eocene exhumation. Similarities between the Priest River complex and other complexes of the northern North American Cordillera suggest shared regional metamorphic and exhumation histories; however, in contrast to complexes to the north, the Priest River contains less partial melt and no evidence for diapiric exhumation. Improved constraints on metamorphism, deformation, anatexis and exhumation provide greater insight into the initiation and evolution of metamorphic core complexes in the northern Cordillera, and in similar tectonic settings elsewhere.     

The significance of Mn‐rich ilmenite and the determination of P–T paths from zoned garnet in metasedimentary rocks from the western Cape Breton Highlands,Nova Scotia

The Jumping Brook Metamorphic Suite in the western Cape Breton Highlands of Nova Scotia is part of an inverted Barrovian sequence that formed during a Late Silurian–Early Devonian promontory–promontory collision in the Canadian Appalachians. In this study, systematic discrepancies between geochemical observations and thermodynamic model predictions led to the discovery of a systematic relationship linking the style of garnet core isopleth intersection (GCII) to the pyrophanite (MnTiO₃) component of co‐existing ilmenite. Samples that yielded tight GCIIs at or near the garnet‐in curve were found to contain ilmenite with negligible pyrophanite components, whereas samples yielding GCIIs far removed (up to 105°C) from the garnet‐in curve were found to contain ilmenite with significant pyrophanite and/or ecandrewsite (ZnTiO₃) components. Based on petrographic and geochemical observations, Mn(±Zn)‐rich ilmenite are interpreted to have sequestered Mn throughout prograde metamorphism due to sluggish intracrystalline diffusion. The amount of reactive Mn input into the thermodynamic models from whole‐rock analyses were, in some cases, overestimated, resulting in garnet‐in curve topologies that extend to erroneously low P–T conditions. Modifications to the whole‐rock chemistry that account for Mn sequestration into ilmenite, however, yielded robust model results. Our results show that, in addition to uncertainties in thermodynamic data sets and phenomenon related to reaction kinetics, Mn‐rich ilmenite may superimpose additional complexities related to the interpretation of predicted equilibria involving garnet. Numerical simulations of garnet crystallization were used to infer P–T paths of metamorphism for one sample from the garnet zone (Mn corrected) and two samples from the staurolite zone (Mn uncorrected) of the inverted sequence. Model results are remarkably similar among the three samples and indicate that garnet crystallization occurred along relatively steep (31–37°C/km) clockwise P–T paths. The peak conditions of garnet crystallization and metamorphism (560–590°C, 7.4–8.0 kbar) are interpreted to have been attained approximately simultaneously, such that the paths are characterized by tight prograde‐to‐retrograde transitions. The hairpin nature of the P–T paths is interpreted to represent the onset of thrust‐related exhumation and isograd inversion along ductile shear zones, consistent with available field and geochronological constraints.     

Thermal evolution of the Tseel terrane,SW Mongolia and its relation to granitoid intrusions in the Central Asian Orogenic Belt

The timing and thermal effects of granitoid intrusions into accreted sedimentary rocks are important for understanding the growth process of continental crust. In this study, the petrology and geochronology of pelitic gneisses in the Tseel area of the Tseel terrane, SW Mongolia, are examined to understand the relationship between igneous activity and metamorphism during crustal evolution in the Central Asian Orogenic Belt (CAOB). Four mineral zones are recognized on the basis of progressive changes in the mineral assemblages in the pelitic gneisses, namely: the garnet, staurolite, sillimanite and cordierite zones. The gneisses with high metamorphic grades (i.e. sillimanite and cordierite zones) occur in the central part of the Tseel area, where granitoids are abundant. To the north and south of these granitoids, the metamorphic grade shows a gradual decrease. The composition of garnet in the pelitic gneisses varies systematically across the mineral zones, from grossular‐rich garnet in the garnet zone to zoned garnet with grossular‐rich cores and pyrope‐rich rims in the staurolite zone, and pyrope‐rich garnet in the sillimanite and cordierite zones. Thermobarometric analyses of individual garnet crystals reveal two main stages of metamorphism: (i) a high‐P and low‐T stage (as recorded by garnet in the garnet zone and garnet cores in the staurolite zone) at 520–580 °C and 4.5–7 kbar in the kyanite stability field and (ii) a low‐P and high‐T stage (garnet rims in the staurolite zone and garnet in the sillimanite and cordierite zones) at 570–680 °C and 3.0–6.0 kbar in the sillimanite stability field. The earlier high‐P metamorphism resulted in the growth of kyanite in quartz veins within the staurolite and sillimanite zones. The U–Pb zircon ages of pelitic gneisses and granitoids reveal that (i) the protolith (igneous) age of the pelitic gneisses is c. 510 Ma; (ii) the low‐P and high‐T metamorphism occurred at 377 ± 30 Ma; and (iii) this metamorphic stage was coeval with granitoid intrusion at 385 ± 7 Ma. The age of the earlier low‐T and high‐P metamorphism is not clearly recorded in the zircon, but probably corresponds to small age peaks at 450–400 Ma. The low‐P and high‐T metamorphism continued for c. 100 Ma, which is longer than the active period of a single granitoid body. These findings indicate that an elevation of geotherm and a transition from high‐P and low‐T to low‐P and high‐T metamorphism occurred, associated with continuous emplacement of several granitoids, during the crustal evolution in the Devonian CAOB.     

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.     

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.     

Anatectic record and contrasting P–T paths of aluminous gneisses from the central Grenville Province

Anatectic aluminous gneisses, some derived from sedimentary rocks of broadly pelitic composition and others from hydrothermally altered felsic volcanic rocks, are exposed in the mid‐P and high‐P segments of the hinterland in the central Grenville Province. These gneisses consist dominantly of garnet, biotite, K‐feldspar, plagioclase and quartz, with sillimanite or kyanite, and display microstructural evidence of anatexis by fluid‐absent reactions consuming muscovite and/or biotite. Melt‐related microstructures, such as inter‐granular films and/or interstitial quartz or feldspar enclosing relict phases, are most abundant in the metasedimentary samples. Despite anatexis at granulite facies conditions, the hydrothermally altered rocks preserve earlier features attributed to the circulation of hydrothermal fluids, such as sillimanite seams, dismembered quartz veins and garnet‐rich aluminous nodules in a K‐feldspar‐dominated matrix. Microstructural and mineral chemical data, integrated with P–T pseudosections calculated with thermocalc for the metasedimentary rocks, permit qualitative constraints on the P–T paths. Data from a high‐P kyanite‐bearing sample are consistent with a steep prograde P–T path up to ~14.5 kbar and 860–900 °C, followed by decompression with minor cooling to the solidus at ~11 kbar and 870 °C. This pressure‐dominated P–T path is similar to those inferred in other parts of the high‐P segment in the central Grenville Province. In contrast, the P–T path predicted from a mid‐P sillimanite‐bearing paragneiss has a strong temperature gradient with P–T of ~9.5 kbar and 850 °C at the thermal peak, and a retrograde portion down to ~8 kbar and 820 °C. In a broad sense, these two contrasting P–T patterns are consistent with predictions of thermo‐mechanical modelling of large hot orogens in which P–T paths with strong pressure gradients exhume deeper rocks in the orogenic flanks, whereas P–T paths with strong temperature gradients in the orogenic core reflect protracted lateral transport of ductile crust beneath a plateau.     

Metamorphic P‒T paths and Zircon U–Pb ages of Archean ultra-high temperature paragneisses from the Qian’an gneiss dome,East Hebei terrane,North China Craton

The East Hebei terrane of North China Craton is characterized by the dome-and-keel structure, a common feature in most Archean cratons, where supracrustal rocks of granulite facies commonly occur as enclaves or rafts in tonalite–trondhjemite–granodiorite (TTG) gneisses. The metamorphic P–T paths of the granulites are significant for addressing the Archean tectonic regimes. Two types of granulite facies paragneiss with pelitic and greywacke compositions from the western margin of Qian'an gneiss dome are documented for their petrography, mineral chemistry, phase equilibria modelling using thermocalc, and zircon dating. Anticlockwise P–T paths involving the pre-peak pressure increase to the ultra-high temperature peak conditions and post-peak cooling and decompression processes were recognized. The pre-peak pressure increase process was constrained for a pelitic granulite mainly based on the spinel and cordierite inclusions in garnet and rutile corona around ilmenite, where the transition from spinel to garnet is modelled at 6–7 kbar at a fixed T = 1,000°C. For greywacke granulite, the pre-peak pressure increase evolution can be ascertained from the textural relation that orthopyroxene is surrounded by garnet, and the outwards increasing grossular (from 0.03 to 0.05) in the core of the atoll-like garnet (Grt-A), to occur from ~7 kbar at ~1,000°C. The peak P–T conditions for pelitic granulite are roughly limited to 7–11 kbar/890–1,050°C on the basis of the stability of the inferred peak assemblage involving garnet, perthite, sillimanite, rutile/ilmenite, and quartz. The peak P–T condition for greywacke granulite can be well constrained as 9–10 kbar/>1,000°C on the basis of the maximum grossular content (X_Grs = 0.045–0.050) in the core of subhedral garnet (Grt-B) and the mantle of Grt-A together with an average re-integrated anorthite content (X_An = 0.07) in K-feldspar. The peak temperature condition is consistent with the ternary feldspar thermometer results mostly of 950–1,020°C for antiperthite and perthite in greywacke granulite, and in accordance with the development of oriented needle-like exsolution of Ti±Fe oxides in garnet from pelitic granulite. The post-peak cooling and decompression process was consistent with the decreasing X_Grs in the mantle of Grt-A and core of Grt-B in greywacke sample, and the final-stage cooling conditions can be well constrained from the stability of final assemblages marked by the later growth of biotite, as 8–9 kbar/820–880°C for pelitic granulite and 6–9 kbar/840–890°C for greywacke granulite. Zircon dating yields provenance ages from 3.34 to 2.57 Ga and metamorphic ages of c. 2.50 Ga for the two types of granulite. The metamorphic ages overlap the final pulse of the Neoarchean magmatic activity of TTGs that ranges from c. 2.56 to c. 2.48 Ga with a peak at c. 2.52 Ga. Combining the development of dome-and-keel structures, the penecontemporaneity between the metamorphism of supracrustal rocks and TTG magmatic activity, and also the unique anticlockwise P–T paths, we prefer a vertical sagduction regime to interpret the tectonic evolution of the East Hebei terrane, which may be also significant for other Archean cratons.     

P–T–t evolution of garnet amphibolites in the Wutai–Hengshan area,North China Craton: insights from phase equilibria and geochronology

Garnet amphibolites can provide valuable insights into geological processes of orogenic belts, but their metamorphic evolution is still poorly constrained. Garnet amphibolites from the Wutai–Hengshan area of the North China Craton mainly consist of garnet, hornblende, plagioclase, quartz, rutile and ilmenite, with or without titanite and epidote. Four samples selected in a south–north profile were studied by the pseudosection approach in order to elucidate the characteristics of their metamorphic evolution, and to better reveal the northwards prograde change in P–T conditions as established previously. For the sample from the lower Wutai Subgroup, garnet exhibits obvious two‐substage growth zoning characteristic of pyrope (X_py) increasing but grossular (X_gr) decreasing outwards in the core, and both X_py and X_gr increasing outwards in the rim. Phase modelling using thermocalc suggests that the garnet cores were formed by chlorite breakdown over 7–9 kbar at 530–600 °C, and rims grew from hornblende and epidote breakdown over 9.5–11.5 kbar at 600–670 °C. The isopleths of the minimum An in plagioclase and maximum X_py in garnet were used to constrain the peak P–T conditions of ~11.5 kbar/670 °C. The modelled peak assemblage garnet + hornblende + epidote+ plagioclase + rutile + quartz matches well the observed one. Plagioclase–hornblende coronae around garnet indicate post‐peak decompression and fluid ingress. For the samples from the south Hengshan Complex, the garnet zoning weaken gradually, reflecting modifications during decompression of the rocks. Using the same approach, the rocks are inferred to have suprasolidus peak conditions, increasing northwards from 11.5 kbar/745 °C, 12.5 kbar/780 °C to 13 kbar/800 °C. Their modelled peak assemblages involve diopside, garnet, hornblende, plagioclase, rutile and quartz, yet diopside is not observed petrographically. The post‐peak decompression is characterized by diopside + garnet + quartz + melt = hornblende + plagioclase, causing the diopside consumption and garnet compositions to be largely modified. Thus, the pesudosection approach is expected to provide better pressure results than conventional thermobarometry, because the later approach cannot be applied with confidence to rocks with multi‐generation assemblages. U–Pb dating of zircon in the Wutai sample records a protolith age of c. 2.50 Ga, and a metamorphic age of c. 1.95 Ga, while zircon in the Hengshan samples records metamorphic ages of c. 1.92 Ga. The c. 1.95 Ga is interpreted to represent the pre‐peak or peak metamorphic stages, and the ages of c. 1.92 Ga are assigned to represent the cooling stages. All rocks in the Wutai–Hengshan area share similar clockwise P–T morphologies. They may represent metamorphic products at different crustal depths in one orogenic event, which included a main thickening stage at c. 1.95 Ga followed by a prolonged uplift and cooling after 1.92 Ga.