The stability of sapphirine + quartz: calculated phase equilibria in FeO–MgO–Al2O3–SiO2–TiO2–O

The presence in rocks of coexisting sapphirine + quartz has been widely used to diagnose conditions of ultra‐high‐temperature (UHT) metamorphism (>900 °C), an inference based on the restriction of this assemblage to temperatures >980 °C in the conventionally considered FeO–MgO–Al₂O₃–SiO₂ (FMAS) chemical system. With a new thermodynamic model for sapphirine that includes Fe₂O₃, phase equilibra modelling using thermocalc software has been undertaken in the FeO–MgO–Al₂O₃–SiO₂–O (FMASO) and FeO–MgO–Al₂O₃–SiO₂– TiO₂–O (FMASTO) chemical systems. Using a variety of calculated phase diagrams for quartz‐saturated systems, the effects of Fe₂O₃ and TiO₂ on FMAS phase relations are shown to be considerable. Importantly, the stability field of sapphirine + quartz assemblages extends down temperature to 850 °C in oxidized systems and thus out of the UHT range.     

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.     

Quantitative P–T–X constraints on orthopyroxene-bearing high-pressure granulites in felsic–metapelitic rocks: evidence from the Huangtuling granulite, Dabieshan Orogen

High‐pressure granulites are generally characterized by the absence of orthopyroxene. However, orthopyroxene is reported in a few high‐pressure, felsic–metapelitic granulites, such as the Huangtuling felsic high‐pressure granulite in the North Dabie metamorphic core complex in east‐central China, which rarely preserves the high‐pressure granulite facies assemblage of garnet + orthopyroxene + biotite + plagioclase + K‐feldspar + quartz. To investigate the effects of bulk‐rock composition on the stability of orthopyroxene‐bearing, high‐pressure granulite facies assemblages in the NCKFMASHTO (Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–Fe₂O₃) system, a series of P–T–X pseudosections based on the melt‐reintegrated composition of the Huangtuling felsic high‐pressure granulite were constructed. Calculations demonstrate that the orthopyroxene‐bearing, high‐pressure granulite facies assemblages are restricted to low X_Al [Al₂O₃/(Na₂O + CaO + K₂O + FeO + MgO + Al₂O₃) < 0.35, mole proportion] or high X_Mg [MgO/(MgO + FeO) > 0.85] felsic–metapelitic rock types. This study also reveals that the X_Al values in the residual felsic–metapelitic, high‐pressure granulites could be significantly reduced by a high proportion of melt loss. We suggest that orthopyroxene‐bearing, high‐pressure granulites occur in residual overthickened crustal basement under continental subduction–collision zones and arc–continent collision belts.     

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.     

Progress relating to calculation of partial melting equilibria for metapelites

Improved activity–composition relationships for biotite, garnet and silicate liquid are used to construct updated P–T grids and pseudosections for high‐grade metapelites. The biotite model involves Ti charge‐balanced by hydrogen deprotonation on the hydroxyl site, following the substitution , where HD represents the hydroxyl site. Relative to equivalent biotite‐breakdown melting reactions in P–T grids in K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O (KFMASH), those in K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O₂ (KFMASHTO) occur at temperatures close to 50 °C higher. A further consequence of the updated activity models is that spinel‐bearing equilibria occur to higher temperature and higher pressure. In contrast, the addition of Na₂O and CaO to KFMASH to make the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O (NCKFMASH) system lowers key biotite‐breakdown melting reactions in P–T space relative to KFMASH. Combination of the KFMASHTO and NCKFMASH systems to make Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O₂ (NCKFMASHTO) results in key biotite‐breakdown melting reactions occurring at temperatures intermediate between those in KFMASHTO and those in NCKFMASH. Given such differences, the choice of model system will be critical to inferred P–T conditions in the application of mineral equilibria modelling to rocks. Further, pseudosections constructed in KFMASH, NCKFMASH and NCKFMASHTO for several representative rock compositions show substantial differences not only in the P–T conditions of key metamorphic assemblages but also overall topology, with the calculations in NCKFMASHTO more reliably reflecting equilibria in rocks. Application of mineral equilibria modelling to rocks should be undertaken in the most comprehensive system possible, if reliable quantitative P–T information is to be derived.     

