Compositional variation and zoning of epidote supergroup minerals in the Campi Flegrei geothermal ﬁeld, Naples, Italy

.


Introduction
The Campi Flegrei area, located west of Naples, Italy, contains numerous hot springs, fumaroles, and other geothermal manifestations, which were well known before the Roman Empire (Giacomelli and Scandone, 2012).Geothermal exploration, in an attempt to locate a productive field similar to Larderello, was intermittently conducted from 1939 to 1954 but was unsuccessful.Carlino et al. (2012) present the history of geothermal exploration in the Neapolitan volcanic region.Extensive drilling and geophysical exploration programs were resumed in 1978 as a joint venture of the Italian national utilities AGIP and ENEL and the Italian Geodynamic Project (Rosi and Sbrana, 1987).Wells were drilled to a depth of ∼ 3 km with concomitant coring, geophysics, and fluid sampling.Although significant geothermal electric power production was never achieved in the Campi Flegrei region, this exploration stimulated widespread research (e.g., Guglielminetti, 1986;Guidi and Antonelli, 1989;De Vivo et al., 1989;Altaner et al., 1991;Zamora et al., 1994;Caprarelli et al., 1997).
As part of the research which evolved from this extensive drilling program, a detailed fluid inclusion study was carried out to better understand the evolution of the fluid composition and the thermal history of the geothermal field (De Vivo et al., 1989).We have used the results of this fluid inclusion study to better understand, in some wells, the variation in epidote group composition.As part of this study, core and chip samples from three Mofete (MF1, MF2, MF5) and two San Vito (SV1, SV3) wells (Fig. 1) were carefully examined for the occurrence of fluid inclusions and to determine their host mineralogy.Members of the epidote group, epidote and clinozoisite, were identified in many samples in the five studied wells and occurred throughout much of the core except at shallow depth.Preliminary studies (Belkin et al., 1990;Belkin and De Vivo, 2004) of their mineral chemistry revealed that their composition varied from epidote to clinozoisite and had different zoning textures.Two members of the allanite group, allanite-(Ce) and ferriallanite-(Ce), were also identified as authigenic phases in the core samples.
The specific aims of this report are (1) to document the composition and petrography of the epidote group in the Mofete and San Vito well samples, (2) to illustrate the various types of zoning observed in the epidote group, and (3) to document the composition and petrography of the allanite group.We know of no other compositional data for epidote supergroup minerals identified in the Mofete or San Vito geothermal fields.

Geologic setting of the Campi Flegrei region
The Campi Flegrei caldera complex (Fig. 1) is part of the active volcanic region developed on the subsiding western margin of the Apennine chain.This region also contains the Somma-Vesuvius, Procida, and Ischia volcanoes, not shown in Fig. 1.Although the region has been intermittently active during the last 600 kyr (Belkin et al., 2016;De Vivo et al., 2020), a major event occurred in 39 ka (De Vivo et al., 2001) with the very large eruption (∼ 310 km 3 ) of the Campanian Ignimbrite (Rolandi et al., 2003).Subsequent to this ignimbrite eruption, the Campi Flegrei region has had much volcanic activity ending with the Monte Nuovo eruption of 1538 CE (Rosi and Sbrana, 1987).For more detail regarding the volcanism of the Campi Flegrei area, the reader is referred to De Vivo et al. (2020) and Orsi et al. (2022) and chapters therein.Bradyseism and fumarolic activity continue, attesting to a still-active hydrothermal system (e.g., De Vivo and Lima, 2006;Bodnar et al., 2007;Lima et al., 2009;Bagnato et al., 2020;Lima et al., 2021;Chiodini et al., 2022).The chemistry of all the igneous products is Siundersaturated, rich in alkalis, and low in MgO.Volumetrically, trachyte and alkali-trachyte are predominant with much lower amounts of trachybasalt, latite, and peralkalinephonolitic-trachyte (Rosi and Sbrana, 1987).The chemistry of the rock penetrated by the Mofete and San Vito wells is generally similar from top to bottom although some wells (Fig. 2) end in hydrothermally metamorphosed marine marls, silts, and sandstones (Rosi and Sbrana, 1987).Piochi et al. (2014) give more stratigraphic detail on the MF5 and SV1 wells, and Mormone et al. (2011) provide more detail for the MF5 well.

