Elsevier

Computers & Geosciences

Volume 89, April 2016, Pages 32-43
Computers & Geosciences

Research paper
Precursors predicted by artificial neural networks for mass balance calculations: Quantifying hydrothermal alteration in volcanic rocks

https://doi.org/10.1016/j.cageo.2016.01.003Get rights and content

Highlights

  • Method for predicting the fresh precursor to hydrothermally altered volcanic rocks.

  • Precursors are predicted using multi-layer perceptron neural networks.

  • Precursor compositions (major oxides) are predicted from immobile elements ratios.

  • Chemical method for quantifying alteration from large datasets of volcanic rocks.

  • New chemical method developed for grassroot or regional mineral exploration.

Abstract

This study proposes an artificial neural networks-based method for predicting the unaltered (precursor) chemical compositions of hydrothermally altered volcanic rock. The method aims at predicting precursor’s major components contents (SiO2, FeOT, MgO, CaO, Na2O, and K2O). The prediction is based on ratios of elements generally immobile during alteration processes; i.e. Zr, TiO2, Al2O3, Y, Nb, Th, and Cr, which are provided as inputs to the neural networks. Multi-layer perceptron neural networks were trained on a large dataset of least-altered volcanic rock samples that document a wide range of volcanic rock types, tectonic settings and ages. The precursors thus predicted are then used to perform mass balance calculations. Various statistics were calculated to validate the predictions of precursors’ major components, which indicate that, overall, the predictions are precise and accurate. For example, rank-based correlation coefficients were calculated to compare predicted and analysed values from a least-altered test dataset that had not been used to train the networks. Coefficients over 0.87 were obtained for all components, except for Na2O (0.77), indicating that predictions for alkali might be less performant. Also, predictions are performant for most volcanic rock compositions, except for ultra-K rocks. The proposed method provides an easy and rapid solution to the often difficult task of determining appropriate volcanic precursor compositions to rocks modified by hydrothermal alteration. It is intended for large volcanic rock databases and is most useful, for example, to mineral exploration performed in complex or poorly known volcanic settings. The method is implemented as a simple C++ console program.

Introduction

Grassroots, or regional mineral exploration, necessitates the use of powerful and simple methods to interpret large amounts of chemical data. The interpretation and quantification of hydrothermal alterations in volcanic terrains, for instance, is particularly important as many economic substances are accumulated by the circulation of hydrothermal fluids in areas with magmatic activity, and because altered rocks are one of the main vectors used to explore for such mineralisations.

Geochemical studies of hydrothermally altered rocks aim at determining the amount of elements gained and lost during alteration of an initially fresh rock (i.e. the precursor), by performing mass balance calculations while assuming that one or more elements is immobile during alteration (see Gresens (1967)). Gresens' (1967) approach has been extensively and successfully used (MacLean and Kranidiotis, 1987, Shriver and MacLean, 1993, Cadieux et al., 2006; to name a few).

A major question that arises when using such mass balance equations concerns the choice of an appropriate precursor for the studied altered sample. Precursor’s compositions are usually obtained by analysing “fresh” samples (i.e. samples that lack alteration minerals) that have been geologically related to the altered sample using detailed knowledge of the local geology (e.g. Grant, 1986, Grant, 2005). Such approaches are possible only for detailed deposit-scale studies, if fresh rocks exist and if our level of knowledge of the area render their identification possible. Even so, usually nothing proves that the fresh sample is an exact match for the precursor to the altered sample, and these difficulties can make mass balance calculations un-reliable.

At a more regional scale, on the other hand, these difficulties can be in-surmountable. Thus, orebody targeting using chemical datasets more often rely on alteration indexes (ratio of major elements), which values are sensitive to the precursor’s compositions and could provide miss-leading indications on the intensity of alteration (see discussion in Trépanier et al. (2015)).

Mass balance calculations would be more reliable if a suitable precursor composition to an altered rock could be determined directly from the chemical characteristics of the altered sample, thus providing a precursor adapted to this sample.

This paper proposes a method for estimating the composition in SiO2, CaO, MgO, FeOT, K2O, and Na2O of fresh precursors to altered volcanic rocks, using ratios of commonly analysed least mobile elements (i.e. Zr, TiO2, Al2O3, Y, Nb, Th, and Cr) and a regression method (i.e. artificial neural network). The method is provided as supplemental material, as a simple C++ program that calculates precursor compositions (this study) prior performing mass balance calculations (using published methods).

