Article body

INTRODUCTION

Intrusion-related gold systems have recently raised great interest among economic geologists, leading to many new discoveries worldwide (McCoy et al. 1997; Thompson et al. 1999; Lang et al. 2000; Lang and Baker 2001; Fan et al. 2003). In southwestern New Brunswick, part of the Canadian Appalachian orogen, several gold deposits and occurrences (e.g., Clarence Stream, Lake George, and Poplar Mountain) have been found to share some similarities with intrusion-related gold systems elsewhere in the world (McLeod and McCutcheon 2000; Thorne and Lentz 2001; Chi 2002; Lentz et al. 2002; McLeod 2002 ; McLeod and Fyffe 2002; Thorne et al. 2002; Yang et al. 2002a, 2002b; Thorne and Lentz 2003). The model of intrusion-related gold systems has stimulated exploration efforts for gold in the area and has highlighted the need for studies to determine the essential controlling factors for formation of these deposits. Clearly, metal behaviour during magmatic evolution is one of the most important factors controlling mineralization.

A fundamental problem with granitoid intrusions is that the precursor magmas have traveled vertically (or laterally) from their source regions to the upper crustal environment, and hence experienced various geological processes that may have greatly modified their original geochemical signatures; therefore, a comprehensive study involving a multidisciplinary approach is required to deal with the complex problem of their petrogenesis (Clarke 1992). In order to better understand the petrogenesis of a granitoid complex and/or batholith, the temporal and spatial framework of the complex must be established from field relations and geochronological data.

The Late Devonian Mount Pleasant Caldera Complex contains two granitic suites (Kooiman et al. 1986; Sinclair et al. 1988; Sinclair and Kooiman 1990; McCutcheon 1990a, 1990b; McCutcheon et al. 1997, 2001). The older McDougall Brook Porphyritic Microgranite of McCutcheon (1990a) is associated with the early stage of caldera development, whereas the younger Mount Pleasant Porphyry of McCutcheon (1990a) is associated with the final stage of the caldera evolution. These two suites are termed the McDougall Brook Granitic Suite (MBG) and the Mount Pleasant Granitic Suite (MPG) in this study (Fig. 1), because both consist of multiple phases with intrusive contacts. The well-known W-Mo-Bi and Sn-Cu-Zn-Pb-In deposits at Mount Pleasant are genetically related to various phases of the MPG (Kooiman et al. 1986; Lentz et al. 1988; Sinclair et al. 1988; McCutcheon 1990b; Sinclair and Kooiman 1990; Inverno 1991; McCutcheon et al. 2001). Interestingly, auriferous greisenized hydrothermal quartz breccias and veins (200 ppb up to 1 g/t Au, and up to 27 g/t Ag; see Johnson 2003) cut the MBG and its wall-rock, both of which were intruded by the MPG (Figs. 1 and 2).

Previous studies mainly dealt with the Mount Pleasant Granitic Suite (i.e. GRI, GRII, and GRIII) and related mineralization characteristics (see Kooiman et al. 1986; Sinclair et al. 1988; Sinclair and Kooiman 1990; Inverno 1991). In these papers, few complete geochemical data sets on the granites were available. Taylor (1992) reported geochemical data for one sample taken from the MPG, and compared it with the Pleasant Ridge topaz granite to the southwest. Whalen et al. (1996) reported geochemical data for one sample from each of the MPG and MBG suites, and noted geochemical differences and similarities between them. No gold data were reported in these studies. Inverno (1991) concentrated on tin mineralization in the MPG, but did not investigate the MBG. McCutcheon (1990a) and McCutcheon et al. (1997, 2001) focused on the stratigraphy and volcanic sequence correlation within the caldera. He and his coworkers analyzed samples from different phases of the MBG for major elements and selected trace elements as a part of a study on the chemical evolution of magmas within the Mount Pleasant caldera (McCutcheon 1990a; McCutcheon et al. 1997). Again, no gold data were reported in these studies, and no work was done on the petrogenetic linkages between the MPG and MBG, and their relations to gold mineralization.

This study focuses on the petrochemical aspects of these two subvolcanic granitic suites within the Mount Pleasant Caldera in order to determine the essential factors that made one suite produce gold mineralization and the other Sn-W-Mo-polymetallic mineralization. On the basis of field relations, petrography, petrochemistry, and trace element modeling, the MBG and MPG are interpreted to be linked through extreme fractional crystallization and supracrustal contamination; fluid fractionation may have played a role in formation of the late MPG. Metal behaviour during the evolution of these magmatic systems indicates that Au behaved as a compatible element, whereas Sn, W, and Mo were incompatible elements. These results can be used to assess the mineralization potential of granitoid intrusions in the area.

Figure 1

Fig. 1 Schematic geological map of the Mount Pleasant Caldera, southwestern New Brunswick (modified from McCutcheon et al. 2001). Field relations indicate that the Mount Pleasant Granitic Suite (MPG) cross cuts the McDougall Brook Granitic Suite (MBG). Locations of the McDougall Brook gold occurrence and Figure 2 are also shown.

-> See the list of figures

GENERAL GEOLOGY

The Mount Pleasant Caldera (Fig. 1), southwestern New Brunswick, is a northeast-trending elliptical feature with minimum dimensions of 13 by 34 km (McCutcheon 1990a, 1990b). The caldera is sited on polydeformed Ordovician to Silurian turbiditic metasedimentary rocks of the Digdeguash and Flume Ridge formations that are thought to form part of the Gander Zone (Fyffe and Fricker 1987). Post-tectonic Middle Mississippian and Pennsylvanian strata overlie the northern part of the caldera, whereas Late Silurian to Devonian granitic rocks of the Saint George Batholith (McLeod et al. 1988; McLeod 1990) form part of its boundary along the southern margin. Rocks within the caldera constitute the Upper Devonian Piskahegan Group that is divided into exocaldera, intracaldera, and late caldera-fill sequences (McCutcheon 1990a, 1990b; McCutcheon et al. 1997, 2001). The Piskahegan Group is largely composed of bimodal volcanic rocks that have been interpreted to have been derived from two parental magmas: mafic with a mantle source and felsic with a lower crustal source (McCutcheon 1990a; McCutcheon et al. 1997). The Mount Pleasant Caldera formed during the transition from a compressional to a transtensional tectonic setting, i.e. during initial development of the Maritimes Basin (McCutcheon et al. 1997). Tectonostratigraphic terrane analysis and geochemi-cal studies (Fyffe and Fricker 1987; Fyffe et al. 1999; Barr and White 1999; Barr et al. 2002) indicate that subduction processes have not taken place in the region since Late Silurian.

Granitic intrusions within the Late Devonian Mount Pleasant Caldera include the McDougall Brook (MBG) and the slightly younger Mount Pleasant (MPG) granitic suites, in the order of intrusive sequence. The MBG is a polyphase suite mainly consisting of marginal feldspar porphyry (MBG1), fine-grained (groundmass) porphyritic monzogranite (MBG2), and minor fine-grained equigranular to subporphyritic quartz monzonite (MBG3); these phases occur with intrusive contacts in the sequence of MBG1 → MBG2 → MBG3 (McCutcheon 1990a; McCutcheon et al. 1997, 2001). The northern half of the MBG intrudes intracaldera rocks, including Seelys Formation (rhyolite and ash flow tuffs) to the south, Scoullar Mountain Formation (sedimentary breccia, andesitic and rhyolitic ash flow tuffs) to the west, and exocaldera Bailey Rock Rhyolite to the east. The southern half of the MBG mostly intrudes the Little Mount Pleasant Formation (rhyolitic ash flow tuffs), although small offshoots intrude into the pre-caldera, Early Devonian Flume Ridge Formation and intracaldera Scoullar Mountain Formation to the west (McCutcheon 1990a; Fig. 1). According to McCutcheon (1990a), the MBG represents the transition between a mafic magma that was contaminated by crystal cumulates from overlying felsic magma in a composi-tionally zoned, high-level subvolcanic magma chamber.

