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.     

Three metamorphic events recorded in a single garnet: Integrated phase modelling,in situ LA‐ICPMS and SIMS geochronology from the Moine Supergroup,NW Scotland

In situ LA‐ICP‐MS monazite geochronology from a garnet‐bearing diatexite within the Moine Supergroup (Glenfinnan Group) NW Scotland records three temporally distinct metamorphic events within a single garnet porphyroblast. The initial growth of garnet occurred in the interval c. 825–780 Ma, as recorded by monazite inclusions located in the garnet core. Modelled P–T conditions based on the preserved garnet core composition indicate an initially comparatively high geothermal gradient regime and peak conditions of ～650 °C and 7 kbar. Monazite within a compositionally distinct second shell of garnet has an age of 724 ± 6 Ma. This is indistinguishable from a SIMS age of 725 ± 4 Ma obtained from metamorphic zircon in the sample, which is interpreted to record the timing of migmatization. This second stage of garnet growth occurred on a P–T path from ～6 kbar and 650 °C rising to ～9 kbar and 700 °C, with the peak conditions associated with partial melting. A third garnet zone which forms the rim contains monazite with an age of 464 ± 3 Ma. Monazite in the surrounding matrix has an age of 462 ± 2 Ma. This corresponds well with a U–Pb SIMS zircon age of 463 ± 4 Ma obtained from a deformed pegmatite that was emplaced during widespread folding and reworking of the migmatite fabric. The P–T conditions associated with the final phase of garnet growth were ～7 kbar and 650 °C. The monazite ages coupled with the phase relations modelled from this multistage garnet indicate that it records two Neoproterozoic tectonothermal events as well as the widespread Ordovician Grampian event associated with Caledonian orogenesis. Thus, this single garnet records much of the Neoproterozoic to Ordovician thermal history in NW Scotland, and highlights the long history of porphyroblast growth that can be revealed by in situ isotopic dating and associated P–T modelling. This approach has the potential to reveal much of the thermal architecture of Neoproterozoic events within the Moine Supergroup, despite intense Caledonian reworking, if suitable textural and mineralogical relationships can be indentified elsewhere.     

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.     

Exhumation of the Dayman dome metamorphic core complex, eastern Papua New Guinea

The ～750 km² Dayman dome of the Late Cretaceous Suckling‐Dayman massif, eastern Papua New Guinea, is a domed landform that rises to an elevation of 2850 m. The northern edge of the dome is a fault scarp >1000 m high that is now part of an active microplate boundary separating continental crust of the New Guinea highlands from continental and oceanic crust of the Woodlark microplate. Previous work has shown that a parallel belt of eclogite‐bearing core complexes north‐east of the Dayman dome were exhumed from up to 24–28 kbar in the last few millions of years. The remarkably fresh and lightly eroded scarp of the Dayman dome exposes shallowly‐dipping mylonitic (S1) metabasite rocks (500 m thick) on the northern flank of Mount Dayman. Field relationships near the base of this scarp show a cross cutting suite of ductile and brittle meso‐structures that includes: (i) rare ductile S2 folia with a shallowly ESE‐plunging mineral elongation lineation defined by sodic‐calcic blue amphibole; (ii) narrow steeply‐dipping ductile D2 shear zones; and (iii) semi‐brittle to brittle fault zones. Pumpellyite‐actinolite facies assemblages reported by previous workers to contain local aragonite, lawsonite and/or glaucophane are found in the core of the complex at elevations greater than 2000 m. These assemblages indicate peak metamorphic pressures of 6–9.5 kbar, demonstrating exhumation of the core of the Dayman dome from depths of 20–30 km. The S1 metamorphic mineral assemblage in metabasite includes actinolite‐chlorite‐epidote‐albite‐quartz‐calcite‐titanite, indicative of greenschist facies conditions for the main deformation. New mineral equilibria modelling suggests that this S1 assemblage evolved at 5.9–7.2 kbar at ～425 °C. Modelling variable Fe³⁺ indicates that the sodic‐calcic blue amphibole (D2) formed under a higher oxidation state compared with the S1 assemblage, probably at <4.5 kbar. A SE‐dipping, Mio‐Pliocene sedimentary sequence (Gwoira Conglomerate) forms a hangingwall block juxtaposed by low‐angle fault contact with the metabasite footwall. Prehnite‐bearing D3 brittle fault zones separate the two blocks and likely accommodated the final exhumation of the S1 greenschist facies assemblage in the footwall. These results indicate that the extensive Mt Dayman fault surface coincides with a domed S1 greenschist facies foliation that was last active at >20 km depth. Exhumation of this foliation must therefore be controlled by brittle faults of the active microplate boundary that are largely not observed in the study area. The structural record of the final exhumation of the Dayman dome to the surface was likely lost as a result of erosion, poor exposure or wide spacing of semi‐brittle to brittle fault zones.     

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.     

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.     

