Rhönite in Cenozoic alkali basalt from Changle, Shandong Province, China, and its significance

Four Cenozoic, rhönite-bearing alkali-olivine basalt samples from the Changle area (Shandong Province, China) show an intracontinental character and were generated in an extensional setting. Petrographic studies document different generations of rhönite. In three samples, rhönite occurs either as a reaction product surrounding coarse-grained corundum, spinel and phlogopite or along cleavage planes in phlogopite. In one sample rhönite forms disseminated crystals in a mantle xenolith, possibly formed by a reaction of coarsegrained orthopyroxene or spinel with a melt. Rhönite exhibits a wide range of compositions: 22.9 wt %–33.0 wt % SiO2, 13.3 wt %–19.0 wt % Al2O3, 9.4 wt %–19.9 wt % MgO and 10.210.2 wt %–24.5 wt % FeO. The derived primary substitutions include (1) SiIV+NaVII= (Al, Fe3+)IV+CaVII, (2) MgVI= (Fe2+, Mn2+)VI and (3) TiVI + (Mg+Fe2+ +Mn2+)VI= 2Fe3+VI. Rare-earth-element (REE) patterns of euhedral rhönite crystals from the mantle xenolith (sample SS17) and those surrounding spinel (sample CL04) have a concave-upward shape for the heavy rare-earth elements (HREEs) and are slightly enriched in the light rare-earth elements (LREEs). These patterns resemble those of kaersutitic amphibole and kaersutite reported from basanite, olivine nephelinite, transitional alkali-olivine basalt and hawaiite. In contrast, REE patterns of the other two samples containing fine-grained, anhedral and acicular rhönite crystals (samples CL01 and EGS03) are relatively steep, with lower HREE and higher LREE abundances, similar to those of ocean island basalts (OIBs). All types of Changle rhönite show positive Nb, Ti and V anomalies in spidergrams normalized to primitive mantle. Mineral assemblages of the studied samples indicate that rhönite crystallized at different stages within a temperature range from about 950 to 1180 C and at pressures below 0.5 kbar, with f O2 below the NNO buffer. The chemical composition of Changle rhönite is interpreted to depend on the composition of the initial silicate melt, the redox conditions during crystallization and the composition of the minerals involved in reactions to form rhönite. Similar to metasomatic mantle amphibole, the compositions of Changle rhönites cover the I-Amph (Iamphibole) and S-Amph (S-amphibole) fields, indicating that they may have formed due to an intraplate metasomatic event, overprinting an older metasomatic subduction episode. Published by Copernicus Publications on behalf of the European mineralogical societies DMG, SEM, SIMP & SFMC. 326 F.-M. Kong et al.: Changle rhönite


