Copyright Statement

Abstract

Rocks in metalliferous ore deposits are typically classified as either ore or waste rock. An alternative distinction, more helpful for discriminating alteration haloes around ore, is protolith rock and altered rock. Metasomatic alteration occurs during hydrothermal ore deposition when mineralising fluids are in disequilibrium with surrounding host rocks. Geologists in brownfields areas commonly seek to interpret ore body geometry and surrounding alteration halos using the spatial distribution of whole rock chemical analyses. However, modelling the relative highs and lows in geochemical data does not give a complete representation of alteration geometry; weight percent values provide no direct information about the nature and magnitude of mass transfer during hydrothermal alteration. The compositional data of rocks are inherently multivariate and the concentrations of all constituents are related to each other by the process of closure. For example, in a simplified three element system, reducing one element will lead to an apparent enrichment of the other two.

Here we present a workflow to process a whole rock geochemical dataset and produce labels for groups of similar least-altered protolith rocks, and then apply these labels to groups of altered rocks. This method involves quick means of exploratory data analysis and lends itself to iterative refinement, relying on a geologist to use their knowledge of a particular site.