Crystals from the Powellite-Scheelite Series at the Nanoscale: A Case Study from the Zhibula Cu Skarn, Gangdese Belt, Tibet

Abstract

Scheelite (CaWO4) and powellite (CaMoO4) are isostructural minerals considered as a non-ideal solid solution series. Micron- to nanoscale investigation of a specimen of skarnoid from Zhibula, Gangdese Belt, Tibet, China, was carried out to assess the identity of the phases within a broad scheelite-powellite (Sch-Pow) compositional range, and to place additional constraints on redox changes during ore formation. An electron probe microanalysis shows that Mo-rich domains within complex oscillatory-zoned single crystals, and as thin sliver-like domains, have a compositional range from 20 mol.% to 80 mol.% Pow. These occur within a matrix of unzoned, close-to-end-member scheelite aggregates (87 mol.%–95 mol.% Sch). Laser-ablation inductively coupled plasma mass spectrometry spot analysis and element mapping reveal systematic partitioning behaviour of trace elements in skarn minerals (grossular50, diopside80, anorthite, and retrograde clinozoisite) and scheelite-powellite aggregates. The Mo-rich domains feature higher concentrations of As, Nb, and light rare earth elements LREE, whereas W-rich domains are comparatively enriched in Y and Sr. Transmission electron microscopy (TEM) was carried out on focused-ion-beam-prepared foils extracted in situ from domains with oscillatory zoning occurring as slivers of 20 mol.%–40 mol.% Pow and 48 mol.%–80 mol.% Pow composition within an unzoned low-Mo matrix (20 mol.% Pow). Electron diffractions, high-angular annular dark field (HAADF) scanning-TEM (STEM) imaging, and energy-dispersive spectroscopy STEM mapping show chemical oscillatory zoning with interfaces that have continuity in crystal orientation throughout each defined structure, zoned grain or sliver. Non-linear thermodynamics likely govern the patterning and presence of compositionally and texturally distinct domains, in agreement with a non-ideal solid solution. We show that the sharpest compositional contrasts are also recognisable by variation in growth direction. Atomic-scale resolution imaging and STEM simulation confirm the presence of scheelite-powellite within the analysed range (20 mol.%–80 mol.% Pow). Xenotime-(Y) inclusions occur as nm-wide needles with epitaxial orientation to the host scheelite-powellite matrix throughout both types of patterns, but no discrete Mo- or W-bearing inclusions are observed. The observed geochemical and petrographic features can be reconciled with a redox model involving prograde deposition of a scheelite+molybdenite assemblage (reduced), followed by interaction with low-T fluids, leading to molybdenite dissolution and reprecipitation of Mo as powellite-rich domains (retrograde stage, oxidised). The observation of nanoscale inclusions of xenotime-(Y) within scheelite carries implications for the meaningful interpretation of petrogenesis based on rare earth element (REE) concentrations and fractionation patterns. This research demonstrates that HAADF-STEM is a versatile technique to address issues of solid solution and compositional heterogeneity.