Next Article in Journal
A Nonlinear Gradient-Coiling Metamaterial for Enhanced Acoustic Signal Sensing
Next Article in Special Issue
Supramolecular Assemblies of Melamine-2-Thiobarbiturate and Melamine-Barbiturate-2-Thiobarbiturate: Experimental and Theoretical Studies
Previous Article in Journal
Study of Microstructure, Crystallographic Phases and Setting Time Evolution over Time of Portland Cement, Coarse Silica Fume, and Limestone (PC-SF-LS) Ternary Portland Cements
Previous Article in Special Issue
Zn(II) Metal–Organic Frameworks with a Long Spacer Ligand and a Tricarboxylate Coligand
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Compositional Variation and Crystal-Chemical Characterization of a Watermelon Variety of Tourmaline from Anjanabonoina, Central Madagascar

1
Department of Earth and Geo-Environmental Sciences, University of Bari “Aldo Moro”, 70125 Bari, Italy
2
Department of Earth Sciences, Sapienza University of Rome, Piazzale Aldo Moro 5, 00185 Rome, Italy
*
Author to whom correspondence should be addressed.
Crystals 2023, 13(8), 1290; https://doi.org/10.3390/cryst13081290
Submission received: 27 July 2023 / Revised: 11 August 2023 / Accepted: 16 August 2023 / Published: 21 August 2023
(This article belongs to the Special Issue Feature Papers in Crystals 2023)

Abstract

:
A polychrome tourmaline crystal from Anjanabonoina pegmatite (Madagascar) was characterized using a multi-analytical approach. The sample showed a complex concentric zoning and a wide range of colors typical of the variety known as “watermelon”. The sample was cut perpendicularly to the c axis. The basal slice exhibits a rim characterized by narrow, differently colored layers parallel to the prism faces and a relatively homogeneous triangular core. Four main pronounced color zones were identified from the rim to core: a dark green rim (M1RVS); a pale green rim (M1RVC); a pale pink rim (M1CR); and a brownish yellow core (M1CG). Compositional variations in the basal slice were studied by scanning electron microscopy and electron microprobe analyses (WDS mode). The Li content was determined via micro-laser-induced breakdown spectroscopy. To deeply characterize the sample, single crystal structure refinement was also performed on fragments extracted from the four zones. The results show that the polychrome tourmaline sample consists of two different species: the three outer zones are Mn-rich fluor-liddicoatite, whereas the inner zone is Mn-rich fluor-elbaite. The structural and compositional characterization of the color zoning shows that each step of the tourmaline growth is related to a change in the geological environment.

1. Introduction

Tourmalines represent the most widespread borosilicate supergroup of minerals in the Earth’s crust. The incredible increase in interest in these minerals in recent years is due to the wide chemical and structural variability that can be found in the same specimen [1]. These minerals typically occur in evolved granitic rocks, especially pegmatites, but also leucogranites and associated veins. Tourmalines are typical primary accessory phases in peraluminous granites, but not in metaluminous granites, where they may occur as a late hydrothermal phase. The crystallization of mainly elbaite-shorl-liddicoatite tourmalines in late-stage igneous and hydrothermal rocks is favored by enrichment in the residual granitic melt of incompatible elements, like boron and lithium. In particular, hydrothermalism can also lead to the crystallization of tourmaline as a metasomatic product in the aureole surrounding granitic plutons. This process, called ‘tourmalinization’, refers to the replacement of ferromagnesian minerals, like biotite, by tourmaline.
Dravite-schorl tourmalines occur as accessory phases mainly in metamorphic rocks as a recrystallization product of detrital tourmaline present in the sedimentary protolith or due to boron metasomatism.
The multiple occurrences of tourmaline in different geologic environments reflect a large stability range in both pressure and temperature. It ranges from diagenetic/epigenetic through high-grade metamorphic conditions up to magmatic conditions. The presence of tourmalines in diagenetic zones and hydrothermal deposits suggests that the lower limit for the formation of these minerals is defined approximately at 150 °C and 100 MPa [2,3]. The upper thermal stability has been observed in high temperature experiments between 725 °C and 950 °C depending on the pressure and composition [4,5,6,7,8,9].
The general formula of the tourmaline supergroup can be written as:
XY3Z6(T6 O18)(BO3)3V3W
where the most common constituents are:
X = Na, Ca, K or □ (vacancy)
Y = Al, Li, Fe, Mg, Mn, V, Cr, Ti, etc.
Z = Al, Mg, Cr, V, Fe
T = Si, Al, B
B = B
V = OH, O
W = F, O, OH [1]
The letters in the formula (X, Y, Z, T, and B) represent groups of constituents at the [9]X, [6]Y, [6]Z, [4]T, and [3]B crystallographic sites. The letters V and W represent groups of anions at the [3]O3 and [3]O1 sites, respectively. The H atoms occupy the H3 and H1 sites, which are related to O3 and O1, respectively [10].
The B site is fully occupied by boron and the T site is predominantly occupied by Si but also can contain Al and B. The Y and Z sites can incorporate many elements with different sizes and charges.
Tourmaline minerals belong to the cyclosilicate subclass. Paraphrasing from Bosi [11], the crystal structure may be described as follows: The nine-coordinated XO9 antiprism, which may also be vacant, and T6O18 ring combine with two sets of three octahedra YO6; a [Y3O15] triplet of octahedra caps the XO9 polyhedron toward the +c axis, with the other [Y3O12] capping the [T6O18] ring of tetrahedra toward the −c axis. The BO3 groups are oriented sub-parallel to (0001) and lie between the rings of tetrahedra. The structural arrangement of [T6O18], XO9, [Y6O18], and (BO3)3 form “islands” that are stacked in columns along the c axis. These islands are attached to one another along the a and b crystallographic axes by spiral chains of ZO6 octahedra, which also extend along the c axis via a 31 triad screw axis. The three-dimensional framework of the tourmaline structure is therefore given by the screw-like arrangement of ZO6.
Due to the chemical composition flexibility, tourmaline is classified as a supergroup that comprises thirty-nine different species approved by the International Mineralogical Association.
Thanks to the great variability of the substitutions occurring in tourmalines, these minerals show very different colors and often appear zoned. The color zoning mainly reflects the growth condition evolution under which the minerals were formed. The changes in chemical composition of the growth medium and variations in temperature and pressure force the mineral to readjust its composition during its formation, creating, in some cases, abrupt modification of color and composition in the same crystal [11].
The stunning array of colors that characterize tourmalines makes them highly appreciated in the gemological trade. The most important varieties used as gems belong to elbaite (indicolite, rubellite, verdelite, etc.) but also, for example, species such as dravite, fluor-liddicoatite, and, rarely, schorl.
A very singular variety of tourmaline heavily zoned perpendicularly to the c axis, with a dark green rim and a usually pink or multicolored core, is named watermelon.
Generally, this variety consists of a liddicoatitic or elbaitic composition, even so, to the best of our knowledge, the literature contains an uncomplete crystal chemical characterization of watermelon tourmaline. Using only chemical analyses, an elbaitic composition was found by Johnson et al. [12] for watermelon tourmalines from Zambian deposits and by Roda-Robles et al. [13] for watermelon tourmalines from Berry–Havey pegmatite (Maine, Oxford, MS, USA). An elbaitic composition was also found by Natkaniec-Nowak et al. [14] and Wesełucha-Birczyńska et al. [15] in samples of watermelon tourmaline from Paprok mine (Afghanistan) using both chemical and spectroscopic Raman analyses. Conversely, a liddicoatitic composition was found only based on optical observations on samples from Anjanabonoina district (Madagascar) and Möng Long and Molo (Myanmar) that showed the typical trigonal stars of this tourmaline species [16,17]. A coexistence between liddicoatitic and elbaitic compositions was observed by Akizuki et al. [18] and Dirlam et al. [19], both of Malagasy origin. In detail, Akizuki et al. [18] observed, by EPMA and X-ray diffraction, growth sectors of liddicoatitic composition and an outer rim zone of elbaitic composition. Similarly, Dirlam et al. [19] observed by X-ray mapping an elbaitic and liddicoatitic portion in different zones of the same sample. In both studies, light elements such as H, Li and B were calculated without any crystal structure information.
In this study, a detailed characterization of a watermelon tourmaline from Anjanabonoina was performed for the first time using chemical and structural information, including measurements of light elements like Li.