Section snippets

Description of the method

In this contribution, a method to predict fresh precursors from the chemical composition (i.e. immobile elements content) of altered volcanic rocks is presented. Predictions are performed for magmatic rocks only because relatively simple relationships, which are controlled by magmatic processes, exist between the trace and major elements content of such rocks. For example, Zr and Si are both incompatible elements that tend to increase during fractional crystallisation, making Zr a proxy for Si

Least-altered volcanic rocks dataset

Relationships between major and immobile elements are evaluated here with neural networks using a large chemical dataset of modern volcanic rocks from around the World. Post-Archean to modern rocks are used because they are well constrained samples for which names, geodynamic settings and freshness can be confidently established; i.e. characteristics that would be harder to determine for more ancient samples. Besides, there might be large differences between ancient (Archean) and modern

Software

This study aims at modelling the chemical composition of fresh precursors to altered rocks, in order to facilitate mass balance calculations. The outputs of the neural networks (i.e. the results of the training process) can however be difficult to use as such, as the software used to train the neural networks (NeuroDimension, 2012) provides outputs files with a “.NSW” extension, which are simple ASCII text files that record the network’s structure, parameters and weights.

To facilitate access to

Neural network validations

Prior to using the method to predict new sample’s precursors, the neural networks’ outputs were validated. The consequences that Zr, Y, Nb, Cr, Th, Al, and Ti mobility might have on the prediction’s accuracy are also tested. Overall, the proposed approach appears accurate, precise, and may tolerate a limited mobility of the immobile elements used by the model.

Method validation using case studies

In this section, several fresh volcanic rocks, which were not included in the Georoc database used to train the neural network, are used to validate the approach (see case studies 1–3). Then, altered rocks from several mineralised areas (i.e. volcanogenic massive sulphide, or VMS, and porphyry deposits) are compiled, their precursors compositions are predicted and mass balance calculations are performed (see case studies 4–6).

Discussion and conclusions

This study proposes a method for predicting the composition in major elements of precursors to hydrothermally altered volcanic rocks, using artificial neural networks. To facilitate its utilisation by the reader, the method is embedded in a C++ console program that also performs mass balance calculations, allowing for the evaluation of alteration, in terms of type and magnitude, from large datasets of volcanic rocks.

Neural networks used here “learned” to predict the composition in major

Computer code

The method presented in this contribution is provided as a computer code (C++ language).

Acknowledgements

The authors wish to address special thanks to the editor and reviewers of the manuscript, Jef Caers and David Lentz, as well as an anonymous reviewer. This project was supported by Canada Economic Development, the Ministère de l'Énergie et des Ressources Naturelles du Québec, the Conférence Régionale des élus Saguenay-Lac-Saint-Jean and companies members of the Consorem. The authors warmly thank their colleague Silvain Rafini for constructive discussions on this project, Geneviève Boudrias for