Quantitative phase petrology of cordierite–orthoamphibole gneisses and related rocks

Cordierite–orthoamphibole gneisses and rocks of similar composition commonly contain low‐variance mineral assemblages that can provide useful information about the metamorphic evolution of a terrane. New calculated petrogenetic grids and pseudosections are presented in the FeO–MgO–Al₂O₃–SiO₂–H₂O (FMASH), Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O (NCKFMASH) and Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–Fe₂O₃ (NCKFMASHTO) chemical systems to investigate quantitatively the phase relations in these rocks. Although the bulk compositions of cordierite–orthoamphibole gneisses are close to FMASH, calculations in this system do not adequately account for the observed range of mineral assemblages. Calculations in NCKFMASH and NCKFMASHTO highlight the role of minor constituents such as Ca, Na and Fe³⁺ in the mineral assemblage evolution of such rocks and these systems are more appropriate for interpreting the evolution of natural examples.     

Sapphirine parageneses from the Caraiba complex, Bahia, Brazil: the influence of Fe2+–Fe3+ distribution on the stability of sapphirine in natural assemblages

Abstract Sapphirine-bearing rocks occur in three conformable, metre-size lenses in intrusive quartzo-feldspathic orthogneisses in the Curaçà valley of the Archaean Caraiba complex of Brazil. In the lenses there are six different sapphirine-bearing rock types, which have the following phases (each containing phlogopite in addition): A: Sapphirine, orthopyroxene; B: Sapphirine, cordierite, orthopyroxene, spinel; C: Sapphirine, cordierite; D: Sapphirine, cordierite, orthopyroxene, quartz; E: Sapphirine, cordierite, orthopyroxene, sillimanite, quartz; F: Sapphirine, cordierite, K-feldspar, quartz. Neither sapphirine and quartz nor orthopyroxene and sillimanite have been found in contact, however. During mylonitization, introduction of silica into the three quartz-free rocks (which represent relict protolith material) gave rise to the three cordierite and quartz-bearing rocks. Stable parageneses in the more magnesian rocks were sapphirine–orthopyroxene and sapphirine–cordierite. In more iron-rich rocks, sapphirine–cordierite, sapphirine-cordierite–sillimanite, cordierite–sillimanite, sapphirine–cordierite–spinel–magnetite and quartz–cordierite–orthopyroxene were stable. The iron oxide content in sapphirine of the six rocks increases from an average of 2.0 to 10.5 wt % (total Fe as FeO) in the order: C,F–A,D–B,E. With increase in Fe there is an increase in recalculated Fe₂O₃ in sapphirine. The four rock types associated with the sapphirine-bearing lenses are: I: Orthopyroxene, cordierite, biotite, quartz, feldspar tonalitic to grandioritic gneiss; II: Biotite, quartz, feldspar gneiss; III: Orthopyroxene, clinopyroxene, hornblende, plagioclase meta-norite; IV: Biotite, orthopyroxene, quartz, feldspar, garnet, cordierite, sillimanite granulite gneiss. The stable parageneses in type IV are orthopyroxene–cordierite–quartz, garnet–sillimanite–quartz and garnet–cordierite–sillimanite. Geothermobarometry suggests that the associated host rocks equilibrated at 720–750°C and 5.5–6.5 kbar. Petrogenetic grids for the FMASH and FMAFSH (FeO–MgO–Al₂O₃–Fe₂O₃–SiO₂–H₂O) model systems indicate that sapphirine-bearing assemblages without garnet were stabilized by a high Fe³⁺ content and a high X_Mg= (Mg/ (Mg+Fe²⁺)) under these P–T conditions.     

A petrogenetic grid for aluminous granulite facies metapelites in the KFMASH system

Due to the retrograde cation exchange problems experienced by conventional geothermobarometers above their closure temperatures, petrogenetic grids are a potentially powerful alternative to unravelling the P–T evolution of ultrahigh‐T granulite terranes. A new qualitative KFMASH (K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O) petrogenetic grid for Mg–Al rich metapelites containing K‐feldspar, sillimanite and quartzofeldspathic melt that successfully accounts for the majority of assemblages composed of variations of sapphirine, spinel, garnet, orthopyroxene, cordierite, biotite and quartz is developed. Univariant reactions are predicted utilizing a newly derived ‘melt projection’ and these reactions are entirely consistent with algebraically calculated reaction coefficients obtained using a set of standard phase compositions. Based upon observations of commonly associated mineral assemblages in natural lithologies the [Spr, Spl], [Qtz, Spl], [Bt, Spl], [Opx, Spr], [Opx, Qtz] and [Bt, Opx] invariant points are assumed to be stable, whilst the [Grt, Spr], [Grt, Qtz], [Spr, Qtz] and [Crd, Qtz] are assumed to be metastable. Biotite‐bearing assemblages are confined to the lowest temperatures, and sapphirine + quartz to the highest temperatures. Orthopyroxene + sillimanite ± quartz assemblages are confined to the highest pressures, whilst spinel‐bearing assemblages are stabilized by lower pressures. The alternative choice of invariant point stability leads to significant differences between this grid and previously proposed topologies. Spinel cannot be stable along with the orthopyroxene and sillimanite assemblage as previously proposed. Further, more subtle differences in topology result from the treatment of H₂O in the chemographic projection used to deduce univariant reactions, and projecting from a water‐bearing quartzofeldspathic melt does not yield the same reaction coefficients as projection from H₂O. The new grid allows reinterpretation of previously proposed evolutionary P–T paths for Mg–Al rich granulites from the Napier Complex and Rauer Group, East Antarctica, and In Ouzzal, Algeria.     

Calculated phase equilibria in K2O-FeO-MgO-Al2O3-SiO2-H2O for sapphirine-quartz-bearing mineral assemblages

Sapphirine, coexisting with quartz, is an indicator mineral for ultrahigh‐temperature metamorphism in aluminous rock compositions. Here a new activity‐composition model for sapphirine is combined with the internally consistent thermodynamic dataset used by THERMOCALC, for calculations primarily in K₂O‐FeO‐MgO‐Al₂O₃‐SiO₂‐H₂O (KFMASH). A discrepancy between published experimentally derived FMAS grids and our calculations is understood with reference to H₂O. Published FMAS grids effectively represent constant a_H2O sections, thereby limiting their detailed use for the interpretation of mineral reaction textures in compositions with differing H₂O. For the calculated KFMASH univariant reaction grid, sapphirine + quartz assemblages occur at P–T in excess of 6–7 kbar and 1005 °C. Sapphirine compositions and composition ranges are consistent with natural examples. However, as many univariant equilibria are typically not ‘seen’ by a specific bulk composition, the univariant reaction grid may reveal little about the detailed topology of multi‐variant equilibria, and therefore is of limited use for interpreting the P–T evolution of mineral assemblages and reaction sequences. Calculated pseudosections, which quantify bulk composition and multi‐variant equilibria, predict experimentally determined KFMASH mineral assemblages with consistent topology, and also indicate that sapphirine stabilizes at increasingly higher pressure and temperature as X_Mg increases. Although coexisting sapphirine and quartz can occur in relatively iron‐rich rocks if the bulk chemistry is sufficiently aluminous, the P–T window of stability shrinks with decreasing X_Mg. An array of mineral assemblages and mineral reaction sequences from natural sapphirine + quartz and other rocks from Enderby Land, Antarctica, are reproducible with calculated pseudosections. That consistent phase diagram calculations involving sapphirine can be performed allows for a more thorough assessment of the metamorphic evolution of high‐temperature granulite facies terranes than was previously possible. The establishment of a a‐x model for sapphirine provides the basis for expansion to larger, more geologically realistic chemical systems (e.g. involving Fe³⁺).     

FeO*-Al2O3-TiO2-Rich Rocks of the Tertiary Bana Igneous Complex, West Cameroon

FeO*‐Al₂O₃‐TiO₂‐rich rocks are found associated with transitional tholeiitic lava flows in the Tertiary Bana plutono‐volcanic complex in the continental sector of the Cameroon Line. These peculiar rocks consist principally of iron‐titanium oxides, aluminosilicates and phosphates, and occur as layers 1–3 m thick occupying the upper part of lava flows on the southwest (site 1) and northwest (site 2) sites of the complex. Mineral constituents of the rocks include magnetite, ilmenite, hematite, rutile, corundum, andalusite, sillimanite, cordierite, quartz, plagioclase, alkali feldspar, apatite, Fe‐Mn phosphate, Al phosphate, micas and fine mixtures of sericite and silica. Texturally and compositionally, the rocks can be subdivided into globular type, banded type, and Al‐rich fine‐gained massive type. The first two types consist of dark globule or band enriched in Fe‐Ti oxides and apatite and lighter colored groundmass or bands enriched in aluminosilicates and quartz, respectively. The occurrence of andalusite and sillimanite and the compositional relations of magnetite and ilmenite in the FeO*‐Al₂O₃‐TiO₂‐rich rocks suggest temperatures of crystallization in a range of 690–830°C at low pressures. The Bana FeO*‐Al₂O₃‐TiO₂‐rich rocks are characterized by low concentrations of SiO₂ (25–54.2 wt%), Na₂O + K₂O (0–1%), CaO (0–2%) and MgO (0–0.5%), and high concentrations of FeO* (total iron as FeO, 20–42%), Al₂O₃ (20–42%), TiO₂ (3–9.2%), and P₂O₅ (0.26–1.30%). TiO₂ is positively correlated with Al₂O₃ and inversely correlated with FeO*. The bulk rock compositions cannot be derived from the associated basaltic magma by crystal fractionation or by partial melting of the mantle or lower crustal materials. In ternary diagrams of (Al₂O₃)?(CaO + Na₂O + K₂O)?(FeO*+ MnO + MgO) and (SiO₂)?(FeO*)?(Al₂O₃), the compositional field of the rocks is close to that of laterite and is distinct from the common volcanic rocks, suggesting that the rocks are derived from lateritic materials by recrystallization when the materials are heated by the basaltic magmas. A hydrothermal origin is discounted because the rocks contain high‐temperature mineral assemblages and lack sulfide minerals. It is proposed that the FeO*‐Al₂O₃‐TiO₂‐rich rocks of the Bana complex were formed by pyrometamorphism of laterite by the heat of basaltic magmas.     

Using calculated chemical potential relationships to account for replacement of kyanite by symplectite in high pressure granulites

Calculated mineral equilibria are used to account for the formation of sapphirine–plagioclase, spinel–plagioclase and corundum–plagioclase symplectites replacing kyanite in quartz–plagioclase–garnet–kyanite granulite facies gneisses from the Southern Domain of the Athabasca granulite terrane, a segment of the Snowbird tectonic zone in northern Saskatchewan, Canada. Metamorphic conditions of >14 kbar and 800 °C are established for the high pressure, garnet–kyanite assemblage using constraints from P–T pseudosections and Zr‐in‐rutile thermometry. Replacement of kyanite by symplectites reflects the reaction of kyanite with the matrix following near‐isothermal decompression to <10 kbar. The chemical potential gradients developed between the kyanite and the matrix led to diffusion that attempted to flatten the gradients, kyanite persisting as a stable phase while it is consumed by symplectite from its edge. In this local equilibrium model, the mineral and mineral compositional spatial relationships are shown to correspond to paths in μ(Na₂O)–μ(CaO)–μ(K₂O)–μ(FeO)–μ(MgO) in the model chemical system, Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂ (NCKFMAS), with SiO₂ and Al₂O₃ taken to be completely immobile. The values of μ(Na₂O) and μ(CaO) are constrained by fixing P–T conditions and choosing appropriate μ(Na₂O) and μ(CaO) values that correspond to the observed plagioclase compositions. μ(FeO)–μ(MgO) diagrams show the corresponding spatial relationships with kyanite and the symplectite phases. These results demonstrate that the replacement of kyanite by sapphirine–plagioclase and spinel–plagioclase appears to be metastable with respect to replacement by corundum–plagioclase. Replacement by corundum–plagioclase does also occur, apparently overprinting pre‐existing symplectite and also kyanite. Ignoring corundum, the resulting diagrams account for the spatial relationships and compositions observed in the spinel–plagioclase and sapphirine–plagioclase symplectites. They are predicted to occur over both a wide range of P–T conditions (6–11 kbar, 650–850 °C) and plagioclase compositions (X_An = 0.5–0.9). The wide range of P–T conditions that may result in identical spatial and compositional relationships suggests that such reaction textures may be of limited use in accurately quantifying the P–T conditions of retrograde metamorphism.     

