Geochemical behaviour of trace elements during fractional crystallization and crustal assimilation of the felsic alkaline magmas of the state of Rio de Janeiro , Brazil

This paper presents geochemical behaviour of trace elements of the felsic alkaline rocks of the state of Rio de Janeiro, Brazil, with special attention of fractional crystallization and continental crust assimilation. Fractionation of leucite and K-feldspar increases Rb/K and decreases K2O/(K2O+Na2O). Primitive nepheline syenite magmas have low Zr/TiO2, Sr, and Ba. On the Nb/Y vs. Zr/TiO2 diagram, these rocks are projected on the field of alkaline basalt, basanite, and nephelinite, instead of phonolite. Well-fractionated peralkaline nepheline syenite has high Zr/TiO2 but there are no zircon. The diagrams of silica saturation index (SSI) distinguish the trends originated form fractional crystallization and crustal assimilation. In the field of SSI<-200, Zr/TiO2 and Ba/Sr have negative correlations to SSI in consequence of fractional crystallization. In the field of SSI>-200, they show positive correlations due to continental crust assimilation. Total REEs (Rare Earth Elements) is nearly 10 times that of granitic rocks, but LaN/SmN and LaN/YbN are similar. REE trend is linear and Eu anomaly is irrelevant. The pegmatitic liquid generated by country rock partial melting is SiO2-oversaturated and peraluminous with high Ba, Sr, Ba/Sr, Zr/TiO2, and SSI, with high content of fluids. This model justifies the peraluminous and SiO2-oversaturated composition of the rocks with relevant effects of continental crust assimilation. key words: nepheline syenite, alkaline syenite, trace elements, rare earth elements, fractional crystallization, continental crust assimilation. Correspondence to: Kenji Freire Motoki E-mail: kenji_dl@hotmail.com *In memoriam intRoDuction Geochemical behaviour of trace elements are sometimes different from major elements. Absolute and relative abundance of determined trace elements, such as Zr, Y, Nb, Ga, and Sc, are stable during alteration processes of metamorphism, hydrothermalism and weathering. The behaviour of immobile elements are sometimes related to those of major elements. For example, Zr/Ti is related to SiO2 and Nb/Y is related to alkaline elements (e.g. Floyd and Winchester 1975, Winchester and Floyd 1977). The content can be used as a proxy for K2O and Co content for SiO2 (Hastie et al. 2007). Immobile elements, especially high field strength elements (HFSE, e.g. Zr, Nb, Hf, Ta, and Ti), are convenient for classification of altered rocks


INTRODUCTION
Geochemical behaviour of trace elements are sometimes different from major elements.Absolute and relative abundance of determined trace elements, such as Zr, Y, Nb, Ga, and Sc, are stable during alteration processes of metamorphism, hydrothermalism and weathering.The behaviour of immobile elements are sometimes related to those of major elements.For example, Zr/Ti is related to SiO 2 and Nb/Y is related to alkaline elements (e.g.Floyd and Winchester 1975, Winchester andFloyd 1977).The content can be used as a proxy for K 2 O and Co content for SiO 2 (Hastie et al. 2007).Immobile elements, especially high field strength elements (HFSE, e.g. Zr, Nb, Hf, Ta, and Ti), are convenient for classification of altered rocks An Acad Bras Cienc (2015) 87 (4)
Original rock composition can be estimated even from strongly altered volcanic rock samples and ash-fall deposits (e.g.Caddah et al. 1994, Alves 2005, Kitsopoulos 2010).On the other hand, large ion lithophile elements (LILE, e.g.K, Rb, Sr, and Ba) and the light rare earth elements, especially Ba, are mobile and easily transferred by fluids.In addition, certain elements, such as REEs, show geochemical behaviour different from major elements, providing unique information on magma generation and evolution.Geochemistry of immobile elements is also useful to estimate tectonic settings of basaltic and granitic rocks (e.g.Hanson 1978, Pearce et al. 1984, Pearce 1996).However, detailed behaviour of trace elements in alkaline rocks are unknown (e.g.Wolff 1984, Bryan 2006).Specifically, geochemical studies on felsic alkaline rocks, such as nepheline syenite, alkaline syenite, phonolite, and trachyte are few.
Recent geochemical researches for the Cretaceous to Early Cenozoic felsic alkaline rocks of the state of Rio de Janeiro, Brazil, have accumulated a significant amount of major and trace element data, which are enough for the genetic discussions (e.g.Motoki et al. 2010, Sichel et al. 2012).This paper presents the behaviour of trace and rare earth elements of these alkaline rocks with additional analytical data.Based on them, the authors discuss fractional crystallization and continental crust assimilation of the nepheline syenite magmas.

