Organization of oscillatory zoning in zircon: analysis, scaling, geochemistry, and model of a zircon from Kipawa, Quebec, Canada

The character of oscillatory zoning within a zircon crystal from the syenite Kipawa Complex, Quebec, varies with scale of observation. Analysis of an scanning electron microscopy (SEM) back-scatter gray-scale traverse at a resolution of one pixel = 2.43 μm revealed 145 zones over 5130 μm, whereas a detailed high-resolution (one pixel = 0.195 μm) section near the crystal rim revealed 225 zones over 795 μm. In order to mathematically characterize the zoning pattern, wavelet, Fourier, and nonlinear analysis techniques were used on profiles of the SEM gray-scale data, and a series constructed was from the zone widths. Results demonstrate that the zircon oscillatory zoning preserves nonlinear and periodic components. Secondary ion mass spectrometry, electron microprobe, and SEM analyses of trace elements show the SEM back-scatter bright zones are enriched in U, Th, and rare earth elements (REE) in comparison to the darker zones. REE patterns are sharply heavy REE enriched and have negative Eu anomalies and prominent positive Ce anomalies. We model the oscillatory zoning, including a measure of its chemical variation, by use of a periodically forced nonlinear system. Results of this data-driven model are quantitatively similar to the natural data. We envisage that the small-scale oscillatory zoning was the result of a nonlinear feedback process wherein the crystal growth modified the adjacent melt, which in turn affected the crystal composition. The large-scale harmonic zones likely reflect changes in the bulk geochemistry of the system from which the zircon grew.