2. Materials and Methods

The studied tourmaline, labelled MADAG 1, comes from the Anjanabonoina aplite–pegmatite located in central Madagascar. The bulk sample (Figure 1a) originally showed a dark-brownish color and an elongated prismatic morphology deeply striated parallel to the c axis, with a basal section with a size of about 3.5 cm in length and 1.4 cm in diameter. From this raw sample, one basal thin section was cut and polished to perform a multi-analytical investigation. The basal slice shows the color zoning typical of the watermelon variety, consisting of a triangular core with a brownish yellowish color, an intermediate pink zone and a green rim (Figure 1b).

2.1. Occurrence

The Anjanabonoina pegmatite (55 km west-southwest of the city of Antsirabe, central Madagascar) is one of the most important sources of multicolor and giant crystal tourmalines. These pegmatites were emplaced 490 million years ago during late-stage granitic plutonism related to the Pan-African event, which occurred from 570 to 455 million years ago [20]. The host rocks consist of quartzites, schists and marbles (Itremo group, Figure 2) locally tourmalinized near the pegmatites [21]. The Anjanabonoina pegmatite shows the characteristics of both lithium–cesium–tantalum (LCT) and niobium–yttrium–fluorine (NYF) pegmatite families according to pegmatite classification by Simmons [21]; these hybrid pegmatites are very rich in rare elements which contributes to the crystals’ formation in mialoritic cavities known as pockets. Gem-bearing pockets usually represent the ultimate concentrations of exotic elements, volatiles and other fluxes present in pegmatites.
During the final stage of solidification, the fluid pressure inside this pocket may cause rupture, and other elements can be introduced from the outside (for example, Ca or Fe), which contributes to donating polychrome endings to this variety of tourmaline [21].

2.2. Experimental Methods