SERRA DO MAR ALKALINE MAGMATIC PROVINCE
Felsic alkaline rocks, such as nepheline syenite, alkaline syenite, phonolite, and trachyte, are rare in field occurrences.In the coastal region of the state of Rio de Janeiro and São Paulo, southeastern Brazil, there are about 20 occurrences of alkaline intrusive bodies, called Serra do Mar alkaline rock province (Ulbrich and Gomes 1981).These bodies intrude into Pan-African metamorphic basement of granitic orthogneiss and pelitic paragneiss (Heilbron et al. 2000, Valladares et al. 2008), post-tectonic granite (Valeriano et al. 2011), silicified tectonic breccia (Motoki et al. 2011), and Early Cretaceous mafic dykes (Guedes et al. 2005, Motoki et al. 2009).
The Serra do Mar alkaline rock province is constituted by two magmatic alignments (Fig. 1a): A) Poços de Caldas-Cabo Frio and B) Monte de Trigo-Vitória Island.In spite of the weathering vulnerability of the constituent minerals, they have strong erosive resistance originated from mechanical strength and weathering passivity (Petrakis et al. 2010).Therefore, the intrusive bodies form massifs of 300 m to 1500 m of relative height, called alkaline massifs (Aires et al. 2012).
Poços de Caldas-Cabo Frio alignment is about 490 km long and has WNW-ESE trend (Fig. 1a, alignment A).Intrusive bodies are composed mainly of nepheline syenite and trachyte, with local occurrences of alkaline syenite on the first and lamprophyre dykes on the second.(e.g.Ulbrich 1984, Brotzu et al. 1997, 2007, Motoki et al. 2007a).They expose bottom level of flattened funnel-shaped plutons (Motoki andSichel 2006, 2008).The fission track datings for apatite indicate that the present exposures correspond to the subvolcanic intrusive bodies of about 3 km of depth from the surface of the intrusive time.The felsic alkaline bodies of the state of Rio de Janeiro belong on Poços de Caldas-Cabo Frio alignment (Fig. 1b).
The Monte de Trigo-Vitória Island alignment is about 60 km long and of WNW-ESE trend (Fig. 1a, alignment B).Different from the Poços de Caldas-Cabo Frio alignment, the intrusions are constituted mainly by alkaline syenite with eventual presence of modal quartz (Alves and Gomes 2001).

