Petrogenetic history and melt inclusion characteristics of mantle plume-derived ijolites from NE India: implications for multistage crystallization and occurrence of “nano-calciocarbonatites”

The Sung Valley ultramafic–alkaline–carbonatite complex (UACC) of Meghalaya, NE, India, is a result of magmatic activity related to the Kerguelen mantle plume spanning from 101 to 115 Ma. In the present study, an integrated crystal size distribution (CSD), mineral chemistry, and melt inclusion analysis are carried out on the ijolites present within this UACC. The CSD analysis shows that these ijolites were formed in multiple stages through changes in the crystallization environment, such as cooling and nucleation rates. Raman spectroscopy of mineral inclusions of rutile, aphthitalite, apatite, carbonate–silicate melt inclusions, and disordered graphite within clinopyroxene and titanite, respectively, indicates a heterogeneous composition of the parental magma. These mineral and melt inclusion phases further suggest localized changes in oxygen fugacity ( f O 2 ) due to redox reactions in the lower crust. SEM–EDX analysis of the exposed melt inclusions reveals the presence of alkali-bearing diopside, phlogopite, and andradite, along with an unidentified carbonated silicate daughter phase. The studied melt inclusions are dominated by carbonate, whereas silicates are subordinate. The presence of this fully crystallized carbonate–silicate melt as calcite, diopside, phlogopite, magnetite, apatite, and andradite suggests the presence of “nano-calciocarbonatites” in these ijolites. Our study provides insights into different mechanisms of the loss of alkalies from initially entrapped alkaline carbonate melt in clinopyroxenes. The predominant occurrence of calcite as the only carbonate phase in the studied melt inclusions is a result of silicate–carbonate melt immiscibility, calcite-normative system in these inclusions, dealkalization of the alkaline carbonates in the presence of external fluid, and/or redistribution of the alkalies to the daughter alkali-bearing silicates.

A complete petrogenetic history of igneous rocks also involves the physical history of crystallization, which is recorded by the textures of igneous rocks and explained through crystal size distribution (CSD) analysis. Quantitative textural measurements such as CSD allow direct assessment of the crystallization history of the minerals constituting the igneous rocks. Processes related to crystallizing magmas, such as magma rheology, fractionation, mixing, and cooling, can be very well explained by the CSD analysis (Jaeger 1968;Marsh 1988;Cashman and Marsh 1988;Armienti et al. 1994;Higgins 2000;Klein et al. 2018). CSD gives deep insights into the crystallization processes through the observed textural changes governed by the magma chamber's mechanical processes (Higgins 2011). The slope of the CSD curve helps to understand different processes during the crystallization of magma batches (Marsh 1988). A quantitative measure of the gain or loss of the crystals over a particular size range can also be explained through these CSD curves, which may be useful in evaluating the importance of physical processes such as crystal accumulation and fractionation in a petrologic system (Cashman and Marsh 1988;Marsh 1988).
The physical and chemical evolution of the crystallizing minerals during magmatic processes can also be understood with the help of mineral chemical studies combined with quantitative measurements such as CSD. Therefore, the chemical composition of minerals and trapped melt inclusions are also equally important in order to construct a thorough petrogenetic record of igneous rocks. Melt inclusions in the crystals are the pockets of parental melt trapped during the crystal growth (Faure and Tissandier 2014). These inclusions not only represent the characteristics of the magma from which their host crystals grew (Anderson 1979;Roedder 1979;Sobolev 1996;Frezzotti 2001;Danyushevsky et al. 2002;Guzmics et al. 2008Guzmics et al. , 2011Mitchell 2009) but also record evolving melt compositions and physicochemical conditions prevailing during crystallization (Roedder 1972). Carbonatesilicate melt inclusions serve as a powerful tool to assess the paragenetic relationship among silicate and carbonate rocks of alkaline-carbonatite complexes (Andreeva et al. 2007). These carbonate melts also work as effective agents for the transportation of rare earth and alkaline earth elements (Guzmics et al. 2009;Mitchell 2009). The alkali-rich composition of carbonate-silicate melt inclusions in alkaline-carbonatite rocks places them very close to the representation of the actual parental magma composition that these rocks crystallize from (Nielsen et al.1997;Sokolov et al. 1999;Yaxley et al. 2022). The investigation of these melt inclusions of alkaline rocks (ijolites) becomes even more crucial while assessing the petrogenetic relationship among alkaline-carbonatite rocks, such as liquid immiscibility during crystallization of their parental magmas (Guzmics et al. 2012;Chayka et al. 2021;Berkesi et al. 2020).
In the northeastern part of India, ijolites are exposed in the Sung Valley UACC, along with ultramafic rocks and carbonatites. Previous geochemical and stable and radiogenic isotopic studies suggest that these different suites of rocks from Sung Valley UACC are not co-genetic (Ray and Pande 2001;Srivastava and Sinha 2004;Srivastava et al. 2005). In the present work, we integrated crystal size distribution, mineral chemistry, and melt inclusions studies to explain a complete sequence of formation of Sung Valley ijolites from their petrogenetic conditions to crystallization history along with the nature of parental magma and history of melts involved. Our results suggest that these ijolites crystallized in multiple stages and had some petrogenetic relation with associated carbonatites. The composition of the trapped melt inclusions in these rocks must have changed significantly during syn-to late-magmatic processes.