The polyphase MPG (Figs. 1 and 2) is composed of fine-grained equigranular granite (GRI), aplitic to porphyritic granite (GRII), and fine- to medium-grained equigranular granite (GRIII), although locally porphyritic and pegmatitc varieties of GRIII are also present (Kooiman et al. 1986; Sinclair et al. 1988; McCutcheon et al. 2001; Figs. 2 and 3). These units have intrusive contacts indicating the intrusive sequence GRI → GRII → GRIII (Fig. 2b). The GRI unit is heavily brecci-ated and altered, and is genetically related to early W-Mo-Bi mineralization (Kooiman et al. 1986; Sinclair et al. 1988; McCutcheon et al. 2001). The GRI unit is not included in the present study. The MPG intrudes intracaldera rocks (Scoullar Mountain and Little Mount Pleasant formations) and cuts the MBG (Fig. 2) along the southwestern margin of the caldera (Figs. 1 and 2; Sinclair et al. 1988; McCutcheon et al. 1997, 2001). McCutcheon (1990a) interpreted the MPG as the top of a zoned, high-level felsic magma chamber that crystallized inward in successive stages.

Radiometric ages are consistent with the field relations mentioned above, indicating that the MPG suite is slightly younger than the MBG suite. An 40Ar/39Ar age of ca. 361 Ma was obtained for GRIII (Sinclair and Kooiman 1990) and a U-Pb zircon age of 363.4 ± 1.8 Ma (Tucker et al. 1998) was obtained for the Bailey Rock rhyolite that is considered to be the extrusive equivalent of the McDougall Brook feldspar porphyry (McCutcheon et al. 1997). Conventional K-Ar data and whole-rock Rb-Sr isotopic data yielded less reliable ages of 340– 330 Ma for the MPG (Kooiman et al. 1986). A Rb-Sr isochron age of 368 Ma was obtained for the Piskahegan Group (Anderson 1992). These data imply that volcanism and subvo-lcanic plutonism in the caldera took place during a relatively narrow time interval.

Figure 2

Fig. 2 (a) Simplified geological map of the Mount Pleasant Mine area and (b) cross section A-B through the Fire Tower and North Zones showing the W-Mo-W-Bi and Sn-Cu-Zn-Pb-In deposits and their relations to units ,GRI, GRII, and GRIII, in the order of intrusive sequence, of the Mount Pleasant Granitic Suite (modified from McCutcheon et al. 2001). The least altered samples were taken from drill core. Drill holes shown on this map: 1 - AM96-3; 2 - PRL95-1; 3 - PRL95-2; 4 -LNZ18; 5- NMP89-1; 6 - E8D1; 7 – NMP90-1.

-> See the list of figures

Figure 3