BEHAVIOUR OF THE MAJOR ELEMENTS
The geochemical data for the present paper are originated from Valença (1980), Motoki et al. (2010Motoki et al. ( , 2013) ) and Sichel et al. (2012).Additional trace element data of Tanguá intrusive complex are presented (Table I).The samples were analysed by atomic absorption (Valença 1980), X-ray fluorescence, ICP-AES (Motoki et al. 2010, Sichel et al. 2012) for major and trace elements, and ICP-MS (Motoki et al. 2010, Sichel et al. 2012;present data) for trace and rare earth elements.The geochemical behaviour of major elements of these alkaline rocks are widely different from non-alkaline rocks, such as granite and granodiorite.However, in fact, they form a negative trend (Fig. 2a) and the fractional crystallization is not a main factor of the geochemical variation.The sample Cf-1 of Cabo Frio Island complex, Rb-1 of Rio Bonito complex, and the rocks of Soarinho complex have characteristically low Na 2 O+K 2 O relative to SiO 2 .These rocks have Norm quartz, being similar to granitic rocks.The Cf-1 was collected from border of the intrusive bodies and the Rb-1 was extracted from a pyroclastic dyke.Bonito complexes are ultrapotassic and those of Cabo Frio Island complex are potassic (Fig. 2b).
The samples from the central part of the intrusive bodies, especially those of Cabo Frio Island, are peralkaline nepheline syenite (Fig. 2c).In contrast, the samples of Soarinho, Rb-1 of Rio Bonito, Tg-1 of Tanguá, and Cf-1 of Cabo Frio Island have remarkably low alkalis and high alumina.The Tg-1 is an alkaline syenite collected from the contact zone of Tanguá intrusive body.
The geochemical data range from strongly SiO 2 -undersaturated alkaline composition to oversaturated non-alkaline composition, forming a continuous trend.Such a distribution, crossing over the thermal divide, is unable to be explained by fractional crystallization.Motoki et al. (2010Motoki et al. ( , 2013) ) and Sichel et al. ( 2012) pointed out the effects of assimilation of continental crust country rocks.The relation between fractional crystallization and crustal assimilation for these alkaline rocks is expressed by the residual diagram of Norm quartz, nepheline, and kaliophilite (Fig. 3).
The relatively primitive nepheline syenite magma has high K 2 O/(Na 2 O+K 2 O) wt%, low (Na+K)/Al mol.and crystallizes leucite.This mineral is observed now as pseudoleucite in Itatiaia, Tanguá, and Morro de São João bodies.According to the fractionation of leucite, the residual liquid changes from potassic to sodic, and the magma composition approaches the cotectic curve, called Stage 1.On the cotectic curve, the residual liquid continues to become more sodic due to K-rich alkaline feldspar fractionation, called Stage 2. The residual magma becomes more SiO 2 -undersaturated, more peralkaline, and more sodic than potassic.It can arrive at the area close to the terminal point (Fig. 3, Ns).At any phase of Stage 2, crustal assimilation events can take place, transforming the magma composition from SiO 2 -undersaturated to oversaturated and form peralkaline to peraluminous, called Stage 3. The samples Rb-1, Cf-1, and those of Soarinho complex, have high grade crustal assimilation, with a ratio of about 50 wt%.

Continental crust granite and orthogneiss are SiO
2 -oversaturated and meta-aluminous, and pelitic paragneiss is SiO 2 -oversaturated and peraluminous.These country rocks have geochemical characteristics opposite to nepheline syenite in alkali-silica and alkali-alumina relations.The geochemical distribution trend of these rocks crosses over the quartz-nepheline and the aegirine- muscovite thermodynamic incompatibilities (Fig. 4).The silica saturation index (SSI) on this figure expresses the degree of its quantitative grade of felsic alkaline rocks (Motoki et al. 2010).The rocks with SSI<0 are SiO 2 -undersaturated with Norm nepheline, and those with (Na+K)/Al mol.
. >1.0 are peralkaline with Norm acmite.On the other hand, the felsic rocks with SSI>0 are SiO 2 -oversaturated with Norm quartz, and those with (Na+K)/Al mol.
. <1.0 are meta-aluminous or peraluminous with eventual occurrence of Norm corundum.The samples Tg-1, Cf-1, Cf-2, and those of Soarinho complex have SSI>0 and (Na+K)/Al mol.<1.0, and are projected on the domain and area of non-alkaline rocks.They are under strong effects of crustal assimilation and occur generally along the border zone of the intrusive complexes.The subvolcanic pyroclastic rocks, such as Rb-1, a strongly SiO 2 -oversaturated and peraluminous composition.
The degree of the fractionation of leucite and K-rich alkaline feldspar in The continuous geochemical distribution crossing over the thermodynamic incompatibilities cannot be originated from fractional crystallization.However, it can be formed by the magma superreheating (Motoki et al. 2010, 2013, Sichel et al. 2012).The SiO 2 -undersaturated hot magma can melt the wall rock generating SiO 2 -oversaturated magma.When the magma temperature is superior to the liquidus, these magmas can mix, forming thermodynamically unstable alkaline syenite.The resorption shape of clinopyroxene in the nepheline syenite and alkaline syenite, supports this model.The volatile materials extracted from the country rocks by the nepheline syenite magma heat, decrease the liquidus temperature, which makes the magma mixture of incompatible composition easier.The reaction rim of amphibole around clinopyroxene in the syenitic rocks corroborates this idea.The R 1 vs. R 2 diagram (Batchelor and Bowden 1985) shows that most of the data are projected in the area of intraplate alkaline rocks (Fig. 6b).The low R 1 values are due to the high alkalis relative to silica.Some of the samples with strong influence of continental crust assimilation, such as Rb-1, Tg-1, and those of Soarinho complex, are projected in the domain of late tectonic granite.