In the present study, a Malagasy tourmaline of the watermelon variety has been studied using a multi-analytical approach: optical microscopy, scanning electron microscopy (SEM) ((Zeiss, Oberkochen, Germany), electron microprobe Cameca, Los Angeles, CA, USA), micro-laser induced breakdown spectroscopy (μ-LIBS) (Marwan, Pisa, Italy) and single-crystal structure refinement (SREF) (Billerica, MA, USA).
Optical observations were made using a Carl Zeiss Axioskop 40 Pol polarizing microscope at the Department of Earth and Geo-environmental Sciences, University of Bari “Aldo Moro”.
The SEM microscope used was a Zeiss LEO, model EVO-50XVP, housed at the Department of Earth and Geo-environmental Sciences, University of Bari “Aldo Moro”, coupled with an X-max (80 mm2) energy dispersive (ED) silicon drift Oxford detector (Oxford Instruments, Oxford, UK) equipped with a Super Atmosphere Thin Window. The accuracy of the analytical data was also checked using several standard minerals manufactured by Micro-Analysis Consultants Ltd. (St. Ives, UK). After a preliminary chemical analysis with SEM-EDS (Zeiss, Oberkochen, Germany), high-resolution images were also acquired under variable pressure (Variable Pressure Scanning Electron Microscopy (VPSEM)) to obtain cathodoluminescence effects.
EMPA data were obtained in wavelength dispersive spectrometry mode (WDS) with a Cameca SX50 instrument available at the Istituto di Geologia Ambientale e Geoingegneria (Rome, Italy), CNR, operating at an accelerating potential of 15 kV and a sample current of 15 nA, with a 10 µm beam diameter. The following minerals and synthetic compounds were used as standards: wollastonite (Si, Ca), magnetite (Fe), rutile (Ti), corundum (Al), vanadinite (V), fluorphlogopite (F), periclase (Mg), jadeite (Na), K-feldspar (K), rhodonite (Mn), sphalerite (Zn) and metallic Cr and Cu. The counting time was 20 s for all elements. The PAP correction procedure for quantitative electron probe micro analysis was applied [24].The relative error of the data was ≤5% and the detection limits were ≤0.03%. The polished slice was investigated performing three traverses: 60 analysis points for the rim (one step point every 12 µm), 12 analysis points for the pink rim zone (one point analysis every 260 µm) and 10 analysis points for the core (one analysis point every 295 µm), which did not show any pronounced oscillatory zoning.
The LIBS method was used to perform traverses on the sample’s surface with an average of five shots for each point in the four selected zones. LIBS was applied using a double pulse Q-switched (Nd-YAG, λ = 1064 nm) laser with a 1 μs delay between the two pulses with an energy of 20 mJ per pulse at the Department of Earth and Geo-environmental Sciences, University of Bari “Aldo Moro”. The small spot size (7–10 μm) was obtained using a petrographic optical microscope (objective lens 10 × NA 0.25 WD 14.75 mm). The LIBS spectra were acquired by an AvaSpec Fiber Optic Spectrometer (390–900 nm with 0.3 nm resolution) with a delay of 3.28 μs after the second pulse and an integration time of 1 ms. The quantitative data were calculated from linear regression using the main lithium emission line intensity (670.706 nm). The calibration was performed using two standard glasses, NIST610 and FLX-SLAG1, with a Li2O content, respectively, of 0.105% and 5.0% and a Tsilaisite sample (NaMn2+3Al6(Si6O18) (BO3)3(OH)3(OH)) from Bosi et al. [25] with a Li2O content of 2.04%, quantified by secondary ion mass spectrometry (SIMS).
Single-crystal structure refinement was performed on four fragments extracted from each of the colored zones of the bulk sample. The fragments were mounted on a glass fiber for X-ray diffraction acquisition. X-ray intensity data were collected with a Bruker KAPPA APEX-II single-crystal diffractometer (Sapienza University of Rome, Earth Science Department) equipped with a charge-coupled device (CCD) area detector and a graphite crystal monochromator using a MoKα X-ray source (0.70930 Ǻ). Crystal structure refinement was performed in SHELXL-2013 software [26]. All structures were refined in the space group R3m.

3. Results

Optical observations of the basal slice (perpendicular to the c axis) show four main pronounced color zones from the rim to the core: a dark green rim (M1RVS); a pale green rim (M1RVC); a pale pink rim (M1CR); and a brownish yellow triangular core (M1CG) (Figure 3). The light yellow zone that surrounds the brownish yellow core indicates the transition between the (M1RVS), (M1RVC) and (M1CR) zones that perpendicularly cross the prismatic growth bands and the (M1CG) core zone that cuts the side of the inner pyramidal growth sectors of the slice.

3.1. Chemical Analyses

To deeply investigate the optically observed color zoning, comparisons were made between optical and high-resolution images acquired at a high depth of field with Variable Pressure Scanning Electron Microscopy (VPSEM) (Figure 4). This correlation shows that the slice is characterized by a succession of narrow tiny layers highlighted by cathodoluminescence and this sequence also characterizes the optically homogeneous pink region.
To investigate the chemical compositional variations, three EPMA traverses were performed (Figure 3): a first traverse investigating the outer rim zone (A-B) (Figure 5), a second traverse was carried out on the pink rim zone (C-D) (Figure 6) and a third traverse was performed on the core zone (E-F) (Figure 7).
The A-B traverse of outer rim was performed for a total length of 724 µm, with 60 analysis points and a step size of 12 µm to fully characterize the entire zone.
The A-B traverse of the outer zones (M1RVS, M1RVC and M1CR) exhibits an important oscillatory zoning; the dark green rim shows a high content of MnO (up to 4.03 wt %) that at the interface with the pink zone decreases (to 0.25 wt %) then slightly increases in the pale green rim (to 2.28 wt %) (Figure 5). MnO shows a negative correlation with Al2O3 and CaO and a positive correlation with Na2O and FeO. The three rim zones show different FeO content, with the dark green outer rim characterized by a higher amount of FeO (maximum value 1.27 wt %). In the pink rim zone, the FeO content decreases (to 0.10 wt %) and slightly increases in the pale green rim (up to 0.74 wt %).
The different content of FeO in the M1RVS, M1CR and M1RVC zones, along with the MnO content, is responsible for the different colors of these zones. The rim is rich in CaO, especially in the pink rim (maximum value 4.07 wt %) and pale green rim zone (maximum value 3.95 wt %). It decreases slightly in the outer dark green rim zone (2.71 wt %), showing a negative correlation with Na2O. The fluorine content is almost constant across the entire sample (about 1.20 wt %) and TiO2 is present in very small amounts in this rim zone (about 0.26 wt %).
The C-D traverse of the pink rim zone (M1CR) exhibits an almost constant trend of the main elements. MnO is present throughout this zone with a slight decrease in the central part of the traverse (minimum value 0.28 wt %) and an increase towards the core (up to 0.69 wt %). A similar trend is observed for Na2O, which tends to increase towards the central zone (from 0.77 wt % to 1.07 wt %), while CaO tends to decrease (from 4.20 wt % to 3.25 wt %). Al2O3 shows some higher values (up to 39.26 wt %), possibly due to spot analyses of inclusions present in the sample. In this zone, the FeO quantity is slightly higher in the outer areas and decreases towards the yellow core (from 0.14 wt % to 0.07 wt %). TiO2 is present in small amounts (about 0.06 wt %) in this pink rim zone. Note that the chemical variations observed in this zone do not correspond to the thin layers observed in VPSEM images. On the other hand, the pink color of this region is almost homogeneous and, consequently, the tiny layers could be due to thin growth bands associated with very little compositional variations of elements under the detection limits of WDS analyses.
The traverse E-F from the pink rim zone to the core zone clearly shows a dramatic variation in the main elements. The light yellow transition zone exhibits a high increase in MnO content (up to 6.42 wt %) and a decrease in Al2O3 (to 37.16 wt %). At the interface between the transition zone and the brownish yellow core, it is also observed that the CaO content decreases (to 1.70 wt %), while the Na2O quantity increases (up to 2.00 wt %). The core zone has the highest amount of TiO2 (about 0.27 wt %); even though it is very small content, it seems to be sufficient to give a brownish coloration with a constant quantity of FeO (about 0.03 wt %) along this traverse.
EPMA traverses performed on the slice were integrated with LIBS traverses to detect the Li2O content in each zone.
The calibration line obtained by LIBS (Figure 8) plotting the counts of the Li emission line at 670.7 nm acquired from the standard materials versus the Li2O wt % shows an excellent fitting (r2 = 0.94). Thus, the quantification of Li2O was obtained using the equation Li2O% = 0.0004x − 3.0543, where x represents the average of the Li emission line counts at 670.7 nm measured from the sample.
The chemical composition of each of the four zones, calculated by an average of different points for each zone (Table 1) using LIBS and EPMA is reported in Table 1.