Introduction
High-Ti minerals discovered in the Cenozoic alkali basalts of the Changle area (Shandong Province, China) were shown to be rhönite by Johnston and Stout (1985). Regrettably, the mineral name is spelled in different ways, such as rhoenite (Johnston and Stout, 1985;Grützner et al., 2013), rhonite (Anan'ev and Selyangin, 2011) or rhönite (Kunzmann, 1989(Kunzmann, , 1999Grapes et al., 2003;Treiman, 2008;Grapes and Keller, 2010;Sharygin et al., 2011). The official name according to the Commission on New Minerals, Nomenclature and Classification (CNMNC) of the International Mineralogical Association (IMA) is rhönite and refers to the Rhön, a low mountain range in Germany found along the border region between Bavaria, Hesse and Thuringia.
A growing number of discoveries of rhönite in alkaline effusives during the last decades show that this mineral plays an important role in the geochemical evolution of crystallizing magma or lava (Kunzmann, 1989(Kunzmann, , 1999Grapes et al., 2003;Treiman, 2008;Grapes and Keller, 2010;Sharygin et al., 2011;Grützner et al., 2013). Rhönite-bearing assemblages can be used for a rough estimate of temperature, pressure and oxygen fugacity during the crystallization of alkali basalt (Sharygin et al., 2011) and can also serve as important general indicators for primary magma evolution (Anan'ev and Selyangin, 2011). However, information on the petrographic interrelationships with other minerals is very scarce (Sharygin et al., 2011). Thus the detailed petrological and geochemical study of rhönite presented in the current work can provide important information on its genesis and the evolution of the host rock in general. Regarding the Cenozoic alkali basalt from Changle, the present study provides new insights into processes of crystallization from siliceous melts and mechanisms of secondary metasomatic formation at high temperatures.

Geological background
The Changle study area is located in Shandong Province in the western-central part of the North China Craton (NCC) and is closely related to the Tancheng-Lujiang (Tan-Lu) fault ( Fig. 1), which separates western Luxi from eastern Jiaodong. The NCC is one of the oldest cratons (3.8 Ga) on earth (Liu et al., 1992). However, it experienced lithospheric rejuvenation in the late Mesozoic and Cenozoic in that the cold, thick (> 180 km) and depleted lithospheric mantle of the eastern NCC was replaced by a hot, thin (80-60 km), younger and fertile lithospheric mantle (Fan et al., 2000;Xu et al., 2004;Zhang et al., 2009). Lithospheric thinning generated intense magmatic activity (e.g. Zheng et al., 1998;Xu et al., 2004;Chen et al., 2006;Dai et al., 2016), of which Cenozoic volcanic rocks are the most typical representatives (Zou et al., 2000;Xu et al., 2017). In addition, the Tan-Lu fault, which is an active strike-slip zone with an extension of more than 5000 km in northeastern Asia (Zheng et al., 1998), is considered to be a major channel for ascending asthenospheric melts (Zheng et al., 1998;Xu et al., 2004;Chen et al., 2006). The erupted melts, which today represent the rhönitebearing Cenozoic alkali-olivine basalts, are a result of the fault systems active in northern Asia.
Despite much work that has been done on Cenozoic alkali basalt in the eastern part of the NCC, the sources of the melts and their characteristics remain a matter of debate (Xu et al., 2012). The Cenozoic alkali basalts typically show an intracontinental background and were generated in an extensional setting (Zou et al., 2000;Xu et al., 2005;Tang et al., 2006;Liu et al., 2008a;Chen et al., 2009). They are interpreted to represent melts that are derived from the asthenosphere with "OIB-type" geochemical features, such as an en- Figure 1. Geological sketch map of the study area (modified after Kong et al., 2017), indicating the sampling localities Fangshan, Ergushan and Suishan. richment both of large-ion-lithophile (e.g. Ba, Rb, Cs, U and Sr) and high-field-strength elements (e.g. Nb and Ta) and a depletion of K (Zeng et al., 2010). The basalts are also reported to have experienced crust-mantle interaction that is reported to have occurred in an old subduction channel (Dai et al., 2016;Li et al., 2017;Xu et al., 2017). The recycled components involved were derived mostly from the Pacific Plate and stored in the transition zone (Meng et al., 2018;Xu et al., 2018;Yu et al., 2018).
On the basis of Sr, Nd, Pb, Hf and Mg isotope studies of Cenozoic basalts in eastern China, some authors have suggested that basalts with OIB-like trace-element patterns could be produced by metasomatism of large portions of the mantle wedge by interaction with partial melts derived from the Pacific slab Xu et al., 2017Xu et al., , 2018Li and Wang, 2018). Three competing petrogenetic models are currently being discussed with regard to the evolution of OIB-like basalts. One model involves mantle metasomatism, where trace-element enrichment similar to OIB was produced by infiltration of small amounts of partial melts (Pilet et al., 2008;Xu et al., 2012). A second model favours a formation of the basalt due to recycling and assimilation of oceanic crustal melts (Hofmann and White, 1982;Dasgupta et al., 2010). Finally, Niu and O'Hara (2003) envisioned a model that considers subducted oceanic crust to become a source of OIBs after a recycling process.
Cenozoic volcanism in Shandong Province occurs in the Changle-Linqu, Yishui, Penglai, Qixia and Wudi areas (Zheng et al., 1998). The Changle-Linqu volcanics studied here occur within the Tan-Lu fault zone of the Luxi part (Xiao et al., 2010); more than 50 Cenozoic volcanoes are known (Fig. 1). The effusive rocks are predominantly represented by alkali-olivine basalt, olivine nephelinite and basanite (Zeng et al., 2010). Note that the Cenozoic magmatic events in Shandong took place during two periods: 24.0-10.3 Ma and 8.7-0.3 Ma (Luo et al., 2009;Zeng et al., 2010). The older of these two periods of magmatism is characterized by the formation of large volcanoes; the major rock type is an alkali-olivine basalt, called "weakly alkaline basalt" by Zeng et al. (2010). The younger period is characterized by the formation of small and isolated volcanoes and dominated by basanite and nephelinite ("strongly alkaline basalt" as defined by Zeng et al., 2010).
The earlier magmatic period of Cenozoic basalt formation is subdivided into two stages, called the Niushan and Yaoshan stages, which have been dated ( 40 Ar / 39 Ar) at 21.0 ± 2.5 Ma and 17.3 ± 1.5 Ma, respectively (He et al., 2011a). The rock types include alkali basalt, olivine basalt and basanite. The samples of the present study belong to the Yaoshan event and originate from Fangshan, about 10 km south of Changle, Ergushan and Suishan (Fig. 1). The Ergushan and Suishan hills are close to each other and just south of the town of Beiyan (Fig. 1); these two localities are next to the Shanwang volcanic region near Linqu (Zheng et al., 1998).
As an interesting feature, megacrysts of corundum, spinel and olivine (Hu et al., 2007;He et al., 2011b;Kong et al., 2017) and mantle xenoliths of low-Mg # (the ratio of Mg / (Mg + Fe 2+ )) peridotite (Zheng et al., 1998;Xiao et al., 2010Xiao et al., , 2013 are found in these basalts; detailed petrological and geochemical studies suggest that they represent fragments of newly accreted lithospheric mantle (Xiao et al., 2010;Su et al., 2014). Examples of both megacrysts and xenoliths are found in the present suite of studied samples.

Analytical techniques
Major and trace elements of rhönite were analysed using a JXA-8800R electron microprobe at the Key Laboratory of Orogenic Belts and Crustal Evolution, School of Earth and Space Sciences, Peking University, operating at an acceleration voltage of 15 kV and a beam current of 20 nA. Natural and synthetic standards from the SPI company (Structure Probe, Inc.) of the United States used for calibration were jadeite (Si), forsterite (Mg), hematite (Fe), albite (Na, Al), rutile (Ti), rhodonite (Mn), diopside (Ca), sanidine (K), nickel silicide (Ni) and chromium oxide (Cr). An approximate error for these analyses is 0.01 wt %; representative microprobe analyses are listed in Tables 1 and 4. Trace-element analyses were performed using an Agilent 7500ce inductively coupled plasma mass spectrometer (ICP-MS), and laser ablation analyses were carried out using a pulsed 193 nm Ar excimer laser with 6 J cm −2 energy at a repetition rate of 5 Hz and with beam sizes between 32 and 60 µm, which were chosen according to the size of the measured mineral. The helium carrier gas flow velocity was 0.651 min. Data reduction was performed using the GLIT-TER 4.4.2 software (Macquarie University). Calibration was performed using USGS reference glasses (Liu et al., 2008b(Liu et al., , 2010. USGS reference material was analysed twice every five analyses. Detection limits were 0.02-0.05 ppm for rareearth elements (REE). Representative ICP-MS analyses are listed in Table 2.

Petrographical and mineralogical characteristics
Mineral abbreviations in this paper are taken from Kretz (1983). Rhönite is abbreviated as Rhoe.

Petrographical features of Changle rhönite
All rhönite-bearing samples originate from alkali-olivine basalt that belongs to the Yaoshan stage. Samples CL01 and CL04 are from Fangshan; samples EGS03 and SS17 are from Ergushan and Suishan, respectively, south of the town of Beiyan (Fig. 1). The modal amount of rhönite in each sample is about 1 %.
In the four studied rock samples, the respective type of rhönite occurrence is different . (1) In sample CL01 rhönite forms anhedral granular grains which are part of a reaction rim around coarse-grained millimetrecentimetre-sized corundum crystals (Fig. 2a). A backscattered electron (BSE) image (Fig. 2b) shows that the reaction rim consists of different zones. Corundum is surrounded by a narrow inner reaction rim of spinel that itself is enveloped by a second rim consisting of fine-grained Tibearing magnetite (Ti-Mag in Figs. 2 and 3), which forms intergrowth textures with rhönite and alkali feldspar, locally also with clinopyroxene (Fig. 2b, c). Rhönite is very fine-grained, with individual crystals not exceeding 90 µm in length.
(2) In sample CL04, euhedral-anhedral grains of rhönite occur in embayments and outer parts of millimetrecentimetre-sized spinel crystals ( Fig. 2d-f) and coexist with fine-grained phlogopite (Fig. 2e), Ti-bearing magnetite and alkali feldspar, locally also with fine-grained clinopyroxene (Fig. 2f). Fragments of mantle xenoliths associated with coarse-grained spinel do not contain rhönite, neither in the xenolith itself nor at its rim in contact with the basalt (Fig. 2d). (3) In sample EGS03, fine-grained acicular rhönite develops in a contact zone between early coarse-grained olivine and titaniferous phlogopite (Fig. 3a, b). Rhönite is also an important constituent of reaction's rim around phlogopite ( Fig. 3a, b) and in addition develops along cleavage planes and cracks in phlogopite (Fig. 3c). Rhönite coexists with Ti-bearing magnetite (Fig. 3b) and alkali feldspar (Fig. 3c), locally also with late fine-grained olivine (Fig. 3c). (4) Sample SS17 represents a distinct occurrence of rhönite because here it forms euhedral grains in broken and locally metasomatically altered parts of mantle xenoliths in basalt (Fig. 3d). It is associated with grey spinel, clinopyroxene, late fine-grained olivine, alkali feldspar and former melt; some rhönite grains are agglomerated to form a synneusis texture (Fig. 3e, f). On the basis of the features described here, all types of rhönite studied are interpreted to represent the reaction product of melt with an early primary magmatic mineral.

Structural formula and chemical composition of rhönite
The structure of rhönite was described by Walenta (1969) and Deer et al. (1978) and discussed together with that of aenigmatite, which is isostructural. Rhönite is triclinic and has the general formula X 2 Y 6 Z 6 O 20 (Deer et al., 1978;Johnston and Stout, 1985;Anan'ev and Selyangin, 2011) (Kunzmann, 1999). The classifi-          cation analysis of Kunzmann (1999) shows that the following three substitutions can be used to describe compositional variability for rhönite if rare substitutions involving B, Be and Sb are neglected (see also discussion in Johnston and Stout, 1985): where the main homovalent substituents are Mg, Fe 2+ and Mn 2+ for M 2+ and Al, Fe 3+ and Cr 3+ for M 3+ . The formulae for rhönite from Changle (Table 1) were calculated on the basis of 14 cations and 20 oxygens per formula unit (Kunzmann, 1999). The ferric/ferrous ratio was determined by charge balance.
The analyses obtained in this study (Table 1) agree with the general compositional definition of Grew et al. (2008) for rhönite within the sapphirine supergroup, and thus this name is used here. Kunzmann (1999) suggested a more detailed classification scheme for the "aenigmatite-rhönite mineral group". In terms of the Kunzmann scheme, rhönite from Changle shows Na < 1 at the X site and 0.5 < Ti < 1.5 at the Y site, conforming to his "rhönite subgroup". However, Kunzmann (1999) restricts rhönite s.s. (sensu stricto) to analyses conforming to 3 < (Si − Na) < 4, whereas the majority of Changle rhönite analyses show 2 < (Si − Na) < 3. Kunzmann (1999) notes that many "rhönite" analyses in the literature are similar and suggests that a new name should be found for such compositions in the future. Generally, the Mg # [Mg/(Mg + Fe 2+ )] of Changle rhönite is in the range of 0.5-0.84 (average of 0.65; Table 1); these values are similar to rhönite reported from enclaves in basanitic volcanics from western Anatolia in Turkey (Grützner et al., 2013), which all belong to the Mg-rich species (Sharygin et al., 2011) (Table 3). In general, rhönites with higher Mg # (0.71-0.84), as in samples SS17 and EGS03, also contain higher amounts of Na and Si than those with lower Mg # (0.5-0.58), as in samples CL01 and CL04 (Table 1; Fig. 5).

Trace-element characteristics
Rhönite typically forms a series of solid solutions (Shchipalkina et al., 2019), and Changle rhönite composition also varies (see Sect. 4) depending on the locality and the petrographical setting in thin section. Thus, apart from its variations in major-element composition, rhönite can also be expected to exhibit distinctive features related to traceelement geochemistry. In order to document these characteristics, they are displayed in chondrite-normalized REE diagrams (Figs. 7, 8) and primitive-mantle normalized diagrams for incompatible elements (Fig. 9). Not much information is available from the literature, but Kogarko et al. (2005) presented trace-element data for rhönite from tephrite-phonolite on the islands of Bravo and Fogu, Cabo Verde archipelago (referred to as BFCV-Rhoes in this paper), which are given for comparison in Figs. 7 and 8. These islands are of volcanic origin and located along the African margin of the central Atlantic Ocean. They are interpreted as resulting from hotspot volcanism (Ramalho et al., 2010). Coltorti et al. (2010) suggested that the lavas are contaminated by old fragments of sub-continental lithospheric mantle. In Fig. 7 the REE patterns of OIBs (Sun and McDonough, 1989) are given for comparison as well.

REE characteristics
The REE patterns of rhönite from the three different localities of the Changle area are presented in Fig. 7. The REE patterns of rhönite from samples CL01 and EGS03 are distinctly different from those of samples CL04 and SS17 in that they are relatively steep (with (La/Yb) N = 5.26-39.11; Table 2; Fig. 7). It is obvious that two of the Fangshan rhönite patterns (CL01-Rhoe 17 and CL01-Rhoe 19), show a pronounced increase towards the light rare-earth elements (LREEs) (Fig. 7a). Towards the heavy rare-earth elements (HREEs) a general decrease from Dy to Tm is observed, followed again by an increase towards Lu (CL01-Rhoe 17) (Fig. 7a). In general, the LREE content in rhönite from sample CL01 (Fig. 7a) is significantly higher than that from sample CL04 (Fig. 7b) from the same locality. The patterns of the two Ergushan rhönites (EGS03; Fig. 7c) show a pronounced enrichment in LREEs and a decrease towards HREEs. Regarding the LREEs, the amounts of La, Ce, Pr, Nd, Sm and Tb in rhönite of sample EGS03-Rhoe 32 are significantly  (Sun and McDonough, 1989) and compared to literature data. BFCV-Rhoes is rhönite from the islands of Bravo and Fogu, Capo Verde archipelago, from Kogarko et al. (2005); OIB is ocean island basalt from Sun and McDonough (1989). higher than those in EGS03-Rhoe 25 (Fig. 7c). On the HREE side, the patterns are similar and marked by a zigzag pattern.
There is essentially no difference in REE chemistry between rhönite within the xenolith of sample SS17 and euhedral rhönite closely associated with spinel megacrysts of sample CL04. The respective REE patterns are slightly enriched in LREEs and show a decrease toward the HREEs, with (Ce/Sm) N = 0.55-1.05 (Fig. 7b, d; Table 2). The rhönite patterns do not show any significant anomalies. Most of the rhönites from Fangshan (CL04) show a convex-upward shape, and, compared to rhönite from Suishan (SS17), the amount of LREEs is somewhat higher (Fig. 7b, d). Furthermore, the patterns of rhönite from Suishan (SS17) are generally relatively flat, with a slight convex-upward shape with a maximum at Sm (Fig. 7d). All 10 measured rhönite patterns from Suishan are similar to each other; only the amount of Lu is distinctively higher in samples SS17-Rhoe 10 and SS17-Rhoe 12 and lower in SS17-Rhoe 9 (Fig. 7d).
The spectrum of REE patterns of Changle rhönite from samples SS17 and CL04 is most similar to BFCV-Rhoes of Kogarko et al. (2005), which is phenocrysts in tephritephonolite (Fig. 7). In general, the REEs of the Changle rhönite are more depleted than those studied by Kogarko et al. (2005;see Fig. 7). On the other hand, rhönite in samples EGS03 and CL01 shows a spectrum of patterns which is similar to that of OIB (Sun and McDonough, 1989) (Fig. 7a, c); however, the EGS03 and CL01 patterns are depleted compared to OIB. Figure 8 shows the distinguishing features of all rhönite REE patterns from this study. It is obvious that the rhönite patterns from samples CL04 and SS17 are quite similar, except for higher amounts of LREEs in CL04. Figure 8 also shows that the spectrum of rhönite REE patterns from samples CL04 and SS17 is similar to that of kaersutitic amphibole megacryst in various types of alkalic basalt such as basanite, olivine nephelinite, transitional alkali-olivine basalt and hawaiite (Irving and Frey, 1984), and kaersutite from Cenozoic alkali basalt from the South China Sea (Yan et al., 2015), which, however, are generally more enriched in REE. Rhönite patterns from samples CL01 and EGS03 are enriched in LREEs and generally follow the trend of wholerock REE patterns of Changle weakly alkali basalt (Zeng et al., 2010), which are old basalts 24.0-10.3 Ma in age (Yu et al., 2010) and represent an important period of volcanism on Shandong peninsula.

Incompatible-element characteristics
In general, the primitive-mantle normalized diagrams for incompatible-element contents of all Changle rhönite samples indicate very similar trends; Nb, Ti and V typically show a pronounced positive anomaly. In addition, the pat-  (Sun and McDonough, 1989) and compared to kaersutitic amphibole (Irving and Frey, 1984) and to weakly alkali basalt (Zeng et al., 2010). "Krs from SCS" is kaersutite from South China Sea, after Yan et al. (2015).
terns of samples CL01 (Fig. 9a) and CL04 (Fig. 9b) show a small positive Zr anomaly. A small positive Sr anomaly is observed only in rhönite from sample CL01 (Fig. 9a). As already observed in their REE patterns (Fig. 7c), the two samples from Ergushan differ significantly from each other, especially with respect to the LREEs, and similarly the incompatible elements (Fig. 9c) also demonstrate considerable differences. Figure 9c shows that sample EGS03-Rhoe 32 has substantially higher amounts of Ba, La, Ce and Pr compared to EGS03-Rhoe 25, in addition to Nd and Sm anomalies.
The spectrum of rhönite patterns from all three localities generally conforms to that of the rhönite sample from Kogarko et al. (2005), the host rock which is interpreted to have been derived from rocks of a hotspot-intraplate, oceanicisland environment. Although the patterns are similar and in part nearly identical, the Changle rhönites (with the exception of Ba, V and Cr) tend to be poorer in incompatible elements than those described by Kogarko et al. (2005) (Fig. 9ad). Note particularly that the rhönite patterns of sample CL04 (Fig. 9b) nearly mimic those from Kogarko et al. (2005), and even the positive Zr anomaly is present. The rhönite from sample CL01, in general agreement with sample EGS03, is somewhat enriched in highly incompatible elements. The rhönite pattern from sample SS17 is again quite similar to those described by Kogarko et al. (2005); nevertheless, the positive Zr anomaly is missing (Fig. 9d). Perhaps the rhönite from Kogarko et al. (2005) is more enriched in incompatible elements because the hosting rock is an alkali tephritephonolite dyke and not of basaltic composition, which then would result in differences in REE and trace elements in general.
In addition, it is noteworthy that in the spidergrams of Fig. 9 the incompatible elements of all studied rhönites demonstrate pronounced positive anomalies of Ti, Nb and V (Fig. 9). Compared to associated minerals (such as spinel, clinopyroxene), rhönite has distinctly more Ti, V and Nb except for the last-crystallized, fine-grained irregular ilmenite (data from unpublished author).

The formation of Changle rhönite in view of other occurrences
The petrographical setting of rhönite, especially with respect to morphology and location in different parts of the magmatic fabric, gives important indications for the order of crystallization. Rhönite grains in samples SS17 and CL04 are euhedral and larger compared to rhönite from sample CL01 and the needle-like anhedral grains surrounding phlogopite or within the cleavage planes of phlogopite in EGS03. These petrographical features indicate that the larger, euhedral rhönite in SS17 and CL04 formed earlier than the smaller, anhedral or acicular grains in EGS03 and CL01. Despite having formed at different stages, all rhönite studied is interpreted to represent the product of a reaction of melt with early or primary minerals (Figs. 2, 3). Sample SS17 is a locally altered mantle xenolith, and rhönite possibly formed by reaction of a melt with coarse-grained orthopyroxene or spinel (Kong et al., 2013 (Fig. 3e-f). The petrographical features and mineral relationships suggest that the following overall rhönite-forming reactions can be constructed. As further confirmation of these suggestions, the chemical analyses of the minerals reacting with rhönite are given in Table 4.
The reaction for the formation of anhedral rhönite as part of a narrow reaction zone surrounding spinel in sample CL01 (Fig. 2a-c) and of euhedral to anhedral rhönite in sample CL04 in embayments surrounding spinel (Fig. 2d-f) is The reaction for the formation of fine-grained acicular rhönite in sample EGS03 that developed between an early olivine megacryst and phlogopite, or around the phlogopite or along cleavage planes and cracks in phlogopite (Fig. 3ac), is (2) phlogopite/olivine megacryst + melt → rhönite + Tibearing magnetite ± late fine-grained olivine + alkali feldspar.
Although the Changle rhönites show considerable variations in major element compositions (Table 1; Figs. 4-6), some of them are similar to specific compositions of rhönite from other localities, for instance to those from tephritephonolite on the islands of Bravo and Fogu, Cabo Verde archipelago (BFCV-Rhoes; Kogarko et al., 2005); to rhönite in olivine phenocrysts from alkali-basaltic fields of Russia, Hungary and Israel (Sharygin et al., 2011); and to rhönite from basanitic volcanic rocks of the Quaternary Kula volcanic province in western Anatolia, Turkey (Grützner et al., 2013) (Table 3). Note that rhönite from samples SS17 and EGS03 contains high amounts of MgO and Cr 2 O 3 and low amounts of FeO and MnO, whereas rhönite from samples CL04 and CL01 shows the opposite trend (Fig. 4). Correspondingly, Grützner et al. (2013) also observed different compositions of rhönite in the same rock. These authors showed that rhönite in enclaves from basanite in Turkey has higher MgO, Al 2 O 3 and Cr 2 O 3 but lower FeO and MnO compared to the rhönite in the matrix of the basanite. Compared to the rhönite in the enclaves of Grützner et al. (2013), those in SS17 and EGS03 from Changle have higher MgO and Cr 2 O 3 but lower FeO and MnO, whereas the rhönite in CL01 and CL04 is similar to the matrix type of Grützner et al. (2013), which, however, is slightly higher in FeO, TiO 2 and Cr 2 O 3 and lower in SiO 2 , MgO and MnO. Thus, the chemical composition of Changle rhönite appears to depend to a large extent on the composition of the initial silicate melt. Sharygin et al. (2011) confirmed that variations in rhönite composition can be correlated with the composition of the initial silicate melt and with the redox conditions during crystallization. Rhönite in EGS03 has a composition of higher Na and Si and lower Ca and Al and is similar to rhönite with NaSi-rich and CaAl-poor compositions in lower-crustal pyroxenite xenoliths within phonolite from the Mount Sidley composite volcano, Marie Byrd Land, Antarctica (Grapes et al., 2003). The latter are reported to indicate a formation at high oxidation conditions at a late growth stage (Grapes et al., 2003). Thus, the chemical composition of Changle rhönite also depends on the redox conditions during crystallization.
Rhönite is also reported to replace Ti-bearing amphibole or kaersutite (Wilshire et al., 1971;Johnston and Stout, 1985;Huckenholz et al., 1988;Kunzmann, 1999;Nédli and Tóth, 2003;Grapes et al., 2003). Note that high amounts of Na + Si IV and low amounts of Ca + Al IV are typical features of rhönite that is assumed to be a secondary phase produced by alteration of kaersutite, whereas primary rhönite is reported to contain low amounts of Na + Si IV and high Ca + Al IV (Nédli and Tóth, 2003;Grapes et al., 2003;Doroozi et al., 2016). In this respect, Changle rhönite from samples SS17, CL01 and CL04 conforms to the chemical characteristics of a primary phase and is low in Na + Si IV and high in Ca + Al IV (Na + Si IV : 3.38-3.62; Ca + Al IV : 4.29-4.59). From the perspective of petrographical observation, Changle rhönites are considered to be products of reactions of melt with spinel, phlogopite, olivine or orthopyroxene (Figs. 2 and 3). Only rhönite from the sample EGS03 shows chemical characteristics that are similar to those reported for secondary rhönite, replacing Ti-bearing amphibole or kaersutite (Fig. 3b, c; Table 1) (Nédli and Tóth, 2003;Grapes et al., 2003;Doroozi et al., 2016), although no Ti-bearing amphibole or kaersutite has been observed in this sample.
The most comprehensive experimental data on rhönite stability relationships have been provided by Kunzmann (1989; see also Kunzmann, 1999, for a summary). On its own composition, rhönite is stable at 1 bar from 850-1000 • C to at least 5 kbar at 900 to 1100 • C (Kunzmann, 1989) at quite variable conditions of oxygen fugacity. For instance, Kunzmann (1989Kunzmann ( , 1999 reports experiments at 1 bar in which rhönite was synthesized under conditions of the HM (hematite-magnetite), NNO (Ni-NiO), QFM (quartzfayalite-magnetite) and IQF (iron-quartz-fayalite) oxygen buffers. However, for bulk compositions reflecting basaltic systems, the P -T stability field of rhönite is much more restricted. Kunzmann (1989) documented the overlapping P -T stability fields of Ti-rich calcium amphibole and rhönite for an alkali-basaltic system under NNO oxygen buffer conditions in the liquidus-solidus P -T region (Fig. 10). Considering that Changle rhönite does not coexist with a Ti-rich amphibole and that Changle rhönite is part of an assemblage that includes Ti-bearing magnetite, feldspar, clinopyroxene and melt, a P -T region of potential formation can be mapped out (Fig. 10). On the basis of Kunzmann's (1989) results, Changle rhönite could form in the temperature range from 950 to 1180 • C at maximum pressures of 0.5 kbar, where the rhönite-out reaction curve is located (Fig. 10). Boivin (1980) studied rhönite stability at 1 bar as a function of temperature and oxygen fugacity (IW is iron-wüstite; NNO) for bulk compositions equivalent to natural kaersutite, basanitoid (a rock intermediate between olivine basalt and basanite) and hawaiite. The results indicate that rhönite sta- Figure 10. Liquidus-solidus relationships of Ti-rich calcium amphibole and rhönite in alkali basalt up to 1 kbar at NNO oxygen buffer conditions after Kunzmann (1989); the stability field of Changle rhönite is shaded in light grey. bility is not very sensitive to oxygen fugacity, as noted above, but depends strongly on bulk-rock composition. Going from the kaersutite to basanitoid and hawaiite bulk compositions, the contents of TiO 2 and the modal foids nepheline and leucite decrease regularly (Boivin, 1980). At 1 bar, rhönite is stable up to ca. 1150 • C (IW buffer) or 1170 • C (NNO buffer) for kaersutite to ca. 1090 • C (both buffers) for basanitoid composition. No rhönite was observed for the hawaiite composition. Thus undersaturated alkali-basaltic melts appear to be a prerequisite for rhönite crystallization, and higher TiO 2 contents favour stability of rhönite to higher temperatures.
For a lack of suitable geobarometers, most studies involving rhönite have accepted the 500-600 bar upper pressure limit for rhönite experimentally derived by Kunzmann (1989). For the Cabo Verde rhönite occurrence, Kogarko et al. (2005) also used the experimental data of Boivin (1980) to conclude that rhönite + clinopyroxene began to crystallize at 1100-1150 • C. However, temperatures for rhönite stability have independently also been estimated from homogenization and crystallization studies on rhönite-bearing melt inclusions in olivine or liquidustemperature determinations on the enclosing basaltic rock. Grapes et al. (2003) concluded that the assemblage rhönite + olivine + clinopyroxene + plagioclase in Mount Sidley xenoliths (Marie Byrd Land, West Antarctica) formed within a narrow temperature range of 1140-1190 • C. Grapes and Keller (2010) inferred that Fe 2+ -dominant rhönite (Mg # < 0.5; see Table 1) coexisting with subsilicic Al-Tibearing augite, olivine and titanian magnetite in tephrite and limburgite from the Kaiserstuhl volcanic complex, SW Germany, formed at temperatures of 1068-1137 • C. Sharygin et al. (2011) suggested that Mg-dominant rhönite in olivinehosted, silicate-melt inclusions in alkali basalt crystallized in the range of 1180 to 1260 • C. Some of the above estimates exceed the temperatures suggested by the work of Kunzmann (1989) in Fig. 10, and the TiO 2 content of rhönite may play a role. In their study, Sharygin et al. (2011) cite quite variable TiO 2 contents, between 8.1 wt % and 12.4 wt % TiO 2 , without correlating their relatively high estimated temperatures directly with TiO 2 contents. Kóthay et al. (2003) indicated in an abstract that TiO 2 -rich rhönite in alkali basalt from the Bakony-Balaton Highland and the Nógrád-Gömör volcanic fields, Pannonian Basin, Hungary, can be stable up to temperatures as high as 1220-1270 • C. Although Kóthay et al. (2003) did not give any explicit information on TiO 2 contents, Sharygin et al. (2011) reported that rhönite from the same locality contains TiO 2 in the range of 9.0 wt %-11.63 wt %. TiO 2 contents in Changle rhönite range from 10.5 wt % to 12.7 wt % (Table 1) and are therefore as high or even higher than those reported by Sharygin et al. (2011). They are also similar to and in part higher than in rhönite (10.7 wt %-11.8 wt %) from the islands of Bravo and Fogu, Cabo Verde archipelago (Kogarko et al., 2005), and significantly higher than in the Kaiserstuhl volcanic complex (8 wt %-10 wt %; Grapes and Keller, 2010). In particular, rhönite in sample SS17 and some in EGS03 contains very high amounts of TiO 2 up to 12.5 wt %.

Comparison of REE patterns between rhönite and kaersutitic amphibole
The composition of the rhönite-bearing basalt of the current study is similar to the weakly alkali basalt of Zeng et al. (2010) that belongs to the generation of old Changle basalts (24.0-10.3 Ma) and that is derived from asthenospheric mantle melts. Thus, in addition to the REE patterns of Changle rhönite from this study, the range of weakly alkali basalt compositions from Zeng et al. (2010) is shown in Fig. 8 for comparison as well. Additional patterns presented in Fig. 8 are from kaersutitic amphibole megacrysts in various types of alkalic basalt such as basanite, olivine nephelinite, transitional alkali-olivine basalt and hawaiite (Irving and Frey, 1984) and from kaersutite of Cenozoic basalt from the South China Sea (Yan et al., 2015). Note that some of the REE patterns of Changle rhönite are very similar to the patterns of mantle amphiboles such as kaersutitic amphibole megacrysts or kaersutite in basalt (Irving and Frey, 1984;Yan et al., 2015;Fig. 8). Especially the REE patterns of euhedral rhönite in sample CL04 which show a maximum at Nd (Fig. 7b) are similar to kaersutitic amphibole megacrysts in basalts, in which only the relative amounts of all the REEs are moderately higher (Irving and Frey, 1984). These results imply that rhönite exhibits a geochemical fractionation behaviour of certain trace or REE elements similar to that of kaersutite or kaersutitic amphibole.
Overall, the REEs of rhönite in samples SS17 and CL04 are more depleted than those of kaersutite from the South China Sea (Yan et al., 2015;Fig. 8 -SCS) and kaersutitic amphibole megacrysts (Irving and Frey, 1984;Fig. 8).
Rhönite from samples EGS03 and CL01 is slightly depleted in MREEs (middle rare-earth elements) compared to kaersutite from the South China Sea and kaersutitic amphibole megacrysts (Yan et al., 2015;Irving and Frey, 1984;Fig. 8). The shape of the REE patterns of rhönite from samples CL01 and EGS03 is similar to that of OIB (Sun and McDonough, 1989;Fig. 7a, c) and the Changle weakly alkali basalt of Zeng et al. (2010) (Fig. 8), although the amounts are slightly lower and rhönite from these two samples shows a marked enrichment of Lu. This phenomenon could also be the result of increasing alkalinity of the melt during the process of magma differentiation, which would lead to an enrichment of light REEs and some other elements such as Lu and Nb (Table 2). Thus not only the petrographical but also the REE features indicate that rhönite in samples EGS03 and CL01 crystallized as a late phase and perhaps contemporaneously with increasing alkalinity of the host magma.

Chemical fingerprints of Changle rhönite with reference to amphibole
By studying suprasubduction and intraplate metasomatic signatures, Coltorti et al. (2007b) concluded that mantle amphiboles are an important indicator that can mirror diverse geochemical characteristics; they thus contribute significantly to a better understanding of the respective geotectonic setting. Coltorti et al. (2007b) distinguished two types of amphibole from mantle xenoliths in basalt, an intraplate type (Iamphibole) and a suprasubduction type (S-amphibole). Geochemical signatures of I-amphibole were shown to be clearly different from those of S-amphibole due to their higher Nb, Ta, Zr, Hf and Ti contents (Coltorti et al., 2007b). By comparing the chemical signatures of rhönite with those of amphibole, it appears possible to derive important petrological and geochemical fingerprints typical for a particular setting or particular melts that originate from the subduction of oceanic crust.
Recently, Dai et al. (2016), Li et al. (2017) and  suggested that the sources of the OIB-like Cenozoic alkali basalts in the North China Craton (such as those from Shandong Peninsula) had experienced a "metasomatizing reaction" between the asthenospheric peridotite with melts derived from deeply subducted oceanic crust. The subducting slab stagnating in the mantle transition zone became a source of melts that carbonated-metasomatized the overlying mantle wedge peridotite (Li and Wang, 2018;. Rhönite-bearing Cenozoic alkali basalts from the Changle area are reported to represent typical intracontinental basalts (Zou et al., 2000;Xu et al., 2005;Tang et al., 2006;Liu et al., 2008a;Chen et al., 2009), and they are tectonically located above the old subduction zone, where the westward subduction of the Pacific slab occurred Li and Wang, 2018). In essence, the Changle rhönitebearing basalts are derived from melts originating in an intracontinental setting where a reaction of the asthenospheric Figure 11. Chondrite-normalized (Sun and McDonough, 1989) incompatible-element spidergram for Changle rhönite. Depleted disseminated amphibole and enriched disseminated amphibole are from the basalt of Kapfenstein, Austria (Coltorti et al., 2007a), and I-Amph (I-amphibole) is from xenoliths from Antarctica (Coltorti et al., 2004(Coltorti et al., , 2007a peridotite with melts that derived from subducted oceanic crust occurred. As discussed above, REE patterns of Changle rhönite are similar to those of kaersutite in Cenozoic basalt (Yan et al., 2015) and kaersutitic amphibole megacrysts in basalt (Irving and Frey, 1984;Fig. 8). Although Changle rhönite does not represent a direct breakdown product of a former amphibole, the geochemical similarity between rhönite and amphibole suggests that rhönite could be able to take over the role of the petrogenetic indicator mineral amphibole by incorporating appropriate typical elements in a similar way. Indeed, the REE patterns of most rhönites resemble those of kaersutitic amphiboles or kaersutites in basalt (Yan et al., 2015;Irving and Frey, 1984;Fig. 8) and can signal suprasubduction and intraplate metasomatic mantle signatures (Coltorti et al., 2004(Coltorti et al., , 2007b. The metasomatic mantle amphibole chosen for comparison originates from xenoliths of peridotite and lherzolite in basalt from Austria and Antarctica (Coltorti et al., 2004(Coltorti et al., , 2007a. Note that the geochemical data of disseminated amphibole in lherzolite xenoliths in basalt from Kapfenstein, Austria, show subduction-related and intraplate metasomatic features (Coltorti et al., 2007a). Note also that LREEs and some incompatible elements of these amphiboles are similar to those of Changle rhönite from sample SS17, which also forms disseminated crystals in xenoliths (Fig. 11).
In a chondrite-normalized spidergram (Fig. 11), the Changle rhönite shows complex geochemical characteristics, revealing a distinct negative K anomaly and positive anomalies of Ta and Ti. Incompatible elements of Changle rhönite show (Fig. 11) that the behaviour of Rb, Ba, Ta and most REE (except Ho and Lu) is similar to that in I-amphibole from xenoliths from Antarctica (Coltorti et al., 2004(Coltorti et al., , 2007b. However, most incompatible elements of Changle rhönite (apart from Ti) are more strongly depleted than those of I-amphibole (Fig. 11). Depleted disseminated amphibole from Kapfenstein yields a similar trend in Th, U, Eu, Ti, Gd and some HREEs (Coltorti et al., 2007b) to most Changle rhönites, which are, however, slightly depleted in these elements (Fig. 11). On the other hand, most Changle rhönites show a trend with respect to Sr and Ho to Lu that is similar to the enriched amphibole from the Kapfenstein locality (Fig. 11).
Importantly, Coltorti et al. (2007b) have shown that Nb is the best discriminator element for distinguishing suprasubduction zone and intraplate settings. I-amphibole domains have relatively high amounts of Nb within a rather restricted range of about 10-1000 ppm compared to S-amphibole domains, where Nb is generally (with very little overlap) below 10 ppm. Assuming that the geochemical behaviour of these elements in amphibole is similar to that in rhönite, we see that although the four studied samples from Changle are located close to each other and also belong to the same geological formation (Yaoshan), they show clear distinguishing features. Whereas rhönite from sample SS17 contains 1.75-5.31 ppm Nb and thus indicates an S-type "amphibole" character, rhönites from the other samples, CL01, CL04 and EGS03, contain distinctively higher amounts of Nb (11.96-29.18 ppm) and point to an I-type character (Fig. 12a). Note that there is also a clear difference in Nb content between samples CL01, CL04 and EGS03; they form distinctive groups.
According to McDonough and Sun (1995), a chondritic character is indicated by Ti/Nb ch = 1823 and Zr/Nb ch = 15.7, where Ti/Nb > 1823 and Zr/Nb > 15.7, Ti/Nb < 1823 and Zr/Nb < 15.7, respectively, indicate a suprachondritic or subchondritic character (Coltorti et al., 2007b). In addition, Coltorti et al. (2007b) showed that Ti/Nb and Zr/Nb ratios of amphibole can be used to distinguish between S-amphibole, which is chondritic to slightly suprachondritic, and I-amphibole, which yields a subchondritic character (Fig. 12b). With respect to rhönite from Changle, Ti/Nb ratios of rhönite from samples SS17 and CL04 are much higher than 1823 and would thus indicate a typical suprachondritic character, whereas rhönite from samples EGS03 and CL01 would indicate a subchondritic character ( Fig. 12b; Table 2). However, with regard to Zr/Nb ratios, most of the Changle rhönites are subchondritic, and only one rhönite crystal in the xenolith from sample SS17 yields higher values than 15.5 and is thus suprachondritic (McDonough and Sun, 1995) (range of 1.65-17.63; Table 2).
In the Zr-Nb diagram (Fig. 12a), all rhönites from sample SS17 are located in the S-amphibole field, whereas rhönite from samples CL01 and EGS03 (with one outlier) plots in the I-amphibole field; all rhönites from sample Cl04 plot in between. In the diagram of Zr/Nb versus Ti/Nb (Fig. 12b), rhönite from samples SS17 and CL04 falls into or close to the S-amphibole field, whereas rhönite from samples CL01 and EGS03 is located outside but close to the I-amphibole field. According to Coltorti et al. (2007b), amphiboles straddling the I-and S-amphibole fields possibly indicate a suprasubduction to intraplate transition stage. Assuming that such data on amphibole are transferable to rhönite, those which document an overlap of I and S character may have experienced an intraplate metasomatic event overprinted by a subduction-related metasomatic episode. Indeed, with respect to the local geological environment of the rhönitebearing Changle intraplate alkali basalt, there is petrological and geochemical evidence of an old Pacific subduction event and further Cenozoic intraplate metasomatism.
Note that Ionov and Hofmann (1995) hypothesized that the amphibole-melt partition coefficient for Nb is lower than that for Ta during mantle metasomatism, based on their study of disseminated and vein-related amphibole separated from mantle xenoliths in alkali basalt (Ionov and Hofmann, 1995). As a consequence of progressively increasing Nb/Ta ratios in the residual melt resulting from the precipitation of amphibole, earlier-formed, vein-related amphibole has a lower Nb/Ta ratio than disseminated amphibole of a later generation (Ionov and Hofmann, 1995). As a characteristic feature of the Changle rhönite, the Nb/Ta ratio of euhedral crystals from sample CL04 is in general quite low (7.2-10.3), whereas the range of ratios of the anhedral crystals from the acicular type CL01 is distinguishably higher and between 16.3 and 17.3 (Table 2; Fig. 2b, c, e, f). The Nb/Ta ratio of euhedral grains from sample SS17 varies widely between 3.3 and 15.9 and reaches the highest values in fine-grained rhönite from sample EGS03 (19.6-23.6) ( Table 2; Fig. 3b, c, e, f). On average, these differences can help to indicate possible growth sequences. In conjunction with the petrographical features, we propose that the larger euhedral rhönite crystals formed earlier than the anhedral smaller ones.

Conclusions
Mineralogical, petrographical and geochemical characteristics of rhönite from three alkali-olivine basalts and one xenolith in basalt from Changle show that rhönite grew at differ-ent stages due to different reactions and formation conditions but generally formed at temperatures exceeding 950 • C and pressures below 0.5 kbar.
Euhedral rhönite from samples SS17 and CL04 is proposed to have crystallized earlier than the acicular and anhedral grains in samples EGS03 and CL01.
Rhönite is documented to have a wide range of compositions in major oxides; the main substitutions derived are Besides providing important local insights into the formation mechanisms of rhönite from the Changle alkali-basalt occurrences, the important objectives of the current paper are to present REE and trace-element data for different rhönite generations as well as a compilation of crystal-chemistry characteristics which may serve as fundamental information on rhönite-bearing rocks in general. Since rhönite and amphibole have structural similarities and also share geochemical characteristics, we anticipate that with future studies of the REE characteristics of amphibole and rhönite from other localities, a more precise "geochemical link" between these minerals can be established, and rhönite may become an important indicator for discriminating different geotectonic settings.
Data availability. All data derived from this research are presented in the enclosed text, tables and figures.
Author contributions. L-QZ and X-HL carried out trace elements of rhönite, F-MK carried out major elements of rhönite and wrote the first draft, H-PS and X-PL contributed to the experimental strategy, interpretations, discussion and manuscript revision.