LARGE ION LITHOPHILE ELEMENTS
The large ion lithophile elements (LILE) are present generally in major silicate minerals, such as feldspars and micas.These elements are highly mobile during weathering and hydrothermal alteration, and are convenient for the studies of partial melting, fractional crystallization, metassomatic alteration, and fluid component behaviour.
Some of the LILEs provide important information about the magmatic evolution process.For example, Rb/K ratio indicates the grade of magma fractionation (e.g.Abbott 1966, Shaw 1968) and Rb/Ba and Rb/Sr are plagioclase fractionation indicators.Figure 7 shows variation diagrams for Rb/K mol., Rb/Ba mol., and Zr/TiO   Batchelor and Bowden (1985).
The Rb/K ratios for the felsic alkaline rocks have a good negative correlation to K 2 O/ (K 2 O+Na 2 O), with R 2 =0.461 (Fig. 7a), which is of the opposite sense to non-alkaline granitic rocks.The samples of Cabo Frio complex have high Rb/K and those of Tanguá and Rio Bonito have low Rb/K.The HFSE-based parameter Zr/TiO 2 also shows a similar characteristic, with R 2 =0.562 (Fig. 7c).The concentration of Sr and Ba, especially Ba, are widely variable, and Rb/Sr and Rb/Ba show no clear correlation to K 2 O/(K 2 O+Na 2 O) (Fig. 7b).These observations suggest that the fractionation of leucite and K-rich alkaline feldspar is a relevant and the plagioclase fractionation is not very expressive.Petrographic observations indicate that plagioclase is quite rare in these rocks (Sichel et al. 2012).
On the other hand, some LILE contents provide information of the crustal assimilation process.The Sr and Ba contents tend to increase according to the SSI.That could explain by, increase of continental crust assimilation.Their upper limits elevate linearly by SSI (Fig. 8a, b, dashed line).However, the concentrations of these elements are distributed randomly in the area beneath the upper limit,  probably due to high mobility and heterogeneous distribution of these elements in these alkaline rocks.
The Ba/Sr diagram distinguishes the trends of fractionation crystallization and continental crust assimilation.In the area of SSI<-200, this ratio increases according to the reduction of SSI due to fractional crystallization, with R 2 =0.372 (Fig. 8c, FC).On the other hand, in the area of SSI>-200, Ba/Sr increases according to elevation of SSI because of continental crust assimilation (Fig. 8c, CA), with R 2 =0.436.The critical point is SSI=-200.The samples Rb-1, Rg-1, Cf-1, Cf-2, and the rocks of Soarinho complex have high Ba, high Ba/Sr, and high SSI because of relevant influence of fluids and continental crust assimilation.
The diagrams of HFSE vs. LILE for basaltic and granitic rocks can discriminate magma source and tectonic condition (e.g.Floyd and Winchester 1975, Whale et al. 1987, Förster et al. 1997).Figure 9 shows the diagrams of Nb and Rb.The studied felsic alkaline rocks are projected in the area between volcanic arcs (Fig. 9a, VA) and within plate magmatisms (WP).Most of them are distributed along the WP line.Exceptionally, the sample Rb-1 is plotted on the VA line, which corresponds to granite of subduction zones and continental collision zones.This was probably due to strong effects of continental crust assimilation.The Nb/Y vs. Rb/K diagram confirms these observations (Fig. 9b).

RARE EARTH ELEMENTS
Rare earth element analyses are available for the selected samples of Tanguá and Rio Bonito complex (Fig. 10a, b).The total REEs is high, in average 669 ppm, which is about 10 times higher than that of granitic rocks.The average La N /Sm N and La N /Yb N (CI chondrite normalized values) are, respectively, 6.26 and 47.6, being similar to granitic rocks.Most of the samples have linear REE variation from La to Yb (Fig. 10c).The difference of the two samples is well expressed on the Gd N /Yb N vs. La N /Sm N diagram (Fig. 10d).The projected data of the analyzed rocks form a well-defined linear trend with R 2 =0.644 in direction perpendicular to the line of La N / Sm N = Gd N /Yb N (dashed line on the figure).This observation indicates that the general enrichment of LREEs is almost the same for all of the samples but the REE pattern ranges widely from concave to convex.
Concave REE pattern is generally attributed to garnet fractionation (e.g.Kay and Gast 1973, Hawkesworth et al. 1979, Terakado 1980).They are generally mafic and ultramafic rocks and show heavy REE depletion represented by Er.However, the Rb-2 has no notable heavy REE depletion.The concave pattern is represented by low Sm, Eu, and Gd.
The Eu anomalies of these samples are generally small, with average Eu/Eu* of 0.95, and therefore, plagioclase fractionation is irrelevant.The diagram of total REEs vs. Eu/Eu* shows that the Eu anomaly has linear upper limit of negative gradient (Fig. 11c).The samples Tg-2, Tg-3, and Tg-4 of Tanguá complex have high total REEs, respectively 1381, 1018, and 934 ppm, and significant negative Eu anomaly, respectively Eu/ Eu*=0.58,0.45, and 0.56.The samples of Tanguá complex present negative correlation on the diagram of SSI vs. Eu/Eu* with R 2 =0.431 (Fig. 11d), indicating that  2 O/(K 2 O+Na 2 O), wt% is 0.61, and the differentiation index is 84.26.These values suggest that this rock is originated from primitive nepheline syenite magma.In general, according to fractional crystallization total REEs increases, La N / Sm N decreases, (Sm/Sm*)* increases, and the REE pattern becomes from concave to linear.In spite of the geochemical indications of less fractionated nepheline syenite magma, the sample Rb-2 is extracted from a syenitic aplite vein.Similar phenomena are found in recent geochemical data of other syenitic aplite veins.This apparently controversial fact, is still in discussion.
In contrast, the sample, Tg-3 has high total REEs, convex pattern, and negative Eu anomaly.This sample has positive SSI of 56 and SiO 2oversaturated with strong effects of continental crust assimilation.According to the assimilation, total REEs and (Sm/Sm*)* increases, Eu/Eu* decreases, and REE pattern becomes from linear to convex.

MULTI-ELEMENTS SPIDER DIAGRAMS
The felsic alkaline rocks of state of Rio de Janeiro are under complex effects of fractional crystallization and continental crust assimilation.The multi-elements spider diagrams (Fig. 12) show peculiar patterns.General concentrations of these incompatible elements are high, being more than 10 to 100 times of average oceanic island basalt.The normalized Sr, Ba, and P are 10 to 100 times lower than the adjacent elements.The content variations of these three elements are very wide.The Ba variation could be due to highly heterogeneous fluid activity in the intrusive bodies, which is indicated by petrographic observations (Motoki et al. 2010).
Some samples show different trace element concentration patterns from the main group.The samples Rb-1 is characterized by very high incompatible elements and high Ba, showing strong influence of fluids.On the other hand, the Rb-2 is marked by low Ba.The major element behaviour indicate that the former sample is under strong effect of continental crust assimilation and the latter, originated from less fractionated nepheline syenite magma.

LOW ZR/TIO 2 OF THE PRIMITIVE NEPHELINE SYENITE
Nepheline syenite samples of Tanguá and Rio Bonito intrusive complexes have unexpectedly low Zr/TiO 2 ratio.Different from the other alkaline complexes, such as Cabo Frio Island, These rocks are projected on the areas of trachy-andesite, alkaline basalt, basanite, and nephelinite, and not of phonolite (Fig. 5a, b).There are two possible geneses for the low Zr/TiO 2 : 1) The nepheline syenite magmas of this region have originally low Zr; 2) Zircon fractionation in the magma chamber decreased Zr content.
According to the fractionation of leucite during Stage 1 and that of K-rich alkaline feldspar during Stage 2, the Zr content and Zr/TiO 2 ratio increase (Fig. 13a).The Zr/TiO 2 increases also according to elevation of peralkalinicity (Fig. 13b).Although the peralkaline syenite has higher Zr content than the meta-alkaline rocks, they contain scarce or no zircon crystals.This phenomenon is commonly observed in the other felsic alkaline complexes (Belousova et al. 2001).Only some alkaline syenite and quartz syenite samples have very small amounts of zircon, and the Zr of nepheline syenite should be present in mafic minerals.
The diagram of SSI vs. Zr/TiO 2 distinguished the trends of fractional crystallization and continental crust assimilation.The fractional crystallization (FC) from a negative correlation trend with R 2 =0.513 in the domain of SSI<-200, and the continental crust assimilation from a positive correlation trend with R 2 =0.123 in the domain of SSI>-200.The behaviour of Zr/TiO 2 on this diagram are similar to those of Ba/Sr (Fig. 13c), but the trend is clearer.In both cases, the critical SSI value is -200 (Fig. 8c, 13c).Motoki et al. (2010, 2013), and Sichel et al. (2012) proposed that the geochemical variation trend of these felsic alkaline rocks, which crosses over the two thermodynamic incompatibilities, is originated from magma super-reheating, country rock melting, and the consequent continental crust assimilation.These papers assumed complete melting of the country rock.This idea can explain the continuous variation from SiO 2 -undersaturated to oversaturated composition.It can also justify the variation from peralkaline to peraluminous composition in the case of the samples of Tanguá, Rio Bonito, and Soarinho complex, whose country rock is paragneiss with abundant modal muscovite and garnet.

PARTIAL MELTING OF THE COUNTRY ROCK
However, Itatiaia, Itaúna, and Cabo Frio Island complexes are intrusive into orthogneissic basement of granitic composition, which is not peraluminous but meta-aluminous.Therefore, some peraluminous alkaline syenite, such as Cf-1 and Cf-2 (Fig. 8), cannot be justified by the abovementioned idea.The SSI vs. (Na+K)/Al diagram indicates that all of the studied rocks, all of the felsic alkaline rocks, from a continuous trend from peralkaline to peraluminous fields regardless of the country rock type, whether it is peraluminous pelitic paragneiss or meta-peraluminous orthogneiss.The authors propose the model of country rock partial melting for the solution.
Along the contact plane, complete melting of country rock can occur because of strong thermal effects of the super-reheated magma.However, at the location of the country rock little distant from the contact between the orthogneiss and the

ACKNOWLEDGMENTS
The present research work was performed under the financial support of Brazilian Petroleum Company, PETROBRAS.The chemical analyses were performed by GEOSOL Ltd., Belo Horizonte, Brazil.Part for the fieldwork instruments, office materials, and the resources of the informatics were supported by Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro (FAPERJ).The authors are grateful to the abovementioned institutions.

Figure 1 -
Figure 1 -Locality map for the Cretaceous to Early Cenozoic felsic alkaline intrusive bodies of the state of Rio de Janeiro, modified from Sichel et al. (2012): a) Alkaline magmatic alignments of Serra do Mar province and Vitória-Trindade Chain; b) Felsic alkaline intrusive bodies of the state of Rio de Janeiro.The subaqueous pyroclastic flows of the São Mateus Volcanic Province follow that specified by Novais et al. (2007).

Figure 3 -
Figure 3 -Projection of the alkaline rocks of state of Rio de Janeiro on the residual diagram of Norm quartz, nepheline, and kaliophilite (Hamilton and MacKenzie 1960) and geochemical evolution of the nepheline syenite magma, modified from Sichel et al. (2012).

FigureFigure 5 -
Figure 6a presents tectonic setting discrimination diagrams based on trace elements.Most of the studied rocks are projected on the within-plate

2 .
The former two ratios are based on LILE and the

Figure 7 -
Figure 7 -Variation diagrams based on trace elements for the felsic alkaline rocks of the state of Rio de Janeiro: a) Rb/K mol.; b) Rb/Ba mol.; c) Zr/TiO 2 wt.The abscissa K 2 O/(K 2 O+Na 2 O) wt% represents fractionation crystallization grade of felsic minerals.
Total REE values increase either by fractional crystallization or by continental crust assimilation.The total REEs and (Sm/Sm*)* have no clear relation either to the magma fractionation indexes, or to the continental crust assimilation indexes, as SSI.

Figure 10 -
Figure 10 -Normalized rare earth elements for the felsic alkaline rocks of state of Rio de Janeiro: a) Tanguá intrusive complex; b) Rio Bonito intrusive complex; c) Unusual REE patterns; d) Gd N /Yb N vs. La N /Sm N diagram.The CI chondrite data are from McDonough and Sun (1995).

Figure 12 -
Figure 12 -Multi-elements spider diagrams for the samples of Tanguá and Rio Bonito intrusive complexes, state of Rio de Janeiro.
Schematic illustration for the contact zone of the nepheline syenite intrusive bodies and its geochemical interpretations on the diagrams of (Na+K) vs. SSI, SSI vs. Ba, SSI vs. Ba/Sr, and SSI vs. Zr/TiO 2 , which are related, respectively, toFig.4,8b,8c,and 13c.In the field of SSI<-200, negative correlation trend is observed because of fractional crystallization.In the domain of positivetrend is found due to continental crustal assimilation.This diagram distinguishes the two trends with critical SSI value of -200.4. The Sr and Ba contents have very wide variation and have no relation to fractional crystallization indexes.The SSI vs. Ba/Sr diagram distinguishes the trends of crustal assimilation and fractional crystallization in the similar way of SSI vs. Zr/TiO 2 diagram. 5. Total REEs are generally high, about 10 times that of granitic rocks.The La N and Yb N are similar to those of granitic rocks.The REE pattern is linear with small Eu anomaly because of irrelevant plagioclase fractionation.The less fractionated nepheline syenite, has low total REEs and concave REE pattern.According to the magma fractionation, total REEs and (Sm/Sm*)* increase.The samples with strong continental assimilation have high total REEs, concave REE pattern, and slight negative Eu anomaly.6.The SiO 2 -oversaturated magma generated by the country rock melting is peraluminous with high Ba, Ba/Sr, Zr/TiO 2 , and total REEs.The low-degree partial melting of country rock can generate fluid-rich peraluminous pegmatitic melt, even from the orthogneiss country rocks of meta-aluminous composition.The hyper-liquidus temperature of the magma enables the mixture between the magmas of thermodynamically incompatible compositions.