Fig. 3 Photographs of representative drill core samples from the Mount Pleasant Granitic Suite and a typical surface sample from the McDougall Brook Granitic Suite. (a) GRIII: fine-grained, equigranular granite (sample NMP89-1-1849; (b) GRIII: fine- to medium-grained, equigranular to porphyritic granite (sample LNZ18-1843); (c) GRII: fine-grained, equigranular to porphyritic granite (sample NMP-1-1759); (d) porphyritic monzogranite from an outcrop of the McDougall Brook Granitic Suite (sample 33A2, see Fig.1 for location).

-> See the list of figures

PETROGRAPHY AND RELATED MINERALIZATION

Based on modal and mesonormative analyses, most of the 27 samples in this study are monzogranite; only three samples from unit MBG1 and one from unit MGB2 are syenogranite (Fig. 4). All of the samples in this study are from subsolvus granite, based on the presence of both plagioclase and K-feldspar.

Figure 4

Fig. 4 Mesonormative Q-A-P classification of granitoid samples from the Mount Pleasant Caldera. The fields are after Le Maitre et al. (1989): 2 - alkali feldspar granite, 3a - syenogranite, 3b - monzogranite, 4 - granodiorite, 5 - tonalite. Chemical data are listed in Table 1.

-> See the list of figures

McDougall Brook Granitic Suite

Unit MBG2 mainly consists of porphyritic monzogranite that locally grades to syenogranite. The monzogranite is composed of plagioclase, K-feldspar, quartz, hornblende, and biotite in order of decreasing abundance, with minor pyroxene. The phenocryst assemblage is plagioclase, K-feldspar, and quartz (~ 40% total abundance), although a few highly altered hornblende (to chlorite) phenocrysts are evident, in which pyroxene cores were noted. Phenocrystic plagioclase forms euhedral to subhedral tabular crystals, ranging from 1 to 12 mm long, with rectangular cross sections and polysynthetic albite twinning. Glomeroporphyritic plagioclase clusters are common. Plagioclase grains are highly altered, giving them a mottled appearance under both plane and cross-polarized light. The alteration products are mainly saussurite (i.e., fine-grained aggregate of epidote group minerals, albite, prehnite, sericite, and possibly clays). K-feldspar (1– 12 mm in grain size) occurs as subhedral, most commonly rectangular, crystals, but rarely euhedral crystals are observed. Coarse patch perthite is common. K-feldspar is also very altered to sericite and clay minerals. Micrographic texture, in which quartz is intergrown in K-feldspar crystals, is present. The quartz forms small, more or less angular grains that are commonly optically continuous over the entire K-feldspar crystal. Phenocrystic quartz generally is less than 4 mm in size; it commonly occurs as subhedral crystals with partly square to hexagonal cross sections, or as anhedral, equant crystals. Groundmass (generally less than 0.5 mm in grain size) constitutes about 60% of the rock and has a similar mineralogy to the phenocryst assemblage. Accessory minerals are apatite, titanite, allanite, zircon, and opaque minerals (mainly magnetite).

The less voluminous MBG1 phase occurs on the margin or border of the MBG2 unit, and consists of porphyritic mon-zogranite (Fig. 3d) and minor equigranular to subporphyritic (containing plagioclase phenocrysts) fine-grained monzogran-ite. The porphyritic monzogranite contains two populations of subhedral to euhedral feldspar phenocrysts (1– 12 mm) with plagioclase generally more abundant than K-feldspar, constituting 30– 35% of the rock. Quartz phenocrysts (1– 4 mm) are minor (< 2%). The groundmass (generally less than 0.1 mm) includes minerals similar to the phenocrysts. Accessory minerals include apatite, zircon, titanite, and magnetite. Most magnetite in the MBG is probably secondary because anhedral magnetite crystals are intimately associated with chlorite that formed by hydrothermal alteration of pyroxene, hornblende, and biotite.

A zone of greisenized quartz breccia and veins about 1500 m long and several hundred meters wide cuts the MBG and its wall-rocks (Seelys Formation ash flow tuffs). This steeply dipping NW-trending zone (Fig. 1) is notably mineralized and contains significant gold (200 ppb up to 1 g/t Au) and silver (up to 27 g/t Ag) (see Johnson 2003). Gold is typically associated with arsenopyrite and pyrite, although the sulphide abundance is low (< 1%) and unevenly distributed in the quartz breccias and veins. These features are common in intrusion-related gold systems elsewhere (Lang et al. 2000; Lang and Baker 2001).

Mount Pleasant Granitic Suite

Units GRII and GRIII of the MPG are similar in mineralogy, and consist mainly of quartz, K-feldspar, plagioclase, and biotite. They may represent two pulses of magma injection from an inward-cooling magma chamber, as indicated by intrusive contacts (Fig. 2b). Unit GRII is mainly composed of porphyritic monzogranite, with local aplitic and seriate textural varieties. Phenocrysts (30– 40%) consist of subhedral to anhedral quartz (0.5 to 1.5 mm), anhedral K-feldspar (0.5 to 3 mm), subhedral to anhedral plagioclase (0.5 to 1 mm), and biotite (0.5 to 2 mm). The groundmass contains minerals similar to the phenocrysts, but the grain size is commonly less than 0.3 mm.

Unit GRIII mainly consists of fine- to medium-grained equigranular monzogranite, although locally porphyritic and pegmatitc varieties are present. Subhedral to anhedral quartz (30– 40%) is dominant, but euhedral quartz is also observed; grain size ranges from 0.5 to 1.5 mm. K-feldspar usually occurs as anhedral to subhedral crystals from 0.5 to 3.5 mm, which constitute 40– 50% of the rock. Plagioclase (20 to 25%) forms subhedral crystals, 1 to 1.5 mm long. Subhedral biotite flakes, 1 to 3 mm in size, contain pleochroic halos caused by radiation damage induced by zircon inclusions, form 1– 3%. Opaque minerals include ilmenite, pyrite, and pyrrhotite. Magnetite is associated with chlorite, suggesting that it formed by hydro-thermal alteration.

Accessory minerals in both GRII and GRIII are monazite, zircon, ilmenite, columbite, and Li mica. Topaz and fluorite in the GRII and GRIII are both interpreted as products of hydro-thermal alteration based on textural relations. However, a few topaz crystals that occur as inclusions in euhedral quartz in GRIII and the presence of primary topaz in an aplitic contact phase of GRII, led Sinclair and Kooiman (1990) to conclude that the MPG can be classified as topaz granite (Manning 1988; Manning and Hill 1990). Magmatic topaz also occurs in leucogranite in the vicinity of Pleasant Ridge and Kedron, as reported by Taylor et al. (1985) and Taylor (1992), who termed the Kedron granite the Bonny River granite. These granite bodies are similar in composition and age to the MPG.

Parts of the aplitic to porphyritic phase of the GRII contain miarolitic cavities, comb quartz layers, and unidirectional solidification textures (UST), suggesting fluid saturation and/ or undercooling of the parental magma (Kirkham and Sinclair 1988). Micrographic intergrowths of quartz and K-feldspar are common in the granite associated with comb quartz layers, both of which reflect a significantly undercooled magma (Fenn 1986; Lentz and Fowler 1992). Detailed drilling and mapping of underground developments indicate that major Sn-Zn-Pb-Cu-In mineralization is closely associated with unit GRII (Kooiman et al. 1986; Sinclair et al. 1988; Sinclair and Kooiman 1990; McCutcheon et al. 2001; Fig. 2b).

In the North Zone at Mount Pleasant, unit GRIII clearly intrudes GRI and GRII (Fig. 2b), which is indicated by sharp contacts that in many places are marked by thin layers of UST (mainly K-feldspar) in the GRIII. Miarolitic cavities (1 mm to 2 cm in size) are locally abundant, in which very fine-grained ser-icite, chlorite, fluorite, and carbonate are present, but arseno-pyrite, molybdenite, and other metallic minerals are relatively rare (Sinclair et al. 1988). Despite well preserved fluid saturation textures in GRIII and evidence of hydrothermal activity, only minor tin mineralization is known to be associated with it (Kooiman et al. 1986; Sinclair et al. 1988; Sinclair and Kooiman 1990; McCutcheon et al. 2001). However, further exploration efforts are needed to confirm this conclusion.

Fluid inclusion studies showed that magmatic fluids were predominant in the early stage of mineralization, but more meteoric fluids were involved in late-stage ore-forming processes (Davis and Williams-Jones 1985; Samson 1990). High salinity (>30 wt.% NaCl equivalent) and high temperature (350– 490°C) fluid is interpreted to have been derived from resurgent boiling, possibly related to magmatic fluids exsolved from the cooling magmas; dilution of this early fluid by con-vecting meteoric water resulted in low- to moderate-salinity fluids that dominated the inclusion population (Davis and Williams-Jones 1985). Low temperature (180– 250°C) and high-salinity fluid may represent late residual fluid from a change in the pressure regime from dominantly lithostatic to hydrostatic conditions (Davis and Williams-Jones 1985). Clearly, hydrothermal breccias, related pervasive alteration and mineralization (silicification, chloritization, and greiseni-zation) in GRI and GRII (Kooiman et al. 1986) reflect pressure fluctuations of fluids associated with the granitic magmas in the Mount Pleasant deposits. These processes are similar to those that operated in the adjacent True Hill biotite porphy-ritic granite intrusions (Lentz et al. 1988; Lentz and McAllister 1990; Lentz 1994; Lentz and Gregoire, 1995), suggesting a comparable genesis for all Sn-W-Mo-Bi mineral deposits in the region.

GEOCHEMISTRY

Analytical procedures

Ten least-altered granite samples collected from seven drill cores at the Mount Pleasant mine and 17 samples from outcrops of the McDougall Brook Granitic Suite (see Figs. 1 and 2 for sample locations) were analyzed for major elements and selected trace elements, including rare earth elements (REE). The methods used in this study include inductively coupled plasma - optical emission spectrometry (ICP-OES), inductively coupled plasma - mass spectrometry (ICP-MS), instrumental neutron activation analysis (INAA), and X-Ray fluorescence spectrometry (XRF).

The ten samples from the MPG were systematically analyzed with ICP-OES (SiO2, TiO2, Al2O3, Fe2O3(total), MnO, CaO, MgO, Na2O, K2O, P2O5, S, As, Ba, Cd, Co, Cr, Cu, Ni, Sc, Sr, V, and Zn) and ICP-MS (Rb, Y, Zr, Nb, Cs, La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, Hf, Ta, Pb, Th, U, Bi, and Sn) at the GSC-Quebec (Geological Survey of Canada-Quebec), with INAA (Au, Ag, As, Ba, Br, Ca, Co, Cr, Cs, Fe, Hf, Hg, Ir, Mo, Na, Ni, Rb, Sb, Sc, Se, Sn, Sr, Ta, Th, U, W, Zn, La, Ce, Nd, Sm, Eu, Tb, Yb, and Lu) at the Activation Labs, Ancaster, Ontario, and with XRF (SiO2, TiO2, Al2O3, Fe2O3(total), MnO, CaO, MgO, Na2O, K2O, P2O5, Rb, Ba, Sr, Ga, Nb, Zr, Y, Th, U, Ce, Cr, Ni, Sc, V, Cu, Pb, Zn, Cl, S, and As) on pressed pellets (Longerich 1995) at the Memorial University of Newfoundland, St. John's. The seventeen samples from the MBG were analyzed with INAA for selected trace elements including REE, also at the Activation Labs; major element data for these samples were taken from the original dataset of McCutcheon (1990a), which were determined by XRF. Certified standard samples SY2 and JG1 (Govindaraju 1994) were used to monitor data quality. Analytical precision of major elements is better than 3 % relative standard deviation (RSD), except for SiO2 (5% RSD); and for most trace elements is better than 10% RSD. Elements for which the analyzed values were beyond 5% of the recommended values in the standards were not used. The analytical results that satisfied the requirements of precision and accuracy are shown in Table 1, and are considered to best represent the chemical compositions of the samples.

Table 1

Table 1. Analyses of rocks from the Mount Pleasant Granitic Suite, Mount Pleasant Caldera, southwestern New Brunswick

Table 1. Analyses of rocks from the Mount Pleasant Granitic Suite, Mount Pleasant Caldera, southwestern New Brunswick

-> See the list of tables

Major elements

Although the MBG has low silica content (SiO2 < 70 wt. %) relative to the MPG (SiO2 > 74 wt. %), both are weakly peralu-minous, with four samples from the MBG being metaluminous (Table 1). Aluminum saturation indices (A/CNK = Al2O3/ (CaO+Na2O+K2O), molar) range from 0.91 to 1.28 (Table 1). Harker diagrams (Figs. 5a, e and i) show that TiO2, MgO, and P2O5 decrease with increasing silica, whereas Al2O3, CaO, K2O, and Na2O show an erratic distribution; Fe2O3 and MnO (Figs. 5c and d) are somewhat less erratic. Samples from the MPG are characterized by high SiO2, (Fe2O3+FeO)/MgO ratio, and alkalis with K2O/Na2O > 1, and low CaO, P2O5, TiO2, and MgO (Table 1; Figs. 5c, e, g, and h), thus sharing many similarities with topaz rhyolite, ongonite, and topaz granite (Christiansen et al. 1983, 1984; Christiansen et al. 1986; Kontak 1990; Taylor 1992; Dostal and Chatterjee 1995; Zhu et al. 2001). The MBG and MPG seem to have two distinctive groups, but each group displays a similar variation trend (Figs. 5a-i) except for those on the plots of A/CNK versus K2O and A/CNK versus CaO (Figs. 5 j and k). These trends may indicate that magmatic evolution was controlled by a similar process such as fractional crystallization, but that minerals involved in the process may have been different in the MBG and MPG.

Figure 5

Fig. 5 Harker diagram for granitoid samples from the Mount Pleasant Caldera. (a) TiO2 vs. SiO2. (b) Al2O3 vs. SiO2. (c) Fe2O3 vs. SiO2. (d) MnO vs. SiO2. (e) MgO vs. SiO2. (f) CaO vs. SiO2. (g) Na2O vs. SiO2. (h) K2O vs. SiO2. (i) P2O5 vs. SiO2 (all data in wt.% in the above plots). (j) A/CNK vs. K2O (wt.%), where A/CNK = Al2O3/( CaO + Na2O + K2O), units in moles. (k) A/CNK vs. CaO (wt.%). (l) Zr (ppm) vs. Nb (ppm); arrows in (l) denote magmatic evolution from early to late stage in the MBG (MBG1 → MBG2) and MPG (GRII → GRIII) suites. See text for detailed discussion.

-> See the list of figures

On the plot of SiO2 versus K2O (Le Maitre et al. 1989), samples from the MBG plot in the fields of high-K and sho-shonite series, whereas those from the MPG fall into the field of medium K (Fig. 6a). Using the MALI (modified alkali-lime index) against SiO2 diagram of Frost et al. (2001), the difference in alkalinity between the MBG and MPG is clear: the former belongs to the alkalic-calcic transitional to the alkalic suite, whereas the latter is calcic-alkalic (Fig. 6b). Interestingly, all samples from the MBG are magnesian granites, defining a trend toward higher FeOtotal/(FeOtotal + MgO) ratios with increasing silica, but the MPG samples range from magnesian to ferroan granite (Fig. 6c), according to criteria established by Frost et al. (2001).

Figure 6

Fig. 6 Petrochemical features of the subvolcanic granitoid samples from the Mount Pleasant Caldera. (a) SiO2 vs. K2O diagram, indicating that these granites are characterized by high K2O (field boundaries are from Le Maitre et al. 1989). (b) Plot of (Na2O+K2O-CaO) (wt.%) vs. SiO2 (wt.%) (field boundaries for alkalic, alkalic-calcic, calcic-alkalic, and calcic-series granitoids are from Frost et al. (2001)). The variable (Na2O+K2O-CaO) is called MALI (modified alkali-lime index) by Frost et al. (2001), based on the alkali-lime classification of Peacock (1931). (c) FeOtot/(FeOtot+MgO) (wt.%) vs. SiO2 (wt.%) diagram. The arrow indicates the fractionation trend, if the magma is evolving under relatively reduced conditions (i.e. no magnetite is involved in the fractional mineral assemblage). The ferroan and magnesian granite boundary is after Frost et al. (2001). (d) ó – ô diagram for granitoids from the Mount Pleasant Caldera. Solid circles are samples from the Mount Pleasant Granitic Suite; solid squares from the McDougall Brook Granitic Suite. Field A represents lavas of volcanoes situated in anorogenic regions, field B volcanoes in orogenic belts and island arcs, and field C includes the alkaline derivatives of both, i.e. trachyte, phonolite, tephrite, etc., among which the sodic types are generally linked to A, and potassic types to B. The tectonic discriminant fields are after Rittmann (1973). The calcic- (ó < 1.2), calc-alkaline- (1.2 – 3.5), alkaline- (3.5 – 8.8), and peralkaline-series (> 8.8) field boundaries are based on empirical studies of a large database of igneous rocks (Qiu 1986; X. Yang, unpublished data).

-> See the list of figures

In terms of the Gottini Index ô (Gottini 1968) and the Rittmann Serial Index ó Rittmann 1957), the lavas of active volcanoes can be divided into two distinctive groups that reflect tectonic setting, i.e. anorogenic and orogenic, as shown in Fig. 6d. The formulations for the parameters ô (Gottini 1968) and ó (Rittmann 1957) are as follows:

