Misinterpretation of zircon ages in layered intrusions

A recent re-interpretation of the Bushveld Complex and other layered intrusions as stacks of randomly emplaced, amalgamated sills is mostly fuelled by finding of zircon ages that are not getting progressively younger from the base upwards, as expected from a classical model for the formation of layered intrusions. Rather, they display several reversals from older to younger ages and vice-versa with moving up-section through the layered intrusions. Here, we show that the reported zircon ages are at odds with the relative ages of rocks as defined by cross-cutting relations in potholes of the Bushveld Complex. This indicates that interpretation of the zircon isotopic data as the emplacement age of the studied rocks/units is incorrect, making a new emplacement model for layered intrusions baseless. This conclusion is further buttressed by the phase equilibria analysis showing that regular cumulate sequences of layered intrusions are not reconcilable with a model of randomly emplaced sills. In this model, the late sills are free to intrude at any stratigraphic position of the pre-existing rocks, producing magmatic bodies with chaotic crystallization sequences and mineral compositional trends that are never observed in layered intrusions. There are thus no valid justifications for the re-evaluation of the current petrological model of the Bushveld Complex and other layered intrusions as large, long-lived and largely molten magma chambers. A fundamental implication of this analysis is that the current high-precision U-Pb TIMS ages from layered intrusions are inherently unreliable on the scale of several million years and cannot therefore be used for rigorous estimations of the timing of crystallization, duration of magmatism, and cooling of these intrusions.