Metamorphic titanite–zircon pseudomorphs after igneous zirconolite

. The formation of metamorphic zircon after baddeleyite is a well-known reaction that can be used to date the metamorphism of igneous silica-undersaturated rocks. By contrast, metamorphic minerals formed after igneous zirconolite have rarely been reported. In this paper, we document metamorphic titanite + zircon pseudomorphs formed from the metamorphic breakdown of igneous zirconolite in syenodiorite and syenite, in the southeastern Sveconorwegian Province, Sweden. Water-rich ﬂuid inﬂux during tectonometamorphism in epidote–amphibolite-facies metamorphic conditions caused the release of silica during a metamorphic reaction involving igneous feldspar and pyroxene and the simultaneous breakdown of igneous Zr-bearing phases. Typical titanite + zircon intergrowths are elongated or platy titanite crystals speckled with tiny inclusions of zircon. Most intergrowths are smaller than 15 µm; some are subrounded in shape. Locally, bead-like grains of titanite and zircon are intergrown with silicate minerals. The precursor igneous zirconolite was found preserved only in a sample of near-pristine igneous syenodiorite, as remnant grains of mainly < 2 µm in size. Two somewhat larger crystals, 8 and 12 µm, allowed semiquantitative conﬁrmation using microprobe analysis. Analogous with zircon pseudomorphs after baddeleyite, titanite + zircon pseudomorphs after zirconolite potentially offer dating of the metamorphic reaction, although the small size of the crystals makes dating with today’s techniques challenging. The scarcity of


Introduction
In continental crust, metamorphic recrystallisation of H 2 Opoor igneous rocks is commonly steered by the influx of hydrous fluid, often guided by deformation (Austrheim, 1987;Engvik et al., 2000Engvik et al., , 2001;;Menegon et al., 2017).Metamorphism is, therefore, typically heterogeneously developed even in the deep parts of orogens (e.g.Austrheim, 2013;Jamtveit et al., 2018;Möller and Andersson, 2018;and, in this special issue, Gyomlai et al., 2023).The unreacted (metastable) igneous rocks and gradual transitions to metamorphosed equivalents allow insights into the progress of metamorphic reactions.In some silica-undersaturated gabbroic rocks, the metamorphic reaction and, by implication, the influx of hydrous fluid and associated deformation can also be dated.This is, for example, possible for metagabbro in which the silica released from the metamorphic reactions has led to the replacement of baddeleyite (ZrO 2 ) by meta-morphic zircon (ZrSiO 4 ; Davidson and van Breemen, 1988;Beckman and Möller, 2018).
While zircon growth after baddeleyite is a comparably well-known reaction (Davidson and van Breemen, 1988), there are to this date only a few reports (Pan, 1997;Pohlner et al., 2020) of metamorphic mineral growth after zirconolite, CaZrTi 2 O 7 .This may be due to zirconolite and pseudomorphs after zirconolite being either extremely rare or, possibly, overlooked.This paper documents, using electron and cathodoluminescence imaging, metamorphic titanite + zircon intergrowths formed from the metamorphic breakdown of igneous zirconolite in metamorphosed varieties of syenodiorite and syenite from the Sveconorwegian Province in southwestern Sweden.Zirconolite was found in the best-preserved igneous syenodiorite only, as remnant igneous grains of mainly very small size (< 2 µm); two larger crystals, 7 and 12 µm large, were analysed by wavelength-The crystal structure of zirconolite is described as a layer structure analogous to the pyrochlore (Mazzi and Munno, 1983;Gieré et al., 1998) and to an anion-deficient fluorite structure type (Pyatenko and Pudovkina, 1964;Gieré et al., 1998).The zirconolite structure consists of alternate layers of TiO 6 polyhedra aligned parallel to (001) with Zr 4+  and Ca 2+ arranged in planes ( 010) and (100), respectively.The complexity of the structure increases with the incorporation of trivalent REEs (e.g.La 3+ and Nd 3+ ) and tetravalent elements U and Th (Nickolsky and Yudintsev, 2021;Rasmussen and Fletcher, 2004).Furthermore, zirconolite may crystallise in three different systems; the varieties are known as polytypes and include zirconolite-3O (three-layered orthorhombic polytype), zirconolite-3T (three-layered trigonal polytype), and zirconolite-2M (two-layered monoclinic polytype).Zirconolite without a suffix denotes non-crystalline (metamict) zirconolite.The cubic variety is named zirkelite (Bayliss et al., 1989;Nickolsky and Yudintsev, 2021;Rasmussen and Fletcher, 2004).

Analytical methods
Microtextural documentation and mineral identification were performed on polished thin sections by using a Tescan MIRA3 field emission scanning electron microscope (FE-SEM) coupled with an Oxford Instruments energydispersive X-ray spectroscopy (EDS) analysis system and a cathodoluminescence (CL) imaging system at the Department of Geology, Lund University.Operation conditions For the mineral chemistry of zirconolite (two ≤ 12 µm large crystals), quantitative analysis was performed in situ in a thin section using the JEOL JXA-8230 electron probe microanalyser (EPMA) at the Universidad Nacional de Colombia.Major and trace elements were collected with a 1 µm diameter beam, 20 nA current, acceleration voltage at 15 kV, and 20 s collection time on the peak with background counting times of 10 s for Al, Nb, Lu, Gd, Eu, Pr, and Ce and 5 s for all other elements.The acquisition time was restricted to minimise the beam damage effects on the small zirconolite crystals, which limited the analytical resolution.Detection limits for Y + REE ranged from 300 to 1000 ppm.The Xray lines and standards used during the WDS analysis counts are summarised in Table 2. Several trial-and-error tests were made to define the interferences at which the analysis gave the best fit.Peak overlap corrections were applied for the interferences Eu-Pr and Er-Tb.For correction of background interferences, the spectrometer was set to measure count rates on both sides of the peak and on one side for Al, Nb, Lu, Gd, Eu, Pr, and Ce.The respective counts were converted to weight percentage oxides applying the CitZAF correction method.Fe 3+ was calculated by stoichiometry following the general equation by Droop (1987).
For the purpose of determining the crystallographic properties, we attempted to analyse the largest (12 µm) zirconolite crystal using electron backscattered diffraction (EBSD).However, the results of the EBSD analysis were not conclusive (data provided in Supplement).

Geological context
The rocks investigated in this study are part of the Vaggeryd syenite complex in southern Sweden.The Vaggeryd syen-ite complex is a N-S elongated body about 50 km long and 7 km wide, which intruded at 1.22 Ga (Pb-Pb ages from single zircon evaporation; Ask, 1996; Fig. 1).The main rock is syenite, but other silica-poor rocks, e.g.quartz syenite, mangerite, charnockite, and leucogabbro, also occur.The Vaggeryd syenite complex is located within the easternmost part of the Eastern Segment of the Sveconorwegian orogen, the remnant of a 1 Ga, large, and hot orogen at the southwestern margin of the Fennoscandian Shield (Fig. 1a).The Vaggeryd syenite complex is heterogeneously deformed by numerous N-S-striking ductile shear zones (Söderlund and Ask, 2006) that have also deformed the surrounding host rocks (Fig. 1b).Similarly aged granitic and syenitoid intrusions are also present further west in the Sveconorwegian Province, within the Eastern Segment (Fig. 1a), in areas that were metamorphosed under higher-grade conditions.
The Eastern Segment records the latest Sveconorwegian collisional stage at 990-960 Ma, during which it was overthrust by Sveconorwegia (Fig. 1).This event resulted in westward-increasing metamorphic grade from epidoteamphibolite facies to high-pressure granulite facies within the Eastern Segment and the exhumation of an eclogite terrane beneath the contact with Sveconorwegia.The easternmost boundary of the Sveconorwegian Province, in which the Vaggeryd syenite complex is situated, is a ca. 25 km wide subvertical N-S-trending deformation zone, progressively changing eastwards into discrete steeply west-dipping shears (the Frontal wedge in Fig. 1a, b; Wahlgren et al., 1994;Möller and Andersson, 2018).The steep structures in the Frontal wedge mark the eastward-waning ductile deformation and recrystallisation in the deeply buried Eastern Segment during the Sveconorwegian orogeny.Lithologies are continuous across this zone, which shows that the Frontal wedge is an intraplate deformation structure (Söderlund and Ask, 2006;Möller et al., 2015;Möller and Andersson, 2018).The ductile deformation within the Frontal wedge was nonpenetrative and associated with metamorphic recrystallisation under epidote-amphibolite-to greenschist-facies conditions.In many places, rocks remain in a pristine igneous or pre-Sveconorwegian state (Möller and Andersson, 2018).https://doi.org/10.5194/ejm-35-773-2023

Sample site
In a quarry at Näsbyholm in the southernmost part of the Vaggeryd intrusion (Fig. 1b), syenitoid rocks are wellexposed in outcrops along a ca.250 m E-W section.The main rock has local compositional variations spanning from reddish feldspar-rich syenite to dark melanocratic syenodiorite (probably incomplete mixing between syenitic and gabbroic melts).This main syenitoid complex also hosts metrewide diabase and pegmatite dykes.All rock types have been heterogeneously overprinted by deformation and associated metamorphism with strain variations over distances of less than a metre.Locally, undeformed rocks are medium-to coarse-grained with crystal sizes up to 1 cm, preserving their igneous fabrics.Where deformed, the rocks are foliated and recrystallised into metamorphic minerals of fine grain size (< 0.5 mm).Although the overall orientation of the foliation is interpreted to align with the regional N-S subvertical trend, the foliation anastomoses in different directions around lens-shaped domains on the scale of a few metres (Fig. 1c).

Sampling
Samples were taken within ca. 100 m of the same quarry outcrop.A well-preserved coarse-grained melanocratic syenodiorite was found in one spot only, in a domain a few metres large, surrounded by recrystallised syenodiorite transi-tioning from weakly deformed to mylonite.Figure 2a and b illustrate the difference between a well-preserved melanocratic syenodiorite, in this paper referred to as near-pristine igneous syenodiorite, and a weakly to moderately deformed and recrystallised variety of the same rock.The samples investigated in this study are mainly of syenodioritic composition (the bulk rock chemistry for three samples is provided in the Supplement, Table S1), but one sample is more reddish, dominated by alkali feldspar, and of syenitic composition.The syenodioritic samples include variations in the degree of strain and associated metamorphic recrystallisation; the reddish syenitic sample is strongly deformed.Five thin sections of the melanocratic syenodiorite, as well as one deformed alkali-feldspar-dominated syenite (09h), were investigated for Zr-bearing accessory minerals.Of these, one sample is near-pristine igneous syenodiorite (09j), one is a metamorphosed and undeformed variety (09a), and the four remaining samples are metamorphosed and deformed (09b, 09c, 09g, and 09h).

Petrography
The near-pristine igneous syenodiorite is dominated by darkgrey plagioclase crystals mantled by mesoperthite and an interstitial aggregate of dark minerals consisting of orthopyrox- ene crystals and lesser amounts of opaques, clinopyroxene, and biotite (Fig. 2a, c).In the deformed varieties, the feldspar occurs as augen or porphyroclasts, wrapped by biotite-rich folia, and pyroxenes are absent (Fig. 2b, d).
Igneous syenodiorite.In the near-pristine igneous syenodiorite, a mesoperthitic mantle (Fig. 2e-f) surrounds coarse plagioclase crystals.These feldspars are oligoclase-andesine (An 32 to An 34 , with an orthoclase component of up to 1 %) and alkali feldspar.The composition of the exsolved alkalifeldspar phase in mesoperthite is Or 68−91 Ab 9−27 .The aggregates of ferromagnesian phases (Fig. 2e-f) consist of coarse-grained orthopyroxene with exsolution lamellae of calcic pyroxene (the latter 15 vol %) and lesser amounts of bi-otite, apatite, and Fe-Ti oxides.Accessory minerals include titanomagnetite (magnetite with ilmenite lamellae), ilmenite, sulfides, baddeleyite, zircon, and zirconolite (detailed below).A few grains of srilankite were also found as minute inclusions in feldspar and baddeleyite grains.Small crystals of metamorphic minerals were also identified: tiny epidoteclinozoisite inclusions (∼ 50 µm in size) in plagioclase and thin rims of very fine-grained biotite on orthopyroxene crystals.
Metamorphosed syenodiorite.The recrystallised and deformed equivalent of the melanocratic syenodiorite (Fig. 2b,  d, g, h) is made up of sodic plagioclase (An 3−17 ), Kfeldspar, biotite, epidote-clinozoisite, quartz, garnet, titanite, https://doi.org/10.5194/ejm-35-773-2023A fine-grained biotite-rich foliation trail is seen in the lower part of the photos.Above is a fine-grained aggregate of recrystallised plagioclase with a euhedral garnet crystal and, at the top of the images, relict igneous mesoperthite.Sub-microscopic (10-30 µm) epidoteclinozoisite grains occur in all three domains (seen in Fig. 4).and calcite.All deformed and metamorphosed varieties of the melanocratic syenitoid are devoid of pyroxene.Sodic plagioclase predominantly forms fine-grained granoblastic aggregates, together with quartz and biotite in lens-shaped and flattened domains (Fig. 2b, d).Locally, coarse-grained mesoperthitic microstructures are still present in the cores of the feldspar domains, as remnants of the igneous feldspar.The ductile deformation fabric is defined by lensoid and flattened feldspar domains, thin and long dark mineral aggregates, and oriented biotite flakes (Fig. 2d, g-h).Fine-grained aggregates of metamorphic epidote-clinozoisite and carbonate are distributed within the felsic domains.Idiomorphic garnet crystals are present locally, containing inclusions of epidoteclinozoisite and quartz.The absence of strain shadows or deflection of the external foliation around garnet porphyroblasts, as well as the inclusion trails being continuous with the matrix foliation, indicate a late-tectonic or post-tectonic formation of the garnet.

Zirconium phases and their microtextural relationships
The accessory Zr phases in the near-pristine igneous syenodiorite are different from those in the metamorphosed varieties.Zirconolite and baddeleyite are present only in the near-pristine igneous variety, while the titanite-zircon pseudomorphs have formed solely in the recrystallised rocks.Below, we report the textural appearance of these different Zr phases.Metamorphic and igneous zircon crystals can be differentiated based on their shape and zoning (briefly described below); they have been dated and chemically characterised by SIMS (secondary ion mass spectrometry; these data are reported in Urueña and Möller, 2023).
Igneous baddeleyite.This is preserved only in the nearpristine igneous syenodiorite (sample 09j).Most baddeleyite crystals range in size from 20 to 100 µm.The baddeleyite blades are anhedral and commonly have up to 5 µm wide zircon rims (Fig. 3a).These zircon rims on baddeleyite are likely metamorphic, but a late igneous origin cannot be excluded based on the textural relationship only.
Metamorphic zircon.This is abundant in all samples of metamorphosed and recrystallised syenodiorite.Unlike igneous zircon, metamorphic zircon typically forms irregular crystals and/or aggregates with a "multi-nuclei" texture (Fig. 3c), commonly referred to as "polycrystalline" (e.g.Beckman and Möller, 2018, and references therein).The in-ner nucleus of each zircon grain is the most luminescent (Fig. 3c).
Zirconolite.This is rare and was found as very small grains only, commonly < 5 µm, in the near-pristine igneous syenodiorite only (sample 09j).The largest two zirconolite crystals are euhedral prismatic to subhedral, up to 12 µm long and 5 µm wide (Fig. 3d).All zirconolite grains are hosted in mesoperthitic feldspar, which commonly also contains small (inferred metamorphic) epidote-clinozoisite crystals.In places, zirconolite contains baddeleyite inclusions, indicating that baddeleyite crystallised before or simultaneously with zirconolite.The zirconolite crystals show low luminescence without internal microstructures (i.e.suggesting no zoning) and lower brightness than baddeleyite in backscattered electron images.The small zirconolite crystals appear pristine without cracks or alteration signs.
Titanite-zircon intergrowths.We interpret these to be pseudomorphs after zirconolite (see Sect. 6), and they are the main focus of this study.Figure 4 shows textural variations among these intergrowths and illustrates their appearances in secondary electron (SE), backscatter electron (BSE), and CL imagery.Titanite-zircon intergrowths are mainly hosted within feldspar domains; only a few were identified next to thin biotite foliae.Elongate or platy titanite crystals "speckled" by abundant tiny inclusions of zircon, < 2 µm large (Fig. 4a), are the most characteristic type of intergrowth.The titanite-zircon intergrowths may also be subrounded in shape (Fig. 4b).In places, the intergrowths have zircon preferentially located within the titanite core and along its rim (Fig. 4b).Titanite-zircon intergrowths were found mainly in metamorphosed syenodiorite but also within metamorphosed alkali-feldspar-dominated syenite (Fig. 4c).Typically, the intergrowths are spatially associated with fine-grained epidoteclinozoisite (Fig. 4a-c; presumably of metamorphic origin).Locally, bead-like, tiny (< 5 µm) grains of titanite and zircon are intergrown with albitic plagioclase and epidoteclinozoisite (Fig. 4d).Other titanite-zircon intergrowths appear to have undergone incipient recrystallisation/ripening, forming two distinct domains of slightly larger size (Fig. 4e); however, the zircon is still of multi-nuclei type.

Zirconolite: mineral and crystal chemistry
The intergrowths of titanite and zircon point to a precursor mineral rich in Ti, Ca, and Zr, suggesting that zirconolite is the most likely candidate.Indeed, tiny crystals with EDS spectra matching zirconolite were found in the nearpristine igneous syenodiorite.We attempted to prove that this phase is zirconolite using WDS and EBSD.Only two small crystals were large enough to analyse (> 5 µm).Three quantitative WDS analyses of these two crystals indicate a relatively limited compositional variation between the REE components: the one crystal large enough (12 µm) to allow the analysis of both core and rim indicates no compositional variation (Fig. 5, Table 3).Even though total sums are below https://doi.org/10.5194/ejm-35-773-2023 100 %, the measured weight percentage of each oxide ranges within the average content reported for natural zirconolites in syenitoid intrusions (cf.Williams and Gieré, 1996).Low totals appear to be common for natural zirconolite, and normalised analyses have been reported (cf.Wang et al., 2021;Williams and Gieré, 1988;Carlier and Lorand, 2008).The low totals in the zirconolite analysed in this study may presumably reflect the presence of elements that were not possible to detect during the elemental identification of the smallsized zirconolite.More importantly, the proportions of constituent elements fit the composition of zirconolite.
Like most natural zirconolite, the analysed phase in this study deviates from the ideal zirconolite composition, Ca 2 Zr 2 Ti  , Al 3+ , and Si 4+ filling vacancies in the structure; and a significant Fe content (most of which substitutes for Ti 4+ in site [V]M).The most abundant REEs are (in decreasing order of abundance) Y, Ce, and Nd, but one crystal has a higher amount of Ce and lower Y and Nd.An empiric formula for analysis Zrc 24_2 (Table 3; chosen for its highest oxide sum and consistent REE pattern; see Supplement, Fig. S1) was calculated based on 8 cations and 14 oxygen atoms: Following the crystal-chemistry models described for the zirconolite endmember CaREEZr 2 Ti 3 Me 3+ O 14 (Haifler et al., 2021), VIII Ca VII 2 Zr VI 2 Ti V 3 TiO 14 (Bayliss et al., 1989), and stefanweissite (Ca, REE) 2 (Zr, Hf, Mn) 2 (Ti, Nb) 2 (Nb, Ti) Fe 2+ O 14 (Chukanov et al., 2019), the components in crystallographic sites for this study were allocated as follows: Ca, REE, and minor bivalent cations for 2 apfu in the [VIII]M site.The [VII]M site is filled with Zr for 2 apfu, and [IV, V, VI]M sites are filled with Ti, Nb, and Fe 3+ (and additional trivalent cations) for 4 apfu.Thus, the simplified formula may be expressed as (Ca, REE) 2 Zr VI 2 Ti IV V 3 (Ti, Nb, Fe 3+ ) O 14 .
EBSD analysis was also performed on the largest (12 µm long) zirconolite crystal (Fig. 5), aiming to constrain the crystal structure.However, the analysis yielded no Kikuchi pattern for indexation of the crystal structure of zirconolite.During analysis, the Kikuchi patterns were non-existent or indistinct (see Supplement, Fig. S2).The EBSD results were thus inconclusive as a majority of the indexing resulted in zero solutions or yielded too high mean angular deviation values of > 1 • .

Discussion
In addition to zircon, the rocks within the Vaggeryd syenite complex host the accessory minerals baddeleyite, zircono-lite, and srilankite, likely due to their relative undersaturation in silica (SiO 2 < 55 % and absence of free quartz) and their Fe-Ti-rich and Zr-rich (> 180 ppm) bulk composition (see Supplement, Table S1).The mineral and textural relations between the near-pristine igneous and the gradually recrystallised syenodiorite varieties indicate that the intergrowths of titanite and zircon are the product of the breakdown of igneous zirconolite.As discussed below (Sect.6.3), the production of metamorphic zircon and titanite depends on the supply of SiO 2 from the consumption of the igneous silicates.

Formation of the igneous Zr minerals
The coexisting magmatic zirconolite, baddeleyite, and zircon, documented in the near-pristine igneous syenodiorite, were likely formed during late-stage crystallisation from Zrsaturated melts (cf.Söderlund et al., 2005;Schaltegger and Davies, 2017).The crystallisation of the magmatic accessory phases was controlled by the relative activities of Si and Zr.Baddeleyite can crystallise before zircon when the magma reaches Zr saturation before silica saturation (Schaltegger Figure 5. Backscatter electron image of one of two > 5 µm zirconolite crystals in near-pristine igneous syenodiorite (Fig. 2a, c, e-f).Closeup view indicating two EPMA spots (blue circles mark Zrc 24_2 and Zrc 24_3 in Table 3).and Davies, 2017).As a consequence, the uptake of Zr by baddeleyite may later be followed by zircon growth.Additionally, the presence of two textural types of igneous zircon in one sample hints at a magmatic system with the recharge of melts of different syenitic compositions.
The coexistence of Zr-bearing phases has been reported from alkaline and mafic intrusions such as gabbro and anorthosite (e.g.Davidson and van Breemen, 1988;Heaman and LeCheminant, 1993;Scoates and Chamberlain, 1995;Rasmussen and Fletcher, 2004;Beckman et al., 2017Beckman et al., , 2014)).In silica-undersaturated syenitic intrusions, baddeleyite is found as an accessory mineral in association with different combinations of zirconolite, titanite, calcite, rutile, and zircon in the presence of a pure CO 2 fluid phase (Lumpkin, 1999).Moreover, Lumpkin (1999) argues that Fe and Ti incorporation into the chemical system broadens the stability field of baddeleyite and may explain the coexistence of baddeleyite and zircon.This may also be the case for coexistence with Fe-rich zirconolite.As the minute inclusions of baddeleyite in zirconolite were observed only in two crystals, without evidence of dissolution and/or replacement, the relationship between both minerals is interpreted as cogenetic.Nevertheless, in contrast to more coarsely sized baddeleyite that is widely distributed in the rock, the zirconolite is typically hosted in the alkali-feldspar domains.This spatial relationship may indicate that zirconolite crystallised in locally supersaturated interstitial liquids near the alkali feldspar.
Alternative mechanisms that describe the formation of Zr oxides associated with zircon involve metasomatic reactions (e.g.Lewerentz et al., 2019;Tropper et al., 2007).During contact metamorphism in marbles, at metamorphic temperatures of 500-600 • C, the introduction of CO 2 -bearing fluids promoted the replacement of baddeleyite by metamor-phic zirconolite, and, as the Si activity increased, metamorphic zircon partly replaced the zirconolite (cf. Tropper et al., 2007).Lewerentz et al. (2019) demonstrate the formation of baddeleyite from the consumption of zircon in amphibolitefacies conditions, in silica-saturated rocks, caused by the ingress of a Ca-bearing fluid, which saturated the system in Ca.Spiridonov et al. (2019) suggested another mechanism involving the formation of zirconolite and baddeleyite via the dissolution of zircon xenocrysts in a shallow-level mafic magma enriched in Zr, REE, and Th or U.Although these three studies include evidence of the mineral relationships and petrogenetic sequences within these rocks, there are no textural indications of the same processes in the rocks investigated here.Except for a lower value for loss on ignition (LOI), the whole-rock composition of the near-pristine igneous syenodiorite matches that of its metamorphosed equivalents, indicating no external supply of Si, Ca, or other elements.Igneous zirconolite and baddeleyite are only preserved in the near-pristine igneous rock, and no relationships suggest that zirconolite replaces baddeleyite or vice versa.