Equation 1

-> See the list of figures

Equation 2

-> See the list of figures

This discrimination diagram was originally proposed by Rittmann (1973), based on plotting of Log ó versus Log ô for 1300 lava samples from active volcanoes worldwide. From mafic through intermediate to felsic igneous rocks, the ó value is a very convenient petrochemical indicator, compared to the Peacock index (Peacock 1931), to portray the alkalinity of igneous rocks. Each whole-rock analysis has a defined alkalinity value, if the Rittmann serial index ó is appropriately applied (X. Yang, unpublished data).

Most of the samples from the MBG and MPG suites fall into field B, except for six samples from the MBG (three are located in field C, and the other three fall on the boundary between the fields B and C) (Fig. 6d). Samples from the MPG have ó values of 1.5 to 2.3 (Table 1), i.e. belonging to the calc-alkaline series, except for one sample (PRL95-2-2019) that has ó value of 0.9, which may be attributed to hydrothermal alteration (i.e. chloritization). Their ô values are extremely high, ranging from 182 to 449, consistent with the character of highly differentiated rocks.

Most of the MBG samples are calc-alkaline (ó ranging from 1.8 to 3.3; Table 1), but six samples are weakly alkaline (ó ranging from 3.5 to 4.3), consistent with their high alkali contents (Fig. 6a and b). Their ô values (13 to 68; Table 1) are much lower than those of the MPG samples, indicating that they are less evolved. Notably, an apparent fractionation trend is exhibited on Figure 6c, suggesting that the MPG could be linked through differentiation of the MBG.

A few samples (30A2, 33A2, 193A; Table 1) from the MBG1 phase of the MBG show very high SiO2 (up to 77.5 wt. %) and K2O (up to 8.41 wt.%) contents, but relatively low Na2O (as low as 0.84 wt.%) and Al2O3 (as low as 11.6 wt.%) contents. These features likely reflect hydrothermal alteration (e.g., silicification, sericitization, and chloritization), as indicated by petrographic observation, and are not the result of mag-matic fractionation. Similarly, two samples (PRL95-2-2019 and NMP89-1-1759) from the GRII phase of MPG contain moderately low SiO2 and Na2O (Table 1), which also is attributed to hydrothermal alteration (e.g., chloritization). Consequently, interpretations based on alkalis must be considered with caution. As shown in Kontak and Clark (2002), unaltered and altered granitoids are distinguishable on plots of A/CNK versus K2O and A/CNK versus CaO. However, except for these partly altered samples mentioned above, no apparent alteration trends are evident on these plots (Figs. 5j & k).

Trace and rare earth element (REE) composition

Trace elements, including REE, in the MBG are distinct from those in the MPG (Table 1). The MBG samples have relatively low incompatible element contents, enriched chon-drite-normalized REE patterns [(La/Yb)N = 3.5– 15.7] with more pronounced negative Nb anomalies, but small negative Eu anomalies (Eu/Eu* = 0.43 – 1.09), and more pronounced negative Ti, Sr, and Ba anomalies (Figs. 7 and 8). The flat "bird-wing shape" REE patterns (Fig. 7b) of the MPG samples, which are characterized by large negative Eu anomalies (Eu/Eu* = 0.014 – 0.002) and low (La/Yb)N (ranging from 1.4 to 2.0) ratios, reflect the highly evolved nature of the MPG. These data are consistent with the major element data (e.g., high SiO2 and alkalis, very low CaO, and high ô values). From the MBG to MPG, negative Eu anomalies become increasingly pronounced, with fractionation based on Th enrichment. The MPG samples plot in the A-type granite field, but the MBG samples fall into the I- and S-type field (Fig. 9). If they are associated with fractionation, feldspar must be one of major phases in the crystallizing assemblage (Eby 1992; Christiansen and Keith 1996). It has to be pointed out that highly evolved granites such as the MPG that display some of the geochemi-cal characteristics of A-type intrusions are not necessarily true A-type granites (Whalen et al. 1987; Eby 1990), as no fayalite or sodic hedenbergite that typify "classic" A-type granite occurs in the MPG. Furthermore, petrographic evidence shows that the MPG was saturated with water during crystallization, which also contrasts with the relatively "dry" nature of anorogenic A-type granitoids (Collins et al. 1982; Whalen et al. 1987; Eby 1990).