Studied samples
All the samples came from the drilling program of AGIP and ENEL and include core and drill cuttings; however, the available samples were not a comprehensive selection of the core, especially for the shallower depths, but emphasized the deeper, hotter parts of the geothermal system.Thus, we do not have data for the shallower samples where epidote was identified by X-ray diffraction (XRD) and optical petrography (De Vivo et al., 1989).The Campi Flegrei geothermal system is liquid dominated with evidence of boiling and fluid reservoirs of different salinity (De Vivo et al., 1989;Caprarelli et al., 1997).The details of the studied samples are described in De Vivo et al. (1989).The studied samples containing epidote supergroups are listed in Table 1 and come from the Mofete 1 (MF1), Mofete 2 (MF2), Mofete 5 (MF5), San Vito 1 (SV1), and San Vito 3 (SV3) exploratory geothermal wells (Fig. 1).On the basis of detailed petrography and XRD determinations, Chelini and Sbrana (1987) defined four descending hydrothermal zones with different characteristic mineralogy: argillic, chlorite-illite, calciumaluminum silicate, and the thermometamorphic zones.The calcium-aluminum silicate zone is characterized by abundant epidote, quartz, feldspars, various sulfides, and layered silicates.The top of this zone was stated by Chelini and Sbrana (1987) to be about 250 • C in the Mofete area and was less well defined in the San Vito area at between 220 and 270 • C (Fig. 2).The bottom of this zone was judged by Chelini and Sbrana (1987) to be at 325 • C in the Mofete area,  (1989).The top of thermometamorphism, total depth of well, and isotherms ( • C) are indicated.Lithologic symbols: 1 -incoherent pyroclastic rocks; 2 -latitic soriae and pumices; 3 -yellow tuffs (post-caldera submarine activity); 4 -chaotic tuffites; 5 -chaotic tuffites (subaerial environment); 6Atrachytic lavas; 6B -latitic lavas and lava domes (pre-caldera period); 7 -interbedded tuffites, tuffs and lavas (submarine environment); and 8 -interbedded shales, marls, siltstones, and sandstones (marine).
whereas in San Vito 1 it corresponds to the current 360 • C isotherm, and in San Vito 3, it corresponds to the current 270 • C isotherm.This lower temperature in San Vito 3 reflects the complexity of recent faulting shifting older hydrothermal zones up and down (Fig. 2).The thermometamorphic zone is characterized by nearly complete recrystallization and/or replacement of the host rock; Chelini and Sbrana (1987) defined two subzones: amphibole-biotite and diopside.All the studied samples come from the calciumaluminum silicate and the thermometamorphic zones as defined by Chelini and Sbrana (1987).
Recently, portions of the Mofete and San Vito cores have been examined in the following studies in order (1) to study the effects of water-rock reactions to aide in geochemical modeling (Piochi et al., 2021), (2) to develop a revised model for the subsurface structure of the Campi Flegrei (Piochi et al., 2014), and (3) to examine the secondary hydrothermal minerals in order to better understand the rock physics and state of the volcanic system (Mormone et al., 2011).Only in the studies by Mormone et al. (2011) and Piochi et al. (2021) is the epidote supergroup described as a component of the hydrothermal mineral assemblage, but no compositional data are given.

Nomenclature
In this report, we use the International Mineralogical Association (IMA)-recommended terminology as outlined by Armbruster et al. (2006) and Mills et al. (2009).The epidote supergroup is subdivided into groups, such as the epidote, allanite, and dollaseite groups.In our study, we identified only two members of the epidote group (clinozoisite, epidote) and two members of the allanite group (allanite-(Ce), ferriallanite-(Ce)).In the recommended nomenclature, the occupancy and dominance of key cation sites determine the group, root name, and individual species (Armbruster et al., 2006).
In the allanite group, A2 is occupied by  (De Vivo et al., 1989).
rare earth element (REE + Y).In the IMA-recommended system (Armbruster et al., 2006), the allanite group is derived from the epidote group by homovalent substitutions and one coupled heterovalent substitution of the type (REE) 3+ A2 + M 2+ M3 → Ca 2+ A2 + M 3+ M3 .Thus, the valences on the key sites, A2 and M3, are the reverse of the epidote group.Gieré and Sorensen (2004), in a comprehensive review of allanite and REE-bearing epidote mineralogy, point out that the distinction of allanite from REE-rich epidote differs by author.In the IMA procedure, REE-bearing epidote and clinozoisite are distinguished from allanite by the dominant cation(s) (> 0.50) in the A2 and M3 sites.The dominance of Fe 3+ in the M1 site distinguishes ferriallanite from allanite.
To normalize the electron microprobe (electron probe microanalyzer, EPMA) data, we used the program WinEpclas (Yavuz and Yildirim, 2018), which follows the IMArecommended naming procedure based on site occupancy and cation dominance.Normalization of EPMA data was on the basis of (T + M + A) = 8 cations.