3.2. Crystal Structure Refinement

Four fragments were extracted for SREF analysis from the four homogeneous zones (Figure 9).
Structural refinements were performed with the SHELXL-2013 program [26]. The starting coordinates were taken from Bosi et al. [27]. The variable parameters were scale factor, extinction coefficient, atom coordinates, site scattering values and atom displacement factors. In detail, the X site was modeled by considering the presence of Na and Ca scattering factors and the Y site was modeled considering Li and Mn. The Z, T, B and anion sites were modeled, respectively, with Al, Si, B and O scattering factors and with a fixed occupancy of 1, because refinement with unconstrained occupancies showed no significant deviations from this value. The position of the H atom bonded to the oxygen at the O3 site in the structure was taken from the Fourier difference map and incorporated into the refinement model; the O3–H3 bond length was restrained (by DFIX command) to be 0.97 Å with the isotropic displacement parameter constrained to be equal to 1.2 times that obtained for the O3 site. Three full matrix refinement cycles with isotropic displacement parameters for all atoms were followed by anisotropic cycles until convergence was attained. Table 2 lists crystal data, data collection information and refinement details. Structural data such as fractional atom coordinates and bond lengths are available as CIFs (Supplementary materials).
It can be observed from the crystal data (Table 2) obtained by refinement and the chemical data (Table 1) that a and c unit cell parameters are positively correlated (r2 = 0.997). This correlation results in an increase in the unit cell volume due to increased Mn + Fe content, i.e., cations in octahedral coordinations of larger sizes (Figure 10).

3.3. Site Populations

In agreement with the structural refinement results, the boron content was assumed to be stoichiometric (B3+ = 3.00 apfu). Both the site scattering results and the bond lengths of B and T are consistent with the B sites being fully occupied by boron and no amount of B3+ at the T site [11], e.g., the mean bond length <T-O> varies between 1.618 and 1.620 Å. Lithium was determined by μ-LIBS. The Fe and Mn oxidation state was assumed as +2 in accordance with the structural information. The OH content and the atoms per formula unit (apfu) were then calculated by charge balance with the assumption that (T + Y + Z) = 15 apfu and 31 anions.
Following the standard site preference suggested for tourmaline [11], the resulting empirical formulae for each of the four zones are as follows:
Dark green rim zone (M1RVS):
X(Ca0.50Na0.410.09)Σ1.00 Y(Li1.39Al1.00Mn0.47Fe0.13Ti0.02)Σ3.00 ZAl6 [T(Si5.99Al0.01)Σ6.00O18] (BO3)3 V(OH)3 W[F0.59(OH)0.35O0.06]Σ1.00
Pale green rim zone (M1RVC):
X(Ca0.64Na0.300.06)Σ1.00 Y(Li1.64Al1.04Mn0.22Fe0.06Ti0.02)Σ2.98 ZAl6(TSi6.02O18)(BO3)3 V(OH)3W[F0.59(OH)0.35O0.06]Σ1.00
Pale pink rim zone (M1CR):
X(Ca 0.52Na0.350.13)Σ1.00 Y(Li1.65Al1.19Mn0.15Ti0.01)Σ3.00 ZAl6 [T(Si5.98Al0.02)Σ6.00O18] (BO3)3 V(OH)3 W[F0.59(OH)0.48]Σ1.07
Brownish yellow core zone (M1CG):
X(Na0.59Ca0.310.10)Σ1.00 Y(Al1.16Li0.93Mn0.87Ti0.03Fe2+0.01)Σ3.00 ZAl6 [T(Si5.99Al0.01)Σ6.00O18] (BO3)3 V[(OH)2.92O0.08]Σ3.00 W(F0.60O0.40)Σ1.00
Note that the occurrence of solely Al in the Z sites is also supported by <Z-O> distances (1.905–1.907 Å), typical of <ZAl-O> in tourmaline [11].

4. Discussion