Relative to I-, S-, A-, and M-type granites (Collins et al. 1982; White and Chappell 1983; Whalen et al. 1987; Eby 1990, 1992), both the MBG and MPG are enriched in F (≤ 0.99 wt.%), Li (≤ 610 ppm) (see Sinclair and Kooiman 1990), Rb (≤ 1210 ppm), Cs (≤ 28 ppm), U (43 ppm), Th (≤ 50 ppm), Nb (≤ 107 ppm), Zr (≤ 580 ppm), and heavy REE (e.g., Yb ≤ 36) (Table 1), which is typical of crustal A-type felsic magmas (Fig. 10), topaz rhyolite (Christiansen et al. 1983, 1984, 1986), ongonite, and topaz granite (Kontak 1990; Taylor 1992; Zhu et al. 2001). Both the MPG and MBG suites are temporally and spatially related, i.e. part of the Mount Pleasant Caldera, and both intrusive suites formed in a within-plate tectonic setting (McCutcheon et al. 1997), but the petrogenetic relationship between them is less clear. It is noteworthy that Nb and Zr behave differently in the MBG and MPG suites. From early to late in time sequence (GRII → GRIII), Nb decreases with decreasing Zr in the MPG, whereas Nb increases with increasing Zr from MBG1 → MBG2 in the MBG (Fig. 5l).

Figure 7

Fig. 7 Chondrite-normalized rare earth element (REE) patterns for granitoids from the Mount Pleasant Caldera. (a) McDougall Brook Granitic Suite samples. (b) Mount Pleasant Granitic Suite samples. Normalizing values for C1 chondrite are from Sun (1982). Symbols are the same as in Figure 4.

-> See the list of figures

Figure 8

Fig. 8 Primitive mantle-normalized spider diagram for granitoids from the Mount Pleasant Caldera. (a) McDougall Brook Granitic Suite. (b) Mount Pleasant Granitic Suite. Primitive mantle values are from Taylor and McLennan (1985). Symbols are the same as in Figure 4.

-> See the list of figures

Figure 9

Fig. 9 Variation in Th (ppm) as a function of the negative Eu anomaly (Eu/Eu*) for granitoid samples from the Mount Pleasant Caldera. The arrow shows the fractional crystallization trend of a magma crystallizing feldspar. The boundary line between the A- and I- and S-type granite fields is from Eby (1992).

-> See the list of figures

Figure 10

Fig. 10 Geochemical discrimination diagrams showing (Zr + Nb + Ce + Y) plotted against FeO*/MgO (a) and (K2O + Na2O)/CaO (b) for granitoid samples from the Mount Pleasant Caldera. A-type granites are located outside of the divided box in each diagram. FG = fractionated felsic granites and OGT = unfractionated M-, I-, and S-type granites; field boundaries are after Whalen et al. (1987).

-> See the list of figures

PETROGENETIC DISCUSSION

The distribution of major and trace elements in the MBG and MPG seems to indicate two distinctive groups with an apparent "gap" (Figs. 5, 11a-c and 12). It would be easy to ascribe this difference to diverse petrogenetic processes for formation of these two granitoid suites within the Mount Pleasant Caldera, without consideration of their field relations, detailed petrographic features, systematic chemical variations, and iso-topic compositions. Brophy (1991) noted that compositional gaps, like the ones above, are not uncommon in calc-alkaline magmatic systems. He proposed that fractionation-generated compositional gaps may result from (1) simultaneous interior crystallization and vigorous chamber-wide convection which leads to progressive crystal suspension, (2) cessation of convection when the percentage of suspended crystals reaches the critical crystallinity of the magma, and (3) eventual buoyancy-driven crystal-melt segregation producing a discrete body of fractionated magma that is separated from the initial magma by a compositional gap. Marsh (1996) pointed out that in a magma chamber with phenocrysts, differentiation is controlled by a balance between capture of crystals by the solidification front and crystal settling; massive crystal separation results in less viscous residual magma with discontinuous compositions. Solidification front instability causes silicic melt segregation and produces a separate magma (Marsh 1996).

Figure 11

