{"title":"Thermal Buffering, Heat Advection and Crustal Thinning in the Ryoke Metamorphic Complex, Yanai, Southwest Japan","authors":"Takeshi Ikeda, Kazuhiro Miyazaki, Yu Sugawara, Risa Hayashi","doi":"10.1111/jmg.70035","DOIUrl":"10.1111/jmg.70035","url":null,"abstract":"<p>The metamorphic zoning and geothermobarometry of the Ryoke metamorphic complex in the Yanai area, southwest Japan, show that its thermobaric structure was buffered by the dehydration melting of biotite. The temperatures over most of the area covered by the three high-grade zones (8.5- to 19.0-km depth) are consistent with those of the dehydration melting reaction and indicate a geothermal gradient of 9.2°C/km. The deepest area (> 19-km depth) of the highest grade zone has a domal structure and records temperatures that are > 50°C higher than the melting reaction. In contrast, the low-grade zones (< 8.5-km depth) record a geothermal gradient of > 100°C/km. The spacing of the isobars indicates that the crust was vertically thinned at all these depths. Thermal modelling suggests that the formation of this thermobaric structure requires thermal buffering by the dehydration melting reactions of muscovite and biotite and also advective heat transfer related to the upward migration of melt produced by these reactions. The melt migration through the thermally buffered middle crust enabled heat to advect efficiently to the shallower crust and modify the conductive geothermal gradient. The spatial–temporal scale of this thermal event in the study area is comparable to that of coeval and present-day volcanism, suggesting a close link between deep crustal and surface processes.</p>","PeriodicalId":16472,"journal":{"name":"Journal of Metamorphic Geology","volume":"44 4","pages":"287-310"},"PeriodicalIF":3.4,"publicationDate":"2026-04-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/jmg.70035","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147715081","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Timothy Chapman, Geoffrey L. Clarke, Sandra Piazolo
{"title":"The Contribution of Volume Change in Preserving High-P Metastable Assemblages","authors":"Timothy Chapman, Geoffrey L. Clarke, Sandra Piazolo","doi":"10.1111/jmg.70028","DOIUrl":"https://doi.org/10.1111/jmg.70028","url":null,"abstract":"<div>\u0000 \u0000 <p>Incomplete retrograde reaction of sodic–calcic pyroxene occurred during the exhumation of Eocene eclogite and blueschist in NE New Caledonia. Sodic–calcic pyroxene is well-preserved in reaction relationships with paragonite and structurally late albite, in conflict with the prediction that it should transform efficiently during decompression to albite, paragonite, glaucophane and quartz. The microstructural distribution of albite accompanying crenulation cleavage development is consistent with albite growth having been restricted by space-filling requirements. In contrast, domains of high strain enable the accommodation of volume and fluid infiltration, increasing the size of the equilibration volume leading to the complete conversion of sodic–calcic clinopyroxene to albite. Despite decompression-driven dehydration, limited volume accommodation presents a plausible mechanism to metastably preserve high-<i>P</i> (up to 18 kbar) sodic–calcic clinopyroxene-bearing assemblages to lower-<i>P</i> conditions (~10 kbar) that should stabilise albite. Efficient retrogression is facilitated by dynamic recrystallisation inducing dilation.</p>\u0000 </div>","PeriodicalId":16472,"journal":{"name":"Journal of Metamorphic Geology","volume":"44 3","pages":"177-199"},"PeriodicalIF":3.4,"publicationDate":"2026-03-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147570153","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Tim Sotelšek, Marian Janák, Sašo Šturm, Sorour Semsari Parapari, Nik Gračanin, Kristian Radan, Mehmet Ali Gülgün, Stanislava Milovská, Mirijam Vrabec
{"title":"Decoding Metamorphic Nanodiamonds Formation Through Their Internal Structures","authors":"Tim Sotelšek, Marian Janák, Sašo Šturm, Sorour Semsari Parapari, Nik Gračanin, Kristian Radan, Mehmet Ali Gülgün, Stanislava Milovská, Mirijam Vrabec","doi":"10.1111/jmg.70032","DOIUrl":"https://doi.org/10.1111/jmg.70032","url":null,"abstract":"<p>Metamorphic diamonds offer insights into Earth's evolution, modulating the global carbon cycle through subduction into the mantle and exhumation. However, the formation and internal structure of minute diamonds remain poorly understood. Here, we study the internal structures of diamond-bearing inclusions from the Eastern Alps using state-of-the-art analytical techniques. Atomically resolved images revealed intricate structures, including well-faceted monocrystalline diamonds, complex nanocrystalline clusters of intergrown diamond and graphite and amorphous matter. Associated fluids contain CO<sub>2</sub>, CH<sub>4</sub> and H<sub>2</sub>O, indicating formation from C–O–H fluids. The proposed model of C–O–H fluid behaviour suggests monocrystalline diamond precipitation at the onset of retrograde metamorphism. Amorphous layers formed through interaction between the trapped fluid and the host garnet during the further exhumation. When the fluid reached the P–T conditions of the transition zone between the diamond and graphite stability fields, polycrystalline clusters formed. For the first time, we can see a detailed record of metamorphic fluids evolution, preserved in the internal structures of microsized diamond inclusions.</p>","PeriodicalId":16472,"journal":{"name":"Journal of Metamorphic Geology","volume":"44 3","pages":"226-245"},"PeriodicalIF":3.4,"publicationDate":"2026-03-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/jmg.70032","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147563170","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Shuang Tang, Jinrui Zhang, Chris Yakymchuk, Hang Chu, Chunjing Wei, Haoda Bian, Yidan Wang, Yuan Li
{"title":"Pressure–Temperature–Time Evolution of a Polymetamorphic Mica Schist With Two-Stage Garnet Growth From the Baoyintu Group of Central Inner Mongolia, China","authors":"Shuang Tang, Jinrui Zhang, Chris Yakymchuk, Hang Chu, Chunjing Wei, Haoda Bian, Yidan Wang, Yuan Li","doi":"10.1111/jmg.70033","DOIUrl":"https://doi.org/10.1111/jmg.70033","url":null,"abstract":"<div>\u0000 \u0000 <p>Understanding metamorphic overprinting is crucial for identifying transitions in tectono-thermal regimes during orogenesis. Here, we report the first definitive example of polymetamorphism in metapelites of the Baoyintu Group, central Inner Mongolia. Our analyses included petrographic observations, phase equilibrium modelling, as well as U–Pb and Rb–Sr geochronology. Garnet–white mica schist exhibits two stages of garnet growth: an early-stage assemblage of chlorite + garnet (I) with epidote + phengite + quartz + rutile (I) inclusions and a late-stage assemblage of chlorite + garnet (II) + muscovite + quartz + ilmenite + staurolite + rutile (II) in the matrix. Phase equilibrium modelling yields pressures of 16–22 kbar and temperatures of 400°C–500°C for the early-stage assemblage corresponding to blueschist facies metamorphism, and 6–7 kbar and ~550°C for the late-stage assemblage corresponding to amphibolite facies metamorphism. A garnet–two mica schist contains a mineral assemblage similar to the late-stage assemblage in garnet–white mica schist, yielding peak pressures of 4–5 kbar and temperatures of 570°C–580°C. Zircon and rutile U–Pb geochronology of the metapelite samples yielded two groups of metamorphic ages (c. 384 and c. 259–236 Ma). In situ mica Rb–Sr geochronology of the garnet-two mica schist yielded an age of 215 Ma interpreted to record cooling after the Permian–Triassic metamorphic peak temperature. These results reveal that medium- and low-dP/dT metamorphism occurred during the late Permian–early Triassic that overprinted Devonian high-dP/dT metamorphism. Based on previous studies of magmatism in the region, we relate the two stages of metamorphism to subduction of Paleo-Asian oceanic lithosphere (Devonian, c. 384 Ma) and back-arc extension (Permian–Triassic, 259–236 Ma). This contribution fills a critical gap in our understanding of the tectono-thermal transition history of subduction within the Paleo-Asian Ocean and provides a robust framework for interpreting polymetamorphism in accretionary orogens.</p>\u0000 </div>","PeriodicalId":16472,"journal":{"name":"Journal of Metamorphic Geology","volume":"44 3","pages":"246-263"},"PeriodicalIF":3.4,"publicationDate":"2026-03-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147563220","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Małgorzata Nowak, Lucie Tajčmanová, Marcin Dąbrowski, Iris Buisman, Jacek Szczepański
{"title":"Coesite Discovery in Eclogites Confirms UHP Metamorphism in the Orlica-Śnieżnik Dome (SW Poland)","authors":"Małgorzata Nowak, Lucie Tajčmanová, Marcin Dąbrowski, Iris Buisman, Jacek Szczepański","doi":"10.1111/jmg.70029","DOIUrl":"https://doi.org/10.1111/jmg.70029","url":null,"abstract":"<p>Eclogite lenses are exposed within the orthogneiss-dominated core of the Orlica-Śnieżnik Dome in the Sudetes, which forms the northeastern margin of the Bohemian Massif (Variscan Belt of Central Europe). The presence of coesite inclusions in garnet and omphacite confirms that these eclogites underwent ultrahigh-pressure metamorphism. This interpretation is further supported by phase equilibria modelling, which indicates peak-pressure metamorphic conditions of 2.9–3.2 GPa and 750°C–830°C. The results are consistent with estimates derived from conventional geothermobarometry and Zr-in-rutile thermometry applied to rutile inclusions in garnet. Based on quartz-in-garnet elastic barometry, a maximum entrapment pressure of approximately 2.0 GPa is obtained. We interpret this discrepancy as a result of viscous relaxation of garnet at high temperature. The first stage of re-equilibration during decompression occurred at a pressure of 2.0–2.2 GPa and a temperature of 680°C–770°C. The observed rock associations exhibit similarities to other UHP occurrences within the Saxo-Thuringian Zone, suggesting a comparable exhumation mechanism. This likely involved initial buoyancy-driven exhumation within a subduction channel, followed by crustal-scale folding. Furthermore, the maximum pressure recorded in the eclogites may partly reflect nonlithostatic components, such as transient pressure variations arising from rheological heterogeneity between the eclogites and their host rocks.</p>","PeriodicalId":16472,"journal":{"name":"Journal of Metamorphic Geology","volume":"44 3","pages":"200-225"},"PeriodicalIF":3.4,"publicationDate":"2026-03-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/jmg.70029","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147562510","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Featured Cover","authors":"T. Mackay-Champion, I. P. Cawood","doi":"10.1111/jmg.70031","DOIUrl":"https://doi.org/10.1111/jmg.70031","url":null,"abstract":"<p>The cover image is based on the article <i>An Inverse Method for Quantifying Petrological Parameters and Uncertainty in Phase Equilibrium Modelling</i> by Ian Cawood and T. Mackay-Champion, https://doi.org/10.1111/jmg.70016.\u0000 \u0000 <figure>\u0000 <div><picture>\u0000 <source></source></picture><p></p>\u0000 </div>\u0000 </figure></p>","PeriodicalId":16472,"journal":{"name":"Journal of Metamorphic Geology","volume":"44 1","pages":""},"PeriodicalIF":3.4,"publicationDate":"2026-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/jmg.70031","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145969912","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Location, Location, Location: Monazite Behaviour During UHT Metamorphism and Melt Crystallization","authors":"Chris Yakymchuk, Jessica Gareau, Michael Williams","doi":"10.1111/jmg.70030","DOIUrl":"https://doi.org/10.1111/jmg.70030","url":null,"abstract":"<p>Monazite is a robust mineral for recording the suprasolidus evolution of migmatites and granulites. However, monazite commonly has diverse compositions and yields variable dates in a single sample; understanding the controls on monazite behaviour and composition during partial melting and melt crystallization are not always straightforward. Here, we integrate in situ monazite petrochronology from sapphirine–quartz granulites in the Arequipa Massif in Peru with an equilibrium model of monazite behaviour. Monazite from leucocratic microdomains—inferred to represent crystallized remnants of melt—is generally older and enriched in Th relative to monazite in melanosome (i.e., residual) microdomains. Enrichment and depletion of Y are unrelated to the presence or absence of garnet in these samples. An equilibrium model of monazite crystallization accounts for the decrease in Th content in monazite with decreasing date but cannot reproduce the wide range of measured Y concentrations and europium anomalies in monazite in the sapphirine–quartz granulites. We suggest that the microstructural setting of monazite is an important control on trace element composition and that whole-rock equilibration of Eu/Eu* and Y (and the HREE) is unlikely. Monazite proximal to the principal sources and sinks of Y (garnet) and Eu (feldspar) may serve as a monitor of the behaviour of these minerals whereas monazite distal from these minerals are a monitor of evolving melt composition during crystallization. Monazite have similar dates to zircon from the same rocks, but monazite in melanocratic domains can be younger than zircon and this monazite may have been affected by late fluid ingress whereas monazite in leucosome was not. We emphasize the importance of in situ analysis of monazite and illustrate the limitations of a whole-rock equilibrium approach to understanding accessory mineral growth and compositions in anatectic systems.</p>","PeriodicalId":16472,"journal":{"name":"Journal of Metamorphic Geology","volume":"44 2","pages":"158-175"},"PeriodicalIF":3.4,"publicationDate":"2025-12-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/jmg.70030","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146148030","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Metamorphic Evolution of Neoarchean Amphibolites in the Hongshanzi Area, Yinshan Block, North China Craton: Insights Into a Drip/Sagduction Regime","authors":"Hexiang Chang, Chunjing Wei, Bin Wang","doi":"10.1111/jmg.70018","DOIUrl":"https://doi.org/10.1111/jmg.70018","url":null,"abstract":"<div>\u0000 \u0000 <p>The metamorphism of amphibolite facies Archean supracrustal rocks has variably been interpreted to relate to subduction–collision or drip/sagduction tectonics. This paper presents a systematic metamorphic study of supracrustal amphibolites from the Hongshanzi area in the Yinshan Block, North China Craton, to decipher their metamorphic evolution and tectonic implications. The amphibolites occur as rafts within host TTGs (trondhjemite–tonalite–granodiorite) and are metamorphosed under lower to upper amphibolite facies conditions. Two samples of lower-amphibolite facies rocks record one episode of metamorphism that exhibits clockwise <i>P–T</i> paths involving pre-<i>P</i><sub>max</sub> burial heating and post-<i>P</i><sub>max</sub> decompressional heating and subsequent cooling. The <i>P</i><sub>max</sub> conditions are constrained to ~595°C, 8.5–9 kbar and ~650°C, 8 kbar, using isopleths of the maximum Ti content in amphibole and minimum <i>X</i><sub>An</sub> (anorthite) in plagioclase in <i>P–T</i> pseudosections. This approach is quite effective for the <i>P–T</i> estimation especially of subsolidus amphibolites. The mAl (Al in M2 site) and tAl (total Al) in amphibole have high potentials to be geobarometers, but they are commonly modified during decompression and tend to underestimate pressures for amphibolites. One sample of upper-amphibolite facies records two episodes of metamorphism. The early episode is characterized by a decompression and cooling evolution from the peak at 790°C–830°C, 7–8 kbar, defined using Ti and tAl in amphibole in <i>P–T</i> and <i>T–M</i> (TiO<sub>2</sub>) pseudosections because the rock is unsaturated with respect to TiO<sub>2</sub>. The later metamorphic event shows a clockwise <i>P–T</i> path and peak condition of 560°C–600°C, 6–7.5 kbar. Zircon dating yields a crystallization age of 2.53–2.52 Ga and three metamorphic ages at 2.52–2.51 Ga, ~2.46 Ga and 1.86–1.82 Ga that are interpreted to date amphibolite facies metamorphism, modification and the overprinting event, respectively. The Neoarchean metamorphism is characterized by variable apparent geothermal gradients from 18°C/km to 29°C/km that occurred coeval with the regional TTG magmatism and that was more likely correlated with a drip or sagduction regime. In addition, the Yinshan Block was involved in the Paleoproterozoic orogeny at 1.86–1.82 Ga.</p>\u0000 </div>","PeriodicalId":16472,"journal":{"name":"Journal of Metamorphic Geology","volume":"44 2","pages":"138-157"},"PeriodicalIF":3.4,"publicationDate":"2025-12-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146139371","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"An Inverse Method for Quantifying Petrological Parameters and Uncertainty in Phase Equilibrium Modelling","authors":"T. Mackay-Champion, I. P. Cawood","doi":"10.1111/jmg.70016","DOIUrl":"https://doi.org/10.1111/jmg.70016","url":null,"abstract":"<p>Phase equilibrium modelling offers a powerful quantitative framework for understanding petrological processes. Yet, many studies still rely on qualitative comparisons between natural datasets and these forward modelled predictions to constrain model parameters, commonly pressure–temperature (<i>P</i>–<i>T</i>) conditions. Compounding this, uncertainties from the observed data or within the modelled predictions are rarely quantified, limiting confidence in the estimated <i>P</i>–<i>T</i> conditions and resulting petrological interpretations. We introduce LinaForma, an inverse modelling workflow that determines best-fit <i>P</i>–<i>T</i> conditions (or other petrological parameters) and their associated uncertainties for a given rock system by minimizing the misfit between observed data (e.g., mineral compositions or modal proportions) and their forward modelled predictions. Uncertainty is quantified by resampling of the observed data with replacement. Diagnostic metrics identify poorly performing variables and assess the sensitivity of the inversion result to variable uncertainty. Applied to an amphibolite-facies pelite and metabasite from the Greater Himalayan Sequence (Zanskar Himalaya, NW India), the approach proves effective across contrasting model systems that require different sets of solution models and variables, and it produces <i>P</i>–<i>T</i> estimates consistent with classical thermobarometry. The workflow offers several advantages: compatibility with outputs from any forward modelling software; flexible variable selection; systematic grid-search inversion in multidimensional space; a robust L1-norm misfit function resistant to outliers; and sensitivity and uncertainty analysis via bootstrap resampling. Limitations include the increasing computational demands for high-dimensional grids (<i>N</i> > 2) and the absence of explicit quantification of uncertainties inherited from the thermodynamic dataset and solution models. Alongside other emerging quantitative methods, LinaForma enables petrologists to make more informed interpretations of complex metamorphic systems and target improvements to thermodynamic datasets.</p>","PeriodicalId":16472,"journal":{"name":"Journal of Metamorphic Geology","volume":"44 1","pages":"70-84"},"PeriodicalIF":3.4,"publicationDate":"2025-12-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/jmg.70016","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145969570","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Jan Schönig, Thomas Zack, Delia Rösel, Raimon Tolosana-Delgado, Horst Marschall, Hilmar von Eynatten, Nils Keno Lünsdorf, Mona Lueder, Matthias Konrad-Schmolke, Jesse B. Walters
{"title":"Chemical Discrimination of Rutile From Different Metamorphic Source Rocks: Dealing With Hierarchical Data Structures in Random Forest Classification","authors":"Jan Schönig, Thomas Zack, Delia Rösel, Raimon Tolosana-Delgado, Horst Marschall, Hilmar von Eynatten, Nils Keno Lünsdorf, Mona Lueder, Matthias Konrad-Schmolke, Jesse B. Walters","doi":"10.1111/jmg.70017","DOIUrl":"https://doi.org/10.1111/jmg.70017","url":null,"abstract":"<p>Rutile provides a wealth of petrochronological information in metamorphic geology and due to its high stability during processes of the sedimentary cycle, rutile takes a special position in sedimentary provenance analysis. Besides being one of the classical minerals datable using the U–Pb system, rutile incorporates a broad range of trace elements, many of those being incompatible in most of the common metamorphic minerals. Although numerous multivariate statistical or machine-learning tools are available, current rutile discrimination schemes suffer from focusing on uni- and bivariate approaches leading to large compositional overlaps. Here we compiled and enlarged a dataset of 2335 rutile trace-element analyses (1646 new analyses) from 110 metamorphic rock samples of 48 localities covering a wide range of pressure–temperature conditions. After showing that the subsampling and testing strategy of the classical random forest algorithm is inappropriate for such hierarchical data structures, we introduce a modified version (random forests for hierarchically structured data) which provides realistic and generalized error estimates, improving hyperparameter tuning and performance. By applying this concept, we present a novel and multivariate rutile discrimination scheme, using the concentrations of 16 elements. The model correctly predicts the source rock composition (felsic versus mafic) in ~89% of the cases and the metamorphic gradient (≤ 350 °C/GPa versus > 350 °C/GPa) in ~84%. Combined with U–Pb dating, this will enable time-resolved insights into the geodynamic evolution of the hinterland, taking rutile provenance analysis to the next level. This approach furthermore highlights specific elements and ratios — previously understudied — that reflect bulk-rock and partitioning effects, making them valuable for future investigation of metamorphic processes. For the ease of application, a user-friendly “rutileRF-HSD” web application is provided.</p>","PeriodicalId":16472,"journal":{"name":"Journal of Metamorphic Geology","volume":"44 2","pages":"112-137"},"PeriodicalIF":3.4,"publicationDate":"2025-11-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/jmg.70017","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146140114","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}