References (52)

  • J.A. Winchester et al.

    Geochemical discrimination of magma series and their differentiation products using immobile elements

    Chem. Geol.

    (1977)
  • D.A. Wood

    The application of a Th-Hf-Ta diagram to problems of tectonomagmatic classification and to establishing the nature of crustal contamination of basaltic lavas of the British Tertiary Volcanic Province

    Earth Planet. Sci. Lett.

    (1980)
  • T.J. Barrett et al.

    Chemostratigraphy and hydrothermal alteration in exploration for VHMS deposits in greenstones and younger volcanic rocks

  • E. Boari et al.

    Transition from ultrapotassic kamafugitic to sub-alkaline magmas: Sr, Nd and Pb isotope, trace element and 40Ar–39Ar Age data from the Middle Latin Valley Volcanic Field, Roman Magmatic Province, Central Italy

    J. Petrol.

    (2009)
  • Cadieux, A.-M., Dubé, B., Williamson, K., Malo, M., Twomey, T., 2006. Characterization of hydrothermal alterations at...
  • J.F. Cohen

    Mineralogy and Geochemistry of Hydrothermal Alteration at the Ann- Mason Porphyry Copper Deposit, Nevada: Comparison of Large-scale Ore Exploration Techniques to Mineral Chemistry (Unpublished Master Thesis)

    (2011)
  • Cohen, J.F., Dilles, J.H., Tosdal, R.M., Halley, S., 2011. Compositional variations in hydrothermal white micas and...
  • J.H. Dilles

    Footprints of Porphyry Cu Deposits: Vectors to the Hydrothermal Center Using Mineral Mapping and Lithogeochemistry

    (2012)
  • Drummond, M.S., Defant, M.J., 1990. A model for trondhjemite‐tonalite‐dacite genesis and crustal growth via slab...
  • T. Finlow-Bates et al.

    The behaviour of so-called immobile elements in hydrothermally altered rocks associated with volcanogenic submarine exhalative ore deposits

    Miner. Deposita

    (1981)
  • GEOROC – Geochemistry of Rocks of the Oceans and Continents, 2011. Max-Planck-Institut für Chemie....
  • J.A. Grant

    The Isocon Diagram – a simple solution to gresens’s rquation for metasomatic alteration

    Econ. Geol.

    (1986)
  • R.L. Gresens

    Composition–volume relationships of metasomatism

    Chem. Geol.

    (1967)
  • A.R. Hastie et al.

    Classification of altered volcanic island arc rocks using immobile trace elements: development of the Th–Co discrimination diagram

    J. Petrol.

    (2007)
  • I.G. Hill et al.

    Yttrium: the immobility–mobility transition during basaltic weathering

    Geology

    (2000)
  • K. Hornik et al.

    Multilayer feedforward networks are universal approximators

    Neural Netw.

    (1989)
  • Cited by (22)

    • Applying neural networks-based modelling to the prediction of mineralization: A case-study using the Western Australian Geochemistry (WACHEM) database

      2023, Ore Geology Reviews
      Citation Excerpt :

      Multiple attempts have been made to apply machine learning (ML) methods to mineral exploration (Caté et al., 2017; Enkhsaikhan et al., 2021; Grunsky and Arne, 2021; Kirkwood et al., 2016; Prado et al., 2020; Rodriguez-Galiano et al., 2015; Sun et al., 2019; Wang et al., 2020; Wedge et al., 2019; Zuo and Xiong, 2017). Some specific examples include prediction of the presence of gold mineralization in drill core from well logs (Caté et al., 2017), deposit estimation from multi-element geochemical data (Hill et al., 2014), prediction of soil geochemistry from high resolution geophysical and remote sensing data (Kirkwood et al., 2016), and estimation of precursor chemical compositions of hydrothermally altered volcanic rocks from the ratios of commonly analyzed least mobile elements (Trépanier et al., 2016). However, all these approaches utilize either shallow neural networks, or methods such as support vector machines (SVMs), regression trees (RTs), and random forests (RFs).

    • Restoring original composition of hydrothermally altered Archean metavolcanic rocks of the Carajás Mineral Province (Brazil): Geodynamic implications for the transition from lid to mobile tectonics

      2020, Lithos
      Citation Excerpt :

      Such a conventional approach has yet to be challenged by new major- and trace-element systematics developed to specifically discriminate Archean volcanic sites. Meanwhile, the ANN-based method of Trépanier et al. (2016) should be developed further in order to propose a full range of model results integrating minor- and trace elements for both predicted precursors and quantified hydrothermal alterations. Precursor contents in high field strength elements (HFSE) and rare earth elements (REE) could then be used to better characterize mantle sources and improve petrogenetic models for Archean metabasalts.

    • Use of semi-volatile metals as a new vectoring tool for VMS exploration: Example from the Zn-rich McLeod deposit, Abitibi, Canada

      2019, Journal of Geochemical Exploration
      Citation Excerpt :

      Trace elements include metals, metalloids, and incompatibles, such as high field strength elements (HFSE) and the large ion lithophile elements (LILE). The HFSE are commonly immobile and mostly used for quantifying the degree of mobility of lithophile elements (MacLean and Kranidiotis, 1987; Rubin et al., 1993; Trépanier et al., 2016). Semi-volatile metals are a suite of trace elements characterized by a low-melting-point and enhanced mobility in hydrothermal systems (fluid or vapor).

    View all citing articles on Scopus
    View full text