The detailed structural and chemical analyses performed on the MADAG1 slice allow us to characterize the optically observed color zoning and show a coexistence of two tourmaline species in the same crystal. The composition of the triangular brownish yellow core (M1CG) is characterized by predominantly XNa > XCa (Figure 11), WF > WO and (Li + Al) at the Y position in the general tourmaline formula, corresponding to a Mn-rich fluor-elbaite, ideally Na(Li1.5Al1.5)Al6(Si6O18)(BO3)3(OH)3F [11].
The outer zones (M1RVS, M1RVC and M1CR) are predominantly XCa > XNa (Figure 8), WF > WO and Y(Li + Al), corresponding to Mn-rich fluor-liddicoatite, Ca(Li2Al)Al6(Si6O18)(BO3)3(OH)3F [11].
For the studied samples, the compositional variation may be explained by the substitution XNa + 3Y(Mn2+ + Fe2+) ⇿ XCa + 2YLi + YAl. The plot of (Na + Mn + Fe) vs. (Ca + Li + Al) shows a negative correlation (r2 = 0.995) described by the equation: (Na + Mn + Fe) = 3.9 − (Ca + Li + Al).
According to Henry and Dutrow [28], tourmaline composition and zoning are generally a result of external influences (temperature, pressure, mineral local assemblages and fluid composition) and internal factors (crystallographic constraints).
The sharp oscillatory zoning observed in the liddicoatitic zone reflects fluctuating external conditions, as well as temperature, recorded as oscillatory compositional variations in adjacent growth stages during the crystallization.
The chemical/color zoning observed in our sample of tourmaline from Anjanabonoina provides useful information about the evolution of the growth environment. As liddicoatitic-elbaitic crystals are the final tourmaline to form, following the crystallization of Ca-bearing schorl and Mn-rich elbaite, they represent geochemically the most evolved tourmaline in the pegmatites [21]. The color zoning observed in the studied sample reflects specific initial conditions for the formation of Mn-rich yellow tourmalines, including a high concentration of Mn, B and Ti in appropriate proportions as suggested by Simmons et al. [29]. Additionally, the studied sample shows that throughout the rim, the chemical composition is totally different from that of the core, especially for the significant enrichment of Ca observed in the outer zones. This growth stage can be related to the final stage of solidification of pegmatite, when the fluid pressure inside the pocket causes rupture, promoting access to external fluids. The host rock of the Itremo group (1.1) adjacent to the hybrid Anjanabonoina pegmatite is interested by hydrothermal alteration, so the Ca necessary to form fluor-liddicoatite comes from the host rocks of this group which contain a unit consisting of marble [30]. The presence of these two coexisting species (fluor-liddicoatite and fluor-elbaite) in the same single crystal proves a melt contamination of this fluid derived from the metasediments of the Itremo group. During the ultimate phase of solidification, the fluid pressure in this cavity could induce rupture, enabling the influx of other elements (such as Ca or Fe) from external sources. These elements donate polychromatic terminations in this particular variety of tourmaline [22].
This conclusion agrees with Lussier et al. [31] who studied an oscillatory zoned liddicoatite tourmaline from Anjanabonoina pegmatite with chemical variation due to more than one input of external fluid contamination. Tindle et al. [32] also explained that the presence of Ca and Fe in the late stage is as a result of low temperature hydrothermal fluids that have reacted with host rocks rich in Ca and ferromagnesian minerals. Tourmaline crystals with a fluor-liddicoatite–fluor-elbaite composition with no potential secondary source of Ca have also been studied by several authors who proposed that Ca can remain in the melt until late fractionation processes [33,34].
All this evidence leads to the conclusion that each step of the growth of the studied watermelon Malagasy tourmaline is deeply correlated with the change in the geological environment, making this sample a brilliant example of the ability of tourmaline to record petrogenetic information of the crystallizing melt.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cryst13081290/s1, File CIF m1cg_a; File m1cr_a; File CIF m1rvc_a; File CIF m1rvs_a.

Author Contributions