Geological background
The oval-shaped ultramafic-alkaline-carbonatite complex (UACC) of Sung Valley is hosted by the Shillong Plateau of Meghalaya in northeastern India (Fig. 1a, b). The Sung Valley UACC is related to the Kerguelen plume (Veena Page 3 of 18 Choudhary et al. Progress in Earth andPlanetary Science (2022) 9:71 et al. 1998;Ray et al. 1999Ray et al. , 2000Srivastava and Sinha 2004;Srivastava et al. 2005;Srivastava 2020). This UACC intrudes the Archean gneiss, schist, and the Proterozoic Shillong Group rocks in the Shillong Plateau (Ray and Pande 2001). The Sung Valley UACC and several other alkaline intrusive bodies are hosted by the N-S trending Um Ngot lineament, which is genetically related to the Ninety-East Ridge in the Indian Ocean (Gupta and Sen 1988). The Sung Valley UACC crystallized within the age range of 101-115 Ma (Ray and Pande 2001;Srivastava and Sinha 2004;Srivastava et al. 2005Srivastava et al. , 2019. This UACC is composed of ultramafic, alkaline, and carbonatite rocks. Ijolites occur as a ring dyke in the Sung Valley UACC, whereas the melilitolite, nepheline syenite, and carbonatites that occur as oval-shaped bodies and small dykes within the pyroxenite and peridotite/serpentinized peridotites (Fig. 1b). Sung Valley pyroxenites are dominantly clinopyroxene rich and show the presence of chalcopyrite of low-temperature origin as the main sulfide mineral . Carbonatites are exposed mostly in the southern part of the Sung Valley UACC (Srivastava and Sinha 2004). Sung Valley ijolites, the third most abundant rock type after pyroxenites and peridotites, are mostly coarse-grained with wide textural variation from porphyritic to poikilitic. The emplacement date of these Sung Valley ijolites is 115.1 Ma (Srivastava et al. 2005), whereas in situ U-Pb SIMS of perovskite from these ijolites yielded an age of 104.0 ± 1.3 Ma (Srivastava et al. 2019). These ijolites occur as a ring dike and the petrological and geochemical studies carried out in the past suggest that these rocks, along with the associated ultramafic (peridotite and pyroxenite) and carbonatite rocks, are formed from the batches of primitive magma with a distinct magmatic affinity such as olivine melilitites, basanites, and carbonatites, and these batches of magma evolved independently (Melluso et al. 2010). These primitive magmas were derived at a pressure greater than 2.5 GPa from a metasomatically enriched carbonated peridotite (Srivastava and Sinha 2004). The whole-rock geochemical characteristics shown by these ijolites are significantly different from nepheline syenite, suggesting that the nepheline syenites are not formed through the fractionation of ijolites (Srivastava and Sinha 2004).
of 300-400 statistically. A high-resolution petrographic microscope was used to take the photomicrographs of each rock-thin section. Crystal margins were outlined using the vector drafting tool of CorelDraw. Software Image-J was used to analyze the mineral outlines and related parameters of the clinopyroxene crystals, such as area, perimeter, and length. Finally, the CSD of the clinopyroxene crystals was calculated with the program CSD-Corrections 1.6 (Higgins 2000).

Electron probe micro-analyses (EPMA)
Polished rock thin sections of 0.03 mm thickness were prepared for petrographic and electron microprobe analyses. For mineral chemistry, four representative samples were analyzed using a CAMECA SXFive Electron Probe Micro Analyzer instrument at SERB-IRHPA National Facility, Department of Geology, Banaras Hindu University, India. For quantitative mineral chemistry analyses, wavelengthdispersive spectrometry and a LaB6 source were used. During the analysis, the instrument was operated at an acceleration voltage of 15 kV and focused beam current of 10 nA having a diameter of 1 µm. The following crystals, such as TAP (thallium acid phthalate), LPET (large pentaerythritol), and LLIF (large lithium fluoride) were used for the measurements. The natural mineral standards such as diopside, forsterite, almandine, albite, and orthoclase supplied by CAMECA-AMETEK were used for calibration and quantification. A precision better than 1% for major element oxides from the repeated analyses of standards was achieved during the analysis. The representative mineral chemical data are given in Tables 2 and 3.

Raman spectroscopy
The compositions of melt and minute crystals in the untreated inclusions from four representative samples were obtained by laser-excited Raman spectrometry using Horiba Jobin Yvan Lab Ram HR Laser Raman Micro Probe in Raman and Fluid inclusion lab at Wadia Institute of Himalayan Geology (WIHG), Dehradun. All the Raman spectra were generated using 100X objectives and with the 514 nm laser of Argon ion (Ar + ) source with a dispersion of a fixed holographic grating 1800 lines/mm. The laser spot size at the time of analysis was ~ 2 µm. To obtain a better signal-to-noise ratio, repeated spectra were recorded in the 100-4000 cm −1 region. Standard silicon was used for the calibration of the instrument. Raman shift at 520.59 cm −1 during the calibration with standard silicon was achieved. Other parameters during the analyses were taken as follows: acquisition time ~ 5-15 s, accumulations ~ 2, and laser power ~ 10 mW.

SEM-EDX
Back-scattered electron (BSE) imaging and energy-dispersive X-ray (EDX) qualitative analyses on exposed melt inclusions were acquired using Carl Zeiss SMT EVO 40 Series-Scanning Electron Microscope (SEM) with EDX, equipped with LaB6 cathode, at Wadia Institute of Himalayan Geology, Dehradun. BSE images were acquired using 20 kV accelerating voltage and beam current of 3-6 nA at different magnifications. The concentration of oxides (wt %) of the elements was analyzed by EDS attachment using the QUANTAS software.

Petrography
A total of ten samples were analyzed in this study. All the fresh ijolite samples were collected from different locations in the Sung Valley UACC (Fig. 1b). Collected ijolite samples are dominantly composed of nepheline and clinopyroxene, mainly aegirine-augite (Aeg-Aug) and diopside (Fig. 2a). These ijolites are medium to coarsegrained and predominantly show hypidiomorphic textures ( Fig. 2a, b, d). The poikilitic texture is also present in sample IJ-9/4, where clinopyroxene crystals are enclosed in large optically continuous nepheline crystals (Fig. 2c). The accessory mineral phases include apatite and titanite (Fig. 2b,c,d). Minor opaque minerals include magnetite, pyrite, pyrrhotite, and chalcopyrite (Fig. 2d). The representative mineral chemistry data for the various phases in these ijolites are given in Tables 2 and 3, and the results are detailed in the following subsection.

Crystal size distribution
Four representative samples (IJ-9/2, IJ-11/2, IJ-9/4, and IJ-11/1) were chosen for CSD analysis. The CSD results related to the clinopyroxene (aegirine-augite and diopside) crystals, such as area, perimeter, and length, were imported into the CSDCorrections 1.6 software (Higgins 2000) to make the stereological corrections (Fig. 3). All the calculations were performed assuming a massive fabric with 0.8 roundness on the blockellipsoid scale. Clinopyroxene crystals contacting the edge of the sample area were not taken into consideration during analysis, as these crystals cannot be representative of complete crystals. The resultant 3D crystal size distributions are given in Table 1. A total of ~ 1475 clinopyroxene grains were analyzed from the representative samples, and crystal size ranged from microns to millimeter scale (Table 1). The crystal size of clinopyroxene grains ranges from 0.08 to 5.65, 0.12 to 7.5, 0.03 to 1.16, and 0.13 to 9.3 mm in the samples IJ-9/2, IJ-11/2, IJ-9/4, and IJ-11/1, respectively (Table 1). The CSD plots between the natural logarithm of crystal population density Ln(n) vs. crystal length (mm) in samples no. IJ-9/2, and IJ-11/1 show a concave upward trend starting with a slight convex upward kink (Fig. 4a, d). Sample no. IJ-11/2 exhibits concave upward CSD (Fig. 4b). The CSD of clinopyroxene in a sample IJ-9/2 shows an inflection point at ~ 2 mm, whereas sample IJ-11/2 shows inflection points at 2 and 7 mm (Fig. 4a, b). Sample no. IJ-9/4, which exhibits poikilitic texture, shows the most distinct CSD profiles with marked differences in the slope nearly a straight line (Fig. 4c). The studied ijolites do not show any signatures of deformation such as annealing, dislocation, or diffusion of clinopyroxenes (Fig. 2).

Mineral composition
The EPMA was carried out on clinopyroxene and titanite of representative ijolite samples. The chemical analysis of clinopyroxene is given in Table 2. Analyzed clinopyroxene grains include both diopside and aegirine-augite (Fig. 3a). The average endmember composition of these clinopyroxene grains is (Di 53-84 Hd 10-32 Ae 4-16 ) (  (Table 2). In terms of International Mineralogical Association (IMA) recommended Q-J parameters (Morimoto et al. 1988), which are based on six oxygens, the analyzed clinopyroxene spots (Table 2) largely belong to jadeite-free quadrilateral pyroxene (Fig. 3b). The Q-J plot shows that these clinopyroxenes are crystallized at low-pressure conditions possibly at crustal depth. Mg# in clinopyroxene shows a positive correlation with Ti content (Fig. 3c).
The clinopyroxene in studied ijolites show enrichment in diopside and to some extent, hedenbergite but are low in aegirine content, which shows that these are unevolved compared to the clinopyroxene in ijolites from Kerimasi and Fen complex ( Fig. 3a) (Church 1996;Mitchell 1980;Káldos et al. 2015). These characteristics, such as almost jadeite-free composition and low aegirine content of clinopyroxene from studied ijolites, suggest that these are early crystalizing phases (Fig. 3a, c) (Káldos et al. 2015). The composition of titanite such as CaO varies from 26.88 to 27.89, TiO 2 from 36.15 to 38.14, FeO t from 1.04 to 1.60, and SiO 2 from 29.51 to 30.15 wt.% (Table 3). A minor amount of trace elements oxides was also observed, such as La 2 O 3 ranging from 0.02 to 0.17 and Nb 2 O 5 from 0.44 to 1.07 wt.% (Table 3). Fe/Al ratios (> 0.5) in the studied titanite (Table 3) point toward their derivation from silicaundersaturated magma forming plutonic rocks such as Page 6 of 18 Choudhary et al. Progress in Earth and Planetary Science (2022) 9:71 Sung Valley ijolites in the present case (Kowallis et al. 2022).

Melt inclusions and mineral inclusions petrography
A diverse population of ~ 90 melt inclusions, hosted by clinopyroxene (diopside and aegirine-augite) and titanite, was observed in the same four samples, which were used for CSD analysis of Sung Valley ijolites. These melt inclusions were recognized from transmitted light microscopy during melt inclusion petrography (Fig. 5).
Most of the melt inclusions are elongated oval in shape ranging between 2 and 15 µm in size (Fig. 5c, e, f ). Some are irregular and rounded to subrounded ranging from 3 to 6 µm in size (Fig. 5g, h). One sample IJ-9/2, also hosts some monophase mineral inclusions along with the melt inclusions in clinopyroxene (Fig. 5d). Melt and mineral inclusions hosted by titanite, were also observed in one sample IJ-9/4 ( Fig. 5a, b). Mineral inclusions in titanite are granular in appearance and of rhombus shape ranging from 2 to 8 µm in size (Fig. 5a). These inclusions occur in a secondary trail, which terminates right at the boundary of titanite crystal (Roedder 1979). Melt inclusions hosted by titanite in the same sample IJ-9/4 are irregular to subrounded with 2-5 µm in size (Fig. 5b). Melt inclusions in clinopyroxene and titanite typically show some distorted/ deformed shrinkage bubbles probably compressed to the wall in the inclusions. This also is a typical behavior of carbonatite-type melt inclusions (Golovin et al. 2020). However, most of the clinopyroxene-hosted melt inclusions contain cavities instead of shrinkage bubbles. The absence of any decrepitation haloes in the melt inclusions suggests that the original compositions of trapped phases are intact and not altered due to any change in ambient pressure conditions after entrapment (Fig. 5). Melt inclusions observed in clinopyroxenes are present in the cores of these grains, which suggests that they are essentially primary. However, many of them could also be secondary in nature as they distinctly form linear and planar alignments within the crystals (Fig. 5). Despite being hosted by the cores of the crystals, these inclusions might have been entrapped at the intermediate stage of the minerals' growth. BSE images of melt inclusions do not show any bubbles (Fig. 6). However, the empty cavities could possibly be representing the places of the former fluid phase, which probably existed earlier in these melt inclusions and escaped during exposure. High-resolution BSE images of the melt inclusions show the presence of subhedral to euhedral daughter crystals that are 2 to 10 µm  (Morimoto et al., 1988) showing that all the points of clinopyroxene (aegirine-augite and diopside) composition fall into the Quad field. c Scatter plot between Mg# and Ti in Clinopyroxene showing the depletion of Ti along the line of fractionation of Clinopyroxene Page 8 of 18 Choudhary et al. Progress in Earth and Planetary Science (2022) 9:71 in size in these melt inclusions (Fig. 6). A significant space of the melt inclusions is occupied by these daughter crystals (Fig. 6). A very limited variation in volume proportion displayed by the crystals in these melt inclusions testifies that these are daughter minerals and not accidentally trapped minerals ( Fig. 6) (Anderson et al. 2003).

Raman spectroscopy of the mineral and melt inclusions
Raman Spectroscopy was used to investigate the phase composition of the trapped melt inclusions in clinopyroxene and titanite. Mineral inclusions hosted by titanite in sample IJ-9/4 were identified as calcite, where the Raman spectrum shows characteristic symmetric stretching vibration v 1 with a very strong band at 1086.8 cm −1 (Frezzotti et al. 2012) (Fig. 5a). Melt inclusions in titanite in sample IJ-9/4 show that they have calcite accompanied by graphite in the shrinkage bubble, where calcite shows its characteristic symmetric stretching vibration v 1 with a very strong band Crystal size (mm)  (Fig. 5b). Special care was taken during the recording of Raman spectra of graphite to make sure that the D and G bands do not appear due to burning of the sample wafer by laser. So all the spectra of graphite were recorded using lesser laser power. Furthermore, to avoid any contamination due to carbon coated on the samples during EPMA, separate wafers of the samples were used during Raman spectroscopy. Asterisks in Raman spectra are assigned to the host peaks of titanite (Fig. 5a, b). The carbonate component in the melt inclusions hosted by clinopyroxene in this sample (IJ-9/4) was also identified as calcite, which again shows a very sharp peak at 1086.5 cm −1 showing symmetric stretching vibration v 1 with a very strong band, strong (s) band at 157 cm −1 , and medium weak (mw) band at 713 cm −1 (Fig. 5c) (Frezzotti et al. 2012). The presence of apatite in the studied melt inclusions was also confirmed where the Raman spectrum shows a very strong band at 964 cm −1 (Fig. 5c). Mineral inclusions hosted by clinopyroxene in sample IJ-9/2 were identified as rutile showing its very strong (vs) bands at 443 cm −1 and 609.5 cm −1 (Frezzotti et al. 2012) (Fig. 5d). The carbonate component in melt inclusions in this sample (IJ-9/2) and two other samples IJ-11/1 and IJ-11/2 also showed symmetric stretching vibration v 1 at, 1085.7, 1086 and 1085.4 cm −1 , respectively, which were identified as calcite (Frezzotti et al. 2012) (Fig. 5e-g). Melt inclusions hosted by clinopyroxene in sample IJ-11/2 were identified as containing aphthitalite (K-Na sulfate), showing its very strong band at 985.6 cm −1 , medium band at 448 cm −1 along with the CO 2 showing two bands of Fermi doublet at 1278.4 and 1384 cm −1 (Frezzotti et al. 2012) (Fig. 5h). Host clinopyroxene in all the samples shows strong Raman bands at ~ 1010 and 667 cm −1 (Thompson et al. 2005).
It is worth noting that during the Raman spectroscopy of these melt inclusions, the spectra of glass were not detected. However, in the present case, it is reasonable because the carbonatitic melts cannot be quenched to a glass (owing to very low viscosity) and inevitably crystallize (Moine et al. 2004;Giuliani et al. 2012;Kamenetsky et al. 2014;Chayka et al. 2021). All the Raman spectra of inclusions were acquired repeatedly, and representative spectra are presented in Fig. 5.

SEM-EDX of the melt inclusions
SEM-EDX analysis was carried out to determine the chemical composition of daughter phases associated with carbonate components in the exposed melt inclusions hosted by clinopyroxene (diopside and aegirine-augite). Points for EDX analysis were selected with the help of BSE images (Fig. 6). Representative oxides data for the daughter phases in studied melt inclusions are given in the Additional file 1: Table S1, and the results are detailed below. The daughter crystals in these carbonate-silicate melt inclusions were identified as phlogopite, andradite garnet, diopside, and magnetite, similar to the phases observed in melt inclusions in ijolites by Andreeva et al. (2007) and Sekisova et al. (2015) (Fig. 6) Similarly, diopside also may not contain ~ 5 to 7% Na 2 O like omphacite, which is certainly not possible to occur in these ijolites. Neither can it have 2 to 3% K 2 O, which is only likely in ultra-high-pressure clinopyroxene. We infer that these signals may also be coming from the phlogopite in these melt inclusions due to contamination from the neighboring phase (Fig. 6)   Page 11 of 18 Choudhary et al. Progress in Earth and Planetary Science (2022) Choudhary et al. Progress in Earth and Planetary Science (2022) 9:71 inclusions are dominated by carbonates, whereas silicates are subordinate (Fig. 5, 6). This further confirms that the parental melt for the inclusions was carbonatitic with silicate component (Golovin et al. 2018(Golovin et al. , 2020. Also, the fact that these inclusions are consistently made up of a variety of crystals of different phases is the key evidence that they are not just mineral inclusions, which are most commonly monophase, but rather represent aliquots of trapped liquid that crystallized within these melt inclusions (Fig. 5, 6). The quantitative EDX data of the analyzed phases are not significantly contaminated by the host as there are almost no alkalies in the host as shown by mineral chemistry data (Table 2), whereas we have a significant amount of alkalies in the daughter phases (Additional file 1: Table S1). Therefore, the studied inclusions are carbonate-silicate (carbonatite-like) melt inclusions (Figs. 5, 6). Based upon the Raman spectroscopy and SEM-EDX data and by virtue of the occurrence of typical liquidus phases of calcio-carbonatitic liquids, i.e., calcite, apatite, diopside, phlogopite, garnet, and magnetite as a fully crystallized melt, these melt inclusions can also be termed as "nano-calciocarbonatites" (Figs. 5, 6).

Crystallization history of Sung Valley ijolites
A complete absence of annealing, dislocation, and diffusion in these clinopyroxenes rules out any possibility of post-magmatic recrystallization (Fig. 2). Therefore, the textural parameters observed during CSD analysis are unaffected by any sub-solidus processes and they represent the primary magmatic character (Fig. 4). All the CSDs of the analyzed samples show curved concaved upward distribution (Fig. 4). Hence, a single population model of crystals cannot be applied. There can be various inferences based upon curved CSDs that can result from the crystallization in igneous systems. We will discuss all possible interpretations in the context of our CSD dataset ( Fig. 4; Table 1) for the Sung Valley ijolites. The first explanation for the curved CSD can be the crystal fractionation process (Marsh 1998). However, this possibility can be negated because curved CSDs that result due to the crystal fractionation process show nearly the same intercepts in CSD plots (Higgins 1996). On the contrary, in the present case, all the CSD plots show different intercepts (Fig. 4). The second explanation can be the change in growth rate as a function of crystal size, where growth rate increases with the size of crystals (Marsh 1988;Eberl et al. 2002). This argument can also be negated because a curved but concave downward CSD pattern can result through this process, as explained by Inanli and Huff (2009), unlike concave upward CSD in the present study (Fig. 4). Additionally, the cases of such variable crystal growth rates with size have not been described yet in the geological environments (Higgins 1996). The third possible explanation for producing curved CSDs can be the magma-mixing hypothesis (Higgins 1996), where melts of different compositions interact physically and chemically, which leads to the disequilibrium condition giving rise to partial resorption e.g., corona texture in clinopyroxene. A complete absence of any such texture (Fig. 2) in the present study further negates this possibility. The fourth possible explanation that can produce the concave upward CSDs is fines destruction (Marsh 1988). In this process, larger crystals are more likely to remain at depth, whereas the small crystals rise at shallower levels in the magma chamber and are resorbed (Higgins 2002). During the fines destruction process, nutrients of the resorbed smaller crystals are fed to the larger crystals (Higgins 1996(Higgins , 2011. The abundance of large crystals of clinopyroxene compared to smaller ones also suggests that crystal growth rate dominated over the nucleation along with the simultaneous dissolution of smaller clinopyroxene crystals (Fig. 4). Such coarsening in the CSD data, which is also substantiated with petrography (Fig. 2), may have a bearing on the high-temperature magmatic storage. Fifth and the last possible explanation to give rise to the curved CSDs with kinks is the abrupt changes in the crystallization environment e.g., cooling rate (Marsh 1988;Armienti et al. 1994), which can also control the nucleation rate (Kamacı and Altunkaynak 2019). This inference is further supported by the similar composition of both the smaller and larger clinopyroxene crystals of Sung Valley ijolites (Table 2), which can be related to the modification of the crystallization environment without any significant changes in the composition of these crystals in the present case (Higgins 2011). The last two possible explanations seem to best fit in the context of observed curved concave upward CSDs in Sung Valley ijolites (Fig. 4). Furthermore, changes in the crystallization environment (cooling rate, nucleation, and growth rate) can incorporate mixed crystals of different generations, which in turn can give rise to concave upward CSD (Morgan et al. 2007). The clinopyroxene CSD in the Sung Valley ijolites can be separated into two categories, where samples IJ-9/2, IJ-11/2, and IJ-11/1 (Fig. 4a, b, d) show curves with inflection points, whereas sample IJ-9/4 ( Fig. 4c) displays nearly a straight line. These observations suggest that clinopyroxene in IJ-9/4 grew in a single stage, and in the rest of the samples IJ-9/2, IJ-11/2 and IJ-11/1, these crystals grew in multiple stages (Wang et al. 2019) (Fig. 4).

Petrogenetic conditions of Sung Valley ijolites: existing models
The Sung Valley UACC rocks are derived from the partial melt of carbonated peridotite at a pressure greater than 2.5 GPa (Srivastava and Sinha 2004;Srivastava et al. 2005). Geochemical studies of the different rock units from the Sung Valley UACC further suggest that these rocks were formed by discrete batches of primitive magmas with different magmatic affinities, and these magmas probably derived from the same source yet evolved independently (Srivastava et al. 2005;Melluso et al. 2010). On the basis of whole-rock geochemical and isotopic studies, Veena et al. (1998) and Ray et al. (2000) concluded that Sung Valley rocks are formed from partial melting of the subcontinental lithospheric mantle, which was previously subjected to metasomatism by Kerguelen mantle plume-derived fluids. Pieces of evidence of carbonate metasomatism of the lithospheric mantle beneath Sung Valley UACC were recently delineated by Choudhary et al. (2021). Therefore, assuming a carbonated peridotite source to the parental melt of ijolite during partial melting would be appropriate, as also suggested by Srivastava et al. (2005). The mineralogical composition of Sung Valley ijolites corroborates the earlier geochemical studies (e.g., Srivastava and Sinha 2004;Melluso et al. 2010) and suggests that these rocks were formed from a magma of nephelinitic affinity, which was Page 14 of 18 Choudhary et al. Progress in Earth and Planetary Science (2022) 9:71 derived from partial melting of a carbonated peridotite. Srivastava et al. (2005) obtained a high concentration of LREE in chondrite-normalized rare-earth patterns in these ijolites, such concentration of LREEs points toward a low degree of partial melting during the derivation of parental melt to the Sung Valley ijolites. Isotopic studies on Sung Valley ijolites carried out by Ray et al. (1999) suggest that parental magma to these ijolites was derived from either a low U/Pb source or interacted with a low U/Pb mantle reservoir. Our mineral chemical data show a jadeite-free composition of the analyzed clinopyroxene in the present study ( Fig. 3b; Table 2), which supports the inference that this melt was crystallized at crustal depth after its origin at greater than 2.5 GPa pressure (Srivastava and Sinha 2004).

Composition of melt inclusions: nature of parental ijolite melt and presence of "nano-calciocarbonatites"
The carbonate-silicate melt inclusions observed in the present study, as well as the models proposed by earlier studies, suggest that these ijolites were formed from a carbonated olivine-nephelinite magma (Melluso et al. 2010), which was probably derived from partial melting of carbonated peridotite (Srivastava and Sinha 2004). The silica-undersaturated nature of the parental magma of these ijolites is further supported by Fe/Al ratios (> 0.5) in the titanite (Table 3) (Kowallis et al. 2022). These nano-calciocarbonatites represent a typical assemblage of liquidus phases that crystallized from calcio-carbonatitic liquid (Figs. 5,6). The predominant occurrence of carbonate component, i.e., calcite in melt inclusions in the core of these crystals suggests that these carbonates represent the pristine magma and are not derived from carbonate metasomatism of ijolite itself as suggested elsewhere (Seifert and Thomas 1995;Jones et al. 2000;Downes et al. 2002;Woolley and Bailey 2012). An abundance of carbonate-silicate melt inclusions in the Sung Valley ijolites (Figs. 5, 6) also testifies that the source rock is a carbonated peridotite, as suggested by Srivastava and Sinha (2004). The presence of disordered graphite with calcite in titanite of a sample IJ-9/4 (Fig. 5b) and mineral inclusions of rutile in clinopyroxene of sample IJ-9/2 ( Fig. 5d), together indicate the fluctuation of fO 2 , invoked by redox reactions in the lower crust during the entrapment of these inclusions (Fig. 5). We opine that during the entrapment of these inclusions, tetravalent Ti 4+ cations from the parental magma consumed enough oxygen provided by carbonate-silicate melt (enriched in CO 2 ) to form rutile (Fig. 5d), which in turn reduced the residual CO 2 to graphite (Fig. 5b). The presence of CO 2 can be justified here because the carbonate components of these melts are enriched in CO 2 and H 2 O (Jones et al. 2013;Choudhary et al. 2021). A positive correlation between Mg# and Ti in clinopyroxene (Fig. 3c) also indicates that Ti is concentrated in the form of rutile inclusions during the early fractionation of the clinopyroxene at high temperature. The Ti content kept decreasing with Mg# along the clinopyroxene fractionation line (Fig. 3c). However, the parental magma of these ijolites got saturated in Ti later, which resulted in the formation of titanite (Fig. 2c,  d). The occurrence of aphthitalite (K-Na sulfate) (Fig. 5h) in inclusions marks the presence of oxidized S-rich fluid in the crystallizing magma (Bataleva et al. 2018). However, the occurrence of these sulfates is common in carbonate-silicate melts (Chayka et al. 2021).
Melt inclusions can record evidence of liquid immiscibility in magmas from a variety of different tectonic settings (Thompson et al. 2007;Panina and Motorina 2008;Mitchell 2009;Kamenetsky and Kamenetsky 2010;Sekisova et al. 2015). It is noteworthy that the carbonate component in all the carbonate-silicate melt inclusions in clinopyroxene and titanite in the studied ijolites are purely calcite (Figs. 5,6), and carbonatites associated with these ijolites in the Sung Valley UACC are also purely calciocarbonatites belonging to the sovitic carbonatite group (Choudhary et al. 2021). Therefore, the first possible explanation of the occurrence of studied melt inclusions can be the separation of carbonate melt from silicate melt or silicate-carbonate melt immiscibility at the waning stages of the formation of the Sung Valley UACC. A number of previous studies on the Sung Valley UACC have also suggested the immiscibility model among silicate and carbonate rocks (Viladkar et al. 1994;Veena et al. 1998;Sen 1999;Ray et al. 2000). However, lower Ba/La ratios in carbonatites in comparison to associated silicate rocks, as well as a paucity of immiscible droplets of carbonate liquid in the associated silicate rocks negated this possibility (Hamilton et al. 1989;Srivastava and Sinha, 2004;Melluso et al., 2010). Moreover, these authors did not completely rule out the immiscibility model and suggested the possibility of the liquid immiscibility of Sung Valley carbonatites from a more primitive liquid. However, a complete absence of any alkali carbonates, such as nyerereite, shortite, and natrite, in the studied melt inclusions questions the possibility of the immiscibility process (Fig. 5). However, calcite being a typical liquidus phase in a carbonatite melt cannot be the only phase resulting from crystallization of such a melt. Although the possibility of silicate-carbonate melt immiscibility cannot be completely ruled out, this calcite discrepancy and loss of alkalies from carbonates need to be addressed through some alternate models.
The second possible explanation for the occurrence of these carbonate and silicate phases can be their accidental entrapment, just as solid crystals in these melt inclusions. However, this process is also unlikely as the ijolites Page 15 of 18 Choudhary et al. Progress in Earth and Planetary Science (2022) 9:71 under investigation are magmatic rocks that are devoid of any post-magmatic recrystallization. Furthermore, a very limited variation in volume proportion in the crystals in these melt inclusions advocates against such a possibility ( Fig. 6) (Anderson et al. 2003). Additionally, the occurrence of a variety of these crystals of different phases suggests that they are the aliquots of trapped crystallizing liquid within these melt inclusions (Figs. 5, 6).
The predominant occurrence of calcite as the only carbonate phase also points toward a calcite-normative system in the studied melt inclusions. This could be the third possible explanation and such calcite-normative system of the studied carbonate melts could result through either extreme fractionation of carbonate-silicate melt or continuous infiltration of alkali-poor carbonate melt during the crystallization of clinopyroxenes and initially present alkali-rich carbonate melt. This infiltration process could take place through magmatic metasomatism during syn to late magmatic processes (Mathez 1995).
Generally, the alkali carbonates may undergo a prompt replacement by calcite even at low-temperature (Zaitsev and Keller 2006). During such processes, the leaching of alkalies (e.g., Na, K) takes place from soluble alkali carbonatite liquids through dealkalization with the formation of stable calcite instead (Le Bas 1981;Chen et al. 2013;Chayka et al. 2021). Therefore, the fourth possible explanation could be the dealkalization of the initially trapped alkaline carbonates to calcite. However, this possibility of such dealkalization appears to be controversial in the present case as this process requires an open system for the trapped melt inclusions and also the presence of some external fluid during late-stage crystallization or sub-solidus conditions of these ijolites. Only a single titanite crystal in one sample (IJ-9/4), hosts secondary mineral inclusions of calcite that points toward occurrence of external fluid (Fig. 5a). However, paucity of such secondary inclusions and absence of any sub-solidus processes indicated by CSD analysis (Fig. 4) argue against this possibility. It is noteworthy that sometimes the melt inclusions, appearing as primary, do not behave as an absolutely closed system, but there could still be micro-/ nano-fissures through which late-or post-magmatic fluids enter the inclusion and modify their composition. These micro-fissures may be healed again due to solidstate diffusion processes in crystals, leaving no evidence of the alteration event, as similarly reported by Chayka et al. (2020Chayka et al. ( & 2021. Based on this argument, the process of dealkalization cannot be ruled out completely. The fifth and last possible scenario for the observed calcite discrepancy in the studied melt inclusions could be the redistribution of alkalies to the coexisting silicate phases within the inclusions, leaving only calcite as a carbonate phase (Fig. 6). The alkali contents of the host clinopyroxene is very low (Table 2). Therefore, the alkalies from the initial alkaline carbonate melt could be redistributed to the daughter silicate phases inside these carbonate-silicate melt inclusions (Fig. 6). This inference is also supported by the presence of a notable amount of alkalies in the daughter diopside, phlogopite, and the unidentified carbonate-silicate phases (Additional file 1: Table S1). Therefore, it may be plausible that the daughter crystals took up the alkalies from the initial carbonate melt, and calcite was left as the final carbonate phase (Figs. 5, 6). Finally, out of all the possible scenarios, first and the last three seem to be responsible for a predominant occurrence of calcite in these carbonate-silicate melt inclusions or nano-calciocarbonatites. These melt inclusions in the studied ijolites clearly indicate the activity of complex alkali-bearing carbonate-silicate-phosphate-sulfate melts during the crystallization of these ijolites (Figs. 5, 6).

Conclusions
Based on the results obtained in this study, the following conclusions can be drawn: • Sung Valley ijolites dominantly show concave upward CSD patterns caused by the fines destruction process and abrupt changes in the crystallization environment; also, these rocks were crystallized in multiple stages. • The presence of calcite and alkali silicates in carbonate-silicate melt inclusions in the clinopyroxene points toward silicate-carbonate immiscibility during the formation of Sung Valley ijolites and carbonatites. • The predominant occurrence of calcite as the only carbonate phase in the studied melt inclusions suggests that these inclusions could also be a result of (i) calcite-normative system in these melts resulting from either extreme fractionation of carbonate-silicate melt or continuous infiltration of alkali-poor carbonate melt, (ii) dealkalization of the initially trapped alkaline carbonates in the presence of external fluid that entered through micro-/nanofissures which subsequently healed, (iii) redistribution of alkalies to the coexisting silicate phases (diopside, garnet, phlogopite and carbonated silicate) within the inclusions, leaving calcite as the only carbonate phase.