Introduction
Dissolution–reprecipitation occurs when a fluid comes into contact with a mineral with which it is undersaturated, resulting in the dissolution of the less stable parent phase and reprecipitation of a more stable daughter phase (Putnis and Putnis, Reference Putnis and Putnis2007; Putnis, Reference Putnis, Carroll and Holloway2009; Ruiz-Agudo et al., Reference Ruiz-Agudo, Putnis and Putnis2014; Altree-Williams et al., Reference Altree-Williams, Pring, Ngothai and Brugger2015; Konrad-Schmolke et al., Reference Konrad-Schmolke, Halama, Wirth, Thomen, Klitscher, Morales, Schreiber and Wilke2018). Spatial coupling of dissolution and reprecipitation commonly occurs, where an interfacial fluid film at the reaction front becomes saturated resulting in precipitation of a stable daughter phase at the reaction surface of the unstable parent phase (Putnis, Reference Putnis, Carroll and Holloway2009; Ruiz-Agudo et al., Reference Ruiz-Agudo, Putnis and Putnis2014). Within regional metamorphic garnet, coupled dissolution–reprecipitation processes may be fingerprinted, producing a compositionally and texturally distinct garnet (Martin et al., Reference Martin, Ballèvre, Boulvais, Halfpenny, Vandergaeghe, Duchêne and Deloule2011; Dempster et al., Reference Dempster, La Piazza, Taylor, Beaudoin and Chung2017; Dempster et al., Reference Dempster, Gilmour and Chung2019). Garnet is characterised by sluggish volume diffusion at amphibolite-facies conditions and can record changing P–T–X (pressure, temperature, composition) during growth and post-growth modification (Jiang and Lasaga, Reference Jiang and Lasaga1990; Spear, Reference Spear1991; Raimondo et al., Reference Raimondo, Payne, Wade, Lanai, Clark and Hand2017).
This study examines the behaviour of zircon in a prograde regional metamorphic sequence of pelitic schists from the Scottish Highlands. Garnet in these schists has been studied previously (Dempster et al., Reference Dempster, La Piazza, Taylor, Beaudoin and Chung2017; Reference Dempster, Gilmour and Chung2019), and the garnet characteristics defined in these studies have been used as a context for understanding the reaction history. Primary clear garnet is inclusion-poor and represents essentially unmodified porphyroblasts, whereas cloudy garnet develops locally, and has been texturally and compositionally modified by coupled dissolution–reprecipitation (Martin et al., Reference Martin, Ballèvre, Boulvais, Halfpenny, Vandergaeghe, Duchêne and Deloule2011; Dempster et al., Reference Dempster, La Piazza, Taylor, Beaudoin and Chung2017). Such garnet contains large, irregularly-shaped quartz inclusions and abundant small, sub-rounded fluid inclusions (Martin et al., Reference Martin, Ballèvre, Boulvais, Halfpenny, Vandergaeghe, Duchêne and Deloule2011; Dempster et al., Reference Dempster, La Piazza, Taylor, Beaudoin and Chung2017). At higher temperatures, the cloudy texture of garnet may be eliminated, and the porphyroblast composition further modified producing secondary clear garnet (Dempster et al., Reference Dempster, Gilmour and Chung2019). Secondary clear garnet forms as rims at porphyroblast margins and surrounding mineral inclusions, these re-equilibrated secondary rims are more texturally and compositionally homogeneous than primary clear and cloudy garnet. Garnet is one of the most studied minerals in metamorphic petrology (Hollister, Reference Hollister1966; Yardley, Reference Yardley1977; Dempster, Reference Dempster1985; Hames and Menard, Reference Hames and Menard1993; Carlson, Reference Carlson2006; Baxter et al., Reference Baxter, Caddick, Dragovic, Kohn, Engi and Lanari2017), making it the perfect canvas to understand the behaviour of zircon during metamorphism.
Zircon is viewed traditionally as a physically and chemically unreactive mineral, highly resistant to metamorphism except at very high temperatures (Poldevaart, Reference Poldevaart1955; Gastil et al., Reference Gastil, Delisle and Morgan1967). A number of studies report evidence of zircon reacting across a range of metamorphic conditions (Hoskin and Black, Reference Hoskin and Black2000; Rubatto et al., Reference Rubatto, Williams and Buick2001; Williams, Reference Williams2001; Dempster et al., Reference Dempster, Hay and Bluck2004; Rasmussen, Reference Rasmussen2005; Rubatto and Hermann, Reference Rubatto and Hermann2007; Dempster et al., Reference Dempster, Hay, Gordon and Kelly2008; Hay and Dempster, Reference Hay and Dempster2009; Vonlanthen et al., Reference Vonlanthen, Fitzgerald, Rubatto and Hermann2012; Dempster and Chung, Reference Dempster and Chung2013; Wang et al., Reference Wang, Dunkley, Clarke and Daczko2014), raising questions about factors that promote zircon reactivity. Zircon growth during metamorphism can occur as outgrowths on detrital grains (e.g. Hay and Dempster, Reference Hay and Dempster2009) or as discrete very fine (<3 μm) zircon termed ‘micro-zircon’ (Dempster et al., Reference Dempster, Hay and Bluck2004; Dempster et al., Reference Dempster, Hay, Gordon and Kelly2008). In this work zircon populations are analysed from schists containing unmodified garnet, cloudy garnet modified during staurolite growth and re-equilibrated clear garnet formed during sillimanite growth to understand the effect each stage of the reaction history has on zircon behaviour.
Geological setting
This study analyses polymetamorphic Leven Schists from the Appin Group in Glen Roy, Scotland. Garnet- and staurolite-zone schists were from NN29861 85668 (UK Ordnance Survey Grid Reference), whereas sillimanite-zone schists were from the north side of Glen Roy NN2671 9289. The samples consist of fine-grained interbedded micaceous and quartzofeldspathic metasedimentary rocks with 1–4 mm garnet porphyroblasts (Dempster et al., Reference Dempster, La Piazza, Taylor, Beaudoin and Chung2017). Such matrix layers range from <0.2 mm up to 3 mm thick. The Appin Group rocks age is determined to be 656 ± 9 Ma, using Re–Os geochronology, representing the age of sedimentation (Rooney et al., Reference Rooney, Chew and Selby2011). The rocks have undergone complex tectonothemal evolution during the Grampian Orogeny in the Late Neoproterozoic and Early Cambrian (Phillips and Key, Reference Phillips and Key1992; Strachan et al., Reference Strachan, Harris, Fettes, Smith and Trewin2002; Vorhies et al., Reference Vorhies, Ague and Schmitt2013). The extent of Precambrian metamorphism is unknown (Dempster et al., Reference Dempster, Rogers, Tanner, Bluck, Muir, Redwood, Ireland and Paterson2002; Hutton and Alsop, Reference Hutton and Alsop2004). Peak metamorphism of the Leven Schists at Glen Roy is estimated at between 500–600°C and 5–8 kbar (Richardson and Powell, Reference Richardson and Powell1976; Wells, Reference Wells1979; Powell and Evans, Reference Powell and Evans1983; Phillips and Key, Reference Phillips and Key1992).
Methods
Analytical methods involved transmitted light microscopy (TLM) and scanning electron microscopy (SEM) of polished sections. TLM was used predominantly for initial petrographic analysis, to highlight the areas of interest for SEM, and when analysing inclusion alignment on multiple planes.
Energy dispersive spectroscopy analysis using the Carl-Zeiss Sigma VP electron microscope at University of Glasgow, operated at 20 kV was undertaken to analyse the composition of garnet. Porphyroblasts were mapped and compositional data extrapolated as line transects using AZtec 3.0 Oxford Instruments software to produce quantitative data. Following the removal of mineral inclusions from transects, running averages were used to smooth compositional profiles and reduce analytical noise. Back-scattered electron (BSE) imaging was used to capture garnet surface features including mineral inclusions, fluid inclusions and fractures. Mineral inclusion abundance in garnet was calculated in vol.% using java analysis on BSE images using ImageJ2, where colour thresholds were adjusted to capture the abundance of different phases.
Zircon abundance and morphology was characterised using an FEI Quanta 200F environmental SEM, at the University of Glasgow, to capture BSE images and map zircon populations. It was operated at 20 kV and moderate beam currents of 26 μA with a dwell time of 100 μs at 1500× magnification to image the smallest observable micro-zircon, ~10 nm wide. Brightness and contrast was adjusted to ensure micro-zircons were identifiable from other high mean atomic number features and spot analyses was used on every suspected zircon to ensure accurate identification. Transects and entire porphyroblasts of garnet were mapped using a series of frames ~100 μm wide that were montaged together on Inkscape producing a map of the analysed area. From this map quantitative data were extracted through Inkscape, producing data on zircon morphology and abundance. To assess micro-zircon abundance relative to fractures, a coordinate grid was produced over a BSE image map of garnet in which the location of all micro-zircons were identified and an equal number of randomly generated coordinates was produced and plotted. The distance from both populations of micro-zircons to the nearest fracture was then measured in μm, allowing a comparison of actual micro-zircon distribution versus randomly generated points.
Where average values are quoted e.g. average size of zircon within a domain, the standard deviation is to 1σ. Average size values produce large errors and a much larger sample size would be required to produce robust average sizes. Hence size distribution is discussed instead throughout this investigation. In total ~26 mm2 of the schists were analysed, 22.3 mm2 of garnet and 3.6 mm2 of the matrix. Within the staurolite-zone schists six porphyroblasts were imaged, three in staurolite-bearing schist GR01 and three in staurolite-absent schists GR02 and GR07. In the sillimanite-zone schists three garnet porphyroblasts were imaged, one in sillimanite-rich schist UGR0 and two in sillimanite-poor schist UGR1.
Results
Petrography
Staurolite zone
Three polished sections were selected for analysis on the basis of the varying proportions of staurolite and cloudy garnet (Table 1). The schists are composed of an interbedded quartzofeldspathic and micaceous matrix with garnet, biotite and/or staurolite porphyroblasts (GR01, GR02 and GR07, Table 1).
– = not detected
The micaceous layers display a strong S2 fabric (Fig. 1a). S1 formed early in the metamorphic history during the development of muscovite + biotite + chlorite, whereas S2 developed during the second stage of deformation accompanying the formation of biotite + muscovite + garnet ± staurolite (Phillips and Key, Reference Phillips and Key1992). The micaceous matrix is fine grained, <0.1 mm, with layers dominated by ~65 vol.% muscovite (Fig. 1a), with ~25 vol.% biotite, and ~5 vol.% quartz. Quartzofeldspathic layers are coarser, ~0.25 mm, and are dominated by ~70 vol.% quartz and ~20 vol.% plagioclase (Fig. 1a).
Garnet is typically ~3 mm in size, occasionally up to 6 mm, and euhedral-to-subhedral. It commonly has one or more irregular, embayed margins (Fig. 1b). Biotite occurs as a ~0.5 mm matrix phase and as larger 2–5 mm porphyroblasts, with retrograde chlorite replacing up to 30 vol.% of biotite (Fig. 1c). Staurolite is developed in a few of the garnet-bearing schists, immediately adjacent to apparently identical schists lacking staurolite (Dempster et al., Reference Dempster, La Piazza, Taylor, Beaudoin and Chung2017). It is typically <2 mm in size and forms adjacent to euhedral garnet margins (Fig. 1b), preferentially within a micaceous matrix (Dempster et al., Reference Dempster, La Piazza, Taylor, Beaudoin and Chung2017). As the proportion of staurolite increases, so too does the proportion of cloudy garnet (Table 1). During staurolite formation, chlorite is invoked typically as a reactant (Philips and Key, Reference Phillips and Key1992), however in the Glen Roy schists there is no evidence of primary chlorite (Dempster et al., Reference Dempster, La Piazza, Taylor, Beaudoin and Chung2017). Instead, in these schists, a small amount of H2O is invoked to balance the hydrous component in the following prograde reaction (Dempster et al., Reference Dempster, La Piazza, Taylor, Beaudoin and Chung2017):
Fe-rich primary clear garnet + muscovite + H2O → staurolite + biotite + quartz + cloudy garnet (Dempster et al., Reference Dempster, La Piazza, Taylor, Beaudoin and Chung2017).
Sillimanite zone
Two polished sections were analysed, a sillimanite-poor schist and a sillimanite-rich schist (UGR1, UGR0, Table 1). The matrix of these Leven Schists is also composed of interbedded quartzofeldspathic and micaceous layers with the latter containing the majority of the garnet, biotite and staurolite porphyroblasts. The matrix mineralogy is similar to the staurolite-zone schists, although the fabric is more strongly folded, probably by a third deformation event, S3 (Phillips and Key, Reference Phillips and Key1992) making original planar compositional layering less prominent (Fig. 1d). The schistose matrix is dominantly micaceous with <1 mm quartzofeldspathic layers.
Garnet porphyroblasts are 2–4 mm in size and subhedral, containing irregular and embayed margins. Most garnet porphyroblasts have at least one planar margin, and the proportion of modified cloudy garnet is higher than in the corresponding staurolite-zone schists (Table 1). Biotite occurs as both matrix and porphyroblast forms and replacement by chlorite is limited, with retrogression associated more typically in biotite included in garnet or at the core of atoll garnet structures. Staurolite occurs as two populations: ~3 mm subhedral–anhedral porphyroblasts forming at margins of garnet; and <1 mm euhedral-to-subhedral porphyroblasts in a micaceous matrix.
Sillimanite typically forms fibrolite mats, 0.25–1.5 mm in diameter, along the margins of garnet in the micaceous matrix (Fig. 1e). Sillimanite also forms mats included within garnet, generally <0.5 mm (Fig. 1f), and more rarely occurs as fracture fill close to embayed garnet margins. Sillimanite inclusions are limited to cloudy garnet and are more common in the sillimanite–rich schist, UGR0.
According to the KFMASH (K2O–FeO–MgO–Al2O3–SiO2–H2O) system, reactions involved in sillimanite formation typically involve the consumption of staurolite (Guidotti, Reference Guidotti1974; McLellan, Reference McLellan1985), the Glen Roy schists show no evidence of staurolite dissolution (Dempster et al., Reference Dempster, Gilmour and Chung2019). Sillimanite is thought to form as a result of the following garnet breakdown reaction:
Cloudy garnet + muscovite → sillimanite + biotite + quartz + secondary clear garnet (Dempster et al., Reference Dempster, Gilmour and Chung2019).
Characterisations of staurolite-zone garnet
Garnet in the staurolite-zone schists consists of unmodified primary clear garnet, cloudy garnet and transitional garnet. Cloudy garnet is characterised by the presence of small, aligned fluid inclusions and large, irregularly-shaped quartz inclusions (Dempster et al., Reference Dempster, La Piazza, Taylor, Beaudoin and Chung2017). Compositional profiles show disruption to normal garnet zoning associated with the cloudy garnet (Dempster et al., Reference Dempster, La Piazza, Taylor, Beaudoin and Chung2017).
Primary clear garnet
Primary clear garnet has mostly planar margins and few fractures. Such garnet contains few fluid inclusions, <0.5 vol.%, though it hosts a variety of mineral inclusions comprising between 6–27 vol.%. Quartz is the most abundant inclusion representing >60 vol.% of the population, followed by ~15 vol.% ilmenite, ~13 vol.% detrital zircon and <10 vol.% allanite. Quartz inclusions display regular margins, are rounded to tabular and are typically <50 μm in size. Ilmenite occurs as two distinct populations; large irregularly-shaped detrital grains (>50 μm) and small (<3 μm) euhedral grains. Allanite is sub-rounded and most inclusions are >30 μm in size with complex internal textures thought to be due to hydrothermal alteration, locally producing different rare earth element (REE) compositions (Smith et al., Reference Smith, Henderson and Jeffries2002; Gros et al., Reference Gros, Slaby, Jokubaukas, Sláma and Kozub-Budzyń2020). Zircon inclusions are discussed separately below.
Mineral inclusions form trails aligned with the matrix fabric, typically curved towards the porphyroblast rim (Fig. 2a) indicating syntectonic growth. Mineral inclusions form either quartz-rich or ilmenite-rich trails. The former can be traced into the quartzofeldspathic matrix and the latter to the micaceous matrix. There is no consistent change in mineral inclusion abundance or morphology from the garnet core to the rim.
Primary clear garnet has concentric divalent compositional zoning from core to rim (Fig. 2a) typical of growth zoning (Atherton, Reference Atherton1968; Dempster, Reference Dempster1985; Chakraborty and Ganguly, Reference Chakraborty, Ganguly and Ganguly1991; Gatewood et al., Reference Gatewood, Dragovic, Stowell, Baxter, Hirsch and Bloom2015). Manganese concentrations have a bell-shaped profile decreasing from ~15 mol.% spessartine in the core to <1 mol.% in the rim, whereas Fe and Mg display complementary bowl-shaped profiles. Calcium shows a less well-defined concentric pattern though generally mimics Mn, decreasing from core to rim (Fig. 2).
Cloudy garnet
Cloudy garnet forms bands in staurolite-zone garnet porphyroblasts (Fig. 3a). The extent of cloudiness in porphyroblasts varies, some are entirely clear, others contain ~75% cloudy garnet. Cloudy garnet contains a high concentration of 10s–100s μm-sized quartz inclusions and very abundant concentrations of fluid inclusions <3 μm (Fig. 3d). Quartz inclusions are irregularly shaped and display no alignment. Fluid inclusions are rounded to subrounded and trails of inclusions show strong planar alignment (see figure 5a in Dempster et al., Reference Dempster, La Piazza, Taylor, Beaudoin and Chung2017) (Fig. 4b), comprising 2–10 vol.% of cloudy garnet. Fluid inclusions show one dominant plane of alignment in garnet with another less well-defined plane perpendicular to it (Fig. 4). The dominant alignment is approximately parallel to σ1 of the S2 mica fabric and the secondary is parallel to σ3 (Fig. 4).
Mineral inclusions include zircon, monazite, quartz, rutile and xenotime. Monazite is the third most abundant mineral inclusion, representing 1.1 vol.% of cloudy garnet. Monazite inclusions can rarely reach up to 20 μm, but are more commonly 3–5 μm in size. Larger monazite crystals are typically hosted in quartz inclusions and fractures, whereas smaller inclusions are hosted by garnet and more rarely in fluid inclusions. Regarding minerals containing heavy and light rare earth elements (HREE, LREE), cloudy garnet contains no HREE-rich allanite, whereas primary clear garnet contains abundant HREE-rich allanite but no LREE-rich monazite. Cloudy garnet is darker in plane polarised light than clear garnet owing to the high abundance of fluid inclusions (Fig. 3b,c). Such garnet contains a high fracture density, displaying long, >200 μm, cross porphyroblast fractures and shorter radial fractures <100 μm, typically surrounding inclusions.
Within cloudy garnet, almandine and pyrope components are greater relative to adjacent primary clear garnet (+7–10 mol.% Fe, and +2–4 mol.% Mg, respectively). There are lower grossular and spessartine components in cloudy garnet than in primary clear garnet (–6–10 mol.% Ca, and –5–6 mol.% Mn (Fig. 2b), respectively). Manganese behaves less consistently than the other cations, occasionally increasing in the outer 100 μm at garnet rims.
Transitional garnet
Transitional garnet forms in areas between cloudy and primary clear garnet and is considered to represent a partially modified intermediate zone ~100 μm in width. Transitional garnet generally contains <1 vol.% fluid inclusions, whereas cloudy garnet contains >1 vol.% fluid inclusions and primary clear garnet contains significantly less, generally <0.1 vol.%. There is no change in the size distribution or morphology of fluid inclusions between transitional garnet and cloudy garnet. Transitional garnet contains few mineral inclusions, only rarely containing monazite, quartz and ilmenite. Where present, quartz inclusions in transitional garnet are small, aligned and regularly shaped, consistent with quartz inclusion morphology in primary clear garnet. The fracture density is also similar to that within primary clear garnet i.e. there is an absence of small-scale, interconnected fractures, which are typically abundant in cloudy garnet.
Whereas the texture of transitional garnet is more consistent with primary clear garnet, the composition of these zones is consistent with cloudy garnet (Fig. 2c). Transitional garnet contains high Mg and Fe, and low Mn and Ca contents. The compositional change at cloudy–clear margins is gradual, the distance over which the composition changes varies for each cation. Manganese shows the most gradual change over a distance of 50 μm, whereas Ca shows the sharpest change with composition adjusting over just 15 μm (Fig. 2c).
Characteristics of sillimanite-zone garnet
Sillimanite-zone garnet is composed of primary clear garnet, cloudy garnet and secondary clear garnet. Dempster et al. (Reference Dempster, Gilmour and Chung2019) constrained the textural and compositional characteristics of clear, cloudy and ‘texturally re-equilibrated garnet’ within the sillimanite-zone schists from Glen Roy; the latter is referred to as ‘secondary clear garnet’ in this study. Transitional garnet is not obvious in the sillimanite-zone schist owing to an uneven distribution of fluid inclusions due to the presence of secondary clear garnet, as such it is not described in the following section.
Primary clear garnet
The sillimanite-bearing schists contain much less primary clear garnet than within the staurolite schists (Table 1). Garnet within the sillimanite-rich schist UGR0 is composed of just 14 vol.% primary clear garnet whereas porphyroblasts in the sillimanite-poor UGR1-1 contain 36 vol.% primary clear garnet. Primary clear garnet is texturally and compositionally identical to that within the staurolite-zone schists.
Cloudy garnet
The geometry of cloudy zones is generally more regular in sillimanite-zone porphyroblasts than those in the staurolite zone, forming preferentially in the core of porphyroblasts. The texture of cloudy garnet is similar to that in the staurolite-zone schists, with abundant irregularly-shaped quartz inclusions and fluid inclusions. The abundance of the latter is variable, with a mixture of quartz-rich cloudy bands (Fig. 3f) and quartz-poor cloudy garnet (Fig. 3g). The quartz inclusion-rich domains contain up to 70 vol.% quartz (Fig. 3f), these quartz-rich bands are typically marginal and are present in ~50% of garnet porphyroblasts in these schists. The quartz-rich domains contain abundant inclusions, such as detrital zircon, monazite and rutile, and increased fracture abundance, predominantly small-scale <200 μm fractures terminating at adjacent inclusions. Fluid-inclusion abundance in quartz-rich domains is typically <1 vol.%. Cloudy garnet outside of quartz-rich bands contain <20 vol.% quartz inclusions. These quartz-poor areas contain a much higher abundance of fluid inclusions (Fig. 3g), more consistent with the staurolite-zone cloudy garnet.
Quartz inclusions are irregularly shaped, typical of cloudy garnet at staurolite-grade conditions, but are relatively small, generally <100 μm. Larger quartz inclusions ~500 μm are rare, monocrystalline and consistently contain sillimanite inclusions. In addition to quartz there are inclusions of ilmenite, rutile, zircon, monazite and sillimanite. Rutile inclusions are more abundant ~1.2 vol.% in sillimanite-zone garnet and the proportion of ilmenite ~1 vol.% is lower than within staurolite-zone garnet, ~1.5 vol.%.
Compositional transects of cloudy garnet are complicated by high inclusion and fracture abundance resulting in analyses of garnet compositions in transects difficult. Spot analyses reveal low Ca and high Mg contents, consistent with staurolite-zone schists.
Secondary clear garnet
Secondary clear garnet is exclusive to the sillimanite-zone schists. It is texturally similar to primary clear garnet containing <0.1 vol.% fluid inclusions. Secondary clear garnet occurs locally surrounding mineral inclusions and infilled fractures, and at most porphyroblast margins (Fig. 3h). Within the porphyroblast interiors, the development of secondary clear garnet is widest adjacent to sillimanite inclusions, locally reaching >500 μm. Here the garnet may contain a high fracture density and some irregularly-shaped inclusions, typical of cloudy garnet, though almost no fluid inclusions. In contrast, secondary clear garnet at the porphyroblast margins generally contains no mineral or fluid inclusions and very limited fractures (Fig. 3h). Such margins are typically 100–200 μm in width, reaching >1 mm adjacent to fibrolite mats in the matrix.
Secondary clear garnet has consistently low Mn and Ca contents, whereas Mg and Fe contents are high and more variable. Most notably, grossular contents in primary clear garnet are ~10 mol.%, decreasing to ~5 mol.% in cloudy garnet, and <5 mol.% in secondary clear garnet (Fig. 5).
Zircon abundance and morphology
Zircon occurs as two distinct populations, micro-zircon (Fig. 6a,b), and detrital zircon (Fig. 6d–f), distinguishable on the basis of their internal structure and size. Detrital zircon displays an internal structure in BSE images, indicative of growth over multiple events (Hoskin and Schaltegger, Reference Hoskin, Schaltegger, Hanchar and Hoskin2003) and are typically subhedral–anhedral. Micro-zircons are very fine grained, typically euhedral-to-subhedral and display no obvious internal structure (Dempster et al., Reference Dempster, Hay, Gordon and Kelly2008). A plot of all zircon sizes analysed produces a trimodal distribution (Fig. 7). There is a large peak for <1 μm2, which is interpreted to represent micro-zircon and a smaller peak for >3 μm2 representing detrital zircon, above this size the internal structure is visible in BSE imaging allowing distinction. Zircon in the 1–3 μm2 range cannot be classified easily, and such zircon is henceforth excluded from further analysis and interpretation, however attention was paid to ensure there were no systematic changes to the proportion of these unclassified grains. Within the analysed sections, 1967 zircons were characterised: 1352 micro-zircons, 351 detrital zircons and 282 unclassified zircons.
Zircon in staurolite-zone schists
Zircon in the matrix
The micaceous matrix in GR01 contains 60 micro-zircons per mm2 whereas the quartzofeldspathic matrix contains 27 micro-zircons per mm2. Zircon represents 5.2 vol.% of the matrix, both detrital and micro-zircon. Outgrowths (cf. Hay and Dempster, Reference Hay and Dempster2009) are rarely present on detrital grains, though are typically <1 μm (Fig. 6e). Most micro-zircons occur as inclusions within grains, 91%, whereas 71% of detrital zircons are located on grain boundaries. Within the micaceous matrix, >50 vol.% of the micro-zircons are hosted by muscovite and >60 vol.% in the quartzofeldspathic matrix are hosted in quartz. The amount of detrital zircon and micro-zircon within each phase is proportional to the modal mineralogy of the matrix, with the notable exception of retrograde chlorite, which contains detrital zircon but no micro-zircon.
Zircon in primary clear garnet
The matrix contains fewer micro-zircons than garnet, GR02 contains 29 micro-zircons per mm2 whereas adjacent garnet GR02-5 contains 72 micro-zircons per mm2. In contrast, garnet contains fewer detrital zircons than the matrix (Fig. 8), matrix GR02 contains 22 detrital zircons per mm2 whereas garnet GR02-5 contains 18 detrital zircons per mm2. Garnet GR02-1 that overgrew quartzofeldspathic matrix contains 30 micro-zircons per mm2, the same porphyroblast that overgrew micaceous matrix contains 73 micro-zircons per mm2. In the matrix there are 15–20% more micro-zircons in micaceous layers than in the adjacent quartzofeldspathic matrix, garnet that overgrows the former contains 25–50% more micro-zircons than garnet which overgrows the quartzofeldspathic matrix (Fig. 8).
Zircon shows no systematic variation in its distribution from the core to the rim of garnet, additionally there is no change in morphology with zircon maintaining a consistent size distribution. The amount of zircon within garnet porphyroblasts in the same polished section varies considerably, GR02-1 contains 46 micro-zircons per mm2 whereas GR02-5 contains almost double that amount, (72 micro-zircons per mm2). Garnet porphyroblasts which contain more detrital zircon also contain more micro-zircons (Fig. 8). Micro-zircons display a spatial link to detrital zircon and ilmenite, with most micro-zircons occurring <100 μm from detrital inclusions.
Zircon in cloudy garnet
Cloudy garnet contains on average 65 vol.% fewer micro-zircons than primary clear garnet in the same porphyroblast (Table 2). In GR01-9 primary clear garnet contains 96 micro-zircons per mm2 and cloudy garnet contains just 31 micro-zircons per mm2. Micro-zircons represent 27 μm2 per mm2 of primary clear garnet, whereas within cloudy garnet micro-zircons comprise 9 μm2 per mm2 (Table 2).
In cloudy garnet 62% of the micro-zircons are <0.3 μm2, whereas in primary clear garnet 45% are <0.3 μm2 in size (Fig. 9). Cloudy garnet contains a smaller proportion of fine detrital zircon <30 per mm2 and more >50 per mm2 zircon. Within cloudy garnet 35% of the micro-zircons are in fluid and mineral inclusions, whereas the remainder are hosted in garnet, compared to <1% of micro-zircon hosted by inclusions in primary clear garnet.
Zircon in sillimanite-zone schists
Zircon in the matrix
The micaceous matrix contains on average 32 micro-zircons per mm2 and 5 detrital zircons per mm2. The quartzofeldspathic matrix contains less zircon, 21 micro-zircons per mm2 and 4 detrital zircons per mm2. Zircon comprises <3% of the matrix of these schists. In the sillimanite-zone matrix 98% of the micro-zircon occurs as inclusions within single grains, whereas only 25% of the detrital zircons are hosted within single grains in the matrix. Quartz predominantly hosts micro-zircon in the quartzofeldspathic matrix and muscovite in the micaceous matrix. Chlorite contains no micro-zircon, although it is present in low modal abundances in these schists (Table 1). Most detrital zircons contain very thin outgrowths typically <1 μm, though occasionally reaching a few μm in width (Fig. 6e), and which probably formed in situ as the outgrowth morphology typically follows grain boundaries in the matrix.
Zircon in primary clear garnet
There is very limited primary clear, unmodified garnet in the sillimanite-zone schists. Primary clear garnet in UGR1-1 contains similar micro-zircon abundances to garnet which overgrew the quartzofeldspathic matrix in the staurolite-zone schists. However, compositional layering in the sillimanite-zone schists is less obvious and thus inclusion banding within garnet that overgrows the more finely interbedded sediments is less obvious. Micro-zircon and detrital zircon show a spatial link to fractures, with detrital zircons commonly occurring within fractures and micro-zircons occurring <10 μm from fractures in UGR1-1. In comparison a random distribution of inclusions produces a theoretical average distance of 200 μm between zircon inclusions and the nearest fracture.
Zircon in cloudy garnet
Cloudy garnet has a similar micro-zircon abundance to that in the staurolite-zone schists, much lower than within primary clear garnet (Table 2). Micro-zircons are especially rare in all sillimanite-zone schists owing to the abundance of cloudy garnet.
Micro-zircons occur predominantly in fluid inclusions, in porphyroblast UGR0-4 82% of micro-zircons are contained in fluid inclusions. Outgrowths on detrital zircon are relatively common in cloudy garnet, though they are generally <10 μm thick. Cloudy garnet in sillimanite-zone schists contain quartz-rich bands, these quartz-rich bands contain even fewer micro-zircons. Micro-zircons are more abundant in cloudy garnet immediately surrounding quartz-rich bands.
Zircon in secondary clear garnet
Secondary clear garnet that forms locally surrounding inclusions and at porphyroblast rims (Dempster et al., Reference Dempster, Gilmour and Chung2019), contains the fewest micro-zircons (15 per mm2) (Table 2). Secondary clear garnet at porphyroblast rims contains 9 vol.% less micro-zircon than secondary clear garnet surrounding mineral inclusions within garnet interiors (Fig. 3). Detrital zircon abundance in secondary clear garnet is consistent with that of cloudy garnet averaging 9 detrital zircons per mm2, only slightly lower than in primary clear garnet. There is no difference in micro-zircon morphology or abundance between cloudy garnet and secondary clear garnet.
Discussion
Staurolite-zone garnet
The variable nature of the inclusion trails within garnet porphyroblasts correlates with the nature of the original sedimentary layering. Quartz inclusion-rich primary clear garnet overgrew the quartzofeldspathic matrix and ilmenite inclusion-rich garnet overgrew the micaceous matrix. Primary clear garnet has undergone limited postgrowth modification, preserving the major divalent cation growth zoning, thus volume diffusion appears to be essentially ineffective at these metamorphic conditions (Hollister, Reference Hollister1966; Caddick et al., Reference Caddick, Konopásek and Thompson2010; Baxter et al., Reference Baxter, Caddick, Dragovic, Kohn, Engi and Lanari2017; Dempster et al., Reference Dempster, La Piazza, Taylor, Beaudoin and Chung2017).
Cloudy garnet forms as a result of coupled dissolution–reprecipitation during staurolite formation (Dempster et al., Reference Dempster, La Piazza, Taylor, Beaudoin and Chung2017). Following this, cloudy garnet contains ~38 vol.% more inclusions. Paired with increased fracture abundance this enables the transport of material away in the fluid phase. Coupled dissolution–reprecipitation may be focussed in quartz inclusion-rich garnet owing to the abundance of microcracks, formed radially around quartz, and grain boundaries (Gillet et al., Reference Gillet, Ingrin and Chopin1984; Engleder, Reference Engleder1987; Whitney et al., Reference Whitney, Cooke and Du Franc2000; Whitney et al., Reference Whitney, Goergen, Ketcham and Kunze2008). Quartz has a low bulk modulus (Wang et al., Reference Wang, Mao, Jiang and Duffy2015) compared to garnet (Leitner et al., Reference Leitner, Wiender and Liebermann1980; O'Neill et al., Reference O'Neill, Bass, Rossman, Geiger and Langer1991) resulting in fracturing during decompression. Microcracks can form interconnected networks, which act as fluid pathways promoting coupled dissolution–reprecipitation within garnet located in quartz-rich layers. Garnet that overgrows the micaceous matrix contains more abundant ilmenite, which has a higher bulk modulus (Wilson et al., Reference Wilson, Muscat, Mkhonto, Ngoene and Harrison2005; Cunha et al., Reference Cunha, Rodrigues, Rodrigues, Sampaio, Moulton, da Costa and Pizani2019) more similar to garnet, resulting in limited fracturing within garnet which has overgrown these quartz-poor layers.
Whereas some quartz inclusions in cloudy garnet might be remnants of quartz from the matrix, there is an increase in the mode of quartz within garnet following coupled dissolution–reprecipitation. The formation of quartz during this process might assist in the transport of alkalis between reactant and product phases (Konrad-Schmolke et al., Reference Konrad-Schmolke, Halama, Wirth, Thomen, Klitscher, Morales, Schreiber and Wilke2018). Dissolution–reprecipitation of pyroxene and amphibole has been suggested to involve the formation of a silica-rich amorphous phase (Keller et al., Reference Keller, Abart, Wirth, Schmid and Kunze2006; Bukovská et al., Reference Bukovská, Wirth and Morales2015; Konrad-Schmolke et al., Reference Konrad-Schmolke, Halama, Wirth, Thomen, Klitscher, Morales, Schreiber and Wilke2018). This amorphous phase allows the repolymerisation of product phases in addition to promoting fluid flow through the creation of grain boundaries (Konrad-Schmolke et al., Reference Konrad-Schmolke, Halama, Wirth, Thomen, Klitscher, Morales, Schreiber and Wilke2018). The quartz inclusions in this study are crystalline, however the complex tectonothermal history of the Glen Roy rocks could potentially produce a crystalline intermediate phase during coupled dissolution–reprecipitation of garnet.
Monazite within garnet is probably metamorphic, forming at the expense of allanite at the staurolite isograd (Wing et al., Reference Wing, Ferry and Harrison2003; Corrie and Kohn, Reference Corrie and Kohn2008; Gregory et al., Reference Gregory, Rubatto, Allen, Williams, Hermann and Ireland2007; Spear, Reference Spear2010). Allanite contains less Ce per mole than monazite (Ercit, Reference Ercit2002), as a result the volume of monazite in cloudy garnet is less. Heavy REE are more compatible with garnet (Rubatto, Reference Rubatto2002; Bui Thi Sinh et al., Reference Bui Thi Sinh, Osanai, Lenz, Nakano, Adachi, Belousova and Kitano2019); thus it may incorporate some HREE from allanite instead producing LREE-enriched monazite. The increased abundance and variety of mineral inclusions in cloudy garnet is further evidence of the ability of fluids involved in coupled dissolution–reprecipitation to transport high-field-strength elements.
Fluid inclusions form due to rapid growth in a fluid-rich environment (Crawford and Hollister, Reference Crawford and Hollister1986). The dominant alignment of fluid inclusions is approximately parallel to extensional fractures, whereas the secondary alignment is parallel to the matrix foliation. Both planes of alignment appear to be deformation-controlled, consistent with fracture orientation; fluid inclusions probably formed during the healing of fractures (Fig. 4). Unaligned inclusions may form asynchronously during coupled dissolution–reprecipitation, produced at an earlier stage when fluids are more available. Alternatively, evidence of alignment of fluid inclusions may be missing from samples where sections are cut >45° from the plane of alignment.
During coupled dissolution–reprecipitation, the garnet composition is modified as a result of partial re-equilibration with the fluid phase (Martin et al., Reference Martin, Ballèvre, Boulvais, Halfpenny, Vandergaeghe, Duchêne and Deloule2011; Ague and Axler, Reference Ague and Axler2016; Dempster et al., Reference Dempster, La Piazza, Taylor, Beaudoin and Chung2017, Reference Dempster, Gilmour and Chung2019). As coupled dissolution–reprecipitation typically forms sharp compositional boundaries (Putnis, Reference Putnis, Carroll and Holloway2009), the gradational boundary between cloudy and primary clear garnet might be indicative of later diffusive modification. Manganese undergoes rapid diffusion (Carlson, Reference Carlson2006) producing a gradual, sloping spessartine distribution profile between primary clear and cloudy garnet, whereas diffusion of Ca is slower (Carlson, Reference Carlson2006) resulting in a more abrupt grossular boundary. Cloudy garnet forms bands, which typically cross-cut compositional zoning, indicating that coupled dissolution–reprecipitation occurs irrespective of the original garnet composition and is instead focussed in fluid-rich areas, such as at crystal margins and along fractures. Most cloudy bands form along tensile planes parallel to the main extensional stress, fractures forming on this plane enable fluid access to the reactive parent garnet (Jamtveit and Austrheim, Reference Jamtveit and Austrheim2010; Dempster et al., Reference Dempster, La Piazza, Taylor, Beaudoin and Chung2017). Reaction fronts propagate parallel to these planes into unaltered parent garnet producing thick cloudy bands. The second common alignment of cloudy bands is parallel to the matrix fabric.
The consistent size distribution of fluid inclusions in cloudy and transitional garnet is indicative that transitional garnet formed during coupled dissolution–reprecipitation and is not a feature of textural re-equilibration. If fluid inclusions were a feature of later re-equilibration, then smaller fluid inclusions would be eradicated preferentially owing to their larger surface energy producing a population of larger fluid inclusions in transitional garnet. The development of coupled dissolution–reprecipitation in the Leven Schists is summarised in Fig. 10. The process begins at the margins of porphyroblasts and propagates towards the interior of the garnet through interconnected fractures. As the reaction interface moves into garnet interiors, fluid becomes less available primarily due to its consumption during reactions and entrapment as inclusions. As a result fewer fluid and mineral inclusions are formed in cloudy garnet with increasing distance from the original reaction interface, continuing until fluid availability falls below a threshold and the process terminates. The partial textural alteration in these zones where composition has been completely modified indicates that at staurolite-grade conditions, textural modification is less effective than compositional modification in garnet.
Sillimanite-zone garnet
In the sillimanite-zone schists the change in composition and texture in cloudy garnet mimics that within staurolite-zone garnet, meaning that all cloudy garnet probably formed synchronously during staurolite formation (Dempster et al., Reference Dempster, Gilmour and Chung2019). Products of fluid-mediated reactions exhibit higher defect densities in their lattice (Spruzeniece et al., Reference Spruzeniece, Piazolo and Maynard-Casely2017) as a result cloudy garnet is more reactive than primary clear garnet. Additionally, during coupled dissolution–reprecipitation Ca is released moving the sillimanite isograd to lower P–T space (Dempster et al., Reference Dempster, Gilmour and Chung2019), thus schists with a large proportion of cloudy garnet are more likely to form sillimanite.
Within the sillimanite-zone schists the original sediments were probably more finely interbedded. These layers produce more channelised structures which can act as conduits (Derry, Reference Derry, Kaufman and Jacobsen1992; Ague, Reference Ague2011). Estimates of fluid flow across layers are an order of magnitude slower than flow parallel to layering (Ferry and Gerdes, Reference Ferry and Gerdes1998). Fluid availability may therefore be higher in these schists, promoting coupled dissolution–reprecipitation and staurolite and cloudy garnet formation, and as a result sillimanite formation. Fluid availability is probably the limiting factor in the reaction process and in the less permeable staurolite-zone schists is limited, whereas in the more permeable sillimanite-zone schists coupled dissolution–reprecipitation and cloudy garnet formation is more effective. The quartz inclusion-rich cloudy garnet in sillimanite-bearing schists shares a similar morphology to cloudy zones in staurolite-zone schists, the quartz-rich bands may represent cloudy garnet forming from the initial propagation of fluid during coupled dissolution–reprecipitation at staurolite grade conditions. During the initial stages of the reaction fluid availability is higher (Putnis and Putnis, Reference Putnis and Putnis2007), increasing the potential for quartz formation and as the process progresses and temperature increases, mineral inclusion-poor cloudy garnet forms.
Secondary clear garnet is exclusive to sillimanite-zone schists, probably forming during sillimanite growth (Dempster et al., Reference Dempster, Gilmour and Chung2019). Secondary clear garnet has a distinctive low Ca composition, indicating it represents an additional stage in garnet evolution. Calcium contents can be used locally within a porphyroblast to monitor the degree of modification of garnet at higher grade where re-equilibration occurs making textural classification complex.
Secondary clear garnet only occurs around mineral inclusions, fractures and at porphyroblast rims, all of which are potentially fluid-rich areas. Fluid availability therefore could be an important factor in the re-equilibration of garnet during sillimanite formation. Both secondary clear garnet margins and secondary clear inclusion rims have the same low Ca, high Mg composition (Fig. 5), indicating that they formed under the same conditions. Porphyroblast margins are more homogeneous, mineral inclusions and fractures have been eliminated together with fluid inclusions, whereas mineral inclusion rims contain more original features with just fluid inclusions eradicated during re-equilibration. When dissolution and reprecipitation are closely coupled the product phase may be a pseudomorph of the dissolved parent phase, preserving original features through epitaxial nucleation (Putnis and Putnis, Reference Putnis and Putnis2007; Putnis and Austrheim, Reference Putnis and Austrheim2010; Ruiz-Agudo et al., Reference Ruiz-Agudo, Putnis and Putnis2014). Temporal decoupling may occur at porphyroblast margins where fluid availability is higher and thus dissolution is more effective than reprecipitation. The result is the eradication of original features and precipitation of a more homogeneous garnet product.
Zircon in staurolite-zone schists
The absence of micro-zircon on grain boundaries in the matrix is probably a facet of preservation. Grain boundaries have high fluid fluxes, whereas zircon contained within grains are sheltered. The abundance of micro-zircon as inclusions in all matrix phases suggests they probably formed at the same time as their host mineral in the matrix. The abundance of micro-zircon and detrital zircon within each matrix phase is proportional to their modal mineralogy, except for retrograde chlorite which contains a high proportion of detrital zircon and no micro-zircon. During prograde metamorphism, fluids are generated by dissolution of hydrous minerals, however during retrograde metamorphism fluid enters the rock from an external source (Etheridge et al., Reference Etheridge, Wall and Vernon1983; Connolly and Thompson, Reference Connolly and Thompson1989; Yardley et al., Reference Yardley, Gleeson, Bruce and Banks2000). The conduits exploited for fluid infiltration during retrogression may be exploited to remove fluid before it has the chance to become saturated in Zr and reprecipitate as zircon, resulting in chlorite with no micro-zircon.
The quartzofeldspathic matrix contains fewer micro-zircons than the micaceous matrix, this might be the result of lower fluid availability in quartzofeldspathic layers and or lower Zr fluxes for new zircon growth (Fig. 11). The micaceous matrix may produce halogen-enriched fluids during dehydration reactions, capable of preferentially dissolving and transporting Zr4+ (Rasmussen, Reference Rasmussen2005). Pelitic host rocks appear to contain a larger population of finer grained and potentially more metamict detrital zircon which could yield more Zr through dissolution (Dempster and Chung, Reference Dempster and Chung2013). In contrast the larger, more structurally intact detrital zircon within quartzofeldspathic host rocks are perceived to be less reactive. Additionally, the micaceous matrix contains abundant ilmenite, which can contain 1000s ppm Zr (Bingen et al., Reference Bingen, Austrheim and Whitehouse2001; Charlier et al., Reference Charlier, Skär, Korneliussen, Duchesne and Auwera2007; Bea et al., Reference Bea, Montero and Ortega2006), and could act as a Zr source. A combination of these three factors produces a micro-zircon-rich micaceous matrix and micro-zircon-poor quartzofeldspathic matrix (Fig. 11).
Garnet contains more micro-zircon than the matrix, indicative that during garnet growth micro-zircons are crystallising at the growing edge in addition to those captured from the matrix, which grew prior to garnet growth. Micro-zircon nucleating at the growing edge of garnet are probably encapsulated quickly resulting in a finer grained population than in the matrix. Micro-zircon formation is constant throughout porphyroblast growth as there is no change in micro-zircon morphology and abundance concentrically within the garnet. The matrix mineralogy controls the amount of Zr available for potential zircon growth, when garnet overgrows the micaceous matrix it accentuates the contrast of micro-zircon abundance across the micaceous and quartzofeldspathic matrix (Fig. 11). In the mica-rich matrix either: (a) ilmenite recrystallises at the garnet isograd releasing Mn for garnet growth (Woodsworth, Reference Woodsworth1977; Jiang et al., Reference Jiang, Palmer and Slack1996) and Zr for zircon; or (b) finer grained fragments of the more metamict detrital zircon are more reactive and dissolve further as fluid is released during garnet forming reactions (Dempster and Chung, Reference Dempster and Chung2013). Formation of abundant micro-zircon in garnet is therefore a result of micro-zircon encapsulation from the matrix but also formation of new micro-zircon where the Zr is derived from dissolution of detrital and micro-zircon at the garnet growing face and potentially the recrystallisation of ilmenite.
Cloudy garnet consistently contains much fewer micro-zircons than primary clear garnet. Coupled dissolution–reprecipitation preferentially occurs in garnet that overgrows the quartzofeldspathic matrix silicate inclusion-rich garnet, due to increased fracturing and hence fluid availability. Cloudy garnet contains fewer micro-zircon than primary clear garnet that overgrew the quartzofeldspathic matrix. The original lithology might in part account for a lower volume of micro-zircon in cloudy garnet however it is not solely responsible indicating that coupled dissolution–reprecipitation results in significant micro-zircon loss by dissolution with metamorphic grade. Cloudy garnet contains finer grained micro-zircon than primary clear garnet; such zircon would be the initial phase eradicated during dissolution owing to their higher surface energy. There is no disparity in micro-zircon size between garnet that overgrew the quartzofeldspathic and micaceous matrix, hence it is possible that during coupled dissolution–reprecipitation a new population of finer grained micro-zircon form at the expense of coarser micro-zircon in the original host primary clear garnet. As coupled dissolution–reprecipitation is a fast process (Putnis and Putnis, Reference Putnis and Putnis2007; Ruiz-Agudo et al., Reference Ruiz-Agudo, Putnis and Putnis2014), the reaction front propagates quickly encapsulating zircon in new garnet preventing further growth and producing finer grained micro-zircon.
Zircon in sillimanite zone
The decrease in the abundance of micro-zircon on grain boundaries in the matrix could be a facet of coarsening of the matrix, however it is also probably influenced by dissolution of micro-zircon at higher temperatures along grain boundaries, which act as fluid pathways.
Primary clear garnet in sillimanite-zone schists contains less zircon than primary clear garnet in staurolite-zone schists (Table 2) and the zircon is distributed unevenly, with micro-zircon and detrital zircon occurring in close proximity to fractures, typically <10 μm. Zircon has a higher bulk modulus, 223–250 GPa (Özkan and Jamieson, Reference Özkan and Jamieson1978; Hazen and Finger, Reference Hazen and Finger1979; Marqués et al., Reference Marqués, Florez, Recio, Gerward and Olsen2006), than garnet, ~170 GPa (Leitner et al., Reference Leitner, Wiender and Liebermann1980). Micro-zircons subsequently form close to detrital zircon during garnet growth indicative that detrital zircons are an important source of Zr. This close spatial relationship may also enable the formation of outgrowths in sillimanite-zone schists where micro-zircon nucleates close to detrital grains as transport ability is increased at higher temperatures.
Cloudy garnet in sillimanite-zone schists contain similar proportions of zircon to staurolite-zone cloudy garnet. Quartz inclusion-rich cloudy bands in garnet form in fluid-rich marginal areas of porphyroblasts and contain abundant fractures, it is unsurprising these quartz-rich areas of cloudy garnet contain fewer micro-zircon than in quartz-poor cloudy garnet owing to their greater dissolution potential. Zircon dissolved in quartz inclusion-rich cloudy garnet produces Zr-rich fluids, which are transported into adjacent quartz-poor cloudy garnet producing abundant micro-zircon in areas proximal to quartz-rich garnet. Micro-zircons are commonly hosted in fluid inclusions, indicating that they have formed in a fluid-rich environment and/or that micro-zircon hosted within inclusions are preferentially protected.
Secondary clear garnet has undergone an additional stage of modification associated with sillimanite formation and contains fewer micro-zircons. Micro-zircon dissolution seems particularly effective during the production of secondary clear garnet at porphyroblast margins, owing to increased fluid availability from the matrix, whereas inclusion rims in porphyroblast interiors are more sheltered from fluids. The consistent size distribution of micro-zircon between cloudy garnet and secondary clear garnet is indicative that no new micro-zircons are precipitated during sillimanite growth. It is instead possible a small proportion of micro-zircons formed in early events are preserved. Dissolution of detrital zircon is ineffective during sillimanite formation, secondary clear garnet contains similar detrital zircon abundances as cloudy garnet. At high temperatures, original metamict detrital zircon begins to undergo annealing producing more stable detrital zircon (Balan et al., Reference Balan, Neuville, Trocellier, Fritsch, Muller and Calas2001; Nasdala et al., Reference Nasdala, Wenzel, Varva, Irmer, Wenzel and Kober2001; Hay and Dempster, Reference Hay and Dempster2009) which is less likely to be dissolved, yielding limited Zr.
Sillimanite-zone schists contain more evidence of zircon outgrowths than within the staurolite-zone schists. The sillimanite-bearing schists have undergone further dehydration reactions producing abundant fluids, which when paired with the higher temperatures, make Zr transportation more effective. Increased transport rates and higher Zr equilibrium volumes promote the formation of outgrowths around existing detrital grains instead of micro-zircon neoblasts.
There is a cumulative 79% loss in the amount of micro-zircon following coupled dissolution–reprecipitation during staurolite formation and dissolution during sillimanite formation. Although it is currently not possible to do an overall assessment of the Zr budget in these rocks, the loss of micro-zircon certainly does not seem to be balanced by an increase in the abundance and size of zircon outgrowths on large detrital grains. Though there is an increase in the abundance of outgrowths within sillimanite-zone schists these are very thin and only present on some detrital zircon. Zirconium is being removed from zircon and potentially reincorporated into another phase for which it is not an essential component. Ilmenite and rutile have been shown to incorporate Zr in relatively large abundances, with rutile containing up to 1360 ppm Zr (Cave et al., Reference Cave, Stepanov, Craw, Large, Halpin and Thompson2015) and ilmenite Zr-contents generally less than rutile but proving slightly more difficult to characterise and highly dependent on silica content (Bea et al., Reference Bea, Montero and Ortega2006). A detectable amount of Zr is redistributed in these schists, if the zircon removed from the average 3 mm garnet following staurolite and sillimanite growth was incorporated solely into ilmenite and rutile, in the abundance they are present, it would represent an easily detectable amount of >1 wt.%. If it were incorporated within the garnet lattice it would represent 185 ppm Zr. It is important to note there is a grade dependence on how Zr is distributed. Sillimanite-zone schists contain evidence of overgrowths and rutile partially replaces ilmenite, potentially accounting for a large amount of Zr lost from micro-zircon dissolution. Whereas staurolite-zone schists have not undergone the ilmenite–rutile transition and contain almost no evidence of outgrowths, the distribution of Zr in these schists following coupled dissolution–reprecipitation is more uncertain.
Conclusions
Analysis of garnet and zircon behaviour during prograde metamorphism reveals important controls on coupled dissolution–reprecipitation processes and the importance of fluid availability on progressive metamorphism.
Coupled dissolution–reprecipitation of garnet is controlled initially by lithological layering and deformation. Garnet porphyroblasts that overgrow more quartzofeldspathic matrix layers generate quartz dominated inclusion suites, subsequently producing more fractures or conduits within the garnet and increasing fluid availability. Once access to the core of garnet is secured porphyroblast composition appears to be the dominant factor controlling metamorphic reaction progress, preferentially dissolving the Mn-rich core producing atoll garnet.
As metamorphic grade increases, dissolution and reprecipitation appear to become decoupled owing to faster dissolution rates. The result is the precipitation of a more homogeneous product phase with the removal of original features, as is the case with secondary clear garnet.
The mineralogy of the matrix impacts the micro-zircon abundance by influencing original detrital zircon populations and fluid availability, and therefore potential Zr fluxes.
Micro-zircon abundance increases during initial garnet formation though decreases with every stage of fluid-mediated modification in the reaction history. A detectable quantity of Zr is removed from zircon and is probably incorporated into other mineral phases – probably rutile and ilmenite, and potentially garnet, owing to lattice relaxation at higher temperatures.
Detrital zircon abundance decreases during garnet formation but remains consistent during any subsequent modification, potentially a result of increased stability caused by annealing of metamict grains at higher temperatures.
Retrogression is effective at dissolving micro-zircon however Zr is transported out of the system before new Zr can precipitate producing chlorite with no micro-zircon.
Much of the metamorphic history of these schists is influenced strongly by fluid availability and controlled by fine scale lithological variation in the matrix promoting or inhibiting channelisation and fluid availability. It has even been suggested the staurolite- and sillimanite-zone schists at Glen Roy have undergone identical P–T–X histories however fluid availability promotes or inhibits the production of index minerals (Dempster et al, Reference Dempster, Gilmour and Chung2019).
Acknowledgements
This work was carried out at the University of Glasgow. The authors thank Robert MacDonald for his assistance with sample preparation. Special thanks are given to Dr Thomas Mueller, Dr Freya George and Dr Roger Mitchell for their thorough and thoughtful reviews, which helped improve and refine this work.
Supplementary material
To view supplementary material for this article, please visit https://doi.org/10.1180/mgm.2021.97