Pub Date : 2024-04-04DOI: 10.1093/petrology/egae038
Matthew P Loocke, Jonathan E Snow
Our understanding of the processes at work in the lower crust/upper mantle transition zone during subduction initiation and early arc development has suffered from a general lack of in situ samples. Here, we present the results of petrographic and geochemical analysis of 34 samples (9 harzburgites, 13 dunites, 2 orthopyroxenites, 3 olivine-gabbros, and 7 wehrlites) collected from the inner trench wall of the Bonin Ridge, Izu-Bonin forearc. The sample suite records three main melt-rock reaction events involving melts with forearc basalt (FAB)-like, boninitic, and transitional compositions. The wehrlitic and gabbroic rocks trend towards more transitional to FAB compositions and the rest towards more boninitic compositions. The crosscutting occurrence of all three events in a single sample (wehrlite D31-106) establishes a relative timing of the events like that reported for the volcanic edifice of the Bonin Ridge, which transitioned from forearc basalt volcanism at subduction initiation (c.a., 51-52 Ma) to boninitic volcanism (c.a., 50-51 Ma) as the subduction system matured. We therefore suggest that the lower crust/upper mantle transition of the Bonin Ridge preserves a record of the transition from FAB melts created by decompression melting at subduction initiation to arc-type flux melting and boninitic volcanism thereafter. Orthopyroxenites and two anomalously fresh harzburgites from the sample suite are suggested to represent the later boninitic melts and possibly the result of hybridization between such melts and residual peridotites, respectively. Diffuse melt-rock reaction between the later boninites and/or subduction-related fluids and the earlier-formed FAB-related crust is recorded by enrichments in fluid mobile elements and depletions in first row transition metals in clinopyroxenes from a metasomatic vein in wehrlite sample D31-106. The chemistry of the wehrlitic and gabbroic clinopyroxenes suggests that they crystallized from hydrous, highly-depleted melts which lack a slab fluid signature. We thus suggest that highly-depleted melt fractions might be created early on during subduction initiation by the introduction of seawater into the proto-mantle wedge. The overall FAB-like nature of the crustal wehrlites and gabbros would suggest that most of the lower arc crust was created by forearc extension during/following subduction initiation and that later, mature arc volcanism may have contributed little or no material to the lower crust/upper mantle record in the outer forearc.
{"title":"Arc foundations and subduction initiation: Insights into the magmatic evolution of the lower crust/upper mantle of the Izu-Bonin forearc during subduction initiation","authors":"Matthew P Loocke, Jonathan E Snow","doi":"10.1093/petrology/egae038","DOIUrl":"https://doi.org/10.1093/petrology/egae038","url":null,"abstract":"Our understanding of the processes at work in the lower crust/upper mantle transition zone during subduction initiation and early arc development has suffered from a general lack of in situ samples. Here, we present the results of petrographic and geochemical analysis of 34 samples (9 harzburgites, 13 dunites, 2 orthopyroxenites, 3 olivine-gabbros, and 7 wehrlites) collected from the inner trench wall of the Bonin Ridge, Izu-Bonin forearc. The sample suite records three main melt-rock reaction events involving melts with forearc basalt (FAB)-like, boninitic, and transitional compositions. The wehrlitic and gabbroic rocks trend towards more transitional to FAB compositions and the rest towards more boninitic compositions. The crosscutting occurrence of all three events in a single sample (wehrlite D31-106) establishes a relative timing of the events like that reported for the volcanic edifice of the Bonin Ridge, which transitioned from forearc basalt volcanism at subduction initiation (c.a., 51-52 Ma) to boninitic volcanism (c.a., 50-51 Ma) as the subduction system matured. We therefore suggest that the lower crust/upper mantle transition of the Bonin Ridge preserves a record of the transition from FAB melts created by decompression melting at subduction initiation to arc-type flux melting and boninitic volcanism thereafter. Orthopyroxenites and two anomalously fresh harzburgites from the sample suite are suggested to represent the later boninitic melts and possibly the result of hybridization between such melts and residual peridotites, respectively. Diffuse melt-rock reaction between the later boninites and/or subduction-related fluids and the earlier-formed FAB-related crust is recorded by enrichments in fluid mobile elements and depletions in first row transition metals in clinopyroxenes from a metasomatic vein in wehrlite sample D31-106. The chemistry of the wehrlitic and gabbroic clinopyroxenes suggests that they crystallized from hydrous, highly-depleted melts which lack a slab fluid signature. We thus suggest that highly-depleted melt fractions might be created early on during subduction initiation by the introduction of seawater into the proto-mantle wedge. The overall FAB-like nature of the crustal wehrlites and gabbros would suggest that most of the lower arc crust was created by forearc extension during/following subduction initiation and that later, mature arc volcanism may have contributed little or no material to the lower crust/upper mantle record in the outer forearc.","PeriodicalId":16751,"journal":{"name":"Journal of Petrology","volume":"53 1","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-04-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140585430","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}
Pub Date : 2024-03-28DOI: 10.1093/petrology/egae036
Maurizio Petrelli
The present manuscript reports on the state-of-the-art and future perspectives of Machine Learning (ML) in petrology. To achieve this goal, it first introduces the basics of ML, including definitions, core concepts, and applications. Then, it starts reviewing the state-of-the-art of ML in petrology. Established applications mainly concern the so-called data-driven discovery and involve specific tasks like clustering, dimensionality reduction, classification, and regression. Among them, clustering and dimensionality reduction have been demonstrated to be valuable for decoding the chemical record stored in igneous and metamorphic phases and to enhance data visualization, respectively. Classification and regression tasks find applications, for example, in petrotectonic discrimination and geo-thermobarometry, respectively. The main core of the manuscript consists of depicting emerging trends and the future directions of ML in petrological investigations. I propose a future scenario where ML methods will progressively integrate and support established petrological methods in automating time-consuming and repetitive tasks, improving current models, and boosting discovery. In this framework, promising applications include (a) the acquisition of new multimodal petrologic data, (b) the development of data fusion techniques, physics-informed ML models, and ML-supported numerical simulations, and (c) the continuous exploration of the ML potential in petrology. To boost the contribution of ML in petrology, our main challenges are: (a) to improve the ability of ML models to capture the complexity of petrologic processes, (b) progressively link machine learning algorithms with the physical and thermodynamic nature of the investigated problems, (c) to start a collaborative effort among researchers coming from different disciplines, both in research and teaching.
本手稿报告了岩石学中机器学习(ML)的最新进展和未来展望。为实现这一目标,本文首先介绍了机器学习的基础知识,包括定义、核心概念和应用。然后,开始回顾机器学习在岩石学中的最新应用。已有的应用主要涉及所谓的数据驱动发现,并涉及聚类、降维、分类和回归等具体任务。其中,聚类和降维已被证明对解码火成岩相和变质岩相中存储的化学记录以及增强数据可视化分别具有重要价值。分类和回归任务分别应用于岩石构造判别和地质测温等领域。手稿的主要核心内容包括描绘 ML 在岩石学研究中的新兴趋势和未来方向。我提出的未来设想是,ML 方法将逐步整合并支持既有的岩石学方法,实现耗时和重复性任务的自动化,改进现有模型并促进发现。在此框架下,有前景的应用包括:(a) 获取新的多模态岩石学数据;(b) 开发数据融合技术、物理信息 ML 模型和 ML 支持的数值模拟;(c) 不断探索 ML 在岩石学中的潜力。为了提高 ML 在岩石学中的贡献,我们面临的主要挑战是(a) 提高 ML 模型捕捉岩石学过程复杂性的能力,(b) 逐步将机器学习算法与所研究问题的物理和热力学性质联系起来,(c) 启动来自不同学科的研究人员在研究和教学方面的合作努力。
{"title":"Machine Learning in Petrology: State-of-the-Art and Future Perspectives","authors":"Maurizio Petrelli","doi":"10.1093/petrology/egae036","DOIUrl":"https://doi.org/10.1093/petrology/egae036","url":null,"abstract":"The present manuscript reports on the state-of-the-art and future perspectives of Machine Learning (ML) in petrology. To achieve this goal, it first introduces the basics of ML, including definitions, core concepts, and applications. Then, it starts reviewing the state-of-the-art of ML in petrology. Established applications mainly concern the so-called data-driven discovery and involve specific tasks like clustering, dimensionality reduction, classification, and regression. Among them, clustering and dimensionality reduction have been demonstrated to be valuable for decoding the chemical record stored in igneous and metamorphic phases and to enhance data visualization, respectively. Classification and regression tasks find applications, for example, in petrotectonic discrimination and geo-thermobarometry, respectively. The main core of the manuscript consists of depicting emerging trends and the future directions of ML in petrological investigations. I propose a future scenario where ML methods will progressively integrate and support established petrological methods in automating time-consuming and repetitive tasks, improving current models, and boosting discovery. In this framework, promising applications include (a) the acquisition of new multimodal petrologic data, (b) the development of data fusion techniques, physics-informed ML models, and ML-supported numerical simulations, and (c) the continuous exploration of the ML potential in petrology. To boost the contribution of ML in petrology, our main challenges are: (a) to improve the ability of ML models to capture the complexity of petrologic processes, (b) progressively link machine learning algorithms with the physical and thermodynamic nature of the investigated problems, (c) to start a collaborative effort among researchers coming from different disciplines, both in research and teaching.","PeriodicalId":16751,"journal":{"name":"Journal of Petrology","volume":"117 1","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-03-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140322199","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}
Pub Date : 2024-03-28DOI: 10.1093/petrology/egae033
Aleksandr S Stepanov, Jovid Aminov, Sharifjon Odinaev, Farukh Sh Iskandarov, Shao-Yong Jiang, Nikolai S Karmanov
Fluorite-dominated rocks are occasionally found in association with carbonatites, but their geologic and petrologic relations are rarely reported. The Dunkeldyk area of the Pamir mountains in south-eastern Tajikistan contains dikes of distinctive rocks composed of calcite, fluorite, celestine-barite, sulfides, apatite, with minor quartz, biotite, and REE-fluorcarbonates. The dikes have sharp contacts with the host (meta-)sedimentary rocks and layering with ribbons, ranging from fluorite-bearing calcite carbonatites to fluoritites (rocks with >50% fluorite). The fluoritites are characterised by high Ca, F, Ba, Sr, REE and S coupled with anomalously low O. The geologic relations and textures suggest a magmatic origin of the dikes from melts close to calcite-fluorite eutectic that experienced nucleation-controlled differentiation during the crystallization of dikes and the formation of fluoritite cumulates in larger intrusions. The Dunkeldyk dikes demonstrate that sizable geological bodies of fluorite-dominated rocks could form from carbonate-fluoride melts originating by differentiation of alkaline silicate magmas.
以萤石为主的岩石偶尔会与碳酸盐岩伴生,但它们之间的地质学和岩石学关系却鲜有报道。塔吉克斯坦东南部帕米尔山脉的敦克尔代克地区含有由方解石、萤石、天青石-重晶石、硫化物、磷灰石以及少量石英、生物岩和稀土元素-氟碳酸盐组成的独特岩石尖峰。尖晶石与主(元)沉积岩有尖锐的接触,并呈带状分层,从含萤石的方解石碳酸盐岩到萤石岩(含>50%萤石的岩石)不等。萤石的特征是高 Ca、F、Ba、Sr、REE 和 S 以及异常低的 O。地质关系和纹理表明,这些岩峰起源于接近方解石-萤石共晶的熔体,在岩峰结晶过程中经历了成核控制的分化,并在较大的侵入体中形成了萤石积层。Dunkeldyk岩钉表明,以萤石为主的岩石的大型地质体可能是由碱性硅酸盐岩浆分异产生的碳酸盐-氟化物熔体形成的。
{"title":"Fluoritites produced by crystallization of carbonate-fluoride magma","authors":"Aleksandr S Stepanov, Jovid Aminov, Sharifjon Odinaev, Farukh Sh Iskandarov, Shao-Yong Jiang, Nikolai S Karmanov","doi":"10.1093/petrology/egae033","DOIUrl":"https://doi.org/10.1093/petrology/egae033","url":null,"abstract":"Fluorite-dominated rocks are occasionally found in association with carbonatites, but their geologic and petrologic relations are rarely reported. The Dunkeldyk area of the Pamir mountains in south-eastern Tajikistan contains dikes of distinctive rocks composed of calcite, fluorite, celestine-barite, sulfides, apatite, with minor quartz, biotite, and REE-fluorcarbonates. The dikes have sharp contacts with the host (meta-)sedimentary rocks and layering with ribbons, ranging from fluorite-bearing calcite carbonatites to fluoritites (rocks with >50% fluorite). The fluoritites are characterised by high Ca, F, Ba, Sr, REE and S coupled with anomalously low O. The geologic relations and textures suggest a magmatic origin of the dikes from melts close to calcite-fluorite eutectic that experienced nucleation-controlled differentiation during the crystallization of dikes and the formation of fluoritite cumulates in larger intrusions. The Dunkeldyk dikes demonstrate that sizable geological bodies of fluorite-dominated rocks could form from carbonate-fluoride melts originating by differentiation of alkaline silicate magmas.","PeriodicalId":16751,"journal":{"name":"Journal of Petrology","volume":"87 1","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-03-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140325884","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}
Pub Date : 2024-03-28DOI: 10.1093/petrology/egae037
J Gregory Shellnutt, Meng-Wan Yeh, Tung-Yi Lee, Yoshiyuki Iizuka, Wei-Yu Chen, M P Manu Prasanth
Mantle xenoliths hosted in volcanic rocks from the island of Lutao offer a glimpse into the nature of the mantle beneath the northern Luzon volcanic arc. The xenoliths are spinel-bearing and composed mostly of harzburgite with one lherzolite and one olivine orthopyroxenite. The olivine (Fo92.5-88.9), orthopyroxene (Mg# = 94.6-89.2), and clinopyroxene (Wo49.1-38.1En57.0-45.4Fs3.0-11.0) compositions are similar to those of abyssal peridotites. The spinel compositions are variable and can be principally divided into high-Al (Cr# < 45) and low-Al (Cr# > 45) groupings. The whole rock compositions are similar to abyssal peridotite (Al2O3 = 0.95-2.07 wt%; Mg# = 88.5-90.9) and have U-shaped chondrite normalized rare earth element patterns. The Sr-Nd isotopes of the xenoliths are broadly chondritic (87Sr/86Sri = 0.704400-0.707908; εNd(t) = 0.0-+1.5). The two-pyroxene equilibrium temperatures range from 900-1200°C with the majority of temperature estimates >1000°C. The olivine-orthopyroxene-spinel oxygen barometry estimates yielded ΔFMQ values from 0 to +2 and correspond to moderately oxidizing to oxidizing conditions. The xenoliths are likely derived from the Philippine Sea Plate lithospheric mantle that was modified by melt extraction and/or fluid enrichment processes. Trace element and isotopic mixing modeling indicate that 1% to 2% contamination by subducted South China Sea sediment can explain the Sr-Nd isotopic enrichment and Th and U elemental variability within the xenoliths assuming an initial composition similar to enriched depleted mid-ocean ridge mantle (E-DMM). The anomalously high two-pyroxene equilibrium temperatures of the Lutao xenoliths relative to other regions of the northern Luzon volcanic arc (Iraya < 1000°C) indicate that they were affected by a high temperature event that was likely a consequence of recent intra-arc rifting that occurred after collision (< 6 Ma) between the Luzon arc and the Eurasian margin.
{"title":"Sub-arc mantle heterogeneity of the Northern Luzon Volcanic Arc: mineral and whole rock compositional variability in mantle xenoliths from Lutao Island","authors":"J Gregory Shellnutt, Meng-Wan Yeh, Tung-Yi Lee, Yoshiyuki Iizuka, Wei-Yu Chen, M P Manu Prasanth","doi":"10.1093/petrology/egae037","DOIUrl":"https://doi.org/10.1093/petrology/egae037","url":null,"abstract":"Mantle xenoliths hosted in volcanic rocks from the island of Lutao offer a glimpse into the nature of the mantle beneath the northern Luzon volcanic arc. The xenoliths are spinel-bearing and composed mostly of harzburgite with one lherzolite and one olivine orthopyroxenite. The olivine (Fo92.5-88.9), orthopyroxene (Mg# = 94.6-89.2), and clinopyroxene (Wo49.1-38.1En57.0-45.4Fs3.0-11.0) compositions are similar to those of abyssal peridotites. The spinel compositions are variable and can be principally divided into high-Al (Cr# &lt; 45) and low-Al (Cr# &gt; 45) groupings. The whole rock compositions are similar to abyssal peridotite (Al2O3 = 0.95-2.07 wt%; Mg# = 88.5-90.9) and have U-shaped chondrite normalized rare earth element patterns. The Sr-Nd isotopes of the xenoliths are broadly chondritic (87Sr/86Sri = 0.704400-0.707908; εNd(t) = 0.0-+1.5). The two-pyroxene equilibrium temperatures range from 900-1200°C with the majority of temperature estimates &gt;1000°C. The olivine-orthopyroxene-spinel oxygen barometry estimates yielded ΔFMQ values from 0 to +2 and correspond to moderately oxidizing to oxidizing conditions. The xenoliths are likely derived from the Philippine Sea Plate lithospheric mantle that was modified by melt extraction and/or fluid enrichment processes. Trace element and isotopic mixing modeling indicate that 1% to 2% contamination by subducted South China Sea sediment can explain the Sr-Nd isotopic enrichment and Th and U elemental variability within the xenoliths assuming an initial composition similar to enriched depleted mid-ocean ridge mantle (E-DMM). The anomalously high two-pyroxene equilibrium temperatures of the Lutao xenoliths relative to other regions of the northern Luzon volcanic arc (Iraya &lt; 1000°C) indicate that they were affected by a high temperature event that was likely a consequence of recent intra-arc rifting that occurred after collision (&lt; 6 Ma) between the Luzon arc and the Eurasian margin.","PeriodicalId":16751,"journal":{"name":"Journal of Petrology","volume":"108 1","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-03-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140322287","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}
Pub Date : 2024-03-28DOI: 10.1093/petrology/egae034
Marian B Holness
Self-organisation in plutonic igneous rocks has been suggested to form by a variety of mechanisms including oscillatory nucleation and growth, competitive particle growth (CPG), and preferential dissolution and re-precipitation during fluid infiltration enhanced by compaction, with driving forces including reduction of the interfacial energy budget by either Ostwald ripening or because the energy of boundaries between two grains of the same mineral is less than that between two grains of different minerals. An investigation of the Stillwater inch-scale layering shows that the CPG patterning mechanism leaves a characteristic microstructural signature preserving evidence for a highly interconnected melt in textural equilibrium, and slow super- and sub-solidus cooling: such a signature is also preserved in chromite-bearing fine-scale layers in the Bushveld intrusion. The cm-scale micro-rhythmic layering of the Skaergaard intrusion, superimposed on single modally-graded layers, does not have these microstructural features. Furthermore, the energy of all relevant inter-phase grain boundaries in the Skaergaard gabbros is less than that of grain boundaries involving only one mineral, viscous compaction was not a significant process in the Skaergaard intrusion, and patterning by oscillatory nucleation and growth is precluded by the fact that the micro-rhythmic layering is superimposed on modally graded layers formed by sedimentation. A new patterning mechanism is proposed, operational only in intrusions in which the interstitial liquid of the crystal mush intersects a binode and splits into two immiscible conjugates. Cm-scale separation of the immiscible conjugate liquids in a compositionally-graded mush, due to both gravity and capillary forces, leads to layering due to differences in their wetting properties. The positive feedback required for pattern formation is due to the two immiscible conjugates predominantly crystallising the minerals which they preferentially wet.
{"title":"Self-organisation in gabbroic cumulates: a new patterning mechanism driven by differential migration of immiscible liquids in a crystal mush?","authors":"Marian B Holness","doi":"10.1093/petrology/egae034","DOIUrl":"https://doi.org/10.1093/petrology/egae034","url":null,"abstract":"Self-organisation in plutonic igneous rocks has been suggested to form by a variety of mechanisms including oscillatory nucleation and growth, competitive particle growth (CPG), and preferential dissolution and re-precipitation during fluid infiltration enhanced by compaction, with driving forces including reduction of the interfacial energy budget by either Ostwald ripening or because the energy of boundaries between two grains of the same mineral is less than that between two grains of different minerals. An investigation of the Stillwater inch-scale layering shows that the CPG patterning mechanism leaves a characteristic microstructural signature preserving evidence for a highly interconnected melt in textural equilibrium, and slow super- and sub-solidus cooling: such a signature is also preserved in chromite-bearing fine-scale layers in the Bushveld intrusion. The cm-scale micro-rhythmic layering of the Skaergaard intrusion, superimposed on single modally-graded layers, does not have these microstructural features. Furthermore, the energy of all relevant inter-phase grain boundaries in the Skaergaard gabbros is less than that of grain boundaries involving only one mineral, viscous compaction was not a significant process in the Skaergaard intrusion, and patterning by oscillatory nucleation and growth is precluded by the fact that the micro-rhythmic layering is superimposed on modally graded layers formed by sedimentation. A new patterning mechanism is proposed, operational only in intrusions in which the interstitial liquid of the crystal mush intersects a binode and splits into two immiscible conjugates. Cm-scale separation of the immiscible conjugate liquids in a compositionally-graded mush, due to both gravity and capillary forces, leads to layering due to differences in their wetting properties. The positive feedback required for pattern formation is due to the two immiscible conjugates predominantly crystallising the minerals which they preferentially wet.","PeriodicalId":16751,"journal":{"name":"Journal of Petrology","volume":"45 1","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-03-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140322193","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}
Pub Date : 2024-03-28DOI: 10.1093/petrology/egae035
Cliff S J Shaw
Clinopyroxenite xenoliths comprised of cumulus clinopyroxene ± amphibole together with intercumulus phlogopite ± olivine ± apatite ± titanite form a large part of the xenolith load in the oldest deposits of the Rockeskyllerkopf Volcanic Complex (RVC) in the West Eifel volcanic field. The xenoliths also contain xenocrysts of olivine and clinopyroxene derived from mantle peridotite and clinopyroxene from lower crustal granulite. The clinopyroxenite xenoliths are divided into five groups on the basis of their modal mineralogy and mineral compositions. Groups 1 to 4 define a continuous compositional trend indicative of fractionation of a mafic alkaline magma. Group 5 xenoliths are compositionally distinct and have been tentatively linked to high pressure crystallization of phonolitic magma within the RVC system. Thermobarometry of the group 1 – 4 xenoliths indicates that they crystallised between 1 and 4 kilobars, equivalent to a depth of 4-14 km. Group 1 -3 xenoliths all crystallised at between 1050 and 1150 °C whereas the amphibole-rich group 4 xenoliths give temperature estimates of ~ 900 °C. The clinopyroxenites share a common parent magma with clinopyroxene – phlogopite veins found in subcontinental lithospheric mantle xenoliths. However, the vein forming mama was richer in incompatible elements, in particular Zr and Hf and is interpreted to be an early formed batch of magma with the clinopyroxenites crystallising from magma derived from the same mantle, which had been depleted by the earlier phase of melting. Intrusion of magma began around 155 ky prior to the eruption of the RVC. Fe-Mg interdiffusion profiles in zoned clinopyroxene show that the magma that formed the xenoliths was present in the crust for up to 28 ky prior to the eruption. However, most samples give interdiffusion times between 1.5 and 9.9 ky. Based on xenocryst residence times and the calculated P-T conditions for clinopyroxene, there were at least seven separate batches of magma emplaced below Rockeskyllerkopf, probably as sills.
{"title":"Clinopyroxenite xenoliths record magma transport and crystallisation in the middle and upper crust: A case study from the Rockeskyllerkopf Volcanic Complex, West Eifel, Germany","authors":"Cliff S J Shaw","doi":"10.1093/petrology/egae035","DOIUrl":"https://doi.org/10.1093/petrology/egae035","url":null,"abstract":"Clinopyroxenite xenoliths comprised of cumulus clinopyroxene ± amphibole together with intercumulus phlogopite ± olivine ± apatite ± titanite form a large part of the xenolith load in the oldest deposits of the Rockeskyllerkopf Volcanic Complex (RVC) in the West Eifel volcanic field. The xenoliths also contain xenocrysts of olivine and clinopyroxene derived from mantle peridotite and clinopyroxene from lower crustal granulite. The clinopyroxenite xenoliths are divided into five groups on the basis of their modal mineralogy and mineral compositions. Groups 1 to 4 define a continuous compositional trend indicative of fractionation of a mafic alkaline magma. Group 5 xenoliths are compositionally distinct and have been tentatively linked to high pressure crystallization of phonolitic magma within the RVC system. Thermobarometry of the group 1 – 4 xenoliths indicates that they crystallised between 1 and 4 kilobars, equivalent to a depth of 4-14 km. Group 1 -3 xenoliths all crystallised at between 1050 and 1150 °C whereas the amphibole-rich group 4 xenoliths give temperature estimates of ~ 900 °C. The clinopyroxenites share a common parent magma with clinopyroxene – phlogopite veins found in subcontinental lithospheric mantle xenoliths. However, the vein forming mama was richer in incompatible elements, in particular Zr and Hf and is interpreted to be an early formed batch of magma with the clinopyroxenites crystallising from magma derived from the same mantle, which had been depleted by the earlier phase of melting. Intrusion of magma began around 155 ky prior to the eruption of the RVC. Fe-Mg interdiffusion profiles in zoned clinopyroxene show that the magma that formed the xenoliths was present in the crust for up to 28 ky prior to the eruption. However, most samples give interdiffusion times between 1.5 and 9.9 ky. Based on xenocryst residence times and the calculated P-T conditions for clinopyroxene, there were at least seven separate batches of magma emplaced below Rockeskyllerkopf, probably as sills.","PeriodicalId":16751,"journal":{"name":"Journal of Petrology","volume":"31 1","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-03-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140322291","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}
Pub Date : 2024-03-25DOI: 10.1093/petrology/egae031
Jackie M Kleinsasser, Adam C Simon, Dean Peterson, Amartya Kattemalavadi, Ian R Goan, Tobias Keller, George J Hudak, Kaitlin Koshurba
The Duluth Complex is a large mafic intrusive system located in northeastern Minnesota emplaced as part of the 1.1 Ga Midcontinent Rift. Several Fe-Ti oxide-bearing ultramafic intrusions are hosted along the Western Margin of the Duluth Complex, and are discordant bodies present in a variety of geometries, hosted in multiple rock types, and dominated by peridotite, pyroxenite, and semi-massive to massive Fe-Ti oxide rock types. Their origin has been debated, and here we present geochemical evidence and modeling that supports a purely magmatic origin for the Titac and Longnose Fe-Ti oxide-bearing ultramafic intrusions. Ilmenite and titanomagnetite textures indicate a protracted cooling process, and δ34S values of sulfides reveal little assimilation of the footwall Virginia Formation, a fine-grained pelitic unit that contains sulfide-rich bands. We model the crystallization of a hypothetical parental magma composition to the host intrusion of Longnose using Rhyolite-MELTS and demonstrate that the accumulation of Fe-Ti oxides in the discordant intrusions cannot be explained by density-driven segregation of crystallized Fe-Ti oxides. Instead, we show that the development of silicate liquid immiscibility, occurring by the unmixing of the silicate melt into conjugate Si- and Fe-rich melts, can result in the effective segregation and transportation of the Fe-rich melt. The Fe-rich melt is ~2 orders of magnitude less viscous than the Si-rich melt, allowing the Fe-rich melt to be more effectively segregated and transported in the mush regime (crystallinities > 50%). This suggests that viscosity, in addition to density, plays a significant role in forming the discordant Fe-Ti oxide-bearing ultramafic intrusions. We propose a genetic model that could also be responsible for the Fe-Ti oxide-rich layers or bands that are hosted within the igneous stratigraphy of mafic intrusions of the Duluth Complex.
德卢斯复合体是位于明尼苏达州东北部的大型岩浆侵入系统,是 1.1 Ga 中大陆裂谷的一部分。沿德卢斯复合体西缘有几处含铁钛氧化物的超基性侵入体,这些侵入体呈不和谐体,几何形态各异,赋存于多种岩石类型中,以橄榄岩、辉绿岩和半块状至块状铁钛氧化物岩石类型为主。它们的成因一直存在争议,在此,我们提出地球化学证据和模型,支持Titac和Longnose含铁钛氧化物超基性侵入体的纯岩浆成因。钛铁矿和榍石的纹理显示了一个漫长的冷却过程,而硫化物的δ34S值则显示了脚墙弗吉尼亚地层(包含富硫化物带的细粒辉长岩单元)的同化程度很低。我们使用 Rhyolite-MELTS 模拟了 Longnose 主侵入体的假定母岩浆成分的结晶过程,并证明不和谐侵入体中铁钛氧化物的累积不能用结晶铁钛氧化物的密度驱动偏析来解释。相反,我们证明了硅酸盐熔体不混合成共轭的富硅和富铁熔体所导致的硅酸盐液不溶解性的发展,可导致富铁熔体的有效偏析和运移。富Fe熔体的粘度比富Si熔体低约2个数量级,这使得富Fe熔体能够在蘑菇云状态(结晶度> 50%)下更有效地分离和运输。这表明,在形成不和谐的含铁钛氧化物超基性侵入体的过程中,除了密度之外,粘度也起着重要作用。我们提出了一个遗传模型,该模型也可能是德卢斯岩群黑云母侵入体火成岩地层中富含铁-钛氧化物层或带的成因。
{"title":"Genesis of Fe-Ti oxide-bearing ultramafic intrusions in the Duluth Complex, Minnesota, USA","authors":"Jackie M Kleinsasser, Adam C Simon, Dean Peterson, Amartya Kattemalavadi, Ian R Goan, Tobias Keller, George J Hudak, Kaitlin Koshurba","doi":"10.1093/petrology/egae031","DOIUrl":"https://doi.org/10.1093/petrology/egae031","url":null,"abstract":"The Duluth Complex is a large mafic intrusive system located in northeastern Minnesota emplaced as part of the 1.1 Ga Midcontinent Rift. Several Fe-Ti oxide-bearing ultramafic intrusions are hosted along the Western Margin of the Duluth Complex, and are discordant bodies present in a variety of geometries, hosted in multiple rock types, and dominated by peridotite, pyroxenite, and semi-massive to massive Fe-Ti oxide rock types. Their origin has been debated, and here we present geochemical evidence and modeling that supports a purely magmatic origin for the Titac and Longnose Fe-Ti oxide-bearing ultramafic intrusions. Ilmenite and titanomagnetite textures indicate a protracted cooling process, and δ34S values of sulfides reveal little assimilation of the footwall Virginia Formation, a fine-grained pelitic unit that contains sulfide-rich bands. We model the crystallization of a hypothetical parental magma composition to the host intrusion of Longnose using Rhyolite-MELTS and demonstrate that the accumulation of Fe-Ti oxides in the discordant intrusions cannot be explained by density-driven segregation of crystallized Fe-Ti oxides. Instead, we show that the development of silicate liquid immiscibility, occurring by the unmixing of the silicate melt into conjugate Si- and Fe-rich melts, can result in the effective segregation and transportation of the Fe-rich melt. The Fe-rich melt is ~2 orders of magnitude less viscous than the Si-rich melt, allowing the Fe-rich melt to be more effectively segregated and transported in the mush regime (crystallinities &gt; 50%). This suggests that viscosity, in addition to density, plays a significant role in forming the discordant Fe-Ti oxide-bearing ultramafic intrusions. We propose a genetic model that could also be responsible for the Fe-Ti oxide-rich layers or bands that are hosted within the igneous stratigraphy of mafic intrusions of the Duluth Complex.","PeriodicalId":16751,"journal":{"name":"Journal of Petrology","volume":"1 1","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-03-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140297599","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}
Pub Date : 2024-03-12DOI: 10.1093/petrology/egae023
Simon J E Large, Chetan L Nathwani, Jamie J Wilkinson, Thomas R Knott, Simon R Tapster, Yannick Buret
Subduction zone magmatism is a major control of volcanism, the generation of modern continental crust and the formation of economically important porphyry Cu-(Mo-Au) deposits. Reading the magmatic record of individual arc segments, and constraining the rates of magmatic changes, are critical in order to fully understand and quantify the processes that drive magma evolution in subduction settings during arc growth. This study focuses on the San Francisco Batholith and the Rio Blanco-Los Bronces porphyry deposit cluster in central Chile, which provides an igneous rock record over ~13.5 Myr of arc evolution. We use whole-rock geochemistry, zircon geochronology and Hf-isotope geochemistry to track changes in the crustal magmatic system of this arc segment during crustal thickening and porphyry Cu deposit formation. By combining the analytical dataset with Monte Carlo fractional crystallisation and assimilation fractional crystallisation modelling, we test a model for significant crustal involvement during magma evolution. Systematic and continuous increases in Dy/Yb, La/Yb, V/Sc and Sr/Y in the magmas over time indicate a transition in the main fractionation assemblage from plagioclase-dominated to amphibole-dominated that reflects deeper crystallisation and/or a higher melt water content. Concomitant decreases in εHf and Th/La as well as increasing Ba/Th are best explained by assimilation of progressively deeper crustal lithologies from low (Chilenia) to high Ba/Th (Cuyania) basement terranes. Our study highlights that an increasingly hydrous magma and a deepening locus of crustal magma differentiation and assimilation, driven by crustal thickening contemporaneous with increased tectonic convergence and ingression of the aseismic Juan Fernandez ridge, can account for all investigated aspects of the multi-Myr magmatic evolution leading up to the formation of the Rio Blanco-Los Bronces porphyry Cu deposits. Our findings corroborate the importance of high-pressure differentiation of hydrous magma for the formation of Andean-style porphyry deposits. Once magmas favourable for porphyry Cu mineralisation were generated in the lower crust, multiple episodes of efficient magma migration into the upper crust fed several, discrete, shallow magmatic-hydrothermal systems over ~3.5 Myr to form the world’s largest known Cu resource at Rio Blanco-Los Bronces.
俯冲带岩浆活动是火山活动、现代大陆地壳的生成以及具有重要经济价值的斑岩型铜(钼-金)矿床形成的主要控制因素。解读单个弧段的岩浆记录,并制约岩浆变化的速率,对于充分理解和量化弧长期间俯冲环境中岩浆演化的驱动过程至关重要。这项研究的重点是智利中部的旧金山岩浆岩和里奥布兰科-洛斯布朗塞斑岩矿床群,它们提供了约 13.5 Myr 的弧演化火成岩记录。我们利用全岩地球化学、锆石地质年代学和Hf-同位素地球化学来追踪该弧段在地壳增厚和斑岩铜矿床形成过程中地壳岩浆系统的变化。通过将分析数据集与蒙特卡洛分型结晶和同化分型结晶模型相结合,我们对岩浆演化过程中地壳的重要参与模型进行了测试。随着时间的推移,岩浆中的Dy/Yb、La/Yb、V/Sc和Sr/Y持续系统地增加,表明主要分馏组合从斜长石为主过渡到闪石为主,这反映了更深的结晶和/或更高的熔融水含量。εHf和Th/La的相应减少以及Ba/Th的相应增加最能解释从低Ba/Th(Chilenia)到高Ba/Th(Cuyania)基底地层岩性的逐渐深入同化。我们的研究强调,在地壳增厚的驱动下,岩浆含水量不断增加,地壳岩浆分化和同化的位置不断加深,与此同时,构造辐合加剧,无地震的胡安-费尔南德斯海脊侵入,这些因素可以解释导致 Rio Blanco-Los Bronces 斑岩铜矿床形成的多年岩浆演化的所有调查方面。我们的研究结果证实了含水岩浆高压分异对安第斯型斑岩矿床形成的重要性。一旦有利于斑岩型铜矿化的岩浆在下地壳产生,多次有效的岩浆迁移到上地壳,在大约3.5百万年的时间里,为几个离散的浅层岩浆-热液系统提供了能量,从而在Rio Blanco-Los Bronces形成了世界上已知的最大铜资源。
{"title":"Tectonic and crustal processes drive multi-million year arc magma evolution leading up to porphyry copper deposit formation in central Chile","authors":"Simon J E Large, Chetan L Nathwani, Jamie J Wilkinson, Thomas R Knott, Simon R Tapster, Yannick Buret","doi":"10.1093/petrology/egae023","DOIUrl":"https://doi.org/10.1093/petrology/egae023","url":null,"abstract":"Subduction zone magmatism is a major control of volcanism, the generation of modern continental crust and the formation of economically important porphyry Cu-(Mo-Au) deposits. Reading the magmatic record of individual arc segments, and constraining the rates of magmatic changes, are critical in order to fully understand and quantify the processes that drive magma evolution in subduction settings during arc growth. This study focuses on the San Francisco Batholith and the Rio Blanco-Los Bronces porphyry deposit cluster in central Chile, which provides an igneous rock record over ~13.5 Myr of arc evolution. We use whole-rock geochemistry, zircon geochronology and Hf-isotope geochemistry to track changes in the crustal magmatic system of this arc segment during crustal thickening and porphyry Cu deposit formation. By combining the analytical dataset with Monte Carlo fractional crystallisation and assimilation fractional crystallisation modelling, we test a model for significant crustal involvement during magma evolution. Systematic and continuous increases in Dy/Yb, La/Yb, V/Sc and Sr/Y in the magmas over time indicate a transition in the main fractionation assemblage from plagioclase-dominated to amphibole-dominated that reflects deeper crystallisation and/or a higher melt water content. Concomitant decreases in εHf and Th/La as well as increasing Ba/Th are best explained by assimilation of progressively deeper crustal lithologies from low (Chilenia) to high Ba/Th (Cuyania) basement terranes. Our study highlights that an increasingly hydrous magma and a deepening locus of crustal magma differentiation and assimilation, driven by crustal thickening contemporaneous with increased tectonic convergence and ingression of the aseismic Juan Fernandez ridge, can account for all investigated aspects of the multi-Myr magmatic evolution leading up to the formation of the Rio Blanco-Los Bronces porphyry Cu deposits. Our findings corroborate the importance of high-pressure differentiation of hydrous magma for the formation of Andean-style porphyry deposits. Once magmas favourable for porphyry Cu mineralisation were generated in the lower crust, multiple episodes of efficient magma migration into the upper crust fed several, discrete, shallow magmatic-hydrothermal systems over ~3.5 Myr to form the world’s largest known Cu resource at Rio Blanco-Los Bronces.","PeriodicalId":16751,"journal":{"name":"Journal of Petrology","volume":"10 1","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-03-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140115808","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}
Pub Date : 2024-03-04DOI: 10.1093/petrology/egae028
Marie K Takach, Wendy A Bohrson, Frank J Spera, Marco Viccaro
The nearly continuous volcanic eruption record at Mt. Etna dating back ~700 years provides an excellent opportunity to investigate the geochemical evolution of a highly active volcano. Of particular interest is elucidating the cause of selective enrichment in alkali elements (K and Rb) and 87Sr/86Sr observed in various episodes of past activity. More recently, this alkali enrichment trend started to manifest in the 17th century and accelerated after 1971, and was accompanied by an increase in the volume, frequency, and explosivity of eruptions. Explanations for this signature include recharge of alkali-enriched magmas and/or crustal contamination from the subvolcanic basement. This study quantitatively examines the role of crustal contamination in post-1971 Etnean magma compositions via hundreds of open-system phase equilibria and trace element calculations based upon whole-rock major oxides, trace elements, 87Sr/86Sr ratios, and mineral compositional data. Available pre-1971 petrochemical data are satisfactorily reproduced by fractional crystallization of a high whole-rock MgO (12–17 wt.%), Ni (135–285 ppm), and Cr (920–1330 ppm) parental magma composition that is documented in Etna’s ~4-ka fall-stratified deposit. Observed post-1971 whole-rock and glass trends and phase equilibria are reproduced via modeled assimilation of a skarn and flysch mixture, lithologies that represent the uppermost 10–15 km of sedimentary rocks beneath Etna. Notably, models show that K2O (wt.%) and Rb (ppm) behave incompatibly during partial melting of skarn/flysch. Additionally, the observed elevation of 87Sr/86Sr in post-1971 samples is consistent with the addition of radiogenic Sr from wallrock partial melts. In best-fit models, which yield observed post-1971 K2O, Rb, and 87Sr/86Sr trends, ~17% anatectic melt is assimilated and there may be a subordinate stoped wallrock component of ≤2% (percentage is relative to the starting mass of pristine magma). Previous work has shown that metasomatized spinel lherzolite and garnet pyroxenite can be melted in different proportions to reproduce long- and short-term changes observed in Etna’s geochemical products. We propose that the alkali enrichment signature observed after 1971 can be fully explained through the combination of mantle heterogeneity and crustal contamination. In particular, up to ~20% crustal input coupled with mantle heterogeneity of primitive melts explains the geochemical signals quite well. The influence of crustal contamination on post-1971 lavas is, in part, the result of frequent recharge of magmas that thermally primed the middle to upper crust and enhanced its partial melting.
{"title":"The Role of Crustal Contamination Throughout the 1329–2005 CE Eruptive Record of Mt. Etna Volcano, Italy","authors":"Marie K Takach, Wendy A Bohrson, Frank J Spera, Marco Viccaro","doi":"10.1093/petrology/egae028","DOIUrl":"https://doi.org/10.1093/petrology/egae028","url":null,"abstract":"The nearly continuous volcanic eruption record at Mt. Etna dating back ~700 years provides an excellent opportunity to investigate the geochemical evolution of a highly active volcano. Of particular interest is elucidating the cause of selective enrichment in alkali elements (K and Rb) and 87Sr/86Sr observed in various episodes of past activity. More recently, this alkali enrichment trend started to manifest in the 17th century and accelerated after 1971, and was accompanied by an increase in the volume, frequency, and explosivity of eruptions. Explanations for this signature include recharge of alkali-enriched magmas and/or crustal contamination from the subvolcanic basement. This study quantitatively examines the role of crustal contamination in post-1971 Etnean magma compositions via hundreds of open-system phase equilibria and trace element calculations based upon whole-rock major oxides, trace elements, 87Sr/86Sr ratios, and mineral compositional data. Available pre-1971 petrochemical data are satisfactorily reproduced by fractional crystallization of a high whole-rock MgO (12–17 wt.%), Ni (135–285 ppm), and Cr (920–1330 ppm) parental magma composition that is documented in Etna’s ~4-ka fall-stratified deposit. Observed post-1971 whole-rock and glass trends and phase equilibria are reproduced via modeled assimilation of a skarn and flysch mixture, lithologies that represent the uppermost 10–15 km of sedimentary rocks beneath Etna. Notably, models show that K2O (wt.%) and Rb (ppm) behave incompatibly during partial melting of skarn/flysch. Additionally, the observed elevation of 87Sr/86Sr in post-1971 samples is consistent with the addition of radiogenic Sr from wallrock partial melts. In best-fit models, which yield observed post-1971 K2O, Rb, and 87Sr/86Sr trends, ~17% anatectic melt is assimilated and there may be a subordinate stoped wallrock component of ≤2% (percentage is relative to the starting mass of pristine magma). Previous work has shown that metasomatized spinel lherzolite and garnet pyroxenite can be melted in different proportions to reproduce long- and short-term changes observed in Etna’s geochemical products. We propose that the alkali enrichment signature observed after 1971 can be fully explained through the combination of mantle heterogeneity and crustal contamination. In particular, up to ~20% crustal input coupled with mantle heterogeneity of primitive melts explains the geochemical signals quite well. The influence of crustal contamination on post-1971 lavas is, in part, the result of frequent recharge of magmas that thermally primed the middle to upper crust and enhanced its partial melting.","PeriodicalId":16751,"journal":{"name":"Journal of Petrology","volume":"63 1","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-03-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140074593","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}
Pub Date : 2024-03-02DOI: 10.1093/petrology/egae020
Gregor Weber, Jon Blundy
Experimentally calibrated models to recover pressures and temperatures of magmas, are widely used in igneous petrology. However, large errors, especially in barometry, limit the capacity of these models to resolve the architecture of crustal igneous systems. Here we apply machine learning to a large experimental database to calibrate new regression models that recover P-T of magmas based on melt composition plus associated phase assemblage. The method is applicable to compositions from basalt to rhyolite, pressures from 0.2 to 15 kbar, and temperatures of 675-1400°C. Testing and optimisation of the model with a filter that removes estimates with standard deviation above the 50th percentile show that pressures can be recovered with root-mean-square-error (RMSE) of 1.1-1.3 kbar and errors on temperature estimates of 21°C. Our findings demonstrate that, given constraints on the coexisting mineral assemblage melt chemistry, is a reliable recorder of magmatic variables. This is a consequence of the relatively low thermodynamic variance of natural magma compositions despite their relatively large number of constituent oxide components. We apply our model to two contrasting cases with well-constrained geophysical information: Mount St. Helens volcano (USA), and Askja caldera in Iceland. Dacite whole-rocks from Mount St Helens erupted 1980-1986, inferred to represent liquids extracted from cpx-hbl-opx-plag-mt-ilm mush, yield melt extraction source pressures of 5.1-6.7 kbar in excellent agreement with geophysical constraints. Melt inclusions and matrix glasses record lower pressures (0.7-3.8 kbar), consistent with magma crystallisation within the upper reaches of the imaged geophysical anomaly and during ascent. Magma reservoir depth estimates for historical eruptions from Askja match the location of seismic wave speed anomalies. Vp/Vs anomalies at 5-10 km depth correspond to hot (~990°C) rhyolite source regions, while basaltic magmas (~1120°C) were stored at 7-17 km depth under the caldera. These examples illustrate how our model can link petrology and geophysics to better constrain the architecture of volcanic feeding systems. Our model (MagMaTaB) is accessible through a user-friendly web application (https://igdrasil.shinyapps.io/MagmaTaBv4/).
{"title":"A machine learning-based thermobarometer for magmatic liquids","authors":"Gregor Weber, Jon Blundy","doi":"10.1093/petrology/egae020","DOIUrl":"https://doi.org/10.1093/petrology/egae020","url":null,"abstract":"Experimentally calibrated models to recover pressures and temperatures of magmas, are widely used in igneous petrology. However, large errors, especially in barometry, limit the capacity of these models to resolve the architecture of crustal igneous systems. Here we apply machine learning to a large experimental database to calibrate new regression models that recover P-T of magmas based on melt composition plus associated phase assemblage. The method is applicable to compositions from basalt to rhyolite, pressures from 0.2 to 15 kbar, and temperatures of 675-1400°C. Testing and optimisation of the model with a filter that removes estimates with standard deviation above the 50th percentile show that pressures can be recovered with root-mean-square-error (RMSE) of 1.1-1.3 kbar and errors on temperature estimates of 21°C. Our findings demonstrate that, given constraints on the coexisting mineral assemblage melt chemistry, is a reliable recorder of magmatic variables. This is a consequence of the relatively low thermodynamic variance of natural magma compositions despite their relatively large number of constituent oxide components. We apply our model to two contrasting cases with well-constrained geophysical information: Mount St. Helens volcano (USA), and Askja caldera in Iceland. Dacite whole-rocks from Mount St Helens erupted 1980-1986, inferred to represent liquids extracted from cpx-hbl-opx-plag-mt-ilm mush, yield melt extraction source pressures of 5.1-6.7 kbar in excellent agreement with geophysical constraints. Melt inclusions and matrix glasses record lower pressures (0.7-3.8 kbar), consistent with magma crystallisation within the upper reaches of the imaged geophysical anomaly and during ascent. Magma reservoir depth estimates for historical eruptions from Askja match the location of seismic wave speed anomalies. Vp/Vs anomalies at 5-10 km depth correspond to hot (~990°C) rhyolite source regions, while basaltic magmas (~1120°C) were stored at 7-17 km depth under the caldera. These examples illustrate how our model can link petrology and geophysics to better constrain the architecture of volcanic feeding systems. Our model (MagMaTaB) is accessible through a user-friendly web application (https://igdrasil.shinyapps.io/MagmaTaBv4/).","PeriodicalId":16751,"journal":{"name":"Journal of Petrology","volume":"37 1","pages":""},"PeriodicalIF":3.9,"publicationDate":"2024-03-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140018042","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}