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

The geochemical evolution of metamorphic rocks during subduction-related metamorphism is described on the basis of multivariate statistical analyses. The studied dataset 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 datasets. 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 endmember decomposition which can be useful for evaluating the spatial distribution of continuous compositional trends. The analyzed dataset, 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 analyzed. The KCA and NMF analyses indicate five clusters and four endmembers, respectively, successfully explaining compositional variations within the dataset. 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 SiO2 and Na2O and increases in other components. On the other hand, the compositional change with metamorphic grade is less obvious in the eastern area. Endmember 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 rock compositions under the same P-T conditions; (2) with higher-grade whole–rock compositions, more H2O is retained. These results provide insight into the mechanism suppressing dehydration under high-pressure metamorphic conditions. This mechanism should be considered in forward modelling of the fluid cycle in subduction zones, although such a quantitative model has yet to be developed. This article is protected by copyright. All rights reserved.

[1]  T. Pettke,et al.  Trace element signature of subduction-zone fluids, melts and supercritical liquids at 120–180 km depth , 2005, Nature.

[2]  R. Powell,et al.  A thermodynamic model for Ca–Na clinoamphiboles in Na2O–CaO–FeO–MgO–Al2O3–SiO2–H2O–O for petrological calculations , 2005 .

[3]  T. Evans.  A method for calculating effective bulk composition modification due to crystal fractionation in garnet‐bearing schist: implications for isopleth thermobarometry , 2004 .

[4]  J. Hermann,et al.  Aqueous fluids and hydrous melts in high-pressure and ultra-high pressure rocks: Implications for element transfer in subduction zones , 2006 .

[5]  G. Wörner,et al.  Compositional variations of ignimbrite magmas in the Central Andes over the past 26 Ma — A multivariate statistical perspective , 2016 .

[6]  K. Kiminami,et al.  The parentage of low-grade metasediments in the Sanbagawa Metamorphic Belt, Shikoku, southwest Japan, based on whole-rock geochemistry , 2003 .

[7]  S. K. Verma,et al.  Multidimensional classification of magma types for altered igneous rocks and application to their tectonomagmatic discrimination and igneous provenance of siliciclastic sediments , 2017 .

[8]  James B. Thompson.  The Graphical Analysis of Mineral Assemblages in Pelitic Schists , 1957 .

[9]  T. Higashino.  The higher grade metamorphic zonation of the Sambagawa metamorphic belt in central Shikoku, Japan , 1990 .

[10]  E. Ghent,et al.  ESTIMATING P-T CONDITIONS OF GARNET GROWTH WITH ISOCHEMICAL PHASE-DIAGRAM SECTIONS AND THE PROBLEM OF EFFECTIVE BULK-COMPOSITION , 2005 .

[11]  R. Batchelor,et al.  Petrogenetic interpretation of granitoid rock series using multicationic parameters , 1985 .

[12]  James A. D. Connolly,et al.  Computation of phase equilibria by linear programming: A tool for geodynamic modeling and its application to subduction zone decarbonation , 2005 .

[13]  A. Putnis,et al.  Replacement Processes in the Earth's Crust , 2010 .

[14]  S. Endo.  Pressure‐temperature history of titanite‐bearing eclogite from the Western Iratsu body, Sanbagawa Metamorphic Belt, Japan , 2010 .

[15]  Roger Powell,et al.  An internally consistent thermodynamic data set for phases of petrological interest , 1998 .

[16]  N. Tsuchiya,et al.  Geochemical behavior of zirconium during Cl–rich fluid or melt infiltration under upper amphibolite facies metamorphism — A case study from Brattnipene, Sør Rondane Mountains, East Antarctica , 2015 .

[17]  J. D. Connolly,et al.  Titanium in phengite: a geobarometer for high temperature eclogites , 2009 .

[18]  C. Graham,et al.  Cycling of B, Li, and LILE (K, Cs, Rb, Ba, Sr) into subduction zones: SIMS evidence from micas in high-P/T metasedimentary rocks , 2007 .

[19]  V. Pawlowsky-Glahn,et al.  Compositional data and their analysis: an introduction , 2006, Geological Society, London, Special Publications.

[20]  W. Leeman,et al.  BBe systematics in subduction-related metamorphic rocks: Characterization of the subducted component , 1993 .

[21]  S. Wallis,et al.  A re‐evaluation of eclogite facies metamorphism in SW Japan: proposal for an eclogite nappe , 2000 .

[22]  H. Iwamori,et al.  Elemental transport upon hydration of basic schists during regional metamorphism : Geochemical evidence from the Sanbagawa metamorphic belt, Japan , 2014 .

[23]  Y. Kouketsu,et al.  Calculated stabilities of sodic phases in the Sambagawa metapelites and their implications , 2011 .

[24]  M. Enami.  Pressure‐temperature path of Sanbagawa prograde metamorphism deduced from grossular zoning of garnet , 1998 .

[25]  T. Hirajima,et al.  Evidence of the lawsonite eclogite facies metamorphism from an epidote-glaucophane eclogite in the Kotsu area of the Sanbagawa belt, Japan , 2013 .

[26]  Michael I. Jordan,et al.  Latent Dirichlet Allocation , 2001, J. Mach. Learn. Res..

[27]  S. Uehara,et al.  Slow subduction and buoyant exhumation of the Sanbagawa eclogite , 2012 .

[28]  P. O'Brien,et al.  Fluid Migration above a Subducted Slab—Constraints on Amount, Pathways and Major Element Mobility from Partially Overprinted Eclogite-facies Rocks (Sesia Zone, Western Alps) , 2011 .

[29]  A. Tsuchiyama,et al.  Combined FIB microsampling and X-ray microtomography: a powerful tool for the study of tiny fluid inclusions , 2016 .