Calculated phase equilibria in K2O-FeO-MgO-Al2O3-SiO2-H2O for silica-undersaturated sapphirine-bearing mineral assemblages

Silica‐undersaturated, sapphirine‐bearing granulites occur in a large number of localities worldwide. Such rocks have historically been under‐utilized for estimating P–T evolution histories because of limited experimental work, and a consequent poor understanding of the topology and P–T location of silica‐undersaturated mineral equilibria. Here, a calculated P–T projection for sapphirine‐bearing, silica‐undersaturated metapelitic rock compositions is constructed using THERMOCALC for the FeO‐MgO‐Al₂O₃‐SiO₂ (FMAS) and KFMASH (+K₂O + H₂O) chemical systems, allowing quantitative analysis of silica‐undersaturated mineral assemblages. This study builds on that for KFMASH sapphirine + quartz equilibria [Kelsey et al. (2004) Journal of Metamorphic Geology, vol. 22, pp. 559–578]. FMAS equilibria are significantly displaced in P–T space from silicate melt‐bearing KFMASH equilibria. The large number of univariant silica‐undersaturated KFMASH equilibria result in a P–T projection that is topologically more complex than could be established on the basis of experiments and/or natural assemblages. Coexisting sapphirine and silicate melt (with or without corundum) occur down to c. 900 °C in KFMASH, some 100 °C lower than in silica‐saturated compositions, and from pressures of c.≤1 to ≥12 kbar. Mineral compositions and composition ranges for the calculated phases are consistent with natural examples. Bulk silica has a significant effect on the stability of sapphirine‐bearing assemblages at a given P–T, resulting in a wide variety of possible granulite facies assemblages in silica‐undersaturated metapelites. Calculated pseudosections are able to reproduce many naturally occurring silica‐undersaturated assemblages, either within a single assemblage field or as the product of a P–T trajectory crossing several fields. With an understanding of the importance of bulk composition on sapphirine stability and textural development, silica‐undersaturated assemblages may be utilized in a quantitative manner in the detailed metamorphic investigation of high‐grade terranes.     

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.     

Petrogenesis of andalusite–kyanite–sillimanite veins and host rocks, Sanandaj-Sirjan metamorphic belt, Hamadan, Iran

Quartz‐rich veins in metapelitic schists of the Sanandaj‐Sirjan belt, Hamadan region, Iran, commonly contain two Al₂SiO₅ polymorphs, and, more rarely, three coexisting Al₂SiO₅ polymorphs. In most andalusite and sillimanite schists, the types of polymorphs in veins correlate with Al₂SiO₅ polymorph(s) in the host rocks, although vein polymorphs are texturally and compositionally distinct from those in adjacent host rocks; e.g. vein andalusite is enriched in Fe₂O₃ relative to host rock andalusite. Low‐grade rocks contain andalusite + quartz veins, medium‐grade rocks contain andalusite + sillimanite + quartz ± plagioclase veins, and high‐grade rocks contain sillimanite + quartz + plagioclase veins/leucosomes. Although most andalusite and sillimanite‐bearing veins occur in host rocks that also contain Al₂SiO₅, kyanite‐quartz veins crosscut rocks that lack Al₂SiO₅ (e.g. staurolite schist, granite). A quartz vein containing andalusite + kyanite + sillimanite + staurolite + muscovite occurs in andalusite–sillimanite host rocks. Textural relationships in this vein indicate the crystallization sequence andalusite to kyanite to sillimanite. This crystallization sequence conflicts with the observation that kyanite‐quartz veins post‐date andalusite–sillimanite veins and at least one intrusive phase of a granite that produced a low‐pressure–high‐temperature contact aureole; these relationships imply a sequence of andalusite to sillimanite to kyanite. Varying crystallization sequences for rocks in a largely coherent metamorphic belt can be explained by P–T paths of different rocks passing near (slightly above, slightly below) the Al₂SiO₅ triple point, and by overprinting of multiple metamorphic events in a terrane that evolved from a continental arc to a collisional orogen.     