Mineral characterisation of the igneous zirconolite
Natural zirconolite is rarely reported because of its dominant association with silica-undersaturated and alkaline rocks (e.g.Williams and Gieré, 1996;Gieré et al., 1998;Rasmussen and Fletcher, 2004).Our findings of pseudomorphs of zirconolite indicate that this mineral may be more common than generally assumed, although rarely identified due to its small size.
The composition of zirconolite in the Vaggeryd syenitoid can be described by the endmember CaREEZr 2 Ti 3 Me 3+ O 14 , with minor substitutions only.For example, the amounts https://doi.org/10.5194/ejm-35-773-2023 of Nb 5+ are too low for two of the other endmembers (Table 1).The mineral chemistry is consistent with Febearing zirconolite being poor in Al and Mg reported for silica-undersaturated syenitic rocks (Fowler and Williams, 1986;Williams and Gieré, 1996;Gieré et al., 1998).The charge-balanced stoichiometric formula calculation (based on 8 cations and 14 O atoms) indicates that most of the total Fe (6.1 wt %-7.9 wt %) corresponds to Fe 2 O 3 (5.6 wt %-6.1 wt %), which would hint at a high oxidation state during the late stage of crystallisation.However, this indication should be taken with caution because the calculation is sensitive to the weight percentage totals.
Although the zirconolite composition in the syenodiorite is closest to zirconolite-3O (cf.Bayliss et al., 1989;Spiridonov et al., 2019;Haifler et al., 2021) and the analogous polytype stefanweissite (Chukanov et al., 2019), the wholegrain EBSD map of the selected zirconolite crystal resulted in an extremely low contrast (black colour) image that indicates no patterns to be indexed for any mineral phase or other zirconolite polytype.The inability to determine a crystalline phase in this study may be a consequence of either immiscible intergrowths of zirconolite-like phases (e.g.zirkelite or zirconolite polytypes) during the igneous crystallisation or loss of crystallinity due to metamictisation.The most reasonable explanation for the lack of electron backscatter diffraction patterns is the amorphisation of the crystal.However, the analytical data in this study do not provide information to determine the cause of amorphisation.Further crystallographic analysis might ascertain the microstructural complexities of this mineral, but that is beyond the scope of this paper.

Metamorphic zircon-forming reactions
The formation of metamorphic zircon by the reaction of baddeleyite with a silica-bearing aqueous fluid has been documented in mafic rocks within a wide range of pressure and temperature (P -T ) conditions (Davidson and van Breemen, 1988;Lumpkin, 1999;Söderlund et al., 2005;Beckman et al., 2017;Beckman and Möller, 2018).Commonly, remnants of baddeleyite crystals are surrounded by polycrystalline zircon coronas.Metamorphic zircon-forming reactions from other Zr-bearing accessory minerals, specifically zirconolite, are less commonly discussed but have been reported by Pan (1997) andPohlner et al. (2020).Under subgreenschist conditions, the breakdown of zirconolite may involve changes by dissolution, alteration, and replacement (Gieré et al., 1998).
In the melanocratic syenodiorite and syenite studied here, tiny epidote-clinozoisite inclusions in the mesoperthite (10-30 µm, seen at high magnification, e.g.Fig. 4) are formed at the beginning of the metamorphic recrystallisation of a ternary feldspar that initially was richer in Ca.The main metamorphic mineral assemblages of the deformed (and undeformed) rocks, in particular the stability of epidoteclinozoisite and albite-rich plagioclase, bear witness to the reactions taking place under epidote-amphibolite conditions (supported by P -T estimates reported in Urueña and Möller, 2023).During metamorphic recrystallisation, the anhydrous mineral paragenesis in the near-pristine igneous syenodiorite was transformed into a mineral assemblage in which most ferromagnesian minerals are hydrous (biotite and epidoteclinozoisite). Thus, the mineral reactions, including the formation of zircon from Zr oxides, represent a specific stage during regional metamorphism when the influx of hydrous fluid associated with deformation triggered the breakdown of the (metastable) igneous minerals to form the metamorphic assemblage.
A general reaction can be balanced, which involves the main minerals in the near-pristine syenodiorite and is dependent on the introduction of hydrous fluid (Urueña and Möller, 2023) The breakdown of the igneous minerals (i.e.dominantly feldspars and orthopyroxene, Reaction 1) released Si (quartz) that was important for the formation of zircon from baddeleyite and zirconolite, and the Ca from calcic plagioclase component was used for the formation of metamorphic clinozoisite, titanite, and garnet.The production of biotite and clinozoisite from the dry pyroxene-bearing igneous assemblage necessarily involved the introduction of H 2 O.This is also corroborated by the higher LOI values for the metamorphosed rocks than for the near-pristine igneous (Supplement, Table S1).
The breakdown of zirconolite into metamorphic zircon and titanite is a mineral reaction that is significantly subordinate in volume to that of the major minerals (Reaction 1) but which was also a consequence of the metamorphism in epidote-amphibolite-facies conditions.It can be described by the following simplified reaction: It is emphasised that igneous zircon crystals in the most deformed varieties lack textures that could suggest metamorphic dissolution and recrystallisation of the igneous zircon crystals, e.g.embayments or overgrowths.The polycrystalline zircon textures and the thin zircon rims around baddeleyite (Figs. 3,4) suggest that metamorphic zircon formed from the consumption of baddeleyite and zirconolite.The result is that metamorphic zircon and igneous zircon are both present in the metamorphosed syenodiorite.
The formation of metamorphic zircon and titanite after zirconolite was first described by Pan (1997) from volcanosedimentary rocks in Manitouwadge, Ontario (Canada), which underwent metamorphism and anatexis in high-grade conditions.Textural relationships within pseudomorphs attest to the breakdown reaction, Reaction (R4), of zirconolite to zircon + titanite + rutile, resulting in high contents of Nb and REE in the product minerals: This transformation of zirconolite to metamorphic zircon + titanite + rutile was observed in rocks that underwent metamorphism in high-pressure upper-amphibolitefacies to granulite-facies conditions (Pan, 1997).In the syenodioritic rocks studied herein, rutile is absent and epidote present, which suggests that Reaction (R2) is representative of epidote-amphibolite-facies conditions.

