Progressive evolution of whole‐rock composition during metamorphism revealed by multivariate statistical analyses

Abstract

The geochemical evolution of metamorphic rocks during subduction-related metamorphism is described on the basis of multivariate statistical analyses. The studied data set comprises a series of mapped metamorphic rocks collected from the Sanbagawa metamorphic belt in central Shikoku, Japan, where metamorphic conditions range from the pumpellyite-actinolite to epidote-amphibolite facies. Recent progress in computational and information science provides a number of algorithms capable of revealing structures in large data sets. This study applies k-means cluster analysis (KCA) and non-negative matrix factorization (NMF) to a series of metapelites, which is the main lithotype of the Sanbagawa metamorphic belt. KCA describes the structures of the high-dimensional data, while NMF provides end-member decomposition which can be useful for evaluating the spatial distribution of continuous compositional trends. The analysed data set, derived from previously published work, contains 296 samples for which 14 elements (Si, Ti, Al, Fe, Mn, Mg, Ca, Na, K, P, Rb, Sr, Zr and Ba) have been analysed. The KCA and NMF analyses indicate five clusters and four end-members, respectively, successfully explaining compositional variations within the data set. KCA indicates that the chemical compositions of metapelite samples from the western (Besshi) part of the sampled area differ significantly from those in the east (Asemigawa). In the west, clusters show a good correlation with the metamorphic grade. With increasing metamorphic grade, there are decreases in SiO 2 and Na 2 O and increases in other components. However, the compositional change with metamorphic grade is less obvious in the eastern area. End-member decomposition using NMF revealed that the evolutional change of whole-rock composition, as correlated with metamorphic grade, approximates a stoichiometric increase of a garnet-like component in the whole-rock composition, possibly due to the precipitation of garnet and effusion of other components during progressive dehydration. Thermodynamic modelling of the evolution of the whole-rock composition yielded the following results: (1) the whole-rock composition at lower metamorphic grade favours the preferential crystallization of garnet under the conditions of the garnet zone, with biotite becoming stable together with garnet in higher-grade

Innovation adoption: Broadcasting versus virality

“…Correspondingly, researchers in other disciplines may read the original article or the citations of the algorithm, and computer scientists may collaborate with researchers from other disciplines and pass different algorithms to them. For example, in a recent study (Yoshida, Kuwatani, Hirajima, Iwamori, & Akaho, 2018) about the evolution of whole‐rock composition during metamorphism, LDA was introduced into geology and used to find endmembers based on the frequencies of elements that make up the rock of interest. Most of the authors, Yoshida, Kuwatani, Hirajima, and Iwamori, are scientists from geology or earth science, while one of them, Shotaro Akaho, is from computer science.…”

Asymptotic Bayesian Generalization Error in Latent Dirichlet Allocation and Stochastic Matrix Factorization

“…The topic model [16] is a ubiquitous learning machine used in many research areas, including text mining [9,17], computer vision [25], marketing research [40], and geology [51]. The topic model is also known as latent Dirichlet allocation (LDA) [9] in the Bayesian terminology.…”

“…Latent Dirichlet allocation (LDA) [7] is one of topic models [10] which is a ubiquitous statistical model used in many research areas. Text mining [7,11], computer vision [18], marketing research [24], and geology [34] are such examples. LDA had been originally proposed for natural language processing and it can extract an essential information from documents by defining the topics of the words.…”

GEOFCM: a new method for statistical classification of geochemical data using spatial contextual information

“…Multivariate statistical analysis can be a powerful tool for detecting the characteristic features of multidimensional geochemical datasets, and it is employed in areas of study including environmental sciences, petrolo-gy, and resource exploration (e.g., Templ et al, 2008;Kuwatani et al, 2014;Ofner et al, 2015;Yasukawa et al, 2016;Ueki and Iwamori, 2017;Kuwatani et al, 2018;Ueki et al, 2018;Yoshida et al, 2018). Geochemical data including chemical composition, petrographic description, and location coordinates are increasingly provided in the literature, making them available for analysis approaches that deal with large-scale, multidimensional data (Igarashi et al, 2016).…”

“…Most geologic samples have their sample location as potential information that is generally provided as two-dimensional longitude/latitude-coordinates. Yoshida et al (2018) performed geocoding of the sampling map from the literature and statistical analyses on the corresponding geochemical data, comparing the clustering results from the geochemical data and its areal distributions. Recently, Haraguchi et al (2017) performed geocoding of published geochemical data since the 1980s and constructed a large location-based geochemical database called 'DODAI'.…”