Inference of a detailed P–T path from P–T pseudosections using metapelitic rocks of variable composition from a single outcrop, Shackleton Range, Antarctica

Generally, P–T pseudosections for reduced compositional systems, such as K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O, Na₂O–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O and MnO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O, are well suited for inferring detailed P–T paths, comparing mineral assemblages observed in natural rocks with those calculated. Examples are provided by P–T paths inferred for four metapelitic samples from a 1 m² wide outcrop of the Herbert Mountains in the Shackleton Range, Antarctica. The method works well if the bulk composition used is reconstituted from average mineral modes and mineral compositions (AMC) or when X‐ray fluorescence (XRF) data are corrected for Al₂O₃ and FeO. A plagioclase correction is suitable for Al₂O₃. Correction for FeO is dependent on additional microscopic observations, e.g. the kind and amount of opaque minerals. In some cases, all iron can be treated as FeO^tot, whereas in others a magnetite or hematite correction yields much better results. Comparison between calculated and observed mineral modes and mineral compositions shows that the AMC bulk composition is best suited to the interpretation of rock textures using P–T pseudosections, whereas corrected XRF data yield good results only when the investigated sample has few opaque minerals. The results indicate that metapelitic rocks from the Herbert Mountains of the Northern Shackleton Range underwent a prograde P–T evolution from about 600 °C/5.5 kbar to 660 °C/7 kbar, followed by nearly adiabatic cooling to about 600 °C at 4.5 kbar.     

Crustal growth during back‐arc closure: Cretaceous exhumation history of Cordillera Darwin,southern Patagonia

The Cordillera Darwin metamorphic complex is unique in the Andes in exposing kyanite–staurolite schist north of the Beagle Channel in southern Patagonia. Garnet in amphibolite facies pelitic schists from Bahía Pia has patchy textures whereby some grains consist of clear, grossular‐rich garnet with fine‐grained S1 inclusion trails truncated by regions of turbid spessartine–pyrope‐rich garnet with biotite, muscovite, plagioclase and quartz inclusions. Micron‐scale aqueous inclusions in turbid garnet are consistent with recrystallization facilitated by fluid ingress; S2 inclusion trails indicate this was broadly contemporary with the growth of kyanite and staurolite in the matrix. Pseudosection modelling in Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–Fe₂O₃ (NCKFMASHTO) is used to infer a P–T path dominated by decompression from 12 to 9 kbar at T ≈ 620 °C, coupled with garnet mode decreasing from ～5% to <1%. U–Th–Pb in situ dating of S2 monazite indicates that staurolite and kyanite growth and thus exhumation was underway before 72.6 ± 1.1 Ma. Contact aureoles developed adjacent to late granite intrusions include sillimanite‐bearing migmatites formed at P ≈ 6 kbar after 72 Ma. Metamorphism of southern Cordillera Darwin induced by continental underthrusting beneath the arc, related to closure of the Rocas Verdes back‐arc basin, was terminated by thrusting‐controlled exhumation, with the rocks at P ≈ 9 kbar by c. 73 Ma and 6 kbar by c. 70 Ma.     

The metamorphic evolution of metapelitic granulites from Radok Lake, northern Prince Charles Mountains, east Antarctica; evidence for an anticlockwise P–T path

Mineral textures in metapelitic granulites from the northern Prince Charles Mountains, coupled with thermodynamic modelling in the K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–Fe₂O₃ (KFMASHTO) model system, point to pressure increasing with increasing temperature on the prograde metamorphic path, followed by retrograde cooling (i.e. an anticlockwise P–T path). Textural evidence for the increasing temperature part of the path is given by the breakdown of garnet and biotite to form orthopyroxene and cordierite in sillimanite‐absent rocks, and through the break‐down of biotite and sillimanite to form spinel, cordierite and garnet in more aluminous assemblages. This is equated to the advective addition of heat from the regional emplacement of granitic and charnockitic magmas dated at c. 980 Ma. A subsequent increase in pressure, inferred from the break‐down of spinel and quartz to sillimanite, cordierite and garnet in aluminous rocks, is attributed to crustal thickening related to upright folding dated at 940–910 Ma. The terrane attained peak metamorphic temperatures of c. 880 °C at pressures of c. 6.0–6.5 kbar during this event. Subsequent cooling is inferred from the localised breakdown of cordierite and garnet to form biotite and sillimanite that developed in the latter stages of the same event. The textural observations described are interpreted via the application of P–T and P–T–X pseudosections. The latter show that most rock compositions preserve only fragments of the overall P–T path; a result of different rock compositions undergoing mineral assemblage changes, or changes in mineral modal abundance, on different sections of the P–T path. The results also suggest that partial melting during granulite facies metamorphism, coupled with melt loss and dehydration, initiated a switch from pervasive ductile, to discrete ductile/brittle deformation, during retrograde cooling.     

Calculation of partial melting equilibria in the system Na2O–CaO–K2O–FeO–MgO–Al2O3–SiO2–H2O (NCKFMASH)

A thermodynamic model for haplogranitic melts in the system Na₂O–CaO–K₂O–Al₂O₃–SiO₂–H₂O (NCKASH) is extended by the addition of FeO and MgO, with the data for the additional end‐members of the liquid incorporated in the Holland & Powell (1998) internally consistent thermodynamic dataset. The resulting dataset, with the software thermocalc , is then used to calculate melting relationships for metapelitic rock compositions. The main forms for this are P–T and T–X pseudosections calculated for particular rock compositions and composition ranges. The relationships in these full‐system pseudosections are controlled by the low‐variance equilibria in subsystems of NCKFMASH. In particular, the solidus relationships are controlled by the solidus relationships in NKASH, and the ferromagnesian mineral relationships are controlled by those in KFMASH. However, calculations in NCKFMASH allow the relationships between the common metapelitic minerals and silicate melt to be determined. In particular, the production of silicate melt and melt loss from such rocks allow observations to be made about the processes involved in producing granulite facies rocks, particularly relating to open‐system behaviour of rocks under high‐grade conditions.     

Experimental study of the stability of cordierite and garnet in pelitic compositions at high pressures and temperatures

The stability of cordierite and garnet relative to their anhydrous breakdown products, i.e. hypersthene, sapphirine, olivine, spinel, sillimanite and quartz, has been studied experimentally in model pelitic compositions (system MgO-FeO-Al₂O₃-CaO-K₂O-SiO₂). Below 1000° C cordierite breaks down according to the divariant reaction cordierite garnet+sillimanite+quartz (1) for most values of the MgO/MgO + FeO ratio (X). At very high values of X (ca. X0.9) garnet in reaction (1) is replaced by hypersthene. The position and width of the divariant field (in terms of pressure and temperature) in which cordierite and garnet coexist, is a function of the MgO/MgO + FeO ratio. If this ratio is increased then the stability field of garnet is reduced and that of cordierite extended towards higher pressure. Compositions of coexisting cordierite and garnet in divariant equilibrium have been analysed by electron probe micro-analyser. These compositions are unique functions of pressure and temperature. Above ca. 1000° C the breakdown of cordierite involves the phases sapphirine and hercynite-rich spinel in Mg-rich and Fe-rich compositions respectively.