Fig. 11 Tectonomagmatic discrimination diagrams for granitoid samples from the Mount Pleasant Caldera. Log– log plots of (a) (Y + Nb) versus Rb and (b) Y versus Nb . Field boundaries are from Pearce et al. (1984) as modified by Christiansen and Keith (1996). (c) Triangle plot of Y– Nb– Ce (A1 Group granites are characterized by element ratios similar to the mantle, whereas A2 Group granites originated from continental crust or arcs; boundary line between groups is from Eby (1992). (d) Diagram showing Yb versus Th/Ta; WPVZ = within plate volcanic zones, field boundary is from Gorton and Schandl (2000). Symbols are the same as in Figure 4.

-> See the list of figures

Figure 12

Fig. 12 Representative logarithmic variation diagrams (ppm) for granitoid samples from the Mount Pleasant Caldera. The bulk partition coefficient (D) for each element can be calculated from the equation D = 1 - m, where m is the slope of the line based on least-square regression analysis. R is the correlation coefficient, which is statistically significant at a confidence level of 95%. Symbols are the same as in Figure 4.

-> See the list of figures

Isotopic compositions

Published isotopic data suggest that the MBG and MPG have a similar source, as they have the same åNd(t) (-0.2) and similar ä18O (7.1 and 8.2 ~ 8.6‰) values (see Taylor 1992; Whalen et al. 1996). The slightly lower ä18O value of the MBG may be attributed to meteoric water interaction (Whalen et al. 1996). These characteristics suggest that the MBG and MPG either were derived from the same parent magma or from different batches of partial melting of the same protolith. The initial 87Sr/86Sr ratio for the MPG is 0.7126 (Kooiman et al. 1986), which allows for a crustal source dominated by recycled crustal materials or a mantle source strongly contaminated by supracrustal rocks. The initial 87Sr/86Sr ratio for the Piskahegan Group is 0.7061 (Anderson 1992), requiring some input from the mantle, consistent with the bimodal nature of the caldera rocks (McCutcheon 1990a; McCutcheon et al. 1997, 2001). A mantle-derived mafic magma could cause partial melting of crustal rocks and generate the felsic magmas responsible for the MBG and MPG. Given low Sr and high initial 87Sr/86Sr ratio together with åNd(t) and ä18O values of the MPG rocks, the most likely scenario is that the magma originated from partial melting of juvenile materials and was subsequently contaminated by supracrustal rocks. Based on a systematic study of geochemistry, Taylor (1992) reached the same conclusion for the source material of the Pleasant Ridge topaz granite that has many chemical similarities to the MPG.

Petrogenetic linkage between the MBG and MPG

Two models are commonly used to explain topaz-bearing granitic rocks: very small degrees of partial melting of continental crustal material that was subjected to an earlier episode of melting (Christiansen et al. 1986; Manning and Hill 1990) and extensive fractionation of felsic magma (Dostal and Chatterjee 1995, 2000; Zhu et al. 2001). Generally, the magmas are generated in an extensional tectonic regime (Christiansen et al. 1986; Whalen et al. 1987, 1996). Moreover, extensive liquid fractionation involving thermogravitational diffusion, liquid immiscibility, and vapour-phase transport has also been suggested to explain the petrochemical variations and the substantial enrichment of incompatible trace elements (Christiansen et al. 1983, 1984, 1986; Dostal and Chatterjee 1995, 2000; Zhu et al. 2001).

Large variations in compatible elements, such as Sr, Ba, and Eu in the MBG and MPG (Table 1, Figs. 7 and 8), rule out the possibility of variable degrees of partial melting of the same source, although small degrees of partial melting can generate high contents of incompatible elements (see Rollinson 1993; Christiansen and Keith 1996). Fractional crystallization, however, not only reduces the abundance of compatible elements, but also increases the concentrations of incompatible elements (e.g., Li, Rb, Cs, U, Th, Nb, and Ta), without any variations in the isotopic ratios if the effect of supracrustal contamination is not considered.

The petrochemical data presented in this study (Table 1; Figs. 5– 9, 11, 12) are modeled assuming that the MBG and MPG were produced by extensive low-T fractional crystallization of the same parental magma, rather than from two batch-melts that evolved separately. The extreme fractional crystallization is interpreted to have taken place mainly at depth, as indicated by two separate compositional groups, but with the same evolutionary trends for the MBG and MPG (Figs. 5, 6c & d, 9, 11, and 12). However, supracrustal contamination of selective components is also reflected in some major elements (Fig. 5) and high initial 87Sr/86Sr ratio, as well as trace elements such as Rb and Cs (Fig. 12; Table 1). The erratic distribution of the alkali elements is interpreted as the manifestation of contamination or alteration, whereas the tightly constrained distribution of TiO2, P2O5, and MgO is interpreted as primary eatures. No alteration trends are apparent on Figures 5j, k (cf. Kontak et al. 2002). Figure 13 illustrates the log-log plot of incompatible (Nb) versus compatible (Sr) elements for the samples from the MBG and MPG, in which the samples are mainly confined in the field between fractional crystallization (FC) and assimilation-fractional crystallization (AFC) (r = 0.5) curves. Assuming a leucogabbro (Thorne and Lentz 2002) magma as the starting point and using upper continental crust (Taylor and McLennan 1985) as contaminant, the FC and AFC paths (Fig. 13) can be constructed with the equations of Rollinson (1993), given the bulk partition coefficients (DSr and DNb) estimated by the method of Allégre et al. (1977) (see details below). The parameters used in the modeling are given in the caption of Fig. 13. General vectors produced by Rayleigh fractionation of single phases (e.g., plagioclase, K-feldspar, bio-tite, hornblende, and clinopyroxene) and their effects on melt composition are also shown for comparison. The compositions of average N-type MORB (Saunders and Tarney 1984), average OIB (Sun 1980), lower continental crust (Weaver and Tarney 1984), and upper continental crust (Taylor and McLennan 1985) are plotted on Fig. 13, which constrains the probable sources for the granitoids from the Mount Pleasant Caldera. All these data suggest that magma evolution was controlled by extreme fractional crystallization plus selective supracrustal contamination. Figure 6c indicates that fractional crystallization may have proceeded under relatively reduced conditions, as high ƒ(O2) would have resulted in crystallization of magnetite and depleted the residual magma in iron.

A simplified calculation was performed to test the probability of petrogenetic connection between the MBG and MPG through fractional crystallization. Preliminary trace-element modeling, together with major-element petrochemical evidence, indicates that the mineral assemblage involved in fractional crystallization was clinopyroxene, amphibole, pla-gioclase, K-feldspar, biotite, zircon, and apatite (Fig. 13). Some minerals rich in light REE (e.g., monazite and allanite) may also occur as separate phases, although only a very small proportion is required. The most evolved sample (NMP89-1-1849) from the MPG is equivalent to 10– 15% residual melt (for most incompatible elements, such as Nb, Ta, and HREE) generated by the most primitive sample (295) from the MBG that is assumed to represent the composition of the parental magma in terms of a fractional crystallization model (Christiansen and Keith 1996) and the estimated bulk partition coefficient (D) as described below. Inconsistency with a lower degree of frac-tionation for incompatible elements (e.g., Rb) is interpreted to imply that the enrichment of the high field strength elements, such as Nb, Ta, and HREE, may require fluid involvement. Dostal and Chatterjee (1995, 2000) showed that fluid fraction-ation could have resulted in differential enrichment of these elements to the level observed in ca. 370 Ma topaz leucogranite from South Mountain Batholith in Nova Scotia. Compared to their data, the present study indicates that fluid fractionation may have played an important role in the extreme enrichment of these elements (e.g., Nb, Ta, and HREE) in the MPG (see discussion below).

This model of fractional crystallization coupled with fluid fractionation plus supracrustal contamination proposed in this study (Fig. 14) differs from the previous model by McCutcheon (1990a), wherein he suggested that two different magmas (i.e. two stages of partial melting of lower crust) evolved separately in high-level chambers by AFC processes to give rise to the different phases. McLeod (1990) proposed that highly evolved granites like the MPG may be derived from partial melting of a F-rich source (crust) at depth, which subsequently fractionated repeatedly in a higher level subvolcanic chamber. In that model, the origin of the MBG and its relation to the Mount Pleasant granites were not discussed.

Figure 13

Fig. 13 Log-log plot of Nb (ppm) versus Sr (ppm) for granitoid samples from the Mount Pleasant Caldera. Symbols are as those in Fig. 4. Start point is assumed to be: Nb = 4 ppm, Sr = 500 ppm, equivalent to Late Silurian leucogabbro from the Saint George Batholith (Thorne and Lentz 2002). The upper continental crust (Taylor and McLennan 1985) is used as the contaminant (Nb =25 ppm, Sr = 350 ppm) in AFC modeling. Bulk coefficients are calculated based on the method of Allégre et al. (1977) for processing the dataset (Table 1): DSr = 1.81, DNb = 0.22; see text for detailed discussion. General trajectories for fractional crystallization (FC) and for assimilation-fractional crystallization (AFC) are shown for comparison, with marks of residual magma fraction. The equations of FC and AFC are taken from Rollinson (1993) for the calculations. Ratio of mass assimilated to mass fractionation (r) was simulated for 0.2 and 0.5, which is shown for the AFC curves. General theoretical vectors produced by Rayleigh fractionation of single phases and their effect on melt composition are also shown for reference: plagioclase (Pl), K-feldspar (K-fld), biotite (Bt), hornblende (Hb), and clinopyroxene (Cpx); mineral-melt partition coefficients are from the compilations of Rollinson (1993). U - upper continental crust (Taylor and McLennan, 1985) (Nb = 25 ppm, Sr = 350 ppm); L - lower continental crust (Weaver and Tarney 1984) (Nb = 5 ppm, Sr = 569 ppm); M - average N-type MORB (Saunders and Tarney 1984) (Nb = 2.5 ppm, Sr = 136 ppm); I - average OIB (Sun 1980) (Nb = 53 ppm, Sr = 800 ppm).

-> See the list of figures

Figure 14

Fig. 14 A hypothetical model for the origin of the granitic suites in the Mount Pleasant Caldera (modified from McCutcheon et al. 1997). Fractional crystallization took place in the magma chamber at depth, cooling and crystallizing inward to the center of the chamber and resulting in zoning in composition and physiochemical variables. The marginal zone of the chamber is more contaminated by supracrustal rocks. Sequential pulses of magma injection were emplaced into the upper crustal rocks (the Piskahegan Group within the caldera). High F, Li and H2O contents in the magma enhanced the protracted crystallization, resulting in extreme enrichment in some incompatible elements (e.g., Rb, Li, Cs, Nb, Ta, HREE, Sn, and W), and depletion in compatible elements such as Eu, Sr, Ba, Ti, and P, Au, and Zr (Hf) in the late MPG.

-> See the list of figures

Metal behaviour during fractional crystallization

In terms of the Rayleigh fractionation law, the partition coefficient D of a hygromagmatophile element such as U can be assumed to be zero in order to evaluate geochemical behaviour of the other elements in a fractionating magmatic system. An equation for this relation is deduced based on this assumption (Allégre et al. 1977):

Equation 3

-> See the list of figures

where CL and Co denote the concentration of an element of interest in residual melt and parental melt, respectively, D is the bulk partition coefficient of the element between the crystallizing mineral assemblage and the residual melt, C*L and C*o represent the content of a hygromagmatophile element in the residual and parental melts, respectively, whose D* is supposed to equal to zero. Therefore, plotting ln CL against lnC*L gives a straight line, given D remains constant during fractional crystallization (Allégre et al. 1977; Christiansen et al. 1983; Fowler and Jensen 1989). Using least-square regression, the slope (m) of the regressed line can be obtained (Fig. 12). Hence, the bulk partition coefficient (D) can be determined for that element (i.e. m = 1 – D).

The D calculation results indicate that Au (1.14) together with Sr (1.81), Ba (2.73), Hf (1.27), Zr (1.53), La (1.20), Ce (1.12), Pr (1.31), Eu (3.29), Sc (1.65), Zn (1.21), and Pb (1.24) were compatible elements during magmatic evolution, suggesting that there is less chance to find Au mineralization in the evolved MPG. Thus, it is better to look for gold within and/or around the early MBG. Based on their D values, Rb (0.34), Cs (0.80), Ga (0.70), Ta (0.16), Nb (0.22), Th (0.56), Y (0.22), Sm (0.83), Gd (0.6), Yb (0.21), Lu (0.23), Sn (-0.18), W (0.96), Cu (0.52), and Mo (0.60) behaved as incompatible elements during magmatic evolution.

The technique of Allégre et al. (1977) applied to granitoid systems, such as those in the Mount Pleasant Caldera, may not be interpreted easily. Firstly, the mobility of U in such subvolca-nic systems must be considered during late- to post-magmatic hydrothermal alteration processes, although a large proportion of U is hosted in zircon in the MBG and MPG. Incompatibility of Th in a granitoid system is fairly similar to that of U, but Th is considered to be immobile during alteration processes (see Fowler and Jensen 1989; Jenner 1996). Therefore, Th is a better fractionation index than U. However, using Th as the frac-tionation index to conduct the same modeling as U described above gave the same tendency of compatibility for all elements listed above. Secondly, the fact that some elements such as Sn, Sb, Cu, and Au (see Jenner 1996) would partition into volatile phases during crystallization of the magma, must be also considered in the trace-element simulation. This consideration suggests that a significant fraction of the metals enters volatile phases that reduce their abundance in the degassed magmas. Thus, compatibility of these elements modeled in this empirical fashion (Allégre et al. 1977) may not really reflect their behaviour in the crystal-melt system alone. It is noticeable that the bulk partition coefficient of Sn is negative (-0.18). This number itself is meaningless, although it is smaller than 1. Thirdly, alteration may elevate or decrease some elements in altered granitoids, which also influences the modeling results. All these effects must be carefully evaluated before a conclusive determination of D for a metal during evolution of a granitic magma system as described here can be achieved. Then, these data for estimated metal bulk partition coefficients can be applied to assess the mineralization potential.

It is notable that the MPG has very low K/Rb ratios (26– 56), requiring fluid involvement in formation of the rocks (Shaw 1970; Dostal and Chatterjee 1995, 2000). Conversely, the MBG has high K/Rb ratios (172– 376), which are typical for magmatic processes. Nb/Ta ratios are lower for the MPG (2.7– 7.6) compared to the MBG (7.7– 28.6), also requiring fluid involvement in the late-stage granites. Intriguingly, the MPG has relatively constant Zr/Hf ratios (12.8-17.1), whereas the MBG shows remarkably higher Zr/Hf ratios (24.3– 48.0), suggesting that Zr/Hf ratios in the residual melts were controlled mainly by fractional crystallization of zircon (Linnen and Keppler 2002), or by fluid fractionation (Dostal and Chatterjee 1995, 2000). The Th/U ratios are lower in the MPG (1.1– 1.8) and relatively high in the MBG (3.0– 7.1), suggesting that U increases more rapidly than Th, consistent with the above estimation of the bulk partition coefficient (D). In other words, U was more incompatible than Th in the magma system. It is known that Nb and Ta are strongly affected by the presence of volatiles and can be considered with Li, B, F, Rb, and Sn as volatile-associated elements (Manning and Hill 1990). In contrast, Zr and Hf are considered to be more refractory, in view of the low solubility of their accessory minerals in peraluminous granitic melts (Watson 1979; Watson and Harrison 1983). P, Ti, Sr, and Ba mainly behave as compatible elements in the granitic system.

Estimation of liquidus temperatures

Liquidus temperatures of granitic rocks can be estimated from accessory mineral saturation temperatures (Watson and Harrison 1983; Harrison and Watson 1984; Montel 1993). These temperatures are capable of reflecting the mechanism of magmatic evolution. In a fractional crystallization sequence, the liquidus temperatures of successive intrusive phases in a granitoid complex decrease from early to late stage crystallization. On the other hand, the liquidus temperatures of a second stage melt would be higher than a first-stage melt during partial melting, but once formed both evolve by FC related to decreasing temperatures.

Monazite-saturation geothermometry (Montel 1993) yields temperatures as high as 825– 950°C for the MPG, similar to those of the MBG (834– 950°C). Abnormally high estimated temperatures for the MPG, assuming it formed through fractional crystallization discussed above, are attributed to LREE contents (Table 1 and Fig. 7) similar to those of the MBG. Possibly other unidentified minerals host a significant part of the LREE in the bulk rock. If this portion of LREE were removed from the calculation, then the estimated temperatures would have been lowered in terms of the monazite-saturation geothermometry. However, zircon-saturation geothermom-etry (Watson and Harrison 1983) produces different results for the MPG (725– 789°C) and MBG (769– 933°C). Apatite-saturation geothermometry (Harrison and Watson 1984) gives results similar to those from zircon-saturation geother-mometry for the MPG (712– 791°C) and MBG (833– 963°C). Decrease of the estimated liquidus temperatures (Fig. 15) from the MBG to MPG is consistent with dropping temperatures in a fractionating magma sequence.

High F concentration (3500– 9900 ppm,) together with high Li (up to 610 ppm) (see Sinclair and Kooiman 1990) and high alkalis in the water-saturated MPG could remarkably reduce solidus temperatures to as low as 525°C, as shown by experimental work on F-rich (1.2 wt.%) topaz-bearing granitic melts (Webster et al. 1987). The large temperature range for crystallization (950– 525°C) allows protracted extreme fractionation of granitic melts such as the MPG, resulting in unusual enrichment of incompatible elements such as F, Li, Rb, Cs, Y, Nb, Ta, Ga, Sn, and W (Fig. 14).

Figure 15

Fig. 15 Diagram of estimated liquidus temperatures Tzircon versus Tapatite (°C) for granitoid samples from Mount Pleasant Caldera. Symbols are the same as in Figure 4. The arrow points the trend of the temperatures decreasing from the MBG to MPG. Zircon and apatite saturation temperatures (°C) were calculated using Watson and Harrison (1983) and Harrison and Watson (1984) equations, respectively, for the samples listed in Table 1.

-> See the list of figures

Discrimination of tectonic setting

The geochemical features of granitoids may reflect their tectonic environment to some extent (Pearce et al. 1984). Various geochemical discrimination diagrams for tectono-magmatic classification have been proposed to distinguish tectonic settings of magmatic rocks in terms of their chemistry (see Rollinson 1993).

Geochemical discrimination diagrams suggest that the MBG could have formed in volcanic-arc to within-plate tectonic settings, whereas the MPG samples plot entirely in the within-plate field (Figs. 11a and 11b). Figure 11d yields a similar result. As shown above, extensive fractionation plus fluid involvement have strongly enriched high field strength elements, which may shift the MPG samples to the within-plate field (Figs. 11b & d). Supracrustal contamination also affects the discriminant result (Fig. 11a). Therefore, caution must be exercised when using these diagrams to discriminate tectonic setting. These classification diagrams are suspect (cf. Barker et al. 1992), because they characterize the tectono-magmatic environment under which the granite protolith formed, not necessarily the environment prevailing during granite emplacement (Whalen et al. 1996; Christiansen and Keith 1996). Thus, the trace element signatures of granitoids directly reflect their protolith, melting, and crystallization histories, rather than their tectonic settings. Granitoids from the Mount Pleasant Caldera appear to be crustal A2-type (Fig. 11c), inferred to have been derived from continental crust or underplated crust that has been through a cycle of continent-continent collision or island-arc magmatism (Eby 1990, 1992). Particularly, the MPG is similar in geochemistry to topaz rhyolite (Christiansen et al. 1983, 1984; Christiansen et al. 1986), ongonite, and topaz granite (Kontak 1990; Taylor 1992; Dostal and Chatterjee 1995; Zhu et al. 2001), and also to the Late Devonian Mount Douglas Granite in the eastern part of the Saint George Batholith (Lentz et al. 1988; McLeod et al. 1988; McLeod 1990).

During Late Devonian, the studied region was in the early stages of basin extension (McCutcheon et al. 1990a) followed the Acadian Orogeny, consistent with the discriminant results in terms of Fig. 6d and the multicationic diagram of Batchelor and Bowden (1985) (not shown).

Metallogenic implications

The highly evolved MPG is associated with granophile element mineralization and characterized by high fluorine, lithium (see Sinclair and Kooiman 1990), and alkali concentrations, which enhance the low-T fractionation. The presence of Li, B, and F have a significant effect on phase equilibria in the water-saturated granite system (Q-Ab-Or-H2O), reducing liquidus and solidus temperatures and affecting minimum melt compositions (Manning 1981; Webster et al. 1987; Manning and Hill 1990) and decreasing viscosity of the melt (Baker and Vaillancourt 1995), thus promoting the duration of crystallization to low temperatures. On the basis of F alone, the shift in minimum position toward the Ab apex in the Q-Ab-Or system has been taken as evidence for origin of topaz granite, and other F-rich granitic rocks, through fractional crystallization of biotite granite magma (Manning 1982). F enrichment generally leads to a concomitant enrichment in Sn, Nb, and Ta, which may possibly be directly complexed with F in the melt (Manning and Hill 1990; Dostal and Chatterjee 1995). A strongly depolymerized F-rich melt is more capable of hosting incompatible elements relative to polymerized volatile-poor granite (Watson 1979; Watson and Harrison 1983; Linnen and Keppler 1997), irrespective of complexing with F (Manning and Hill 1990). Based on the estimated D value of Sn, the late GRIII appears to have high potential for producing Sn mineralization.

Geochemically, the less differentiated MBG is broadly comparable to the granitoids associated with intrusion-related gold systems elsewhere (McCoy et al. 1997; Thompson et al. 1999; Lang et al. 2000; Lang and Baker 2001; Yang et al. 2002a, b). Granitoid intrusions related to gold mineralization are mostly metaluminous, subalkalic (calc-alkalic) to weakly alkalic with intermediate to felsic compositions, forming simple to multiphase complexes (Lang et al. 2000). Locally, peraluminous intrusions are present, but they are ascribed to late-stage frac-tionation rather than different magmas with distinctive sources (Duncan et al. 1998). According to the summary of Lang et al. (2000), the main-stage, subalkalic (calc-alkalic) intrusions have silica contents between 62 and 72 wt. % SiO2, and late differentiates reach 78 wt. % SiO2. REE patterns display LREE enrichment relative to HREE, with La and Yb concentrations about 100 and 10 times chondrite values, respectively. Eu anomalies are absent to slightly negative. Most intrusions are characterized as reduced I-type granitoids.

Most of the analyzed samples from the MGB and MPG (Table 1) have low Au contents, although the MBG samples appear to contain higher gold (up to 8 ppb). Compared to continental crust (2.5 ppb) and normal sedimentary rocks, such as greywacke (4.8 ppb) and shale (3 ppb) (see Lentz 2003), Au does not appear to be elevated significantly in the granitoids from the Mount Pleasant Caldera. Regionally, the granitoids in the area also contain relatively low Au compared to continental crust (McLeod 1990; Whalen 1993). Available studies of other intrusion-related Au deposits have shown that the content of Au in unaltered granitoids is not significantly different from those of barren granitoid rocks (McCoy et al. 1997). These data suggest that background Au abundance in the granitoids is not an essential parameter controlling associated gold mineralization. This observation is in accord with the conclusion by Tilling et al. (1973) that it is impossible to relate gold abundance in igneous rocks and gold mineralization. The geochemical behaviour of gold during the magmatic-hydro-thermal evolution of granitoid systems may be a key factor responsible for gold mineralization, which is controlled by the physiochemical conditions and composition of the magmas. Local geological setting (e.g. structure and lithology) and interaction between the intrusions and country rocks affect the magmatic-hydrothermal systems, which ultimately places controls on the emplacement of gold deposits.

Gold is not expected to concentrate in the residual melt during protracted fractionation processes, because it occurs as a compatible metal. Gold is preferentially partitioned into sulphides in equilibrium with a reduced magma, based on empirical observation (McCoy et al. 1997) and experimental studies (Jugo et al. 1999; Kesler et al. 2002). However, early crystallization of magnetite did not appear to severely influence the ability of a calc-alkaline granitic magma to yield a gold-rich ore fluid, according to experimental study of Au partitioning between magnetite and granitic magma (Simon et al. 2003). Consequently, any sulphide saturation in magma systems must greatly deplete Au in the silicate melts. Inverno (1991) reported much higher Au values in his granite samples from the Mount Pleasant Mine, e.g., Au in the GRII ranging from 8 to 85 ppb (average = 28.7 ± 2.6 ppb, n = 6), and the GRIII ranging from < 5 to 27 ppb (average = 13.6 ± 2.5 ppb, n = 6) (see his appendix 7). These high Au contents may be ascribed to hydrothermal alteration, although the Au concentrations decrease from the GRII to GRIII. Our greisenized granite samples taken from True Hill also show higher Au contents (up to 100 ppb), suggesting that late-stage hydrothermal fluids could scavenge gold that is incorporated in early immiscible sulphides in granitoids (Yang and Lentz 2003). At the present stage, it is not easy to evaluate timing of sulphide saturation in felsic magmatic systems, in particular in granite intrusions that were usually influenced by subsolidus processes. These processes may alter the primary sulphides on one hand, and form hydrothermal sulphides on the other. To date, few convincing methods exist for distinguishing primary and hydrothermal sulphides in granitoids, although studies on petrography, chemical, and sulphur isotopic compositions of the sulphides could provide some clues as to their origin. More work needs to be done to better understand the origin of intrusion-related Au ore deposits.

CONCLUSIONS

Based on the present study, it is suggested that subvolcanic granitoid rocks in the Late Devonian Mount Pleasant Caldera were formed by extreme fractional crystallization of magmas that may have been derived from recycled crustal materials, although the early MBG may have been influenced by local contamination of wall-rocks through selective assimilation. The fractional crystallization is interpreted to have taken place mainly at depth, as indicated by two separate compositional trends for the MBG and MPG. Extensive liquid fractionation may have played an important role in magmatic-hydrother-mal evolution, especially facilitating enrichment of high field strength elements, such as Nb, Ta, Y, HREE, Sn, and W, as well as large ion lithophile elements (e.g. Li, Cs, Rb) in late-stage melts (i.e. the MPG). High volatile (i.e. F, Li and H2O) concentrations in the residual melts enhanced the protracted low-T fractional crystallization. Metal elements behaved distinctively during magmatic evolution. Gold occurred as a compatible element, suggesting that gold and metal mineralization are associated with the early granitoid phases such as the MBG, whereas tin behaved as a very incompatible element so that highly evolved granite phases such as the MPG are more likely to generate tin and related mineralization. Hence, it is suggested that exploration for gold related to granitoid complexes such as those in the Mount Pleasant Caldera should be concentrated around the early phase granitoids. In contrast, the late, highly evolved granite may be a good target for exploration for Sn-W-Mo-Bi deposits. It is noteworthy that the youngest phase, GRIII, of the MPG seems to have highest potential for forming Sn deposits.

Some questions remain to be answered, in particular the relationship of the granitoids in the Mount Pleasant Caldera to those in the Saint George Batholith and the Pomeroy Intrusive Suite (McLeod 1990). More detailed petrochemical investigations and reliable age determinations (e.g. U-Pb zircon or monazite) are required to characterize the petroge-netic linkages among these granitoids, although they generally have been thought to be related temporally, spatially, and thus genetically (Lentz et al. 1988; Lentz and McAllister 1990; McLeod 1990).