Conceptualization, F.R. and G.A.; Methodology, F.B., G.T. and G.A.; Formal analysis, F.R., F.B. and G.T.; Writing—original draft, F.R.; Writing—review & editing, G.A. and F.B.; Supervision, G.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are grateful to Marcello Serracino and Nicola Mongelli for the essential assistance during the EMPA and SEM-EDS chemical analyses, respectively. The authors also thanks Giuseppe Secondo Elettivo for providing the sample.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Henry, D.J.; Dutrow, B.L. Tourmaline crystallography, crystal chemistry and nomenclature: Current status. Natura 2021, 111, 47–48. [Google Scholar]
  2. Hawthorne, F.C.; Henry, D.J. Classification of the Minerals of the Tourmaline Group. Eur. J. Mineral. Czech Republic 1999, 11, 201–215. [Google Scholar] [CrossRef]
  3. Moore, J.N.; Christenson, B.W.; Allis, R.G.; Browne, P.R.; Lutz, S.J. The mineralogical consequences and behavior of descending acid-sulfate waters: An example from the Karaha–Telaga Bodas geothermal system, Indonesia. Can. Mineral. 2004, 42, 1483–1499. [Google Scholar] [CrossRef]
  4. Manning, D.A.C.; Pichavant, M. The role of fluorine and boron in the generation of granitic melts. In High Grade Metamorphism, Migmatites and Melting, Meeting of the Geochemical Group of the Mineralogical Society; Shiva Publ: Nantwich, UK, 1983; pp. 94–109. [Google Scholar]
  5. Marschall, H.R.; Meyer, C.; Wunder, B.; Ludwig, T.; Heinrich, W. Experimental boron isotope fractionation between tourmaline and fluid: Confirmation from in situ analyses by secondary ion mass spectrometry and from Rayleigh fractionation modelling. Contrib. Mineral. Petrol. 2009, 158, 675–681. [Google Scholar] [CrossRef]
  6. van Hinsberg, V.J.; Henry, D.J.; Dutrow, B.L. Tourmaline as a petrologic forensic mineral: A unique recorder of its geologic past. Elements 2011, 7, 327–332. [Google Scholar] [CrossRef]
  7. Ballirano, P.; Celata, B.; Bosi, F. In situ high-temperature behavior and breakdown conditions of uvite at room pressure. Phys. Chem. Miner. 2022, 49, 40. [Google Scholar] [CrossRef]
  8. Ballirano, P.; Celata, B.; Bosi, F.; Beckett-Brown, C.E.; Andreozzi, G.B. Thermal behavior of schorl up to breakdown temperature at room pressure. Period. Di Mineral. 2023, 92, 23–32. [Google Scholar]
  9. Celata, B.; Stagno, V.; Capizzi, L.S.; Bosi, F.; Ballirano, P.; D’Arco, A.; Stopponi, V.; Lupi, S.; Scarlato, P.; Skogby, H.; et al. Schorl breakdown at upper mantle conditions: Insights from an experimental study at 3.5 GPa. Lithos 2023, 438, 106999. [Google Scholar] [CrossRef]
  10. Bosi, F. Tourmaline crystal chemistry. Am. Mineral. J. Earth Planet. Mater. 2018, 103, 298–306. [Google Scholar] [CrossRef]
  11. Agrosì, G.; Bosi, F.; Lucchesi, S.; Melchiorre, G.; Scandale, E. Mn-tourmaline crystals from island of Elba (Italy): Growth history and growth marks. Am. Mineral. 2006, 91, 944–952. [Google Scholar] [CrossRef]
  12. Johnson, M.L.; Wentzell, C.Y.; Elen, S. Multicolored bismuth-bearing tourmaline from Lundazi, Zambia. Gems Gemol. 1997, 33, 204–211. [Google Scholar] [CrossRef]
  13. Roda-Robles, E.; Simmons, W.; Pesquera, A.; Gil-Crespo, P.P.; Nizamoff, J.; Torres-Ruiz, J. Tourmaline as a petrogenetic monitor of the origin and evolution of the Berry-Havey pegmatite (Maine, USA). Am. Mineral. 2015, 100, 95–109. [Google Scholar] [CrossRef]
  14. Natkaniec-Nowak, L.; Dumanska-Slowik, M.; Ertl, A. “ Watermelon” tourmaline from the Paprok mine (Nuristan, Afghanistan). Neues Jahrb. Fur Mineral. Abh. 2009, 186, 185. [Google Scholar] [CrossRef]
  15. Wesełucha-Birczyńska, A.; Natkaniec-Nowak, L. A Raman microspectroscopic study of organic inclusions in “watermelon” tourmaline from the Paprok mine (Nuristan, Afghanistan). Vib. Spectrosc. 2011, 57, 248–253. [Google Scholar] [CrossRef]
  16. Cook, R.B. Connoisseur’s choice: Liddicoatite, fluor-liddicoatite, and liddicoatitic tourmalines, Anjanabonoina district, Madagascar. Rocks Miner. 2013, 88, 346–353. [Google Scholar] [CrossRef]
  17. Laurs, B.M. Tourmaline slices from Myanmar. J. Gemmol. 2015, 34, 668–670. [Google Scholar]
  18. Akizuki, M.; Kuribayashi, T.; Nagase, T.; Kitakaze, A. Triclinic liddicoatite and elbaite in growth sectors of tourmaline from Madagascar. Am. Mineral. 2001, 86, 364–369. [Google Scholar] [CrossRef]
  19. Dirlam, D.M.; Laurs, B.M.; Pezzotta, F.; Simmons, W.B. Liddicoatite tourmaline from Anjanabonoina, Madagascar. Gems Gemol. 2002, 38, 28–53. [Google Scholar] [CrossRef]
  20. Paquette, J.L.; Nédélec, A. A new insight into Pan-African tectonics in the East–West Gondwana collision zone by U–Pb zircon dating of granites from central Madagascar. Earth Planet. Sci. Lett. 1998, 155, 45–56. [Google Scholar] [CrossRef]
  21. Simmons, W.B.; Groat, L.A. Gem-bearing pegmatites. Geol. Gem. Depos. 2007, 37, 169–206. [Google Scholar]
  22. Collins, A.S. Madagascar and the amalgamation of Central Gondwana. Gondwana Res. 2006, 9, 3–16. [Google Scholar] [CrossRef]
  23. Raharimahefa, T.; Kusky, T.M. Structural and remote sensing studies of the southern Betsimisaraka Suture, Madagascar. Gondwana Res. 2006, 10, 186–197. [Google Scholar] [CrossRef]
  24. Pouchou, J.L.; Pichoir, F. Quantitative analysis of homogeneous or stratified microvolumes applying the model “PAP”. In Electron Probe Quantitation; Springer: Boston, MA, USA, 1991; pp. 31–75. [Google Scholar]
  25. Bosi, F.; Skogby, H.; Agrosì, G.; Scandale, E. Tsilaisite, NaMn3Al6(Si6O18)(BO3)3(OH)3OH, a new mineral species of the tourmaline supergroup from Grotta d’Oggi, San Pietro in Campo, island of Elba, Italy. Am. Mineral. 2012, 97, 989–994. [Google Scholar] [CrossRef]
  26. Sheldrick, G.M. Crystal structure refinement with SHELXL. Acta Crystallogr. Sect. C Struct. Chem. 2015, 71, 3–8. [Google Scholar] [CrossRef] [PubMed]
  27. Bosi, F.; Celata, B.; Skogby, H.; Hålenius, U.; Tempesta, G.; Ciriotti, M.E.; Bittarello, E.; Marengo, A. Mn-bearing purplish-red tourmaline from the Anjanabonoina pegmatite, Madagascar. Mineral. Mag. 2021, 85, 242–253. [Google Scholar] [CrossRef]
  28. Henry, D.J.; Novák, M.; Hawthorne, F.C.; Ertl, A.; Dutrow, B.L.; Uher, P.; Pezzotta, F. Nomenclature of the tourmaline-supergroup minerals. Am. Mineral. 2011, 96, 895–913. [Google Scholar] [CrossRef]
  29. Simmons, W.B.; Falster, A.U.; Laurs, B.M. A survey of Mn-rich yellow tourmaline from worldwide localities and implications for the petrogenesis of granitic pegmatites. Can. Mineral. 2011, 49, 301–319. [Google Scholar] [CrossRef]
  30. Pezzotta, F.; Praszkier, T. Giant liddicoatites from Madagascar. Minerals 2013, 6, 13–17. [Google Scholar]
  31. Lussier, A.J.; Hawthorne, F.C. Oscillatory zoned liddicoatite from Anjanabonoina, central Madagascar. II. Compositional variation and mechanisms of substitution. Can. Mineral. 2011, 49, 89–104. [Google Scholar] [CrossRef]
  32. Tindle, A.G.; Selway, J.B.; Breaks, F.W. Liddicoatite and associated species from the McCombe spodumene-subtype rare-element granitic pegmatite, northwestern Ontario, Canada. Can. Mineral. 2005, 43, 769–793. [Google Scholar] [CrossRef]
  33. Novak, M.; Selway, J.B.; Cerny, P.; Hawthorne, F.C.; Ottolini, L. Tourmaline of the elbaite-dravite series from an elbaite-subtype pegmatite at Blizna, southern Bohemia, Czech Republic. Eur. J. Mineral. 1999, 11, 557–568. [Google Scholar] [CrossRef]
  34. Selway, J.B.; Smeds, S.A.; Černý, P.; Hawthorne, F.C. Compositional evolution of tourmaline in the petalite-subtype Nyköpingsgruvan pegmatites, Utö, Stockholm Archipelago, Sweden. GFF 2002, 124, 93–102. [Google Scholar] [CrossRef]