[30]  C. Sakai,et al.  XRF analyses of Sanbagawa pelitic schists in central Shikoku, Japan , 1996 .

[31]  M. Enami,et al.  Coexistence of jadeite and quartz in garnet of the Sanbagawa metapelite from the Asemi–gawa region, central Shikoku, Japan , 2014 .

[32]  S. Banno,et al.  Pressure-temperature trajectory of the Sanbagawa metamorphism deduced from garnet zoning , 1986 .

[33]  M. Aoya.  P–T–D Path of Eclogite from the Sambagawa Belt Deduced from Combination of Petrological and Microstructural Analyses , 2001 .

[34]  Kenta Yoshida,et al.  3D chemical mapping of ‘Mn–caldera shaped zoning’ garnet found from the Sanbagawa metamorphic belt of the Besshi district, SW Japan , 2015 .

[35]  S. Banno,et al.  Low-grade Progressive Metamorphism of Pelitic Schists of the Sazare area, Sanbagawa Metamorphic Terrain in central Sikoku, Japan , 1974 .

[36]  R. Powell,et al.  Calculation of partial melting equilibria in the system Na2O–CaO–K2O–FeO–MgO–Al2O3–SiO2–H2O (NCKFMASH) , 2001 .

[37]  正樹 榎並.  四国中央部別子地域・三波川帯の灰曹長石-黒雲母帯 , 1982 .

[38]  T. Kuwatani,et al.  Thermodynamic forward modeling of progressive dehydration reactions during subduction of oceanic crust under greenschist facies conditions , 2011 .

[39]  B. Yardley,et al.  Quartz, albite and diopside solubilities in H2O–NaCl and H2O–CO2 fluids at 0.5–0.9 GPa , 2001 .

[40]  A. Takasu.  Prograde and Retrograde Eclogites in the Sambagawa Metamorphic Belt, Besshi District, Japan , 1984 .

[41]  T. Ohta,et al.  Problems in compositional data analysis and their solutions , 2006 .

[42]  R. Arculus,et al.  Redistribution of trace elements during prograde metamorphism from lawsonite blueschist to eclogite facies; implications for deep subduction-zone processes , 2003 .

[43]  Yoshinori Nakanishi-Ohno,et al.  Three levels of data-driven science , 2016 .

[44]  Kenta Yoshida,et al.  Geochemical features and relative B–Li–Cl compositions of deep-origin fluids trapped in high-pressure metamorphic rocks , 2015 .

[45]  S. Banno.  Brief history of petrotectonic research on the Sanbagawa Belt, Japan , 2004 .

[46]  Donna L. Whitney,et al.  Abbreviations for names of rock-forming minerals , 2010 .

[47]  A. Takasu,et al.  Prograde eclogites from the Tonaru epidote amphibolite mass in the Sambagawa Metamorphic Belt, central Shikoku, southwest Japan , 2005 .

[48]  John Aitchison,et al.  The Statistical Analysis of Compositional Data , 1986 .

[49]  T. Holland,et al.  Mixing properties of phengitic micas and revised garnet‐phengite thermobarometers , 2002 .

[50]  H. Keppler,et al.  Solubility of rutile in subduction zone fluids, as determined by experiments in the hydrothermal diamond anvil cell , 2005 .

[51]  H. Iwamori,et al.  Transition from dehydration to hydration during exhumation of the Sanbagawa metamorphic belt, Japan, revealed by the continuous P–T path recorded in garnet and amphibole zoning , 2015, Contributions to Mineralogy and Petrology.

[52]  J. Platt,et al.  Age and early metamorphic history of the Sanbagawa belt: Lu–Hf and P–T constraints from the Western Iratsu eclogite , 2009 .

[53]  Charu C. Aggarwal,et al.  Mining Text Data , 2012, Springer US.

[54]  T. Mouri,et al.  Composite metamorphic history recorded in garnet porphyroblasts of Sambagawa metasediments in the Besshi region, central Shikoku, Southwest Japan , 2014 .

[55]  Y. Podladchikov,et al.  Volcanic arcs fed by rapid pulsed fluid flow through subducting slabs , 2012 .

[56]  J. Aitchison,et al.  Biplots of Compositional Data , 2002 .

[57]  T. Kuwatani,et al.  Classification of geochemical data based on multivariate statistical analyses: Complementary roles of cluster, principal component, and independent component analyses , 2017 .

[58]  R. Tolosana-Delgado,et al.  Some Basic Concepts of Compositional Geometry , 2005 .

[59]  N. Arndt,et al.  Progressive crustal contamination of the Bushveld Complex: evidence from Nd isotopic analyses of the cumulate rocks , 2000 .

[60]  W. Leeman,et al.  FRACTIONATION OF TRACE ELEMENTS BY SUBDUCTION-ZONE METAMORPHISM : EFFECT OF CONVERGENT-MARGIN THERMAL EVOLUTION , 1999 .

[61]  S. Wallis,et al.  Paragenesis of sodic pyroxene-bearing quartz schists: implications for the P-T history of the Sanbagawa belt , 1994 .

[62]  Kenta Yoshida,et al.  Grain size reduction due to fracturing and subsequent grain-size-sensitive creep in a lower crustal shear zone in the presence of a CO2-bearing fluid , 2017 .

[63]  Kenta Yoshida,et al.  Phase relations of lawsonite-blueschists and their role as a water-budget monitor: a case study from the Hakoishi sub-unit of the Kurosegawa belt, SW Japan , 2016 .

[64]  Worley,et al.  The effect of TiO2 and Fe2O3 on metapelitic assemblages at greenschist and amphibolite facies conditions: mineral equilibria calculations in the system K2O–FeO–MgO–Al2O3–SiO2–H2O–TiO2–Fe2O3 , 2000 .