Analytical methods
Polished thin sections were prepared either from pieces of core or from individual hand-picked crystals.After routine optical petrography, the samples were examined with a JEOL JSM-840 scanning electron microscope (SEM) equipped with a Princeton Gamma-Tech Si(Li) energy-dispersive Xray analyzer or a Hitachi SU5000 Schottky thermally as-sisted field-emission SEM equipped with an Oxford Ultima 100 mm 2 energy-dispersive X-ray spectroscopy (EDS) silicon drift detector to select areas for analysis.Quantitative electron microprobe analyses of major and minor elements were obtained with a JEOL JXA-8900R five-spectrometer, fully automated EPMA using wavelength-dispersive X-ray spectrometry (WDS).Analyses were made at 15 keV accelerating voltage and 20 or 30 nA probe current measured with a Faraday cup, and counting times on both the peak and the background varied from 20 to 120 s.

Results -epidote group
Epidote and clinozoisite, members of the monoclinic epidote group, occur in a wide diversity of parageneses within a great range of pressures, temperatures, and host rock compositions (e.g., Zen and Hammarstrom, 1984;Deer et al., 1986;Bird and Spieler, 2004;Chen et al., 2019;Ahmed et al., 2020).The general formula for the epidote group miner- Other cations can substitute in minor amounts; A sites can contain large, high-coordination cations such as Sr, Pb, REE, Y, Th, Na, and K, and octahedrally coordinated M sites can contain Mn 3+ , Mn 2+ , Fe 2+ , Ti, Cr 3+ , and V cations.REEs have been measured in the A2 site of epidote group minerals from a wide variety of paragenetic environments such as magmatic (Dawes, 1990), ultrahigh pressure (Rolfo et al., 2000), various metamorphic facies (Smulikowski and Kozłowski, 1994; see also Gieré and Sorensen, 2004), and geothermal (Aggarwal, 1999;Bird and Spieler, 2004).

Petrography
In the studied samples, the epidote group occurs in a wide variety of habits, shapes, and sizes.Very common throughout the cores were vugs (Fig. 3a) and veins containing hundreds to just a few crystals (Fig. 3b).Also, we observe crystal vugs completely or partially filled by later deposition of K-feldspar (Fig. 3c), quartz, and calcite.Limited optical petrographic study suggests that epidote crystals appear to be elongated in the b-crystallographic direction but have various aspect ratios.Some crystals were relatively stubby (Fig. 3c), and some had length-to-width ratios of as much as 10 : 1 (Fig. 3a).No systematic distribution of sizes as reported by Patrier et al. (1990) was found, but it appeared that the mode of occurrence influenced the size.Larger crystals are observed in larger veins or vugs, and smaller crystals are noted in the thinner veins and smaller pores.Another very common epidote group occurrence is in masses of crystals of up to 1.5 cm usually cored by another earlier phase (Fig. 4a, b).Sulfides, allanite-(Ce), ferriallanite-(Ce), cerite-(Ce), aluminocerite-(Ce), titanite, britholite-(Ce), and a TiO 2 phase were observed in the cores of epidote group masses.Epidote was observed to fill probable vesicles with a welldefined outer edge (Fig. 4c) and partially filled cavities with a more irregular boundary (Fig. 4d).Veins, commonly transecting the thin section, were usually accompanied by pyrite or pyrrhotite (Fig. 4e), and the individual crystals comprising the veins were invariably zoned (Fig. 4f).

Epidote group composition
The clinozoisite (czo)-epidote (ep) members of the epidote group are related by the exchange of Fe 3+ ↔ Al in the M3 site; we used X Fe [Fe 3+ / (Fe 3+ + Al)] apfu as a measure of this exchange.

San Vito samples
In the San Vito core samples, X Fe ranged from 0.10 to 0.32, where SV1 varied for ep from 0.25-0.32 and for czo from 0.11-0.16and SV3 ranged from czo 0.10-0.16 to ep 0.17-0.28(X Fe ).The San Vito epidote group data show a good linear co-variation in the Fe 3+ -Al substitution (Fig. 5b).Aluminum varies from 2.03-2.68apfu, and Fe 3+ varies from 0.939-0.295apfu.Minor elements analyzed in San Vito clinozoisite and epidote were Ti, Mn, Mg, Sr, K, Na, Pb, Cr, Y, and Th.Manganese, Ti, Mg, and Sr were the most abundant, with MnO ranging (in wt %; dl denotes detection limit) from dl-0.85,TiO 2 ranging from dl-0.50, MgO ranging from dl-0.36, and SrO ranging from dl-0.48.Potassium, Na, and Y were very rarely detected with most analyses just above detection; Pb, Cr, Th, and Y were not detected.The minor elements in the San Vito epidote group were, in general, less abundant than in the Mofete epidote group.No major gap in X Fe (Fig. 5a, b) that would suggest immiscibility (cf.Holdaway, 1972;Brunsmann et al., 2002) was observed in either the Mofete or the San Vito samples.Table 3 gives representative compositions of the San Vito epidote group.(De Vivo et al., 1989).Fe 3+ / Fe 2+ is calculated on the basis of stoichiometry; Z, average atomic number.

Table 3.
Representative compositions of epidote and clinozoisite from the San Vito samples.

Compositional zoning in the epidote group
Various styles of zoning in epidote group minerals were observed in all the studied Mofete and San Vito well samples.
Zoning involving a smooth gradation either from a Fe-rich core to an Al-rich rim or vice versa was not observed.All the compositional zones were separated by relatively sharp boundaries, at least within the resolution limits of the SEM, although in some complex zoning areas, these were difficult to clearly resolve.We observed oscillatory, complex, and sector zoning.Oscillatory zoning was most clearly observed when the crystal happened to lie lengthwise in the plane of the thin section.Figure 7a shows a portion of an 850 µm crystal oriented so that the oscillatory zoning is clearly displayed.The core (right arrow) is a partially resorbed epidote (spot A).In the same vug, Fig. 7b shows another partially resorbed core in a crystal oriented perpendicular to Fig. 7b.The core of Fig. 7b is iron-rich (spot C), but the crystal growth following the resorption was relatively iron-poor (spot B).Two traverses are conducted on the crystal in Fig. 7a: (1) WDS line scan recording the Fe Kα and Al Kα counts per second and (2) WDS analytical analysis every ∼ 25 µm.Figure 7c shows the line scan across the variable-width oscillatory zones.Two zones, spots D (ep) and E (czo), are separated by a very sharp boundary, although because the microprobe beam has finite volume, the line scan transitions are somewhat sloped.Inset   3 gives their complete analysis.
Complex zoning was the most common type observed in both Mofete and San Vito samples and was usually mixed with other zoning types.It was recognized in crystals growing in vugs (Fig. 8a), in masses of crystals, and in veins (Fig. 4f). Figure 8a shows a portion of a vug (MF2-1824) with some of the crystals with britholite-(Ce) centers.The enlarged portion (Fig. 8b) of a euhedral crystal shows a chaotic assemblage of zones in the center that are terminated by os-cillatory zoning.Spots analyzed in Fig. 8b show that the center is iron-poor relative to the peripheral zones.Figure 8e gives the X Fe data for selected spots, and their EPMA data are given in Table 2.
Sector zoning refers to a compositional difference between coeval growth sectors in a crystal.This type was less common than either complex or oscillatory zoning in the samples and was typically seen in crystals in vugs (Fig. 8c   shown in Fig. 8d also illustrates the mixtures of zoning types: oscillatory and complex zoning with sector zoning in some of the late-formed crystals.The X Fe values are shown in Fig. 8e and the complete EPMA data in Table 2. Sector zoning refers to a compositional difference between coeval growth sectors in a crystal.This type was less common than either complex or oscillatory zoning in the samples and was typically seen in crystals in vugs (Fig. 8c).Where measured by the EPMA, the different sectors had significantly different X Fe values and BSE responses.The zoning shown in Fig. 8d also illustrates the mixtures of zoning types: oscillatory and complex zoning with sector zoning in some of the late-formed crystals.The X Fe values are shown in Fig. 8e and the complete EPMA data in Table 2.

Results -REE-bearing epidote supergroup
We have separated, for description, members of the allanite group from those of the epidote group that had detectable REE + Y.The critical distinction that separates the allanite group from other members of the epidote supergroup containing REE + Y is the dominance (> 50 % apfu) of REE + Y in the A2 site (Armbruster et al., 2006).

Allanite group
We identified two members of the allanite group, allanite- 1 gives the distribution of the allanite group in the studied samples.The allanite group was much more abundant in the Mofete + Al)] apfu for the EPMA analytical spots A-F for epidote (ep) and clinozoisite (czo) shown in Fig. 7a, b, and c; analyses are given in Table 3. wells than in the San Vito wells, especially in samples MF1-1597 and MF2-1824.

Petrography
We observed two general textures of the allanite group; (1) radiating sprays of crystals, usually at the centers of epidote group masses, and (2) groups of crystals, many euhedral, in the center of epidote group masses.These textures are authigenic, precipitated from geothermal fluids, as opposed to uncommon igneous allanite fragments we identified in the igneous products.The starburst pattern (Fig. 9a, b) was also seen without any accompanying epidote.An uncommon occurrence, only observed in sample MF1-1597, was allanite groups with cores of aluminocerite-(Ce) (Al 2 O 3 = 2.0 wt % to 2.3 wt %) (Fig. 9c, f).Oscillatory zoned crystals (Fig. 9d) and other textures (Fig. 9e) indicate precipitation of allanite in pulses.Anenburg et al. (2015, Fig. 6a) illustrate similar oscillatory zoned allanite.Ferriallanite was only observed in one sample (MF1-1597) and forms euhedral crystals in the centers of larger epidote masses (Fig. 9f).Minor zoning of https://doi.org/10.5194/ejm-35-25-2023+ Al) apfu] for the EPMA analytical spots G-K for epidote (ep) and clinozoisite (czo) shown in Fig. 8b and d; analyses are given in Table 3. light rare earth elements (LREEs, La-Gd) was observed in both allanite and ferriallanite.

Composition
The allanite and ferriallanite are enriched in LREEs and have variable compositions, mostly due to minor zoning.We did not observe any compositional relationship with temperature and/or depth.Both minerals are Ce dominant and hence are allanite-(Ce) and ferriallanite-(Ce).Minor elements analyzed in the allanite group were Ti, Mn, Mg, Sr, K, Na, Pb, Cr, Y, and Th.Manganese, Sr, Mg, and Ti are the most abundant with MnO ranging (in wt %; dl denotes detection limit) from dl-1.70,SrO ranging from dl-1.41,MgO ranging from dl-1.19, and TiO 2 ranging from dl-0.82.Potassium, Na, Y, and Th were infrequently detected with all analyses < 0.50 wt %; Pb and Cr were not detected.The LREE abundance (in wt %) is variable: La 2 O 3 varies from 3.98-8.30;Ce 2 O 3 varies from 7.11-16.17;Pr 2 O 3 varies from dl-1.41;Nd 2 O 3 varies from 1.05-4.40.When measured, Sm 2 O 3 ranges from dl-0.61,Eu 2 O 3 ranges from dl-0.35, and Gd 2 O 3 ranges from dl-0.37.Representative compositions of the allanite group are given in Table 4.

REE-bearing epidote group
During routine BSE observation, epidote group crystals, especially those bordering masses of allanite and ferriallanite (Fig. 9f, g), had bright rims or zones.Microprobe analysis yielded LREEs but in concentrations such that either Al or  4.
Fe 3+ was dominant (> 50% apfu) in the A2 site (Fig. 9g points N and O).In some cases, it appears that the transition from the REE-bearing epidote group to the allanite group is gradual (Fig. 9d), and in other cases, the transition seems sharp within the resolution limits of the SEM (Fig. 9e).Bonazzi and Menchetti (1995) suggest that the increasing content of Fe 2+ in epidote distorts the structure to accommodate more REE.We observe this relationship (Table 4), but our data set is too limited to be definitive.Table 4 gives representative compositions of REE-bearing clinozoisite and epidote; all LREE compositions are Ce dominant.

Geothermal epidote
Although epidote was described as a mineral species more than 180 years ago (Franz and Liebscher, 2004), the world energy "crises" of the 1970s and 1980s produced much exploration for high-enthalpy geothermal fluids for electric power production, and mineralogical studies associated with this exploration generated numerous literature descriptions of epidote from geothermal fields (e.g., Cavarretta et al., 1980Cavarretta et al., , 1982;;Shikazono, 1984;Bird et al., 1988).Bird and Spieler (2004)     dote in geothermal systems and emphasize that the epidote group is common and abundant in geothermal systems at temperatures usually > 200 • C, exhibits a wide range of octahedral substitution (Al ↔ Fe 3+ ), and displays frequent zoning.For example, Absar (1991) uses the presence of epidote to model the epidote-formation regime at T ≥ 250 • C in the Ohaaki-Broadlands, New Zealand, geothermal field.Bird and Spieler (2004) and Caruso et al. (1988) emphasize the sensitive dependence of epidote composition (coupled Al ↔ Fe 3+ substitution) on temporal and spatial variations in fluid composition, redox conditions, pH, CO 2 concentrations, aqueous speciation of Fe and Al complexes, etc. Bird and Spieler (2004, p. 284) conclude their review by warning the epidote researcher that "interpretation of the chemical and physical evolution of magma-hydrothermal systems based on chemically zoned epidotes must be considered with caution".

Epidote composition and well temperature
Figure 6a and b show the co-variation in X Fe with temperature for the Mofete and San Vito wells.The San Vito system has been disturbed by faulting, and this is evident by the distribution of epidote compositions, structural relations, and drill-hole stratigraphy (Chelini and Sbrana, 1987;De Vivo et al., 1989).In contrast, the Mofete system has been described as a simple, generally prograde, active, and structurally undisturbed area.Although there is scatter, we found decreasing X Fe with increasing well temperature and more clinozoisite below 320 • C. De Vivo et al. (1989) assessed the fossil temperature variation in the boreholes by fluid inclusion microthermometry for comparison with the measured well temperature.The fluid inclusion data indicate relatively stable thermal conditions for MF1, MF2, and SV1, but the data suggest that the thermal regime in well MF5 has cooled subsequent to fluid inclusion trapping and host mineral crystallization.Overlapping mineral assemblage zones in the SV3 well indicate that minerals that hosted the measured fluid inclusions are relatively recent and that a significantly hotter regime existed in the past (De Vivo et al., 1989).Guglielminetti (1986) and De Vivo et al. (1989) describe multiple fluid reservoirs in the Mofete field characterized by differences in the salinity and total dissolved solids (TDSs), but we did not observe any effect of these differences on the chemistry of distribution of the epidote supergroup.Liou (1973) suggests that Al increases in epidote with increasing temperature on the basis of experimental studies.Poli and Schmidt (2004) remark that generally epidote becomes more aluminous as oxygen fugacity decreases, but this relationship was equivocal.Potel et al. (2002) also describe epidote composition in New Caledonia metamorphics that varies as a function of oxygen fugacity and that rock composition is the main factor that controls both mineral assemblages and their composition.Liou et al. (1985), Caruso et al. (1988), and Arnason and Bird (1992) report decreas-ing X Fe with increasing temperature and/or depth, whereas Shearer et al. (1988) report no systematic variation with depth.In fact, many descriptions of geothermal epidote report no relationship between X Fe and temperature and/or depth and describe complex epidote histories, such as multiple generations, and disequilibrium mineral assemblages (Bird and Spieler, 2004, and references therein).In particular, Bird and Spieler (2004) illustrate the complexity of the influence of these variables by pointing out that the Al content of epidote will both increase and decrease with increasing temperature as a function of CO 2 fugacity (e.g., Arnason et al., 1993).
Iron in the Mofete and San Vito geothermal fluids was a mixture of Fe 2+ and Fe 3+ as there was common simultaneous precipitation of various divalent iron sulfides with epidote group minerals.We also identified iron oxides in all the samples.Hematite, intergrown with epidote (Fig. 4d), was observed in some samples, and this suggests mineralization from more oxidizing fluids as also suggested by McKibben et al. (1988) for the Salton Sea geothermal system.The sulfur fugacity also varied although from the restricted occurrence of anhydrite, most sulfur was relatively reduced.However, the fact that a general decrease in X Fe with increasing temperature is observed at all, especially in the MF1 well, suggests that the most important factors that determined X Fe distribution have remained relatively constant during the formation of the measured epidote group.