Figure 1. Optical micrographs of the tourmaline sample from the Anjanabonoina pegmatite (Madagascar): (a) raw crystal with the red lines corresponding to slice cut perpendicular to the c axis; (b) basal polished slice, labelled MADAG1, showing color zoning typical of the “watermelon” variety.
Figure 1. Optical micrographs of the tourmaline sample from the Anjanabonoina pegmatite (Madagascar): (a) raw crystal with the red lines corresponding to slice cut perpendicular to the c axis; (b) basal polished slice, labelled MADAG1, showing color zoning typical of the “watermelon” variety.
Crystals 13 01290 g001
Figure 2. Simplified geology map of Madagascar (modified from [22,23]).
Figure 2. Simplified geology map of Madagascar (modified from [22,23]).
Crystals 13 01290 g002
Figure 3. Basal slice of MADAG1 (perpendicular to the c axis): optical image (left) and schematic sketch (right). Note the four main homogeneous colored zones. The red arrows indicate the EPMA traverses performed across the rim (A-B), the pink zone (C-D) and the core (E-F). Light yellow represents the transition zone around the M1CG.
Figure 3. Basal slice of MADAG1 (perpendicular to the c axis): optical image (left) and schematic sketch (right). Note the four main homogeneous colored zones. The red arrows indicate the EPMA traverses performed across the rim (A-B), the pink zone (C-D) and the core (E-F). Light yellow represents the transition zone around the M1CG.
Crystals 13 01290 g003
Figure 4. Outer zone of the MADAG1 slice: (a) optical image; (b) VPSEM image showing the complex oscillatory zoning in the pink region.
Figure 4. Outer zone of the MADAG1 slice: (a) optical image; (b) VPSEM image showing the complex oscillatory zoning in the pink region.
Crystals 13 01290 g004
Figure 5. EPMA traverse across the outer zone of slice (length 724 µm) with 60 analysis points and a step size of 12 µm with the variation in the main elements (wt %).
Figure 5. EPMA traverse across the outer zone of slice (length 724 µm) with 60 analysis points and a step size of 12 µm with the variation in the main elements (wt %).
Crystals 13 01290 g005
Figure 6. EPMA traverse across the pink rim zone of slice (length 2850 µm) with 12 analysis points and a step size of 260 µm with the variation in the main elements (wt %).
Figure 6. EPMA traverse across the pink rim zone of slice (length 2850 µm) with 12 analysis points and a step size of 260 µm with the variation in the main elements (wt %).
Crystals 13 01290 g006
Figure 7. EPMA traverse at the interface between the M1CR and the core zone of slice (length 4230 µm) with 20 analysis points and a step size of 211 µm with the variation in the main elements (oxide wt %). Note the sharp variation in MnO in the light yellow transition zone that surrounds the core.
Figure 7. EPMA traverse at the interface between the M1CR and the core zone of slice (length 4230 µm) with 20 analysis points and a step size of 211 µm with the variation in the main elements (oxide wt %). Note the sharp variation in MnO in the light yellow transition zone that surrounds the core.
Crystals 13 01290 g007
Figure 8. Calibration of the Li emission line at 670.7 nm counts vs. Li2O% of the standards. The blue circles represent the standards; The colored circles represent the four different zones of the MADAG1 tourmaline sample.
Figure 8. Calibration of the Li emission line at 670.7 nm counts vs. Li2O% of the standards. The blue circles represent the standards; The colored circles represent the four different zones of the MADAG1 tourmaline sample.
Crystals 13 01290 g008
Figure 9. Basal slice of MADAG1; the red circles show the zones of the slice where fragments were extracted for XRD analysis.
Figure 9. Basal slice of MADAG1; the red circles show the zones of the slice where fragments were extracted for XRD analysis.
Crystals 13 01290 g009
Figure 10. Correlation between the unit cell volume (Å3) and Mn + Fe content (apfu) in the present samples. The blue dashed line represents regression line.
Figure 10. Correlation between the unit cell volume (Å3) and Mn + Fe content (apfu) in the present samples. The blue dashed line represents regression line.
Crystals 13 01290 g010
Figure 11. Ternary system for the primary tourmaline groups based on the dominant occupancy of the X site. The colored circles represent the average chemical composition obtained with EPMA and SREF data for each zone (M1RVS, M1RVC, M1CR and M1CG).
Figure 11. Ternary system for the primary tourmaline groups based on the dominant occupancy of the X site. The colored circles represent the average chemical composition obtained with EPMA and SREF data for each zone (M1RVS, M1RVC, M1CR and M1CG).
Crystals 13 01290 g011
Table 1. Chemical composition obtained with EPMA analyses of MADAG1 slice for the four selected zones (M1RVS, M1RVC, M1CR and M1CG).
Table 1. Chemical composition obtained with EPMA analyses of MADAG1 slice for the four selected zones (M1RVS, M1RVC, M1CR and M1CG).
Average of 14 Spots (wt %)Average of 15 Spots (wt %)Average of 20 Spots (wt %)Average of 10 Spots (wt %)
M1RVSM1RVCM1CRM1CG
SiO238.10 (17)38.44 (23)38.58 (16)37.22 (21)
TiO20.18 (6)0.15 (4)0.09 (2)0.27 (4)
B2O3 (a)11.0111.0911.2110.81
Al2O337.82 (23)38.12 (22)39.43 (8)37.84 (24)
FeO tot0.95 (20)0.45 (17)-0.03 (3)
MnO tot3.53 (50)1.64 (31)1.16 (26)6.42 (12)
CaO2.99 (22)3.80 (31)3.16 (34)1.81 (8)
Na2O1.35 (13)0.99 (6)1.16 (12)1.90 (4)
Li2O (b)2.19 (10)2.60 (6)2.65 (37)1.44 (9)
F1.18 (16)1.19 (18)1.20 (15)1.17 (17)
H2O (a)3.193.203.372.72
O = F−0.50−0.50−0.50−0.49
Total102.03101.17101.49101.13
Atoms normalized to 31 anions
Si (apfu)5.996.025.985.98
Ti4+0.0210.020.010.03
B3.003.003.003.00
Al7.017.047.217.17
Fe2+0.130.060.000.00
Mn2+0.470.220.150.87
Ca0.500.640.520.31
Na0.410.300.350.59
Li1.391.641.650.93
F0.590.590.590.59
OH3.353.353.482.92
(a) Calculated by stoichiometry and SREF information. (b) Determined by μ-LIBS. Errors for oxides and fluorine are standard deviations (in brackets).
Table 2. Single-crystal X-ray diffraction data details for MADAG1.
Table 2. Single-crystal X-ray diffraction data details for MADAG1.
M1RVSM1RVCM1CRM1CG
Crystal data
Size (mm)0.30 × 0.30 × 0.300.29 × 0.35 × 0.350.15 × 0.20 × 0.200.25 × 0.27 × 0.35
Space group; ZR3m, 3
a (Å)15.8740 (3)15.8345 (2)15.8313 (2)15.8839 (2)
c (Å)7.1138 (1)7.1039 (1)7.1040 (1)7.1160 (1)
V3)1552.40 (6)1542.54 (4)1541.94 (4)1554.82 (4)
Data collection
Data collection temperature (K)293
Range for data collection, 2θ (°)75
Radiation, wavelength (Å)MoKα, 0.71073
Reciprocal space range, hkl−26 ≤ h ≤ 24−26 ≤ h ≤ 26−26 ≤ h ≤ 25−24 ≤ h ≤ 21
−26 ≤ k ≤ 27−26 ≤ k ≤ 25−26 ≤ k ≤ 27−26 ≤ k ≤ 27
−11 ≤ l ≤ 11−11 ≤ l ≤ 11−11 ≤ l ≤ 10−11 ≤ l ≤ 12
Measured reflections11,55811,47211,59811,591
Unique reflections1876185117331891
Redundancy12
Absorption correction methodSADABS
Refinement
Refinement methodFull matrixleast squares on F2
Structural refinement programSHELXL-2013
Extinction coefficient0.0023 (3)0.0055 (3)0.0020 (3)0.0054 (4)
Flack parameter0.23 (3)0.25 (4)0.30 (4)0.23 (3)
wR20.04210.04030.04090.0437
R1 all data0.01680.01620.01700.0176
R1 for I > 2σI0.01680.01610.01660.0175
GooF1.1241.1291.1211.129
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rizzo, F.; Bosi, F.; Tempesta, G.; Agrosì, G. Compositional Variation and Crystal-Chemical Characterization of a Watermelon Variety of Tourmaline from Anjanabonoina, Central Madagascar. Crystals 2023, 13, 1290. https://doi.org/10.3390/cryst13081290

AMA Style

Rizzo F, Bosi F, Tempesta G, Agrosì G. Compositional Variation and Crystal-Chemical Characterization of a Watermelon Variety of Tourmaline from Anjanabonoina, Central Madagascar. Crystals. 2023; 13(8):1290. https://doi.org/10.3390/cryst13081290

Chicago/Turabian Style

Rizzo, Floriana, Ferdinando Bosi, Gioacchino Tempesta, and Giovanna Agrosì. 2023. "Compositional Variation and Crystal-Chemical Characterization of a Watermelon Variety of Tourmaline from Anjanabonoina, Central Madagascar" Crystals 13, no. 8: 1290. https://doi.org/10.3390/cryst13081290

APA Style

Rizzo, F., Bosi, F., Tempesta, G., & Agrosì, G. (2023). Compositional Variation and Crystal-Chemical Characterization of a Watermelon Variety of Tourmaline from Anjanabonoina, Central Madagascar. Crystals, 13(8), 1290. https://doi.org/10.3390/cryst13081290

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop