To constrain the rate of cooling of lower-crustal rocks from an ultrahot orogen, we determined both the age and equilibration temperature of metamorphic zircon from six widely spaced samples of metasedimentary garnet-sillimanite gneiss from the Eastern Ghats Province in eastern India. For the combined data set of metamorphic zircon, concordant dates decrease continuously within 2 sigma uncertainty from around 950 Ma to 800 Ma, consistent with similar to 150 m.y. of zircon crystallization. Ti-in-zircon temperatures for each dated spot during this period decrease with age, corresponding to linear cooling rates ranging from 0.26 to 0.90 degrees C/m.y. We propose that retention of heat-producing elements in the lower crust of the Eastern Ghats Province and a low net erosion rate were responsible for similar to 150 m.y. of ultraslow cooling.

The rate of cooling of metamorphic rocks provides a first-order constraint on the tectonic processes controlling heat flow and exhumation. For example, for small crustal terranes that were subducted to ultrahigh pressure conditions during the early stages of collisional orogenesis, exhumation is generally fast with rates similar to plate velocities, such that cooling is also rapid. Similarly, rates of cooling are commonly fast (generally ~20–30°C/Myr) during exhumation of metamorphic core complexes or due to transpression. By contrast, cooling in some granulite terranes can be slow and close-to-isobaric, leading to time-integrated cooling rates of <5°C/Myr. The implication of such slow rates of cooling is that these granulite terranes were close to isostatic equilibrium as a result of sustained high mantle heat flow that limited exhumation by erosion. However, constraining initial cooling rates in granulite terranes can be difficult, particularly where the rocks reached ultrahigh temperatures (>900 °C) that exceed the closure temperature of many geochronometers. In order to overcome this difficulty, we combine U–Pb zircon geochronology with Ti-in-zircon thermometry to investigate the thermal history of metapelitic rocks from the Eastern Ghats Province of eastern India. For the combined dataset of metamorphic zircon from the samples, concordant dates decrease continuously within 2σ uncertainty from around 950 Ma to 800 Ma, consistent with c. 150 Ma of zircon crystallization. Ti-in-zircon temperatures for each dated spot during this period decrease with age, corresponding to linear cooling rates ranging from 0.26 to 0.90°C/Myr. We propose that retention of heat producing elements in the lower crust of the Eastern Ghats Province and a low net erosion rate were responsible for c. 150 Myr of ultra-slow cooling. The location of the Eastern Ghats Province on the margin of the supercontinent Rodinia may have been a contributing factor enabling the region to remain relatively undisturbed until it was exhumed during the formation of Gondwana.

Phase equilibrium modelling is commonly employed to constrain the pressure–temperature ( P–T ) evolution of granulite facies rocks, from which their geodynamic setting may be inferred. However, defining a suitable equilibrium volume in such rocks is non‐trivial due to heterogeneities in protolith composition and open system behaviour, including melt loss or gain and deformation. Consequently, equilibrium volumes and the mineral assemblages they contain may vary both temporally and spatially within a single rock. Additionally, despite the extreme temperatures they attained, granulites commonly contain microtextures indicating gradients in chemical potential after the metamorphic peak. This study assesses the processes that control compositional heterogeneity between equilibrium volumes in a suprasolidus granulite and the consequences for the development of a range of microstructures. An approach combining phase equilibrium forward modelling and chemical potential diagrams is used to simulate the suprasolidus evolution of two adjacent equilibrium volumes in an Mg‐rich metapelite. Assuming equilibrium within the two compositional domains, evolving granulite facies assemblages vary along a high T/P (125°C/kbar) clockwise P–T path representative of high T/P metamorphic terranes. Many retrograde microtextures in metapelitic rocks can be reproduced by considering chemical potential (µ) gradients in µ MgO , µ FeO and µ CaO between phases, assuming (a) the presence of melt on grain boundaries, (b) that Al 2 O 3 is perfectly immobile and (c) that K 2 O and Na 2 O are perfectly mobile. Additionally, documented microtextures in granulites imply local SiO 2 ‐undersaturation, despite the presence of matrix quartz in the rock. Preserving such chemical potential gradients requires that solid phases are chemically and physically isolated from melt. Efficient prograde melt drainage may lead to a loss in melt interconnectivity at high temperature (>900°C), meaning that equilibration of centimetre‐scale compositional domains is controlled by solid‐state diffusion. The presence of isolated pockets of viscous melt allows the formation of discontinuous intergranular reaction microtextures.

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) H2O‐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.

The relative validity and closure temperature of the Zr-in-rutile thermometer for recording UHT metamorphism are process dependent and hotly debated. We present an integrated petrological approach to Zr-in-rutile thermometry including phase equilibrium (pseudosection) modelling in complex chemical systems with updated mineral a-X models and systematic in-situ microanalysis of rutile. This study is centred on high-pressure rutile bearing UHT granulites from Mt. Charles, Napier Complex, Antarctica. P–T phase equilibrium modelling of two garnet bearing granulites (samples 49677, 49701) constrains an overall post-peak near isobaric cooling (IBC) evolution for the Napier Complex at Mt. Charles; from ~14kbar, 1100°C with moderate decompression to ~11kbar, 800–900°C. Local hydration on cooling over this temperature range is recorded in a kyanite bearing granulite (sample 49688) with an inferred injection of aqueous fluid equivalent to up to 9mol% H2O from T–MH2O modelling. Further late stage cooling to <740°C is recorded by voluminous retrograde mica growth and partial preservation of a ky-pl-kfs-bt-liq bearing equilibrium assemblage. Overall, Zr-in-rutile temperatures at 11kbar (Tomkins et al., 2007) are reset to between 606°C and 780°C across all samples, with flat core–rim Zr concentration profiles in all rutiles. However, zircon precipitates as inclusions, needle exsolutions, or rods along rutile grain boundaries are recrystallised from rutiles in qz/fsp domains. Reintegrating the Zr-in-rutile concentration ‘lost’ via the recrystallisation of these zircon precipitates (e.g. Pape et al., 2016) can recover maximum concentrations of up to ~2.2wt% and thus maximum peak temperatures of ~1149°C at 11kbar. Rutile Nb–Ta signatures and rounded rutile grains without zircon precipitates in hydrated mica domains in sample 49688 provide evidence for fluid-mediated mobility of Zr and Nb during retrograde cooling in hydrated lithologies. Aqueous fluid supplemented excess H2O liberated by melt crystallisation, interacting with rutile on cooling of sample 49688 to reset Zr-in-rutile temperatures (606–780°C) at the Mt. Charles locality. The wide range of geochemical and petrological characteristics of Napier Complex rutile highlights that Zr-in-rutile reintegration and a broad petrological approach are required for successful interpretation of Zr-in-rutile geothermometry for long-lived regional UHT metamorphism.