The reaction history of kyanite in high‐P aluminous granulites

Cathodoluminescence (CL) mapping of kyanite in high pressure, aluminous granulites from the central Grenville Province reveals internal structures that are linked to their metamorphic reaction history. In two samples, individual kyanite crystals are shown to be composite porphyroblasts comprising three distinct generations, defined by their CL intensity and Cr (±V, Ti, Fe and Ga) content, and each separated by resorbed interfaces. In contrast, a sub‐aluminous sample contains two types of kyanite, one as resorbed inclusions in garnet and another in the groundmass or replacing garnet. These textural variants of kyanite are interpreted within the framework of phase equilibria modelling. In P–T pseudosections, a first generation of kyanite, which is only present in the most aluminous samples, is potentially linked to staurolite breakdown, and its resorption is consistent with a subsequent increase in pressure. This kyanite represents the earliest remnant of prograde metamorphism identifiable in these rocks. The second generation, present in the porphyroblasts in the same samples and as inclusions in garnet in the sub‐aluminous sample, is interpreted to be the peritectic product of muscovite dehydration melting. Resorption of this kyanite is consistent with subsequent continuous dehydration melting of biotite, which is also inferred based on microstructural considerations. The final generation of kyanite, present as rims on the prograde kyanite porphyroblasts in aluminous samples and as part of the groundmass or replacing garnet in the sub‐aluminous rock, is interpreted to have grown during melt crystallization upon retrogression. The presence of retrograde kyanite implies that the melt crystallized over a wide range of temperatures, and provides an important constraint on the P–T conditions of the metamorphic peak and on the retrograde P–T path. CL mapping is crucial for identifying retrograde kyanite in aluminous samples, as it preferentially overgrows existing kyanite rather than replacing other prograde phases. The scarcity of kyanite in sub‐aluminous rocks allows retrograde kyanite to grow as discrete crystals that can be identified by optical microscopy. This work attests to the potential of unconventional tools such as CL imaging for deciphering the metamorphic history of rocks.     

On the roles of deformation and fluid during rejuvenation of a polymetamorphic terrane: inferences on the geodynamic evolution of the Ruker Province, East Antarctica

Evaluating pressure–temperature (P–T) conditions through mineral equilibria modelling within an amphibolite facies polymetamorphic terrane requires knowledge of the fluid content of the rocks. The Archean‐Palaeoproterozoic basement rocks of the Ruker Province, East Antarctica, preserve evidence of three metamorphic events (M1–M3). Of particular interest is the M3 event, which is constrained to the early Palaeozoic (c. 550–480 Ma). Evaluation of the tectonic setting during this time is important because the Ruker Province is located within a critical region with respect to models of Gondwana assembly. Structural evidence of the early Palaeozoic event is preserved as large (up to ～500 m wide) high strain zones that cut the orthogneiss‐metasedimentary basement (Tingey Complex) of the Ruker Province. Rocks within these zones have been thoroughly recrystallized and preserve a dominant shear fabric and M3 mineral assemblages that formed at P–T conditions of 4.0–5.2 kbar and 565–640 °C. Distal to these zones, rocks preserve more complex petrographic relationships with S1 and S2 foliations, being incompletely overgrown by M3 retrograde assemblages. We show that the mineral assemblages preserved during the M3 event are highly dependent on the availability of fluid H₂O, which is strongly influenced by the structural setting (i.e. proximity to the high‐strain zones). P–T structural and fluid flow constraints support a model of basin inversion during early Palaeozoic crustal rejuvenation in the Ruker Province.     

