Stability of sapphirine + quartz in the oxidized rocks of the Wilson Lake terrane,Labrador: calculated equilibria in NCKFMASHTO

The equilibrium coexistence of sapphirine + quartz is inferred to record temperatures in excess of 980 °C, based on the stability of this assemblage in the simplified chemical system FeO–MgO–Al₂O₃–SiO₂ (FMAS) system. However, the potential for sapphirine to contain significant Fe³⁺ suggests that the stability of sapphirine + quartz could extend to lower temperatures than those constrained in this ideal system. The Wilson Lake terrane in the Grenville Province of central Labrador preserves sapphirine + quartz‐bearing assemblages in highly oxidized bulk compositions, and provides an opportunity to explore the stability of sapphirine + quartz in such rock compositions within the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O (NCKFMASHTO) chemical system. Starting with the phase equilibria in FeO–MgO–Al₂O₃–SiO₂–TiO₂–O (FMASTO), expansion into K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O (KFMASHTO) allows the effect of the stability of the additional phases, biotite, K‐feldspar and melt, on the stability of sapphirine + quartz to be assessed. These phase relations are evaluated generally using P–T projections, and the ultimate extension into NCKFMASHTO is done with pseudosections. Conditions of peak metamorphism in the Wilson Lake terrane are constrained using P–T pseudosections, and the appropriate H₂O and O contents to use in the modelled compositions are investigated using T–M_H2O and T–M_O pseudosections. The peak P–T estimates from a sapphirine + quartz‐bearing sample are ～960 to 935 °C at ～10 to 8.6 kbar, similar to estimates from orthopyroxene + sillimanite + quartz ± garnet‐bearing samples. Whereas the sapphirine + quartz‐bearing sample is more Fe‐rich than the orthopyroxene + sillimanite‐bearing sample on an all‐Fe‐as‐FeO basis, once the oxidation state is taken into account, the former is effectively more magnesian than the latter, accounting for the sapphirine occurrence.     

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.     

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.     

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.     

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.     

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.     

New mineral activity–composition relations for thermodynamic calculations in metapelitic systems

New activity–composition (a–x) relations for minerals commonly occurring in metapelites are presented for use with the internally consistent thermodynamic dataset of Holland & Powell ( 2011 , Journal of Metamorphic Geology, 29 , 333–383). The a–x relations include a broader consideration of Fe₂O₃ in minerals, changes to the formalism of several phases and order–disorder in all ferromagnesian minerals where Fe–Mg mixing occurs on multiple sites. The a–x relations for chlorite, biotite, garnet, chloritoid, staurolite, cordierite, orthopyroxene, muscovite, paragonite and margarite have been substantially reparameterized using the approach outlined in the companion paper in this issue. For the first time, the entire set of a–x relations for the common ferromagnesian minerals in metapelitic rocks is parameterized simultaneously, with attention paid to ensuring that they can be used together to calculate phase diagrams of geologically appropriate topology. The a–x relations developed are for use in the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O₂ (NCKFMASHTO) system for both subsolidus and suprasolidus conditions. Petrogenetic grids in KFMASH and KFMASHTO are similar in topology to those produced with earlier end‐member datasets and a–x relations, but with some notable differences. In particular, in subsolidus equilibria, the FeO/(FeO + MgO) of garnet is now greater than in coexisting staurolite, bringing a number of key staurolite‐bearing equilibria into better agreement with inferences from field and petrographic observations. Furthermore, the addition of Fe³⁺ and Ti to a number of silicate phases allows more plausible equilibria to be calculated in relevant systems. Pseudosections calculated with the new a–x relations are also topologically similar to equivalent diagrams using earlier a–x relations, although with many low variance fields shifting in P–T space to somewhat lower pressure conditions.     

Phase equilibria constraints on melting of stromatic migmatites from Ronda (S. Spain): insights on the formation of peritectic garnet

Stromatic metatexites occurring structurally below the contact with the Ronda peridotite (Ojén nappe, Betic Cordillera, S Spain) are characterized by the mineral assemblage Qtz+Pl+Kfs+Bt+Sil+Grt+Ap+Gr+Ilm. Garnet occurs in low modal amount (2–5 vol.%). Very rare muscovite is present as armoured inclusions, indicating prograde exhaustion. Microstructural evidence of melting in the migmatites includes pseudomorphs after melt films and nanogranite and glassy inclusions hosted in garnet cores. The latter microstructure demonstrates that garnet crystallized in the presence of melt. Re‐melted nanogranites and preserved glassy inclusions show leucogranitic compositions. Phase equilibria modelling of the stromatic migmatite in the MnO–Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂–O₂–C (MnNCaKFMASHOC) system with graphite‐saturated fluid shows P–T conditions of equilibration of 4.5–5 kbar, 660–700 °C. These results are consistent with the complete experimental re‐melting of nanogranites at 700 °C and indicate that nanogranites represent the anatectic melt generated immediately after entering supersolidus conditions. The P–T estimate for garnet and melt development does not, however, overlap with the low‐temperature tip of the pure melt field in the phase diagram calculated for the composition of preserved glassy inclusions in garnet in the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O (NCKFMASH) system. A comparison of measured melt compositions formed immediately beyond the solidus with results of phase equilibria modelling points to the systematic underestimation of FeO, MgO and CaO in the calculated melt. These discrepancies are present also when calculated melts are compared with low‐T natural and experimental melts from the literature. Under such conditions, the available melt model does not perform well. Given the presence of melt inclusions in garnet cores and the P–T estimates for their formation, we argue that small amounts (<5 vol.%) of peritectic garnet may grow at low temperatures (≤700 °C), as a result of continuous melting reactions consuming biotite.     

Phase equilibrium modelling of the amphibolite to granulite facies transition in metabasic rocks (Ivrea Zone,NW Italy)

The development of thermodynamic models for tonalitic melt and the updated clinopyroxene and amphibole models now allow the use of phase equilibrium modelling to estimate P–T conditions and melt production for anatectic mafic and intermediate rock types at high‐T conditions. The Permian mid‐lower crustal section of the Ivrea Zone preserves a metamorphic field gradient from mid amphibolite facies to granulite facies, and thus records the onset of partial melting in metabasic rocks. Interlayered metabasic and metapelitic rocks allows the direct comparison of P–T estimates and partial melting between both rock types with the same metamorphic evolution. Pseudosections for metabasic compositions calculated in the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O (NCKFMASHTO) system are presented and compared with those of metapelitic rocks calculated with consistent end‐member data and a–x models. The results presented in this study show that P–T conditions obtained by phase equilibria modelling of both metabasic and metapelitic rocks give consistent results within uncertainties, allowing integration of results obtained for both rock types. In combination, the calculations for both metabasic and metapelitic rocks allows an updated and more precisely constrained metamorphic field gradient for Val Strona di Omegna to be defined. The new field gradient has a slightly lower dP/dT which is in better agreement with the onset of crustal thinning of the Adriatic margin during the Permian inferred in recent studies.     

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.     

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.     

Age and origin of sillimanite schist from the Chinese Altai metamorphic belt: implications for late Palaeozoic tectonic evolution of the Central Asian Orogenic Belt

Reconstructing late Palaeozoic metamorphism of the Central Asian Orogenic Belt (CAOB) can provide a better understanding of how the CAOB formed. The petrology of sillimanite-bearing metapelitic schists from high-grade portions of the Permian Chinese Altai metamorphic belt (andalusite-type) reflects the effects of poorly understood high-T, low-P metamorphism. Phase equilibria modelling in the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–TiO₂–O (NCKFMASHTO) system restricts P–T conditions of the sillimanite schists to approximately 635–670°C at approximately 5.8–6.8 kbar. SHRIMP U–Pb analyses of zircon from the rocks yield a concordant age of 299.2 ± 3.4 Ma. Combined with the slightly younger (292.8 ± 2.3 Ma) areally restricted pelitic granulite with peak P?T conditions of approximately 780–800°C at approximately 5–6 kbar and high-T granulite with P?T conditions of approximately 860°C at approximately 6 kbar, these metamorphic rocks reflect prograde heating at relatively low pressure in early Permian time. Together with contemporary and widespread magmatic activities, they are best explained in the context of a post-orogenic extensional environment related to a mantle plume.     

The effect of TiO2 and Fe2O3 on metapelitic assemblages at greenschist and amphibolite facies conditions: mineral equilibria calculations in the system K2O–FeO–MgO–Al2O3–SiO2–H2O–TiO2–Fe2O3

Mineral equilibria calculations in the system K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–Fe₂O₃ (KFMASHTO) using thermocalc and its internally consistent thermodynamic dataset constrain the effect of TiO₂ and Fe₂O₃ on greenschist and amphibolite facies mineral equilibria in metapelites. The end‐member data and activity–composition relationships for biotite and chloritoid, calibrated with natural rock data, and activity–composition data for garnet, calibrated using experimental data, provide new constraints on the effects of TiO₂ and Fe₂O₃ on the stability of these minerals. Thermodynamic models for ilmenite–hematite and magnetite–ulvospinel solid solutions accounting for order–disorder in these phases allow the distribution of TiO₂ and Fe₂O₃ between oxide minerals and silicate minerals to be calculated. The calculations indicate that small to moderate amounts of TiO₂ and Fe₂O₃ in typical metapelitic bulk compositions have little effect on silicate mineral equilibria in metapelites at greenschist to amphibolite facies, compared with those calculated in KFMASH. The addition of large amounts of TiO₂ to typical pelitic bulk compositions has little effect on the stability of silicate assemblages; in contrast, rocks rich in Fe₂O₃ develop a markedly different metamorphic succession from that of common Barrovian sequences. In particular, Fe₂O₃‐rich metapelites show a marked reduction in the stability fields of staurolite and garnet to higher pressures, in comparison to those predicted by KFMASH grids.     

Calculated phase equilibria for low- and medium-pressure metapelites in the KFMASH and KMnFMASH systems

Petrogenetic grids in the KFMASH and KMnFMASH model systems calculated with the software thermocalc 3.1 are presented for the P–T range 0.5–12 kbar and 450–900 °C, for assemblages involving garnet, muscovite, chloritoid, biotite, chlorite, staurolite, cordierite, spinel, orthopyroxene, K‐feldspar, Al₂SiO₅ phases, quartz, water and melt. Based on calculated compatibility diagrams and P–T and T–M_Mn [Mn/(Mg + Fe + Mn)] pseudosections for different metapelitic bulk compositions, the principal conclusions are that the addition of Mn to the KFMASH system: (i) enhances the stability of garnet, and, to a lesser extent, aluminosilicates; (ii) reduces the stability of staurolite, cordierite and, to a lesser extent, chlorite; and (iii) extends the medium pressure stability of muscovite and the low‐P stability field of K‐feldspar. The influence of Mn on individual mineral stabilities is strongly related to rock composition, in particular, to the relative contents of Al₂O₃ and K₂O. For metapelites of a range of compositions and M_Mn values, P–T pseudosections in the KFMASH system, in most cases, do not adequately predict the mineral assemblages observed in natural assemblages under medium and low‐pressure conditions. In contrast, the P–T pseudosections in the KMnFMASH system generally provide more satisfactory results, suggesting that MnO is one of the non‐KFMASH components that should not be neglected in documenting the phase equilibria of medium‐ and low‐P metapelites.     

Phase relations in high-pressure metapelites in the system KFMASH (K2O–FeO–MgO–Al2O3–SiO2–H2O) with application to natural rocks

Using a previously published, internally consistent thermodynamic dataset and updated models of activity–composition relations for solid solutions, petrogenetic grids in the model system KFMASH (K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O) and the subsystems KMASH and KFASH have been calculated with the software THERMOCALC 3.1 in the P–T range 5–36 kbar and 400–810 °C, involving garnet, chloritoid, biotite, carpholite, talc, chlorite, staurolite and kyanite/sillimanite with phengite, quartz/coesite and H₂O in excess. These grids, together with calculated AFM compatibility diagrams and pseudosections, are shown to be powerful tools for delineating the phase equilibria and P–T conditions of pelitic high-P assemblages for a variety of bulk compositions. The calculated equilibria and mineral compositions are in good agreement with petrological observation. The calculation indicates that the typical whiteschist assemblage kyanite–talc is restricted to the rocks with extremely high X_Mg values, decreasing X_Mg in a bulk composition favoring the stability of chloritoid and garnet. Also, the chloritoid–talc paragenesis is stable over 19–20 kbar in a temperature range of ca. 520–620 °C, being more petrologically important than the previously highlighted assemblage talc–phengite. Moreover, contours of the calculated Si isopleths in phengite in P–T and P–X pseudosections for different bulk compositions extend the experimentally derived phengite geobarometers to various KFMASH assemblages.     

The lawsonite paradox: a comparison of field evidence and mineral equilibria modelling

Lawsonite equilibria are predicted to occur over a broad P–T spectrum developed during subduction, yet lawsonite‐bearing assemblages are rare. In the context of mafic mineral equilibria modelled for the range of common crustal metamorphism (4–23 kbar, 400–750 °C) using the system Na₂O‐CaO‐K₂O‐FeO‐MgO‐Al₂O₃‐SiO₂‐H₂O and the software thermocalc , unusually high water contents are demanded by lawsonite assemblages. As a consequence, lawsonite assemblages are predicted to have difficulty forming and lawsonite equilibria to be uncommon. Metabasalt undergoing cooler subduction may experience substantial periods involving the metastable persistence of mineral assemblages because of water under‐saturation with non‐occurrence of recrystallization. If formed, lawsonite‐bearing assemblages are observed to be highly unstable; their preservation requires that exhumation be accompanied by substantial cooling. The amount of structurally bound H₂O in minerals plays a critical role in the formation and preservation of mineral assemblages, controlling key changes in rocks undergoing subduction.     

Thermodynamic modelling of the reaction muscovite+cordierite→Al2SiO5+biotite+quartz+ H2O: constraints from natural assemblages and implications for the metapelitic petrogenetic grid

The reaction muscovite+cordierite→biotite+Al₂SiO₅ +quartz+H₂O is of considerable importance in the low pressure metamorphism of pelitic rocks: (1) its operation is implied in the widespread assemblage Ms + Crd +And± Sil + Bt + Qtz, a common mineral assemblage in contact aureoles and low pressure regional terranes; (2) it is potentially an important equilibrium for pressure estimation in low pressure assemblages lacking garnet; and (3) it has been used to distinguish between clockwise and anticlockwise P–T paths in low pressure metamorphic settings. Experiments and thermodynamic databases provide conflicting constraints on the slope and position of the reaction, with most thermodynamic databases predicting a positive slope for the reaction. Evidence from mineral assemblages and microtextures from a large number of natural prograde sequences, in particular contact aureoles, is most consistent with a negative slope (andalusite and/or sillimanite occurs upgrade of, and may show evidence for replacement of, cordierite). Mineral compositional trends as a function of grade are variable but taken as a whole are more consistent with a negative slope than a positive slope. Thermodynamic modelling of reaction 1 and associated equilibria results in a low pressure metapelitic petrogenetic grid in the system K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O (KFMASH) which satisfies most of the natural and experimental constraints. Contouring of the Fe–Mg divariant interval represented by reaction 1 allows for pressure estimation in garnet‐absent andalusite+cordierite‐bearing schists and hornfelses. The revised topology of reaction 1 allows for improved analysis of P–T paths from mineral assemblage sequences and microtextures in the same rocks.     

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.     

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 for a morb composition in a P–T range, 450–650 °C and 18–28 kbar: the stability of eclogite

Lower temperature eclogite (with T = 600 °C) represents a significant part of the occurrences of eclogite in orogenic belts. ‘True’ eclogite, with, for example, garnet + omphacite >70%, is well represented in such an occurrence. Calculated phase equilibria in Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O (NCKFMASHTO), for just one rock composition – that of a representative mid‐ocean ridge basalt, morb – are used to see under what circumstances ‘true’ eclogite is predicted to occur. The variables considered are not only pressure (P) and temperature (T) but also water content and oxidation state. The latter two variables are known to exert a significant control on mineral assemblage but are difficult to establish retrospectively from the observed rocks themselves. It is found that whereas oxidation state does have a strong effect on mineral assemblage, the key control on developing ‘true’ eclogite is shown to be temperature and water content. If temperature is established to be <600 °C, water content has to be low (less or much less than that for H₂O saturation) in order for ‘true’ eclogite to form. Moreover, unless pressure is at the high end in the range considered, lawsonite eclogite and ‘true’ eclogite will tend to be mutually exclusive, with the former requiring high water content at the lower temperature where it occurs, but the latter requiring low water content.