Zoning
Zoning in rock-forming minerals was described in the 19th century from some of the first petrographic studies of rocks using polarized light microscopy (e.g., Sorby, 1858).Different types of zoning have been subsequently recognized, and various theories have been put forward to explain the various types of zoning.Some reported epidote zoning seems related to the host rock composition where, for example, more iron-rich host rocks produce more iron-rich epidotes (Shikazono, 1984), whereas in other epidote studies (Santaguida, 1999) where host rock composition widely varies, no relationship was observed.De Vivo et al. (1989) and Chelini and Sbrana (1987) indicate no significant differences in host rock composition with regard to iron content in the Campi Flegrei lithologies from which the studied samples were selected.
Relatively simple zoning, where epidote crystals have either Al-rich or Fe-rich rims, has been explained by Arnason and Bird (1992) as the reaction path for plagioclase and calcite dissolution, respectively.Arnason and Bird (1992) also emphasize the sensitivity of epidote composition to such factors as pH and Fe 3+ and Al speciation.The complex and widely varied epidote zoning observed in the Campi Flegrei samples suggests complicated and relatively rapid fluctuation in intensive parameters to produce a great variability in zoning.Furthermore, epidote crystal chemical sites, such as M1, https://doi.org/10.5194/ejm-35-25-2023M2, and M3, seem to be especially sensitive to these condition changes.
Oscillatory zoning can arise from various oscillations in the growth medium (Loomis, 1983) or by autonomous, time-independent oscillations that are self-organized such as those in reaction-diffusion systems, like the Belousov-Zhabotinsky reaction type (Ortoleva, 1994).Figure 7a shows a typical epidote crystal with well-developed oscillatory zoning.A detailed, high-resolution SEM-BSE examination of the oscillatory zones in the Campi Flegrei epidotes did not show any sawtooth periodic decay, nor did it show any other pattern typical of those described in self-organized zoning.The oscillatory zoning was separated by sharp interfaces, and in any single band, the composition was usually uniform within the limits of the available instruments (Fig. 7a,  c).The alternation of Fe-rich and Fe-poor bands suggests oscillations in the fluid media such as rapid changes in pH, f CO 2 , f O 2 , or composition.The sharpness of the individual zones, at least on the scale of the SEM resolution, suggests that the fluid diffusion rates for the local conditions surrounding the growing crystals were rapid, but the solid diffusion rates were essentially zero with respect to this environment.No resorbed edges among the oscillatory bands were seen in the studied samples, suggesting a relatively rapid prograde environment.
Oscillatory zoning in epidote is common (e.g., Janeczek and Sachanbinski, 1992;Choo, 2002;Inoue and Utada, 2017;Bird and Spieler, 2004, and references therein).Choo (2002) describes oscillatory zoning in the Bobae sericite deposits, and microprobe experiments across oscillatory bands yield amplitudes which range from 2 to 13 X Fe .This is the general range measured in the oscillatory zoned epidote group in the present study (4-10 X Fe ).
Sector zoning is zoning localized at particular crystallographic orientations with differences in composition between sectors.Dowty (1976) discusses this zoning type with respect to epidote formation and indicates that the M1 and M2 sites prefer Al very strongly over Fe 3+ , whereas the M3 site prefers Fe 3+ , although Fehr and Heuss-Aßbichler (1999) suggest that the M2 is occupied by Al exclusively.Sector zoning in epidote can be quite complex: Yoshizawa (1984) describes an example in a low-grade schist where the epidote is zoned in 10 different crystallographic sectors.Thus, epidote growth in a fluid with both Fe and Al species available may produce zones controlled by crystallographic compositional differences.However, other factors need to be favorable to develop sector zoning; probably a relatively stagnant fluid is needed for the site-selective compositional differences to be operative.For example, Fig. 8d shows sector zoning in some sample areas although interfering crystal growth areas have complicated the pattern.
Complex zoning, the most common type observed in the Campi Flegrei epidotes, is a mixture of oscillatory banding, sector zoning, and patchy growth.Figure 8a and b show typical complexly zoned epidotes in vugs, and Fig. 4f shows this zoning type in a vein.Zoning of any type was more common in the deeper, hotter parts of the Campi Flegrei geothermal field, but this also may, in part, be influenced by the fact that epidote was also more abundant.
Thus, there are two sources of compositional variability identified in this study: one related to oscillatory zoning in response to intensive parameter and or fluid composition change and the other related to crystallographic features such as sector zoning.Unfortunately, we know of no analyses of REEs in the geothermal fluids sampled during the drilling of the Mofete and San Vito wells (De Vivo et al., 1989) or any other nearby geothermal wells (Somma et al., 2021) or thermal springs and fumaroles (Valentino et al., 1999).However, REEs have been analyzed in high-temperature fluids from the Larderello-Travale and Monte Amiata geothermal fields, Tuscany, Italy (Möller, 2002), and from the Reykjanes geothermal system, Iceland (Fowler and Zierenberg, 2015).In these two studies, the REEs were LREE enriched with some variation in the degree of enrichment.Williams-Jones (2015) indicates that LREEs form more stable complexes in chloride and fluoride solutions than in heavy rare earth elements (HREEs, Tb-Lu).This is compatible with the San Vito and Mofete geothermal saline fluids (De Vivo et al., 1989) and the common presence of fluorite and fluorapatite.