Metastable persistence of pelitic metamorphic assemblages at the root of a Cretaceous magmatic arc – Fiordland, New Zealand

Four aluminosilicate-bearing, amphibolite facies pelitic schists sampled from the root of the long-lived eastern Gondwana continental magmatic arc now exposed in southwest Fiordland, New Zealand, record remarkably different P–T–t histories. The four samples were collected from within 20 km of each other within the Fanny Bay Group and Deep Cove Gneiss near Dusky Sound. Integrated petrography, mineral chemistry, mineral equilibria modelling and in situ electron microprobe chemical dating of monazite shows that the sample of the Fanny Bay Group south of the Dusky Fault records a Carboniferous history with peak conditions of 4–4.5 kbar at 570–590 °C, while one sample of the Deep Cove Gneiss from Long Island records a Cretaceous history with apparent peak conditions of 7.5 kbar at 650 °C. Two other samples of the Deep Cove Gneiss from Resolution Island record mixed Carboniferous and Cretaceous histories with apparent peak conditions of 7 kbar at 650 °C and 3–7 kbar at 640–720 °C. The metapelitic schists on Resolution Island were intruded by arc magmas including the voluminous high- P Western Fiordland Orthogneiss, yet they lack mineralogical evidence of the Cretaceous high- P (>12 kbar) event. Analysis of water isopleths in a model system shows that the amount of water accommodated in the rock mineral assemblage increases with pressure. With the exhaustion of all free water, and without the addition of external water, these rocks persisted metastably within the deep arc during the high- P event. The emplacement of large volumes of diorite (i.e. the Western Fiordland Orthogneiss) into the root of the Early Cretaceous continental magmatic arc did not lead to regional granulite facies metamorphism of the country rock schists, as large volumes of amphibolite facies rock metamorphosed under medium- P conditions persisted metastably in the deep arc crust.     

Retrograde melt–residue interaction and the formation of near‐anhydrous leucosomes in migmatites

Considering physical segregation of melt from its residue, the chemical potentials of the components (oxides) are the same in both when segregation occurs. Then, as P–T conditions change, gradients in chemical potential are established between the melt‐rich domains and residue permitting diffusional interaction to occur. In particular, on cooling, the chemical potential of H₂O becomes higher in the melt segregation than in the residue, particularly when biotite becomes stable in the residue assemblage. Diffusion of water from the melt to the residue promotes crystallization of anhydrous products from the melt and hydrous products in the residue. This diffusive process, when coupled with melt loss from the rocks subsequent to some degree of crystallization, can result in a significant degree of anhydrous leucosome being preserved in a migmatite with only minor retrogression of the residue. If H₂O can diffuse between the melt segregation and all of the residue, then no apparent selvedge between the two will be observed. Alternatively, if H₂O can diffuse between the melt segregation and only part of the residue, then a distinct selvedge may be produced. Diffusion of H₂O into the residue may be in part responsible for the commonly anhydrous nature of leucosomes, especially in granulite facies migmatites. Diffusion of other relatively mobile species such as Na₂O and K₂O has a lesser effect on overall melt crystallization but can change the proportion of quartz, plagioclase and K‐feldspar in the resultant leucosome. The diffusion of H₂O out of the melt results in the enhanced crystallization of the melt in the segregation and increases the amount of resulting anhydrous leucosome relative to the amount produced if melt crystallized in chemical isolation from the residue. For high residue:melt ratios, the proportion of resulting near‐anhydrous leucosome can approach that of the proportion of melt present at the onset of cooling with only minor loss of melt from a given segregation required. Crystallization of melt segregations via the diffusion of H₂O out of them into the host may also play a major role in driving melt‐rich segregations across key rheological transitions that would allow the expulsion of remaining melt from the system.