Conclusions
In the syenitoid rocks described in this study, the formation of titanite + zircon after zirconolite took place during metamorphism under epidote-amphibolite-facies conditions and was associated with the breakdown of igneous ternary feldspar and orthopyroxene.The field relationships and the formation of hydrous ferromagnesian phases (and increase in bulk LOI) suggest that this metamorphic transformation was triggered by hydrous fluid infiltration associated with ductile deformation.Titanite + zircon intergrowths, typically titanite speckled with micrometre-sized zircon and, less commonly, polycrystalline aggregates, are characteristic pseudomorphs of zirconolite.The precise reaction for the formation of metamorphic titanite + zircon intergrowths will depend on the wholerock composition and the pressure-temperature conditions during metamorphism.However, in all cases, analogous to the formation of metamorphic zircon after baddeleyite, it depends on the release of silica from the breakdown of silicate minerals in the rock.
The formation of metamorphic zircon in silicaundersaturated rocks takes place during the breakdown of precursor igneous Zr oxides such as baddeleyite and/or zirconolite.Metamorphic titanite-zircon intergrowths represent pseudomorphs after igneous zirconolite, which are analogous with polycrystalline zircon pseudomorphs after igneous baddeleyite.These textural relationships allow the linkage of the metamorphic zircon formation and a regional metamorphic event.
Financial support.This research received financial support from the Royal Physiographic Society of Lund on behalf of Endowments for the Natural Sciences, Medicine and Technology -Geoscience 2019 through a grant to the project "Tracking mountain-building processes by the coupled records of garnet and zircon -towards an increased understanding of deep processes in collisional orogens".This study is part of a PhD project funded by the Geological Survey of Sweden via grant no.36-1957/2017 to Charlotte Möller.
Review statement.This paper was edited by Samuel Angiboust and reviewed by Jan De Hoog and Daniel Harlov.

Figure 1 .
Figure 1.Tectonic and geological setting.(a) Simplified tectonic map of the Sveconorwegian Province with the internal subdivision of the Eastern Segment(Möller et al., 2015).The Vaggeryd syenite complex is located within the Frontal wedge and shown in reddish-brown.(b) Overview bedrock map of the southern part of the Vaggeryd syenite complex and surroundings (based on a 1 : 250 000 scale database), including previously unpublished structural data from Rikard Ask.The studied locality north of Näsbyholm is marked with a green spot (source: map viewer Geological Survey of Sweden, SGU; coordinate system: SWEREF 99; meta-TIB: metamorphosed Transscandinavian Igneous Belt rocks).(c) Sketch illustrating progressive changes in deformation and metamorphic recrystallisation in the studied syenodiorite, with the near-pristine igneous rock preserved in the core of a lens with weakly to moderately deformed and metamorphosed margins.

Table 2 .
Analytical setup and calibration standards used for the electron microprobe analyses of zirconolite.

Table 3 .
Chemical compositions of zirconolite and stoichiometry for the empirical formula.The analysed zirconolite crystals are shown in Fig.3d(analysis Zrc 64_1) and Fig.5(analyses Zrc 24_2 and Zrc 24_3).