Conclusions
We conclude the following.(1) Authigenic epidote group minerals, epidote and clinozoisite, are common to abundant in the calcium-aluminum silicate and thermometamorphic zones of the Mofete and San Vito geothermal fields, Campi Flegrei, Italy.(2) Their compositions encompass nearly the complete range of coupled Al ↔ Fe 3+ substitution from X Fe 0.06 to 0.33 (Fe 3+ = 0.185-0.967apfu), and in the Mofete wells, X Fe tended to decrease with increasing temperature and depth.(3) The measured well temperatures ranged from 275-350 • C and 300-390 • C in the Mofete and San Vito wells, respectively, although the fluid inclusion study by De Vivo et al. (1989) indicates some differences in past temperatures.(4) Oscillatory, sector, and complex zoning are ubiquitous and indicate rapid changes in one or more intensive parameters.(5) Two members of the allanite group, allanite-(Ce) and ferriallanite-(Ce), were identified; their texture indicates that they are authigenic and that the geothermal fluids were enriched in, at least, LREEs.REE-bearing epidote and clinozoisite were observed, usually in proximity to the allanite group.
Data availability.All data are provided in the tables and in the Supplement.
Author contributions.BDV obtained the samples and contributed to the discussion.HEB examined and described petrographically the samples with the SEM, analyzed the phases with electron microprobe, interpreted the data, and wrote the text.
Competing interests.The contact author has declared that neither of the authors has any competing interests.
Disclaimer.Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Depth
Figure and spot Figure and spot

Fig
Fig.7dgives the X Fe values measured by an EPMA at selected points; Table3gives their complete analysis.Complex zoning was the most common type observed in both Mofete and San Vito samples and was usually mixed with other zoning types.It was recognized in crystals growing in vugs (Fig.8a), in masses of crystals, and in veins (Fig.4f).Figure8ashows a portion of a vug (MF2-1824) with some of the crystals with britholite-(Ce) centers.The enlarged portion (Fig.8b) of a euhedral crystal shows a chaotic assemblage of zones in the center that are terminated by os- ).Where measured by the EPMA, the different sectors had significantly different X Fe values and BSE responses.The zoning https://doi.org/10.5194/ejm-35-25-2023Eur.J. Mineral., 35, 25-44, 2023

Figure 5 .
Figure 5. Co-variation in Al and Fe 3+ (apfu) in the epidote group for the three Mofete cores (a) and the two San Vito cores (b).Epidote (ep) is shown as a solid circle and clinozoisite (czo) as a solid square.

Figure 7 .
Figure 7. BSE images: (a) several millimeter-long, complexly zoned epidote crystals from sample SV3-2353.The left arrow indicates oscillatory growth in the b-crystallographic direction of the euhedral elongate crystal.The core of the crystal (right arrow) has an irregularly shaped iron-rich core.The line X-Y is the traverse shown in Fig. 7c.(b) A basal section of a complexly zoned epidote crystal showing an iron-rich core (arrow) that has been partially resorbed.The dark areas in the core are biotite crystals.(c) Enlarged detail of the EPMA traverse shown in Fig. 7a.The EPMA line traverses for the elements Fe and Al are shown below registered in position.The Fe-rich bands are epidote (ep), and Al-rich bands are clinozoisite (czo).The conditions of the EPMA line traverse were 15 keV and 30 nA, with a pixel size of ∼ 0.1 µm.The excitation of the sample with the electron beam is a tear-drop-shaped volume so that as the beam traverses from one section to another, the boundaries are sloped.(d) Values of X Fe [(Fe 3+ ) / (Fe 3++ Al)] apfu for the EPMA analytical spots A-F for epidote (ep) and clinozoisite (czo) shown in Fig.7a, b, and c; analyses are given in Table3.

Depth
Figure and spot

8. 4
REE-bearing epidote supergroup and REEs in geothermal fluidsAllanite-(Ce) was commonly observed in the Mofete samples but was less abundant in the San Vito samples.Ferriallanite was only identified in one Mofete sample, MF1-1597.The distribution and abundance of the REE-bearing epidote group (epidote and clinozoisite) follow those of the allanite group.During SEM-EDS petrographic study we identified the following geothermal authigenic REE-bearing minerals: zircon, zirconolite, synchysite-(Ce), bastnäsite-(Ce), monazite-(Ce), cerite-(Ce), aluminocerite-(Ce), britholite-(Ce), and probable hellandite-(Y); all except zircon and zirconolite were LREE enriched.We also noted the presence of authigenic fluorite and fluorapatite in the samples, attesting to F as a component of the geothermal fluids.Mormone et al. (2011) report allanite in MF1 and a REE phosphate in MF5, and Piochi et al. (2021) identified an La-and Ce-bearing epidote in SV1 at 1713 m depth; no analyses of any of these REE phases are given.

Table 1 .
Distribution of epidote supergroup minerals in the Mofete and San Vito samples.
* Measured well temperature after stabilization

Table 2 .
Representative compositions of epidote and clinozoisite from the Mofete samples.